[gdal-dev] Georeference EUMETSAT's NDVI (hdf5)

[Simon Richard Proud]

Hi,
I use the following to warp SEVIRI images:

ulx_msg=-5570248.832537
uly_msg=5570248.832537
lrx_msg=5567248.429179
lry_msg=-5567248.429179
gdal_translate -a_srs  "+proj=geos +a=6378169 +b=6356583.8 +lon_0=0 +h=35785831" -a_ullr $ulx_msg $uly_msg $lrx_msg $lry_msg (infile) (outfile)
gdalwarp -s_srs '+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0' -t_srs '+proj=latlong +datum=WGS84' -tr 0.04 0.04 -te -40 -40 40 40 -order 3 (infile) (outfile)

The correct parameters for a, b, lon_0 and h can be found from the weekly MSG position reports, the above values were valid for MSG-2 in late December. Change -tr and -te as needed for whatever geographical area you're interested in.
Hope that helps,
Simon
________________________________________
From: gdal-dev-bounces at lists.osgeo.org [gdal-dev-bounces at lists.osgeo.org] on behalf of Johannes Landmann [j_landmann at gmx.de]
Sent: 25 February 2013 20:27
To: gdal-dev at lists.osgeo.org
Subject: [gdal-dev] Georeference EUMETSAT's NDVI (hdf5)

Dear GDAL-friends,

I am quite new to GDAL and I've got a substantial question:

For my bachelor thesis I work with EUMETSAT's NDVI data published in HDF5 format. These data are "produced" in the geostationary satellite projection but do not contain geospatial reference in the HDF5 dataset.

I have already tried out to assign the projection to the images with  "gdalwarp -s_srs SR-ORG:81" as well as "gdalwarp -t_srs SR-ORG:81" but with no result.

What am I doing wrong? Is there another (better) way to "reassign" the projection?

Thanks in advance,

Johannes
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev