posts:mutt-ldap: Mention Wade and Niels and add tags

[blog.git] / posts / pyproj.mdwn

diff --git a/posts/pyproj.mdwn b/posts/pyproj.mdwn
index 9a10e4072aec9a2ef90836918eda25bc9c145f59..ce328484b84b1e4a5b741cdf68152367ca25cb54 100644 (file)
--- a/posts/pyproj.mdwn
+++ b/posts/pyproj.mdwn
@@ -58,13 +58,14 @@ To plot the above New York to Vancouver route on a flat map, we need a
     #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]
 and [AWIPS 226][awips226] projections respectively (`lat_1`
 corresponds to `Latin1`, `lat_2` corresponds to `Latin2`, and `lon_0`
-corresponds to `Lov`).
+corresponds to `Lov`; see [this description][lcc-param] of the
+two-standard-parallel LCC and its PROJ.4 parameters).
 
 Convert our lat/lng pairs into grid points:
 
@@ -126,7 +127,12 @@ lat/lng pairs stored in a text file:
     >>> lines = plot(xs, ys, 'r.-')
     >>> show()
 
-Note that you can easily get lat/lng pairs using [geopy][]:
+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]]):
 
     >>> import geopy
     >>> g = geopy.geocoders.Google()
@@ -149,8 +155,9 @@ conversions, [GeographicLib][] looks promising.
 [LCC]: http://en.wikipedia.org/wiki/Lambert_conformal_conic_projection
 [awips221]: http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html#GRID221
 [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/