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

ฉันพยายามใช้ ogrinfo เพื่อรับรายละเอียดบางอย่างเกี่ยวกับรูปร่างไฟล์ที่ฉันดาวน์โหลด ขณะนี้วิธีเดียวที่ฉันรู้วิธีการทำเช่นนี้คือการโหลดลงใน QGIS และคลิกด้วยตนเองเพื่อค้นหาข้อมูลใด ๆ เช่นเปิดตารางแอตทริบิวต์ ฉันแค่ต้องการให้สามารถดูข้อมูลเมตาใด ๆ ที่ติดแท็กพร้อมกับคุณสมบัติ ถ้าฉันทำ: ogrinfo -al USA_adm0.shp ฉันเห็นได้ตั้งแต่ต้นว่ามีข้อมูลที่มีประโยชน์มากมาย แต่แล้วมันก็บินผ่านไปพร้อมกับข้อมูลคุณสมบัติทั้งหมด มีคนช่วยฉันได้ไหม แก้ไข นี่คือสิ่งที่ฉันได้รับบน mac ของฉันโดยใช้แฟล็ก -ro และ -so ดูเหมือนจะไม่ช่วยอะไรมาก ->ogrinfo -ro -so USA_adm0.shp INFO: Open of `USA_adm0.shp' using driver `ESRI Shapefile' successful. 1: USA_adm0 (Polygon)

12 shapefile  gdal  ogr 

ฉันมีชุดแรสเตอร์ (.vrt) พร้อมข้อมูลความชื้นในดินทุกวัน ฉันต้องการที่จะรวมพิกเซลของแรสเตอร์ทั้งหมดเพื่อให้มีการวัดตามเดือน อย่างไรก็ตามโลกทั้งใบไม่ได้ครอบคลุมในแต่ละวันซึ่งส่งผลให้เกิดค่านิยมในสถานที่ที่มีการวัดสำหรับวันอื่น ๆ สิ่งที่ฉันต้องการทำคือผลรวมของค่าแรสเตอร์แต่ละอัน อย่างไรก็ตามดูเหมือนว่าในแต่ละครั้งที่พิกเซล nodata อยู่ในผลรวมพิกเซลที่ได้จะถูกจัดประเภทเป็นโนดาต้าโดยตรง ฉันอยากจะมีสิ่งที่ตรงกันข้าม: ไม่สนใจคุณค่าของสิ่งที่มีอยู่ทั้งหมดและรวมส่วนที่เหลือ ฉันคิดถึงวิธีแก้ปัญหา 2 วิธี: การรวม rasters โดยไม่สนใจค่าของโหนด การแปลงพิกเซลโนดาต้าให้เป็นค่า 0 จากนั้นรวมแรสเตอร์ทั้งหมด น่าเสียดายที่ฉันไม่พบเครื่องมือใด ๆ ในการทำเช่นนี้ มีใครช่วยฉันบ้าง

12 raster  gdal  qgis-2 

มีวิธีในการทำงานเช่นเดียวกับยูทิลิตี gdalbuildvrt โดยใช้การผูก GDAL Python หรือไม่ จนถึงตอนนี้ฉันยังไม่พบวิธีการอื่นนอกเหนือจากการสร้าง vrt ของชุดข้อมูลเดียวและแก้ไข xml ด้วยตนเอง ฉันต้องการสร้าง vrt จากหลายภาพแรสเตอร์ เป็นไปได้ไหมถ้าใช้ Python แท้? ตัวเลือกอื่น ๆ ของฉันคือใช้ subprocess เพื่อเรียก gdalbuildvrt

12 python  raster  gdal  vrt  gdalbuildvrt 

ฉันกำลังมองหาฟังก์ชั่นที่คล้ายกับ GRASS r.resample ฉันมีโฟลเดอร์ที่เต็มไปด้วยไฟล์ TIF และต้องการลดความละเอียด ("ลดขนาดไฟล์ลง") เพื่อลดขนาดไฟล์

12 raster  gdal  ogr 

