Hi R-user,
I am new with R and I have a problem with the map and mapproj fonctions :
with the following code :
>test_proj=mapproject(LONG_d01_vec_t1, LAT_d01_vec_t1, projection="lambert",
parameters=c(30,60))
>longitude_vec=test_proj$x
>latitude_vec=test_proj$y
>longitude_mat=matrix(longitude_vec,ncol=e_sn_d01)
>latitude_mat=matrix(latitude_vec,ncol=e_sn_d01)
>coordonnees=data.frame(longitude_vec,latitude_vec)
>image.plot(as.image(MOD_VENT_d01_t1,coordonnees,nrow=e_we_d01,ncol=e_sn_d01
),zlim=c(0,20),xlab="LONG",ylab="LAT",col=positiv_color(20),axes=FALSE)
>map("world",proj="lambert",pa=c(30,60),
lwd=1,col="slategrey",add=TRUE)
>map.grid(c(-20,30,35,55),labels=TRUE,nx=5,ny=5,pretty=TRUE)
>map("france",proj="lambert",pa=c(30,60),
lwd=1,col="slategrey",add=TRUE)
>map.grid(c(-20,30,35,55),labels=TRUE,nx=5,ny=5,pretty=TRUE)
>map("italy",proj="lambert",pa=c(30,60),
lwd=1,col="slategrey",add=TRUE)
>map.grid(c(-20,30,35,55),labels=TRUE,nx=5,ny=5,pretty=TRUE)
France and Italy are not co-localised with France and Italy from the world
mapÂ…
Now, I suppose that any of the three maps are well situatedÂ…
Do you have any idea for my problem?
Thanks very much for any help,
Regards,
Amandine
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
