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

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


13
ติดตั้ง GDAL ด้วย Python บน windows หรือไม่
ทุกคนสามารถอธิบายวิธีการติดตั้ง GDAL / OGR ด้วย Python บน Windows ได้ไหม ฉันมี Windows Vista และฉันได้ลองทำตามข้อมูลในเว็บไซต์และดูเหมือนจะไม่ผูกไฟล์ exe ให้ฉัน ใครสามารถอธิบายกระบวนการรวมถึงลิงค์ไปยังไฟล์ / โฟลเดอร์ที่ฉันต้องการ ตอนนี้ฉันพยายามเรียกใช้การตั้งค่า gdal ด้วย minGW แต่สิ่งนี้ก็ล้มเหลวด้วย:

11
ทางเลือกอื่นในการใช้ Arcpy
ฉันดูเหมือนจะใช้แพคเกจเว็บไซต์ Arcpy ของ ESRI สำหรับการประมวลผลหลามทางภูมิศาสตร์ของฉัน สำหรับเครดิตของ ESRI เหล่านี้เป็นชุดเครื่องมือที่ยอดเยี่ยมที่สามารถช่วยให้บรรลุข้อตกลงอย่างมาก อย่างไรก็ตามฉันต้องการสร้างสคริปต์การประมวลผลทางภูมิศาสตร์นอกโดเมน ESRI Arcpy ตัวอย่างเช่นถ้าฉันต้องการคลิปแรสเตอร์เป็นรูปหลายเหลี่ยมฉันจะเริ่มต้นด้วยสคริปต์ต่อไปนี้จากESRI : # Import system modules import arcpy from arcpy import env from arcpy.sa import * # Set environment settings env.workspace = "C:/sapyexamples/data" # Set local variables inRaster = "elevation" inMaskData = "mask.shp" # Check out the ArcGIS Spatial …

23
Python GDAL แพ็คเกจไฟล์ส่วนหัวหายไปเมื่อติดตั้งผ่าน pip
ฉันถามคำถามนี้เกี่ยวกับ Stack Overflowก่อนที่ฉันจะรู้ว่ามีเว็บไซต์นี้อยู่ เนื่องจากไซต์นี้เหมาะสมกว่าฉันจะถามที่นี่ด้วย ฉันกำลังพยายามติดตั้ง gdal จาก pip pip install gdalในสภาพแวดล้อมเสมือนจริง (Ubuntu) มันล้มเหลวเพราะหาไม่เจอcpl_port.h extensions/gdal_wrap.cpp:2853:22: fatal error: cpl_port.h: No such file or directory compilation terminated อย่างไรก็ตาม GDAL /usr/include/gdal/cpl_port.hมีการติดตั้งอย่างถูกต้องและไฟล์ส่วนหัวจะอยู่ที่ มีตัวแปรสภาพแวดล้อมสำหรับ GDAL ที่ต้องตั้งค่าเพื่อให้ pip ค้นหาไฟล์ส่วนหัวหรือไม่

