座標を使って密度プロットを作成する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_sp <- SpatialPoints(Farm_points, proj4string = proj4string(NZ)) 
counts <- poly.counts(Farms_points_sp, NZ)