#llcrnrlon,llcrnrlat are lon and lat (in degrees) of lower
# left hand corner of projection region.
-where `proj='lcc` selects the [Lambert conformal conic][lcc]
+where `proj='lcc'` selects the [Lambert conformal conic][lcc]
projection for the x/y points, and `ellps='clrk66'` selects the
reference ellipsoid for the lat/lng coordinates. The other
coordinates are LCC parameters that select the [AWIPS 221][awips221]
>>> lines = plot(xs, ys, 'r.-')
>>> show()
+I've written up a simple script using this approach:
+[[geoscript/maproute.py]]. I've also written up a simple script to
+draw a map with labeled points: [[geoscript/maplabel.py]].
+
Note that you can easily get lat/lng pairs using [geopy][] (ebuild in
my [[Gentoo overlay]]):
[awips226]: http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html#GRID218
[lcc-param]: http://www.remotesensing.org/geotiff/proj_list/lambert_conic_conformal_2sp.html
[UTM]: http://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system
-[zone]: http://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system#UTM_zone
+[zones]: http://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system#UTM_zone
[geopy]: http://code.google.com/p/geopy/
[GeographicLib]: http://geographiclib.sourceforge.net/html/