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

Python เป็นภาษาโปรแกรมโอเพ่นซอร์สที่ใช้ในโปรแกรม GIS จำนวนมาก

2
การซ้อนทับรูปหลายเหลี่ยมด้วยหุ่นดี
ฉันพยายามจับรูปหลายเหลี่ยมที่ไม่ทับซ้อนกันทั้งหมดที่ระบุด้านล่างโดยใช้ Shapely (รูปหลายเหลี่ยม A, B & C) ยิ่งกว่านั้นฉันหวังว่าจะทำเช่นนั้นโดยไม่ต้องทำซ้ำทดสอบจุดตัดเป็นต้นคำตอบที่ยอมรับสำหรับคำถามนี้เป็นการแสดงออกถึงวิธีการของ PostGIS แต่ดูเหมือนว่า 'สหภาพ' หมายถึงสิ่งที่แตกต่างกันสำหรับผู้คนที่แตกต่างกัน

6
Pure Python Library สำหรับปฏิบัติการทางเรขาคณิต
ไม่มีใครรู้เกี่ยวกับห้องสมุด Python บริสุทธิ์ (ไม่มีการพึ่งพา C หรือ Java) สำหรับการดำเนินการทางเรขาคณิตเช่นบัฟเฟอร์, ยูเนี่ยน, ปริภูมิ, ฯลฯ ? ฉันกำลังมองหาบางอย่างที่เข้ากันได้กับ Google AppEngine

2
รับชื่อฟิลด์ของ shapefiles โดยใช้ GDAL
ฉันใช้ GDAL ใน Python สำหรับการนำเข้า shapefile ฉันต้องการทราบชื่อฟิลด์สำหรับไฟล์ด้วยวิธีปัจจุบันของฉันคือ: fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) แต่ด้วยวิธีนี้ฉันได้รับคุณสมบัติสำหรับเลเยอร์แรก มันหมายความว่าเป็นไปได้หรือไม่ที่เลเยอร์ที่ต่างกันสามารถมีคุณสมบัติที่แตกต่างกันได้? ถ้าไม่เป็นไปได้ที่จะรับชื่อฟิลด์พร้อมกันแทนที่จะเข้าสู่ส่วนลึกนี้หรือไม่? ถ้าใช่มีวิธีใดในการเรียกชื่อฟิลด์ที่ง่ายกว่านี้ไหม
15 python  shapefile  gdal 

3
การแปลงไฟล์ LAS ไปเป็นอาร์เรย์ที่ไม่มีค่า?
ฉันเริ่มเรียนรู้วิธีจัดการข้อมูล LAS ในหลามและต้องการดูว่าคนอื่นจัดการไฟล์ LAS อย่างไร ฉันต้องการอ่านคะแนน (ฉันใช้อาร์เรย์ numpy) และกรองคลาสที่ 1 และ 2 (ไม่จัดประเภทและกราวด์) ไปยังอาร์เรย์ที่แยกต่างหาก ฉันมีรหัสต่อไปนี้ แต่ดูเหมือนจะไม่สามารถรับคะแนนที่กรองได้ # Import modules from liblas import file import numpy as np if __name__=="__main__": '''Read LAS file and create an array to hold X, Y, Z values''' # Get file las_file = r"E:\Testing\ground_filtered.las" # Read …
15 python  lidar  numpy  array  liblas 

