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

ฉันมีรหัสทางภูมิศาสตร์กับคะแนนและข้อมูลที่เกี่ยวข้องบางส่วน ฉันต้องการลงจุดบนแผนที่โดยใช้ geopandas และมีขนาดของจุดที่สอดคล้องกับหนึ่งในคอลัมน์ใน geodataframe จนถึงตอนนี้ฉันมีรหัสต่อไปนี้: base = world.plot(color='white', figsize=(20,10)) geo_df.plot(ax=base, marker='.', color='red', markersize = geo_df['Pop_2005']) plt.xlim([-85, -60]) plt.ylim([-5, 12.5]); แต่ฉันได้รับข้อผิดพลาดต่อไปนี้: TypeError: cannot convert the series to <class 'float'> ความคิดใด ๆ

10 python  geopandas  matplotlib  pandas 

ฉันยังใหม่กับงูหลามและตั้งใจที่จะใช้ gdal_merge ในสคริปต์ไพ ธ อน (สายลับ) ฉันใช้ Windows 10, Python 3.6 และมีเครื่องมือ gdal ติดตั้งจาก osgeo4w ฉันรู้ว่าโพสต์อื่น ๆ อธิบายถึงปัญหานี้แต่ไม่มีใครสามารถช่วยฉันแก้ไขปัญหานี้ได้ เมื่อฉันเรียกโมดูล gdal จากคำสั่งมันทำงานเหมือนมีเสน่ห์: python "C:\OSGeo4W64\bin\gdal_merge.py" -o merged.tif input1.tif input2.tif อย่างไรก็ตามฉันไม่สามารถทำงานได้อย่างถูกต้องในสคริปต์ไพ ธ อน (Spyder) วิธีแรกสร้างผลลัพธ์ แต่ไม่ใช่ชื่อที่ถูกต้อง: สร้างไฟล์ 'out.tif' และไม่ใช่ไฟล์ 'merged.tif' ตามที่ฉันร้องขอ: import sys sys.path.append('C:\\OSGeo4W64\\bin') import gdal_merge as gm gm.main(['-o', 'merged.tif', 'N00W078.tif', 'N00W079.tif']) วิธีที่สองไม่สร้างผลลัพธ์: …

10 python  dem  windows  gdal-merge  spyder 

มีวิธีรับรายการรูปแบบ gdal ที่ติดตั้งบนระบบที่มี python หรือไม่? บนบรรทัดคำสั่งสามารถทำได้อย่างง่ายดายด้วย: gdalinfo --formats ... แต่มันสามารถทำได้ในหลาม?

10 python  raster  gdal 

ฉันพยายามที่จะสร้างรูปหลายเหลี่ยมสำหรับ swaths วงโคจรของดาวเทียม จนถึงตอนนี้ฉันมีวิธีในการสร้างสองบรรทัดซึ่งแสดงถึงขอบของแต่ละ swath ใน [lat, long] บางส่วนของแนวเขตข้อมูลข้ามประเทศและล้อมรอบ: ฉันสามารถแก้ปัญหานี้ด้วยogr2ogr -wrapdateline: ogr2ogr -wrapdateline -f "ESRI Shapefile" test.shp orbits.shp ซึ่งจะแบ่งเส้นอาจ ตอนนี้ฉันต้องการที่จะสามารถสร้างรูปหลายเหลี่ยมในการตกแต่งภายในของทั้งสองสาย ตัวอย่างเช่นในกรณีที่ขอบหนึ่งของ swath ข้ามเส้นข้อมูลรูปหลายเหลี่ยมเติมเมื่อเกิดขึ้นในอีกด้านหนึ่งเช่น: ฉันต้องการวิธีการที่เป็นไปโดยอัตโนมัติตามที่ฉันต้องการเพื่อทำงานซ้ำมาก โดยเฉพาะอย่างยิ่งในงูหลามเนื่องจากเป็นวิธีที่ฉันสร้างบรรทัด นี่คือสองไฟล์รูปร่างที่มีเส้น: wraparound ; datelinefixed

10 python  ogr  wgs84  shapely  antimeridian 

ฉันใช้ GDAL กับ python เพื่อจัดการรูปร่างของไฟล์ เมื่อเปิดไฟล์ต้นฉบับฉันจะตั้งค่าไดรเวอร์เป็น“ ESRI Shapefile” ก่อนเรียกวิธีการ "เปิด" ของไดรเวอร์ด้วยรูปร่างที่ต้องการ ฉันไม่มีปัญหากับสิ่งนี้ อย่างไรก็ตามฉันอยากรู้ว่าอะไรคือความแตกต่างระหว่างวิธีการใช้งานอินสแตนซ์ไดร์เวอร์“ เปิด” และฟังก์ชั่น ogr“ เปิด” ตัวอย่างเช่นอะไรคือความแตกต่างระหว่าง: driver = ogr.GetDriverByName("ESRI Shapefile") source = driver.Open(“O:\myshp.shp”,0) และ source = ogr.Open(“O:\myshp.shp”,0) มีสถานการณ์เฉพาะที่คุณจะใช้อีกสถานการณ์หนึ่งหรือไม่?

