Skip to content
Prev 1978 / 29559 Next

Good projection for N/S America?

On Mon, 9 Apr 2007 White.Denis at epamail.epa.gov wrote:

            
The attached script shows how to do the interrupted sinusoidal projection 
using spTransform in rgdal, for the whemi.lin data posted with the 
free-standing functions by Denis White a couple of days ago. Once the 
lines are converted into SpatialLines objects, the rest is robust and 
simple, as is the use of gridlines() in sp. The one catch is calculating 
the offset, here in an x_0= offset along the Equator in metres between the 
two central longitude values. The output is attached as a PNG image. The 
point about the sp objects is that they contain enough metadata (here a 
PROJ.4 projection description) to let them be moved to other R packages or 
external software. 

The half-dozen basic projections are easy to specify in PROJ.4, for 
example from the geotiff list:

http://www.remotesensing.org/geotiff/proj_list/

which is what I used here. The other projections mentioned are:

Lambert Cylindrical Equal Area "+proj=cea +lon_0=-80"

Lambert Azimuthal Equal Area "+proj=laea +lat_0=0 +lon_0=-80"

while the Northern hemisphere sinusoidal is:

Sinusoidal "+proj=sinu +lon_0=-100"

So I'd argue that PROJ.4 projection descriptions are not difficult to use, 
and with sp objects, do stay stuck to the data (has anyone else ever 
forgotten what projection was used when revisiting data, not just me?).

Using the maptools map2SpatialLines() interface function, or the Rgshhs() 
interface to GSHHS shorelines, even getting the lines is quite easy, 
qualified by clipping and bounding box issues in extremities for 
projection from geographical coordinates.

Of course, it would help to have MacOS X and selected Linux binaries of 
rgdal, we're very lucky that Uwe Ligges is so helpful with the Windows 
binaries.

Roger