[gdal-dev] TransformPoint returning Infinity
Frank Warmerdam
warmerdam at pobox.com
Tue Oct 20 10:35:59 EDT 2009
Henneke, Amanda M wrote:
> Hello all-
>
> I've been running into a problem where I'm trying to transform points into UTM NAD83 and occasionally (does not seem to be repeatable, more random) the TransformPoint method returns Infinity for the x and y coordinate. Just for your info, the coordinate system I'm converting from is UTM NAD27. Anyone else seen this? I've included my code below...am I doing something wrong or is there some sort of an issue with the TransformPont method?
>
> // Get the input layer's spatial reference information
> OSGeo.OSR.SpatialReference sr = lyrOriginal.GetSpatialRef();
> string sSpRefWKT = "";
>
> // If this is a valid spatial reference, then export the data to a string and return it
> if (Validator.IsValid_SpatialRef(sr))
> {
> sr.ExportToPrettyWkt(out sSpRefWKT, 1);
> }
> sr.Dispose();
>
> OSGeo.OSR.SpatialReference srOld = new OSGeo.OSR.SpatialReference(sSpRefParams);
> OSGeo.OSR.SpatialReference srNew = new OSGeo.OSR.SpatialReference("");
>
> srNew.SetUTM(11, 1);
> srNew.SetWellKnownGeogCS("NAD83");
>
> ctTransform = new OSGeo.OSR.CoordinateTransformation(srOld, srNew);
>
> double[] dLnSegPt1 = { dX1, dY1 };
> ctTransform.TransformPoint(dLnSegPt1);
>
Amanda,
TransformPoint() takes a three entry array, try it with:
double [] dLnSegPt1 = { dX1, dY1, 0.0 };
ctTransform.TransformPoint(dLnSegPt1);
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 gdal-dev
mailing list