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

ฉันพยายามสร้างปลั๊กอินเพื่อโหลดผู้แต่งเพลงจากไฟล์สร้างแอตลาสและส่งออกไปยังรูปภาพ จนถึงตอนนี้ฉันประสบความสำเร็จในการโหลดเทมเพลตและส่งออกไปยังรูปภาพ ฉันไม่สามารถเพิ่มเลเยอร์ใด ๆ ในตำนาน (ซึ่งอยู่ใน toc) ลงในแผนที่ที่ส่งออกซึ่งส่งผลให้แผนที่ว่างเปล่าและไม่มีการแสดงออกของฟิลด์ใดทำงาน # Get layers in the legend and append, must be a cleaner way to do this? layers = self.iface.legendInterface().layers() layerStringList = [] for layer in layers: layerID = layer.id() layerStringList.append(layerID) # Add layer to map render myMapRenderer = QgsMapRenderer() myMapRenderer.setLayerSet(layerStringList) myMapRenderer.setProjectionsEnabled(False) # …

10 python  qgis-plugins  pyqgis  qgis-2  print-composer 

เป็นไปได้ไหมที่จะเปลี่ยนสีในขณะที่เพิ่มเลเยอร์เวกเตอร์ ฉันมีเวกเตอร์เลเยอร์ที่มีพื้นที่น้ำและต้องการโหลดด้วยสไตล์สีฟ้า เพิ่มเลเยอร์ด้วย: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) ฉันสามารถเปลี่ยนสีหลังจากหรือขณะโหลดได้หรือไม่?

10 qgis  python  pyqgis  color 

ฉันพบข้อมูลการยกระดับความแม่นยำสูง (1/3 - 1/9 ความละเอียดอาร์ซีวินาที) จากชุดข้อมูลความสูงแห่งชาติที่จัดทำโดย USGS มันมาพร้อมรูปแบบ IMG, Arcgrid และ GridFloat ฉันรู้ว่า ArcGrid นั้นสอดคล้องกับแพคเกจซอฟต์แวร์แบบชำระเงิน แต่ฉันพยายามที่จะใช้เครื่องมือที่มีให้อย่างอิสระ ฉันมีข้อมูล GPS ที่ฉันพยายามจะสัมพันธ์กับระดับพื้นดิน มีห้องสมุดไพ ธ อนใดบ้างที่จะให้ฉันแปลงข้อมูล gps เป็นระดับความสูงโดยใช้ดาต้าไทล์ในรูปแบบ IMG, ArcGrid หรือ GridFloat

10 python  elevation  esri-grid-format 

ฉันต้องการสร้างฟังก์ชั่นที่เลือกคุณสมบัติและซูมเข้าไป (คล้ายใน QGIS) ดังนั้นจึงมีฟังก์ชั่นดังต่อไปนี้: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) def zoomTo(self): layer = self.vlayer atable = self.ui.table selectList=[] for i in atable.selectionModel().selectedRows(): ID = atable.item(i.row(),0).text() selectList.append(int(ID)) layer.setSelectedFeatures(selectList) คุณสมบัติที่เลือกจะถูกเน้นบนแผนที่ แต่ฉันไม่รู้ว่าจะทำ "ซูม" กับคุณสมบัติที่เลือกหรือโฟกัสบางอย่างในแผนที่ได้อย่างไร

10 qgis  python  pyqgis 

ฉันใช้สาย folowing เพื่อแจ้งผู้ใช้เกี่ยวกับสถานะ: iface.mainWindow().statusBar().showMessage("Status:" + str(i)) ปลั๊กอินใช้เวลาประมาณ 2 นาทีในการรันบนชุดข้อมูลของฉัน แต่ windows ตรวจพบว่าเป็น "ไม่ตอบสนอง" และหยุดแสดงการอัปเดตสถานะ สำหรับผู้ใช้ใหม่มันไม่ค่อยดีนักเนื่องจากดูเหมือนว่าโปรแกรมจะพัง มีวิธีแก้ไขหรือไม่เพื่อให้ผู้ใช้ไม่ถูกทิ้งในที่มืดเกี่ยวกับสถานะของปลั๊กอิน?

10 qgis  python 

เป็นครั้งคราว แต่น่าเสียดายที่บ่อยครั้งบริการ AGS ของฉันก็หยุดโดยไม่มีเหตุผล ฉันมักจะไม่สังเกตเห็นสิ่งนี้จนกว่าผู้ใช้จะโทรหาฉันและถามว่าทำไมแผนที่ของพวกเขาจึงไม่มีข้อมูล มีวิธีใดบ้าง (ฉันกำลังคิด Python) ที่ฉันสามารถรับการตอบสนองจาก AGS เมื่อบริการหยุดทำงาน ผู้ใช้ AGS รายอื่นจัดการกับปัญหานี้อย่างไร