3
gdal ImportError ใน python บน Windows
กำลังพยายามนำเข้า gdal สำหรับ Python 2.7.3 บน Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <mo _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swit_helper _mod = imp.load_module('_gdal', …

2
การแยกรูปร่างไฟล์ต่อคุณสมบัติใน Python โดยใช้ GDAL?
เป็นไปได้ไหมที่จะแบ่ง shapefile ต่อคุณสมบัติใน python (ดีที่สุดจะเป็นวิธีการแก้ปัญหาที่ฉันสามารถบันทึกวัตถุเวกเตอร์ผลลัพธ์ลงในหน่วยความจำแทนลงดิสก์ได้ชั่วคราว) เหตุผล: ฉันต้องการใช้ฟังก์ชั่น gdal rasterizeLayer กับหลาย ๆ เซตย่อยของ shapefile ฟังก์ชั่นต้องการวัตถุ osgeo.ogr.Layer mkay ฉันลองไปรอบ ๆ แล้วมันสามารถทำงานได้ดังต่อไปนี้ คุณสามารถรับรูปทรงเรขาคณิตของวัตถุเลเยอร์ gdal ต่อคุณสมบัติดังต่อไปนี้ # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() ตอนนี้ฉันแค่ต้องรู้วิธีการสร้างวัตถุ osgeo.ogr.layer ตามรูปทรงเรขาคณิตนี้ สำหรับการชี้แจง ฉันต้องการฟังก์ชั่นในรหัส ogr …

1
การจัดกลุ่มข้อมูลทางภูมิศาสตร์ตามตำแหน่งของจุดและค่าจุดที่เกี่ยวข้อง
กำหนดจุดข้อมูลด้วยลองจิจูดละติจูดและค่าคุณสมบัติที่สามของจุดนี้ ฉันจะจัดกลุ่มคะแนนเป็นกลุ่ม (ภูมิภาคย่อยทางภูมิศาสตร์) ตามมูลค่าทรัพย์สินได้อย่างไร ฉันค้นหาโดย google และพบว่าปัญหานี้ดูเหมือนว่าจะเรียกว่า "การจัดกลุ่มแบบ จำกัด เชิงพื้นที่" หรือ "การกำหนดภูมิภาค" อย่างไรก็ตามฉันไม่คุ้นเคยกับการจัดการข้อมูลทางภูมิศาสตร์และยังไม่ทราบว่าอัลกอริทึมชนิดใดที่ดีและแพ็คเกจ python / R ใดที่ดีสำหรับงานนี้ หากต้องการให้แนวคิดที่ง่ายขึ้นเกี่ยวกับสิ่งที่ฉันต้องการสมมติว่าแผนการกระจายข้อมูลของฉันมีดังนี้ ดังนั้นแต่ละจุดคือจุด, x คือลองจิจูด, y คือละติจูด, และตารางสีแสดงว่าค่านั้นใหญ่หรือเล็ก ฉันต้องการแบ่งจุดเหล่านั้นออกเป็นภูมิภาคย่อย / กลุ่ม / กลุ่มตามสถานที่และความคล้ายคลึงกันของค่า ชอบสิ่งต่อไปนี้ (ไม่ใช่สิ่งที่ฉันต้องการเพียงแค่แสดงความคิดที่เป็นธรรมชาติ): ดังนั้นฉันจะบรรลุสิ่งนี้ได้อย่างไร

3
การละลายรูปหลายเหลี่ยมขึ้นอยู่กับคุณลักษณะด้วย Python (หุ่นดีฟิโอน่า)?
ฉันพยายามสร้างฟังก์ชั่นที่ทำสิ่งเดียวกันกับที่ QGIS "ละลาย" ฟังก์ชั่น ฉันคิดว่ามันจะง่ายสุด ๆ แต่ก็ดูเหมือนจะไม่ดี ดังนั้นจากสิ่งที่ฉันรวบรวมได้การใช้ฟิโอน่ากับหุ่นดีควรเป็นตัวเลือกที่ดีที่สุดที่นี่ ฉันเพิ่งเริ่มยุ่งเกี่ยวกับไฟล์เวคเตอร์ดังนั้นโลกนี้ก็ค่อนข้างใหม่สำหรับฉันและไพ ธ อนเช่นกัน สำหรับตัวอย่างเหล่านี้ฉันทำงานกับ County shapefile ก่อตั้งขึ้นที่นี่http://tinyurl.com/odfbanu ดังนั้นนี่คือรหัสบางส่วนที่ฉันรวบรวม แต่ไม่สามารถหาวิธีที่จะทำให้พวกเขาทำงานร่วมกันได้ สำหรับตอนนี้วิธีที่ดีที่สุดของฉันคือต่อไปนี้ขึ้นอยู่กับ: https://sgillies.net/2009/01/27/a-more-perfect-union-continued.html มันใช้งานได้ดีและฉันได้รับรายชื่อของ 52 รัฐเป็นรูปทรงเรขาคณิตที่ดี โปรดแสดงความคิดเห็นหากมีวิธีที่ตรงไปข้างหน้ามากขึ้นในการทำส่วนนี้ from osgeo import ogr from shapely.wkb import loads from numpy import asarray from shapely.ops import cascaded_union ds = ogr.Open('counties.shp') layer = ds.GetLayer(0) #create a list of unique …

2
กำลังเข้าร่วม (concat) รายการของ dataframes ที่คล้ายกันใน geopandas หรือไม่
ฉันมีรายการของ dataframes โครงสร้างเดียวกันทั้งหมดและทั้งหมดเป็นรูปหลายเหลี่ยม ฉันจะเพิ่มพวกเขาทั้งหมดเข้าด้วยกันโดยใช้ geopandas ได้อย่างไร
15 python  geopandas 

1
รับค่าพิกเซลที่จุดเดียวโดยใช้ rasterio
ในการรับค่าพิกเซลเพียงจุดเดียวใน raster โดยใช้ rasterio มีตัวอย่างอยู่ที่นี่: https://github.com/mapbox/rasterio/pull/275 อย่างไรก็ตามมี API โดยตรงภายใน rasterio (ไม่ใช่ cli) ซึ่งสามารถใช้เพื่อแยกค่าที่จุดเดียวใน raster หรือไม่? - แก้ไข with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src: x = (src.bounds.left + src.bounds.right) / 2.0 y = (src.bounds.bottom + src.bounds.top) / 2.0 vals = src.sample((x, y)) …
14 python  numpy  rasterio 

2
วิธีการคำนวณกล่องขอบเขตของวัตถุ geojson โดยใช้ python หรือ javascript
ฉันมีแอพ django ที่มีแผ่นพับที่ส่วนหน้าและฉันจำเป็นต้องคำนวณ bbox ของวัตถุ geojson เพื่อให้ฉันสามารถส่งพิกัดไปยัง map.fitBounds() ฉันลองใช้ libs แล้วลองใช้getBoundsวิธีการของวัตถุบางอย่างใน Leaflet เช่นFeatureGroupแต่มันก็บ่นว่าขอบเขตไม่ได้ถูกกำหนดไว้ ใครช่วยชี้ให้ฉันเห็นวิธีแก้ปัญหาง่ายๆ

1
กำลังแปลงทวีตของ Twitter เป็นจุดหรือไม่
ฉันได้ทำการวิจัยเกี่ยวกับการใช้ R เพื่อ datamine twitter แต่ฉันไม่พบคำตอบหรือแบบฝึกหัดที่ดีสำหรับคำถามของฉัน ฉันสนใจดึงทวีตจากทวิตเตอร์ด้วยแฮชแท็กภายในกรอบเวลาที่กำหนดและวางแผนตำแหน่งของทวีตเหล่านั้นบนแผนที่ใน QGIS หรือ ArcMap ฉันรู้ว่าทวีตสามารถมีตำแหน่งทางภูมิศาสตร์เชื่อมโยงกับพวกเขา แต่ฉันจะดึงข้อมูลนี้ได้อย่างไรในตอนแรก?
14 python  point  twitter 

2
วิธีการวาดรูปหลายเหลี่ยมจากคอนโซลหลาม?
ฉันได้เยี่ยมชมและเยี่ยมชมหน้าใหม่เกี่ยวกับการจัดการรูปทรงเรขาคณิตในตำราอาหาร PyQGIS: http://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.htmlแต่ดูเหมือนว่าจะไม่เข้าใจวิธีการรับ รูปหลายเหลี่ยมที่วาดจากคอนโซล Python ใครช่วยได้บ้าง

5
ภาษาการเขียนโปรแกรมที่แนะนำสำหรับการสำรวจระยะไกล?
ฉันเริ่มการศึกษาที่หวังว่าจะนำไปสู่อาชีพที่ยาวนานในฐานะผู้เชี่ยวชาญด้านการรับรู้ระยะไกล ขณะนี้ฉันกำลังทำงานกับ ArcGIS สำหรับบางแอปพลิเคชันและเรียนรู้ ENVI สำหรับผู้อื่น ฉันได้ตระหนักว่ามีความจำเป็นที่ฉันต้องเรียนรู้ภาษาการเขียนโปรแกรมและติดอยู่กับการเลือกระหว่าง IDL และ Python ฉันชอบที่จะได้ยินภาษาการเขียนโปรแกรมที่ชุมชนแนะนำสำหรับการสำรวจระยะไกลของกระบวนการพื้นผิว

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

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