This already in your question:
"+init=epsg:4326"
That is equivalent to "+proj=longlat +ellps=WGS84" just looked up via the
EPSG code. A good site to explore this topic is spatialreference.org
On Saturday, May 19, 2012, Belay Gebregiorgis wrote:
> Hello everyone,
>
> I want to get a 1km by lkm grid raster image using my csv data. If I call
> latitude=a, longitude=b and preciptation=c.
>
> a<-(1,2,3,4,5)
> b<-(6,7,8,9,10)
> c<-(10,20, 30,40, 50)
>
> Then I found an example in r help which goes like
>
> pts = read.table("file.csv",......)
> library(sp)
> library(rgdal)
>
> proj4string(pts)=CRS("+init=epsg:4326") # set it to lat-long
> pts = spTransform(pts,CRS("insert your proj4 string here"))
> gridded(pts) = TRUE
> r = raster(pts)
> projection(r) = CRS("insert your proj4 string here")
>
> Because I am new to R, I have no idea what to put into the proj4 string and
> all that. Can anyone help me on this? If there is a different way of doing
> this, that is also fine.
>
> Kind Regards,
>
> Belay
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org <javascript:;> 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.
>
--
Michael Sumner
Institute for Marine and Antarctic Studies, University of Tasmania
Hobart, Australia
e-mail: mdsum...@gmail.com
[[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.
