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

ฉันมีไฟล์ geojson ที่แสดงผลไม่ถูกต้องและฉันเชื่อว่าเกี่ยวข้องกับไฟล์ที่มีมัลติโพลิกอนที่ไม่ปฏิบัติตามกฎทางขวา เครื่องมือใดที่สามารถช่วยฉันแก้ไขปัญหานี้ได้ ฉันได้ลองนำเข้าและส่งออกกับ QGIS และแปลงไปมาระหว่างรูปแบบด้วยorg2orgแต่นั่นไม่ได้พาฉันไปทุกที่

9 geojson 

ฉันพยายามที่จะดึงข้อมูลข้อความจากเขตข้อมูลสตริงที่มีข้อความคั่นด้วยจุดและวางไว้ในสาขาใหม่โดยใช้เครื่องคิดเลขฟิลด์ ฉันใช้ฟังก์ชัน python นี้ (นำมาจากวิธีการแยกข้อความก่อน / ใน QGIS? ): from qgis.core import * from qgis.gui import * @qgsfunction(args='auto', group='Custom') def func(value1, feature, parent): return value1.split('.')[0] ส่งคืนข้อความทั้งหมดก่อนจุดแรก ตอนนี้ฉันสงสัยว่าจะเขียนฟังก์ชันเพื่อวางแถวข้อความที่คั่นด้วยแต่ละตัวในฟิลด์แยก ก่อน: หลังจาก:

9 qgis  pyqgis  field-calculator  attribute-table  qgis-custom-function 

ฉันจะยุติสคริปต์ใด ๆ ในหน้าต่างคอนโซล Python ของ QGIS ได้อย่างไร ctrl + c ไม่ทำงาน

9 pyqgis  qgis-python-console 

