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

ฉันพยายามที่จะสร้างแผนที่ฐานไทล์ ฉันมีภาพที่ฉันต้องการสร้างไทล์และฉันใช้ GDAL ฉันจัดการเพื่อสร้างไทล์ แต่เมื่อฉันดูผลลัพธ์ดูเหมือนว่าแผนที่แบบเรียงต่อกันคือ: มันพลิก สัดส่วนไม่ถูกต้อง นี่คือภาพซ้อนทับแผนที่: นี่คือผลของการเรียงต่อกัน: ฉันตรวจสอบหลาย ๆ ครั้งว่าพิกัดนั้นดีโดยการเปรียบเทียบมุมภาพกับภาพออร์โธกราฟ อีกอย่างคือใน google earth ที่มีพิกัดเดียวกันภาพแผนที่จะพอดีกับแผนที่ google แก้ไข gdalinfo map_overlay.png ด้วยผลลัพธ์นี้: ไดร์เวอร์: PNG / ไฟล์กราฟิกเครือข่ายพกพา: map_overlay.png ขนาดคือ 3527, 2494 ระบบพิกัดคือ `'โครงสร้างรูปภาพข้อมูลเมตา: INTERLEAVE = พิกเซลพิกัดมุม: ด้านซ้าย (0.0, 0.0) มุมซ้ายล่าง (0.0, 2494.0) มุมขวาบน (3527.0 , 0.0) มุมขวาล่าง (3527.0, 2494.0) ศูนย์ (1763.5, 1247.0) …

10 gdal  tiles 

ฉันยังใหม่กับการเขียนโปรแกรมเชิงพื้นที่ของหลาม ฉันรันสคริปต์ต่อไปนี้และได้รับข้อความแสดงข้อผิดพลาดที่เกี่ยวข้อง >>> import osgeo >>> import osgeo.ogr >>> shapefile = osgeo.ogr.Open("tl_2009_us_state.shp") >>> numLayers = shapefile.GetLayerCount() Traceback (most recent call last): File "<pyshell#5>", line 1, in <module> numLayers = shapefile.GetLayerCount() AttributeError: 'NoneType' object has no attribute 'GetLayerCount'

10 python  gdal 

ฉันพยายามแสดงไฟล์แรสเตอร์ความละเอียดขนาดเล็กในเครื่องมือ QGIS Print Manager ฉันต้องการแสดงป้ายกำกับขนาดเล็กที่มีค่าแต่ละเซลล์ (ตำแหน่ง - ประมาณกึ่งกลางของแต่ละเซลล์) ฉันจะทำมันได้อย่างไร มีปลั๊กอินใดบ้างที่สามารถช่วยให้ฉันทำงานนี้ได้สำเร็จ ความคิดของฉันคือการใช้ gdal2xyz โหลดไฟล์ดังกล่าวโดยใช้ปลั๊กอิน CSV แล้วบันทึกเป็นไฟล์เวกเตอร์ น่าเสียดายที่มันใช้งานได้ดีมากเพราะฉันมีไฟล์แรสเตอร์มากมาย อีกแนวคิดหนึ่งคือการใช้ gdal2xyz.py แล้ว ogr2ogr เพื่อบันทึกเป็นไฟล์เวกเตอร์ มีวิธีการชุดกระบวนการนี้หรือไม่ มีโปรแกรมโอเพนซอร์ซใดที่ฉันสามารถทำได้ง่ายกว่านี้อีกไหม?

10 qgis  raster  gdal  ogr 

ฉันโหลดรูปร่างไฟล์ไปยังฐานข้อมูล PostGIS ได้สำเร็จ ฉันต้องการที่จะดึงรูปทรงเรขาคณิตและเลเยอร์ที่เฉพาะเจาะจงเป็นไฟล์ dxf ตามแบบสอบถาม SQL ฉันใช้ PostGIS, GDAL, PostgreSQL และ Python จากชุด OSGeo4W บน Windows ฉันคิดว่าฉันควรใช้ ogr2ogr เล็กน้อยเช่นนี้: $ ogr2ogr -sql "SELECT * FROM some_table" -f "dxf" outputfile.dxf -overwrite "PG:dbname=mydb user=dbusername password=dbpassword" แต่เมื่อฉันทำเช่นนั้นด้วยค่าที่เหมาะสมฉันได้รับข้อผิดพลาดแบบนี้: ERROR 1: DXF layer does not support arbitrary field creation, field 'objectid' not created. ส่งผลให้เกิดข้อผิดพลาดสำหรับทุกคอลัมน์คุณลักษณะในชั้นsome_tableแต่ฉันจะได้รับไฟล์ …

10 postgis  gdal  ogr  dxf 

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

10 python  raster  gdal 

ฉันพยายามอย่างหนักในการสร้างไทล์สำหรับภาพความละเอียดสูงที่เรามี ภาพปัจจุบันที่เรามีคือภาพที่มีขนาดใหญ่มาก (+ 20GB) ซึ่งบันทึกเป็นไฟล์ GeoTiff ฉันต้องการสร้างไทล์โดยใช้ยูทิลิตีบรรทัดคำสั่ง gdal2tiles จากนั้นเปิดและดูใน Cesium โดยใช้ผู้ให้บริการภาพ TMS เพื่อจัดเรียงไทล์ ใช้ gdalinfo นี่คือรายละเอียดบางส่วนของภาพ: Driver: GTiff/GeoTIFF Files: image.tif Size is 52250, 56119 Coordinate System is: PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",27], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",10000000], UNIT["meters",1], …

10 gdal  gdalwarp  cesium  gdal2tiles 

ฉันใช้ 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 

ฉันจำเป็นต้องสร้างรูปหลายเหลี่ยมของภาพแรสเตอร์เลเยอร์เดียวหลายรูปไม่ใช่กล่องขอบเขต / ขอบ แต่พื้นที่ที่ไม่มีค่า Nodata ดังแสดงที่นี่: การสร้างรูปร่างไฟล์แสดงรอยเท้าของ Rasters? . ในการตอบคำถามข้างต้นปลั๊กอิน Image Boundaryได้ถูกกล่าวถึง แต่ฉันไม่พบมันใน QGIS 1.8.0 Lisboa ที่ทำงานบน Ubuntu เครื่องมือยังคงมีอยู่หรือไม่? ถ้าไม่มีมีวิธีการทำเช่นนี้โดยใช้เครื่องมือ opensource ไม่ว่าจะเป็น gdal, R, QGIS, GRASS หรือคล้ายกันไม่ใช่กับ ArcMap?

10 gdal  pyqgis  qgis-plugins  footprint 

ฉันต้องการสร้าง.tiffไฟล์อ้างอิงทางภูมิศาสตร์จาก.pngไฟล์ ดังนั้นฉันแปลง. png เป็น.tiffไฟล์ปกติโดยใช้บรรทัดคำสั่งต่อไปนี้: convert image.png image.tiff ที่นี่ทำงานได้ดีฉันมี.tiffไฟล์ที่ถูกต้องซึ่งแสดงผลได้ดี จากนั้นฉันต้องการเพิ่มข้อมูลการอ้างอิงทางภูมิศาสตร์ลงไป ในการทำเช่นนั้นฉันพบเครื่องมือชื่อtiff2geotiffฉันใช้มันเช่นนี้: tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff สิ่งนี้ทำให้ฉันมีภาพสีดำซึ่งนำเข้ามาใน QGis แสดงได้ดีในสถานที่ที่ดียกเว้นว่าเป็นสีดำทุกที่ (พิสูจน์ได้ว่ามีการเพิ่มข้อมูลทางภูมิศาสตร์ที่ได้รับการพิสูจน์มาอย่างดี ) ดังนั้นเมื่อฉัน: gdalinfo georeferenced.tiff ฉันได้รับ: Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength Driver: GTiff/GeoTIFF Files: araujo_tmp.tiff Size is …

10 qgis  gdal  php  geotiff-tiff 

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

ฉันพยายามจัดทำสคริปต์ Python ที่ใช้ GDAL สำหรับ PyPI ฉันเริ่มต้นด้วยการอ้างอิงโดยตรงในsetup.py: install_requires=['GDAL==1.11.2'], วิธีนี้แพ็คเกจไม่สามารถติดตั้งในสภาพแวดล้อมการทดสอบเสมือนจริงของฉัน: extensions/gdal_wrap.cpp:2855:22: fatal error: cpl_port.h: No such file or directory #include "cpl_port.h" ^ compilation terminated. error: Setup script exited with error: command 'x86_64-linux-gnu-gcc' failed with exit status 1 ฉันลองแล้วมีการอ้างอิงถึงpygdalเนื่องจากมีการทำเครื่องหมายเป็นรุ่นที่เป็นมิตร virtualenv: install_requires=['pygdal'], วิธีนี้การติดตั้งจะเสร็จสิ้นโดยไม่มีข้อผิดพลาด (แต่ด้วยการโหลดคำเตือนการคอมไพล์ตามปกติ) อย่างไรก็ตามเมื่อนั้นฉันจึงเรียกใช้สคริปต์ฉันได้รับข้อผิดพลาดนี้กลับมา: Traceback (most recent call last): File "/home/desouslu/.virtualenvs/test_p3/bin/hasc2gml", line …

9 python  gdal 

มีเอกสารอย่างเป็นทางการสำหรับโมดูล python GDAL หรือไม่ มีตำราหลายเล่ม แต่มีข้อมูลไม่เพียงพอในแต่ละฟังก์ชันหรืออาร์กิวเมนต์สำหรับฟังก์ชันพื้นฐานที่สุด

9 python  gdal  documentation 

ฉันพยายามฉายภาพแรสเตอร์ ใน R มีprojectRaster()ฟังก์ชั่นสำหรับสิ่งนี้ (ด้านล่างเป็นตัวอย่างที่ทำซ้ำได้อย่างสมบูรณ์): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') ซึ่งใช้งานได้ดี อย่างไรก็ตามมันค่อนข้างช้า เพื่อเพิ่มความเร็วฉันแม้ว่าจะใช้gdalwarpแทน (ด้วย SSD ค่าใช้จ่ายในการอ่านและเขียนจาก / …

9 raster  gdal  r  coordinate-system  gdalwarp 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันกำลังสร้างสคริปต์ใน python โดยใช้ OGR / GDAL ฉันมีชุดไฟล์รูปร่างและชุดไฟล์แรสเตอร์ GeoTiff ฉันต้องการให้สคริปต์ของฉันเพิกเฉยกับรูปร่างของไฟล์หากพวกเขาไม่ตัดกันกับพื้นที่แรสเตอร์ Shapefile ไม่ได้เป็นรูปสี่เหลี่ยมผืนผ้าดังนั้นฉันไม่สามารถเปรียบเทียบค่า xmin / xmax, ymin / ymax ที่ส่งกลับโดย layer เท่านั้น GetGetExtent () ฉันต้องการรูปหลายเหลี่ยมที่เกิดขึ้นจริงซึ่งเป็นรูปร่างโดยรวมแล้ววิธีพิจารณาว่ารูปหลายเหลี่ยมนั้นตัดกับสี่เหลี่ยมแรสเตอร์หรือไม่ ฉันคิดว่าฉันสามารถรวมรูปหลายเหลี่ยมทั้งหมดใน shapefile เข้ากับคุณลักษณะหนึ่งแล้วอ่านเรขาคณิตในคุณลักษณะนั้นแล้วเปรียบเทียบข้อมูลนั้นกับขอบเขตแรสเตอร์ อย่างไรก็ตามฉันไม่แน่ใจเกี่ยวกับวิธีการดำเนินการนี้โดยเฉพาะ วิธีดึงข้อมูลรูปหลายเหลี่ยมเส้นขอบจาก shapefile? วิธีการตรวจสอบว่ารูปหลายเหลี่ยมนั้นตัดกันพื้นที่สี่เหลี่ยมที่กำหนดหรือไม่

9 python  raster  shapefile  gdal  ogr 

ฉันรวบรวมและติดตั้งต้นแบบ QGIS ปัจจุบันตามปกติและฉันนำเข้าข้อมูลแรสเตอร์ลงในฐานข้อมูล 9.3 ของฉันที่มีส่วนขยาย postgis 2.0.2 เมื่อใดก็ตามที่ฉันพยายามนำเข้าข้อมูลแรสเตอร์เป็น qgis โดยใช้ตัวจัดการฐานข้อมูลฉันได้รับข้อความนี้: GDAL provider: Cannot get GDAL raster band: Raster layer: Provider is not valid (provider: gdal, URI: PG: dbname=gis host=127.0.0.1 user=xxx password=xxx port=5432 mode=2 schema=public table=elevation ในบรรทัดคำสั่งฉันได้รับ ERROR 1: Computed PostGIS Raster dimension is invalid. You've probably specified unappropriate resolution. ฉันไม่รู้ว่าจะค้นหาข้อผิดพลาดได้ที่ไหน …

9 qgis  postgis  raster  gdal  raster2pgsql