3
ปรับรูปหลายเหลี่ยมให้เรียบลงในแผนที่รูปร่าง
นี่คือแผนที่รูปร่างที่มีรูปหลายเหลี่ยมหลายระดับ ให้ถามวิธีทำให้รูปหลายเหลี่ยมทำให้จุดยอดทั้งหมดคงอยู่ในตำแหน่งที่แน่นอนได้อย่างไร จริง ๆ แล้วรูปร่างนั้นถูกสร้างขึ้นบนข้อมูลกริดคุณอาจแนะนำให้ปรับข้อมูลกริดให้เรียบและดังนั้นผลลัพธ์ที่ได้จะราบรื่นขึ้น โปรดทราบว่านี่ไม่ทำงานตามที่ฉันต้องการเนื่องจากฟังก์ชั่นการปรับให้เรียบเช่นตัวกรองแบบเกาส์จะลบชุดข้อมูลขนาดเล็กและจะเปลี่ยนช่วงของตัวแปรที่สามเช่นความสูงที่ไม่ได้รับอนุญาตในแอปพลิเคชันของฉัน อันที่จริงฉันกำลังมองหาชิ้นส่วนของรหัส (โดยเฉพาะอย่างยิ่งในPython ) ซึ่งสามารถทำการปรับให้เรียบของรูปหลายเหลี่ยม 2 มิติ (ชนิดใดก็ได้: นูนเว้าการตัดกันด้วยตนเองเป็นต้น) ไม่เจ็บปวดพอสมควร (ลืมหน้ารหัส) และแม่นยำ FYI มีฟังก์ชั่นใน ArcGISที่ทำสิ่งนี้ได้อย่างสมบูรณ์แบบ แต่การใช้แอพพลิเคชั่นเพื่อการค้าของบุคคลที่สามไม่ใช่ตัวเลือกของฉันสำหรับคำถามนี้ 1) Scipy.interpolate: ในขณะที่คุณเห็นเส้นโค้งที่เกิด (สีแดง) ไม่เป็นที่พอใจ! 2) นี่คือผลลัพธ์โดยใช้รหัสที่ให้ไว้ในที่นี้ มันใช้งานไม่ได้! 3) สำหรับฉันทางออกที่ดีที่สุดควรเป็นสิ่งต่อไปนี้ซึ่งเป็นรูปสี่เหลี่ยมจัตุรัสที่ถูกปรับให้เรียบโดยค่อยๆเปลี่ยนค่าเพียงค่าเดียว ฉันหวังว่าแนวคิดที่คล้ายกันเพื่อทำให้รูปหลายเหลี่ยมเรียบขึ้น ตอบสนองเงื่อนไขที่เส้นโค้งผ่านจุด: 4) นี่คือการดำเนินการของฉัน "ความคิดของคนจรจัด" ทีละบรรทัดในPythonกับข้อมูลของเขา อาจมีข้อบกพร่องบางอย่างเนื่องจากผลลัพธ์ไม่ดี K = 2 เป็นหายนะดังนั้น k> = 4 5) ฉันลบจุดหนึ่งไปยังตำแหน่งที่มีปัญหาและตอนนี้เดือยที่เกิดขึ้นนั้นเหมือนกับของ whuber's แต่ก็ยังคงมีคำถามว่าทำไมวิธีการนี้จึงใช้ไม่ได้กับทุกกรณี 6) …

8
การแปลงพิกัดที่ฉายเป็น lat / lon โดยใช้ Python หรือไม่
ฉันถือว่านี่เป็นคำถามพื้นฐาน แต่ดูเหมือนไม่สามารถค้นหาหรือรู้จักวิธีแก้ปัญหา ไซต์นี้ส่งคืน Point: X: -11705274.6374 Y: 4826473.6922 เมื่อคุณค้นหาด้วยค่าคีย์แรกของ 000090 เป็นตัวอย่าง ฉันเดาว่านี่เป็นการอ้างอิงเชิงพื้นที่และฉันก็เข้าใจว่ามันคืออะไร ฉันกำลังมองหาคำแนะนำหรือตัวอย่างวิธีแปลงเป็น Latitude และลองจิจูดโดยใช้ Python

17
เรียนรู้การเขียนโปรแกรมหลามด้วยเป้าหมาย GIS ทั่วไปในใจ?
ในความเห็นของคุณหนังสือ / ไซต์ที่ดีที่สุดในการเรียนรู้ Python โดยมีเป้าหมายใน GIS คืออะไร โดย "ดีที่สุด" ฉันหมายถึง: ไม่นานมาก (หนังสือ) ง่ายต่อการเข้าใจ (หนังสือ / เว็บไซต์) ตัวอย่างการปฏิบัติที่ดี (หนังสือ / เว็บไซต์) สำหรับคำตอบเฉพาะสำหรับการเรียนรู้วิธีการใช้ ArcPy site-package สำหรับ Python ในการปรับแต่ง ArcGIS สำหรับเดสก์ท็อปมีคำถามและคำตอบแยกต่างหาก: อะไรคือทรัพยากรสำหรับการเรียนรู้ ArcPy

