Dear Simon,
I understand better now and it is running.
//Belay
On Sat, May 19, 2012 at 1:22 PM, Simon Knapp sleepingw...@gmail.com wrote:
I'm not sure what line you mean. The thing you pass to raster should
be a matrix of the dimensions of your image data.
You should have have xmin, xmax,
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)
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
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
Thanks Simon. I have other vectors in the CSV file so after naming my
latitude and longitude x and y, I tried this. The first two lines run fine
but I get error in the second line.
r - raster(as.matrix(F1))
extent(r) - extent(xmin, xmax, ymin, ymax)
F1-data.frame(x, y,z)
Error in function
5 matches
Mail list logo