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

ฉันมีไฟล์รูปร่างที่มีรูปหลายเหลี่ยม และฉันมีไฟล์แรสเตอร์ระดับโลก ฉันต้องการซ้อนรูปหลายเหลี่ยมของ shapefile ลงบนตารางแรสเตอร์และคำนวณค่าเฉลี่ย raster สำหรับรูปหลายเหลี่ยมแต่ละรูป ฉันจะทำสิ่งนี้โดยใช้ GDAL เขียนผลลัพธ์ไปยังรูปร่างไฟล์ได้อย่างไร

17 python  raster  shapefile  gdal  spatial-statistics 

ฉันมี GeoPandas Dataframe ง่าย ๆ : ฉันต้องการอัปโหลด GeoDataframe นี้ไปยังตาราง PostGIS ฉันมีการตั้งค่าฐานข้อมูลพร้อมส่วนขยาย PostGIS แล้ว แต่ไม่สามารถเพิ่ม Dataframe นี้เป็นตารางได้ ฉันได้ลองทำสิ่งต่อไปนี้แล้ว: engine = <> meta = MetaData(engine) eld_test = Table('eld_test', meta, Column('id', Integer, primary_key=True), Column('key_comb_drvr', Text), Column('geometry', Geometry('Point', srid=4326))) eld_test.create(engine) conn = engine.connect() conn.execute(eld_test.insert(), df.to_dict('records'))

17 python  postgis  geopandas 

ขณะนี้ฉันกำลังทำงานเพื่อนำเข้าข้อมูลสภาพอากาศ CANGRID (จัดให้เป็น Surfer Grid ascii, ไฟล์ ".grd") ไปยัง ArcGIS ตารางมีขนาด 95 แถวคูณ 125 คอลัมน์ เมทาดาทาให้ละติจูด / ลองจิจูดของแหล่งกำเนิด (มุมซ้ายล่าง), ขนาดเซลล์ (50 กม.) เช่นเดียวกับการบันทึกภาพเป็นขั้วภาพสามมิติที่มีจุดศูนย์กลางกลาง (110 องศา W) และละติจูดของจุดกำเนิด (60 องศา N) หลังจากครั้งแรกที่พยายามแปลง. grd เป็น rasc as .ascii และ. flt ไม่สำเร็จฉันได้จัดการใช้ GDAL เพื่อกำหนดขอบเขตและการฉายภาพอย่างไรก็ตามชุดข้อมูลไม่ถูกต้องสอดคล้องกับขอบเขตของพื้นที่ที่ต้องการ ดูภาพด้านล่าง มีการเปลี่ยนรูปทางภูมิศาสตร์ที่ยอมรับได้สำหรับ stereographic ขั้วโลกที่สามารถอธิบายการขาดการจัดตำแหน่งนี้ได้หรือไม่? ตัวอย่างเช่นมีปัจจัยการแปลงหรือการหมุนเฉพาะที่ฉันควรใช้หรือไม่ ไฟล์ตัวอย่างจากชุดข้อมูลอยู่ที่นี่: "t201113.grd" นี่คือรหัสที่ฉันใช้ใน GDAL …

16 arcgis-desktop  python  gdal  coordinate-system 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน9 เดือนที่ผ่านมา หากไม่รวม ArcPy จะมีห้องสมุดไพ ธ อนใดบ้างที่สามารถทำการประมวลผลทางภูมิศาสตร์เช่นบัฟเฟอร์ / ตัดกับรูปร่างของไฟล์ได้?

16 python  buffer  software-recommendations 

ฉันสนใจที่จะเรียนรู้วิธีการใช้อาร์เรย์ NumPy เพื่อเพิ่มประสิทธิภาพการประมวลผลทางภูมิศาสตร์ งานส่วนใหญ่ของฉันเกี่ยวข้องกับ "ข้อมูลขนาดใหญ่" ซึ่งการประมวลผลทางภูมิศาสตร์มักใช้เวลาหลายวันเพื่อทำงานบางอย่างให้สำเร็จ จำเป็นต้องพูดฉันสนใจมากในการเพิ่มประสิทธิภาพการปฏิบัติเหล่านี้ ArcGIS 10.1 มีฟังก์ชั่น NumPy จำนวนมากที่สามารถเข้าถึงได้ผ่านทาง arcpy รวมไปถึง: NumPyArrayToFeatureClass (arcpy.da) RasterToNumPyArray (arcpy) TableToNumPyArray (arcpy.da) ตัวอย่างเช่นสมมติว่าฉันต้องการเพิ่มประสิทธิภาพการประมวลผลเวิร์กโฟลว์เข้มข้นต่อไปนี้โดยใช้อาร์เรย์ NumPy: แนวคิดทั่วไปที่นี่คือมีจุดจำนวนมากที่ใช้เวกเตอร์ที่เคลื่อนที่ผ่านทั้งการดำเนินการตามเวกเตอร์และแรสเตอร์ทำให้เกิดชุดข้อมูลแรสเตอร์เลขฐานสองจำนวนเต็ม ฉันจะรวมอาร์เรย์ NumPy เพื่อเพิ่มประสิทธิภาพเวิร์กโฟลว์ประเภทนี้ได้อย่างไร

