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

ฉันพยายามเขียนหลายเลเยอร์ไปยัง GeoPackage เดียวกันใน R แต่ฉันได้รับข้อผิดพลาดCreation of output file failed. ฉันพยายามค้นหาเอกสารการอ่านและเขียนไปยังไฟล์. gpkg ด้วย RGDAL เช่นเพื่อค้นหาว่าwriteOGR()รองรับหลายเลเยอร์จริงหรือไม่ เป็นไปได้ไหมถ้าเป็นเช่นนั้นจะทำอย่างไร? ตัวอย่างการทำงานที่น้อยที่สุด: library(sp) library(maptools) library(rgdal) data(wrld_simpl) norway <- wrld_simpl[wrld_simpl$NAME == "Norway", ] sweden <- wrld_simpl[wrld_simpl$NAME == "Sweden", ] file <- tempfile("scandinavia", fileext = c(".gpkg")) writeOGR(norway, dsn = file, layer = "norway", driver = "GPKG") writeOGR(sweden, dsn …

11 gdal  r  rgdal  geopackage  sf 

ฉันพยายามสร้างรูปหลายเหลี่ยมรูปร่างจากรายการพิกัดโดยใช้เครื่องมือหลามโอเพนซอร์ส ต่อไปนี้คือสิ่งที่ฉันได้เพื่อให้ห่างไกลที่ถูกแฮ็กด้วยกันจากงูหลาม GDAL / OGR ตำราและGIS นี้ตอบ มีคำถามที่คล้ายกันPython: วิธีการสร้างรูปหลายเหลี่ยม Shapefile จากรายการของพิกัด X, Y pyshpแม้ว่าคำถามนี้จะหมายถึงการใช้ อย่างไรก็ตามฉันสนใจที่จะสร้างไฟล์รูปหลายเหลี่ยมโดยใช้เครื่องมือหลาม gdal / ogr เท่านั้น import ogr def create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in coords: ring.AddPoint(coord[0], coord[1]) # Create polygon poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) return poly.ExportToWkt() def write_shapefile(poly, out_shp): """ https://gis.stackexchange.com/a/52708/8104 """ # Now convert …

10 python  gdal  polygon  ogr 

แนวคิดของคอมโพสิตฮิลล์เชดคือการรวมหลายฮิลส์เชดกับการวางแนวแสงแดดที่แตกต่างกันเพื่อหลีกเลี่ยงจุดบอด มันต้องรวมภาพฮิลส์ 3 ภาพด้วย315⁰ (แสงแดด NW, เลเยอร์เริ่มต้น), + 355⁰เลเยอร์, ​​+ 275⁰แสงตะวันตามลำดับ นี้รายละเอียดเพิ่มขึ้นและความสง่างามเต็มรูปแบบที่ไม่ใช่ GIS ดังแสดง below.A คำอธิบายอยู่ที่นี่ การใช้ไฟล์ srtm ใด ๆ ฮิลล์เชดเริ่มต้น 3 ตัวสามารถทำได้ผ่าน: gdaldem hillshade input.tif hillshades_A.tmp.tif -s 111120 -z 5 -az 315 -alt 60 -compute_edges gdaldem hillshade input.tif hillshades_B.tmp.tif -s 111120 -z 5 -az 355 -alt 60 -compute_edges gdaldem …

10 gdal  dem  hillshade  gdal-calc 

ฉันจะได้รับพิกัดที่ฉายรวมถึงค่าพิกเซลจริงที่พิกัดเหล่านั้นจากไฟล์ GeoTiff จากนั้นบันทึกลงในอาร์เรย์แบบ numpy ได้อย่างไร ฉันมีไฟล์ arsenci020l.tif และพิกัดอยู่ในหน่วยเมตร ด้านล่างคือผลลัพธ์ย่อของ gdalinfo ที่ฉันวิ่งไป ~$ gdalinfo arsenci020l.tif Driver: GTiff/GeoTIFF Files: arsenci020l.tif arsenci020l.tfw Size is 10366, 7273 Coordinate System is: PROJCS["Lambert Azimuthal Equal Area projection with arbitrary plane grid; projection center 100.0 degrees W, 45.0 degrees N", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], …

10 python  gdal  geotiff-tiff  numpy 

เป็นไปได้ไหมที่จะส่งข้อมูลโดยตรงไปที่ogr2ogr? ฉันพยายามcurlไฟล์ GeoJSON และไพพ์ไปที่ogr2ogrเพื่อแปลงเป็นรูปแบบอื่น โดยไม่ต้องใช้ไพพ์ฉันสามารถทำสิ่งนี้ได้โดยการเขียนไฟล์ก่อนแล้วแปลงมันโดยใช้คำสั่งต่อไปนี้ แต่ฉันอยากรู้ว่ามันเป็นไปได้ไหมที่จะข้ามขั้นตอนนั้น curl "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson" -o "countries.geojson" && ogr2ogr -f "KML" countries.kml countries.geojson ฉันยังพยายามสร้างตัวแปรด้วยเอาต์พุต curl แทนที่จะเป็น piping แต่ogr2ogrไม่สามารถอ่านได้: DATA=$(curl -s "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson") && ogr2ogr -f "KML" countries.kml $DATA

10 gdal  ogr2ogr 

ฉันต้องการที่จะ georeference แรสเตอร์ที่ใช้และpython GDALแนวทางปัจจุบันของฉันคือการโทรgdal_translateและgdalwarpใช้os.systemและรายการจุดควบคุมภาคพื้นดินที่น่าเกลียด ฉันต้องการวิธีการทำสิ่งนี้ภายในpythonจริงๆ นี่เป็นกระบวนการปัจจุบันที่ฉันใช้: import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"') os.system('gdalwarp -r bilinear -tps -co COMPRESS=NONE "outraster2.tif" "outraster3.tif"') มีก่อนหน้านี้คือคำถามและคำตอบจาก 2012 ซึ่งระบุว่าสามารถเข้าถึงได้หลังจากการนำเข้าgdal_translate gdalฉันไม่แน่ใจว่าล้าสมัยหรือว่าผิด แต่เมื่อฉันรันfrom osgeo import gdalฉันไม่เห็นgdal.gdal_translateตัวเลือก ฉันไม่รู้ว่ามันมีอยู่หรือไม่ แต่ฉันจะรักถ้าฉันสามารถแปลและปฏิเสธ rasters ในแบบ pythonic ตัวอย่างเช่น: # translate gcp_points = [(1251.92, …

10 python  raster  gdal  coordinate-system  gdalwarp 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันพยายามรวม 14 geotiff แบบนี้: geotiff แต่ละอันมีค่าประมาณ 50Mb ฉันต้องการ geotiff ที่เอาต์พุต เวิร์กโฟลว์ของฉัน: gdalbuildvrt -input_file_list list.txt test.vrt (ที่รายการของฉันมีชื่อของ tifs) จากนั้น: gdal_translate -of Gtiff test.vrt test.tif Input file size is 79841, 59955 มันใช้งานได้ แต่ผลที่ได้คือ geotiff ที่ 13,3 Gb! สำหรับไฟล์ 14 ไฟล์แต่ละ 50 Mb ฉันพยายาม Geotiff ที่ …

10 raster  gdal  geotiff-tiff  gdal-translate  gdalbuildvrt 

ฉันมีข้อมูล GML บางส่วนที่มีคอลัมน์ที่ไม่จำเป็นในนั้น ฉันต้องการลบคอลัมน์เหล่านี้เพื่อลดขนาดไฟล์ (เพื่อชี้แจง: ฉันหมายความว่าไฟล์มีหลายจุดในนั้นและแต่ละจุดเหล่านี้มีข้อมูลเมตาบางส่วนติดอยู่และฉันไม่ต้องการข้อมูลเมตาทั้งหมด) ฉันต้องการทำสิ่งนี้จากบรรทัดคำสั่งเนื่องจากฉันมีไฟล์จำนวนมาก ฉันสามารถแปลงไฟล์ GML เป็นรูปแบบต่าง ๆ ได้โดยใช้ogr2ogrดังนั้นคำถามของฉันคือ: วิธีที่ง่ายที่สุดในการลบคอลัมน์ออกจากไฟล์ข้อมูลทางภูมิศาสตร์จากบรรทัดคำสั่งคืออะไร ฉันเป็นโปรแกรมเมอร์ Python ดังนั้นสัญชาตญาณของฉันคือการแปลงไฟล์เป็น GeoJSON จากนั้นใช้วิธีการแยกวิเคราะห์ JSON ของ Python เพื่ออ่านและเขียนไฟล์ JSON อีกครั้ง แต่อาจมีวิธีที่ตรงกว่านี้โดยใช้เครื่องมือ GDAL หรือคล้ายกัน

10 gdal  data  geojson  gml  file-size 

ฉันมีไฟล์ NetCDF ซึ่งมีข้อมูลตามจำนวนวันที่ฝนตกในสหราชอาณาจักร ฉันต้องการแปลงเป็น GeoTIFF แต่ไม่สามารถรับตำแหน่งทางภูมิศาสตร์ของ GDAL ได้อย่างถูกต้อง ดูการถ่ายโอนข้อมูลจาก gdalinfo ด้านล่าง ฉันต้องการ TIFF หนึ่งวงที่แสดงตัวแปร "raindays_tmean_abs" ตัวแปร "lat" และ "lon" ดูเหมือนจะมีค่า min / max ที่สมเหตุสมผลเพื่อนำข้อมูลไปทั่วสหราชอาณาจักร แต่เมื่อฉันแปลงด้วย gdal_translate มันจะจบลงที่ชายฝั่งตะวันตกของแอฟริกา มีคำเตือนเกี่ยวกับมิติแรกที่ไม่ใช่มิติ "เวลาหรือแนวตั้ง" ฉันอ่านว่า gdal คาดหวังว่าข้อมูลจะอยู่ในลำดับที่แน่นอนดังนั้นนี่อาจเป็นปัญหา แต่ฉันไม่ทราบวิธีการแก้ไข ฉันจะแปลง / ปฏิเสธข้อมูลนี้ได้อย่างไร gdalinfo NETCDF: RainDays1_1961-1990_LTA_25km.nc Warning 1: dimension #0 (meaning_period) is not a Time or Vertical …

10 gdal  coordinate-system  netcdf 

คำถามต่อไปนี้เกี่ยวกับการจัดการ Aerials ด้วย ArcGIS: รูปแบบที่มีประสิทธิภาพที่สุดในการจัดการภาพถ่ายทางอากาศเพื่อการดูเท่านั้น ดูเหมือนว่ามี 2 ตัวเลือกที่สำคัญสำหรับการจัดเก็บ / resampling / reprojecting ฯลฯ เสาอากาศ: JP2000 / JP2 / JPEG 2000 (เร็ว ๆ นี้ 5 รหัสสำหรับการจัดการ GDAL) ECW (ERDAS Wavelets ที่ถูกบีบอัด (.ecw)) อื่น ๆ ที่ฉันพลาด? สิ่งที่ฉันเข้าใจขึ้นอยู่กับรุ่น QGIS สำหรับทั้งสองอย่างนั้นมักจะต้องมีการติดตั้งไลบรารีเพิ่มเติมบางอย่าง ECW มีข้อ จำกัด บางประการ - สำหรับการบีบอัดข้อมูลจำเป็นต้องซื้อใบอนุญาตหรือไม่? ฉันทดสอบ jpeg ที่ฉันไม่สามารถใช้สำหรับไฟล์ขนาดใหญ่ (ข้อ จำกัด ด้านขนาดสูงสุด) และมันก็ช้าด้วยขนาดที่ใหญ่กว่า …

10 qgis  gdal  ecw  aerial-photography  jpeg-2000 

ฉันพยายามเรียกใช้รหัส Python โดยใช้หน้าต่าง Python ใน ArcGIS 10.1 ที่ใช้ทั้งโมดูลarcpyและ gdalอย่างไรก็ตามเมื่อฉันลองและนำเข้าgdalโมดูลฉันได้รับข้อผิดพลาด: ImportError: No module named osgeo เห็นได้ชัดว่ามันหาโมดูลไม่ได้ดังนั้นฉันจึงเพิ่มไดเรกทอรีหลักของแพ็คเกจเว็บไซต์หลามลงในsys.pathรายการ: sys.path.append(r"C:\Python27\lib\site-packages") อย่างไรก็ตามเมื่อฉันลองและเรียกใช้import osgeoหรือfrom osgeo import gdalฉันได้รับข้อผิดพลาดต่อไปนี้: Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <module> _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swig_import_helper _mod …

10 arcgis-desktop  arcgis-10.1  gdal  importerror  python-window 

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

ฉันจะถ่ายภาพและแก้ไขภาพโดยอัตโนมัติด้วยซอฟต์แวร์โอเพ่นซอร์สได้อย่างไร ข้อมูลที่เกี่ยวข้องกับแต่ละภาพคือ: รายละเอียดกล้อง ละติจูด, ลองจิจูด, ระดับความสูงของกล้อง การวางแนวกล้อง: หันหน้าไปทาง UAV, ระยะพิทซ์และหันเหจะได้รับ ฉันต้องการส่งออกไฟล์ geotiff และหวังว่าจะได้ภาพที่ประกอบเข้าด้วยกันเพื่อสร้างแผนที่

10 raster  gdal  osgeo  photogrammetry  unmanned-aerial-vehicle 

ฉันต้องการสร้างภาพ multispectral จาก cero เพื่อทำการทดสอบบางอย่างกับมัน บางสิ่งที่ง่ายจริงๆเช่น 5 วงที่มีลักษณะเหมือนกันอย่างสมบูรณ์ด้วยเสียงเกลือและพริกไทยบนพวกเขาหรือค่าที่แตกต่างกันที่จุดศูนย์กลาง เห็นได้ชัดว่านี่จะเป็นสแต็กของเมทริกซ์ซึ่งเป็นอาร์เรย์หลายมิติซึ่งค่อนข้างตรงไปตรงมาเพื่อสร้าง ฉันต้องการประสบความสำเร็จในการใช้ python และ gdal แต่ gdal นั้นค่อนข้างลึกลับฉันไม่เข้าใจเลย ฉันต้องการสร้างไฟล์พิกัด มีใครช่วยฉันได้บ้าง พอยน์เตอร์หรือแบบฝึกหัด gdal บางอันอ่อนโยนมาก ขอบคุณทุกคน.

10 python  raster  gdal  imagery 

ฉันได้ทำตามคำแนะนำความช่วยเหลือออนไลน์ต่าง ๆ เกี่ยวกับวิธีทำให้ ecw สำหรับ QGIS 1.8 ทำงานใน Ubuntu 12.04 ฉันใช้งานได้ทุกอย่างใน Ubuntu 10.04 แต่เนื่องจากการโยกย้ายไปยัง 12.04 ไม่มีโชคเช่นนั้น มีคำแนะนำที่ใช้ได้กับ Ubuntu 12.04 หรือไม่? ในขณะนี้ดูเหมือนว่าจะง่ายกว่าในการติดตั้งรุ่น windows ที่มีการสนับสนุน ecw ในการติดตั้ง แต่ฉันต้องการอยู่กับ Ubuntu

10 qgis  gdal  ubuntu  references  ecw