«gdalwarp» etiketlenmiş sorular

0'dan 360'a kadar longtitudes olan bir koordinat sistemine sahip bir geotiff raster resme sahibim. Görüntünün yatay merkezi 180 boyunda. Aşağıdaki resme bakınız: -180 180 longtitude aralığında EPSG: 4326 SRS'ye dönüştürmek istiyorum. Ve görüntünün merkezinin Greenwich meridyeninde (0) olmasını istiyorum. Sanırım bu srs çokça kullanılıyor. Sonuç şöyle görünmesini bekliyorum: Bu nedenle, …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

Raster dosyamı bir shapefile geometrisine göre kırpmaya çalışıyorum. Aşağıdaki kodu kullanıyorum gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Ancak şekil dosyası geometri kapsamı dışında siyah bir renk ortaya çıkıyor. Ben burada örnek veriyorum. Klip yapmak istediğim ilk görüntüde. İkinci görüntü elde edilen rasterdir ancak ben geometrinin dışının boş olacağı 3. Raster'ı istiyorum. …

26 gdal  gdalwarp 

Yeniden proje / GDAL python bağlamaları ile yeniden örneklemeye çalışıyorum, ama komut satırı yardımcı programına göre biraz farklı sonuçlar alıyorum gdalwarp. Daha kısa örnek için aşağıdaki güncellemeye bakın Bu komut dosyası Python yaklaşımını göstermektedir: from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair of …

20 python  gdal  gdalwarp  numpy 

Birkaç gdalwarpraster projelendirmek ve hizalamak için -tapkullandıktan sonra çıktı rasterlerinin orijinal rasterlerden önemli ölçüde daha büyük olduğunu fark ettim. Oldukça kapsamlı bir web araması bu Trac sorununu ortaya çıkardı : Frank Warmerdam bunun nedenini açıkladı: "Dikkatli bir incelemede, söz konusu dosyadaki fark, gdal_translate'in çıktı dosyasını GTiffDataset :: CreateCopy? () İçinden …

19 gdal  coordinate-system  geotiff-tiff  gdalwarp  file-size 

Her biri 100MB civarında 20 dem tiff'i birleştirmek istedim ve gdalwarp kullandım (sanal dosya seçeneğini biliyorum): gdalwarp $(list_of_tiffs) merged.tiff Ancak bu komut çok uzun zaman aldı ve 15 dakika sonra setten 5. görüntüyü işliyordu. Sonlandırıyorum. Daha sonra gdal_merge.py betiğinin mozaikleme için de kullanılabileceğini buldum ve denedim: gdal_merge.py $(list_of_tiffs) hangi daha …

16 gdal  gdalwarp  gdal-merge 

Ben aynı (daha küçük) alana çok bantlı geotiff raster dosyaları toplu kırpmak için basit bir yardımcı program yazıyorum. Gdalwarp kullanarak, tek bir çokgen kırpma şekil dosyası kullanarak bir dosyayı kolayca kırpabilirim: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif Ancak, kırpmak istediğim gerçek alan her zaman başlangıçta şekil dosyası değil, başka bir …

14 python  raster  gdal  clip  gdalwarp 

QGIS 2.XX'te, açılır menü çubuğu araçlarının altında, gui'de bulunmayan komutların manuel olarak girilebileceği (bir kurşun kalem düğmesine tıklayarak) düzenlenebilir bir arama kutusu vardı. 3.0'da, menü çubuğu araçları, işleme araç kutusunda bulunan GDAL araçlarıyla aynıdır, bu da daha mantıklıdır, ancak şimdi bir araç işlemine basit komutlar eklemenin bir yolunu görmüyorum. warp …

13 gdalwarp  qgis-3 

QGis Save Raster kodunu Rendered Image olarak kullanmaya çalıştım . ibir raster katmanıdır: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs = pipelayer.crs().toWkt() pipe = QgsRasterPipe() pipe.set(pipeprovider.clone()) pipe.set(piperenderer.clone()) pipedFile = os.path.join(tempfile.gettempdir(), safeLayerName + '_pipe.tif') print pipedFile file_writer = QgsRasterFileWriter(pipedFile) file_writer.writeRaster(pipe, …

12 qgis  gdal  pyqgis  gdalwarp  gdal-translate 

Bazı görüntüleri gdal ile mozaikliyorum ve mozaiğin ortasındaki keskin kenarları kaldırmak için her görüntünün kenarına doğru solma / kademeli bir alfa katmanı kullanarak nihai sonucu iyileştirmek istiyorum. Yaşadığım sorun, her bir görüntünün aşamalı alfa katmanıyla bölümünün, aşağıda gösterildiği gibi yarı saydam olmak yerine, altındaki mozaikin altındaki görüntüleri maskelemesi: İdeal olarak, …

12 gdal  geotiff-tiff  gdalwarp  gdal-translate  gdalbuildvrt 

Ben dateline yakınındaki SRTM fayans manipüle etmek için gdalwarp kullanıyorum (yani 180 °, yani antimeridian). SRTM karolarının meridyenle çok hafif (1/2 piksel) çakışması vardır. Bunu gdalinfo kullanarak görebilirsiniz: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size is 1201, 1201 [...] Lower Left (-180.0004167, -16.0004167) (180d 0' 1.50"W, 16d 0' …

11 srtm  gdalwarp 

Gdalwarp kullanarak Geostationary'den Lambert'e uyumlu olmaya çalışıyorum. Giriş verilerim netcdf ve coğrafi koordinatlarda (derece) ve yeniden eşlenen verileri netcdf'ye çıkarmak istiyorum. Giriş netcdf verileri için karşılık gelen bir vrt dosyası oluşturduk. Gdalwarp netcdf dosyasını çıktılar, ancak çıktı verilerinin tümü sıfırdır ve aşağıdaki hatayı alıyorum: Creating output file that is 5120P …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

Ben kullanarak raster coğrafi olarak tanımlanmasına istiyorum pythonve GDAL. Benim şu anki yaklaşımı çağrısına olduğunu gdal_translateve gdalwarpkullanma os.systemve yer kontrol noktalarının çirkin listesinde. Bunu yerel olarak yapmanın bir yolunu gerçekten istiyorum python. Bu şu anda kullandığım işlem: import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 …

10 python  raster  gdal  coordinate-system  gdalwarp 

Sahip olduğumuz yüksek çözünürlüklü bir görüntü için fayans üretmek için biraz uğraşıyorum. Şu anki görüntümüz GeoTiff dosyası olarak kaydedilen çok büyük (+ 20 GB) bir görüntüdür. Gdal2tiles komut satırı yardımcı programını kullanarak fayans oluşturmak ve daha sonra fayans sağlamak için TMS görüntü sağlayıcı kullanarak Cesium, açmak ve görüntülemek istiyorum. Gdalinfo …

10 gdal  gdalwarp  cesium  gdal2tiles 

Bir Raster yansıtmaya çalışıyorum. R'de bunun projectRaster()işlevi vardır (tamamen tekrarlanabilir bir örneğin altında): # 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 = …

9 raster  gdal  r  coordinate-system  gdalwarp 

Python ve QGIS 2.0 kullanıyorum. Bir çokgen özelliği ile bir klasörde rasters klip çalışıyorum. "PyQGIS" i (ilk diyelim) kullanarak benim için ilk kez, daha önce kemirmeye alışkındım. Neyse, benim basit senaryomun çalışmasını alamıyorum, herhangi bir öneri çok takdir edilecektir! import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() CLIP= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER="C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00" OUTPUT= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/foscagno_pyqgis/" …

9 gdal  pyqgis  gdalwarp