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

องค์กรของเรากำลังพิจารณาที่จะย้ายเวิร์กโฟลว์การประมวลผลทางภูมิศาสตร์ไปยัง PostGIS ขณะนี้เรากำลังใช้ ArcGIS ด้วยเครื่องมือ Python ที่กำหนดเองมากมายที่ใช้ใน ModelBuilder เรากำลังย้ายข้อมูลส่วนใหญ่ของเราไปยัง PostGIS เพื่อการใช้งานแอพที่หลากหลายและตอนนี้เรากำลังถามว่ามันเหมาะสมที่จะทำการประมวลผลข้อมูลที่นั่นหรือไม่ เราประมวลผลข้อมูลเพื่อให้เข้ากันได้กับซอฟต์แวร์ของเรา ลูกค้าซื้อซอฟต์แวร์ของเราให้ข้อมูลกับเราและเราดำเนินการเพื่อให้ได้รับประโยชน์สูงสุดสำหรับการใช้งานในซอฟต์แวร์ของเรา สิ่งนี้ต้องการให้เราสร้างเครื่องมือที่หลากหลายเพื่อจัดการคุณภาพการป้อนข้อมูลที่หลากหลาย เราไม่สามารถคาดหวังว่าจะได้รับข้อมูลในรูปแบบเฉพาะหรือสคีมาดังนั้นเราจึงสร้างเครื่องมือในการแมปฟิลด์อินพุตไปยังฟิลด์เอาต์พุตแยกวิเคราะห์ฟิลด์เดียวเป็นหลายฟิลด์ผสานหลายชุดข้อมูล ฯลฯ นอกจากนี้เรายังทำการรวมเชิงพื้นที่ และเชื่อมต่อฟิลด์และการดำเนินการทั่วไปอื่น ๆ อีกมากมาย PostGIS นั้นมีความสามารถในการดำเนินการตามความต้องการในการประมวลผลทั้งหมดของเราอย่างสมบูรณ์แบบ สำหรับผู้ที่ใช้ PostGIS เพื่อทำการประมวลผลข้อมูลคุณมีคำแนะนำสำหรับองค์กรเครื่องมือในการใช้ ฯลฯ หรือไม่? คุณใช้ร่วมกับการประมวลผลหลาม QGIS หรือไม่? คนใช้ Python ORM สำหรับการประมวลผลที่ไม่ใช่เว็บ ฉันโน้มตัวไปใช้GeoDjangoเพราะมันมี Python ORM สำหรับ PostGIS การทดสอบเริ่มต้นของเราในการใช้ PostGIS เพื่อประมวลผลข้อมูลมีบล็อคข้อความ SQL ขนาดใหญ่จำนวนมากในรหัส Python และเราคิดว่า GeoDjango ORM อาจช่วยในการสร้างโค้ดที่จัดการและอ่านได้ง่ายขึ้น นอกจากนี้ยังมีGeoAlchemy …

12 python  postgis  geoprocessing  geodjango  geoalchemy 

ฉันกำลังเขียนรหัสหลามเพื่อแปลง X, Y คะแนนใน excel เป็น shapefile ในกระบวนการนี้ฉันมีคะแนน tyad xy จาก shhet1 ขั้นตอนที่ฉันต้องรวมในกระบวนการของฉันคือการอ่านข้อมูลจากแผ่นงาน 1 ของสมุดงาน excel (97-2003) ฉันเขียนโค้ดดังนี้ .. import arcpy from arcpy import env import Xlrd env.workspace="E:\" input_table="123.xls\Sheet1" แต่มันไม่ทำงาน

12 python  excel 

ฉันถูกขอให้ผู้จัดการจัดทำมาตรการวัดระดับความสามารถเพื่อให้พนักงานเข้าใจการเขียนสคริปต์ / การเขียนโปรแกรม เนื่องจากฉันเป็นโปรแกรมเมอร์ที่สอนตัวเองฉันจึงไม่มีความเข้าใจอย่างเป็นทางการเกี่ยวกับวิธีการวัดระดับความเข้าใจเมื่อพูดถึงการเขียนบท เราใช้ Python ใน ArcGIS และ QGIS แล้ว เราใช้ VB.Net กับผลิตภัณฑ์ MS Office ของเราและเราได้เริ่มตะลุยใน Java สำหรับการทำแผนที่เว็บ ฉันสงสัยว่าผู้จัดการ GIS คนใดในพวกคุณมีเมทริกซ์ / วิธีการบางอย่างเพื่อกำหนดระดับความสามารถของพนักงานของคุณเมื่อพูดถึงสคริปต์

