CGI script to convert lat/long to cartesian coordinates for the Lambert Conformal Conic projection

get_cart.py

#!/usr/bin/env python

# This file gets raw cartesian coordinates for the Lambert Conformal Conic projection given a latitude/longitude pair.
# It is used to cross reference a user's location with an ESRI shapefile database.

import os
from urlparse import parse_qs
print "Content-type: text/html\n\n"

from pyproj import Proj

p = Proj(proj='lcc',ellps='WGS84')

qs = parse_qs(os.environ['QUERY_STRING'])
cdt = p(qs['lon'][0], qs['lat'][0])
print '{0},{1}'.format(cdt[0], cdt[1])