คำถามติดแท็ก ogr

ฉันมี shapefile ของเขตการปกครองของประเทศและฉันต้องการแยกอาร์เรย์ของละติจูดและลองจิจูดสำหรับแต่ละแผนก .. มันเป็นไปได้ไหมที่จะทำ?

12 python  shapefile  ogr 

ฉันกำลังมองหาฟังก์ชั่นที่คล้ายกับ GRASS r.resample ฉันมีโฟลเดอร์ที่เต็มไปด้วยไฟล์ TIF และต้องการลดความละเอียด ("ลดขนาดไฟล์ลง") เพื่อลดขนาดไฟล์

12 raster  gdal  ogr 

การใช้ Gdal / OGR มีวิธีแปลงเลเยอร์บรรทัด shp, kml หรือ PostGIS เป็นเครือข่ายของโหนดและเซ็กเมนต์ดังนั้นจึงสามารถใช้ในแพ็คเกจเช่น networkx?

12 python  topology  ogr  networkx 

มีวิธีแปลงไฟล์เชิงเส้นเป็นไฟล์ CSV ในรูปแบบต่อไปนี้หรือไม่: Node_1, Node_2, attribute_1, ...., attribute_x 1 2 "hello" 567845.334 ฉันดูที่GDAL / OGR2OGRและPostGISและฉันไม่คิดว่าจะแปลงข้อมูลเป็นรูปแบบที่ฉันต้องการ ฉันไม่รังเกียจที่จะเขียนสคริปต์ให้ทำถ้าฉันถูกชี้นำไปในทิศทางที่ถูกต้อง

12 shapefile  convert  ogr  csv 

ฉันมีไฟล์ในการฉายภาพ British National Grid: Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936", SPHEROID["Airy_1830",6377563.396,299.3249646]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",49], PARAMETER["central_meridian",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["Meter",1]] cat: Integer (9.0) ฉันสามารถใช้ GDAL / OGR เพื่อให้ได้พื้นที่รวมของรูปหลายเหลี่ยมทั้งหมดใน shapefile เป็นเฮกตาร์หรือไม่ ฉันสงสัยว่ามันเป็นไปได้ด้วย-sqlหรือเปล่า: ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp …

11 gdal  ogr  area  ogrinfo 

ฉันมี shapefile มีคุณลักษณะหลายประการเช่นYEAR, และCOUNTY AREAฉันยังมีไฟล์ CSV กับเขตอื่น ๆ ที่ฉันต้องการใน shapefile POPULATIONเช่น ทั้ง shapefile และไฟล์ CSV GISJOINมีฟิลด์ ฉันรู้วิธีเข้าร่วมใน QGIS แต่ฉันจะทำการเข้าร่วมแบบถาวรและเขียนไปยัง Shapefile โดยใช้ogr2ogrหรือหนึ่งในเครื่องมืออื่น ๆ ใน GDAL / OGR ได้อย่างไร

11 shapefile  ogr  csv  attribute-joins 

