一种综合人口密度与景观结构的城市功能区的识别方法与流程

本方法属于景观生态学、城市地理学、地理信息技术交叉的
技术领域：
。
背景技术：
：伴随着城市化的进程，沿城市拓展轴必然会出现不同的城市功能分区，如城市中心区、外围区、（隔离带）、（卫星城）、远郊区等。目前使用土地利用数据多用于区分城市建成区与非建成区，对这些功能分区之间的界限研究较少。此外，有关城市内部功能格局的研究主要是利用各区县分行业人口或社会经济相关指标对城市内部的功能分工特征进行分析，较少从景观角度进行研究的。综合人口密度和土地利用数据来识别城市功能区的应用研究相对较少。技术实现要素：本发明的目的是利用人口密度和景观结构，在城市拓展轴上识别城市功能区。为城市规划和城市功能的优化配置提供科学依据。具体方法包括以下步骤：a.沿城市拓展轴设立样带与样方城市拓展轴浓缩了城市化的发展历程，代表了城市发展的方向和广度。因此沿城市拓展轴设立相应的样带和样方可以反映城市功能的分布和城市发展的程度。样方的面积大小尽可能相似，沿城市拓展轴方向分别给样方赋予不同的ID号。根据土地利用数据建立城市功能景观体系，如图2。提取样方的人口密度与土地利用类型信息将样带的矢量图层与村级行政区划图做叠加分析，提取出样带范围内的村级行政区范围，同时基于各个行政村的人口数据和面积，得到各个行政村的人口密度，在此基础上计算各行政村与每个样方重叠部分的人口，从而相加得到各样方内的总人口，除以样方面积进而得到各样方的人口密度。计算公式如下：Pi代表样方i的人口密度，其中Aj表示第j个行政村的人口，Mj表示第j个行政村的面积，Bij表示第j个行政村落在第i个样方的面积，Ti表示样方i的面积。将样带的矢量图层与土地利用类型图层做叠加分析，提取出样带范围内的土地利用信息，将居住用地赋予居住功能景观属性，对应代码为1；将工业用地和物流仓持用地赋予工业功能景观属性，对应代码为2；将商服用地赋予商业功能景观属性，对应代码为3；将公共设施用地、公共管理与公共服务用地、科教用地、文教体卫用地和各类交通运输用地赋予公共服务功能景观属性，对应代码为4；将绿地赋予生态服务功能景观属性，对应代码为5。将耕地、园地、农村宅基地、空闲地和裸地赋予潜在城市景观属性，对应代码为6；将林地、水体和滩涂赋予其他景观属性，对应代码为7。并且将每个样方的城市功能景观（代码1-5）分别生成一个图层。计算城市功能景观的百分比结构和景观格局指数对每个样方的城市功能景观和非城市功能景观分别统计面积，计算城市功能景观的百分比结构，包括城市景观和非城市景观分别占样方的比例，城市景观各类型占城市景观的比例以及非城市景观各类型占非城市景观的比例。对每个样方的城市功能景观矢量图按照转化为1m*1m的栅格图，计算景观破碎度指数、香农多样性指数和景观分离度指数。1、景观破碎度指数：FN代表景观破碎度指数，其中NP表示景观的图斑个数，TA是景观面积。2、香农多样性指数：SHDI代表香农多样性指数，其中pi表示景观中斑块类型i的面积比重。3、景观分离度指数：LDI代表景观分离度指数，其中aij为斑块ij的面积，A为整个景观的面积。综合人口密度与景观结构来识别城市功能分区将人口密度、景观百分比结构和景观格局指数作为聚类分析的因子，采用系统聚类分析，自动对样方进行分类，识别出沿城市拓展轴的城市功能区。附图说明图1是样带与样方分布图图2是城市功能景观与土地利用类型对应表图3是人口密度空间分布图图4是聚类结果与城市功能分区示意图。具体实施方式本发明一种基于人口密度与景观结构的城市功能分区的方法可以通过以下实施例进行说明：a.本发明选择厦门市为实施案例城市，沿城市拓展轴（厦禾路-嘉禾路-厦门大桥-同集南路-同集中路-同集北路）做半径1km的缓冲区作为样方，在厦禾路起点布设第一个点，之后每2km布点，过点做垂线分割缓冲区，依次划分为17个样方。结果如图1。b.基于土地利用对城市功能景观进行分类，结果如图2。c.计算样方的人口密度基于厦门市各个行政村的人口数据和面积，得到各个行政村的人口密度，将各行政村与样方进行叠加，计算各行政村与样方重叠部分的人口，从而相加得到各样方内的总人口，除以样方面积进而得到各样方的人口密度。结果如图3。d.计算城市功能景观的百分比结构和景观格局指数计算城市功能景观的百分比结构，包括城市景观和非城市景观分别占样方的比例，城市景观各类型（居住功能景观、工业功能景观、商业功能景观、公共服务功能景观和生态服务功能景观）占城市景观的比例以及非城市景观各类型（潜在城市景观和其他景观）占非城市景观的比例。结果如下表：道路ID号城市景观占样方的比例非城市景观占样方的比例居住功能景观占城市景观的比例工业功能景观占城市景观的比例商业功能景观占城市景观的比例公共服务功能景观占城市景观的比例生态功能景观占城市景观的比例潜在城市景观占非城市景观的比例其他景观占非城市景观的比例192.98%7.02%55.23%0.00%9.92%26.26%8.59%47.98%52.02%283.08%16.92%60.94%1.14%7.03%25.35%5.53%48.61%51.39%378.86%21.14%55.00%4.88%18.01%21.31%0.80%30.79%69.21%493.22%6.78%53.45%1.70%12.62%26.00%6.24%36.52%63.48%577.55%22.45%56.92%3.40%11.88%23.88%3.93%49.34%50.66%679.93%20.07%16.44%38.01%11.55%27.66%6.34%96.03%3.97%777.78%22.22%5.68%39.08%6.39%43.17%5.69%99.67%0.33%854.96%45.04%9.17%32.16%0.04%55.77%2.85%72.73%27.27%952.33%47.67%16.15%0.00%0.00%83.85%0.00%25.27%74.73%1056.27%43.73%17.04%3.32%5.07%74.22%0.35%82.82%17.18%1155.10%44.90%15.38%46.11%0.00%38.45%0.06%83.94%16.06%1242.58%57.42%23.78%31.68%3.40%30.85%10.28%38.00%62.00%1326.35%73.65%0.00%41.02%0.69%44.62%13.67%45.09%54.91%1453.45%46.55%0.00%49.24%0.00%46.82%3.93%30.02%69.98%1569.74%30.26%0.67%58.43%1.27%37.96%1.68%88.31%11.69%1627.24%72.76%0.00%17.63%2.05%74.44%5.89%88.40%11.60%1742.70%57.30%0.44%50.67%1.22%45.09%2.58%92.04%7.96%计算景观格局指数（香农多样性指数、景观破碎度指数和景观分离度指数），结果如下表：道路ID号香农多样性指数景观破碎度指数景观分离度指数road11.120.200.92road21.050.220.92road31.150.210.93road41.190.250.92road51.160.270.93road61.440.220.88road71.230.210.80road81.010.250.77road90.440.090.39road100.810.150.46road111.020.230.84road121.420.250.86road131.030.240.75road140.830.110.79road150.840.210.89road160.770.340.68road170.880.360.79road181.070.480.60e.综合人口密度与景观结构来识别城市功能分区将人口密度、景观百分比结构和景观格局指数作为聚类分析的因子，采用系统聚类分析，自动对样方进行分类，分类结果与城市化演变规律相一致，根据城市化演变规律将聚类结果分别对应为城市中心区、外围区、隔离带、卫星城和远郊区。结果如图4。当前第1页1 2 3