密度 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 ではありません。では、密度を正方形のピクセルにする方法を教えてください。