10 python  arcgis-server  administration 

ปิด. คำถามนี้เป็นคำถามปิดหัวข้อ ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้เป็นไปตามหัวข้อสำหรับระบบแลกเปลี่ยนข้อมูลทางภูมิศาสตร์ของกองซ้อน ปิดเมื่อเดือนที่แล้ว ฉันมีกลุ่มสตรีมที่ยาว 1,000 กม. ฉันต้องการค้นหาความแตกต่างของระดับความสูงระหว่างจุดสองจุดติดต่อกันระยะทาง 1 กม. เริ่มต้นจากต้นน้ำถึงปลายน้ำ ฉันจะรับความแตกต่างจากระดับความสูงจาก DEM ได้อย่างไร ฉันมีส่วนของสตรีมในรูปแบบแรสเตอร์และในรูปแบบเวกเตอร์ มันจะดีกว่านี้ถ้าฉันมีความคิดเกี่ยวกับสคริปต์ Python

10 python  dem  elevation 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน3 ปีที่ผ่านมา ไม่มีใครมีตัวอย่างที่ดีของการใช้ Python ในการทำภารกิจสำรวจโดยอัตโนมัติหรือทำให้งานง่ายขึ้น? ฉันกำลังสอนแล็บ Python / ArcPy ในสัปดาห์นี้เพื่อเรียน GIS มีการสำรวจนักเรียนจำนวนมากในชั้นเรียนที่คิดว่า GIS เป็นส่วนใหญ่เสียเวลาของพวกเขาและใช้เวลาเรียนเพียงเพราะมันจำเป็น ฉันต้องการให้สิ่งที่พวกเขาตื่นเต้น บริบทเล็ก ๆ น้อย ๆ : นักเรียนส่วนใหญ่อาจจะทำงานให้กับ บริษัท ป่าไม้หรือ บริษัท ชื่อเมื่อพวกเขาสำเร็จการศึกษาและเราอยู่ในอาร์คันซอดังนั้นเกือบทุกอย่างที่นี่อยู่บนพื้นฐานของ PLSS จนถึงตอนนี้ความคิดของฉันแสดงให้เห็นถึงฟังก์ชันคณิตศาสตร์เชิงมุมและกล่าวถึงโปรแกรม CAD ที่ใช้ Python ฟรี ฉันยังได้เรียนรู้ว่ามีโครงการที่จะเขียนปลั๊กอินสำหรับ AutoCAD ที่ให้อินเทอร์เฟซ Python ฉันได้สอนชั้นเรียนแล้วและตามที่คาดไว้ผู้สำรวจไม่ได้ตื่นเต้นเกินไป ฉันยังชอบที่จะได้ยินตัวอย่างที่ดีของการใช้ Python เพื่ออำนวยความสะดวกในการสำรวจ

10 python  arcpy  education  land-survey  teaching 

ฉันพยายามเรียนรู้วิธีใช้ 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 

ฉันมีไฟล์รูปร่างเป็นจำนวนมากที่แสดงถึงพื้นที่ที่น่าสนใจสำหรับการวิเคราะห์ที่จะดำเนินการโดยใช้แหล่งที่มาของภาพถ่ายดาวเทียม (IKONOS, RapidEye ฯลฯ ) น่าเสียดายที่ภาพไม่ได้ใช้ระบบ Pathrow เช่น Landsat เป็นต้นขอบเขตจึงแตกต่างกันอย่างมาก ฉันมีรูปร่างที่ตัดกับ AOI แต่ละตัวซึ่งแสดงถึงขอบเขตของการได้รับภาพที่แตกต่างกันซึ่งทั้งหมดนี้ได้รับการยอมรับแล้ว บางรูปร่างเหล่านี้มีรูปหลายเหลี่ยมมากกว่า 500 รูป ฉันต้องการค้นหาวิธีที่ควรเป็นแบบอัตโนมัติ (Python และ ArcInfo 10 โดยเฉพาะอย่างยิ่ง FOSS ก็เป็นที่ยอมรับเช่นกัน) เพื่อกำหนดรูปหลายเหลี่ยมจำนวนน้อยที่สุดเพื่อครอบคลุมพื้นที่ที่ฉันสนใจแต่ละแห่ง

10 arcgis-10.0  python  shapefile  geoprocessing 