12 python 

ฉันประหลาดใจจริง ๆ ว่าขั้นตอนใหญ่ ๆ ที่การพัฒนา QGIS เกิดขึ้นในช่วงหลายเดือนที่ผ่านมา ชุดไอคอนใหม่ตัวแต่งการพิมพ์ที่ได้รับการปรับปรุงและการเปลี่ยนแปลงอื่น ๆ อีกมากมายใน QGIS 2.0 นั้นยอดเยี่ยมมาก! ท่ามกลางการปรับปรุงที่ยอดเยี่ยมอื่น ๆ อีกมากมายนอกจากนี้ยังมีแถบงานสีน้ำเงินที่ดีซึ่งฉันพบว่าน่าดึงดูดมาก เป็นไปได้หรือไม่ที่จะกำหนดแถบนี้เป็นไพ ธ อนสำหรับสคริปต์และปลั๊กอิน (เช่นเพื่อแสดงข้อความเล็ก ๆ หลังจากการคำนวณเสร็จสิ้น)? ฉันถามคำถามนี้ที่นี่เพราะฉันเชื่อว่านักพัฒนา QGIS จำนวนมากยังเยี่ยมชมเว็บไซต์ถามตอบนี้ PS: แถบจะปรากฏเฉพาะใน QGIS dev และสำหรับการดำเนินการบางอย่างเท่านั้น

12 qgis  python  pyqgis 

ฉันมีArc / Info Binary Grid --- โดยเฉพาะแรสเตอร์การไหลของ ArcGIS --- และฉันต้องการระบุเซลล์ทั้งหมดที่มีค่าเฉพาะ (หรืออยู่ในช่วงของค่า) ในที่สุดฉันต้องการไฟล์รูปร่างของคะแนนที่แสดงถึงเซลล์เหล่านี้ ฉันสามารถใช้ QGIS เพื่อเปิด hdr.adf และรับผลลัพธ์นี้เวิร์กโฟลว์คือ: QGIS> เมนู Raster> Raster Calculator (ทำเครื่องหมายคะแนนทั้งหมดด้วยมูลค่าเป้าหมาย) QGIS> เมนู Raster> รูปหลายเหลี่ยม QGIS> เมนูเวกเตอร์> เมนูย่อยเรขาคณิต> รูปหลายเหลี่ยม centroids แก้ไขเซนทรอยด์เพื่อลบโพลีเซนทรอยด์ที่ไม่ต้องการ (เหล่านั้น = 0) วิธีการนี้ "ทำงานได้" แต่มันไม่ได้ดึงดูดความสนใจฉันเพราะมันสร้างไฟล์ 2 ไฟล์ที่ฉันต้องลบแล้วฉันต้องลบบันทึกที่ไม่ต้องการออกจาก shapefile ของ centroids (นั่นคือ = 0) คำถามที่มีอยู่วิธีการเรื่องนี้ แต่มันเหมาะสำหรับ ArcGIS …

12 python  raster  gdal 

เมื่อพยายามติดตั้ง Shapely บนคอมพิวเตอร์ Windows 64 บิตของฉันฉันไม่สามารถทำให้ไลบรารี GEOS ทำงานได้ จนถึงตอนนี้ฉันได้เรียกใช้ตัวติดตั้ง OSGeo4W ซึ่งฉันติดตั้ง GDAL (ฉันเชื่อว่ามีไลบรารีทางภูมิศาสตร์รวมอยู่ในแพ็คเกจนั้น) หลังจากนั้นฉันตรวจสอบและฉันมี geos_c.dll ในC:\OSGeo4W\binไดเรกทอรีของฉันแต่ฉันไม่ได้ทำตามขั้นตอนการกำหนดค่าบางอย่างหรือไลบรารีไม่ทำงาน ฉันต้องการหุ่นดีในการทำงานดังนั้นฉันจึงวิ่งpip install shapelyหลังจากติดตั้ง GDAL และดูเหมือนว่าจะใช้งานได้ (แม้ว่าจะไม่สามารถหาไลบรารี C สำหรับ GEOS) ในรหัสของฉันฉันสามารถนำเข้า Shapely แต่เมื่อฉันพยายามที่จะใช้มันฉันได้รับข้อผิดพลาดบอกฉัน"geos.dll"ไม่พบ

