[R] mapproject Lambert in metres

From: Nevil Amos <nevil.amos_at_gmail.com>
Date: Mon, 21 Jun 2010 14:07:51 +1000


I have the following script to convert Longitude and Latitude into lambert projection ( adapted from Geneland manual). How do I get the proejction with units in metres rather than whatever they are currently? alos is it possible to add a fualse XY so that all values are positive ( all my lattidudes are south ie negative)

Many thanks

Nevil Amos

Monash University

LONGLAT is a headerless two column file of longitude and latitude such as

141 -32
142 -35
136 -20

## load packages

require(mapproj)
require(maps)
## read raw data

coord.lonlat <- LONGLAT
## convert (Lon,Lat) coordinates into Lambert
mapproj.res <- mapproject(x=coord.lonlat[,1], y=coord.lonlat[,2],
projection="lambert",
param=c(min(coord.lonlat[,2]),max(coord.lonlat[,2])))
## save planar coordinates as a two-column matrix
coord.lamb <- cbind(mapproj.res$x,mapproj.res$y)



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon 21 Jun 2010 - 04:11:09 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 21 Jun 2010 - 04:50:33 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive