last day (15 days later) » 

3:40 PM
2
A: What is the projection of my raster?

Billy34Your raster is defined with an unusual transverse mercator projection using a classic WGS84 geoid (to be quick, shape of earth). As universal transverse mercator (UTM) are defined using equator as latitude of origin and equally distributed meridians along equator, your's is using 42° north as lat...

 
how do I get latitude and longitude for every single point starting from these mercator coordinates?
 
your raster coordinates are projected that means they are on a flat surface instead of being on a spherical surface. latitude/longitude as you request are for sphere so your raster should be "unprojected" to return to spherical coordinates (lat/lon).
edited my answer with project to lat/lon
 
here is the output of the gdalwarp command: ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,1200,1400, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done.
 
I built a raster very similar to yours that displays the same information through gdalinfo and applied gdal_warp command as above and got a result raster without any errors and displaying good informations with gdalinfo. Please double check your parameters or is it your gdal version (mine is 2.2.3) or the weird character in scale_factor
try adding -s_srs="+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" along other parameters to force input spatial reference
 
This is my command (gdal_warp not found, I used gdalwarp, but I think it's the same thing): gdalwarp -s_srs='+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' C:\Users\a.frabotta\Desktop\ProtCiv\VMI\01-04-2018-00-00.tif C:\Users\a.frabotta\pippo.tif but I obtain ERROR 6 unknown option name ... My GDAL library version is 2.2.4
Hi
 
3:41 PM
you must provide other options to so same command but with added parameter -s_srs '+ ...
 
what are other options?
 
gdal_warp -s_srs='+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' -t_srs "EPSG:4326" your_raster result.tif
-s_srs is source projection, -t_srs is result projection
sorry remove space between -s_srs and '+proj...
gdal_warp -s_srs '+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' -t_srs "EPSG:4326" your_raster result.tif
 
gdalwarp -s_srs '+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' -t_srs 'EPSG:4326' C:\Users\a.frabotta\Desktop\ProtCiv\VMI\01-04-2018-00-00.tif C:\Users\a.frabotta\pippo.tif

ERROR 1: Translating source or target SRS failed:
I give up :(
 
sorry, it work for me. maybe a bogus information in your file
 
4:05 PM
I made another attempt and I get it

gdalwarp -s_srs "+proj=tmerc +lat_0=42 +lon_0=12.5 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" -t_srs "EPSG:4326" C:\Users\a.frabotta\Desktop\ProtCiv\VMI\01-04-2018-00-10.tif C:\Users\a.frabotta\pippo.tif

ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,1200,1400, skipping.
0...10...20...30...40...50...60...70...80...90...100 - done.
 

last day (15 days later) »