1
รับค่าพิกเซลของแรสเตอร์ GDAL ภายใต้จุด OGR โดยไม่ต้อง NumPy?
ฉันกำลังทำงานกับแบบจำลองการคำนวณของความหลากหลายของเรณูป่าในภูมิประเทศ ตัวแบบเสร็จสมบูรณ์แล้วและตอนนี้ฉันกำลังดิ้นรนกับขั้นตอนหลังการประมวลผล ฉันมีเรณูซัพพลาย GDAL ของฉันที่มีลักษณะคล้ายนี้ (สีที่เบากว่าหมายถึงการเยี่ยมชมของเรณูที่พิกเซล): และฉันมี OGR รูปร่างไฟล์ของจุดที่เป็นตัวแทนของตัวอย่างสถานที่ในแนวนอน: ฉันกำลังพยายามวิเคราะห์จุดพิกเซลภายใต้จุดเหล่านี้ แต่ต้องทำเช่นนั้นฉันต้องสามารถแยกค่าของพิกเซลที่อยู่ใต้จุดใดจุดหนึ่งได้ เป็นไปได้หรือไม่ที่จะแยกค่าของพิกเซลภายใต้จุดโดยใช้ OGR และ GDAL ผ่าน Python เท่านั้น ฉันต้องการหลีกเลี่ยงการอ่านแรสเตอร์ทั้งหมดลงในหน่วยความจำReadAsArray()เพราะเอาท์พุทแรสเตอร์ของฉันมีขนาดใหญ่มาก (ใหญ่เกินไปที่จะใส่ลงในหน่วยความจำ) ฉันสังเกตเห็นโพสต์นี้ซึ่งคล้ายกัน แต่ต้องมีการเรียกบรรทัดคำสั่ง
45 python  gdal  ogr  pixel 

5
วิธีการโทร gdal_translate จากรหัสหลาม?
เป็นไปได้ไหมที่จะใช้ gdal API โทรgdal_translateจากรหัส Python? ฉันไม่ได้หมายถึงเพียงแค่เรียกใช้ gdal_translate.exe จากระบบไฟล์ แต่จะเรียกมันว่าอย่างใดในรหัสดังนั้นฉันไม่จำเป็นต้องรู้ไดเรกทอรีที่แน่นอนว่าไฟล์ gdal_translate นั้นอยู่ในอะไร
40 python  gdal 

3
เป็นไปได้หรือไม่ที่จะดูเนื้อหาของ Shapefile โดยใช้ Python ที่ไม่มีลิขสิทธิ์ ArcMap?
ฉันสงสัยว่ามันเป็นไปได้ที่จะดูเนื้อหาของ shapefile โดยใช้ Python โดยไม่ต้องมีและสิทธิ์การใช้งาน ArcMap สถานการณ์คือคุณสามารถสร้าง shapefiles จากหลาย ๆ แอปพลิเคชั่นไม่เพียง แต่จากซอฟต์แวร์ ESRI ฉันต้องการสร้างสคริปต์ Python ที่ตรวจสอบการอ้างอิงเชิงพื้นที่ประเภทสถานที่ชื่อคุณลักษณะและคำจำกัดความและเนื้อหาของเขตข้อมูลในรูปแบบไฟล์และเปรียบเทียบกับชุดของค่าที่ยอมรับได้ ฉันต้องการให้สคริปต์นี้ทำงานแม้ว่าองค์กรจะไม่มีใบอนุญาต ESRI ในการทำสิ่งนี้คุณต้องใช้ ArcPy หรือคุณสามารถขุดลงในไฟล์รูปร่างโดยไม่ต้องใช้ ArcPy ได้หรือไม่?
40 python  shapefile 

4
การแม็พใจความอย่างง่ายของ Shapefile โดยใช้ Python?
ฉันต้องการแสดงข้อมูลทางภูมิศาสตร์ใน Python โดยไม่ต้องใช้ ArcGIS / ArcPy และสร้างแผนที่ บนอินเทอร์เน็ตฉันพบวิธีทำแผนที่เฉพาะเรื่องโดยใช้ Python : นี่คือรหัสที่ฉันลอง: import shpUtils import matplotlib.pyplot as plt shpRecords = shpUtils.loadShapefile('C:\\Users\\shapefile.shp') for i in range(0,len(shpRecords)): x = [] y = [] for j in range(0,len(shpRecords[i]['shp_data']['parts'][0]['points'])): tempx = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['x']) tempy = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['y']) x.append(tempx) y.append(tempy) plt.fill(x,y) plt.axis('equal') plt.title("Testing") plt.show() อย่างไรก็ตามเมื่อฉันเรียกใช้มันจะให้สีแบบสุ่มแก่ฉัน ถ้าฉันต้องการเห็นคอลัมน์ของแฟ้มข้อมูลรูปร่างฉันจะใช้สิ่งนี้โดยใช้รหัสที่คล้ายกันได้อย่างไร สิ่งนี้ไม่ชัดเจนในลิงก์ที่ให้ไว้ข้างต้นซึ่งเขากล่าวถึงการใช้สีเท่านั้น ... …

