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

ฉันมีภาพ. TIFF ที่ฉันพยายามซ้อนทับใน Google Earth เพื่อที่ฉันจะต้อง georeference มันก่อน การฉายภาพคือ UTM ที่มี NAD83 เป็นตัวเลข ฉันรู้ว่าละติจูด / ลองจิจูดของมุมทั้งสี่ของรูปภาพ ฉันกำลังทำวิธีต่อไปนี้: gdalinfo เพื่อรับพิกัดทั้งสี่ gdal_translate เพื่อระบุข้อมูล lat / long gdal2tiles.py เพื่อเรียงภาพ มีเพียงสิ่งเดียวที่ดักฟังฉันคือฉันไม่เคยระบุ UTM และ NAD83 จากนั้นภาพซ้อนทับของฉันจะมาได้อย่างไร ใครทราบวิธีการระบุข้อมูลการฉายนี้

16 gdal  georeferencing 

ฉันมีแหล่งข้อมูลหลายแหล่งที่สร้างขึ้นโดยใช้ gdal2tiles ที่ฉันต้องการนำเสนอบนแผนที่เดียวกัน ดังนั้นเมื่อนำเสนอไทล์ฉันจำเป็นต้องกำหนดแหล่งที่มาที่จะให้บริการจากการตรวจสอบกับขอบเขต ใครรู้วิธีการคำนวณ lat lon bounds สำหรับไทล์เดี่ยวโดยอ้างอิงจากการซูม, x และ y (จากโครงสร้างไฟล์ที่สร้างโดย gdal2tiles)? BTW: ฉันใช้ Google Maps API v3 ในกรณีที่มีความจำเป็นต้องเรียกใช้ฟังก์ชันบางอย่างจาก API เพื่อช่วยในการคำนวณ เหตุผลที่ฉัน จำกัด เพียงการซูม x และ y ล้วนเป็นเพราะไทล์ไม่ใช่เพียงการซ้อนทับบนแผนที่ แต่ยังมีฟังก์ชั่นการพิมพ์แบบกำหนดเองที่อนุญาตให้พิมพ์นอกแผนที่ได้

16 gdal  google-maps  gdal2tiles 