ฉันกำลังทำงานกับเครื่องมือที่แปลงเลเยอร์แรสเตอร์เป็นอาร์เรย์สำหรับการประมวลผลด้วย NumPy และในอุดมคติฉันต้องการที่จะทำงานกับ rasters ที่มาในแพ็คเกจด้วย. gdb โดยไม่ต้องส่งออกทั้งหมด (โดยเฉพาะถ้าต้องการ ArcGIS หรือ ArcPy ) เป็นไปได้กับไดรเวอร์ OpenFileGDB หรือไม่ จากสิ่งที่ฉันบอกได้ว่าไดรเวอร์นี้ดูเหมือนว่าจะรักษาเลเยอร์แรสเตอร์เหมือนกับเลเยอร์เวกเตอร์ซึ่งให้คุณเข้าถึงข้อมูลบางอย่างเกี่ยวกับเลเยอร์ แต่ไม่ได้ให้ฟังก์ชั่น ReadAsArray

12 python  gdal  file-geodatabase 

ฉันรัน gdal2tiles บนเครื่องเดเบียนที่มี 4 คอร์: sudo /usr/bin/gdal2tiles.py -r cubic -a 0,0,0 -z 10-15 /home/adm/topo/ekb.vrt /var/www/tiles ในhtopฉันเห็น: ทำไม gdal ใช้แกนเดียวในการทำงาน รุ่น Gdal คือ 2.1.2

11 gdal 

ฉันต้องการเร่งความเร็วกระบวนการที่สอดคล้องกับคำสั่งนี้: gdalwarp -t_srs EPSG:4326 -overwrite input.ntf output.tif มีวิธีใช้การประมวลผลแบบขนานใน GDAL หรือไม่ กึ่งนอกหัวข้อ: ถ้าไม่คุณแนะนำวิธีแก้ปัญหาที่ไม่ใช่ GDAL เพื่อเร่งกระบวนการ gdalwarp เหมือน? นี่คือหน้าเว็บที่ฉันได้ดู: เป็นไปได้ไหมที่จะทำการประมวลผลแบบขนานใน GDAL และ QGIS? [gdal-dev] การประมวลผลแรสเตอร์ของ GDAL: การคำนวณแบบขนาน แก้ไข: ฉันถามคำถามนี้เพราะฉันคิดว่าฉันเห็น CPU% อยู่ต่ำกว่า 100% เมื่อประมวลผล gdalwarp อย่างไรก็ตามในการค้นหาครั้งที่สองก็มาถึง 555.5%

11 gdal  parallel-processing 

ฉันมีไฟล์ในการฉายภาพ British National Grid: Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936", SPHEROID["Airy_1830",6377563.396,299.3249646]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",49], PARAMETER["central_meridian",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["Meter",1]] cat: Integer (9.0) ฉันสามารถใช้ GDAL / OGR เพื่อให้ได้พื้นที่รวมของรูปหลายเหลี่ยมทั้งหมดใน shapefile เป็นเฮกตาร์หรือไม่ ฉันสงสัยว่ามันเป็นไปได้ด้วย-sqlหรือเปล่า: ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp …

11 gdal  ogr  area  ogrinfo 

ฉันมีข้อมูลจำนวนมากในรูปแบบของ CSV ที่มีพิกัด UTM XและYและค่าระดับความสูงเป็นZข้อมูล ฉันจำเป็นต้องแปลงข้อมูลเหล่านี้เป็น DEM เป็น GeoTIFF เพื่อการวิเคราะห์เพิ่มเติม ในกรณีนี้จำนวนมากหมายถึง 16 เมตร เส้นที่มีในจุดหนึ่งX, YและZต่อบรรทัด คะแนนมีการกระจายอย่างเท่าเทียมกันดังนั้นจึงไม่จำเป็นต้องมีการแก้ไข แต่ละจุดจะต้องถูกแปลงเป็นเซลล์แรสเตอร์ ข้อมูลต้นฉบับมาโดยไม่มีตัวคั่นพร้อมความกว้างคอลัมน์คงที่ ฉันคิดแล้วออกวิธีการแปลงไฟล์ไวยากรณ์ที่จะใช้คั่นแทนความกว้างคงที่และกำจัดตัวละครทุกพื้นที่โดยใช้โปรแกรมแก้ไขข้อความกระแสsed จากที่นี่ที่ปกติขั้นตอนการทำงานของฉันจะเป็นที่จะนำเข้าข้อมูลลงใน ArcGIS โดยการสร้างระดับคุณลักษณะจากX, YและZข้อมูลและเป็นขั้นตอนที่สองแปลง shapefile จุดเป็น GeoTIFF ใช้ชี้ไปที่ Rasterเครื่องมือ อย่างไรก็ตามไฟล์ที่ฉันมีในขณะนี้มีขนาดใหญ่เกินไปสำหรับกระบวนการนี้ แทนที่จะเป็นกระบวนการทำงานที่อธิบายไว้ข้างต้นฉันกำลังมองหาทางเลือกที่มีประสิทธิภาพและค้นพบ GDAL อย่างไรก็ตามในgdal_translateรูปแบบที่สนับสนุนที่ใกล้เคียงที่สุดที่ฉันสามารถหาได้ในรายการประเภทไฟล์ที่รองรับคือตาราง ASCII แต่ไม่มี XYZ คั่นด้วยเครื่องหมายจุลภาค ปัญหาอีกประการหนึ่งคือฉันมีพิกัด UTMในขณะที่ตัวอย่างส่วนใหญ่ดูเหมือนจะใช้พิกัดองศาแบบทศนิยม อย่างไรก็ตามฉันต้องอยู่ในระบบ UTM (หรืออย่างน้อย GeoTIFF ของฉันจะต้องอยู่ในระบบพิกัด UTM) ดังนั้นฉันกำลังมองหาวิธีแปลง CSV XYZ เป็น GeoTIFF โดยใช้ …

