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

ฉันกำลังพยายามสร้างรูปหลายเหลี่ยม voronoi ใน QGIS ที่จะพิจารณา "หลุม" ในโดเมนทั่วไป ตัวอย่างจะเป็น: ฉันสร้าง Voronois ในภาพนี้โดยใช้ QGIS ผ่านคำสั่ง GRASS จากนั้นใช้เครื่องมือ "ความแตกต่าง" เพื่อสร้างรู มีการใช้รูปร่างไฟล์รูปหลายเหลี่ยมแบบแยกซึ่งประกอบด้วยส่วนต่าง ๆ ของหลุมเป็นชั้น "ส่วนต่าง" แอปพลิเคชันตัวอย่างจะสร้างรูปหลายเหลี่ยมรอบจุดสุ่มตัวอย่างที่รวบรวมระหว่างโครงสร้างที่ควรแยกออกจากการวิเคราะห์ ปัญหาสองประการเกิดขึ้นที่นี่: ฟังก์ชั่น "ความแตกต่าง" ไม่ทำงานอย่างถูกต้อง 100% โดยมีขอบเขตรูปหลายเหลี่ยมบางส่วนขยายเข้าไปใน "หลุม" สิ่งนี้สามารถแก้ไขได้โดยการค้นหาแถวในตารางแอตทริบิวต์ซึ่งไม่มีหมายเลขโพลีกอน (หรือ ID ของ "0") "การเจาะรู" หลังจากนั้นความจริงประเภทนี้อาจส่งผลให้เกิดรูปหลายเหลี่ยมที่ไม่ต่อเนื่องดังที่แสดงโดยลูกศรสีแดงในภาพ คำถามของฉันคือ: มีเครื่องมือหรือปลั๊กอิน Voronoi ที่สามารถพิจารณาการปรากฏตัวของ "หลุม" ในใจกลางของโดเมนเป็นกระบวนการขั้นตอนเดียวและยังกำจัดการสร้างรูปหลายเหลี่ยมไม่ต่อเนื่อง? ฉันจินตนาการว่าเครื่องมือดังกล่าวจะขยายขอบเขตของรูปหลายเหลี่ยมไปยังจุดตัดที่ใกล้ที่สุดพร้อมกับขอบเขตอื่นยกเว้นว่าจะมีการต่อยกับเขต "หลุม" เป็นครั้งแรก

12 qgis  masking  voronoi-thiessen 

ดังนั้นฉันจึงเขียนสคริปต์ที่ทำสิ่งที่ฉันต้องการซ้ำแล้วซ้ำอีกโดยใช้ "ในขณะที่ True:" ในช่วงเวลาหนึ่ง (ทุก ๆ 5 วินาทีโดยใช้ time.sleep (5)) จนถึงตอนนี้ดีมาก แต่เมื่อฉันต้องการหยุดมันฉันก็ทำไม่ได้ ฉันได้ลอง Control + C, Control + Break, Escape แล้วก็ข้ามคีย์บอร์ดของฉันไป วิธีเดียวที่จะหยุดมันได้คือปิด QGIS ความคิดใด ๆ นอกจากนี้เมื่อสคริปต์เปิดใช้งาน time.sleep (5) ชนิด QGIS ของ lags และค้างเป็นเวลา 5 วินาทีและฉันไม่สามารถยกเลเยอร์ได้ แต่ฉันถือว่านี่เป็นเรื่องปกติ นี่คือสคริปต์ของฉัน: from PyQt4.QtGui import * from PyQt4.QtCore import * from qgis.core import * from …

12 pyqgis  loop 

มีแพ็คเกจ Python ที่ให้การทำงานของ Straight Skeleton algorithm หรือไม่? ฉันทราบว่าCGALโครงการโอเพ่นซอร์ส (C ++) ได้รับการนำไปใช้แต่ดูเหมือนว่าการผูกมัดของ cgal จะไม่รวมแพ็คเกจ CGALนี้ ไม่ว่าในกรณีใดฉันต้องการใช้ Python อย่างแท้จริงซึ่งฉันสามารถแก้ไข / ขยายได้ตามความต้องการ ในขณะที่การใช้งานที่สามารถจัดการรูปหลายเหลี่ยมกับหลุมจะดีกว่าก็ไม่จำเป็นอย่างเคร่งครัด

