I am trying to convert a pair of coordinates specified in US survey feet on a custom transverse mercator projection to lat/long using PyProj. I must be doing something wrong because I cannot get the correct decimal degree values:
import pyproj
# Coordinates in ft and false easting/northing in ft
x1,y1 = [1383566.534, 1627687.047]
sourceProjection = "+proj=tmerc +lat_0=42 +lon_0=-115.5 +k=0.9996 +x_0=1640416 +y_0=328083 +datum=NAD83 +units=us-ft +no_defs"
x2,y2 = pyproj.transform(pyproj.Proj(sourceProjection), pyproj.Proj(init='epsg:4326'), x1, y1)
print "Original: {0}, {1}".format(x1, y1)
print "Transformed: {0}, {1}".format(x2, y2)
The output I get is:
Original: 1383566.534, 1627687.047
Transformed: -119.384930994, 53.6300677515
but the correct result is -116.503170 45.562201
What am I doing wrong please?