ระบบสารสนเทศทางภูมิศาสตร์

มีวิธีคำนวณระยะทางไปยังจุดที่ใกล้ที่สุดใน QGIS หรือไม่? ใน ArcGIS มีฟังก์ชั่นที่เรียกว่าEuclidean Distanceซึ่งสามารถนำเวกเตอร์อินพุตหรือแรสเตอร์และสร้างแรสเตอร์ของระยะทาง มีฟังก์ชั่นที่คล้ายกันใน QGIS หรือไม่? ฉันหาไม่พบในเอกสารประกอบ ฉันสำรวจRasterฟังก์ชั่นใน QGIS และHub DistanceคำนวณจากMMGIS (สำหรับโซลูชันที่ไม่ใช่แรสเตอร์) แต่ฉันก็ไม่สามารถคิดได้ว่าจะใช้วิธีนี้อย่างไร หากวิธีเวกเตอร์เป็นวิธีแก้ปัญหาเดียวนั่นก็อาจใช้ได้เช่นกัน

16 qgis  raster  distance  point 

ฟังก์ชั่นสำหรับรับ SRID จากฟิลด์เรขาคณิตคืออะไร ฉันได้แทรก shapefile ลงในตาราง SQL Server ด้วยฟิลด์ Geometry ที่กำหนดไว้ ฉันรู้ว่าการฉายภาพของ shapefile คือ Swreff99_1800 (EPSG: 3011) ตอนนี้คำถามของฉันคือฉันจะรับ SRID นี้ได้อย่างไร: 3011 โดยใช้วิธี OGC ของฟิลด์เรขาคณิต

16 geometry  sql-server  srid 

ฉันใหม่จริง ๆ กับ Python และฉันต้องการทราบว่ามีวิธีที่รวดเร็วในการรับค่าเซลล์ของพิกเซลแรสเตอร์ทีละพิกเซลและพิกัด (พิกัด XY แผนที่ของศูนย์กลางของแต่ละพิกเซล) โดยใช้ Python ใน ArcGIS 10 หรือไม่ เพื่ออธิบายเพิ่มเติมนี้ฉันต้องได้รับ map X, map Y และค่าเซลล์ของพิกเซลแรกและกำหนดค่าทั้งสามให้กับตัวแปรสามตัวและทำซ้ำขั้นตอนนี้สำหรับส่วนที่เหลือของพิกเซลอื่น ๆ (วนรอบแรสเตอร์ทั้งหมด) ฉันคิดว่าฉันต้องอธิบายคำถามของฉันให้มากขึ้น ปัญหาคือฉันต้องได้รับตำแหน่ง XY ของพิกเซลของแรสเตอร์แรกและรับค่าเซลล์ของแรสเตอร์อื่น ๆ ที่สอดคล้องกับตำแหน่ง XY นั้น กระบวนการนี้ควรวนซ้ำทุกพิกเซลของแรสเตอร์แรกโดยไม่ต้องสร้างเชพไฟล์จุดกึ่งกลางใด ๆ เนื่องจากมันต้องใช้เวลามากจริงๆเพราะฉันต้องจัดการแรสเตอร์ที่มีเกือบ 8 พันล้านพิกเซล นอกจากนี้ฉันต้องทำสิ่งนี้โดยใช้ Python ใน ArcGIS 10 @ James: ขอบคุณมากสำหรับคำแนะนำของคุณ ใช่สิ่งนี้จะใช้ได้กับแรสเตอร์เดียว แต่ฉันต้องเก็บค่าเซลล์สำหรับแรสเตอร์อื่น ๆ ด้วยเช่นกัน ปัญหาคือหลังจากได้รับพิกัด X และ …

16 arcgis-10.0  python  raster 

คำกล่าวนี้เกี่ยวกับว่า 80% ของข้อมูลทั้งหมดมีองค์ประกอบเชิงพื้นที่ได้รับการรวบรวมทั่วทุกแห่ง รู้เบื้องต้นเกี่ยวกับคลาส GIS, คำอธิบายผลิตภัณฑ์, การพูดคุยและอื่น ๆ ใครจะรู้บ้างว่าใครเป็นคนพูดก่อน หรือลิงค์ไปยังบทความต้นฉบับจะดี

16 history 

ขออภัยหากคำถามต่อไปนี้ค่อนข้างงี่เง่า แต่ฉันเพิ่งใหม่กับสิ่ง GIS ทั้งหมดนี้ ฉันพยายามแปลงภาพ geoTiff ที่คาดการณ์ไว้เป็น WGS84 โดยใช้ gdal ใน python ฉันได้พบโพสต์ที่แสดงกระบวนการในการแปลงคะแนนภายใน GeoTiff ที่ฉายโดยใช้บางอย่างที่คล้ายกับต่อไปนี้: from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], …

16 python  gdal  coordinate-system  wgs84 

มือใหม่ที่จะแลกเปลี่ยนสแต็คที่นี่ ขณะนี้ฉันกำลังใช้ตัวตัดกระเบื้องอัตโนมัติที่นี่เพื่อสร้างแผ่นกระเบื้อง 256x256 มาตรฐานที่จำเป็นในการสร้างประเภทแผนที่ที่กำหนดเองสำหรับแผนที่ของฉัน สคริปต์ตัวตัดไทล์ของแผนที่รู้สึกค่อนข้างงุ่มง่ามในการใช้งานและฉันใช้ RAM ขนาด 4GB ของเครื่องของฉันทั้งหมดเพื่อเปิด photoshop และใช้มัน นั่นและมันน่ารำคาญมากที่จะวางรูปภาพของฉันด้วยพื้นที่ว่างด้านบนและด้านล่างเพื่อให้ได้ตำแหน่งที่ถูกต้อง ไม่มีใครรู้เกี่ยวกับเครื่องมืออื่น ๆ ที่สามารถใช้สร้างไทล์สำหรับ google maps ได้บ้าง?

16 google-maps  tiles 

ฉันกำลังโยกย้ายรหัสจาก sp แพ็คเกจไปยังแพ็คเกจ sf ที่ใหม่กว่า รหัสก่อนหน้าของฉันฉันมีรูปหลายเหลี่ยม SpatialDataFrame (censimentoMap) และ SpatialPointDataFrame (indirizzi.sp) และฉันได้รับรูปหลายเหลี่ยมเซลล์ id ("Cell110") สำหรับแต่ละจุดที่วางไว้ภายในด้วยคำแนะนำด้านล่าง: points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) ที่จริงฉันสร้างวัตถุ sf สองอัน: shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, crs=crs_string) และ indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string) และฉันกำลังมองหาเอสเอฟเทียบเท่าคำสั่งข้างต้น ... ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else …

16 r  overlay  sp  sf 

เป็นไปได้ไหมที่จะแปลงเลเยอร์ฟีเจอร์ใน ArcGIS เป็น GeoPackage โดยตรง ใน ArcTISbox ของเดสก์ท็อป 10.5 ArcTIS ภายใต้ "เครื่องมือการแปลง / เพื่อ Geopackage" ตัวเลือกเดียวที่ปรากฏสำหรับฉันคือ "เพิ่มแรสเตอร์ลงใน GeoPackage" ไม่มี "เพิ่มคุณสมบัติ" หรืออะไรทำนองนั้น ในโพสต์บล็อกนี้เมื่อเกือบ 3 ปีที่แล้ว ESRI อ้างว่า "ที่ 10.2.1 หรือกับเดสก์ท็อป ArcGIS 10.2.2 คุณสามารถสร้าง GeoPackage ที่ว่างเปล่าและเติม GeoPackage โดยการคัดลอกข้อมูลคุณสมบัติลงไป" อย่างไรก็ตามฉันไม่สามารถดูได้ว่า ใน ArcCatalog ไม่มีตัวเลือกในการสร้างไฟล์ GeoPackage

16 arcgis-desktop  convert  arccatalog  geopackage  arcgis-10.5 

ฉันกำลังล้างข้อมูลรูปหลายเหลี่ยม (ขอบเขต) และฉันสงสัยว่าเป็นไปได้หรือไม่ที่จะเพิ่มโหนด / จุดยอดให้กับชุดรูปหลายเหลี่ยมที่เลือกหรือแม้แต่รูปหลายเหลี่ยมเดียว อย่างไรก็ตามเครื่องมือแก้ไขโทโพโลยี QGIS นั้นยอดเยี่ยมมาก - เรียบง่าย แต่แม่นยำมากและไม่ต้องทำงานมากนัก QGIS ให้ความสามารถในการลบโหนดที่เลือก - น่าทึ่งที่ทอพอโลยีถูกเก็บรักษาไว้และปรับปรุงให้ดีขึ้นเมื่อทำการล้างโหนดส่วนเกิน สิ่งที่ยอดเยี่ยม!

16 qgis  editing  vertices 