16 arcgis-desktop  python  arcpy  numpy  big-data 

ฉันต้องการสร้างไฟล์ GeoJSON โดยทางโปรแกรมโดยใช้รูปหลายเหลี่ยมจาก shapefile แต่เพิ่มคุณสมบัติจากแอปพลิเคชันของฉันเอง สิ่งนี้สามารถทำได้อย่างง่ายดายสำหรับไฟล์รูปร่าง: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the list of fields names associated with the nodes. data is a list of lists whose first element is the geocode and the remaining elements are values of the …

16 python  geojson  ogr  fiona 

ฉันมีซอฟต์แวร์บางส่วน (ไม่ใช่ Arc) ที่ทำงานทุกคืนบนเครื่องที่ไม่ได้ใช้งานจริงซึ่งอัพเดตคุณสมบัติ shapefile จากฐานข้อมูลภายนอก ฉันต้องการลบฟีเจอร์ทั้งหมดใน shapefile เป็นระยะ ๆ (ไม่ใช่ไฟล์เองซึ่งจะต้องคงอยู่) และให้ซอฟต์แวร์ "สร้าง" shapefile ใหม่ขึ้นมาใหม่ ฉันต้องการทำให้กระบวนการนี้เป็นแบบอัตโนมัติ ฉันไม่มีซอฟต์แวร์ GIS ใด ๆ ติดตั้งอยู่ในเครื่องนั้น ฉันหวังว่าฉันจะเขียนสคริปต์เป็นประจำใน Python ซึ่งจะลบคุณลักษณะดังกล่าวโดยอัตโนมัติเหมือนกับเครื่องมือลบคุณสมบัติทางภูมิศาสตร์ของ Arc ในการประมวลผล มีโมดูล Python ใดบ้างที่จะอนุญาตให้ฉันทำสิ่งนี้? โดยเฉพาะอย่างยิ่งโอเพนซอร์ส? ผมมองไปที่หุ่นดีและPyShpแต่ไม่เห็นอะไรที่จะช่วยให้ฉันลบคุณลักษณะในมวลหรือที่ตรงกับประโยคที่ พวกเขาสามารถเขียนคุณสมบัติและวิเคราะห์ แต่ไม่เห็นคุณลักษณะการลบที่ใดก็ได้ แน่นอนฉันต้องคิดถึงบางสิ่งที่ง่าย ... แก้ไข:ฉันมี 35 โฟลเดอร์ (ขอบเขตทางภูมิศาสตร์ที่แตกต่างกันทั้งหมดในการฉายของพวกเขาเอง) ด้วย 35-65 รูปร่างไฟล์ที่มีประมาณ 1,000 รูปร่างไฟล์ที่จะจัดการกับ

16 python  shapefile  open-source-gis 

ฉันกำลังทำงานกับ Add-in ของ ArcMap ใน C # จากรหัส C # ฉันได้เรียกใช้งานสคริปต์ Python ตอนนี้เพื่อเรียกใช้สคริปต์เหล่านี้ฉันมีเส้นทางหลามแบบตายตัว แต่นี่ไม่ใช่พกพา ดังนั้นฉันต้องการได้รับเส้นทางของการปฏิบัติการ Python จากรหัสและใช้งาน คำถาม: ฉันจะได้รับเส้นทางของการปฏิบัติการของ Python ที่ใช้โดย ArcMap จากรหัส C # ได้อย่างไร? แก้ไข: จากคำแนะนำของคุณตอนนี้ฉันใช้ "สภาพแวดล้อมของเส้นทาง" เพื่อรับเส้นทาง Python //get python path from environtment variable string GetPythonPath() { IDictionary environmentVariables = Environment.GetEnvironmentVariables(); string pathVariable = environmentVariables["Path"] as string; …

16 arcgis-10.0  python  c# 