เราเพิ่งได้รับชุด DEM จำนวนมากในที่ทำงานและฉันต้องการสร้างรูปทรงจากพวกเขา DEM นั้นมีความละเอียด 1 ม. และขนาด 1kmx1km ผลลัพธ์จาก gdalinfo: Driver: AAIGrid/Arc/Info ASCII Grid Files: 380000_6888000_1k_1m_DEM_ESRI.asc Size is 1000, 1000 Coordinate System is `' Origin = (380000.000000000000000,6888000.000000000000000) Pixel Size = (1.000000000000000,-1.000000000000000) Corner Coordinates: Upper Left ( 380000.000, 6888000.000) Lower Left ( 380000.000, 6887000.000) Upper Right ( 381000.000, 6888000.000) Lower …

16 raster  gdal  dem  contour 

ตัวอย่างฉันจะใช้ไทล์ต่อไปนี้http://a.tile.openstreetmap.org/3/4/2.pngและบันทึกเป็น "4_2.png" พิกัด WGS84 ของไทล์นี้สามารถคำนวณหรืออ่านได้ที่นั่นโดยคลิกที่ไทล์ที่สอดคล้องกัน: 0 66.51326044311185 45 40.97989806962013 (West North East South) วิธีการระบุตำแหน่งทางภูมิศาสตร์ของกระเบื้องอย่างถูกต้อง (โดยใช้ gdal เพื่อสร้างรูปแบบ geotiff หรือรูปแบบทางภูมิศาสตร์อื่น ๆ ) เพื่อให้: ไม่จำเป็นต้องยืดบิตแมป (= พิกเซลใน geotiff ตรงกับในบิตแมปดั้งเดิม) ภาพที่ได้จะเปิดขึ้นในตำแหน่งที่ถูกต้องใน GIS viewer / editor (เช่นในTatukGIS Free Viewer ) (แก้ไขเมื่อวันที่ 19 กันยายน 2011 เพื่อทำให้คำถามของฉันชัดเจนขึ้นและรวมถึงข้อสรุปของฉัน) บทสรุปของฉัน: ฉันก่อนว่าแนวคิดที่สาม (ดูด้านล่าง) เป็นแนวคิดที่ถูกต้อง ฉันเปิด geotiff ใน GIS Viewer …

16 web-mapping  gdal  epsg  tiles 

ขออภัยหากคำถามต่อไปนี้ค่อนข้างงี่เง่า แต่ฉันเพิ่งใหม่กับสิ่ง GIS ทั้งหมดนี้ ฉันพยายามแปลงภาพ geoTiff ที่คาดการณ์ไว้เป็น WGS84 โดยใช้ gdal ใน python ฉันได้พบโพสต์ที่แสดงกระบวนการในการแปลงคะแนนภายใน GeoTiff ที่ฉายโดยใช้บางอย่างที่คล้ายกับต่อไปนี้: from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], …

16 python  gdal  coordinate-system  wgs84 

ดูเหมือนว่ามันจะตรงไปตรงมา แต่ฉันไม่สามารถใช้ชีวิตของฉันได้ ด้วย GeoTiff ที่มีอยู่แล้วที่มีวงดนตรีเดียวคุณจะเพิ่มวงอื่นใน GeoTiff นั้นได้อย่างไร เป็นวิธีเดียวที่จะทำเช่นนี้เพื่อใช้ GDALCreate () เพื่อสร้างไฟล์ใหม่ที่มีวงดนตรีเพิ่มเติมและจากนั้นคัดลอกวงดนตรีที่มีอยู่จากไฟล์ต้นฉบับไปยังไฟล์ใหม่หรือไม่? นั่นเป็นทางออกเดียวที่ฉันเกิดขึ้นจนถึงตอนนี้ แต่ดูเหมือนว่าควรมีวิธีง่ายๆในการเพิ่มแบนด์ไปยังชุดข้อมูลที่มีอยู่

15 raster  gdal  band 

โอเคดังนั้นฉันจึงรู้สึกเบื่อหน่ายกับโซลูชัน raster-in-database ที่แตกต่างกันซึ่งฉันพยายามในช่วง 2 สัปดาห์ที่ผ่านมาซึ่งไม่สามารถใช้งานได้ในปัจจุบัน (PostGIS raster) หรือไม่สามารถจัดการปริมาณข้อมูลของฉันได้ (Rasterlite) ฉันกำลังอ่านเกี่ยวกับรูปแบบ VRT อดทนกับฉันฉันต้องการความคิดเห็นของคุณเกี่ยวกับเรื่องนี้และถ้าสิ่งที่ฉันเสนอทำให้รู้สึกทั้งหมด: VRT นั้นเป็นไฟล์ XML เมตาดาต้าที่อธิบายคุณสมบัติต่างๆของไฟล์แรสเตอร์จริงเช่นขนาดพิกเซลการระบุตำแหน่งทางภูมิศาสตร์ ฯลฯ ขั้นตอนที่ 1: ฉันแปลง (ไม่ใช่เพียงสร้าง VRTs) ไฟล์ GeoTIFF ทั้งหมดของฉันเป็น VRT ตอนนี้เนื่องจากฉันเห็นว่าฉันสามารถผสาน VRT ได้ฉันจึงคิดว่าเป็นไปได้ที่จะสร้างโมเสกของแฟ้มข้อมูลเมตา VRT ซึ่งเป็น VRT ขนาดใหญ่หนึ่งตัวที่จะเก็บเมทาดาทาทั้งหมด ยอดเยี่ยมนั่นหมายถึงเมตาดาต้าอาจถูกนำมาใช้เพื่อปรับการเรนเดอร์ให้เหมาะสมเนื่องจาก QGIS จะตัดสินใจ (สมมติฐาน) ว่าไทล์ใดที่อยู่ในขอบเขตของผืนผ้าใบปัจจุบันและที่ไม่ ... ขั้นตอนที่ 2: สร้างโมเสก VRT ตกลงตอนนี้ผลประโยชน์ต่อไปที่ฉันเห็นคือฉันต้องการเพิ่มหนึ่งรายการลงในสารบัญของ QGIS และฉันสามารถตั้งค่าสิ่งต่าง ๆ ได้เช่นการแสดงผลตามขนาดและความโปร่งใสของค่าข้อมูลโมฆะถ้าจำเป็น ใคร ๆ …

15 qgis  gdal  vrt 

สั้น: ฉันต้องการเปลี่ยนความละเอียดของแรสเตอร์และปรับทางลาดสีเทาตามที่แสดงในภาพตะโกน การตั้งค่าคือการใช้ GDAL, PIL หรือ Numpy คำอธิบาย: ฉัน Kriging จุดลงไปในเชิงภาพที่มีความละเอียดการส่งออกของ 20 เมตรกับGeostatistical ห้องสมุดที่มีประสิทธิภาพสูง ฉันไม่ต้องการเปลี่ยนความละเอียดเอาต์พุตเนื่องจากเวลาการแก้ไขเพิ่มขึ้นเป็นทวีคูณ ด้วยความละเอียดนี้ภาพที่ส่งออกจะน่าเกลียด (pixelated และ aliased) ฉันไม่รู้ว่ามันถูกต้องในแนวความคิดหรือไม่ แต่ฉันต้องการภาพที่นุ่มนวลเหมือนในตัวอย่างร้อง มันเหมือนกับการ 'ปรับกำลังภาพ' ใหม่ให้เป็นภาพที่มีความละเอียดมากกว่า ฉันใช้ python ดังนั้นการตั้งค่าของฉันคือ GDAL, Python Imaging Library หรือ Numpy คำตอบอาจเป็นเชิงทฤษฎีเช่นชี้ให้เห็นชื่ออัลกอริทึมหรือแนวคิดของการดำเนินการชนิดนี้ ที่มา: ปลายทาง: แก้ไขผลลัพธ์ด้วย gdalwarp cubic spline:

15 python  raster  gdal 

ฉันใช้ GDAL ใน Python สำหรับการนำเข้า shapefile ฉันต้องการทราบชื่อฟิลด์สำหรับไฟล์ด้วยวิธีปัจจุบันของฉันคือ: fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) แต่ด้วยวิธีนี้ฉันได้รับคุณสมบัติสำหรับเลเยอร์แรก มันหมายความว่าเป็นไปได้หรือไม่ที่เลเยอร์ที่ต่างกันสามารถมีคุณสมบัติที่แตกต่างกันได้? ถ้าไม่เป็นไปได้ที่จะรับชื่อฟิลด์พร้อมกันแทนที่จะเข้าสู่ส่วนลึกนี้หรือไม่? ถ้าใช่มีวิธีใดในการเรียกชื่อฟิลด์ที่ง่ายกว่านี้ไหม

15 python  shapefile  gdal 

กำลังพยายามนำเข้า gdal สำหรับ Python 2.7.3 บน Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <mo _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swit_helper _mod = imp.load_module('_gdal', …

15 python  gdal  windows-xp 

เป็นไปได้ไหมที่จะแบ่ง shapefile ต่อคุณสมบัติใน python (ดีที่สุดจะเป็นวิธีการแก้ปัญหาที่ฉันสามารถบันทึกวัตถุเวกเตอร์ผลลัพธ์ลงในหน่วยความจำแทนลงดิสก์ได้ชั่วคราว) เหตุผล: ฉันต้องการใช้ฟังก์ชั่น gdal rasterizeLayer กับหลาย ๆ เซตย่อยของ shapefile ฟังก์ชั่นต้องการวัตถุ osgeo.ogr.Layer mkay ฉันลองไปรอบ ๆ แล้วมันสามารถทำงานได้ดังต่อไปนี้ คุณสามารถรับรูปทรงเรขาคณิตของวัตถุเลเยอร์ gdal ต่อคุณสมบัติดังต่อไปนี้ # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() ตอนนี้ฉันแค่ต้องรู้วิธีการสร้างวัตถุ osgeo.ogr.layer ตามรูปทรงเรขาคณิตนี้ สำหรับการชี้แจง ฉันต้องการฟังก์ชั่นในรหัส ogr …

15 python  gdal  rasterization  osgeo  split-by-attribute 

ฉันต้องการพัฒนาเว็บพอร์ทัลที่จะแสดงภาพที่สร้างขึ้นโดย QGIS แบบไดนามิก แต่ QGIS ที่ให้รูปแบบ TIFF มีขนาดใหญ่เกินไปและไม่สามารถแสดงบนเว็บไซต์โดยใช้โปรแกรมดูรูปภาพ นั่นเป็นเหตุผลที่ฉันตัดสินใจแปลงรูปภาพ tiff เหล่านั้นเป็น jpg / jpeg เพื่อลดขนาดเช่นเดียวกับเบราว์เซอร์ที่เป็นมิตรเพื่อให้ฉันสามารถแสดงภาพเหล่านั้นบนเว็บไซต์ของฉันได้อย่างง่ายดาย แต่ฉันไม่ทราบวิธีการแปลงภาพ TIFF เป็น jpg / jpeg โดยใช้คำสั่ง gdal โปรดช่วยฉันในเรื่องนี้

15 qgis  gdal  geotiff-tiff 

ฉันกำลังพยายามติดตั้ง GDAL 2.X บน CentOS7 โดยไม่ต้องสร้างจากแหล่งที่มาเพราะฉันต้องการกำหนดค่าในเซิร์ฟเวอร์อื่นและการติดตั้งโดยซอร์สโค้ดจะเป็นงานที่เหนื่อยล้าและใช้เวลานาน ฉันพบที่เก็บข้อมูลบางอย่างเช่น opengeo และ ELGIS แต่สำหรับ CentOS7 ใน EPEL repo ไม่มี GDAL 2.x [root@localhost ~]# yum repo-pkgs epel list gdal* (...) Available Packages gdal.x86_64 1.11.4-1.el7 @epel gdal-libs.x86_64 1.11.4-1.el7 @epel gdal-python.x86_64 1.11.4-1.el7 @epel gdal-devel.x86_64 1.11.4-1.el7 @epel gdal-doc.noarch 1.11.4-1.el7 @epel gdal-java.x86_64 1.11.4-1.el7 @epel gdal-javadoc.noarch 1.11.4-1.el7 @epel gdal-perl.x86_64 …

15 gdal  centos  repository 

ฉันต้องการรวมประมาณ 20 tiff dem ต่อแต่ละประมาณ 100MB และฉันใช้ gdalwarp (ฉันรู้เกี่ยวกับตัวเลือกไฟล์เสมือน): gdalwarp $(list_of_tiffs) merged.tiff อย่างไรก็ตามคำสั่งนี้ใช้เวลานานมากและหลังจากผ่านไป 15 นาทีก็มีการประมวลผลภาพที่ 5 จากชุด ฉันยุติมัน จากนั้นฉันก็พบว่าสคริปต์ gdal_merge.py ยังสามารถใช้สำหรับโมเสกได้เช่นกันและลองใช้งาน: gdal_merge.py $(list_of_tiffs) ซึ่งเสร็จสิ้นภายในเวลาไม่เกิน 3 นาที อย่างที่ฉันคาดไว้ว่าทั้งสองคำสั่งจะให้ผลลัพธ์เดียวกันฉันสงสัยว่าความแตกต่างระหว่างทั้งสองเป็นอย่างไรทำไม gdalwarp ใช้เวลานานมากถ้าผลลัพธ์เหมือนกัน

14 gdal  gdalwarp  gdal-merge 

ฉันพยายามที่จะเพิ่มประสิทธิภาพของฉันรหัสหลามใช้pypy ในขณะที่ฉันสามารถนำเข้าห้องสมุดส่วนใหญ่และทำงานได้สำเร็จ แต่ฉันไม่สามารถนำเข้าห้องสมุด GDAL ได้ ไลบรารีเดียวกันเมื่อนำเข้าใน Python 2.7 ทำงานได้ดี ฉันได้ลองเพิ่มตำแหน่ง GDAL ไปยัง PYTHONPATH แล้วและลองดาวน์โหลดการรวม GDAL Python ไปยังไดเรกทอรี PyPy แต่อย่างใด pypy ไม่ต้องการรู้จัก GDAL และมาพร้อมกับข้อผิดพลาดต่อไปนี้ Traceback (most recent call last): File "app_main.py", line 51, in run_toplevel File "ImportENC-P.py", line 50, in <module> from osgeo import gdal, ogr, osr File "B:\PYPY\site-packages\osgeo\__init__.py", line 21, …

14 python  gdal