[R-sig-Geo] Plotting an image on a stereographic grid orientation

Dag Johan Steinskog dag.johan.steinskog at nersc.no
Fri Aug 12 09:45:14 CEST 2005


Dear all

First of all, I am sorry for all inconvinience I am causing you, but I 
have a problem on plotting some data on stereographic grids.

Say I have data given to be (for simplicity):
lon<-seq(-180,180,2.5)
lat<-seq(50,90,2.5)
data<-matrix(rnorm(145*17,sd=10),ncol=17,nrow=145)

And this I want to plot on a map like this (maps and mapproj package 
needed):
map("world",proj="stereographic",xlim=c(-180,180),ylim=c(50,90),interior=FALSE,lwd=1)
map.grid(lim=c(-180,180,50,90),nx=5,ny=3,pretty=TRUE,col=1,label=FALSE)

I would like to have a colour presentation of the data. Before I have 
used the image command to plot a datamatrix when there is no specific 
projection, and I would like to use the image command on this problem too.

I hope someone have time to help me.

Best Wishes
Dag Johan Steinskog




More information about the R-sig-Geo mailing list