$> echo "16E 45N"| proj -v +init="epsg:32633"
$> echo "16E 45N"| proj -v +proj=utm +zone=33 +datum=WGS84 +units=m +no_defs
Nota: è possibile sostituire il comando invproj con il seguente: proj -i
$> echo "578815.30 4983436.77"| invproj +init="epsg:32633"
$> echo "578815.30 4983436.77"| proj -I +init="epsg:32633"
$> echo "16E 45N"| proj -v -f "%.1f" +init="epsg:32633"
$> echo "578814 4983432"| proj -I -f "%.5f" +init="epsg:32633"
km 1000. Kilometer
m 1. Meter
kmi 1852.0 International Nautical Mile
in 0.0254 International Inch
ft 0.3048 International Foot
us-in 1./39.37 U.S. Surveyor's Inch
us-ft 0.304800609601219 U.S. Surveyor's Foot
[---]
Opzione | Significato |
+a | definizione del semiasse maggiore |
+b | definizione del semiasse minore |
+datum | nome del datum (proj -ld) |
+ellps | nome dell'ellissoide (proj -le) |
+k / +k_0 | fattore di scala |
+lat_0 | latitudine dell'origine |
+lon_0 | meridiano centrale |
Opzione | Significato |
+towgs84 | parametri (3 o 7) per il datum shifting (trasformazione di Helmert) |
+nadgrids | file in formato ntv1/ntv2, per la trasformazione mediante reticolo |
+no_defs | inibisce l'uso del file /usr/share/proj/proj_def.dat per i parametri di default |
+proj | nome della proiezione da utilizzare (proj -lp) |
+zone | definizione della zona UTM |
+x_0 | falso est |
+y_0 | falso nord |
+units | fattore di conversione da unità cartesiane a unità metriche |
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
+proj=longlat +ellps=intl +towgs=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +no_defs
+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
+proj=tmerc +towgs=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +lat_0=0 +lon_0=15 +k=0.9996 +x_0=2520000 +y_0=0 +ellps=intl +units=m +no_defs
Nota: i sistemi di riferimento specificati possono essere entrambi geografici, entrambi cartografici o misti
echo “2520000 456789” | cs2cs +init=”epsg:3004”
echo “16 45” | cs2cs +init=”epsg:4230” +to +proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
echo “16 45” | cs2cs +init=”epsg:4230” +to +init=”epsg:32633”
echo “16 45” | cs2cs +proj=latlong +ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +no_defs +to +init=”epsg:32633”
gdalinfo –formats
gdalinfo –format GTiff
gdalinfo file_raster
-noct: sopprime l’outpput delle color tables
-stats: calcola e mostra le statistiche
gdal_translate [opzioni] src_dataset dst_dataset
Opzione | Significato |
-of | formato di output |
-ot | data type della banda esportata |
-b | banda/e da esportare |
-a_srs | sovrascrive le informazioni sul datum |
-outsize | rescaling dell'output |
-scale | resampling dell'output |
-srcwin | subsetting |
-mo | aggiunta di metatag all'output |
-co | opzioni di creazione specifiche del formato (v.di gdalinfo -format format_name) |
gdal_translate -of JPEG data.tiff data.jpg
gdal_translate -of JPEG -co “PROGRESSIVE=TRUE” -co “QUALITY=90” -mo “PIPPO=DATA” -a_srs epsg:3004 data.tif data.jpg
gdalwarp [opzioni] src_dataset* dst_dataset
-s_srs | source SRS |
-t_srs | destination SRS |
-order | ordine del polinomio interpolatore (da 1 a 3) |
-te | estensione geografica dell'output (extension) |
-tr | risoluzione dell'output (in unità di misura del SRS target) |
-r | metodo di ricampionamento (near, bilinear, cubic, cubicspline, lanczos) |
-multi | modalità di processamento multicore (parallela) |
-of | formato dell'output |
-co | opzioni di creazione dell'output (specifiche del formato) |
gdalwarp -t_srs epsg:3004 -multi -order cubic srtm_40_04.tif srtm_40_05.tif dem.tif
gdaladdo [opzioni] src_data livelli
gdaladdo -r cubic data.tif 2 4 6 8 16
gdaladdo -r cubic -ro –config COMPRESS_OVERVIEW JPEG data.tif 2 4 6 8 16
Nota: ogni modalità presenta una serie di opzioni specifiche per modificare il comportamento del programma
gdaldem hillshade srcfile dstfile
-z | fattore di scala per le altitudini (asse z) |
-s | rapporto tra unità verticali e orizzontali (utile se le latitudini sono espresse in gradi e le longitudini in metri) |
-az | azimut della fonte luce (in gradi). Il valore 0 corrisponde al lato superiore del raster, 90 alla direzione EST |
-alt | altitudine della fonte luce, in gradi |
gdaldem color-relief srcfile colorsfile dstfile
La sintassi del file di configurazione dei colori è estremamente semplice, e mutuata dal comando grass r.colors
Si tratta di un file di testo in cui ciascuna riga è costituita da quattro valori separati da spazio: (altitudine R G B). Un quinto valore può essere aggiunto per caratterizzare la componente di opacità del pixel in uscita
Il valore speciale nv matcha le celle del raster con valore NULL
Nelle ultime versioni di gdal, le curve di livello calcolate sono orientate, per garantire una consistenza topologica
gdal_contour [opzioni] srcfile dstfile
-a | imposta il nome della colonna dei valori altimetrici |
-f | formato di output (di default shp) |
-i | inposta l'intervallo delle isoipse |
-fl | estrae solo le curve all'altitudine selezionata |
gdal_rasterize [opzioni] srcfile* dstfile
-b | definizione delle bande su cui salvare i valori numerici |
-burn | valori fissi da utilizzare nel processo di rasterizzazione |
-a | nome dell'attributo del file vettoriale da cui prendere i valori numerici |
-where | clausula WHERE per la selezione delle feature |
-sql | stringa SQL per la selezione delle feature |
-of | formato del file di output |
-a_srs | sovrascrive l'SRS del file in output |
ogrinfo [opzioni] srcfile
ogrinfo [opzioni] srcfile layername
-al | nella prima modalità di utilizzo, mostra tutte le geometrie di tutti i layers |
-so | nella seconda modalità di utilizzo, mostra solo i metadati e non le geometrie del layer |
-spat | seleziona solo le geometrie che intersecano il rettangolo specificato |
-where | seleziona le geometrie da visualizzare in base alla clausola SQL WHERE specificata |
-sql | seleziona le geometrie da visualizzare in base ad un'espressione SQL |
ogr2ogr [opzioni] dstfile srcfile [opzioni layer] [layer [layer …]]
-f | specifica il formato di output |
-a_srs | sovrascrive l'SRS del file in output |
-s_srs | specifica l'SRS del file in input |
-t_srs | specifica l'SRS del file in output |
-where | seleziona le geometrie in base ad una clausola WHERE SQL |
-sql | seleziona le geometrie in base ad una espressione SQL |
-clipsrc | effettua il clipping delle geometrie da selezionare in base a diversi parametri (ad es. un rettangolo) |
-append | aggiunge geometrie ad un layer/datasource già esistente |
-segmentize | effettua la segmentazione delle geometrie rispetto ad una distanza massima |
Maggiori informazioni sono reperibili all’URL: https://lpdaac.usgs.gov/lpdaac/products