土地利用空间格局是地球表层系统最重要的景观标志,能直接反映出城市化过程中土地利用的演化特征,为城市规划、城市建设和区域可持续发展提供有效的决策支持。本文从土地利用的空间特征入手,以2010年和2014年两期区域土地利用数据为实例数据,通过计算图斑的平均最近邻指数,并采用核密度估计进行空间可视化表达与对比分析,研究城市土地利用时空格局的变化特征。结果显示,2010~2014年实验数据土地利用各地类的空间距离在缩小,整体呈现空间聚集趋势;而空间聚集的原因主要有图斑增量明显、原空间密集区向外扩展以及图斑空间位置改变等。 Spatial pattern of land use is the most important landscape mark of the earth surface, which can reflect the evolution characteristics of land use in urbanization process, and provide effective de-cision support for city planning, urban construction and regional sustainable development. From the perspective of spatial characteristics of land use, taking urban land use data in the year of 2010 and 2014 for examples, by means of calculating the average nearest neighbor index, and using kernel density estimation for spatial visualization and comparative analysis, urban land use evolution patterns are studied. As a result, spatial distance between the patches of experimental land use data is narrowed around 2010-2014, which presents the trend of space gathering. The main reasons of the spatial aggregation are the obviously increase of certain land use type, the outward expansion of original space’s intensive area, the change of space position and so on.
庄澜1,孟雯1*,杨光2
1上海市不动产登记事务中心,上海
2同济大学测绘与地理信息学院,上海
收稿日期:2016年5月24日;录用日期:2016年6月17日;发布日期:2016年6月20日
土地利用空间格局是地球表层系统最重要的景观标志,能直接反映出城市化过程中土地利用的演化特征,为城市规划、城市建设和区域可持续发展提供有效的决策支持。本文从土地利用的空间特征入手,以2010年和2014年两期区域土地利用数据为实例数据,通过计算图斑的平均最近邻指数,并采用核密度估计进行空间可视化表达与对比分析,研究城市土地利用时空格局的变化特征。结果显示,2010~2014年实验数据土地利用各地类的空间距离在缩小,整体呈现空间聚集趋势;而空间聚集的原因主要有图斑增量明显、原空间密集区向外扩展以及图斑空间位置改变等。
关键词 :土地利用,空间特征,空间分析,空间距离,空间格局
土地利用空间格局是地球表层系统最重要的景观标志,土地利用直接反映了人类与自然的相互影响与交互作用 [
从20世纪80年代中后期以来,随着全球化与可持续发展领域问题的日益严峻,土地利用/覆被变化研究受到国际社会的重视 [
本文从土地利用数据的空间特征入手,通过空间分析的方法分析土地利用数据的空间距离、空间格局的变化特点,研究城市土地利用的变化特征、变化趋势,并为城市土地利用、城市建设和区域可持续发展提供有效决策支持 [
空间分析是基于地理对象的位置和形态特征的空间数据分析技术,它将空间信息分为空间位置、空间分布、空间形态、空间距离、空间方位、拓扑、相似和相关 [
空间距离分析主要针对空间对象之间的远近程度进行测量,并可从多个尺度进行对比分析。本文采3用平均最近邻指数(Average Nearest Neighbor, ANN)来定量计算土地利用各地类图斑之间的空间距离及整体分布变化情况。
平均最近邻指数ANN算法公式如下:
其中,
而
公式(2)~(3)中,
平均最近邻指数表示的是实际图斑之间最小距离均值与随机分布下期望距离均值之间的比值,如果0 < ANN < 1,则表示土地利用图斑对象之间的实际最小距离小于随机分布状态,在空间上呈现聚集分布;若ANN > 1,则表示土地利用图斑对象之间的最小距离大于随机状态,呈现离散分布。通过该指数可以获取土地利用各地类之间的空间距离和空间分布信息。
空间格局分析主要定量描述图斑对象在具体空间上的分布模式,并基于分布情况进一步分析空间聚集或离散区域的土地利用结构及变化特征。
本文采用非参数核密度估计对土地利用数据进行空间分布分析,核密度估计可将数据对象的特征在任何尺度下进行空间可视化表达,具体采用Rosenblatt-Parzen全局核密度估计,公式如下:
其中
实例数据采用上海市某区域2010年和2014年两期土地利用现状数据进行对比分析,如图1所示。实例数据根据《全国土地分类(试行)》进行分类,每期数据都包含三大类所有三级类。为了便于计算与分析,本文分类截至到二级类,即一级类包含农用地、建设用地、未利用地;二级类包含耕地、园地、林地、牧草地、其他农用地、商服用地、工矿仓储用地、公共设施用地、公共建筑用地、住宅用地、交通用地、水利设施、特殊用地、未利用土地和其他土地。
针对2010年和2014年两期实例土地利用数据,采用平均最近邻指数统计各二级类图斑之间的最小空间距离特征及空间距离的变化趋势,对空间特征进行基本的数值统计;基于空间距离分析结果,再利
图1. 2010年(上图)和2014年(下图)土地利用实例数据
用核密度估计进行空间可视化表达,获取具体的空间格局特征与变化趋势,结合数值统计与空间统计的方法进一步研究土地利用变化特征及趋势。
根据公式(1)~(3),分别计算实例区域内15个二级类图斑数据的平均最近邻指数,结果如表1所示。
根据表1结果可知,2010年该区域土地利用中图斑之间平均空间距离最短的是住宅用地、耕地和其他农用地,分别为47.424米、65.398米和71.236米,而平均距离最大的是特殊用地、未利用土地和商服用类型,均超过200米;2014年该区域图斑之间平均距离最小的同样也是住宅用地、耕地和公用设施这三个类型,距离最大的是未利用土地,达到328.413米。在平均最近邻指数方面,2010年数据平均最近邻值最小的是园地、水利设施和住宅用地,均小于0.45,地类呈现空间聚集分布模式,而其他用地、耕地和其他农用地的ANN值较高,表明这些地类图斑间距较大,且分布在区域外围;2014年计算结果与2010年基本保持一致。
图2表示2010~2014年土地利用空间距离变化情况,结果显示,图斑之间平均最近距离减少最大的是特殊用地、商服用地和公共建筑用地类型,图斑距离缩小了187.27米、93米和79米,而住宅用地、林地和耕地之间的距离变化差异小,均在1.5米以内;在平均最近邻指数差值方面,变化较大的特殊用地、公共建筑用地和交通运输用地,而住宅用地差值大于0。综述表1和图2,结果表明2010~2014年间区域土地利用数据整体呈现图斑距离减少且ANN指数下降,在空间上呈现聚集趋势,在地类方面,特殊用地、商服、公共建筑用地以及交通运输在过去四年间因新增用地或位置变化而在空间上呈现聚集趋势。
土地利用地类 | 平均最近邻指数 | ||||||
---|---|---|---|---|---|---|---|
Do_2010 (米) | Do_2014 (米) | 差值 (米) | ANN_2010 | ANN_2014 | 差值 | ||
1. 农用地 | 耕地 | 65.398 | 63.908 | −1.49 | 0.684 | 0.679 | −0.005 |
园地 | 101.521 | 96.128 | −5.393 | 0.314 | 0.298 | −0.016 | |
林地 | 77.002 | 75.969 | −1.033 | 0.573 | 0.562 | −0.011 | |
牧草地 | 0 | 0 | 0 | 0 | 0 | 0 | |
其他农用地 | 71.236 | 68.241 | −2.995 | 0.663 | 0.641 | −0.022 | |
2. 建设用地 | 商服用地 | 282.976 | 189.417 | −93.559 | 0.519 | 0.485 | −0.034 |
工矿仓储用地 | 119.733 | 107.581 | −12.152 | 0.587 | 0.575 | −0.012 | |
公用设施用地 | 92.518 | 63.917 | −28.601 | 0.502 | 0.447 | −0.055 | |
公共建筑用地 | 239.827 | 160.782 | −79.045 | 0.472 | 0.383 | −0.089 | |
住宅用地 | 47.424 | 46.594 | −0.83 | 0.449 | 0.451 | 0.002 | |
交通运输 | 116.965 | 80.588 | −36.377 | 0.528 | 0.468 | −0.06 | |
水利设施 | 185.231 | 140.829 | −44.402 | 0.418 | 0.361 | −0.057 | |
特殊用地 | 504.387 | 317.117 | −187.27 | 0.526 | 0.399 | −0.127 | |
3. 未利用地 | 未利用土地 | 345.36 | 328.413 | −16.947 | 0.55 | 0.535 | −0.015 |
其他用地 | 138.316 | 130.771 | −7.545 | 0.75 | 0.72 | −0.03 |
表1. 平均最近邻指数结果
图2. 2010~2014年实例数据平均最近邻差值统计图
根据基于平均最近邻的空间距离分析,可以获取实验区域土地利用地类在2010~2014年间的空间距离变化信息。然而仅仅通过数值表达的一维信息无法全面分析土地利用的变化特征,如商服、公共建筑、交通和特殊用地在2010~2014年间的ANN指数波动较大,这四个地类只能通过空间距离变化分析整体趋势,无法对不同地类内部的空间分布模式进行变化分析。因此,本节将通过核密度估计进行空间可视化表达,对土地利用的空间格局变化进行进一步的对比分析,全面了解土地利用的变化情况。
1) 商服用地
图3为2010年和2014年两期商服用地空间分布对比,结合空间距离分析结果可知:2010~2014年实验区域内商服用地平均最近邻距离下降了93米,下降幅度达33%,主要由于其图斑数量较2010年有明显的新增,导致密度高值点数也有明显的增加,使得商服用地空间格局更加密集,主要的增量分布于实验区域的东北角,同时在四年间呈现以区域中心向周围扩散的变化特征。
2) 公共建筑用地
图4为2010年和2014年两期公共建筑用地空间分布对比,结合空间距离分析结果可知:2010~2014年公共建筑用地图斑最近距离减少了79米,ANN指数两期下降超过18%,如图4所示,2014年在空间上新增了多个密度高值点且原有高值点面积扩大,与商服用地相似,空间格局的变化特征以中心区域向外扩展,但新增密集点数或新增图斑面积不明显。
3) 交通运输用地
图5为2010年和2014年两期交通用地空间分布对比,根据表1可知2010~2014年交通用地图斑之间最近距离减少了36米,ANN下降幅度超过11%;在空间格局方面,四年间实验区域内主要变化在于东北角的密集区域增大且成片,而南部及西部交通用地密集程度较低,空间分布较稀疏离散。
4) 特殊用地
根据表1结果,特殊用地在2010~2014年间的图斑空间最近距离减少量最大为187米,下降幅度超过37%,是所有土地利用二级类中距离减少量和幅度最大的类型。如图6所示,在空间上,特殊用地的空间格局主要分布于实验区域边缘,因此造成该地类图斑之间的平均距离绝对值大,而2014年特殊用地整体没有明显的增量,即空间加密过程,在北部有新增高密集值点,且与原有密度高值点距离相近,使其产生局部高度聚集的效果。特殊用地类型四年间空间分布模式不变,增量不明显,两期空间聚集趋势明显主要由于图斑空间位置的变化形成区域高度聚集,整体空间格局较为稳定。
图3. 商服用地空间分布(上图2010年,下图2014年)
图4. 公共建筑用地空间分布(上图2010年,下图2014年)
图5. 交通用地空间分布(上图2010年,下图2014年)
图6. 特殊用地空间分布(上图2010年,下图2014年)
本文基于土地利用的空间特征,通过对图斑数据进行空间距离分析,并采用核密度估计进行空间可视化表达,对比分析2010年和2014年两期实验数据的土地利用变化情况,并针对变化较大的商服、公共建筑、交通和特殊用地进行空间分布模式分析,获得以下结论:
1) 2010~2014年实验数据土地利用各地类间的空间距离在缩小,整体呈现空间聚集趋势;
2) 各地类空间聚集原因不同,其中商服用地面积有明显新增,在空间上形成明显的加密点,空间格局以东北角和中心区域密集分布;公共建筑和交通用地没有明显的空间加密过程,空间格局主要以原密集区面积增大且向外扩展形成;特殊用地空间分布模式不变,由于原始分布于实验区外围,空间位置的变化使得区域尺度下产生明显高度聚集,但空间格局较为稳定。
上海市规划和国土资源管理局科研项目。
庄澜,孟雯,杨光. 基于空间分析的城市土地利用变化特征研究 The Study on Urban Land Use Evolution Based on Spatial Analysis Method[J]. 测绘科学技术, 2016, 04(03): 69-79. http://dx.doi.org/10.12677/GST.2016.43009