Pytania oznaczone «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)....

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ąć....

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...

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...

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...

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...

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...

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...

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 =...

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:...

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...

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...

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ż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,...

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 <-...