[Liblas-devel] Precision of reprojection transforms

Frank Warmerdam warmerdam at pobox.com
Fri Apr 29 10:03:53 EDT 2011


On 11-04-29 08:35 AM, Edwards, David wrote:
> I am trying to use libLAS to read an LAS file and perform a coordinate
> transformation.  Points in the file are stored in UTM coordinates.  I would
> like to transform them into geocentric coordinates.  As best I can tell, the
> version of libLAS I am using (1.6.1) does not support conversion to
> geocentric coordinates, so my backup plan is to convert to WGS84 geodetic
> coordinates (lat., long. , height above ellipsoid), and then perform the
> geodetic-to-geocentric conversion myself.
>
> The problem is that when I set up a ReprojectionTransform to convert the
> points from UTM to geodetic coordinates, the resulting latitude and
> longitude of each point has nonzero digits only to the hundredth of a
> degree, even though the UTM coordinates of my input points have data to the
> millimeter.  For example, here are the results of two such conversions:
>
> UTM (Easting, Northing) --->  Geodetic (Longitude, Latitude) 831805.86
> 4296378.26 --->  -77.18 38.75 831806.31 4296377.64 --->  -77.18 38.75
>
> Even though the points have different UTM easting and northing coordinates,
> they are transformed to the same latitude and longitude.  Is this limited
> precision inherent in the way libLAS does its coordinate conversions, or am
> I setting up the transform object(s) incorrectly?
>
> I apologize if this is a duplicate post.  I tried to send this yesterday
> before joining the mailing list.

Dave,

Reprojection is generally accurate to a small fraction of a millimeter.
I believe the problem you are encountering is how you are formatting the
coordinates for reporting - not the nature of the reprojection itself.

Best regards,
-- 
---------------------------------------+--------------------------------------
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    | Geospatial Programmer for Rent



More information about the Liblas-devel mailing list