Displaying 1 result from an estimated 1 matches for "longicrs".
Did you mean:
logics
2008 Jul 09
0
"Rotated Lat-Lon" projection in mapproj
...ot other outputs in a lambert projection would be:
##Read lat/lon data
lat<-as.numeric(system(paste(pycmd,'control',"200108150000","LATITCRS",1,0,sep="
"),intern=TRUE))
lon<-as.numeric(system(paste(pycmd,'control',"200108150000","LONGICRS",1,0,sep="
"),intern=TRUE))
## Projected coordinates
coords_lamb<-mapproject(lon,lat,projection="lambert",parameters=c(60.,30.))
xs<-t(matrix(coords_lamb$x,ncol=NCC,nrow=NRC))[,1]
ys<-t(matrix(coords_lamb$y,ncol=NCC,nrow=NRC))[1,]}
## Plot map
filled.contour(x...