基于建筑物与POI数据的人口空间化研究.pdf

书书书第 卷第 期（）年 月中 国 地 震 李金香，谭明，古力孜帕木拉提，等 基于建筑物与 数据的人口空间化研究 中国地震，（）：基于建筑物与 数据的人口空间化研究李金香谭明古力孜帕 木拉提李波张金燕新疆维吾尔自治区地震局，乌鲁木齐 摘要人口数据是地震灾害人员伤亡评估的基础，准确的人口空间分布信息对提高应急灾情评估的准确度和开展高效的应急决策具有重要意义。本文结合建筑物空间分布、（兴趣点）、路网与人口统计等数据，基于 相关分析以及多元线性回归分析方法，构建人口空间化模型，实现了人口统计数据基于规则格网的更为精确的可视化表达。结果表明：建筑物空间分布能较好地反映人口的宏观分布特征，但是对于人口分布细节特征的刻画，则是以 数据来反映人口的微观分布特征效果更为理想；受区域经济发展影响，人口分布具有明显的空间差异，由城镇向乡村区域递减的趋势明显，中心城区、周边城镇、偏远乡镇之间的人口密度差异巨大。基于建筑物与 数据的人口模拟值与实际值的偏差较小，数据精度符合实际情况，可为灾情研判提供可靠的基础数据。关键词：人口空间化 多元线性回归千米网格 文章编号 （）中图分类号 文献标识码 收稿日期 修定日期 项目类别 新疆维吾尔自治区自然科学基金面上项目（）、新疆地震科学基金项目（）共同资助 作者简介 李金香，女，年生，高级工程师，主要从事地震应急与遥感震害监测研究。：引言新疆地区地震频发，给新疆人民造成了巨大的人员伤亡和经济损失。在地震应急救援工作中，以行政区划为单元的人口统计数据在使用中存在的空间分辨率低、直观性差、不支持空间运算和分析等不足亟待解决（高义等，）。无论是震前预测还是震后快速评估，准确的人口空间分布信息是人员伤亡评估的基础（徐敬海等，；谢江丽等，；周中红等，；朱鹏宇等，），以往地震应急工作中以乡镇级行政区划的人口平均密度或面积权重法计算的人口数量，与实际调查结果间存在较大误差。尤其是新疆人口的不均匀分布，使得这种方法获取的数据精度不高。而人口数据空间化能够更为客观地展示人口分布情况，是解决该问题的有效手段。近年来，得益于遥感（）和地理信息系统（）技术的迅速发展，统计数据格网化得到了快速发展。人口数据空间化一直是统计数据空间化的主要研究方向（董南等，）。王雪梅等（）在回顾了国内外人口空间化研究的主要方法后，认为国外的研究主要包括从遥感解译信息反演人口数据、从 夜间灯光数据反演人口数据和从遥感直接获取的光谱特征直接反演人口数据几个方面。而国内的研究主要是根据土地利用数据和其他地理因子（如高程、道路、居民区等）建立回归模型来反演人口数据（柏中强等，；王珂靖 期李金香等：基于建筑物与 数据的人口空间化研究等，；崔晓临等，；杜培培等，）。随着人口数据空间化研究的深入，形成了一系列具有代表性的模型和方法，主要有核密度估计法（孙艳萍等，）、面积内插法（吕安民等，）、土地利用影响模型（江东等，；黄河清等，；唐奇等，）以及多源数据融合模型（董春等，；王春菊等，；淳锦等，；赵鑫等，；王晓洁等，；王芳等，）。已有研究大多基于遥感夜间灯光数据、土地利用数据等（高倩等，；王明明等，；于婷婷，），并采用模型算法进行估算，是全局性的建模分析，未能完全地考虑微观尺度上人口分布的随机性，不能反映内部差异，缺少对人口空间分布异质性、非平稳性的研究，特别是针对以县级区划为基本研究单元的理论研究。随着时代发展，目前房屋建筑空间分布数据越来越精细，（）点数据越来越详实，基于房屋真实空间分布数据进行人口统计数据空间化会大大提升千米格网数据的准确性，为新疆地区地震应急辅助决策提供更加科学准确的数据支撑。图 研究区概况本研究以新疆“基于遥感影像和经验估计的区域房屋震害风险初判”项目的示范区 新疆巴楚县、库车市和乌鲁木齐县为研究区，以规则格网为研究单元，利用建筑物空间分布数据、数据、路网数据、人口统计数据等数据源，进行人口空间化算法模型研究，以期获得空间精度和写实程度能够满足地震应急应用要求且使用方便灵活的格网化空间数据，提升地震灾情评估速度和准确度。研究区概况与数据处理 研究区概况巴楚县位于天山南麓，塔里木盆地和塔克拉玛干沙漠边缘，辖 镇 乡，处于巴楚隆起大地构造单元（图 （）。库车市位于天山中部南麓，塔里木盆地北缘，辖 镇 乡 街道，地势北高 南 低、自 西 北 向 东 南 倾 斜，地 貌 分 为 北 部 山 地、中 部 戈 壁 和 南 部 冲 积 平 原（图 （）。乌鲁木齐县南依天山，北与准噶尔盆地相连，地势东南高，西北低，辖 镇 乡（图 （）。巴楚县、库车市与乌鲁木齐县历史上地震频发，人口分布却极不均匀，三县人口均集中分布在绿洲区域，中心城区人口密集，经济发展迅速，交通设施、购物、教育医疗等服务机构完善，山地、戈壁及沙漠地区人烟稀少，各乡之间人口差异明显。本文选取三县为研究区域，具有一定的研究意义。中国地震 卷 数据处理与技术路线数据中的人口统计数据包含乡镇人口统计数据；行政区划数据来自新疆地震应急基础数据库，包含乡镇行政区划面矢量数据，存储格式为 格式，字段属性包含行政区划代码、行政区划名称等；交通路网空间数据来自新疆地震应急基础数据库，包括国道、省道、县道、乡道、专用道、城市内部道路等；房屋空间分布数据通过“基于遥感影像和经验估计的区域房屋震害风险初判”项目获取；电子地图兴趣点来自高德导航数据，兴趣点包括学校、超市、医院等点位数据。数据预处理主要包括：空间数据配准，完成各类空间数据位置校正；人口统计数据与乡镇街道行政区划数据完成关联；对 数据进行清洗与剪切，进而进行分类整理。数据是一种代表真实地理实体的点状地理空间大数据，数据量十分庞大，涉及人们生产生活的各个方面，发掘其中所蕴含的内在信息关系，提取有用的内容，是当前的研究热点之一。目前常见的电子地图均包含 ，且都具有名称、类别和位置等主要属性，为了方便查询，需对其进行分类。的分类体系不可能将包罗万象的所有信息进行精准分类，只能将基础和普遍的信息进行分类，以满足大众的基本需求。在此条件下，分类代码体系的编制应遵循一定的原则。高德软件公司与天地图有限公司、北京四维图新科技股份有限公司、中国物品编码中心等多家单位共同参与完成了中华人民共和国国家标准 地理信息兴趣点分类与编码 的起草。高德软件 编码规则符合国标标准，采用线分类法将 分为大类、中类、小类三个层次，其中依据 使用的普遍性和社会公众对于 的关注程度划分大类，每一大类按照其不同特点和相互之间的内在联系划分中类和小类。本文的 数据来自高德导航，数据在获取的同时配有高德数据分类代码表，根据分类代码，对 数据进行分类整理（表 ）。最后根据研究区域经济发展情况进行分区。在对研究区进行分区后，计算各区域内建筑物面积、道路密度、各类型 数量等，与人口数据进行相关性分析，选择与人口数据存在显著相关的因子进行多元线性回归，建立各分区的人口空间数据集。技术路线如图 所示。研究方法本文以建筑物空间分布数据、数据、路网数据、人口统计数据等数据源为基础，进行人口空间化算法模型研究（图 ）。建筑物是人们生产生活的房屋及其附属设施，是最直接反映人口分布的因子（刘焕金，）。道路对人口分布具有指示作用，然而新疆地广人稀，人口在绿洲区集中分布，很多乡镇面积大、道路长度长而人口稀少，因此采用道路长度进行人口分布研究存在偏差，故本文引进路网密度与居民区道路长度两个因子与人口分布进行相关性研究。是指具有地理标识的空间特征物，包含名称、类别、经纬度等信息，具有易获取、现势性强、数据量丰富、定位精度高、更能反映微观细节信息等特点（，），其能够在一定程度上直观反映人口的空间分布（，）。本文将建筑物数据（）、路网密度（）、居民区道路长度（）、大厦小区（个）、政府机关（个）、餐饮住宿（个）、购物（个）、交通服务（个）、教育医疗（个）、公共服务（个）、商业机构（个）、文体休闲（个）、观光旅游（个）、农牧场点数（个）与人口统计（人）数据进行 相关性分析，利用与人口分布相关性较高的因子进行人口空间化建模。在统计学 期李金香等：基于建筑物与 数据的人口空间化研究表 部分 数据分类大类中类小类大类中类小类大厦小区大厦大厦购物商场 购物商场 购物中心超市超市小区住宅小区便利店便利店市场 商店专营市场 商店政府机关政府机关 基层群众自治组织中央机关公共服务电信邮政电信省级机关邮政市级 区县级机关传媒机构广播电视乡镇级机关 街道办事处新闻出版村级单位福利机构福利机构其他机关单位事务所事务所消防机关消防机关鉴定所司法机关公安局 派出所公共厕所公共厕所检察院 法院其他专业服务美容美发交警队旅行社司法局 司法处中介机构公证机关公证局 公证处保洁 家政 搬家海关海关装潢 设计 广告涉外机构大使馆 领事馆婚纱 摄影 冲印国际组织机构物流 货运 快递驻地机构驻地机构殡仪社会团体社会团体其他服务中，相关性分析是评价 个统计变量相关性的一种指标，其对原始变量的分布不做要求，适 用 范 围 较 广。其 主 要 思 想 是：分 别 对 个 变 量 、做 等 级 变 换（），用等级 和 表示；然后按 相关性分析的方法计算 和 的相关性。相关性分析公式为 （）（）（）（）式中，为第 个乡镇街道的统计人口数，为第 个乡镇街道第 类指标因子的数值，（）与 （）分别为对应的元素 和 在各自列向量中的排名，为乡镇街道总数，为第 类因子的数值与乡镇街道人口数的 相关系数。经过计算，选择相关性显著的因子作为建模因子。确定建模因子后，采用多元线性回归分析方法进行人口空间化模型构建。因经济发展不同，不同行政单元内单位建筑面积的人口数各不相同，为提高建模精度，将人口分布特征较为接近的地区划分为一类进行分区建模，因模型中的系数具有物理意义，需保证系数为非负数。中国地震 卷图 技术路线图 研究框架示意图采用多元线性回归模型进行数据拟合，以各建模因子作为回归分析的自变量，以各乡镇统计人口数据作为因变量，将回归置信度设置为 ，建立回归方程组公式，即 ，（）式中，为第 类建模因子的回归系数，为方程的截距。在建模时将常数项设置为 ，表明人口分布在与其强相关的因子分布的区域，且建立的模型回归系数应满足变量的显著性检验（、检验）。期李金香等：基于建筑物与 数据的人口空间化研究 实例验证 建模区划分及相关性分析建筑物作为人类生产和生活的载体，在很大程度上影响着人口的分布和迁移，可以说人口的分布都是围绕着建筑物展开的，尤其是居住建筑。因此，利用地理空间技术，以建筑物为离散的载体研究人口的空间化具有重要意义。本文更换了众多研究中采用的高程、坡度等影响因素，以人们生活的建筑物为载体，实现人口数据的精细空间化。受地理位置、经济发展等多因素影响，不同地区单位建筑物占地面积、人口密度各不相同，为提高模型构建的精度，对影响因素差异大的区域进行分区，以凸显影响因子的差异性。本研究以乡镇街道为最小研究单元，依据建筑物空间分布、兴趣点密集程度、道路密度等进行分区，将研究区分为个分区。第一类分区为靠近中心城区的街道及城镇，该类分区兴趣点密集，经济发达。第二类分区以农村居民点为主，农村居民点、兴趣点分散在各个乡镇。城镇建筑物存在较多的多层建筑，且各街道经济发展水平不同，多层建筑物的数量和楼层数差异也较大。假设人口数量总是对应一定的建筑面积（刘正廉等，），以各县乡镇人均建筑物占地面积为基准，对城镇区域建筑物进行平均楼层数设定，设城镇街道建筑物占地面积为 ，平均楼层数设定系数为，则城镇区域建筑物面积为 ，该方法可以使得城镇人均建筑面积更接近该县人均建筑面积，满足人口数量总是对应一定的建筑面积的假设。进而进行基于建筑物空间分布的人口空间化研究。对第一类分区和第二类分区的样本数据进行定量分析，统计各乡镇街道中各类型兴趣点的数量。分别计算第一类、第二类分区各乡镇建筑物数据、路网密度、居民区道路长度、大厦小区、政府机关、餐饮住宿、购物、交通服务、教育医疗、公共服务、商业机构、文体休闲、观光旅游、农牧场点的数量，并与人口统计数据进行 相关性分析。相关性分析结果（表 ）显示，第一类分区中建筑物建筑面积、路网密度、购物、公共服务与人口存在显著正相关，相关性系数分别为 、；第二类分区中建筑物建筑面积、居民区道路长度、政府机关与人口存在显著正相关，相关性系数分别为 、。将存在显著相关的因子作为各分区建模因子，对 个分区进行模型构建。基于建模因子的多元线性回归分析根据各分区建模因子与人口的相关性，采用多元线性回归分析方法建立不同分区回归模型。以各分区建模因子数值为自变量，各乡镇街道人口数值为因变量，利用公式（）重新构建多元线性回归模型。多元线性回归要求因子间相互独立，经计算第一分区中购物与公共服务间的 相关系数结果为 ，存在显著正相关，购物与公共服务不是相互独立的因子，均参与回归会使模型产生多重共线性，故两者间选择一个因子进行回归分析，通常选择相关性高的因子。本研究中购物、公共服务与人口的相关程度相当，相关系数差仅为 ，对比两类 点的数据量及分布范围，购物 数量远多于公共服务，且分布更广，能够更好地指示人口分布，故最终选择购物因子进行回归分析。因此，选择购物、建筑物建筑面积、路网密度三个因子进行第一分区的回归分析，选择建筑物建筑面积、居民区道路长度、政府机关三个因子进行第二分区的回归分析，得到的回归方程为 第一分区：，（）中国地震 卷表 第一、第二分区人口与各因子的 相关系数因子类别第一分区 相关系数第二分区 相关系数建筑物建筑面积 路网密度 居民区道路长度 大厦小区 政府机关 餐饮住宿 购物 交通服务 教育医疗 公共服务 商业机构 文体休闲 农牧场 观光旅游 注：代表检验统计量在 水平上显著相关（双侧检验），代表检验统计量在 水平上显著相关（双侧检验）。第二分区：，（）式中，为第 个乡镇建筑物占地面积，为第 个乡镇路网密度，为第 个乡镇居民区道路长度，为第 个乡镇政府机关数量，为第 个乡镇购物点的数量，为第一、第二分区中第 个城镇街道的建筑物建筑面积拟合系数，其中第二分区中设定 为 。在 置信度下，第一分区和第二分区模型具有统计显著性，且建立的模型回归系数均满足变量的显著性检验（、检验）。人口空间化及精度评价 人口空间化人口空间化是将统计数据分布到其对应地理空间上的一个过程。根据分区结果，利用栅格计算方式，将统计人口数据转换为更能反映人口空间分布的栅格数据，完成人口数据空间化。栅格大小选择 规则格网，首先利用 软件中的 工具创建研究区域内规则格网的矢量数据，用规则格网裁剪研究区建筑物数据、路网数据，统计每个格网内的建筑物占地面积、路网长度、各分类 点数量等。根据分区方式，将格网数据分为第一分区格网和第二分区格网，分别计算格网内建筑物面积及路网密度，基于式（）、式（），计算各分区内每个格网的人口数。对分区人口的初始栅格数据进行拼接，得到整个研究区的人口栅格数据。精度评价将每个格网的模拟人口数汇总到乡镇街道行政单元上，采用各乡镇街道的空间回归结果与人口统计数据之间的相对误差进行精度评价，计算公式为 期李金香等：基于建筑物与 数据的人口空间化研究 （）式中，为第 个乡镇街道模拟人口的相对误差，为模拟人口数据。计算每个乡镇街道模拟人口的相对误差，结果如表 和图 所示。表 模拟人口相对误差序号乡镇街道相对误差 序号乡镇街道相对误差 序号乡镇街道相对误差 序号乡镇街道相对误差 巴楚镇 阿拉格尔乡 齐满镇 牙哈镇 色力布亚镇 琼库尔恰克乡 墩阔坦镇 乌尊镇 阿瓦提镇 英吾斯塘乡 雅克拉镇 伊西哈拉镇 三岔口镇 热斯坦街道 玉奇吾斯塘乡 水西沟镇 恰尔巴格乡 萨克萨克街道 比西巴格乡 永丰镇 多来提巴格乡 新城街道 哈尼喀塔木乡 板房沟镇 阿纳库勒乡 东城街道 阿克吾斯塘乡 萨尔达坂乡 夏玛勒乡 乌恰镇 阿格乡 甘沟乡 阿克萨克马热勒乡 阿拉哈格镇 塔里木乡 托里乡 图 模拟人口相对误差对比通过计算，比较各乡镇街道的人口模拟值与实际值的偏差情况。由表 和图 可知，经回归分析得到的 个行政单元人口相对误差值中，相对误差大于 的乡镇为 ，大于 的有 个乡镇，其中最大误差为 ；以上的乡镇其模拟人口相对误差范围在 以内。对人口模拟值和实际值结果进行拟合，得到两者之间的线性拟合率为 ，因此本文的模型具有较高的模拟精度。基于规则格网的可视化表达人口空间化将统计数据分布到规则的地理空间格网上，避免了原始统计单元面积大小不一造成的尺度混杂问题。通过小尺度规则格网的可视化表达，可以展示研究区更加精确的人口空间分布情况，结果如图 所示。本研究成果能较好地反映人口的空间分布特征，对于人口分布细节特征的刻画较为理想。由图 可知，受地形地貌、交通便利度、经济等因素的影响，研究区人口空间分布具有明显的差异性，表现为人口密度由城镇向乡村逐渐递减，城市中心人口密集，沿路网向外辐射，中国地震 卷城镇区域点状聚集，偏远乡村人口较少，沙漠、戈壁、山地等区域人烟稀少，人口密度差异巨大。基于建筑物与 数据的千米格网人口分布符合实际情况，可为灾情研判提供可靠的基础数据。图 研究区 格网人口分布 结论本文基于建筑物空间数据、路网数据、数据及人口统计数据，开展人口空间化方法研究，基于 相关分析选取相关性显著的模型构建因子，采用多元线性回归分析方法构建研究区人口空间化模型，实现了人口统计数据基于规则格网的更为精确的可视化表达。结果表明：（）受区域经济发展影响，不同区域模型构建的影响因子存在差异，通过对影响因子差异性大的区域进行分区建模，有助于提高模型构建的精度。（）相关分析方法可以快速提取模型构建的影响因子，并保证影响因子与因变量之间的相关性具有显著性；第一类分区中提取的模型构建影响因子为建筑物建筑面积、路网密度、购物，其 相关性系数分别为 、，第二类分区中提取的模型构建影响因子为建筑物建筑面积、居民区道路长度、政府机关，其 相关性系数为 、，建模因子均与人口存在显著正相关。（）多元线性回归分析方法可以对多个影响因子进行回归建模，将统计人口分布在建模因子分布的地区，建立的模型具有显著性，精度较高，计算的各乡镇街道的人口模拟值与实际值的偏差较小。（）通过小尺度规则格网的可视化表达可以展示研究区更加精确的人口空间分布情况。研究区人口分布具有明显的空间差异，由城镇向乡村区域递减的趋势明显，中心城区、周边城镇、偏远乡镇之间的人口密度差异巨大。总的来看，基于建筑物与 数据的人口空间化方法适用性强，精度较高，人口空间化成果能较好地反映实际人口的空间分布特征，且细节刻画较为准确，可以为地震应急救援提供决策依据，有利于在救灾初期帮助决策者对灾情做出正确判断和评估，提高地震应急救援的时效性。期李金香等：基于建筑物与 数据的人口空间化研究参考文献柏中强，王卷乐，姜浩，等 基于多源信息的人口分布格网化方法研究 地球信息科学学报，（）：淳锦，张新长，黄健锋，等 基于 数据的人口分布格网化方法研究 地理与地理信息科学，（）：，崔晓临，张佳蓓，吴锋，等 基于多源数据融合的北京市人口时空动态分析 地球信息科学学报，（）：董春，刘纪平，赵荣，等 地理因子与空间人口分布的相关性研究 遥感信息，（）：董南，杨小唤，蔡红艳 人口数据空间化研究进展 地球信息科学学报，（）：杜培培，侯西勇 基于多源数据的中国海岸带地区人口空间化模拟 地球信息科学学报，（）：高倩，阿里木江卡斯木 基于 夜间灯光数据的天山北坡城市群人口分布空间模拟 西北人口，（）：高义，王 辉，王 培涛，等 基于人口普查 与 多源夜 间灯 光 数 据 的 海岸 带 人 口 空 间 化分 析资 源 科 学，（）：黄河清，王有亮，胡宝清，等 基于神经网络和 的广西都安县人口数据空间化研究 测绘与空间地理信息，（）：江东，杨小唤，王乃斌，等 基于 、的人口空间分布研究 地球科学进展，（）：刘焕金 基于多源数据的太湖流域人口空间化研究 硕士学位论文 南京：南京农业大学刘正廉，桂志鹏，吴华意，等 融合建筑物与 数据的精细人口空间化研究 测绘地理信息，（）：吕安民，刘海启，李成名，等 人口密度的面积内插算法研究 中国农业资源与区划，（）：孙艳萍，陈文凯，张苏平，等 基于居民点密度的人口空间化及在四川九寨沟 级地震中的应用 地震工程学报，（增刊）：唐奇，许新宜，俞淞，等 基于 的人口空间离散化方法及其应用 以北方地区为例 北京师范大学学报（自然科学版），（）：王春菊，汤小华 基于 的福建省人口统计数据空间化 地理与地理信息科学，（）：王芳，朱明稳，陈崇旺，等 基于步行指数与人口空间格局的城市健康资源公平性 以广州市中心城区为例 资源科学，（）：王珂靖，蔡红艳，杨小唤 多元统计回归及地理加权回归方法在多尺度人口空间化研究中的应用 地理科学进展，（）：王明明，王卷乐 基于夜间灯光与土地利用数据的山东省乡镇级人口数据空间化 地球信息科学学报，（）：王晓洁，王卷乐，薛润生 基于普查和手机定位数据的乡镇尺度人口空间化方法研究 地球信息科学学报，（）：王雪梅，李新，马明国 基于遥感和 的人口数据空间化研究进展及案例分析 遥感技术与应用，（）：谢江丽，李帅，姚远 地震应急数据库中人口数据预测 以乌鲁木齐市为例 中国地震，（）：徐敬海，安基文，聂高众 基于千米格网的地震应急灾情预评估数据开发 地震地质，（）：于婷婷 基于珞珈一号夜间灯光数据的大连市甘井子区人口空间化研究 科学技术创新，（）：赵鑫，宋英强，刘轶伦，等 基于卫星遥感和 数据的人口空间化研究 以广州市为例 热带地理，（）：周中红，陈文凯，何少林，等 基于人口千米格网的地震压埋人员分布评估方法的应用与评价 以甘肃岷县漳县 级地震为例 地震研究，（）：朱鹏宇，程佳，杜龙 基于深度学习神经网络方法的中国大陆地震伤亡预测模型 中国地震，（）：，（）：，（）：中国地震 卷 ，（），：，：；