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

mett mettbest at gmail.com
Wed Jul 6 16:54:27 EDT 2011


Jean-Claude Repetto wrote:
> 
> Probably because your projection is not cylindrical. What projection are 
> you using ?
> Jean-Claude
> 
> The projection is
> `PROJCS["User-Defined WGS84/Lambert Conformal Conic (IS.)",GEOGCS["WGS
> 84",DATUM["WGS_1984",SPHEROID["WGS
> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",64.25],PARAMETER["standard_parallel_2",65.75],PARAMETER["latitude_of_origin",65],PARAMETER["central_meridian",-19],PARAMETER["false_easting",500000],PARAMETER["false_northing",500000],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]'
> 
> I did the follow step:
> - 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();
> 
> So, what should I do to have correct values from the conversion? I cannot
> find a good solution.
> 
> Thank you very much for any help,
> Matt
> 
> 
> 
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
> 


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


More information about the gdal-dev mailing list