Hi All,
I am trying to create a grid but something getting wrong 
here is the code 

rm(list=ls())
library(maptools)
library(rgdal)

options("digits"=22)
lat=seq(from= -123.4, to= -121.4, length=20)
long=seq(from=48.95, to=49.5, length=10)
#generates the grid
grd = expand.grid(lat,long)
sp_grd <- SpatialPoints(cbind(grd[,2],grd[,1]), proj4string=CRS("+proj=longlat +datum=WGS84"))

#do projection

spTransform(sp_grd,CRS("+proj=merc +datum=WGS84"))



any hint is appreciated

thx

	[[alternative HTML version deleted]]


