R-spatstat im与方形像素?(用于导出到弧形光栅)

Question

将spatstat中的密度im导出为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。那么…如何使密度为平方像素？

Answer 1

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)