ฉันกำลังพยายามสร้างดัชนีเชิงพื้นที่สำหรับรูปร่างบางไฟล์โดยใช้ Qgis API ฉันได้ทำตามขั้นตอนที่อธิบายไว้ในบล็อกของ Nathan Woodrow ( https://nathanw.net/2013/01/04/using-a-qgis-spatial-index-to-speed-up-your-code/ ) และใช้งานได้ดี : layer = QgsVectorLayer(path, name, 'ogr') idx = QgsSpatialIndex() all_features = layer.getFeatures() map(idx.insertFeature, all_features) ปัญหาของฉันคือฉันไม่มีไฟล์ที่เป็นผลลัพธ์ (.qix? .sbn? .sbx?) ฉันจะบันทึกดัชนีเชิงพื้นที่นี้เพื่อให้ผู้ใช้ในอนาคตของรูปร่างไฟล์ได้ประโยชน์จากมันได้อย่างไร

9 qgis  pyqgis  spatial-index 

ฉันต้องการตรวจสอบไฟล์เวคเตอร์ถ้ามีคอลัมน์ "แหล่งที่มา" ถ้าใช่ไม่มีอะไรให้ทำ หากไม่มีคอลัมน์ "แหล่งที่มา" ควรสร้างและเติมด้วย '0' เป็นไปได้ไหมที่จะทำผ่านตัวสร้างโมเดล

9 qgis  field-calculator  qgis-processing  qgis-modeler 

ฉันกำลังพยายามแยกชั้นเวกเตอร์ที่ซ้อนกันสองชั้นและค้นหาความแตกต่างใน QGIS 2.18 ฉันพยายามค้นหาความแตกต่างกับเลเยอร์ที่มีคะแนนและบัฟเฟอร์ด้วย ทั้งคู่ให้ข้อผิดพลาดต่อไปนี้: global name 'QgsWKBTypes' is not defined See log for more details ฉันค้นหาสิ่งนี้แล้ว แต่ไม่สามารถหาทางแก้ไขได้

9 qgis  error  qgis-processing 

ฉันกำลังพยายามโหลดโครงการ. qgs ที่มีอยู่และรับรายการเลเยอร์ภายใน ฉันดูเหมือนจะสามารถเปิดโครงการได้ แต่ฉันไม่สามารถรับรายการเลเยอร์ได้ ฉันสับสนเล็กน้อยเกี่ยวกับไวยากรณ์ PyQGIS พื้นฐานและตำราอาหารไม่ครอบคลุมตัวอย่างเฉพาะนี้ from qgis.core import * from PyQt4.QtCore import QFileInfo QgsApplication.setPrefixPath(r"C:\OSGeo4W\apps\qgis-ltr", True) qgs = QgsApplication([], False) qgs.initQgis() # Get the project instance project = QgsProject.instance() # Open the project project.read(QFileInfo( r'C:\path\to\project\project.qgs')) print project.fileName() # Get the layers in the project layers = QgsMapLayerRegistry.instance().mapLayers() print …

9 pyqgis 

ฉันพยายามเปลี่ยนเฉพาะการเข้ารหัสแหล่งข้อมูล shapefile ผ่านคุณสมบัติและวิธีการเข้ารหัส shapefiles จาก LATIN1 เป็น UTF-8 แต่มันไม่เปลี่ยน หากไม่มีการบันทึกโครงการ QGIS ฉันไม่สามารถเปลี่ยนการเข้ารหัสแหล่งข้อมูล shapefile โดยตรงได้หรือไม่ ฉันยังลองใช้รหัสนี้มันแสดงการเปลี่ยนการเข้ารหัส แต่เมื่อฉันเปิด shapefile อีกครั้งมันจะแสดงการเข้ารหัสระบบเดิม: for layer in QgsMapLayerRegistry.instance().mapLayers().values(): layer.setProviderEncoding(u'UTF-8') layer.dataProvider().setEncoding(u'UTF-8') print layer.name(), layer.dataProvider().encoding() (ฉันไม่ได้บันทึกโครงการ QGIS เพียงแค่เข้าถึงและบันทึก shapefile โดยตรงเมื่อฉันบันทึกโครงการและทำงานนี้มันจะทำงานอย่างถูกต้อง )

9 qgis  pyqgis  encoding 

ฉันมีไฟล์ ASCII พร้อมละติจูด, ลองจิจูด, และ data_val ในรูปแบบต่อไปนี้ 35-13.643782N, 080-57.190157W, 118.6 ... ฉันมีไฟล์ภาพ GeoTiff และฉันสามารถดูได้อย่างง่ายดาย ฉันต้องการวาง "พิน" (อาจเป็นจุด / ธง / ดาวหรืออะไรก็ตามที่ง่ายที่สุด) ในภาพที่ตำแหน่งละติจูด / ลองจิจูดเฉพาะที่พบในไฟล์ ASCII นี่คือสิ่งที่ฉันสามารถทำได้: ภาพต้นฉบับของฉันดูเหมือนว่า: Driver: GTiff/GeoTIFF Files: /tmp/Charlotte SEC 100.tif Size is 16867, 12358 Coordinate System is: PROJCS["Lambert Conformal Conic", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010042, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], …

9 python  gdal  latitude-longitude  geotiff-tiff  ascii 

ฉันได้เพิ่มเลเยอร์ข้อความที่คั่นด้วย QGIS มาก่อนด้วยคอลัมน์ที่ประกอบด้วยละติจูดและลองจิจูดแยกกัน แต่ตอนนี้ฉันมีไฟล์ Excel เช่น: station no parcel no parcel coordinates 63 215 X1,Y1; X2,Y2; X3,Y3; X4,Y4 ดังนั้นฉันจึงมีพิกัดพัสดุทั้งหมดในหนึ่งคอลัมน์สำหรับทุกพัสดุ ฉันต้องการเพิ่มข้อมูลนี้ด้วยพัสดุลงใน Shapefile ก่อนหน้าของฉันโดยไม่มีสถานีข้อมูล แต่ฉันไม่สามารถหาวิธีนำเข้าพิกัดนี้ไปยัง QGIS ได้

9 qgis  coordinates  import  excel 

ฉันกำลังแปลงเวกเตอร์เป็นแรสเตอร์ใน R อย่างไรก็ตามกระบวนการนี้ยาวเกินไป มีความเป็นไปได้ไหมที่จะนำสคริปต์ไปใช้ในการประมวลผลแบบมัลติเธรดหรือ GPU เพื่อให้เร็วขึ้น? สคริปต์ของฉันเป็นเวกเตอร์ที่แรสเตอร์ r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') r.raster คลาส: RasterLayer ขนาด: 9636, 11476, 110582736 (nrow, ncol, ncell) ความละเอียด: 10, 10 (x, y) ขอบเขต: 505755, 620515, 8555432, …

9 raster  r  rasterization  rgdal  parallel 

ฉันได้ลองเพิ่มชุดข้อมูลใน QGIS โดยใช้ 'qgis.utils.iface.addRasterLayer ()' แต่มันจะพ่นข้อผิดพลาดนี้ เลเยอร์ไม่ถูกต้อง: ผู้ให้บริการ GDAL ไม่สามารถเปิดชุดข้อมูล GDAL http://precisionedgeag.com:8080/arcgis/rest/services/ImageServices/RapidEye_R_NDVI?f=json&pretty=true : `C: \ DOCUME ~ 1 \ 15937.DSM \ LOCALS ~ 1 \ Temp \ file.dat 'ไม่รู้จักเป็นรูปแบบไฟล์ที่รองรับ ผู้ให้บริการเลเยอร์แรสเตอร์ไม่ถูกต้อง (ผู้ให้บริการ: gdal, URI: http://precisionedgeag.com:8080/arcgis/rest/services/ImageServices/RapidEye_R_NDVI?f=json&pretty=true จะเพิ่มเซิร์ฟเวอร์ ArcGIS นี้ (ชุดข้อมูลอิมเมจเซอร์วิส) ใน QGIS ได้อย่างไร?

9 qgis  arcgis-server  image-service 

ในขณะที่พยายามที่จะทำให้รูปหลายเหลี่ยมชั้น Polyline ใน PyQGIS ใช้ processing.runandload("qgis:linestopolygons",explode_path,polygon_path) คำสั่งในคอนโซลหลามมันไม่ได้สร้างรูปหลายเหลี่ยมที่สอดคล้องกับเส้นขอบ แต่มันจะสร้างรูปหลายเหลี่ยมเศษไม้ตามขอบโค้ง ทำไมรูปหลายเหลี่ยมไม่ถูกสร้างตามแนวเส้น? ในรูปหลายเหลี่ยมเส้นสีดำควรเป็นขอบของรูปหลายเหลี่ยมที่สร้างขึ้น มีเพียงไม่กี่ตัวเท่านั้นที่สร้างขึ้นตามขอบโค้ง (รูปหลายเหลี่ยมสีน้ำเงิน)

9 pyqgis  qgis-processing  polygonize 

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

9 qgis  line 

จากแบบจำลองฉันต้องการบันทึกผลลัพธ์การประมวลผลในไฟล์ ในแบบจำลองของฉันฉันมีExtract by locationอัลกอริทึม (ที่มีอินพุตแบบเวกเตอร์สองตัว) ให้ผลลัพธ์ที่สมเหตุสมผล แต่เป็นเลเยอร์ชั่วคราวซึ่งฉันต้องทำSave as...เอง ดีที่สุดฉันต้องการส่งเส้นทางไปยังแบบจำลองของฉันซึ่งผลลัพธ์ของฉันจะถูกบันทึกเป็น shapefiles (ในขั้นตอนแรก) นี่คือรูปแบบ โปรดทราบว่าเอาต์พุตเวกเตอร์ (ไฟล์) ไม่สามารถเชื่อมต่อกับผลลัพธ์การประมวลผลใด ๆ : ดังที่เห็นในภาพหน้าจอ a) เราสามารถส่งชื่อสำหรับผลลัพธ์สุดท้าย แต่ไม่ใช่ไฟล์แม้แต่ในSave selected featuresอัลกอริทึม b) ก) ข) ฉันมีความคิดที่คลุมเครือว่าจะบรรลุผลอย่างไรด้วยสคริปต์ไพ ธ อน แต่สิ่งนี้สามารถทำได้ด้วยการรวมกันของเครื่องมือที่มีอยู่หรือไม่

9 qgis  pyqgis  qgis-processing  qgis-modeler  save