ฉันมีรูปร่างไฟล์ที่มีรายละเอียดอย่างเหมาะสมพร้อมด้วยคุณสมบัติรูปหลายเหลี่ยม / รูปหลายเหลี่ยม (ไฟล์ประมาณ 500mb) จริงๆแล้วมันเป็นรูปร่างของโลกทั้งใบด้วยคุณสมบัติที่เป็นตัวแทนของแนวชายฝั่ง ฉันต้องการหารข้อมูลนี้โดยใช้กริด เพื่อความชัดเจนฉันไม่ต้องการ 'เรียงลำดับ' ข้อมูล แต่จริง ๆ แล้วตัดรูปหลายเหลี่ยมขึ้นเป็นแบบเรียงต่อกัน ฉันรู้ว่าคำถามนี้ถูกถามมาก่อน แต่วิธีแก้ปัญหาที่พบนั้นใช้ไม่ได้สำหรับฉัน ฉันได้พยายาม: การใช้ QGIS และตัดเนื้อหา Shapefile ของฉันด้วยเวกเตอร์กริด - ผลลัพธ์แย่มาก เจ้าของที่ดินรายใหญ่ส่วนใหญ่หายตัวไปอย่างน่าอัศจรรย์แม้ว่าบางครั้งมันอาจดูเหมือนว่าเป็นที่ดินขนาดเล็กก็ตาม ฉันควรทราบว่าวิธีนี้ใช้ได้ดีกับข้อมูลที่ง่ายกว่า (เช่น. คะแนนน้อยกว่า) การใช้เครื่องมือ Intersection ของ OGR ฉันลองทั้งผ่าน ogr2ogr และแม้กระทั่งโดยการหมุนเครื่องมือ C ++ ของฉันเอง ทั้งคู่มีปัญหาเดียวกันกับ QGIS พวกเขายังไม่แสดงปัญหานี้สำหรับไฟล์แบบง่าย ๆ แต่จะไม่ทำให้ไฟล์ซับซ้อนขึ้น สำหรับการอ้างอิงฉันใช้ shapefile ของออสเตรเลียและนิวซีแลนด์ต่ำกว่า 20mb ในขนาดและทั้ง QGIS และ OGR …

11 qgis  shapefile  ogr  qgis-processing  vector-grid 

ฉันพยายามสร้างรูปหลายเหลี่ยมรูปร่างจากรายการพิกัดโดยใช้เครื่องมือหลามโอเพนซอร์ส ต่อไปนี้คือสิ่งที่ฉันได้เพื่อให้ห่างไกลที่ถูกแฮ็กด้วยกันจากงูหลาม GDAL / OGR ตำราและGIS นี้ตอบ มีคำถามที่คล้ายกันPython: วิธีการสร้างรูปหลายเหลี่ยม Shapefile จากรายการของพิกัด X, Y pyshpแม้ว่าคำถามนี้จะหมายถึงการใช้ อย่างไรก็ตามฉันสนใจที่จะสร้างไฟล์รูปหลายเหลี่ยมโดยใช้เครื่องมือหลาม gdal / ogr เท่านั้น import ogr def create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in coords: ring.AddPoint(coord[0], coord[1]) # Create polygon poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) return poly.ExportToWkt() def write_shapefile(poly, out_shp): """ https://gis.stackexchange.com/a/52708/8104 """ # Now convert …

10 python  gdal  polygon  ogr 

ฉันพยายามเรียนรู้วิธีใช้ ogr ใน python โดยใช้ประเทศและชุดข้อมูลสถานที่ที่มีประชากรจากhttp://www.naturalearthdata.com/downloads/50m-cultural-vectors/. ฉันกำลังพยายามใช้ตัวกรองและบัฟเฟอร์เพื่อค้นหาคะแนน (ne_50m_populated_places.shp) ภายในบัฟเฟอร์ที่ระบุของประเทศที่ตั้งชื่อ (กรองจาก ADMIN คลาสคุณลักษณะใน ne_50m_admin_0_countries.shp) ปัญหาดูเหมือนว่าฉันไม่เข้าใจหน่วยที่จะใช้สำหรับ buffer () ในสคริปท์ฉันใช้ค่า 10 โดยพลการเพื่อทดสอบว่าสคริปต์ทำงานหรือไม่ สคริปต์ทำงาน แต่ส่งคืนสถานที่ที่มีประชากรจากทั่วภูมิภาค Carribean สำหรับประเทศที่มีชื่อ 'Angola' เป็นการดีที่ฉันต้องการที่จะสามารถระบุระยะทางบัฟเฟอร์พูด 500km แต่ไม่สามารถหาวิธีการทำเช่นนี้เป็นความเข้าใจของฉันคือบัฟเฟอร์ () กำลังใช้หน่วยของ countries.shp ที่จะอยู่ในรูปแบบ wgs84 lat / long . คำแนะนำเกี่ยวกับวิธีการเพื่อให้บรรลุนี้จะได้รับการชื่นชม # import modules import ogr, os, sys ## data source os.chdir('C:/data/naturalearth/50m_cultural') # get the …

10 python  ogr  buffer 

ฉันใช้ qgis และหญ้าสำหรับงานของฉัน ฉันได้สร้าง shapefile ซึ่งสูงกว่า 100mbs (ผลการวิเคราะห์ทางสถิติบางอย่าง) หากฉันแปลงเป็นรูปแบบ sqlite (เปิดใช้งานเชิงพื้นที่) ข้อมูลจะลดขนาดลงเหลือ 15 mbs สิ่งที่ยอมรับได้สำหรับการแบ่งปันกับอีเมลหรือดรอปบ็อกซ์ แต่เพื่อนร่วมงานของฉันใช้ arcgis 9.3 - เขาจะสามารถอ่านข้อมูลได้โดยไม่มีปัญหาหรือไม่?

10 arcgis-desktop  open-source-gis  arcgis-9.3  ogr  sqlite 

เป็นไปได้หรือไม่ที่จะจับข้อความแสดงข้อผิดพลาดเอาต์พุต PostgreSQL ใน OGR หลังจากแบบสอบถามล้มเหลว ตัวอย่างเช่น conn.ExecuteSQL('SELECT * FROM non_existing_table;') ใน Postgres สิ่งนี้จะส่งคืนข้อความแสดงข้อผิดพลาด: "Error: relation 'non_existing_table' does not exist." ใน ogr เราได้มาตรฐาน: "Error: current transaction is aborted, commands ignored until end of transaction block" มันจะเป็นโบนัสถ้าเราสามารถจับ / แยกวิเคราะห์ข้อความแสดงข้อผิดพลาดของ Postgres ได้ ฉันสงสัยว่านี่เป็นช็อตเล็กน้อยที่ได้รับจากสิ่งที่ฉันได้อ่านเกี่ยวกับ SWIG / CPL แต่ก็คุ้มค่าที่จะถามว่าใครมีไอเดียไหม เป็นการดีที่จะเป็น Python แต่ฉันจะแสดงความคิดเห็น C / C …

10 python  postgis  gdal  postgresql  ogr 

มือใหม่ที่นี่สับสนจากประสบการณ์ครั้งแรกของฉันใน QGIS ฉันมีเส้นจำนวนหนึ่งที่กำหนดเส้นทางจักรยาน แต่เดิมพวกเขาถูกบันทึกไว้ในไฟล์. gml และแปลเป็น. kml ในภายหลังเพื่อให้สามารถดูได้ใน GE ฉันสามารถเปิดทั้ง. gml และ. kml เป็นเลเยอร์ใน QGIS แต่ไม่อนุญาตให้ฉันแก้ไขบรรทัด ถ้าฉันบันทึกเลเยอร์. gml หรือ. kml เป็น Shapefile จากนั้นให้โหลด shapefile นั้นเป็นเลเยอร์ QGIS จากนั้นฉันสามารถแก้ไขบรรทัดได้ ฉันได้ลองโหลดไฟล์. kml แล้วสร้างเลเยอร์ shapefile ว่างใหม่การเลือกและคัดลอกบรรทัดบนเลเยอร์. kml และพยายามวางลงในเลเยอร์ shapefile ที่ว่างเปล่าซึ่งไม่สามารถทำได้เช่นกัน ฉันไม่เข้าใจว่า QGIS เข้าใจเส้นบนเลเยอร์ kml ดีพอที่จะบันทึกเป็น Shapefile หรือไม่ทำไมพวกเขาถึงไม่สามารถแก้ไขได้ ฉันรู้สึกว่าฉันต้องคิดถึงบางสิ่งพื้นฐานที่อธิบายสิ่งนี้และบางทีมันอาจมีความหมายอื่นที่ฉันควรเข้าใจ

10 qgis  kml  ogr  editing  gml 