12 python  polygon  geometry  skeleton  computational-geometry 

นี่คือสิ่งที่ฉันกำลังทำ: สร้างเลเยอร์โดยอัปโหลดรูปร่างไฟล์ลงใน QGIS เลือกคุณสมบัติของไฟล์นั้น ส่งออกแอตทริบิวต์ที่เลือกด้วยช่องทำเครื่องหมาย "บันทึกคุณสมบัติที่เลือกเท่านั้น" เมื่อฉันบันทึกไฟล์ ฉันต้องการส่งออกชื่อที่สอดคล้องกับรูปหลายเหลี่ยมแต่ละอัน น่าเสียดายที่เมื่อฉันบันทึกเป็นไฟล์ KML และเปิดใน Google แผนที่ของฉันฉันสามารถเห็นขอบเขตรูปหลายเหลี่ยมเท่านั้นดังนั้นฉันต้องตั้งชื่อพวกเขา ข้อมูลอยู่ในคอลัมน์ในตารางแอตทริบิวต์

12 qgis  shapefile  google-maps  kml  attribute-table 

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

12 qgis  polygon  vector  overlapping-features 

ฉันมีหลายเลเยอร์ฐานและจำเป็นต้องเปลี่ยนเป็นโปรแกรมขึ้นอยู่กับการกระทำและสถานะต่าง ๆ ในแอปพลิเคชันของฉัน ไม่สามารถหาได้อย่างไรเพื่อชีวิตของฉัน

12 leaflet  mapbox 

เป็นไปได้ไหมที่จะเปิดฐานข้อมูลส่วนบุคคลของ Esri (* .mdb; ArcGIS 10.2) ใน QGIS 2.6? รุ่นก่อนหน้านี้ถูกต้อง

12 qgis  esri-geodatabase  personal-geodatabase 

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

12 qgis  polygon  buffer 

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

12 qgis  gps  line 

ปิด คำถามนี้เป็นคำถามความคิดเห็นตาม ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้สามารถตอบข้อเท็จจริงและการอ้างอิงได้โดยแก้ไขโพสต์นี้ ปิดให้บริการใน4 ปีที่แล้ว ฉันกำลังมองหาการติดตั้ง GNSS / GPS ราคาประหยัด (<US $ 500) ด้วยความแม่นยำในระดับเซนติเมตรโดยไม่ต้องแฮ็ค HW มาก ฉันไม่สามารถผลิต PCB ของฉันหรือบัดกรีใด ๆ (แม้ว่าฉันจะทำอย่างนั้นถ้าไม่มีวิธีอื่น) ดังนั้นการตั้งค่าที่ง่ายต่อการประกอบจะได้รับการต้อนรับ ฉันรู้เรื่องของพิกซี่ 900 เหรียญแต่มันก็แพงเกินไปสำหรับฉัน ดูเหมือนว่าความแม่นยำซม. ควรจะเป็นไปได้น้อยกว่ามากเช่นใช้เซ็นเซอร์ GPS ดิบ 50 USD พร้อมเสาอากาศและพีซีทั่วไปที่มีซอฟต์แวร์ RTKLIB ฉันไม่แน่ใจว่าจะเป็นการดีกว่าถ้าใช้การตั้งค่าเซ็นเซอร์ GPS สองตัวสำหรับ RTK (หนึ่งสถานีฐานและอีกหนึ่งสำหรับ rover) หรือว่าฉันจะได้รับข้อมูล DGPS ที่ถูกต้องแก้ไขที่อื่น (ภูมิภาคของฉันคือสาธารณรัฐเช็ก - ดูเหมือนว่ามีกริดแห่งชาติที่นี่ เพื่อสตรีมข้อมูลการแก้ไขในราคาที่สมเหตุสมผล) แอปพลิเคชันของฉันจะอยู่ในรถโดยสารดังนั้นฉันจะไม่ถูก จำกัด …

12 gps  accuracy  precision  rtklib 