11 gdal  geotiff-tiff  csv  gdal-translate  xyz 

สคริปต์พื้นฐานเพื่อวนซ้ำแบบซ้ำในโฟลเดอร์ย่อยและรวมรูปร่างไฟล์ทั้งหมดไว้ในไฟล์เดียวคือ: #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f "ESRI Shapefile" $consolidated_file $i else # update the output file with new file content ogr2ogr -f "ESRI Shapefile" -update -append $consolidated_file …

11 gdal  ogr2ogr  merge  bash  mergeshapes 

ฉันกำลังมองหาการปรับขนาดภาพ TIFF ที่ไม่ได้ลงชื่อ 16 บิตเป็นภาพ 8 บิตโดยใช้ GDAL (OSGeo4W) อย่างไรก็ตามฉันต้องการ จำกัด ค่าพิกเซลเพื่อให้อยู่ในช่วงของภาพ 8 บิต สิ่งที่ฉันหมายถึงคือฉันต้องการให้แน่ใจว่าสัดส่วนภาพในภาพ 16 บิตนั้นแสดงในรูป 8 บิต (0-255) การใช้ gdal_translate เพื่อแปลงจาก 16 บิตเป็น 8 บิตจะตัดค่าพิกเซลและไม่สร้างภาพ RGB ฉันไม่แน่ใจในสิ่งที่ฟังก์ชั่นและการตั้งค่าที่จะโทรจากภายใน gdal_translate หรือถ้ามีตัวเลือกที่ดีกว่า GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 …

11 gdal  imagery  gdal-translate 

ฉันพยายามเรียนรู้วิธีการประมวลผลภาพจากการรับรู้ระยะไกลโดยใช้การเชื่อมโยง Python GDAL และการถ่ายภาพจำนวนมาก เป็นครั้งแรกที่ฉันอ่านไฟล์ geotiff Landsat8 ทำการจัดการอย่างง่ายและเขียนผลลัพธ์ลงในไฟล์ใหม่ รหัสด้านล่างดูเหมือนว่าจะทำงานได้ดียกเว้นว่าแรสเตอร์ดั้งเดิมถูกทิ้งในไฟล์เอาต์พุตแทนที่จะเป็นแรสเตอร์ที่ได้รับการจัดการ ความคิดเห็นหรือข้อเสนอแนะใด ๆ ยินดีต้อนรับ แต่โดยเฉพาะอย่างยิ่งบันทึกเกี่ยวกับสาเหตุที่แรสเตอร์ที่จัดการไม่แสดงในผลลัพธ์ import os import gdal gdal.AllRegister() file = "c:\~\LC81980242015071LGN00.tiff" (fileRoot, fileExt) = os.path.splitext(file) outFileName = fileRoot + "_mod" + fileExt ds = gdal.Open(file) band = ds.GetRasterBand(1) arr = band.ReadAsArray() [cols, rows] = arr.shape arr_min = arr.Min() arr_max = …

