Pytania otagowane jako gdalwarp

Mam obraz rastrowy geotiff, który ma układ współrzędnych o długości od 0 do 360. Poziomy środek obrazu to 180 długości geograficznej. Zobacz zdjęcie poniżej: Chcę go przekształcić w EPSG: 4326 SRS o zakresie długości -180 180. Chcę, aby środek obrazu znajdował się na południku Greenwich (0). Myślę, że ten pan …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

Próbuję wyciąć mój plik Raster zgodnie z geometrią pliku kształtu. Korzystam z następującego kodu gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Ale powoduje to powstanie czarnego koloru poza zakresem geometrii pliku kształtu. Podaję tutaj przykład. Na pierwszym obrazie, który chcę wyciąć. Drugi obraz to wynikowy raster, ale chcę trzeciego rastra, gdzie zewnętrzna …

26 gdal  gdalwarp 

Próbuję ponownie rzutować / ponownie próbkować za pomocą powiązań Pythona GDAL, ale otrzymuję nieco inne wyniki w porównaniu do tych z narzędzia wiersza poleceń gdalwarp. Zobacz krótszy przykład aktualizacji Ten skrypt ilustruje podejście do Pythona: from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair …

20 python  gdal  gdalwarp  numpy 

Po użyciu gdalwarpdo wyświetlania i wyrównywania do rastra (przez -tap) pewnej liczby rastrów zauważyłem, że wyjściowe rastry były znacznie większe niż oryginalne rastry. Dość dokładne wyszukiwanie w sieci ujawniło ten problem z Trac : Frank Warmerdam wyjaśnił przyczynę: „Po dokładnym sprawdzeniu różnica w tym pliku jest taka, że ​​gdal_translate używa …

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

Chciałem połączyć około 20 dem tiffów po około 100 MB i użyłem gdalwarp (wiem o opcji pliku wirtualnego): gdalwarp $(list_of_tiffs) merged.tiff Jednak polecenie to trwało bardzo długo i po około 15 minutach przetwarzało 5. obraz z zestawu. Kończę to. Następnie odkryłem, że skrypt gdal_merge.py może być również używany do mozaikowania …

15 gdal  gdalwarp  gdal-merge 

Piszę proste narzędzie do przycinania partii wielopasmowych plików rastrowych geotiff do tego samego (mniejszego) obszaru. Korzystając z gdalwarp, mogę z łatwością przyciąć plik za pomocą pliku kształtowego przycinającego jeden wielokąt: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif Jednak rzeczywisty obszar, do którego chcę przyciąć, będzie zawsze początkowo definiowany przez inny plik …

14 python  raster  gdal  clip  gdalwarp 

W QGIS 2.XX narzędzia rozwijanego paska menu miały na dole edytowalne okno wywoływania, w którym polecenia niedostępne w GUI można było wprowadzić ręcznie (klikając przycisk ołówka). W 3.0 narzędzia paska menu są takie same jak narzędzia GDAL znajdujące się w przyborniku przetwarzania, co ma większy sens, ale teraz nie widzę …

13 gdalwarp  qgis-3 

Próbowałem użyć kodu z QGis Save Raster jako renderowanego obrazu . ijest warstwą rastrową: 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) …

12 qgis  gdal  pyqgis  gdalwarp  gdal-translate 

Mozaikuję niektóre obrazy za pomocą gdal i chciałbym poprawić końcowy wynik, używając zanikającej / stopniowej warstwy alfa w kierunku krawędzi każdego obrazu, aby usunąć ostre krawędzie na środku mozaiki. Problem polega na tym, że część każdego pojedynczego obrazu ze stopniową warstwą alfa maskuje obrazy pod nim w końcowej mozaice, a …

12 gdal  geotiff-tiff  gdalwarp  gdal-translate  gdalbuildvrt 

Używam gdalwarp do manipulowania kafelkami SRTM w pobliżu linii danych (tj. 180 °, inaczej antimeridian). Płytki SRTM mają bardzo niewielkie (1/2 piksela) nakładanie się na południk. Możesz to zobaczyć za pomocą gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size is 1201, 1201 [...] Lower Left (-180.0004167, -16.0004167) (180d …

11 srtm  gdalwarp 

Próbuję przekształcić z Geostationary na konformizm Lambert za pomocą gdalwarp. Moje dane wejściowe są w netcdf i są we współrzędnych geograficznych (stopniach) i chciałbym wyrenderować odwzorowane dane do netcdf. Utworzyłem odpowiedni plik vrt dla wejściowych danych netcdf. Gdalwarp wyświetli plik netcdf, ale dane wyjściowe są zerami i pojawia się następujący …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

Chcę georeferencję rastra za pomocą pythoni GDAL. Moje obecne podejście polega na dzwonieniu gdal_translatei gdalwarpkorzystaniu os.systemz brzydkiej listy punktów kontroli naziemnej. Naprawdę chciałbym zrobić to w sposób natywny python. To jest obecny proces, którego używam: 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 …

10 python  raster  gdal  coordinate-system  gdalwarp 

Trochę staram się wygenerować kafelki dla obrazu w wysokiej rozdzielczości, który mamy. Obecny obraz, który mamy, to bardzo duży obraz (+ 20 GB), zapisany jako plik GeoTiff. Chciałbym wygenerować kafelki za pomocą narzędzia wiersza polecenia gdal2tiles, a następnie otworzyć i wyświetlić w Cesium, używając dostawcy obrazów TMS do dostarczenia kafelków. …

10 gdal  gdalwarp  cesium  gdal2tiles 

Próbuję wyświetlić raster. W R jest do tego projectRaster()funkcja (poniżej w pełni odtwarzalnego przykładu): # 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 

Używam Pythona i QGIS 2.0. Próbuję przyciąć rastry do folderu według jednej funkcji wielokąta. To pierwszy raz, gdy używam (powiedzmy) „PyQGIS”, wcześniej byłem przyzwyczajony do gry. Tak czy inaczej, mój prosty skrypt nie działa, każda sugestia byłaby bardzo mile widziana! 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= …

9 gdal  pyqgis  gdalwarp