7
เพื่อนบ้านที่ใกล้ที่สุดระหว่างเลเยอร์จุดและเลเยอร์บรรทัด [ปิด]
ฉันได้ถามคำถามนี้หลายครั้งใน stackoverflow และ irc ระหว่าง #qgis และ #postgis และฉันยังพยายามที่จะเขียนโค้ดหรือใช้มันด้วยตนเองใน postgis โดยไม่มีคำตอบจริง โดยใช้การเขียนโปรแกรม (ส่วนใหญ่เป็นงูใหญ่) ฉันต้องการวาดเส้นจากจุดเลเยอร์ไปยังเส้นโครงของมันในบรรทัดที่ใกล้ที่สุดของเส้นหรือชั้นรูปหลายเหลี่ยม ณ ตอนนี้ข้อมูลของฉันส่วนใหญ่อยู่ในรูปแบบของ ESRI และรูปแบบโพสต์จิส; อย่างไรก็ตามฉันควรจะอยู่ห่างจากโซลูชัน postgis เนื่องจากฉันเป็นผู้ใช้ shp + qgis ทางออกที่ดีคือการใช้ GDAL / OGR กับ python หรือไลบรารี่ที่คล้ายกัน การใช้ห้องสมุด GDAL / OGR จะเริ่มได้ที่ไหน เป็นไปได้ไหมที่จะให้แผนโซลูชัน ฉันสามารถใช้ NetworkX เพื่อทำการวิเคราะห์เพื่อนบ้านที่ใกล้ที่สุดได้หรือไม่ เป็นไปได้จริงหรือ หากทำได้ง่ายขึ้นคะแนนสามารถเชื่อมต่อกับจุดสิ้นสุดของกลุ่มแทนที่จะเป็นจุดที่คาดการณ์ไว้

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

3
ไวยากรณ์และการใช้งานที่เหมาะสมสำหรับเวิร์กสเปซ ArcGIS“ in_memory” คืออะไร?
ฉันไม่ชัดเจนว่าจะใช้in_memoryพื้นที่ทำงานใน ArcGIS อย่างไรดีที่สุด ทรัพยากรที่ดีที่สุดที่ฉันได้พบเกี่ยวกับเรื่องนี้จาก ESRI ในการใช้พื้นที่ทำงานในหน่วยความจำ น่าเสียดายที่ส่วนช่วยเหลือไม่ได้ให้ข้อมูลมากนัก นอกจากนี้ยังมีการโพสต์บนเว็บไซต์นี้ว่าในเวลาสั้น ๆ สัมผัสกับเรื่องที่นี่และที่นี่ คำถามเฉพาะของฉัน: การใช้งานและไวยากรณ์ที่เหมาะสมสำหรับการใช้in_memoryพื้นที่ทำงานในสคริปต์ที่ใช้ArcGIS / arcpy คืออะไร? เป็นin_memoryพื้นที่ทำงานเช่นเดียวกับยกตัวอย่างเช่นการสร้างชั้นใช้arcpy.MakeFeatureLayer_management()? มีมาตรฐานใด ๆ เช่นการลบin_memoryพื้นที่ทำงานในตอนท้ายของสคริปต์หรือไม่

10
จะเปลี่ยนคุณสมบัติทั้งหมดในชุดข้อมูลเวกเตอร์ได้อย่างไร
สมมุติว่าฉันรวม Shapefile เข้าด้วยกันและฟีเจอร์ทั้งหมดจะเปลี่ยนจุดยอดเป็นจำนวนคงที่ อะไรคือวิธีที่ง่ายที่สุดในการเปลี่ยนฟีเจอร์ทั้งหมด (x, y)ของจุดยอด) โดยการเลื่อนตามอำเภอใจ? ฉันมีไฟล์มากมายที่ฉันจะใช้การแก้ไขนี้กับดังนั้นคำตอบ Bash / OGR จะเป็นที่ต้องการ :) สุดท้ายผมสิ้นสุดการใช้ Spatialite ShiftCoordsนี้มันมีฟังก์ชั่นที่ดี อย่างไรก็ตามด้ายมีข้อมูลมาก! ขอบคุณทุกคน!

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