一种暴露与非暴露绿地识别方法、终端设备及存储介质技术

本发明专利技术涉及一种暴露与非暴露绿地识别方法、终端设备及存储介质，该方法中包括：S1：获取研究区内的原始手机信令数据；S2：将原始手机信令数据转换成具有空间坐标的手机信令数据；S3：获取研究区内的城市绿地数据；S4：将具有空间坐标的手机信令数据转换成与城市绿地数据对应的图层坐标系一致的数据；S5：将转换后的坐标系一致手机信令数据与城市绿地数据进行空间叠加分析，得到叠加图层，基于叠加图层得到研究区内的暴露绿地和非暴露绿地。本发明专利技术改进了现有技术中未能充分反映人群实际绿地暴露的不足，实现了对人群绿地暴露的科学定量测度。量测度。量测度。

【技术实现步骤摘要】
tool”工具箱，找到“Projection and transformation”功能下的“Define projection”，在“Define projection”对话框中选择“Geographic coordinate system”，再选择“World”，最后选择“WGS 1984”，即可得到具有空点坐标的手机信令数据。
[0014]进一步的，步骤S4中将具有空间坐标的手机信令数据转换成与城市绿地数据对应的图层坐标系一致的数据的方法为：在ArcGIS软件中运用Arctoolbox中的“Data management tool”工具箱，找到“Projection and transformation”功能下的“Projection”，在“Projection”对话框中选择与城市绿地数据对应的图层一致的坐标系，即可得到与城市绿地数据对应的图层坐标系一致的数据。
[0015]进一步的，还包括S6：统计研究区内的暴露绿地和非暴露绿地的面积。
[0016]一种暴露与非暴露绿地识别终端设备，包括处理器、存储器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本专利技术实施例上述的方法的步骤。
[0017]一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现本专利技术实施例上述的方法的步骤。
[0018]本专利技术采用如上技术方案，改进了现有技术中未能充分反映人群实际绿地暴露的不足，实现了对人群绿地暴露的科学定量测度。
附图说明
[0019]图1所示为本专利技术实施例一的流程图。
[0020]图2所示为本专利技术实施例一中基站位置图。
[0021]图3所示为本专利技术实施例一中研究区城市绿地图。
[0022]图4所示为本专利技术实施例一中研究区城市绿地和手机信令数据叠加过程图。
[0023]图5所示为本专利技术实施例一中暴露绿地和非暴露绿地提取过程图。
[0024]图6所示为本专利技术实施例一中研究区暴露绿地和非暴露绿地分布图。
具体实施方式
[0025]为进一步说明各实施例，本专利技术提供有附图。这些附图为本专利技术揭露内容的一部分，其主要用以说明实施例，并可配合说明书的相关描述来解释实施例的运作原理。配合参考这些内容，本领域普通技术人员应能理解其他可能的实施方式以及本专利技术的优点。
[0026]现结合附图和具体实施方式对本专利技术进一步说明。
[0027]实施例一：
[0028]本专利技术实施例提供了一种暴露与非暴露绿地识别方法，如图1所示，所述方法包括以下步骤：
[0029]S1：获取研究区内的原始手机信令数据。
[0030]手机信令数据可以从三大电信运营商(移动，电信和联通)获取。手机信令数据可以用于探索居民的出行行为特征，如挖掘人群的出行规律和出行特征、研究居民的通勤行为以及预测人群的出行路线等。
[0031]该实施例从某运营商获取了2020年6月的研究区范围内用户的手机信令脱敏数据(手机信令数据)。借助于该运营商的基站设备，对用户的手机信令数据进行采集，以此反映
对用户的个人行为信息(如图2所示)。从图2中可以看出，手机信令数据是基于每个格网为250m
×
250m采集的点数据，该数据包含了每个基站2020年6月1日
‑
2020年6月30日的逐小时数据。每条手机信令数据仅包括时间、性别、年龄、位置以及人口扩样等信息，个人身份等涉及隐私的信息已隐去。
[0032]S2：将原始手机信令数据转换成具有空间坐标的手机信令数据。
[0033]该实施例中采用的转换方法为：将原始手机信令数据导入ArcGIS软件，运用Arctoolbox中的“Data management tool”工具箱，找到“Projection and transformation”功能下的“Define projection”，在“Define projection”对话框中选择“Geographic coordinate system”，再选择“World”，最后选择“WGS 1984”，即可得到具有空点坐标的手机信令数据。
[0034]S3：获取研究区内的城市绿地数据。
[0035]城市绿地数据可以基于卫星遥感数据(Landsat 8或者高分一号等)通过植被指数法提取。归一化植被指数(Normalized Difference Vegetation Index,NDVI)能够反映土地植被覆盖的情况，是一种常用的遥感指标。该实施例中基于高分一号(GF
‑
1)的遥感影像数据进行提取城市绿地数据，运用植被指数法提取了研究区的绿地信息。基于归一化植被指数的结果，结合经验值设置阈值以提取研究区的植被覆盖信息，即研究区城市绿地，如图3所示。
[0036]S4：将具有空间坐标的手机信令数据转换成与城市绿地数据对应的图层坐标系一致的数据。
[0037]该实施例中采用的转换方法为：在ArcGIS软件中运用Arctoolbox中的“Data management tool”工具箱，找到“Projection and transformation”功能下的“Projection”，在“Projection”对话框中选择与城市绿地数据对应的图层一致的坐标系，即可得到与城市绿地数据对应的图层坐标系一致的数据。
[0038]S5：将转换后的坐标系一致手机信令数据与城市绿地数据进行空间叠加分析，得到叠加图层，基于叠加图层得到研究区内的暴露绿地和非暴露绿地。
[0039]本实施例中空间叠加分析基于ArcGIS软件实现，即在ArcGIS里运用intersect功能将基站矢量数据与研究区城市绿地矢量数据进行叠加，得到“手机信令intersect城市绿地”图层(即叠加图层)，如图4所示。对叠加图层的属性表进行属性选择，其中“pop＝NULL”的绿地为无人暴露于其中的绿地，是非暴露绿地；反之，有人群暴露于其中的绿地为暴露绿地，如图5所示。
[0040]为了得到暴露绿地和非暴露绿地的准确数据，本实施例中还包括S6：统计研究区内的暴露绿地和非暴露绿地的面积。
[0041]具体统计方法为：在ArcGIS软件中，打开叠加图层的属性表，选择“Add Field”，“Name”填“Area”,“Type”选择“Float”，即在属性表里建立了area字段。对叠加图层的属性表进行属性选择，运用“pop＝NULL”选择条件，选择出非暴露绿地，在选择后的属性表中，找到“Area”字段，右键单击，选择“Calculate Geometry”功能，“Property”选择“Area”，“Units”选择“Square Kilometers[sq km]”，即可计算出各个城市非暴露绿地斑块的面积。对“Area”列属性右键单击，选择“Statistics”，其中“Sum”项为整个城市非暴露绿地的面积。同理，对叠加图层的属性表进行属性选择，运用“pop本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种暴露与非暴露绿地识别方法，其特征在于，包括以下步骤：S1：获取研究区内的原始手机信令数据；S2：将原始手机信令数据转换成具有空间坐标的手机信令数据；S3：获取研究区内的城市绿地数据；S4：将具有空间坐标的手机信令数据转换成与城市绿地数据对应的图层坐标系一致的数据；S5：将转换后的坐标系一致手机信令数据与城市绿地数据进行空间叠加分析，得到叠加图层，基于叠加图层得到研究区内的暴露绿地和非暴露绿地。2.根据权利要求1所述的暴露与非暴露绿地识别方法，其特征在于：步骤S2中将原始手机信令数据转换成具有空间坐标的手机信令数据的方法为：将原始手机信令数据导入ArcGIS软件，运用Arctoolbox中的“Data management tool”工具箱，找到“Projection and transformation”功能下的“Define projection”，在“Define projection”对话框中选择“Geographic coordinate system”，再选择“World”，最后选择“WGS 1984”，即可得到具有空点坐标的手机信...

【专利技术属性】
技术研发人员：吝涛，姚霞，张浚茂，陈媛，李方芳，
申请(专利权)人：中国科学院城市环境研究所，
类型：发明
国别省市：