ฉันใหม่จริง ๆ กับ Python และฉันต้องการทราบว่ามีวิธีที่รวดเร็วในการรับค่าเซลล์ของพิกเซลแรสเตอร์ทีละพิกเซลและพิกัด (พิกัด XY แผนที่ของศูนย์กลางของแต่ละพิกเซล) โดยใช้ Python ใน ArcGIS 10 หรือไม่ เพื่ออธิบายเพิ่มเติมนี้ฉันต้องได้รับ map X, map Y และค่าเซลล์ของพิกเซลแรกและกำหนดค่าทั้งสามให้กับตัวแปรสามตัวและทำซ้ำขั้นตอนนี้สำหรับส่วนที่เหลือของพิกเซลอื่น ๆ (วนรอบแรสเตอร์ทั้งหมด) ฉันคิดว่าฉันต้องอธิบายคำถามของฉันให้มากขึ้น ปัญหาคือฉันต้องได้รับตำแหน่ง XY ของพิกเซลของแรสเตอร์แรกและรับค่าเซลล์ของแรสเตอร์อื่น ๆ ที่สอดคล้องกับตำแหน่ง XY นั้น กระบวนการนี้ควรวนซ้ำทุกพิกเซลของแรสเตอร์แรกโดยไม่ต้องสร้างเชพไฟล์จุดกึ่งกลางใด ๆ เนื่องจากมันต้องใช้เวลามากจริงๆเพราะฉันต้องจัดการแรสเตอร์ที่มีเกือบ 8 พันล้านพิกเซล นอกจากนี้ฉันต้องทำสิ่งนี้โดยใช้ Python ใน ArcGIS 10 @ James: ขอบคุณมากสำหรับคำแนะนำของคุณ ใช่สิ่งนี้จะใช้ได้กับแรสเตอร์เดียว แต่ฉันต้องเก็บค่าเซลล์สำหรับแรสเตอร์อื่น ๆ ด้วยเช่นกัน ปัญหาคือหลังจากได้รับพิกัด X และ …

16 arcgis-10.0  python  raster 

ขออภัยหากคำถามต่อไปนี้ค่อนข้างงี่เง่า แต่ฉันเพิ่งใหม่กับสิ่ง GIS ทั้งหมดนี้ ฉันพยายามแปลงภาพ geoTiff ที่คาดการณ์ไว้เป็น WGS84 โดยใช้ gdal ใน python ฉันได้พบโพสต์ที่แสดงกระบวนการในการแปลงคะแนนภายใน GeoTiff ที่ฉายโดยใช้บางอย่างที่คล้ายกับต่อไปนี้: from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], …

16 python  gdal  coordinate-system  wgs84 

งูใหญ่แพคเกจ geopyมีสองวัดระยะทางเทคนิค: เกรทเซอร์เคิและVincenty ของสูตร >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd >>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg >>> vincenty(p1, p2).meters 429.16765838976664 >>> great_circle(p3, p4).meters 428.4088367903001 อะไรคือความแตกต่าง? ต้องการการวัดระยะทางแบบใด

16 python  distance  great-circle  geopy  vincenty-formulae 

ด้วยgeopandas GeoDataFrameรูปหลายเหลี่ยมที่มีชุดของฉันต้องการที่จะได้รับพื้นที่ในกม. sq ของแต่ละคุณสมบัติในรายการของฉัน นี่เป็นปัญหาที่พบได้บ่อยและวิธีแก้ปัญหาที่แนะนำตามปกติในอดีตเคยใช้shapelyและpyprojโดยตรง (เช่นที่นี่และที่นี่ ) มีวิธีการทำสิ่งนี้อย่างบริสุทธิ์geopandasหรือไม่?

16 python  area  shapely  geopandas  pyproj 

ฉันใช้ข้อมูล GIS ของเคาน์ตีเพื่ออัพเดทถนนใน OpenStreetMap ฉันมีสำเนาที่พวกเขาตีพิมพ์เมื่อปีที่แล้วและฉบับที่เพิ่งเผยแพร่และฉันต้องการค้นหา LineStrings ทั้งหมดที่มีคุณสมบัติหรือรูปทรงเรขาคณิตของพวกเขาเปลี่ยนไป สิ่งนี้จะช่วยฉันในการทำให้แน่ใจว่าถนนสายใหม่และถนนที่ได้รับการดัดแปลงนั้นได้รับการอัพเดตใน OpenStreetMap ฉันต้องการทำสิ่งนี้โดยใช้ซอฟต์แวร์ FOSS เช่น QGIS หรือ Python / OGR กลุ่มถนนควรมีตัวระบุที่ไม่ซ้ำกันดังนั้นความคิดเดียวของฉันคือการเขียนสคริปต์ Python ที่เปิดทั้งสอง Shapefiles ค้นหากลุ่มใด ๆ ในใหม่ที่ไม่ได้อยู่ในกลุ่มเก่า (เพิ่มกลุ่ม) และกลุ่มย้อนกลับ (ลบกลุ่ม) แล้ววนซ้ำคู่ที่ตรงกันเพื่อเปรียบเทียบพิกัดที่เป็นส่วนประกอบเพื่อดูว่ามีอะไรเปลี่ยนแปลงหรือไม่ นี่เป็นวิธีที่ดีหรือไม่? มีวิธีที่ง่ายกว่า?

