Hi:
I have some hdf5 data and I want to change the subdataset to geotiff, and I
use gdalinfo to get the subdataset name and then I use "gdal_translate
subdatasetname result.tif".
However I found that the crs of the generated tif is not recognized by the
arcgis,then I want to know why so I create a new tif by arcgis(here it is
named arcgis.tif,and the coordinate system is wgs84). Then I use gdalinfo
command to get the metadata of the arcgis.tif and compare it with the tif
generated by gdal(here I call it gdal.tif). Then I found that the Coordinate
are the same as the following:
--------------------------
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
------------------------------
the difference is the *Corner Coordinates*:
the arcgis.tif is :
-------------------
Corner Coordinates:
Upper Left ( 103.8438587, 28.6769083) (103d50'37.89"E, 28d40'36.87"N)
Lower Left ( 103.8438587, 28.6668077) (103d50'37.89"E, 28d40'0.51"N)
Upper Right ( 103.8546032, 28.6769083) (103d51'16.57"E, 28d40'36.87"N)
Lower Right ( 103.8546032, 28.6668077) (103d51'16.57"E, 28d40'0.51"N)
Center ( 103.8492310, 28.6718580) (103d50'57.23"E, 28d40'18.69"N)
------------------
the gdal.tif is :
----------------
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 720.0)
Upper Right ( 1440.0, 0.0)
Lower Right ( 1440.0, 720.0)
Center ( 720.0, 360.0)
------------------------
It seems that the bbox is defined by different manner , the arcgis.tif is by
Degrees minutes and seconds and the gdal.tif is by grid of some other?
I wonder how to solve this problem?
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
