我试图在纽约市的城市地图上绘制犯罪数据密度（再次！）。作为一个众所周知的问题，有很多例子（http://www.obscureanalytics.com/2012/12/07/visualizing-baltimore-with-r-and-ggplot2-crime-data/）。这些方法通过等倾线绘制犯罪密度，而我需要通过预定区域的相同密度方块来表示它（并且区域/边长可以从一次迭代变为另一次）。这实际上是在商业上可用的COTS包中完成的，比如PredPol（参见http://www.predpol.com）。通过正方形表示犯罪密度的原因是广场是要巡逻的热点区域。规模将影响所需警察的总数。
这就是我想要实现的目标：
我希望能够创建具有预定尺寸的相同方形多边形以覆盖到地图（这是一个栅格吗？道歉但我刚刚开始学习拼写GIS！） 我想使用上面的方块作为颜色的项目，如在等值区域图中（即，与该区域中的犯罪频率相关的不同着色），可能使用ggplot2或类似的。 这应该让我看到每平方公里的犯罪密度如何变化，改变广场的大小（即面积），提出不同的巡逻区域。 我不知道是否可以使用R来创建预定义大小的预定形状的正方形多边形（如下面的代码所示）。欢迎任何帮助或链接到示例。
我很乐意得到关于计算密度的替代方法的一些指示。我使用了stat_density2（ggplot2的一部分），但也许有更好/更快的方法？ （ 事后来看，我是否需要密度功能？我只需要计算一个单元格中的罪行并对其进行颜色绘制......）
这就是我要去的地方：
library(rgdal)
library(raster)
library(sp)
#NY boroughs shapefile downloaded from NY website
shp <- readOGR(dsn = "nybb_14a_av", layer = "nybb")
r <- raster(extent(shp))
res(r)=0.05
# using BoroCode as an experiment...
r <- rasterize(shp, field="BoroCode", r)
plot(r)
plot(shp,lwd=10,add=TRUE)
#don't know the result of the above: the laptop basically hangs processing
#plot(r) :)