[Proj] Cassini coordinate reprojection problem
Chris Faulkner
christopher.faulkner at ntlworld.com
Tue Mar 15 18:12:47 EST 2005
Hello
I have been sent some coordinates reprojected by another package that do
not match what I am getting out of proj4.
The projection is from Cassini to Transverse Mercator, details below.
I think the other package is genamap. The TM coordinate reported by
genamap from the reprojection is 535543,175269. I have downloaded
geotrans and it gives exactly the same result.
Proj4 gives this - 535596.86 175379.03. I am on RedHat Linux 9.0,
running 4.4.9
Most likely, I am driving proj4 incorrectly. Here is what I do
# Convert to lat /long (it doesn't matter if I remove the es and a
parameters)
chris at dinky /cygdrive/c/mapmarketing/software/tabrunner
$ invproj +ellps=airy +proj=cass +a=6377563.3959999 +es=0.0066705398
+lat_0=51.2695503 +lon_0=0.665434
-49630.956 21390.870
0d2'55.474"W 51d27'34.779"N
# Convert the lat / long to Transverse Mercator
chris at dinky /cygdrive/c/mapmarketing/software/tabrunner
$ proj +ellps=airy +proj=tmerc +lat_0=49 +lon_0=-2 +x_0=400000
+y_0=-100000
0d2'55.474"W 51d27'34.779"N
535596.86 175379.03
proj +proj=list gives an error message. Projection initialization
failure
Cause: unknown projection id
Any help would be most welcome !
Thanks
Chris
More information about the Proj
mailing list