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

ถามตอบสำหรับนักทำแผนที่นักภูมิศาสตร์และผู้เชี่ยวชาญด้าน GIS

4
พยายามแยกรายการของค่าที่ไม่ซ้ำออกจากเขตข้อมูลโดยใช้ python
ฉันมีจำนวนคอลัมน์ในตารางจำนวนหนึ่งภายใน FGDB ที่ฉันต้องการแยกค่าที่ไม่ซ้ำกันสำหรับแต่ละคอลัมน์ ตัวอย่างเช่น: ค่าอาจเป็น [1,2,2,2,3,4] และฉันพยายามคืน [1,2,3,4] ฉันสามารถทำงานนี้ได้หลายวิธีใน ARCGIS แต่ฉันพยายามขยายตัว ฉันพบชิ้นส่วนของหลามบนเว็บที่ฉันคิดว่าจะทำงานได้ แต่ฉันพยายามที่จะให้มันทำงาน (ฉันได้รับข้อผิดพลาดทางไวยากรณ์ที่ไม่ถูกต้องเนื่องจากฉันได้รับข้อผิดพลาดทางไวยากรณ์ในบรรทัดที่ 3) ข้อผิดพลาดของผู้ใช้ที่ง่ายมาก ตัวอย่างโค้ดด้านล่าง import arcpy def unique_values(r'N:\GISProjects\Landuse\Plant_Biosecurity_Project\ArcGIS_Online.gdb\Holdings_Property_Merge' , 'LU_ALUMMaj'): with arcpy.da.SearchCursor(table, [field]) as cursor: return sorted({row[0] for row in cursor}) นี่คือข้อความแสดงข้อผิดพลาดที่ฉันได้รับจากข้อความประเสริฐ: File "C:\Users\hawkinle\Desktop\STDTAS\Unique_Data.py", line 3 def unique_values(r'N:\GISProjects\Steve_Eastwood_Landuse\Plant_Biosecurity_Project\ArcGIS_Online.gdb\Holdings_Property_Merge' , 'LU_ALUMMaj'): ^ SyntaxError: invalid syntax [Finished in 0.1s …

4
บันทึกเลเยอร์ WMS ไว้ในฮาร์ดไดรฟ์ของคุณผ่าน QGIS หรือไม่
นี่คือคำถามอื่น ฉันกำลังทำงานกับชุดข้อมูล WMS ในตอนนี้และฉันต้องการทราบว่าคุณสามารถบันทึกลงในฮาร์ดไดรฟ์ได้หรือไม่ ฉันรู้ว่าคุณสามารถบันทึกมุมมองปัจจุบันเป็น Tif ที่อ้างอิงทางภูมิศาสตร์ได้ แต่เมื่อพิจารณาถึงความละเอียดที่ฉันต้องการและพื้นที่ที่ครอบคลุมโดยพื้นที่ทำงานของฉันฉันต้องการสิ่งที่เร็วกว่ามากหรือฉันจะทำสิ่งนี้เป็นเวลาหลายสัปดาห์
11 qgis  wms  export  save 

1
แปลงลองจิจูด / ละติจูดเป็นพิกเซลบนแผนที่
ฉันมีแผนที่จากที่นี่ ฉันต้องการที่จะสามารถแปลงคู่ lon / lat โดยพลการเป็นพิกเซลบนแผนที่ (เช่นเดียวกับความสามารถในการแปลงกลับ) แผนที่มาพร้อมกับไฟล์. tfw และข้อมูลการฉายภาพนี่คือ: 0.02222222222222 0.00000000000000 0.00000000000000 -0.02222222222222 -180.00000000000000 90.00000000000000 และฉายข้อมูล: Projection: Plate Carree aka Geographic or "LatLong" Earth ellipsoid: Sphere, radius 6370997 m Datum: WGS84 Extent: 180 West to 180 East, 90 North to 90 South Size: 16,200 height samples wide x 8,100 …

2
ดึงการเชื่อมต่อ PostGIS ที่มีอยู่ใน PyQGIS
ฉันสามารถดึงการเชื่อมต่อที่มีอยู่ไปยังฐานข้อมูล PostGIS ใน PyQGIS ได้หรือไม่ ฉันต้องการให้รายการของการเชื่อมต่อฐานข้อมูลที่มีอยู่และต่อมาก็คือรายการตารางภายใน UI ของปลั๊กอินของฉัน ฉันตรวจสอบตำราอาหาร แต่ไม่สามารถหาทางที่จะไปให้ไกลกว่านี้ได้
11 qgis  postgis  pyqgis 

1
แผนที่แบบโต้ตอบและแบบไดนามิกที่มีข้อมูลลมในแผ่นพับ R?
ฉันมีข้อมูลบางอย่างเช่นนี้: longitude latitude speed direction 1 6.10722222 46.23639 4 360 2 6.95416667 43.54694 4 360 3 7.21472222 43.66556 13 330 4 4.01666667 48.32167 7 290 5 2.30833333 43.21611 14 290 6 2.48305556 44.40806 13 320 7 5.21500000 43.43694 19 330 8 4.92361111 43.52278 32 320 9 5.10805556 43.60306 26 330 …

3
ใช้ ST_Difference เพื่อลบคุณลักษณะที่ทับซ้อนกันหรือไม่
ฉันกำลังพยายามใช้ ST_Difference เพื่อสร้างชุดรูปหลายเหลี่ยม (การประมวลผล trimmedparcelsnew) ที่ไม่มีพื้นที่ใด ๆ ที่ครอบคลุมโดยรูปหลายเหลี่ยมชุดอื่น (test.single_geometry_1) โดยใช้ PostGis 2.1 (และ Postgres SQL 9.3) นี่คือคำถามของฉัน: CREATE TABLE processing.trimmedparcelsnew AS SELECT orig.id, ST_Difference(orig.geom, cont.geom) AS difference FROM test.single_geometry_1 cont, test.multi_geometry_1 orig; แต่รูปหลายเหลี่ยมที่เกิดขึ้นยังไม่ได้ถูกตัดแต่ง แต่ดูเหมือนว่าจะถูกแยกออกจากกันโดยที่พวกมันตัดกับเลเยอร์อื่น ฉันลองใช้ตัวเลือกโดยไม่ใส่ผลลัพธ์ลงในตารางและทุกอย่างอื่นที่ฉันสามารถนึกได้ แต่ดูเหมือนว่าฉันจะไม่สามารถใช้ฟังก์ชันนี้ได้ ฉันแนบรูปภาพผลลัพธ์แล้ว หลังจากความคิดเห็นฉันได้ลองเพิ่มส่วนคำสั่ง WHERE แล้ว ฉันต้องการพัสดุที่ไม่มีทางแยกและบริเวณที่ตัดกันของพัสดุอื่นถูกลบ (layer test.single_geometry แสดงถึงการปนเปื้อนที่ฉันต้องการลบออกจากพัสดุของฉัน) ฉันลองทางแยก แต่แน่นอนว่าฉันต้องการทางแยกที่ไม่ใช่จริงดังนั้นตอนนี้ฉันจึงลองแยกกัน ฉันได้ลองเพิ่มต้นกำเนิดลงในตารางของฉันแล้ว แต่เอกสารสำหรับ ST_Difference ( …

3
กำลังแปลงละติจูด, ลองจิจูด (epsg: 4326) เป็น EPSG: 3857?
ฉันต้องการศึกษาคณิตศาสตร์ที่เกี่ยวข้องเพื่อดูว่าละติจูดและลองจิจูดที่กำหนดเป็นองศาทศนิยมสามารถแปลงเป็น EPSG: 3857 ได้อย่างไร ทุกคนสามารถชี้ไปที่การอ้างอิงที่ดีหรือเป็นไปได้รหัส Ppython / C / C ++ / Java ห้องสมุดที่เป็นไปได้สำหรับการตรวจสอบคณิตศาสตร์?

3
ฉันจะดูข้อมูล GTFS ใน QGIS ได้อย่างไร?
ฉันต้องการใช้ QGIS เพื่อสร้างแผนที่ของฉัน ฟีดข้อมูลเดียวคือ gtfs ที่เกี่ยวข้องกับเส้นทางรถบัส ฉันสามารถดูได้บน Google Earth (แม้ว่าจะล้าหลัง) แต่เมื่อฉันพยายามนำเข้า. kml ลงใน QGIS มันเพิ่งหมด RAM (เป็นเวลา 30 นาที) และทำให้แผนที่ของฉันเสียหาย หากฉันนำเข้า 10 ในเลเยอร์ 18xx เส้นทางจะซ้ำกัน 10 ครั้ง (Google Earth สามารถแสดงเส้นทางทั้งหมดแม้ว่าจะล่าช้า) ขั้นตอนของฉัน ฉันกำลังลองวิธีการใหม่ดังนั้น ... สิ่งที่ฉันทำคือดาวน์โหลด. txt บนเว็บไซต์ TfNSW (ข้อมูล Google และการเข้าสู่ระบบ tfnsw) จากนั้นใช้สคริปต์. py ชื่อTransitfeedโดย bdferris เพื่อแปลง. txt เป็น. kml ในที่สุดนำเข้า. …

2
การระบุรูปหลายเหลี่ยมที่ทับซ้อนกันในชั้นเดียวโดยใช้ QGIS?
ฉันมีรูปร่างของคะแนนบัฟเฟอร์ใน QGIS ฉันต้องการที่จะแสดงรูปหลายเหลี่ยมบัฟเฟอร์ทั้งหมดในชั้นนี้ซึ่งทับซ้อนกัน ฉันได้ทดลองใช้เครื่องมือตัดกัน แต่ดูเหมือนว่าจะใช้งานได้ถ้าฉันดูที่จุดตัดของ 2 ชั้นแยกกันมิฉะนั้นจะบอกว่ารูปหลายเหลี่ยมทั้งหมดตัดกัน ใครบ้างมีโซลูชัน QGIS

5
ข้อมูลดาวเทียมเรดาร์ (SAR) ฟรี [ปิด]
ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน3 ปีที่ผ่านมา ฉันกำลังมองหาแหล่งที่มาของภาพถ่ายดาวเทียมเรดาร์ (SAR) ฟรี (ไม่ใช่เพื่อการค้า) ไม่มีใครรู้แหล่งที่มา?

3
ปัญหากับเรขาคณิต SRID ใน PostGIS
ฉันเลิกเรียน (Poste) ใน PostGIS id_borne serial NOT NULL, num_borne character varying, shape_borne geometry ฉันเปลี่ยน SRID ของตารางของฉันสำเร็จแล้วโดยใช้ select UpdateGeometrySRID('public', 'borne', 'shape_borne', 26191) ; SRID เก่าคือ 0 ปัญหาคือเมื่อฉันพยายามเพิ่มองค์ประกอบใหม่ในตารางของborneฉันฉันมีข้อผิดพลาดข้อความนี้: geometry srid (0) does not match column srid (26191) นั่นคือปัญหาของฉัน! SRID ได้รับการปรับปรุง แต่เมื่อฉันพยายามเพิ่มจุดใหม่ฉันได้รับข้อความนี้

3
เป็นไปได้ไหมที่จะหมุน Canvas ใน QGIS
ผ้าใบหมุนได้หรือไม่เพื่อให้ทิศเหนือชี้ไปในทิศทางอื่นนอกเหนือจากที่ฉายไว้? คำถามนี้ตอบด้วยความเคารพนักแต่งเพลง แต่ไม่ใช่ผ้าใบ มีรายการอายุ 5 ปีที่เกี่ยวข้องในฐานข้อมูลการติดตามบั๊กที่ปิดไม่จำเป็น นั่นคือสถานะปัจจุบันหรือมีวิธีหมุนผืนผ้าใบในเวอร์ชันที่ใหม่กว่าหรือไม่? การสร้างการฉายภาพใหม่ไม่ใช่คำตอบ
11 qgis  qgis-2 

1
วิธีวาง North Arrow และ Scale บนแผนที่ใน COMPOSER ใน QGIS 2.2.0
ฉันยังใหม่กับ QGIS และคำถามที่ถามและตอบที่นี่เกี่ยวกับ North Arrow ดูเหมือนจะไม่ช่วยฉันแก้ปัญหาด้านล่าง ในบานหน้าต่างการแมปหลักของ QGIS 2.2.0 Valmiera ฉันสามารถเพิ่มมาตราส่วนและลูกศรทิศเหนือ แต่ทันทีที่ฉันไปที่นักแต่งเพลงเพื่อพิมพ์แผนที่สิ่งทั้งสองนั้นหายไปจากแผนที่ที่ฉันเปิดในผู้แต่ง นักแต่งเพลงมีบางสิ่งที่ทำให้ฉันสามารถเพิ่มสเกลสิ่งหนึ่งได้ แต่มันแสดงสเกลจริงบนกระดาษที่ฉันจะพิมพ์และไม่ใช่สเกลจริงของพื้นผิวโลกของแผนที่ที่แสดง ฉันได้ลองรีเฟรชแผนที่แล้ว แต่มันไม่ได้ผล บรรทัดล่างคือ: ฉันจะได้รับขนาด (ขออภัย!) (ขยายและหดตัวในค่าเมื่อฉันซูมเข้าและออก) ในนักแต่งเพลงและฉันจะได้ลูกศรทิศเหนือที่นั่นได้อย่างไร

7
Geocode เป็น Google SpreadSheet ได้อย่างไร
ฉันมีสเปรดชีตใน Google เอกสารประมาณ 20-30 แถวโดยมี 6 คอลัมน์ดังต่อไปนี้: ID, Name, Address, Cell_No, Latitude, Longitude ที่สองคอลัมน์สุดท้ายว่างเปล่า อะไรคือวิธีที่ง่ายสำหรับผู้ใช้มือใหม่ที่ไม่ใช่โปรแกรมเมอร์เพื่อกำหนดที่อยู่ทางภูมิศาสตร์และเติมข้อมูลในช่องละติจูดและลองจิจูดด้วยผลลัพธ์ของการเข้ารหัสทางภูมิศาสตร์ ฉันแค่ต้องเติมฟิลด์ละติจูดและลองจิจูด ฉันไม่จำเป็นต้องเห็นข้อมูลใด ๆ บนแผนที่

3
วิธีการโหลดเอาท์พุทหน่วยความจำจากการประมวลผล QGIS?
อัลกอริทึมการประมวลผลจำนวนมากมีตัวเลือกในการบันทึกผลลัพธ์เป็นหน่วยความจำ หากวิ่งจากกล่องเครื่องมือสิ่งนี้จะใช้ได้เพราะคุณสามารถทำเครื่องหมายในช่อง "เปิดไฟล์เอาต์พุตหลังจากเรียกใช้อัลกอริทึม" มองไปที่ซอร์สโค้ดสำหรับ QGIS Postprocessing.handleAlgorithmResultsรหัสสำหรับการโหลดชั้นหน่วยความจำที่ดูเหมือนว่าจะมีการกำหนดไว้ในการทำงาน ฟังก์ชันใช้การเข้าถึงalg.outputsรายการและทำสิ่งต่อไปนี้: for out in alg.outputs: progress.setPercentage(100 * i / float(len(alg.outputs))) if out.hidden or not out.open: continue if isinstance(out, (OutputRaster, OutputVector, OutputTable)): try: if out.value.startswith('memory:'): layer = out.memoryLayer # NOTE!! QgsMapLayerRegistry.instance().addMapLayers([layer]) # NOTE!! else: # ... เมื่อคุณเรียกใช้อัลกอริทึมการประมวลผลจากคอนโซลมีวิธีโหลดเลเยอร์โดยไม่ต้องเข้าถึงวัตถุนี้หรือไม่? ฉันวิ่งได้ processing.runalg("qgis:intersection", layer1, layer2, "memory:") หรือแม้กระทั่ง processing.runalg("qgis:intersection", layer1, …

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.