An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120519/103b08f4/attachment.pl>
converting csv to image file
4 messages · Belay Gebregiorgis, Simon Knapp, Michael Sumner
provided you get the call to read.table (or perhaps read.csv) right
and presuming that file contains only the image data, you should
be able to say:
r <- raster(as.matrix(read.csv("file.csv")))
extent(r) <- extent(xmin, xmax, ymin, ymax)
and not worry about the projection (if it is plain old decimal degrees)
Regards,
Simon
On Sat, May 19, 2012 at 7:49 PM, Belay Gebregiorgis <belayg12 at gmail.com> 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 at 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.
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120519/b2ec3cc7/attachment.pl>
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120519/00c4d97a/attachment.pl>