[gdal-dev] Re: The GEOS projection and gdalwarp

[Marcel Derrien]

Hi,
First I suggest you use gdal_translate and provide -a_ullr ulx uly lrx 
lry so that your input file is georeferenced.
Assuming that your input file is the full earth disk 3712x371,  the 
corner coordinates (metre) of your image are
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 input.tif 
input_geo.tif
Using input_geo.tif with gdalwarp should allow your desired reprojection.
Marcel
---
Marcel Derrien,
http://www.meteo-spatiale.fr

Date: Mon, 26 Mar 2012 13:31:24 +0200
From: "Simon R. Proud" <Srp at geo.ku.dk>
Subject: [gdal-dev] The GEOS projection and gdalwarp
To: <gdal-dev at lists.osgeo.org>
Message-ID: <4F706FAC.9690.0040.1 at geo.ku.dk>
Content-Type: text/plain; charset=US-ASCII

Hi,
I'm working with meteosat data in the GEOS projection and have been having some problems with gdalwarp.The data I'm using contains no geo-information, so I call gdalwarp with this:

gdalwarp -s_srs '+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0' -t_srs '+proj=latlong +datum=WGS84' -ts 2000 2200 -te -20.5 -40.99992 52.2272 39.0 /media/Data/ISOP.tif /home/processing/GEO.tif

This results in:
"ERROR1: Unable to compute a transformation between pixel/line and georeferenced coordinates for /media/Data/ISOP.tif.
There is no affine transformation and no GCPs."

Do you know how I could overcome this? When I convert *to* GEOS I needed to include some additional info (such as subsatellite pixel size as an argument with -tr) but I can't find how to include that info in the source projection. Any ideas/hints greatly appreciated!

Thanks,
Simon.

~~~~~~~~~~~~~~~~~~~~~~~~
Dr Simon Proud
Post-doc, Satellite Remote Sensing
University of Copenhagen

Tel: +45 35 32 25 84
~~~~~~~~~~~~~~~~~~~~~~~~




-------------- next part --------------
A non-text attachment was scrubbed...
Name: marcel_derrien.vcf
Type: text/x-vcard
Size: 375 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120327/28fdd783/marcel_derrien.vcf