ฉันรู้ว่าต้องมีวิธีการทำเช่นนี้ แต่ฉันมีปัญหาในการหาวิธีการทำเช่นนั้น ฉันต้องการสร้างรายการของฟีเจอร์คลาสเพื่อทำการดำเนินการอย่างไรก็ตามฉันต้องการให้ลิสต์นั้นมีคลาสของฟีเจอร์ที่มีตัวระบุชื่อไฟล์ที่แตกต่างกัน 2 รายการ รหัสจะมีลักษณะดังนี้: FCS = arcpy.ListFeatureClasses("X_*" OR "*_Y") ฉันรู้ว่า OR ไม่ถูกต้อง แต่มีวิธีการบางอย่างนี้โดยไม่ต้องเขียนโค้ดเพื่อทำงานกับ wildcard ตัวที่สอง ขอบคุณที่ช่วยเกาต์ฉันรู้ว่ามีวิธีการเลือกใน Select By Attributes ซึ่งอนุญาตให้ใช้สัญลักษณ์ตัวแทนหลายรายการโดยใช้คำสั่ง OR

12 arcpy  list  wildcard 

ฉันต้องการอัปเดตคุณลักษณะหนึ่งคุณลักษณะ อย่างไรก็ตามฉันไม่ทราบว่าจะอัปเดตโดยใช้คุณสมบัติของวัตถุ ฉันต้องใช้ผู้ให้บริการข้อมูลเพื่ออัปเดต layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it = dp.getFeatures() for i in range(0, dp.featureCount()): feat = it.next() attrs = { 2 : 30 } layer.dataProvider().changeAttributeValues({ feat.id() : attrs }) ฉันสามารถเปลี่ยนค่าของคุณสมบัติโดยใช้QgsFeatureวัตถุได้หรือไม่? ยิ่งไปกว่านั้นมันเป็นไปได้ที่จะวนซ้ำโดยใช้วัตถุตัววนซ้ำ?

12 qgis  python  pyqgis  features 

ฉันมีปัญหาในการใช้ฟังก์ชัน Polygonize ใน python ตัวอย่างตำราสำหรับการนี้สามารถพบได้ที่นี่ ส่วนที่เกี่ยวข้องของรหัสของฉันคือ: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = driver.CreateDataSource(outShapefile+ ".shp") outLayer = outDatasource.CreateLayer("polygonized", srs=None) gdal.Polygonize( band, None, outLayer, -1, [], callback=None ) outDatasource.Destroy() sourceRaster = None ฉันรู้ว่าวงดนตรีมีข้อมูลที่เกี่ยวข้องนี่เป็นตัวอย่างของbandArray: array([[ 4., 4., 3., …

12 python  gdal 

ฉันพยายามใช้ ogrinfo เพื่อรับรายละเอียดบางอย่างเกี่ยวกับรูปร่างไฟล์ที่ฉันดาวน์โหลด ขณะนี้วิธีเดียวที่ฉันรู้วิธีการทำเช่นนี้คือการโหลดลงใน QGIS และคลิกด้วยตนเองเพื่อค้นหาข้อมูลใด ๆ เช่นเปิดตารางแอตทริบิวต์ ฉันแค่ต้องการให้สามารถดูข้อมูลเมตาใด ๆ ที่ติดแท็กพร้อมกับคุณสมบัติ ถ้าฉันทำ: ogrinfo -al USA_adm0.shp ฉันเห็นได้ตั้งแต่ต้นว่ามีข้อมูลที่มีประโยชน์มากมาย แต่แล้วมันก็บินผ่านไปพร้อมกับข้อมูลคุณสมบัติทั้งหมด มีคนช่วยฉันได้ไหม แก้ไข นี่คือสิ่งที่ฉันได้รับบน mac ของฉันโดยใช้แฟล็ก -ro และ -so ดูเหมือนจะไม่ช่วยอะไรมาก ->ogrinfo -ro -so USA_adm0.shp INFO: Open of `USA_adm0.shp' using driver `ESRI Shapefile' successful. 1: USA_adm0 (Polygon)

12 shapefile  gdal  ogr 

เป็นเรื่องปกติหรือไม่ที่จะใช้รูปทรงเรขาคณิต 'หลายจุด' แทน 'จุด' รู้ว่าฉันใช้มันเพื่อเก็บจุดเดียวหรือไม่? ฉันจะแปลงจากหลายจุดเป็นจุดได้อย่างไร

12 postgis  postgresql  geometry