[Proj] Re: To handle geographic info from HDF5 formated Satellite Data

[Frank Warmerdam]

Aydin wrote:
> Frank,
> 
> Thank you for response.
> if I use your suggestion, for example pixel (1000,1000);
> 
> longitude = -25 + 0.560648 * 1000 = 5581.48
> latitude = 75 - 0.560648 * 1000 = -5531.48
> 
> So, it is not possible.
> 
> I do not have problem converting projected space to geographic space using 
> proj lib. My problem is how can calculate the distance from the corner using 
> below parameters than I can convert those to lat/lon.

Aydin,

I stand corrected - the pixel size turns out to be in km not decimal degrees.

So you will have to offset from the top left corner using:

x (easting) = -2782987.27 + 560.648 * pixel
y (northing) = 8348961.81 - 560.648 * line

Then reproject that point back to long/lat using

proj -I +proj=eqc +lat_ts=0 +lon_0=0.000000 +ellps=WGS84

Good luck,
-- 
---------------------------------------+--------------------------------------
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    | President OSGeo, http://osgeo.org