ฉันพยายามแสดงไฟล์แรสเตอร์ความละเอียดขนาดเล็กในเครื่องมือ QGIS Print Manager ฉันต้องการแสดงป้ายกำกับขนาดเล็กที่มีค่าแต่ละเซลล์ (ตำแหน่ง - ประมาณกึ่งกลางของแต่ละเซลล์) ฉันจะทำมันได้อย่างไร มีปลั๊กอินใดบ้างที่สามารถช่วยให้ฉันทำงานนี้ได้สำเร็จ ความคิดของฉันคือการใช้ gdal2xyz โหลดไฟล์ดังกล่าวโดยใช้ปลั๊กอิน CSV แล้วบันทึกเป็นไฟล์เวกเตอร์ น่าเสียดายที่มันใช้งานได้ดีมากเพราะฉันมีไฟล์แรสเตอร์มากมาย อีกแนวคิดหนึ่งคือการใช้ gdal2xyz.py แล้ว ogr2ogr เพื่อบันทึกเป็นไฟล์เวกเตอร์ มีวิธีการชุดกระบวนการนี้หรือไม่ มีโปรแกรมโอเพนซอร์ซใดที่ฉันสามารถทำได้ง่ายกว่านี้อีกไหม?

10 qgis  raster  gdal  ogr 

เราได้รับ DXF ที่มีแอตทริบิวต์ Elevation แต่เราต้องการสร้าง 3D DXF ฉันรู้ว่าฉันสามารถสร้าง 3D dxf จาก shapefile โดยใช้: ogr2ogr -f "DXF" out.dxf in.shp -zfield Elev แต่ไดรเวอร์ ogr ดูเหมือนจะไม่รับแอตทริบิวต์เป็นเขตข้อมูลใน dxf ไม่มีใครรู้ว่าสิ่งนี้เป็นไปได้ ฉันต้องการโซลูชันฟรีหากเป็นไปได้ เรามี FME แต่ฉันกำลังมองหาตัวเลือกอื่น ๆ

10 3d  ogr  fme  dxf 

ฉันได้รับสำเนาหนังสือสุดยอด 'Python Geospatial Development' โดย Erik Westra ( ลิงก์ Amazon ) และฉันกำลังทำงานอยู่ ปัจจุบันสอนให้ฉันโหลดข้อมูลชายฝั่ง GSHHS จากรูปร่างไฟล์ลงในฐานข้อมูล PostGIS เพื่อเตรียมสร้างแอปพลิเคชันเว็บเชิงพื้นที่ ปัญหาของฉันคือ: เมื่อฉันพยายามนำเข้าข้อมูล GSHHS ไปยัง PostGIS นั้นจะถูกปฏิเสธเนื่องจากรูปหลายเหลี่ยมชายฝั่งไม่ได้รับการพิจารณาว่า 'ถูกต้อง' โดยเฉพาะฉันได้รับข้อความแสดงข้อผิดพลาดที่อธิบายถึงรูปหลายเหลี่ยม (แต่ไม่ทั้งหมด) ของแนวชายฝั่งว่าเป็น 'วงไม่ปิด' ฉันเข้าใจว่าข้อผิดพลาดนี้พยายามบอกฉันว่าจุดแรกและจุดสุดท้ายของรูปหลายเหลี่ยมนั้นไม่เหมือนกัน อย่างไรก็ตามสิ่งนี้ไม่เป็นความจริง ฉันได้ตรวจสอบการเป็นตัวแทนของ WKT ของรูปหลายเหลี่ยมหลายอันและถูกต้องแล้ว แน่นอนพวกเขาเริ่มต้นและสิ้นสุดด้วยการประสานงานเดียวกัน รูปหลายเหลี่ยมจะถูกแยกออกจาก shapefiles โดยใช้ไลบรารี OGR และส่งออกแต่ละรูปหลายเหลี่ยมไปยัง WKT ฉันลองสร้างรูปหลายเหลี่ยมใหม่ผ่าน Shapely และทดลองกับ WKB แต่ก็ไม่มีประโยชน์ ฉันได้รับสามารถที่จะโหลดข้อมูลเดียวกันใน PostGIS เป็นตาราง MULTIPOLYGON โดยใช้รถตักดิน …

10 python  postgis  ogr  wkt