11 python  gdal  geotiff-tiff  landsat  numpy 

ฉันกำลังดิ้นรนกับการประมวลผลและหวังว่าฉันจะสามารถแก้ไขได้ที่นี่ ฉันทำงานกับ Remote Sensing ที่ใช้กับการป่าไม้โดยเฉพาะการทำงานกับข้อมูล LiDAR แนวคิดคือการใช้ Scikit-image สำหรับการตรวจจับต้นไม้ด้านบน เนื่องจากฉันใหม่ใน Python ฉันถือว่าชัยชนะส่วนตัวที่ยอดเยี่ยมในการทำสิ่งต่อไปนี้: นำเข้า CHM (พร้อม matplotlib); เรียกใช้ตัวกรอง Gaussian (พร้อมแพ็คเกจ scikit-image); เรียกใช้ตัวกรอง maxima (พร้อมแพ็คเกจ scikit-image); รัน peak_local_max (พร้อมแพ็กเกจ scikit-image); แสดง CHM ​​ด้วย maxima ท้องถิ่น (พร้อม matplotlib); ตอนนี้ปัญหาของฉัน เมื่อฉันนำเข้าด้วย matplot รูปภาพจะสูญเสียพิกัดทางภูมิศาสตร์ ดังนั้นพิกัดที่ฉันมีจึงเป็นเพียงพิกัดรูปภาพพื้นฐาน (เช่น 250,312) สิ่งที่ฉันต้องการคือการรับค่าของพิกเซลภายใต้จุดสูงสุดในภาพ (จุดสีแดงในภาพ) ที่นี่ในฟอรัมฉันเห็นชายคนหนึ่งถามสิ่งเดียวกัน ( รับค่าพิกเซลของแรสเตอร์ GDAL ภายใต้จุด OGR …

11 python  gdal  remote-sensing 

ฉันต้องการเรียกใช้ GDAL QGIS บนRaspberry Piเพื่อทำหน้าที่เป็นเครื่องมือประมวลผล GIS ที่ราคาไม่แพง แผนการจะซื้อประมาณ 50 Pis ใครมีลิงค์หรือข้อมูลเกี่ยวกับการติดตั้งและใช้งาน GDAL บน Raspberry Pi หรือเพียงแค่ติดตั้ง GDAL ในสภาพแวดล้อม RASPBIAN อย่างดีที่สุดที่ฉันเห็น Raspbian เป็นตัวเลือกอันดับต้น ๆ ของฉันและเป็น Debian ที่ใช้ แต่ ARCH, RASPBMC, PIDORA, RISC OS และ OPENELEC เป็นรสชาติของ linux ที่เปิดให้ฉันออกจากกล่อง หวังว่าด้วย Python

11 gdal  linux  raspberry-pi 

ฉันมีรูปภาพ RapidEye จำนวนมาก (+300 Gb) ที่ฉันต้องการใช้เป็นแผนที่ฐานสำหรับแอปพลิเคชันออนไลน์ จนถึงตอนนี้ฉันได้ทำหน้าที่เป็น ImagePyramid ใน geoserver โดยทำตามขั้นตอนที่แสดงในงานนำเสนอที่เรียกว่า "Geoserver on steroids" ปัญหาคือฉันไม่สามารถทำสมดุลสีที่เหมาะสมดังนั้นเมื่อดูชุดข้อมูลทั้งหมดโมเสคดูไม่ดีมาก แต่เมื่อฉันซูมเข้าจะดีขึ้น (ภาพเบลล์) กระบวนการของฉันคือ แปลงภาพทั้งหมดเป็น 8 บิต, epsg 4326, การแก้ไขลูกบาศก์โดยใช้ gdal สร้างแรสเตอร์เสมือน gdal ด้วยภาพทั้งหมด (gdalbuildvrt) สร้างปิรามิดของภาพ (gdal_retile พร้อมการบีบอัดและการไถแบบ geotiff) และเผยแพร่บน geoserver (โดยใช้การยืดฮิสโตแกรมในสไตล์เลเยอร์) เคล็ดลับใด ๆ เกี่ยวกับวิธีปรับปรุงกระบวนการและให้ได้สมดุลสีที่ดีขึ้น

11 geoserver  gdal  basemap  image-pyramids