ฉันมีตารางย้อนหลังและตารางคะแนนในโพสต์จิส ฉันรู้ว่าใกล้ที่สุดถึงจุดใดก็ตาม สิ่งที่ฉันต้องรู้คือจุดที่ "ด้าน" ของเส้นนั้นคืออะไร ฉันเดาว่าฉันต้องทำอย่างนั้นโดยการสร้างเส้นตั้งฉากจากจุดที่กำหนดไปยังบรรทัด (จุดที่ใกล้ที่สุดในบรรทัด) แล้วเปรียบเทียบพิกัด แต่ฉันไม่รู้วิธีการทำเช่นนั้นและถ้ามันเป็นวิธีที่เหมาะสม เนื่องจากการเปลี่ยนแปลงของเส้นเป็นทิศทาง ฉันทำรูปภาพเพื่ออธิบายงานของฉัน เส้นนั้นเป็นสีดำทิศทางของมันจะแสดงด้วยลูกศรสีเขียว ฉันต้องเพิ่มคอลัมน์ "ด้าน" ลงในตารางจุดดังนั้นจุดสีแดงควรมีค่า "ถูกต้อง" และจุดสีน้ำเงินควรมีค่า "ซ้าย" ใครสามารถให้ตัวอย่างรหัส SQL ของการคำนวณค่า "ด้าน" จุด?

16 postgis  point  linestring 

เราจะได้รับการสรุปการนับจำนวนของเซลล์ภายในแรสเตอร์ใน QGIS อย่างไร ฉันรู้ว่า QGIS กำลังทำสิ่งนี้เพื่อจัดทำสถิติสรุปและฮิสโตแกรมในคุณสมบัติของเลเยอร์ ฉันไม่ต้องการสถิติสรุปหรือภาพฮิสโตแกรม แต่ควรนับค่าแต่ละค่า ดูเหมือนว่ามีอยู่ช่วงหนึ่งที่ความพยายามที่จะเพิ่มคุณสมบัติพิเศษสำหรับสิ่งนี้ (ฟีเจอร์ # 4321, http://hub.qgis.org/issues/4321 ) ฉันขาดอะไรบางอย่างที่ชัดเจนจริง ๆ ?

16 qgis  raster 

งูใหญ่แพคเกจ geopyมีสองวัดระยะทางเทคนิค: เกรทเซอร์เคิและVincenty ของสูตร >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd >>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg >>> vincenty(p1, p2).meters 429.16765838976664 >>> great_circle(p3, p4).meters 428.4088367903001 อะไรคือความแตกต่าง? ต้องการการวัดระยะทางแบบใด

16 python  distance  great-circle  geopy  vincenty-formulae 

มีวิธีการใดบ้างใน ArcGIS 10.2 เพื่อเลือกคะแนนย่อยแบบสุ่ม ตัวอย่างเช่นในภาพหน้าจอที่แนบมาฉันสนใจที่จะเก็บคะแนนที่เลือกไว้ 20% และลบที่เหลือ

16 arcgis-desktop  arcpy  editing  select  random 

ฉันกำลังมองหาคำแนะนำในการเขียนการทดสอบอัตโนมัติสำหรับปลั๊กอิน QGIS ที่เขียนใน Python ฉันได้เขียนการทดสอบสำหรับสคริปต์ Python ในอดีตโดยใช้ PyUnit ( unittestโมดูล) แต่ไม่เคยทำเช่นนั้นสำหรับแอปพลิเคชันที่มี GUI ฉันพบหน้าอธิบายวิธีใช้ PyQt4.QTest เพื่อทำการทดสอบหน่วยในวิดเจ็ต Qt ( http://www.voom.net/pyqt-qtest-example ) แต่ฉันพยายามดิ้นรนเพื่อดูว่าฉันสามารถใช้สิ่งนี้ได้อย่างไร ด้วยวิดเจ็ตที่ได้รับการออกแบบให้เรียกใช้จากภายใน QGIS ส่วนที่เกี่ยวกับ "การทดสอบ" ในเอกสารประกอบ PyQGISนั้นไม่มีอยู่ สิ่งที่ฉันมีอยู่คือ: เก็บการประมวลผลข้อมูลจริงในโมดูลหรือฟังก์ชั่นแยกและเขียนการทดสอบหน่วยสำหรับสิ่งเหล่านั้น ทำการทดสอบพื้นฐานของ UI โดยใช้ QTest โปรดอธิษฐานให้ทุกอย่างอยู่ด้วยกันเมื่อใช้ปลั๊กอินจากภายใน QGIS มีวิธีที่ดีกว่า?

16 pyqgis  qgis-plugins  unit-testing 

มันถูกต้องหรือไม่ที่จะมีองค์ประกอบคุณสมบัติที่มี featureCollection-element เป็นผู้ปกครอง? นี่คือตาม geojson.org ที่ถูกต้อง: { "type": "FeatureCollection", "features": [ { "type": "Feature", "geometry": {"type": "Point", "coordinates": [102.0, 0.5]}, "properties": {"prop0": "value0"} } ] } แต่ฉันไม่พบว่ามันถูกต้องหรือไม่ถูกต้องหากมีสิ่งนี้: { "type": "FeatureCollection", "properties" : { "description" : "This is the geometry for..." } "features": [ { "type": "Feature", "geometry": {"type": "Point", "coordinates": …

16 geojson