Why not plot data and map in original coordinates and then overlay
with grid lines showing magnetic coordinates?
--
Clint Bowman INTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600 FAX: (360) 407-7534
Olympia, WA 98504-7600
On Fri, 8 Oct 2010, Barry Rowlingson wrote:
On Fri, Oct 8, 2010 at 7:21 PM, ogbos okike <ogbos.ok...@gmail.com> wrote:
library(maps)
data(worldMapEnv)
map("world")
The three steps above in R gives you the world map I mean.
if you do:
w = map("world")
then w will be a simple list with $x and $y components, as well as
$range and $names. If you write the $x and $y to a file, convert, then
read back into $x and $y (and then convert the $range, which is the
bounding box) you should have something you can plot in transformed
coordinates.
Note that polygons are separated by an NA in $x and $y, so you will
have to make sure you preserve those too.
Barry
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo