So, as far as I understand:<br><br>GDAL can reference image with geos projection like<br> +proj=geos +lon_0=42 +h=35785831 +datum=WGS84<br><br>GDAL can put GeoTransformation information like<br> -a_ullr -924124.2 5421728.5 4179561.6 3468466.1<br>
I guess this is related to maxium / minimum latitude and longitude seen from satellite position.<br><br><br><b>But I'm still wondering if someone success to reproject an image from geostationnary projection to regular lat/lon ?</b><br>
<br>Thanks,<br>Stéphane<br><br><div class="gmail_quote">2010/3/2 Frank Warmerdam <span dir="ltr"><<a href="mailto:warmerdam@pobox.com">warmerdam@pobox.com</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">Stéphane wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
Does anyone know if Gdal, current version or development version, is able to translate from Geostationnary projection to other regular projection ?<br>
<br>
Is it planned in future Gdal release ?<br>
</blockquote>
<br></div>
Stéphane,<br>
<br>
I can confirm that GDAL/OGR supports the geostationary projection<br>
(+proj=geos in PROJ.4). You can specify the +proj=geos projection<br>
strings with GDAL where it would normally accept proj.4 strings.<br>
I see it translates:<br>
<br>
+proj=geos +lon_0=42 +h=35785831 +datum=WGS84<br>
<br>
as:<br>
<br>
PROJCS["unnamed",<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
TOWGS84[0,0,0,0,0,0,0],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.0174532925199433,<br>
AUTHORITY["EPSG","9108"]],<br>
AUTHORITY["EPSG","4326"]],<br>
PROJECTION["Geostationary_Satellite"],<br>
PARAMETER["central_meridian",42],<br>
PARAMETER["satellite_height",35785831],<br>
PARAMETER["false_easting",0],<br>
PARAMETER["false_northing",0]]<br>
<br>
I do not know how you would construct the geotransform<br>
for an arbitrary satellite image though.<br>
<br>
Best regards,<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="h5">
<br>
2010/3/2 António Rocha <<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a> <mailto:<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a>>><br>
<br>
Hi Stéphane<br>
<br>
According to <a href="http://www.gdal.org/gdal_translate.html" target="_blank">http://www.gdal.org/gdal_translate.html</a><br>
<br>
*-a_ullr* /ulx uly lrx lry/:<br>
Assign/override the georeferenced bounds of the output file. This<br>
assigns georeferenced bounds to the output file, ignoring what would<br>
have been derived from the source file.<br>
Yes... Regarding extracting sub-domains, I didn't tried<br>
<br>
Yes, please keep me in touch in case there are better methods<br>
<br>
Antonio<br>
<br>
<br>
Stéphane wrote:<br>
<br>
Thanks for feedback !<br>
<br>
Do you have any ideas of the meaning of -a_ullr arguments they<br>
gave you ?<br>
<br>
Futhermore, once you get your GeoTiff file, did you success to<br>
translate it to other projection, or extracting sub domains ?<br>
<br>
I've heard that GEOS projection is not so really supported<br>
currently by GDAL, and re-projection / sub domain extraction<br>
doesn't work correctly ...<br>
<br>
I'll keep you in touch if any news !<br>
Stéphane<br>
<br>
2010/3/2 António Rocha <<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a><br>
<mailto:<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a>><br>
<mailto:<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a><br>
<mailto:<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a>>>><br>
<br>
<br>
Greetings<br>
<br>
I had exacly the same problem. I have been using LSASAF data<br>
(from<br>
MSG-2nd Generation satellites) that provides LST and FVC in HDF5<br>
files. For some technical reason they are not goerreferenced and<br>
their information.<br>
<br>
I talked with a few ILWIS software developers and they gave the<br>
following Commands (adapted to my spatial coverage images):<br>
<br>
gdal_translate -a_srs "+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0<br>
+y_0=0.0" -a_ullr -924124.2 5421728.5 4179561.6 3468466.1<br>
"HDF5:"HDF5_LSASAF_MSG_LST_Euro_200802170100"://LST"<br>
"euro_area_georeferenced.tif"<br>
<br>
<br>
gdalwarp -s_srs "+proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0<br>
+y_0=0.0" -t_srs "+proj=latlong" "euro_area_georeferenced.tif"<br>
"euro_area_latlong.tif"<br>
<br>
My HDF5 file doesn't contain any infoprmation regarding GEOS<br>
projection. So, I had to do this<br>
<br>
<br>
Hope this helps<br>
<br>
If you get any news/information/method regarding this, please let<br>
me know<br>
<br>
Best regards<br>
Antonio<br>
<br>
<br>
__________ Information from ESET NOD32 Antivirus, version of<br>
virus<br>
signature database 4907 (20100302) __________<br>
<br>
The message was checked by ESET NOD32 Antivirus.<br>
<br>
<a href="http://www.eset.com" target="_blank">http://www.eset.com</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
__________ Information from ESET NOD32 Antivirus, version of virus<br>
signature database 4908 (20100302) __________<br>
<br>
<br>
The message was checked by ESET NOD32 Antivirus.<br>
<br>
<a href="http://www.eset.com" target="_blank">http://www.eset.com</a><br>
<br>
<br>
<br>
<br></div></div>
------------------------------------------------------------------------<div class="im"><br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div></blockquote><font color="#888888">
<br>
<br>
-- <br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/%7Ewarmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush | Geospatial Programmer for Rent<br>
<br>
</font></blockquote></div><br>