12 qgis  python  gdal  shapely 

ฉันพยายามเบิร์นไฟล์ไปยังแรสเตอร์โดยใช้ RasterizeLayer ของ GDAL ฉันสร้างพื้นที่แรสเตอร์ที่น่าสนใจจากรูปร่างที่แตกต่างกันก่อนกำหนดขนาดพิกเซลที่เฉพาะเจาะจง AOI นี้ทำหน้าที่เป็นฐานสำหรับการแรสเตอร์ต่อไปนี้ทั้งหมด (จำนวน collumns และแถวเดียวกัน, การฉายภาพและ geotransform เดียวกัน) ปัญหาเกิดขึ้น แต่เมื่อฉันไปเผารูปร่างเพื่อแรสเตอร์ของตัวเองตามขนาดพิกเซลและประมาณการเดียวกัน ลิงค์ด้านล่าง (มีตัวแทนไม่เพียงพอที่จะโพสต์ภาพ) แสดงไฟล์รูปร่างดั้งเดิมในผิวสีแทนและสีชมพูเข้มที่ RasterizeLayer เขียนข้อมูล สีชมพูอ่อนเป็นค่า nodata สำหรับข้อมูลแรสเตอร์สีชมพูเข้ม สีเทาคือ AOI ตามการเบิร์นไฟล์รูปร่างเสร็จสมบูรณ์ เมื่อพิจารณาถึงขอบเขตของรูปหลายเหลี่ยม shapefile ฉันคาดว่าจะเห็นค่าการเบิร์นที่มุมด้านล่างทั้งสองรวมทั้งพิกเซลสองพิกเซลที่อยู่ใต้ข้อมูลที่แสดง เห็นได้ชัดว่านี่ไม่ใช่กรณี ดังต่อไปนี้เป็นรหัสที่ฉันใช้ในการสร้างสิ่งเหล่านี้ รูปร่างทั้งหมดถูกสร้างขึ้นโดยใช้ QGIS และถูกสร้างขึ้นทั้งหมดในการฉายภาพเดียวกัน (ควรสังเกตว่าการ gridding ในภาพที่แสดงเป็นเพียงเพื่อให้ความคิดเกี่ยวกับขนาดพิกเซลที่ฉันใช้) from osgeo import ogr from osgeo import gdal aoi_uri = 'AOI_Raster.tif' aoi_raster = …

12 python  raster  gdal  rasterization 

ฉันทำงานกับ GDAL ใน Python และต้องการให้ GDAL เพิ่มข้อยกเว้นเพื่อให้ฉันสามารถทำสิ่งที่มีประโยชน์กับพวกเขาได้ เช่นหากไฟล์ไม่สามารถเปิดด้วย gdal.Open () ฉันต้องการที่จะโยนข้อยกเว้นแทน gdal กลับไม่มี ฉันต้องการวิธีจับข้อผิดพลาดแทนที่จะถูกทิ้งไปยัง sys.stdout ขณะนี้ฉันได้ลองใช้gdal.UseExceptions()ฟังก์ชั่น แต่พบปัญหา เมื่อgdal.UseExceptions()มีการตั้งค่าดูเหมือนว่าจะซ่อนข้อผิดพลาดบางอย่างและไม่ยกข้อยกเว้นพวกเขา ตัวอย่างเช่นERROR 5: Access window out of range in RasterIO() Requested (1,15) of size 25x3 on raster of 26x17หายไปและไม่มีข้อยกเว้นยก มันดูอันตรายทีเดียว ขอบคุณสำหรับความคิดใด ๆ

12 python  gdal  ogr  osgeo 

ฉันสนใจที่จะเรียนรู้วิธีการใช้ประโยชน์จากพลังการประมวลผลแบบมัลติคอร์ที่มีอยู่ในคอมพิวเตอร์เดสก์ท็อป อาร์คระบุว่าพื้นหลังการประมวลผลทางภูมิศาสตร์ช่วยให้ผู้ใช้สามารถใช้หลายคอร์อย่างไรก็ตามงานจำเป็นต้องรอให้งานก่อนหน้าเสร็จสิ้น มีใครพัฒนาวิธีการประมวลผลทางภูมิศาสตร์แบบขนานหรือแบบมัลติเธรดใน Arc / Python หรือไม่? มีปัญหาคอขวดของฮาร์ดแวร์ที่ป้องกันการประมวลผลแบบมัลติคอร์ในแต่ละงานหรือไม่? ฉันพบตัวอย่างที่น่าสนใจใน Stackoverflow ที่ดึงดูดความสนใจของฉันแม้ว่าจะไม่ใช่ตัวอย่างของกระบวนการทางภูมิศาสตร์: from multiprocessing import Pool import numpy numToFactor = 976 def isFactor(x): result = None div = (numToFactor / x) if div*x == numToFactor: result = (x,div) return result if __name__ == '__main__': pool = Pool(processes=4) possibleFactors = range(1,int(numpy.floor(numpy.sqrt(numToFactor)))+1) print …

12 python  arcpy  parallel-processing  big-data 

รับการเรียนรู้กระบวนการทางภูมิศาสตร์ด้วย Python โดยใช้หลักสูตรGEOG 825 (การประมวลผลทางภูมิศาสตร์ของหลาม) ของ MSU น่าเสียดายที่ลิงค์เสียชีวิตแล้ว ฉันรู้ว่านี่เป็นช็อตที่ยาวนาน แต่ไม่มีใครมีสำเนาของเนื้อหาจากหลักสูตรนั้นหรือไม่? ฉันอยากจะจบมันจริงๆ ตรวจสอบหลักสูตรการประมวลผล Python ทางออนไลน์อื่นแล้ว หนึ่งที่พบในรัฐเพนน์และในรัฐยูทาห์ แม้ว่าผู้ใช้จะเน้นที่การใช้เครื่องมือไม่ว่าจะเป็น ogr / gdal หรือ ArcGIS 'arcpy / arcgisscripting มีเพียงหลักสูตรของ Michigan State U เท่านั้นที่พูดถึงอัลกอริธึมและโครงสร้างข้อมูลเท่าที่ฉันรู้

12 python  geoprocessing  education 

เมื่อไม่นานมานี้ฉันได้ทำการสำรวจการวางแนวหลังคาบ้านใกล้เรือนเคียงสำหรับการติดตั้งแผงโซลาร์เซลล์โดยการตรวจสอบแต่ละหลังคาและทำเครื่องหมายทีละภาพไม่ว่าจะเป็นทิศใต้หรือเกือบหันไปทางทิศใต้ เนื่องจากใช้เวลานานมากฉันจึงต้องการหาวิธีทำให้เป็นอัตโนมัติ แต่ไม่รู้ว่าจะเป็นไปได้อย่างไรโดยการวิเคราะห์เงาหลังคา ฉันสามารถเขียนสคริปต์เล็กน้อย แต่ไม่เคยเขียนปลั๊กอินใด ๆ ฉันต้องรวมองค์ประกอบอื่นใดอีก เราอยู่ในงบประมาณที่ต่ำมากดังนั้นขอคำแนะนำใด ๆ โดยใช้เครื่องมือโอเพนซอร์ส

12 python  open-source-gis  satellite 

ฉันมีชุดข้อมูลแรสเตอร์ในรูปแบบ ASCII ใช้ Python ฉันจำเป็นต้องกำหนดminและmaxค่าภายในชุดข้อมูล ฉันได้รับแจ้งว่าข้อมูลส่วนหัวเป็นกุญแจซึ่งมีหลายอย่างเช่นจำนวนแถว / คอลัมน์ขนาดของเซลล์และอื่น ๆ คุณไม่สามารถข้ามข้อมูลส่วนหัวและอ่านชุดข้อมูลทั้งหมดเพื่อกำหนดminและmaxค่าได้หรือไม่ นี่คือสิ่งที่ฉันพยายามจะทำ ฉันข้ามคู่แรกของบรรทัดที่มีข้อมูลส่วนหัวและพยายามกำหนดค่าจากที่นั่น ต่อไปนี้เป็นสิ่งที่ฉันมี แต่ต้องการคำแนะนำบางอย่างเนื่องจากฉันยังใหม่กับ Python raster_file = open('data.asc', 'r') # Open the file data = raster_file.readlines()[4:] # Read the lines in the file, and skip the first six lines for lines in data: print max(data) # Find the max value …

12 python  statistics  esri-ascii-raster 

ฉันต้องใช้การดำเนินการทางเรขาคณิตที่ง่ายมาก (อ่าน SHP หรือรูปแบบมาตรฐานอื่น ๆ และการค้นหาง่าย ๆ ) ภายในอุปกรณ์ฝังตัวที่ใช้ระบบปฏิบัติการ Linux ที่ใช้งานง่าย API ที่เล็กที่สุดและง่ายที่สุดที่ฉันสามารถใช้คืออะไร มันสามารถเขียนใน C / C ++ หรือ Python เพื่ออนุญาตให้รวมกับระบบของฉัน ขอบคุณ ซามูเอล

12 python  linux 

ฉันได้รับการเข้ารหัสกับ Python มาหลายเดือนแล้วและได้พัฒนาสคริปต์ที่ซับซ้อนพอสมควรสำหรับงานการประมวลผลทางภูมิศาสตร์เป็นหลัก ที่ถูกกล่าวว่าฉันยังคงเรียนรู้มากขณะที่ฉันมาจากพื้นหลัง SQL / VBA / VBScript ฉันรู้ว่ารหัสที่คอมไพล์ทำงานได้เร็วกว่าโค้ดที่ต้องใช้ตัวแปลภาษาดังนั้นฉันจึงสนใจที่จะรวบรวมสคริปต์ Python สำหรับการประมวลผลทางภูมิศาสตร์กับไฟล์. EXE สำหรับการทำงานกับข้อมูลขนาดใหญ่ เป็นไปได้ไหม ถ้าเป็นวิธีที่ดีที่สุดในการรวบรวมสคริปต์ Python (.py) ที่นำเข้าโมดูล arcgisscripting หรือ arcpy คืออะไร ฉันใช้เวลาสองสามนาทีเพื่อค้นหาสิ่งที่ฉันต้องการจะทำและการค้นหาส่งคืนบทความนี้ให้คนอื่น ๆ : http://www.ehow.com/how_2091641_compile-python-code.html คอมไพเลอร์ดูเหมือนจะใช้งานได้ แต่เมื่อทำการประมวลผลไฟล์. EXE ที่เกิดขึ้นมันทำให้เกิดข้อผิดพลาดที่ซ่อนเร้นว่าบางไฟล์ไม่มี สคริปต์ Python ทำงานในสิ่งที่ดูเหมือนจะดีพอสมควรจากบรรทัดคำสั่ง แต่ฉันสงสัยว่าฉันจะเห็นการปรับปรุงเล็กน้อยถ้าฉันสามารถรวบรวมไฟล์. py อีกครั้งฉันกำลังทำงานกับชุดข้อมูลขนาดใหญ่ที่ใช้เวลาในการประมวลผลประมาณ 20 ชั่วโมง (การแยกแหล่งต้นกำเนิดจากแหล่งตัวอย่างคุณภาพน้ำที่ป้อน) ฉันจะทำทุกอย่างที่จะทำให้ดีขึ้น สคริปต์รันเร็วกว่า 10% นอก ArcGIS จากบรรทัดคำสั่งโดยใช้ชุดทดสอบของเว็บไซต์เทียบกับการตั้งค่าสคริปต์เป็นเครื่องมือสคริปต์ในกล่องเครื่องมือใหม่ใน ArcCatalog ฉันใช้งานสคริปต์จากบรรทัดคำสั่งโดยไม่มีอินสแตนซ์ใด ๆ ของ …

12 arcpy  python  performance  compile  cython 

ฉันมีแรสเตอร์ที่ฉันต้องการจะแก้ไขจุดด้วย ที่นี่ฉันอยู่ที่ไหน: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array = source.GetRasterBand(1).ReadAsArray() # Close raster source = None # Compute mid-point grid spacings ax = array([gt[0] + ix*gt[1] + gt[1]/2.0 for ix in range(nx)]) ay …

12 python  algorithm  interpolation