10 python  gdal 

ฉันสงสัยว่ามีวิธีการเปิดคอนโซล Python โดยตรงที่เปิดตัวโปรแกรมหรือไม่ มีความเป็นไปได้ในการตั้งค่าทางลัดสำหรับคอนโซล Python แต่ฉันไม่พบตัวเลือกดังกล่าวสำหรับการเปิดตัวโปรแกรม QGIS

10 qgis  python  pyqgis  qgis-2.6  qgis-python-console 

เท่านี้การเข้าถึง ArcObjects จาก Python เป็นคำถามและคำตอบที่อ่านมากที่สุดของฉันใน GIS Stack Exchange แม้ว่าความสำเร็จนั้นอาจเป็นหนึ่งในจุดอ่อนที่สุดของฉันเมื่อพูดถึงการใช้งานจริง ส่วนใหญ่ของที่แสดงยากจนมาจากความยากจนของฉันที่จะอ่านและทำความเข้าใจArcObjects เอกสาร ดังนั้นสำหรับงานที่ได้รับอะไรบ้างแนวทางสำหรับการแปล. net / c ++ / java / ... เอกสารและตัวอย่างในการเทียบเท่าของหลาม? (ภาษาใดเป็นภาษาที่ดีที่สุดในการทำงานจากเรื่องนั้น) และดัชนีที่ดีที่สุดหรือหน้า Landing Page เริ่มต้นจากอะไร สิ่งใดที่ควรเพ่งความสนใจไปและอย่างน้อยก็สำคัญเหมือนกันสิ่งที่สามารถเพิกเฉยได้อย่างอิสระ? สมมติว่าผู้ชมของคุณใช้ภาษาไพ ธ อนค่อนข้างน้อยและไม่รู้หนังสือในภาษาที่ใช้ในการพัฒนาอื่น ๆ นำเราไปสู่การฝึกเขียนโปรแกรมขนาดเล็กตั้งแต่แนวคิดเริ่มต้นและการวิจัยไปจนถึงผลลัพธ์การทำงานของหลาม

10 python  arcobjects 

เป็นไปได้ไหมที่จะตั้งค่า CRS โครงการของโครงการใน QGIS ผ่านคำสั่ง Python? ในกรณีของฉันผู้ใช้ควรเปิด QGIS เริ่มปลั๊กอินที่เขียนของฉันและปลั๊กอินนี้จะโหลดข้อมูลเชิงพื้นที่ซึ่งมีรหัส epsg 31467 (Gauss Krueger Zone 3) น่าเสียดายที่ CRS เริ่มต้นของ QGIS คือ WGS84 และฉันก็เลยมีปัญหากับยูนิตและสเกล วิธีแรกของฉันคือการใช้รหัสนี้: my_crs = core.QgsCoordinateReferenceSystem(31467, core.QgsCoordinateReferenceSystem.EpsgCrsId) self.iface.mapCanvas().mapRenderer().setDestinationCrs(my_crs) การเปลี่ยนแปลงนี้ (หลังจากข้อมูลถูกแสดงให้เห็นแล้ว) CRS เป็น Gauss Krueger Zone 3 อย่างน้อยโค้ด espg ที่มุมขวาล่างได้เปลี่ยนไป หน่วยยังคงอยู่ในระดับและขนาดผิดทั้งหมด 1: 5000 ใน Gauss Krügerคือ 1: 528822376 ใน WGS 84 …

10 qgis  python 

ฉันมีปัญหา / คำถามเกี่ยวกับหลามและ QGIS เป็นไปได้หรือไม่ที่จะได้รับข้อมูลของฟีเจอร์หนึ่งที่เลือกใน QGIS โดยใช้ python และถ้าเป็นเช่นนั้น ในตัวอย่างของฉันฉันมีเลเยอร์ซึ่งแสดงพัสดุบางส่วนสำหรับการวางผังเมือง ถ้าฉันเลือกหนึ่งในพัสดุที่ฉันต้องการรับข้อมูลสำหรับพัสดุตรงนี้ ข้อมูลนี้ควรอ่านโดยสคริปต์ไพ ธ อนจากนั้นเพิ่มลงในไฟล์ PDF ที่สร้างขึ้นโดยอัตโนมัติ ขั้นตอนสุดท้ายนั้นง่าย แต่ฉันไม่ได้รับข้อมูลจากคุณสมบัติที่เลือกนี้ เป็นไปได้ที่จะรับข้อมูลจากชั้นที่ใช้งานอยู่ แต่มันเป็นไปได้หรือไม่ที่จะได้รับข้อมูลจากคุณลักษณะที่เลือกหนึ่งของชั้นนี้ หวังว่าคุณจะช่วยฉัน ความนับถือ

10 qgis  python  pyqgis 

