I am trying to convert a set of lat long to state plane coordinates using python and pyproj. The location is Tennessee and this gave me information about the central meridian, lat origin and the proj type which are 35.2500,34.333 and lambert conformal conic respectively. When i use National Geodetic survey's tool to do the conversion for a sample long lat(-86.911,36.316), i get this result:
NORTH(Y) EAST(X) AREA
220358.963 518192.337 TN
When i try doing this by code, the results are different:
Here is the piece of code:
from pyproj import Proj p = Proj(proj='lcc',zone=4100,lon_0='35.35w',lat_1='34.33n') x,y = p(-86.911,36.316) print x print y
Result x,y is: -4367838.06379 5640715.01601
I cannot understand where i am going wrong. Please advice.