[R-sig-Geo] kde output from ks reversed

[Robert Hijmans]

> Here is a self contained sample with a different set of points, same issues
though. 
>
> curset <- read.delim("sample.csv", header = TRUE, sep = ",") 

Thanks for trying, but this is not self-contained, as it refers to a file
that only you have. 


> When I use your routines they fail, saying that the cells unequal in
> horizontal and vertical dimensions. 

That is a limitation of the ascii format, you can write to a different
format. 


Here is a self contained example that *looks* OK to me. 

library(ks) 
library(raster) 
set.seed=1
curset <- cbind(rnorm(100), rnorm(100))
Hpi1 <- Hpi(x = curset, pre="scale") 
ud_hpi <- kde(x = curset, H=Hpi1) 
r <- raster(ud_hpi) 
plot(r) 
points(curset)
writeRaster(r, 'test_r.tif') 

-- 
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/kde-output-from-ks-reversed-tp6037090p6041894.html
Sent from the R-sig-geo mailing list archive at Nabble.com.