[gdal-dev] Re: [gdal-java-1.8.0]TransformPoint() - different values for same inputs, why?

mett mettbest at gmail.com
Wed Jul 6 17:35:06 EDT 2011


Jean-Claude Repetto wrote:
> 
> On 07/06/11 22:54, mett wrote:
> 
> I think your results are correct, because PROJ.4 gives exactly the same 
> ones :
> 
> proj -I -f "%.12f" +proj=lcc +lat_1=64.25 +lat_2=65.75 +lat_0=65 
> +lon_0=-19 +x_0=500000 +y_0=500000 +datum=WGS84 +units=m +no_defs
> 258321.25943477 655225.28062628
> -24.391630217716	66.299640570965
> 258321.25943477 574605.28062628
> -24.242820669640	65.579229590843
> 
> Jean-Claude
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
> 

Thank you again.
Just to close definitely this discussion, do you think the following steps
are correct in this case? 
- Got the Upper Left corner of the image with
GetGeoTransform(adfGeoTransform) 
- Found out the others corners in this way: 
double miny = adfGeoTransform[3] + width*adfGeoTransform[4] +
height*adfGeoTransform[5]; 
double maxx = adfGeoTransform[0] + width*adfGeoTransform[1] +
height*adfGeoTransform[2]; 
- Convert all these 4 points related to the corners with TransformPoint(); 

Especially the point where I found the other corners, by moving in the
raster space.

Thank you,
have a nice day
MATT


--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/gdal-java-1-8-0-TransformPoint-different-values-for-same-inputs-why-tp6555787p6556106.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.


More information about the gdal-dev mailing list