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: [email protected] [[email protected]] on behalf of Johannes Landmann [[email protected]] Sent: 25 February 2013 20:27 To: [email protected] 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 [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