ฉันพยายามที่จะเข้าใจว่าฉันสามารถใช้โมดูล csv ใน python เพื่อเปิดไฟล์ csv ในโฟลเดอร์เดียวกับสคริปต์ python ได้อย่างไรจากนั้นสร้าง shapefile โดยใช้ shapefile module pyshp ไฟล์ csv มีลักษณะเช่นนี้ แต่สามารถมีเรคคอร์ดสองสามพันแถว: id_nr;date;target;start_lat;start_lon 1;2012-05-21;navpoint 25x;55.123654;13.456954 1;2012-05-23;navpoint 11f;55.143654;12.456954

10 python  shapefile  csv  pyshp 

ฉันมีเครื่องมือสคริปต์ที่ใช้งานง่ายที่วนรอบพื้นที่ทำงานและเปลี่ยนชื่อและคัดลอกรูปร่างไฟล์ไปยังชุดข้อมูลคุณสมบัติ อย่างไรก็ตามหากมีไฟล์รูปร่างที่เสียหายอยู่ในพื้นที่ทำงานสคริปต์จะล้มเหลวและหยุดการประมวลผล คุณจัดการข้อผิดพลาดเช่นนี้ได้อย่างไร มีวิธีพิมพ์ไฟล์ข้อผิดพลาดและดำเนินการประมวลผล shapefile ถัดไปใน for for loop to complete หรือไม่ import arcpy from arcpy import env # Allow overwriting of output env.overwriteOutput = True # Parameters env.workspace = arcpy.GetParameterAsText(0) state = arcpy.GetParameterAsText(1) gdb = arcpy.GetParameterAsText(2) # Get a list of shapefiles in folder fcs = arcpy.ListFeatureClasses() # Find …

10 python  arcpy 

ฉันพยายามทำให้ร่างแสดงข้อมูล azimuthal ด้วยช่วงความไม่แน่นอนที่แตกต่างกันในแต่ละจุด รูปเก่าแก่ของโรงเรียนจากกระดาษ 1991 จับแนวคิด "bowtie plot" ที่ฉันตั้งเป้าไว้: ข้อเสนอแนะใด ๆ เกี่ยวกับวิธีที่ฉันจะทำให้ร่างที่คล้ายกัน? ฉันเป็นมือใหม่ที่เกี่ยวข้องกับ GIS แต่ฉันสามารถเข้าถึง ArcGIS ผ่านมหาวิทยาลัยของฉันได้ ประสบการณ์ Arc ของฉันถูก จำกัด ในการสร้างแผนที่ทางธรณีวิทยาดังนั้นฉันจึงไม่ต้องทำอะไรแปลกใหม่เกินไป ฉันแหย่ไปในตัวเลือกสัญลักษณ์ใน Arc และ QGIS แต่ไม่เห็นการตั้งค่าใด ๆ ที่ฉันคิดว่าจะทำงาน โปรดทราบว่านี่ไม่ได้เป็นเพียงเรื่องของการหมุนสัญลักษณ์รูปโบว์โดย azimuth; ช่วงเชิงมุมของ "bowtie" แต่ละอันจะต้องแตกต่างกัน ฉันให้คะแนนทักษะ Python ของฉันเป็น 'strong ตัวกลาง' และทักษะ R ของฉันเป็น 'low middle' ดังนั้นฉันไม่รังเกียจที่จะแฮ็คข้อมูลบางอย่างร่วมกับmatplotlibและmpl_toolkits.basemapหรือไลบรารีที่คล้ายคลึงกันหากจำเป็น แต่ฉันคิดว่าฉันขอคำแนะนำที่นี่ก่อนที่จะลงไปที่ถนนในกรณีที่มีวิธีแก้ปัญหาที่ง่ายกว่าจาก GIS-land ที่ฉันไม่เคยได้ยินมาก่อน

10 python  algorithm  r  visualisation  spherical-geometry 

ในขณะที่วนซ้ำเวกเตอร์เลเยอร์โดยใช้รหัสต่อไปนี้ (ทำให้ร้อนขึ้นจากตัวอย่างในpyqgis cookbook ) มีวิธีตรวจสอบว่ามีการเลือกสถานที่หรือไม่? provider = vlayer.dataProvider() feat = QgsFeature() allAttrs = provider.attributeIndexes() provider.select(allAttrs) while provider.nextFeature(feat): geom = feat.geometry() attrs = feat.attributeMap() for (k,attr) in attrs.iteritems(): print "%d: %s" % (k, attr.toString()) อีกทางหนึ่งฉันสามารถสร้างรายการคุณลักษณะที่เลือกโดยใช้vlayer.selectedFeatures()แต่ฉันหวังว่าจะมีวิธีการตรวจสอบแต่ละคุณสมบัติโดยตรง

10 qgis  python  pyqgis  select 

เป็นไปได้หรือไม่ที่จะจับข้อความแสดงข้อผิดพลาดเอาต์พุต 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