15 qgis  python  linestring  change-detection 

สั้น: ฉันต้องการเปลี่ยนความละเอียดของแรสเตอร์และปรับทางลาดสีเทาตามที่แสดงในภาพตะโกน การตั้งค่าคือการใช้ GDAL, PIL หรือ Numpy คำอธิบาย: ฉัน Kriging จุดลงไปในเชิงภาพที่มีความละเอียดการส่งออกของ 20 เมตรกับGeostatistical ห้องสมุดที่มีประสิทธิภาพสูง ฉันไม่ต้องการเปลี่ยนความละเอียดเอาต์พุตเนื่องจากเวลาการแก้ไขเพิ่มขึ้นเป็นทวีคูณ ด้วยความละเอียดนี้ภาพที่ส่งออกจะน่าเกลียด (pixelated และ aliased) ฉันไม่รู้ว่ามันถูกต้องในแนวความคิดหรือไม่ แต่ฉันต้องการภาพที่นุ่มนวลเหมือนในตัวอย่างร้อง มันเหมือนกับการ 'ปรับกำลังภาพ' ใหม่ให้เป็นภาพที่มีความละเอียดมากกว่า ฉันใช้ python ดังนั้นการตั้งค่าของฉันคือ GDAL, Python Imaging Library หรือ Numpy คำตอบอาจเป็นเชิงทฤษฎีเช่นชี้ให้เห็นชื่ออัลกอริทึมหรือแนวคิดของการดำเนินการชนิดนี้ ที่มา: ปลายทาง: แก้ไขผลลัพธ์ด้วย gdalwarp cubic spline:

15 python  raster  gdal 

ฉันกำลังมองหาทางออกที่แตกต่างและสง่างามกว่าสำหรับปัญหาสถิติเชิงพื้นที่ ข้อมูลดิบประกอบด้วยพิกัด xy สำหรับต้นไม้แต่ละต้น (เช่นแปลงเป็นไฟล์ point .shp) แม้ว่าจะไม่ได้ใช้ในตัวอย่างนี้ต้นไม้ทุกต้นยังมีรูปหลายเหลี่ยมที่สอดคล้องกัน (เช่น. shp) ซึ่งแสดงถึงเส้นผ่านศูนย์กลางมงกุฎ รูปภาพสองภาพทางด้านซ้ายแสดงการประมาณความหนาแน่นเคอร์เนลแนวนอน (KDEs) ที่ได้มาจากไฟล์ point .shp ของตำแหน่งต้นไม้แต่ละแห่ง - หนึ่งแห่งจาก 1989 และอีกแห่งจากปี 2009 กราฟิกทางด้านขวาแสดงความแตกต่างระหว่าง KDE ทั้งสอง โดยที่จะแสดงเฉพาะค่า +/- 2 ส่วนเบี่ยงเบนมาตรฐานของค่าเฉลี่ย เครื่องคำนวณแรสเตอร์ของอาร์คถูกนำมาใช้เพื่อทำการคำนวณอย่างง่าย (2009 KDE - 1989 KDE) ที่จำเป็นในการสร้างภาพซ้อนแรสเตอร์บนภาพขวามือ มีวิธีการที่เหมาะสมกว่าสำหรับการวิเคราะห์ความหนาแน่นของต้นไม้หรือพื้นที่หลังคาเปลี่ยนแปลงตลอดเวลาทั้งทางสถิติหรือแบบกราฟิกหรือไม่? ด้วยข้อมูลเหล่านี้คุณจะประเมินการเปลี่ยนแปลงระหว่างข้อมูลต้นไม้ในปี 1989 และ 2009 ในสภาพแวดล้อมเชิงพื้นที่ได้อย่างไร โซลูชั่นใน ArcGIS, Python, R, Erdas และ ENVI ได้รับการสนับสนุน

15 arcgis-desktop  python  r  remote-sensing  spatial-statistics