[R-sig-Geo] World map in magnetic coordinates

Clint Bowman clint at ecy.wa.gov
Sat Oct 9 00:24:08 CEST 2010


Why not plot data and map in original coordinates and then overlay 
with grid lines showing magnetic coordinates?

-- 
Clint Bowman			INTERNET:	clint at ecy.wa.gov
Air Quality Modeler		INTERNET:	clint at 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.okike at 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 at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>


More information about the R-sig-Geo mailing list