使用R中的坐標創建密度圖

Question

先感謝您

我正在嘗試創建一個農場的密度圖，我有其坐標。 我有農場長/緯度坐標，並已使用以下代碼將其轉換為空間數據並繪制在新西蘭地圖上；

# Read in New Zealand shape file 
NZ <-readShapeLines("ESRI_Shapefile_2016_Digital_Boundaries_High_Def_Clipped/TA2016_HD_Clipped.shp")

# Change projection 
proj4string(NZ)<-CRS("+proj=tmerc +lat_0=0.0 +lon_0=173.0 +k=0.9996 +x_0=1600000.0 +y_0=10000000.0 +datum=WGS84 +units=m")

# Plot map and point data 
plot(NZ)
points(Farm_points) 

有什么簡單的方法可以將點數據轉換為密度數據並在地圖上繪制？

我嘗試使用GISTools軟件包和以下代碼計算每個多邊形中的點數；

counts <- poly.counts(Farms_points, NZ)

但是我收到以下錯誤消息；

Error in identical(spgeom1@proj4string, spgeom2@proj4string) : 
  trying to get slot "proj4string" from an object of a basic class ("matrix") with no slots 

任何建議歡迎...

謝謝

Answer 1

沒有您的數據，我們無法測試為什么它不起作用，但是，錯誤消息很明顯：如果您希望poly.counts起作用，則需要將矩陣Farms_points轉換為SpatialPoints ：

Farms_points_sp <- SpatialPoints(Farm_points, proj4string = proj4string(NZ))
counts <- poly.counts(Farms_points_sp, NZ)