2. 华东师范大学 地理科学学院, 上海 200241
2. School of Geographic Sciences, East China Normal University, Shanghai 200241, China
偷盗"三车"包括偷盗摩托车、电动车和自行车, 这样的犯罪案件普遍与地理环境密切相关且具有时空集聚特性[1-2].通过研究犯罪案件时空特征与预测方法, 例如热点制图[3-4], 时间序列预测[5-6]等能帮助警务人员更合理地安排警力部署[7]、平衡工作重心[8].
现有犯罪案件的分析与预测方法多样.有些研究关注时间序列分析预测, 例如使用探索性数据分析寻找现有案件的时间序列特征, 并据此预测犯罪率[5];有些研究关注案件的空间分布, 例如使用核密度[8]、聚类[9]等方法绘制犯罪热点预测图; 还有些研究结合了时空两方面同时考虑, 例如在时间热点基础上绘制空间热点预测图[8]、时空分析制预测图[3]等.但是少有研究同时关注时空尺度问题, 如果找到时空最优尺度, 将使绘制的热点图更准确可信, 具备预测功能.本文主要关注这种时空尺度问题.
本文在前人研究的基础上, 分别在时间和空间尺度对已有方法做了扩展: ①探查时间热点时, 提出尺度法, 并结合传统的严格法来整合数据; ②绘制空间热点预测图时, 改进DBSCAN (Density-Based Spatial Clustering of Applications with Noise)算法, 使其自适应选择最优聚类尺度.下面分别对主要技术的现有进展和本文贡献作出说明.
在确定案件发生具体时间时, 有平均划分及严格划分两种常用方式[10-11], 前者充分利用数据但探测热点有偏差, 而后者以舍弃部分数据为代价探测较准热点.针对这个问题, 本文提出一种尺度法整合数据, 先以严格法探查大尺度时间周期热点, 经降维压缩到一个周期, 继而在更小尺度探查热点, 可在缺失数据的情况下完成周期的提取.
在绘制财产类案件的空间热点预测图时, 有学者在聚类时以凸包绘制热点, 使热点划分更"软", 同时具备更高的预测准确率[4].但它难以确定最优聚类尺度, 热点预测图的绘制与准确率检验仍然分离, 有较多不便.本文针对财产类案件提出改进的DBSCAN聚类法(Modified DBSCAN, MDBSCAN), 即在DBSCAN聚类中嵌入凸包分析.同时, 以凸包案件密度中值(Median of Case Density in Convex Hull, MDCH)为评估函数, 来自适应选择最优尺度, 使热点预测图的绘制与准确率检验同时进行.并使用预测准确率指数(Prediction Accuracy Index, PAI)[8-12]验证MDBSCAN方法的有效性.
1 研究区与数据 1.1 研究区研究区位于某市中心城区的西部(由于数据涉密, 研究区隐去经纬度).其西南角为一机场, 故在西南角道路分布稀疏.境内道路分布密集, 由东至西分别被三个高架、一个高速贯穿, 并且4条地铁线路在区内形成井字形.
1.2 研究数据如图 1, 本文中涉及到的数据: ①研究区行政界线图和主要道路图; ② 2015年1——10月研究区偷盗"三车"案件数据, 由该区公安分局提供.案件原始数据共2 805个记录, 每个记录包括案件发生地址、时间等信息.
以上原始数据通过三步数据预处理进行后期的分析: ①结合百度API进行地理编码[13], 得到本文所使用的地理坐标经纬度; ②筛选出研究区范围内属性不为空的事件点作为研究对象; ③因地理格网体系形状规则尺度灵活, 信息集成和其他单元结合方便[14], 故本文在获得各月份最优热点图后, 使用格网进行最后的空间叠置分析.
2 方法本节将说明两种多尺度时空分析方法, 研究过程如图 2所示.具体涉及:以尺度法整合犯罪时间数据, 利用小波分析获得犯罪时间热点(见2.1节); 改进DBSCAN算法, 将其尺度变量多尺度化, 以绘制具有预测能力的犯罪空间热点(见2.2节).
案件记录中, 通常包含"三车"拥有者最后见到自己的车以及发现车丢失的时间, 这使得其中任意时刻都有可能成为真正的犯罪时间.选择合适的划分方法将使时间热点探测更准确.本文提出的尺度法如下: ①以主流的严格法划分时间数据, 并用小波提取案件的日尺度热点周期; ②由于小波图显示案件在日尺度上发生的频率相似, 故可以将其在小时尺度合并, 补充小时尺度的数据量; ③重复①步, 以严格法划分时间数据, 以小波提取小时尺度的热点周期.这解决了原本在小时尺度上案件量过少, 无法探测小时尺度热点的问题.
上文中提到的连续小波变换是一种处理非线性时间序列的常用方法, 它在处理不同频率的时间序列中效果更优[15], 且无论信号自身频率的高低, 小波变换都能够提取该信号的特征信息[16].所以, 本文借此提取偷盗"三车"的犯罪周期.
对于连续信号的连续小波变换定义如式(1)和(2)所示[17-18], 其中,
$ {C_X}(a,b) = \int_{ - \infty }^{ + \infty } {x(t){\psi _{a,b}}(t){\rm{d}}t} ,{\rm{ }} $ | (1) |
$ \begin{align} \psi_{a, b}(t)=|a|^{-\frac{1}{2}}\psi\Big(\frac{t-b}{a}\Big), a, b\in {\mathbf R}, a\neq0, \end{align} $ | (2) |
$ \begin{align} \psi (t)=\pi ^{-\frac{1}{4}}{\rm{e}} ^{{\rm{i}}, \omega_0t}{\rm{e}}^{-\frac{t^2}{2}}. \end{align} $ | (3) |
在此基础上需要计算小波功率谱以获得时间热点, 其被定义为小波系数的模, 而小波系数模的求法类似于傅里叶变换, 如式(4)所示.其中
$ \begin{align} W_X(a, b)=C_X(a, b)C_X^*(a, b)=|C_X(a, b)|^2. \end{align} $ | (4) |
DBSCAN[19]核心是使空间中一定距离(Eps)范围内的点对象数目达到最小指定数目(MinPt)的所有对象, 成为一个独立的分类; 以凸包为辅助绘制热点区域.将凸包分析嵌入DBSCAN聚类, 针对财产类犯罪提出一种改进的DBSCAN聚类方法(Modified DBSCAN, MDBSCAN), 在用户指定MinPt后自适应确定最优聚类尺度(Eps).
MDBSCAN法的主要步骤: ①根据指定MinPt使用不同尺度聚类并对每个尺度下的每类单独绘制凸包;
②分尺度遍历所有凸包案件密度, 并对各尺度的密度序列取中值(MDCH), 若只有一个凸包, 则舍弃该尺度, 进入下一次计算; ③最大MDCH值对应尺度为最优尺度.计算方法如式(5)所示, 其中,
$ \begin{align} {\rm {MDCH}}={\rm {Med}}\Big(\frac{C_{CH_i}(Eps)}{S_{CH_i}(Eps)}\Big)(i=1, 2, \cdots, n). \end{align} $ | (5) |
根据前述, 最优尺度聚类结果从热点探测角度看, 应当具备最好的热点预测准确率.目前, PAI被广泛采纳为热点预测精度指标[8, 20], 所以本文以PAI验证自适应选择最优尺度的效果. PAI由两部分组成, 具体公式如式(6)所示, 其中HHR为测准案件数与总案件数之比, AP为预测面积与研究区面积之比.在MDBSCAN算法中,
$ \begin{align} {\rm {PAI}}=\frac{{\rm {HHR}}}{\rm {AP}}=\frac{n_i/N}{a_i/A}(i=1, 2, \cdots, m). \end{align} $ | (6) |
偷盗"三车"类犯罪案件的时间数据从1月记录至10月, 将其整理为两个尺度:日尺度和小时尺度, 并且提取两种不同尺度下的周期.
首先我们将原始数据整理为日尺度数据, 小波功率越大, 相应时间间隔出现的频率越高.通过0.05显著度(粗黑线内部)的偷盗"三车"案件小波功率谱高值如图 3(a)所示. 8 d以下显示的周期可能是由于随机因素导致的, 相较之下30 d周期与65 d周期有更高的波谱能量.此处周期较明显可直接用作空间热点探查的周期.下面以求更小的时间尺度为例, 说明如何以尺度法整理时间数据, 获得更小尺度的周期信息.
本文分析小时尺度热点时, 由于每月(30 d为一个周期)案件发生频率具备周期性, 故而可以将每月的案件数按照相同的日期时刻进行次数统计, 得到每小时的案件数量, 从而尽可能地弥补单月造成的空值, 也使周期更具普遍意义.如图 3(b)的小波功率谱所示, 16
偷盗"三车"案件在日尺度的30 d热点与小时尺度的20 h热点可以很容易地理解为偷盗"三车"事件发生的周期. 20 h尺度, 值得警务实务工作重点关注.后文探查空间热点时, 考虑到以20 h周期计算, 其耗时较多, 所以以30 d为周期探查空间热点.
3.2 空间热点分析 3.2.1 MDBSCAN聚类寻找最优尺度下的空间热点根据2.2节提出的MDBSCAN, 本例中MinPt=10是有良好聚类变化的区间, 现通过迭代法计算各个尺度下的MDCH, 这里随机挑选4月为例, 以式(5)计算的结果如表 1所示, 图 4为各尺度下的DBSCAN聚类与凸包分析的结果.结合图表可知, 当Eps为0.21 km和0.95 km时凸包(类别)个数为1.前者所绘凸包过小, 则所围区域为极大热点, 虽有价值但不具备指导意义.而后者所绘凸包过大, 可能将凸包内部的聚类点模式作为均匀分布或随机分布处理, 所以其结果不具备参考价值.由此得出去除凸包个数为1的结论, 所以Eps为0.32 km时为最优尺度.
根据式(6)易知, PAI指数越高, 预测的准确率越高.同以4月预测5月为例.表 2显示4月预测5月的PAI值计算列表, 可知, 当去除凸包个数为1对应的
图 5(a-c)展示了4月预测5月、6月预测7月、9月预测10月的三个例子, 说明以MDCH和PAI挑选最优Eps的结果一致性; 再根据表 1与表 2的MDCH、PAI验证, 可知MDCH与PAI的值变化趋势相同.由此可见, 使用MDBSCAN选取的最优Eps的聚类结果具备优秀的热点制图预测能力.除了4月外, 我们还以MDBSCAN算法预测了其余8个月的案件, 最终结果如表 3所示.
使用MDCH作为评估函数, 自适应选择最优聚类尺度的有效性已被PAI验证.其有效的原因在于: ①根据式(5)与(6)可以发现, MDCH与PAI都以分子衡量热点预测的精准度, 而以分母(面积)作为惩罚项, 两者相除以衡量热点的准确率.所以MDCH与PAI实际衡量的是相似的指标; ② MDCH与PAI最大的区别在于, 前者在同一周期内衡量准确率, 后者以下一周期衡量上一周期准确率.研究已证明[20-21], 财产类犯罪存在临近重复犯罪的现象, 所以以3.1时间热点的结果作为预测的周期可以弱化时间因素, 使MDCH有效替代PAI.这种替代将在聚类的同时选择最优尺度绘制犯罪热点预测图, 简化过程的同时确保了准确率.
3.2.3 热点区域制图根据每月热点预测图可知下月易发生偷车案件的空间范围.本文使用2015年10个月的最优Eps聚类结果(最优尺度DBSCAN聚类图), 结合案件数量格网分布来绘制犯罪热点预测图(见图 6), 供警务参考使用.具体热点区域为研究区西北部的A地铁站周边, 以及东北部三块区域, B地铁站, C地铁站, 及D周边.上述区域为案件发生最为频繁而密集的地方, 也即, 它们从时间角度看, 是2015年10个月里频数最高的地区, 从空间角度看, 是MDBSCAN聚类密度最高的地区.地铁站因通勤人流量大, 来往地铁站与居住地的交通方式主要为自行车, 故有大量"三车"停放.而D周边有较多居住区, 除此之外, 共有四所高校在此周围, 商贩密集、大学生防护意识较差等多项因素综合形成案件热点区域.
本文针对以偷盗"三车"为例的犯罪案件进行热点分析与预测, 主要内容与方法包括: ①提出尺度法整合时间数据, 并以此法提取时间热点周期; ②将DBSCAN的尺度变量多尺度化, 设计针对性目标函数改进DBSCAN算法, 单独或叠置使用后探查最优空间热点.
结果表明: ①尺度法与严格法的结合, 可以在数据量不足的情况下探测到犯罪周期, 弥补了严格法的缺陷. ②依据探测到的时间热点(周期), 以MDBSCAN算法绘制的热点区域, 能够缓解划分过硬问题, 且以MDCH为评估函数时可以将空间热点探测与预测评估同时进行, 提高效率, 保持准确率. ③尺度法整合数据后小波分析可知, 研究区内的偷盗"三车"案件具有较明显的30 d周期和极明显的20 h周期.这一结论可能与警务工作人员的轮班岗位时间等有关, 应当引起警务工作者注意. ④ MDBSCAN获取各周期最优空间热点图并叠置分析后, 发现个别地铁站与居住区属于高发区, 应在后期的实际警务行动中重点预防.本文将多尺度分析分别应用于时间与空间热点探查, 未来还需进一步研究将其用于时空结合的分析方法中.
[1] | BRANTINGHAM P L, BRANTINGHAM P J. Patterns in Crime[M]. New York: Macmillan, 1984: 1-403. |
[2] | GRUBESIC T H, MACK E A. Spatio-temporal interaction of urban crime[J]. Journal of Quantitative Criminology, 2008, 24(3): 285-306. DOI:10.1007/s10940-008-9047-5 |
[3] | BRUNSDON C, CORCORAN J, HIGGS G. Visualising space and time in crime patterns:A comparison of methods[J]. Computers, Environment and Urban Systems, 2007, 31(1): 52-75. DOI:10.1016/j.compenvurbsys.2005.07.009 |
[4] | GRUBESIC T H. On The Application of fuzzy clustering for crime hot spot detection[J]. Journal of Quantitative Criminology, 2006, 22(1): 77-105. |
[5] | ISMAIL S, RAMLI N. Short-term crime forecasting in Kedah[J]. Procedia-Social and Behavioral Sciences, 2013, 91: 654-660. DOI:10.1016/j.sbspro.2013.08.466 |
[6] | BUNN D W, VASSILOPOULOS A I. Using group seasonal indices in multi-item short-term forecasting[J]. International Journal of Forecasting, 1993, 9(4): 517-526. DOI:10.1016/0169-2070(93)90078-2 |
[7] | GORR W, HARRIES R. Introduction to crime forecasting[J]. International Journal of Forecasting, 2003, 19(4): 551-555. DOI:10.1016/S0169-2070(03)00089-X |
[8] | TOMPSON L, TOWNSLEY M. (Looking) Back to the future:Using space——time patterns to better predict the location of street crime[J]. International Journal of Police Science & Management, 2010, 12(1): 23-40. |
[9] | CRAGLIA M, HAINING R, WILES P. A comparative evaluation of approaches to urban crime pattern analysis[J]. Urban Studies, 2000, 37(4): 711-729. DOI:10.1080/00420980050003982 |
[10] | RATCLIFFE J H. Aoristic analysis:the spatial interpretation of unspecific temporal events[J]. International Journal of Geographical Information Science, 2000, 14(7): 669-679. DOI:10.1080/136588100424963 |
[11] | RATCLIFFE J H, MCCULLAGH M J. Aoristic crime analysis[J]. International Journal of Geographical Information Science, 1998, 12(7): 751-764. DOI:10.1080/136588198241644 |
[12] | PATTEN I T V, MCKELDIN-CONER J, COX D. A microspatial analysis of robbery:Prospective hot spotting in a small city[J]. Crime Mapping:A Journal of Research and Practice, 2009, 1(1): 7-32. |
[13] | 杜传明. 百度地图API在小型地理信息系统中的应用[J]. 测绘与空间地理信息, 2011, 34(2): 152-153. |
[14] | 兀伟, 邓国庆, 武晓莉, 等. 地理格网及在统计分析中应用的探讨[J]. 测绘标准化, 2015(3): 22-26. |
[15] | JEVREJEVA S, MOORE J C, GRINSTED A. Influence of the arctic oscillation and El Niñ-Southern Oscillation (ENSO) on ice conditions in the Baltic Sea:The wavelet approach[J]. Journal of Geophysical Research:Atmospheres, 2003, 108(D21): 1-11. |
[16] | MIAO J, LIU G, CAO B, et al. Identification of strong karst groundwater runoff belt by cross wavelet transform[J]. Water Resources Management, 2014, 28(10): 2903-2916. DOI:10.1007/s11269-014-0645-8 |
[17] | WANG W, HU S, LI Y. Wavelet transform method for synthetic generation of daily streamflow[J]. Water Resources Management, 2010, 25(1): 41-57. |
[18] | LABAT D. Cross wavelet analyses of annual continental freshwater discharge and selected climate indices[J]. Journal of Hydrology, 2010, 385(1): 269-278. |
[19] | ESTER M, KRIEGEL H, SANDER J, et al. A density-based algorithm for discovering clusters in large spatial databases with noise[C]. Knowledge Discovery and Data Mining, 1996: 226-231. |
[20] | 徐冲, 柳林, 周素红. DP半岛街头抢劫案件的临近重复发生模式[J]. 地理研究, 2015, 34(4): 384-394. |
[21] | SHORT M B, D'ORSOGNA M R, BRANTINGHAM P J, et al. Measuring and modeling repeat and near-repeat burglary effects[J]. Journal of Quantitative Criminology, 2009, 25(3): 325-339. |