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

[Stéphane]

So, as far as I understand:

GDAL can reference image with geos projection like
       +proj=geos +lon_0=42 +h=35785831 +datum=WGS84

GDAL can put GeoTransformation information like
       -a_ullr -924124.2 5421728.5 4179561.6 3468466.1
I guess this is related to maxium / minimum latitude and longitude seen from
satellite position.


*But I'm still wondering if someone success to reproject an image from
geostationnary projection to regular lat/lon ?*

Thanks,
Stéphane

2010/3/2 Frank Warmerdam <warmerdam at pobox.com>

> 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<http://pobox.com/%7Ewarmerdam>
> and watch the world go round - Rush    | Geospatial Programmer for Rent
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20100305/5bd1005d/attachment.html