[gdal-dev] Require Help Transforming into Equirectangular

[ssobieraj]

Hello,

I'm trying to use FWTools to transform from an arbitrary Datum/Projection
pair in an arbitrary format into a Geographic (lat/lon) projection with the
WGS84 ellipsoid in GeoTiff format.

To do this, I've tried the following:

gdalwarp -s_srs Input.prj -t_srs Output.prj Input.xyz Output.tif

Where Input.prj is formatted in Well-Known-Text format and contains
information specifically for the Input.xyz file.
Output.prj is also in WKT format but contains the following:

GEOGCS["Geographic Coordinate
System",DATUM["WGS84",SPHEROID["WGS84",6378137,298.257223560493]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]

***
I've been successful in using this set-up for transforming from UTM
projection in WGS84 datum, but I am having trouble transforming from
Sinusoidal projection in RT90 datum.  The output file is no longer skewed in
the Sinusoidal projection, but I cannot seem to get the datum information to
coincide, resulting in a 10+ km discrepancy.

A copy of the WKT-formatted PRJ file is seen below:

PROJCS["Sinusoidal",GEOGCS["Geographic Coordinate
System",DATUM["RT90",SPHEROID["Bessel
1841",6377397.155,299.1528131060786]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Sinusoidal"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Any ideas as to why I would be experiencing this offset?



Regards,
Scott


-- 
View this message in context: http://osgeo-org.1803224.n2.nabble.com/Require-Help-Transforming-into-Equirectangular-tp5789691p5789691.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.