I have downloaded the GMTED2010 elevation data from
http://topotools.cr.usgs.gov/gmted_viewer/viewer.htm
I am merging all the files/tiles (to create a global dataset) with:
gdal_merge.py
There is one band in the tif files and I would like to convert it to meters!
I am converting the resulting file from tif format to netCDF format for my software with:
gdal_translate -of netCDF 1.tif 1.nc
The netCDF-file contains elevation in Band 1. I am unsure how to convert this to meters. Currently it seems (eyeballs norm) a factor 2 smaller than the actual elevation.
ncump -h on the netCDF file yields:
netcdf mea {dimensions: lon = 43200 ; lat = 21600 ;variables: char crs ; crs:grid_mapping_name = "latitude_longitude" ; crs:longitude_of_prime_meridian = 0. ; crs:semi_major_axis = 6378137. ; crs:inverse_flattening = 298.257223563 ; crs:spatial_ref = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]]" ; crs:GeoTransform = "-180.0001388888889 0.008333333333333333 0 89.99986111111112 0 -0.008333333333333333 " ; double lat(lat) ; lat:standard_name = "latitude" ; lat:long_name = "latitude" ; lat:units = "degrees_north" ; double lon(lon) ; lon:standard_name = "longitude" ; lon:long_name = "longitude" ; lon:units = "degrees_east" ; short Band1(lat, lon) ; Band1:long_name = "GDAL Band Number 1" ; Band1:_FillValue = -32767s ; Band1:grid_mapping = "crs" ;
// global attributes: :GDAL_AREA_OR_POINT = "Area" ; :Conventions = "CF-1.5" ; :GDAL = "GDAL 2.0.0, released 2015/06/14" ; :history = "Fri Jul 24 15:33:07 2015: GDAL CreateCopy( mea.nc, ... )" ;}
Any help would be much appreciated.
أكثر...
http://topotools.cr.usgs.gov/gmted_viewer/viewer.htm
I am merging all the files/tiles (to create a global dataset) with:
gdal_merge.py
There is one band in the tif files and I would like to convert it to meters!
I am converting the resulting file from tif format to netCDF format for my software with:
gdal_translate -of netCDF 1.tif 1.nc
The netCDF-file contains elevation in Band 1. I am unsure how to convert this to meters. Currently it seems (eyeballs norm) a factor 2 smaller than the actual elevation.
ncump -h on the netCDF file yields:
netcdf mea {dimensions: lon = 43200 ; lat = 21600 ;variables: char crs ; crs:grid_mapping_name = "latitude_longitude" ; crs:longitude_of_prime_meridian = 0. ; crs:semi_major_axis = 6378137. ; crs:inverse_flattening = 298.257223563 ; crs:spatial_ref = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]]" ; crs:GeoTransform = "-180.0001388888889 0.008333333333333333 0 89.99986111111112 0 -0.008333333333333333 " ; double lat(lat) ; lat:standard_name = "latitude" ; lat:long_name = "latitude" ; lat:units = "degrees_north" ; double lon(lon) ; lon:standard_name = "longitude" ; lon:long_name = "longitude" ; lon:units = "degrees_east" ; short Band1(lat, lon) ; Band1:long_name = "GDAL Band Number 1" ; Band1:_FillValue = -32767s ; Band1:grid_mapping = "crs" ;
// global attributes: :GDAL_AREA_OR_POINT = "Area" ; :Conventions = "CF-1.5" ; :GDAL = "GDAL 2.0.0, released 2015/06/14" ; :history = "Fri Jul 24 15:33:07 2015: GDAL CreateCopy( mea.nc, ... )" ;}
Any help would be much appreciated.
أكثر...