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. _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