มีใครประสบความสำเร็จในการอัปเดตข้อมูลเมตาในโปรแกรม ArcGIS 10 หรือไม่ การพิจารณาใช้ Python / arcpy แต่ ArcObjects (C # หรือ Python / comtypes) ก็เป็นไปได้เช่นกัน ฉันต้องการอัปเดตทั้งเมตาดาต้ารูปแบบ FGDC และ ArcGIS-ISO และวิธีการใดก็ตามที่ใช้จะต้องสามารถรักษาองค์ประกอบที่มีอยู่ (ไม่ว่างเปล่า) พร้อมกับองค์ประกอบที่เพิ่มเข้ามาได้ยกเว้นในกรณีที่ขัดแย้งกัน องค์ประกอบที่เพิ่มเข้ามาเขียนทับองค์ประกอบที่มีอยู่

10 arcobjects  arcgis-10.0  python  c#  metadata 

ฉันกำลังพยายามหาปริมาณของไอน้ำที่ตกตะกอน (PWV), โอโซนและละอองลอยเป็นหน้าที่ของเวลาเหนือจุดเฉพาะบนโลกนั่นคือหอสังเกตการณ์ทางดาราศาสตร์ของเรา ในการทำเช่นนี้ฉันได้รับรหัส Python แล้วmodapsclientซึ่งจะดาวน์โหลดผลิตภัณฑ์ MODIS Aqua และ Terra MYDATML2 และ MODATML2 สองครั้งต่อวันซึ่งครอบคลุมละติจูดและลองจิจูดเฉพาะที่ฉันสนใจ สิ่งที่ฉันไม่แน่ใจเกี่ยวกับคือวิธีการดึงปริมาณเฉพาะที่ฉันต้องการเช่นเวลาที่ใช้ข้อมูล MODIS และ PWV สำหรับตำแหน่งละติจูดและลองจิจูดเฉพาะของหอดูดาวของฉันเพื่อทำให้เป็นค่าอนุกรมเวลา ผลิตภัณฑ์ MYDATML2 ดูเหมือนจะมี 2D ละติจูดและลองจิจูดกริดของCell_Along_Swath_5kmและCell_Across_Swath_5kmดังนั้นฉันเดานี้ทำให้แนวข้อมูลเมื่อเทียบกับข้อมูลตารางกระเบื้องหรือ? ปริมาณที่ฉันต้องการเช่นนั้นPrecipitable_Water_Infrared_ClearSkyดูเหมือนว่าจะขัดกับCell_Along_Swath_5kmและCell_Across_Swath_5kmฉันไม่แน่ใจว่าจะรับค่า PWV ที่ lat ที่เฉพาะเจาะจงได้นานแค่ไหนที่ฉันสนใจโปรดช่วยด้วย?

9 python  remote-sensing  modis 

ฉันต้องดาวน์โหลดภาพ Landsat โดยอัตโนมัติเพื่อที่จะทำป่าชุมชน ฉันกำลังใช้แพ็กเกจ landat-utils อย่างไรก็ตามเมื่อฉันตั้งฉากของฉันจากภูมิภาคความสนใจของฉันรหัสไม่ทำงานตามที่ฉันคาด ดังนั้นนี่คือตัวอย่างหนึ่งของ sceneID (LC82150682015350LGN01) ซึ่งฉันต้องการได้รับเพื่อทดสอบสคริปต์ของฉัน ฉันพยายามใส่ชื่อผู้ใช้และรหัสผ่านจาก USGS แต่ไม่สามารถแก้ไขปัญหาได้ แต่ปัญหายังคงอยู่ ที่นี่รหัส: $ landsat download LC82150682015350LGN01 AttributeError: 'RemoteFileDoesntExist' object has no attribute 'message' เมื่อฉันตั้งชื่อผู้ใช้และรหัสผ่าน USGS ของฉันข้อผิดพลาดได้: AttributeError: 'USGSError' object has no attribute 'message' เหตุใดฉันจึงไม่ได้รับไฟล์พร้อมฉากนี้ ID ฉันเปลี่ยนมันและฉันมีปัญหาเดียวกัน

9 python  landsat  download  usgs 

ฉันมีรูปหลายเหลี่ยมในไฟล์ shp ที่ฉันต้องเปลี่ยนความกว้างของเส้น ฉันหัวเข้าร่วมพวกเขาในกลุ่มเดียว มีวิธีในการเปลี่ยนบรรทัดด้วยสำหรับเลเยอร์ทั้งหมดในกลุ่มพร้อมกันหรือไม่

9 qgis  python  pyqgis  layers  style 

ฉันจะบันทึกฟังก์ชัน Python ที่กำหนดเองของฉันใน QGIS ได้อย่างไร ฉันต้องการให้ผู้ใช้ช่วยเหลือสำหรับฟังก์ชั่นที่ฉันเขียนให้พวกเขา: ฉันได้ลองเพิ่ม docstring แล้ว: แต่เมื่อใดก็ตามที่ฉันคลิกโหลดข้อผิดพลาดจะปรากฏขึ้น: ข้อความแสดงข้อผิดพลาดแสดงข้อผิดพลาดการเยื้องอย่างชัดเจนซึ่งไม่สอดคล้องกับรหัสจริงของฉัน

9 qgis  python  field-calculator  qgis-custom-function 

ฉันมีไฟล์ 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