Daha çox

GDAL istifadə edərək TIFF-dən GeoTIFF yaradılır?


Mənim bir tiff obrazım var və onu daha çox istəmək istərdim, buna görə də bir Geotiff yaradın. Bunu etmək üçün GDAL istifadə etməliyəm. Xüsusilə QGIS sənədindəki yürütülebilir mağazadan istifadə edirəm (gdalwarp.exe, gdaltranslate.exe ...).

Hələlik bu kodu bir DOS əmrində (komanda sətri) istifadə edirəm və bu işə yaramır.

gdal_translate -of GTiff -gcp 33.9983 669.668 -833589 -692.339 -gcp 124.661 1377.97 -833520 -1379.2 -gcp 6788.34 -10.2993 -826876 6.06375 -gcp 6839.33 748.997 -833363336336336336336336636833683368336335 -1021.34 "D: / Users / (…) / Desktop/image_projetee_UTM.tiff" "C: / Users / (…) /AppData/Local/Temp/image_projetee_UTM.tiff" gdalwarp -r-order yaxın 1 -co COMPRESS = YOX -dstalpha "C: / Users / (…) /AppData/Local/Temp/image_projetee_UTM.tiff" "D: / Users / (…) / Desktop/image_projetee_UTM_georef.tif"

1 - Bu kodu işlədikdən sonra çıxış görüntüsü GTiff deyil, giriş şəkli kimi bir Tiffdir.

2 - Məsələn, Google Earth-də çıxış görüntüsünü açmağa çalışsam, yerdən istinad edilmir

Düşünürəm ki, bir Tiff sənədinə yer ayırmaq üçün bir şey buraxmışam ...

Kimsə mənə kömək edə bilər?


Beləliklə, 4 gps / s (yer nəzarət nöqtələri) istifadə edərək GDAL ilə bir tif faylını georeferanslaşdırmağı bacarıram. Bu qınağı etmək üçün komanda sətirində gdal istifadə edirəm.

Əvvəlcə belə tərcümə ilə gdal istifadə edin:

gdal_translate -of GTiff -gcp 0 0 -6.848326 45.501053 -gcp 6862 0 -6.490975 45.501503 -gcp 0 1379 -6.762872 45.377363 -gcp 6862 1379 -6.545354 45.382523 "Inputimage.tif" "Çıxış

Torpaq nəzarət nöqtələri belə qurulur: görüntüdəki piksel koordinatı (x, y) və sonra coğrafi yerləşmə (uzunluq, enlik).

Gdaltranslate-dən sonra çıxan faylın gcps ilə düzgün bir şəkildə bəsləndiyini gdalinfo ilə yoxlamağı təklif edirəm. Gdalinfo ilə dos əmrində belə bir şey görməlisiniz:

Koordinat Sistemi "GCP Proyeksiyası = GCP [0]: Id = 1, Məlumat = (0,0) -> (45.501053, -6.848326,0) GCP [1]: Id = 2, Məlumat = (6862,0) - > (45.501503, -6.490975,0) GCP [2]: Id = 3, Məlumat = (0,1379) -> (45.377363, -6.762872,0) GCP [3]: Id = 4, Məlumat = (6862,1379) ) -> (45.382523, -6.545354,0)

Bu nöqtədə görüntü hələ yerdən istinad edilmir. Gdalwarp'u belə istifadə etməlisən (qınamağın bir neçə yolu var, burada "yaxın" dan istifadə edirəm:

gdalwarp -r sifariş yaxınlığında 1 -co COMPRESS = NONE -dstalpha "OutputImageF_From_gdal_translate.tif" "FinalImage.tif"


Videoya baxın: Read and write raster files with GDAL in Python (Oktyabr 2021).