Thanks Agus, I fixed this in raster version 1.3-0. Robert
On Mon, Jul 19, 2010 at 11:46 PM, Agustin Lobo <[hidden email]> wrote:
While
test <-rasterToPoints(Br025,fun=function(x){x>0},asSpatialPoints=T)
does not work, the following works:
test <-rasterToPoints(Br025,fun=function(x){x>0},asSpatialPoints=F)
test <- as.data.frame(test)
coordinates(test) <- c("x","y")
class(test)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
test at proj4string <- CRS(projection(Br025))
projection(test)
[1] "+proj=longlat +ellps=WGS84"
Object of class SpatialPointsDataFrame
Coordinates:
? ?min ? max
x -8.75 28.00
y 41.50 59.75
Is projected: FALSE
proj4string : [+proj=longlat +ellps=WGS84]
Number of points: 2576
Data attributes:
? Min. 1st Qu. ?Median ? ?Mean 3rd Qu. ? ?Max.
? 0.50 ? ?7.00 ? 33.75 ?178.00 ?152.10 3638.00
although the value of slot "Is projected:" indicates some problem.
Agus
--
View this message in context:
http://r-sig-geo.2731867.n2.nabble.com/Bug-in-rasterToPoints-tp5294219p5315475.html
Sent from the R-sig-geo mailing list archive at Nabble.com.