Сейчас я занимаюсь исследовательским анализом, и цель состоит в том, чтобы построить карту квадрата, определить, существует ли полная пространственная случайность (визуально и с помощью теста хи-квадрат) и т. Д. Однако у меня возникают проблемы с построением карты квадрата. Ранее я задавал этот вопрос, который заставил меня перепроецировать мои данные.
Вот мой код:
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
Здесь мои точки данных ограничены границами почтового индекса Нью-Йорка, но они не отображаются на моей карте. Что я могу сделать, чтобы это исправить?
arrestsраспечатайте его, чтобы увидеть, сколько точек было включено, и прочтите описание пространственной области, чтобы убедиться, что это то, что вы ожидали. Вы, вероятно, обнаружите, что забыли пространственно преобразовать некоторые данные. - person Adrian Baddeley schedule 19.12.2020