有正方形像素的R-spatstat im？(出口至弧栅)
有正方形像素的R-spatstat im？(出口至弧栅)
提问于 2015-04-03 15:46:21
我在将密度im从spatstat导出到ArcGIS可读的文件格式时遇到了问题。这是我的密码
library(raster)
library(spatstat)
library(maptools)
library(sp)
# make a spatstat ppp with California boundary as window
ca <- readShapePoly("ca.shp")
o3 <- readShapePoints("o3.shp")
o3 <- as(o3, "SpatialPoints")
o3p <- as.ppp(o3)
o3p$window <- as.owin(ca)
# calculate density
d.o3p <- density.ppp(o3p)一切都很好。但是当我试图导出到ascii光栅文件时
writeRaster(raster(d.o3p), filename="grid.asc", format="ascii", NAflag=-9999)我知道这个问题
Error in .startAsciiWriting(x, filename, ...) :
x has unequal horizontal and vertical resolutions. Such data cannot be
stored in arc-ascii format我用的是加州的数据，所以高宽比不是1。所以……我如何使密度我有方形像素？
回答 1
Stack Overflow用户
回答已采纳
发布于 2015-04-03 20:58:24
ArcGIS必须使用ascii格式吗？使用任意窗口可能很难获得完美的方形像素。像素分辨率由参数eps或参数dimyx (而不是两者都控制)控制，它们被发送到as.mask，您可以在该帮助文件中读取更多信息。基本上，eps是像素大小，所以最简单的方法是将其设置为给定的值。然而，如果窗口的大小不完全符合窗口在一个方向上，像素可能是轻微的非正方形。例如，这一方法非常有效：
X <- runifpoint(100, win = owin(c(0,9),c(0,10)))
d <- density.ppp(X, eps = 1)但这提供了非方形像素：
d2 <- density.ppp(X, eps = 2)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/29435390
复制相关文章
相似问题
腾讯云开发者
