[gdal-dev] Georeferencing Geostationnary (+proj=geos)

[Frank Warmerdam]

Stéphane wrote:
> Hi,
> 
> Does anyone know if Gdal, current version or development version, is 
> able to translate from Geostationnary projection to other regular 
> projection ?
> 
> Is it planned in future Gdal release ?

Stéphane,

I can confirm that GDAL/OGR supports the geostationary projection
(+proj=geos in PROJ.4).  You can specify the +proj=geos projection
strings with GDAL where it would normally accept proj.4 strings.
I see it translates:

+proj=geos +lon_0=42 +h=35785831 +datum=WGS84

as:

PROJCS["unnamed",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.257223563,
                 AUTHORITY["EPSG","7030"]],
             TOWGS84[0,0,0,0,0,0,0],
             AUTHORITY["EPSG","6326"]],
         PRIMEM["Greenwich",0,
             AUTHORITY["EPSG","8901"]],
         UNIT["degree",0.0174532925199433,
             AUTHORITY["EPSG","9108"]],
         AUTHORITY["EPSG","4326"]],
     PROJECTION["Geostationary_Satellite"],
     PARAMETER["central_meridian",42],
     PARAMETER["satellite_height",35785831],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0]]

I do not know how you would construct the geotransform
for an arbitrary satellite image though.

Best regards,

> 
> 2010/3/2 António Rocha <antonio.rocha at deimos.com.pt 
> <mailto:antonio.rocha at deimos.com.pt>>
> 
>     Hi Stéphane
> 
>     According to http://www.gdal.org/gdal_translate.html
> 
>     *-a_ullr* /ulx uly lrx lry/:
>       Assign/override the georeferenced bounds of the output file. This
>       assigns georeferenced bounds to the output file, ignoring what would
>       have been derived from the source file.
>     Yes... Regarding extracting sub-domains, I didn't tried
> 
>     Yes, please keep me in touch in case there are better methods
> 
>     Antonio
> 
> 
>     Stéphane wrote:
> 
>         Thanks for feedback !
> 
>         Do you have any ideas of the meaning of -a_ullr arguments they
>         gave you ?
> 
>         Futhermore, once you get your GeoTiff file, did you success to
>         translate it to other projection, or extracting sub domains ?
> 
>         I've heard that GEOS projection is not so really supported
>         currently by GDAL, and re-projection / sub domain extraction
>         doesn't work correctly ...
> 
>         I'll keep you in touch if any news !
>         Stéphane
> 
>         2010/3/2 António Rocha <antonio.rocha at deimos.com.pt
>         <mailto:antonio.rocha at deimos.com.pt>
>         <mailto:antonio.rocha at deimos.com.pt
>         <mailto:antonio.rocha at deimos.com.pt>>>
> 
> 
>            Greetings
> 
>            I had exacly the same problem. I have been using LSASAF data
>         (from
>            MSG-2nd Generation satellites) that provides LST and FVC in HDF5
>            files. For some technical reason they are not goerreferenced and
>            their  information.
> 
>            I talked with a few ILWIS software developers and they gave the
>            following Commands (adapted to my spatial coverage images):
> 
>            gdal_translate -a_srs "+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0
>            +y_0=0.0" -a_ullr -924124.2 5421728.5 4179561.6 3468466.1
>            "HDF5:"HDF5_LSASAF_MSG_LST_Euro_200802170100"://LST"
>            "euro_area_georeferenced.tif"
> 
> 
>            gdalwarp -s_srs "+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0
>            +y_0=0.0" -t_srs "+proj=latlong" "euro_area_georeferenced.tif"
>            "euro_area_latlong.tif"
> 
>            My HDF5 file doesn't contain any infoprmation regarding GEOS
>            projection. So, I had to do this
> 
> 
>            Hope this helps
> 
>            If you get any news/information/method regarding this, please let
>            me know
> 
>            Best regards
>            Antonio
> 
> 
>            __________ Information from ESET NOD32 Antivirus, version of
>         virus
>            signature database 4907 (20100302) __________
> 
>            The message was checked by ESET NOD32 Antivirus.
> 
>            http://www.eset.com
> 
> 
> 
> 
> 
> 
>     __________ Information from ESET NOD32 Antivirus, version of virus
>     signature database 4908 (20100302) __________
> 
> 
>     The message was checked by ESET NOD32 Antivirus.
> 
>     http://www.eset.com
> 
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev


-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent