Spatstat中的样方

Question

我现在正在做探索性分析，目标是绘制一个样本地图，确定是否存在完全的空间随机性(视觉上和使用卡方检验)等。然而，我在绘制样本地图时遇到了问题。我之前问过this question，这导致我重新投影了我的数据。

下面是我的代码：

library(rgdal) #Brings Spatial Data in R
library(spatstat) # Spatial Statistics
library(maptools)
library(raster)

# Load nyc zip code boundary polygon shapefile 
s <- readOGR("/Users/my_name/Documents/fproject/zip","zip")
nyc <- as(s,"owin")

### OGR data source with driver: ESRI Shapefile 
#Source: "/Users/my_name/Documents/project/zip", layer: "zip"
#with 263 features

# Load nyc arrests point feature shapefile
s <- readOGR("/Users/my_name/Documents/project/nycarrests/","geo1")

### OGR data source with driver: ESRI Shapefile 
#Source: "/Users/my_name/Documents/project/nycarrests", layer: "geo1"
#with 103376 features
#It has 19 fields

#Converting the dataset into a point pattern
utm <- "+proj=utm +zone=32 +datum=WGS84"
x <- spTransform(s,utm)
arrests <- as(x,"ppp")

marks(arrests) <- NULL
Window(arrests) <- nyc 
plot(arrests, main="2020", cols="dark green",pch=20)

### Quadrat density 2020 
Q <- quadratcount(arrests, nx= 3, ny=3)
plot(arrests, pch=20, cols="grey70", main="2020")  # Plot points
plot(Q, main="2020", add=TRUE)  # Add quadrat grid

在这里，我的数据点受到纽约市邮政编码边界的约束，但它们不会显示在我的地图上。我能做些什么来解决这个问题呢？

Answer 1

您没有提供太多信息。例如，可以显示(打印) s (每次)或crs(s) (两者都加载了栅格包)

无论哪种方式，几乎可以肯定的是，这是因为您使用的坐标参考系不同。您可以将其中一个更改为"UTM zone 21"，但不更改其他分区。它们可能是不同的。此外，21区不适合纽约。