Ho un'immagine raster geotiff che ha un sistema di coordinate con lunghezze da 0 a 360. Il centro orizzontale dell'immagine è di 180 longitudini. Vedi l'immagine qui sotto: Voglio trasformarlo in EPSG: 4326 SRS con intervallo di longitudine -180 180. E voglio che il centro dell'immagine sia al meridiano di …
Sto provando a ritagliare il mio file Raster secondo una geometria shapefile. Sto usando il seguente codice gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Ma sta ottenendo un colore nero al di fuori dell'estensione della geometria dello shapefile. Sto dando l'esempio qui. Nella prima immagine che voglio ritagliare. La seconda immagine è …
Sto cercando di riproiettare / ricampionare con i collegamenti python GDAL, ma sto ottenendo risultati leggermente diversi rispetto a quelli dell'utilità della riga di comando gdalwarp. Vedi aggiornamento di seguito per un esempio più breve Questo script illustra l'approccio Python: from osgeo import osr, gdal import numpy def reproject_point(point, srs, …
Dopo aver utilizzato gdalwarpper proiettare e allineare alla griglia (via -tap) un numero di raster ho notato che i raster di output erano significativamente più grandi dei raster originali. Una ricerca web abbastanza approfondita ha rivelato questo problema di Trac : Frank Warmerdam ha spiegato il motivo: "Su un'attenta revisione, …
Volevo unire circa 20 dem tiffs ciascuno intorno a 100 MB e ho usato gdalwarp (conosco l'opzione del file virtuale): gdalwarp $(list_of_tiffs) merged.tiff Tuttavia, questo comando ha richiesto molto tempo e dopo circa 15 minuti stava elaborando la quinta immagine dal set. Lo interrompo. Ho quindi scoperto che lo script …
Sto scrivendo una semplice utility per ritagliare batch di file raster geotiff multibanda nella stessa area (più piccola). Usando gdalwarp, posso facilmente ritagliare un file usando uno shapefile di ritaglio a singolo poligono: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif Tuttavia, l'area effettiva a cui voglio tagliare sarà sempre inizialmente definita …
In QGIS 2.XX, gli strumenti della barra dei menu a discesa avevano una casella di chiamata modificabile nella parte inferiore, in cui i comandi non disponibili nella GUI potevano essere inseriti manualmente (facendo clic su un pulsante a forma di matita). In 3.0, gli strumenti della barra dei menu sono …
Sto creando un mosaico di alcune immagini con gdal e vorrei migliorare il risultato finale usando uno strato alfa sfumato / graduale verso il bordo di ogni immagine per rimuovere i bordi nitidi nel mezzo del mosaico. Il problema che sto riscontrando è che la porzione di ogni singola immagine …
Sto usando gdalwarp per manipolare le tessere SRTM vicino alla linea di dati (ovvero 180 °, ovvero l'antimeridiano). Le tessere SRTM hanno una leggera sovrapposizione (1/2 pixel) con il meridiano. Puoi vederlo usando gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size is 1201, 1201 [...] Lower Left (-180.0004167, …
Sto cercando di rimappare da Geostationary a Lambert conforme usando gdalwarp. I miei dati di input sono in netcdf, ed è in coordinate geografiche (gradi) e vorrei inviare i dati rimappati su netcdf. Ho creato un file vrt corrispondente per i dati di input netcdf. Gdalwarp produrrà il file netcdf, …
Voglio georeferenziare un raster usando pythone GDAL. Il mio attuale approccio è chiamare gdal_translatee gdalwarpusare os.systeme un brutto elenco di punti di controllo a terra. Mi piacerebbe davvero un modo per farlo nativamente all'interno python. Questo è l'attuale processo che sto usando: import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 …
Ho faticato un po 'a generare tessere per un'immagine ad alta risoluzione che abbiamo. L'immagine attuale che abbiamo è un'immagine molto grande (+ 20 GB), salvata come file GeoTiff. Vorrei generare i riquadri utilizzando l'utilità della riga di comando gdal2tiles e quindi aprirli e visualizzarli in Cesio, utilizzando il provider …
Sto cercando di proiettare un raster. In R c'è la projectRaster()funzione per fare ciò (sotto un esempio completamente riproducibile): # 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 <- …
Sto usando Python e QGIS 2.0. Sto cercando di ritagliare i raster in una cartella con una funzione poligonale. È la prima volta che uso (diciamo) "PyQGIS", prima ero abituato ad arcipare. Comunque, non riesco a far funzionare la mia semplice sceneggiatura, qualsiasi suggerimento sarebbe molto apprezzato! import qgis.core, qgis,utils …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.