[R] map - mapproj : problem of states localisation

Ray Brownrigg Ray.Brownrigg at mcs.vuw.ac.nz
Wed Nov 14 00:38:56 CET 2007


Unfortunately, your code is not self-contained.

Please provide a completely self-contained set of commands which illustrate 
the problem (as described in the posting guide).

If you mean a problem that the boundaries don't quite line up, as illustrated 
by:
map("france")
map("world", col=2, add=T)

then try using the "worldHires" database from the mapdata package.

HTH,
Ray Brownrigg

On Wed, 14 Nov 2007, Amandine Chevalier wrote:
> 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_d0
> >1
>
> ),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]]



More information about the R-help mailing list