Characteristics and Formation Mechanism of Intra-Urban Employment Flows Based on Mobile Phone Data-Taking Wuhan City as an Example
-
摘要: 城市内部就业人口流动作为城市群体的主要移动形式，分析其特征及形成机理对城市规划、交通预测等具有重要意义。基于武汉市手机信令数据，识别职住人口分布与流动，构建城市内部就业流动网络。运用网络分析、可达性计算、逻辑回归等方法，分析城市内部就业流动的特征及其形成机制。研究表明，武汉市内部就业流动在数量上分布不均衡，大量就业流动集中于少数街道间。在空间上，就业流动随距离、可达时间增加而减少，并依地形、文化形成若干联系紧密的就业社区；以就业流出地居住人口、流入地工作人口度量的就业势能是驱动就业流动的最主要因素，而文化差异、空间不邻近、可达性差阻碍就业流动的发生。此外，不同产业特色对就业流动影响不同，商业、科教阻碍就业外流，工业吸引外来就业。Abstract: Human mobility is a cross-disciplinary research hotspot which reflects the complex man-land relationship. Understandings of the intra-urban employment population flow, as an important part of urban group mobility, are crucial for urban planning and traffic forecasting. The common use of location-awareness devices such as mobile phones enables to capture human behavioral data for analyzing intra-urban employment flow. In this paper, we attempt to characterize the employment flow and interpret its formation mechanism using mobile phone data recorded during 30 days (from June 1, 2016 to June 30, 2016) in Wuhan. Considering the spatial distribution of cellular base stations, we divide our study area into grids sized 250 m×250 m. We then adopt approaches to extract daily individual mobility information by identifying positon and duration of stagnations. Next, we propose several rules to infer in which grids people work and reside from individual mobility information we extracted and estimate the job/housing distribution as well as the employment flows between sub-districts. We then build the employment flows network, and methods such as network centrality analysis, network density analysis and community detection are applied to discover the pattern and characteristic of employment flows. Finally, we examine the influence of employment potential, traffic accessibility, spatial proximity, cultural difference and major industry on the size of employment flow between sub-districts and interpret its formation mechanism using logistic regression. The results firstly show that the number of intra-urban employment flow in Wuhan is unevenly distributed as large amounts of employment flow are concentrated in a few sub-districts. Secondly, the employment flows mainly occur in the sub-districts with good spatial proximity. It shows the distribution of the outflow employments decrease gradually as distance and accessibility time increase, and forms several employment communities (such as Hanyang, Wuchang North, Wuchang South, et al.) in the influence of geographical barriers and cultural differences. Thirdly, the employment potential measured by the resident population of the origin sub-districts and the working population of the destination sub-districts is the most important factor to increase the size of employment flow. Rather, cultural differences and poor accessibility hinder employment flow. The result also verifies that differences in major industries affect the employment flow, that commerce and education have a negative effect on the outflow of employments, but the industry attracts employments.
-
高精度数字高程模型(digital elevation model，DEM)是国家重要的战略信息资源，是国家经济建设、社会发展和国防安全必不可少的基础图件。然而，在植被覆盖区，传统航测及光学遥感手段只能获取植被层表面高程信息，无法测绘“裸地球”真实形状。真实地形往往需要通过人工实地调绘估计平均植被高，从光学手段获得的地形数据中扣除，这样的测图手段在植被覆盖区精度差、效率低，难以真实表达林下地形信息。因此，如何准确获得森林高度并将其从遥感或航测手段获得的高程中扣除，或直接得到不包含植被高的真实地形高程，对获取大范围高精度DEM极为重要。
合成孔径雷达干涉测量(interferometric synthetic aperture radar，InSAR)因其大范围、高分辨率、高精度、全天候监测的特点已被列为局部及全球尺度地形测绘的有效技术途径之一[1-3]。此外，与传统光学遥感只能进行“表面监测”不同，InSAR采用微波频段工作模式，在森林覆盖区具备穿透植被层并记录植被垂直结构、物理属性信息的能力，为森林高度反演及林下地形测绘提供了可能[4]。但是，传统InSAR技术在森林区测量的地形高度位于森林冠层顶部与地表之间，获取的DEM产品中包含森林信号，仍未实现真实林下地形测绘。其根本原因在于传统InSAR技术受观测信息不足的限制, 无法区分同一分辨单元内不同散射体及其垂直分布，无法分离森林覆盖区冠层及地表回波信号[4]。
为了解决上述难题，研究人员提出通过增加先验信息或观测量信息使得传统InSAR技术具备区分植被冠层及地表回波信号的能力。基于以上策略，发展了多种林下地形反演方法，主要可以归为3类。①基于InSAR技术的反演方法，主要包括基于相干散射模型的反演方法[5-6]、森林高扣除法[7-8]、时-频分析方法(time-frequency analysis, TFA)[9-10]；②基于极化InSAR(polarimetric InSAR, PolInSAR)技术的林下地形反演方法，主要包括基于相干最优理论的相位分离法[4, 11-13]、基于相干散射模型的反演方法[14-15]、基于相干散射矩阵的散射机理分解方法[16-17]SAR；③基于多基线InSAR/PolInSAR数据的层析SAR(tomographic SAR, TomoSAR)技术林下地形反演方法[18]。
本文旨在通过分析上述基于InSAR技术的林下地形反演方法及其研究现状，使读者能够对林下地形测绘进行系统性了解。
1 面向林下地形测绘的InSAR测高技术与现状
1.1 InSAR测高原理及发展现状
InSAR技术集SAR技术与干涉测量技术于一体。SAR是一种主动式微波遥感，单景SAR影像记录了SAR回波信号的强度信息及相位信息。微波信号为极化波，其散射过程具有对散射目标的几何信息(如形状、大小、空间排布规律)及物理属性信息(如介电常数)敏感的特点，这些信息对散射过程的影响综合表现为SAR影像的强度及相位信息[15]。当相同SAR传感器对同一目标以不同入射角进行观测时，对获取的两景SAR影像进行干涉处理即可用于获取地表高程信息，如图 1(a)所示。图 1中H为SAR平台高度，R为SAR传感器到地表目标点P之间的距离，称为斜距；h为目标点P的高程；B为基线；θ为入射角；α为基线倾角；S表示获取的SAR影像。通过引入干涉测量技术，微波信号具备测高能力，但是受观测信息不足的限制，InSAR只能测量同一分辨单元内多个散射体的综合高度，无法捕捉不同散射体垂直向分布规律。
为了解决该问题，需引入外部先验信息加以约束或增加观测信息。鉴于此，学者们提出在InSAR测高技术基础之上引入极化测量技术，即PolInSAR[4]。从数据获取角度来看(图 1(b))，PolInSAR将传统单极化InSAR工作模式(HH或VV)扩展为多极化InSAR工作模式(HH，HV，VV，VH)，增加了观测信息。此外，不同极化方式的微波信号对散射目标的几何信息及物理属性信息具有不同的敏感程度，具备区分同一分辨单元内不同散射体的能力。以上特点使得PolInSAR为区分同一分辨单元内不同散射体及其垂直分布提供了可能。
除了通过变换极化方式使传统InSAR技术具有分辨同一分辨单元内不同散射体的能力之外，TomoSAR技术通过增加基线数量扩充观测空间(图 1(c))，结合InSAR测高原理与谱估计技术，实现了LOS(line of sight)法向合成孔径，将传统二维成像扩展为三维成像，具备测量同一分辨单元内不同散射体的能力[18]。
1.2 面向林下地形测绘的SAR系统发展现状
文献[19]较为全面系统地总结了星载SAR系统的发展历史，可知总体上现有星载SAR系统特点为多波段、多极化、高时空分辨率。但是，森林几何形态特征及物理属性具有高时变性，要在森林区开展InSAR测量要求星载系统重返周期需足够短，以抵抗时间去相干对干涉质量的影响。对于地形测绘而言，最理想的数据获取模式为“零”时间基线。为此，2000年，美国实施SRTM(Shuttle Radar Topography Mission)计划，利用航天飞机搭载双天线SAR传感器，获取了覆盖全球80%面积的30 m及90 m分辨率DEM[2]。之后，德国分别于2007年、2010年先后发射了TerraSAR-X及TanDEM-X姊妹星组成分布式协同工作模式，获取了全球12 m分辨率的DEM[3]。尽管上述两种DEM产品仍然包含植被高信号，制约了某些应用，但是作为全球地形基础数据，在众多领域发挥了无法替代的作用。
相比于星载系统，机载系统由于研制成本低、调度灵活等优点已被广泛用于星载系统的仿真实验。表 1为现有主要机载系统及其参数，利用这些机载系统在森林覆盖区开展了广泛研究。通过对植被覆盖区参数反演的研究，研究人员意识到L、P波段的SAR信号由于强穿透能力，可采集更为丰富的植被垂直结构信息及林下地形信息，在森林参数反演方面比短波长更具优势。鉴于此，欧空局规划在2020年发射P波段BIOMASS卫星[20]，德国航天局预计在2022年发射L波段Tandem-L卫星[21]。这两个卫星均将林下地形测绘纳为主要监测任务之一。
表 1 现有机载SAR系统及参数Table 1. The Airborne SAR Systems and Their ParametersSAR传感器 国家 开始年限 波段 极化方式 AIRSAR 美国 1988 P/L/C 多极化 EMISAR 丹麦 1995 X/L 多极化 CCRS-SAR 加拿大 1995 C 全极化 E-SAR 德国 1996 X/C/L/P 多极化 PI-SAR 日本 1996 X/L 多极化 GEO-SAR 美国 1999 X/P 多极化 RAMSES-SAR 法国 2004 W-P8种波段 多极化 F-SAR 德国 2007 X/C/L/P 多极化 UAV-SAR 美国 2007 L 全极化 PI-SAR2 日本 2008 X 全极化 CASM-SAR 中国 2009 X/P 多极化 CET38-SAR 中国 2009 X/P 多极化 SETHI-SAR 法国 2009 X/C/L/P 多极化 2 InSAR林下地形测绘方法与进展
2.1 基于InSAR技术的林下地形测绘方法
1) 森林高扣除法。该方法核心思想为借助外部手段估计森林高度，并将其从InSAR反演的DEM中扣除[7-8]。该方法存在的主要问题为：①如何快速、大范围、高精度获取森林高度；②由于InSAR DEM采用微波信号获取，在植被覆盖区具有一定的穿透深度，该穿透深度在空间上受森林密度、含水量及地形等方面的影响呈现显著差异，增加了树高扣除法的不确定性。
2) TFA方法。TFA方法主要原理为通过子孔径分解技术，获取具有不同方位向入射角的子视图像[9-10]。微波信号对散射目标的大小、形状及空间排布规律等几何信息敏感。由于森林目标的密度分布呈现各向异性，当SAR传感器从不同角度观测同一森林场景时，由于密度的差异导致不同视角条件下的后向散射信息存在差异，从而导致不同子视图像具有不同的相位中心。甄选含有最强地表散射贡献的子视回波信号，并认为其相位中心位于地表，进而实现林下地形估计。该方法无需借助任何外部数据、先验信息及散射模型即可实现林下地形反演。但是该方法精度受限于微波信号在森林区的穿透深度，只适用于P波段数据且稀疏森林覆盖条件下的林下地形反演。
3) 基于相干散射模型反演的方法。该方法的主要思路为建立InSAR观测量、微波散射过程及森林生物物理参数之间的关联。为此，将森林覆盖场景抽象概括为植被层及地表层两层，其中植被层由各向同性随机分布的散射粒子构成，该模型被称为随机地体二层散射模型(random vo- lume over ground, RVoG)[5-6]。但是，由于InSAR数据无法提供足够的观测量用于该模型求解，需通过其他手段获取部分模型参数作为已知量。实际应用中，高分辨率、高精度先验信息的难以获取极大限制了该方法的适用范围。
2.2 基于PolInSAR技术的林下地形测绘方法
1) 基于相干最优理论的相位分离法。不同极化方式的微波信号对散射目标的几何信息及物理属性信息具有不同的敏感程度，如HV极化微波信号对植被冠层散射敏感，而HH极化微波信号对地表散射回波信号敏感。根本原因在于不同极化方式的微波信号对应不同的散射机理。基于此，在极化域内，可以甄选对地表散射(表面散射、二面角散射等)最为敏感的极化方式，认为其相位中心位于地表，实现林下地形估计[4, 11-13]。该方法没有考虑森林生物物理参数与PolInSAR观测量之间的关联，可靠性受限于微波信号的穿透深度。当微波具有较小穿透深度时，该方法无法探测到真实地表散射信息。
2) 基于极化相干散射模型反演方法。在单极化InSAR模式下，由于观测信息不足，无法对RVoG模型进行直接解算。在PolInSAR测量模式建立之后，文献[14]首次将RVoG模型引入PolInSAR测量，建立了森林生物物理参数与PolInSAR观测量之间的关联，无需借助任何外部数据，实现了林下地形、森林高度等参数的反演。至此，基于PolInSAR数据的森林参数反演得到国内外研究者的广泛关注与研究。
然而，随着研究的深入，研究人员意识到该反演方法的关键在于RVoG模型能否准确表达微波信号与森林场景的相互作用过程。为此，提出一系列改进模型。如RVoG+VTD(volume temporal decorrelation)用于考虑介电常数变化、风动引起时间去相干对植被冠层干涉信号的影响[22-23]；RMoG(random motion over ground)模型用于考虑风动对植被冠层、地表层干涉信号的影响[24-25]；S-RVoG(slope-RVoG)用于考虑地形坡度对微波信号在植被层穿透深度的影响[26-28]；三层RVoG模型[17]、高斯后向散射模型[29-30]及可变消光系数模型[29, 31]用于考虑植被垂直结构不均匀造成的影响；双站干涉模式下的RVoG模型用于考虑表面散射与二面角散射相位中心不一致的问题[32]。通过上述改进模型，显著提升了经典RVoG模型应对复杂森林场景的能力。
在解算方法方面，文献[14]将RVoG模型参数解算概括为六维非线性方程求解问题；之后，文献[33]针对非线性求解中初值难以确定，容易陷入局部最优等问题提出了三阶段反演方法，该方法直观、简单、稳健，提出后得到广泛关注。尽管三阶段算法大大降低了RVoG模型的解算难度，但是该方法只适用于单基线PolInSAR数据，无法扩展到多基线数据进行联合解算，即便存在多基线数据时，也只是对每条基线数据反演结果进行简单加权平均[34-37]。针对该问题，文献[38-39]分别提出了多基线非线性迭代算法及复数平差法。这些方法不仅通过不同几何条件下观测数据的约束关系提升了参数求解的稳健性，而且能够充分利用多余观测信息考虑更为复杂、更能表征森林散射特征的散射模型。
综上所述，基于极化相干散射模型的反演方法已成为研究人员研究林下地形及森林高度提取的关注焦点。但仍存在一些问题需要考虑。首先，同频率SAR信号视野下，散射体几何属性存在显著差异，如X波段视野下散射体主要为小枝叶，而P波段视野下散射体为大型枝干。这种条件下，散射体空间排布规律存在显著差异，RVoG模型与森林场景散射过程匹配程度发生变化，导致反演结果存在较大偏差、甚至失败。其次，森林几何形态特征、物理属性具有显著的时变性，现有时间去相干模型尚不能准确刻画外界因子对微波散射过程的影响，时间去相干仍是森林参数反演的主要误差源之一。
3) 基于极化相干散射矩阵的散射机理分解方法。该方法思想来源于Freeman三分量分解，假定微波信号与森林场景作用过程中主要发生表面散射(地表)、二面角散射(树干-地表)及体散射(冠层)[16]3种散射机制。在此基础之上，建立3种散射机理的数学表达模型，求取不同散射机理的相位中心。之后，文献[17]以类似思想构建了基于Neumann极化分解模型的反演框架，与前一种分解方法相比，该方法通过有效结合极化散射分解与相干散射模型，建立了极化相干矩阵与生物物理参数之间的关联。该方法提出后，由于模型解算复杂程度高，并没有得到广泛应用。但是，该反演框架为全面考虑森林场景异质性，建立更为精准的散射模型提供了思路，该问题还有待深入研究。
2.3 基于TomoSAR技术的林下地形测绘方法
文献[18]首次用机载L波段数据对森林区进行了三维成像。但是由于采用传统的FFT(fast Fourier transformation)成像方法，垂直向分辨率较低。因此，后续大量研究工作主要关注层析谱估计方法，一系列超分辨率算法相继提出，主要归为非参数化频谱估计法(包括Beamfor- ming、Capon、SVD(sigular value decomposition)等)、参数化频谱估计方法(包括MUSIC(multiple signal classification algorithm)、WSF(weighted subspace fitting)、NLS(non-linear least squares)等)及压缩感知成像方法[40-43] 3类。这些研究工作重点关注如何提升高度向分辨率，尚未对TomoSAR林下地形反演性能进行综合评价。此外，这些方法通常假定散射体为点目标，但在森林区，散射目标为分布式目标，导致这些方法估计的层析谱不够理想。为了解决该问题，文献[44]提出基于COMET(covariance matching estimation technique)的TomoSAR方法，该方法可适用于多基线InSAR/PolInSAR数据，对散射体空间几何特性及散射噪声水平均有良好的容忍度。之后，研究人员发现，单纯估计不同高度位置的后向散射强度信息不足以支撑深入分析不同高度位置的散射机理与回波信号之间的关联。针对该问题，文献[45-46]提出采用SKP(sum of Kronecker product)分解技术对多基线PolInSAR数据的极化矩阵及结构矩阵进行分离，并通过代数合成方法提取体散射及地面散射信息。通过结合极化信息，使得传统TomoSAR不仅能反映不同高度的后向散射强度信息，而且能提供散射机理信息，为精细化分析植被垂直结构与微波信号作用过程提供重要技术支撑。
3 InSAR林下地形测绘应用与进展
3.1 基于InSAR技术的林下地形测绘
文献[5]首次利用C波段机载InSAR数据结合相干散射模型实现了林下地形提取，反演精度为6.5 m。文献[8]利用ICESat (ice, cloud and land elevation satellite)/MODIS(moderate-resolution imaging spectroradiomete)联合估计的全球植被高(空间分辨率1 km)剔除C波段STRM DEM(空间分辨率90 m)中包含的森林高信号，修正后的森林区SRTM DEM精度由14.0 m提高到6.0 m。文献[10]利用P波段机载InSAR数据采用TFA法提取了稀疏针叶林覆盖区林下地形，林下地形精度为2.0 m。
3.2 基于PolInSAR技术的林下地形测绘
尽管PolInSAR技术提出后已广泛用于森林高度提取，但是林下地表干涉相位只作为中间过程参数并没有深入展开应用研究。文献[4]利用SIR-C/X-SAR L波段PolInSAR数据，采用相干最优相位分离法进行森林区DEM反演，与传统InSAR DEM相比，精度由12.9 m提高到7.8 m。文献[47]利用L波段PolInSAR数据结合相干散射模型进行针叶林区林下地形反演，反演结果与LiDAR结果相比偏差均值为2~3 m，并发现夏季与冬季反演的林下地形精度有差异，冬季结果精度更高。文献[35]利用L、P波段机载多基线PolInSAR针叶林区数据结合相干散射模型进行林下地形反演，结果表明，与LiDAR DEM相比，L波段误差主要分布在-10~10 m、P波段误差主要分布在-7~7 m，由于P波段具有强穿透性，更利于林下地形反演。文献[39]利用L、P波段机载多基线PolInSAR复数平差法对针叶林区进行林下地形反演，结果表明，与LiDAR DEM相比，L波段DEM误差为3.5 m, P波段DEM误差为3.2 m。基于极化相干散射矩阵的散射机理分解方法虽具备提取林下地形的能力，尚未见相关文献对其在林下地形反演做定量评价。
3.3 基于TomoSAR技术的林下地形测绘
文献[44]采用9景机载P波段数据结合COMET成像方法对针叶林覆盖区(平地)进行林下地形提取，与LiDAR DEM相比，反演精度为0.9 m。文献[43]采用11景P波段、16景L波段机载数据结合TDBP(time-domain back-projection)成像方法对针叶林覆盖区(平地)进行林下地形提取，与LiDAR相比，L波段DEM反演精度为1.37 m, P波段DEM精度为0.83 m。文献[48]采用6景机载L波段数据结合M-RELAX成像方法对针叶林覆盖区(丘陵)进行林下地形提取，与LiDAR相比，反演的DEM精度为2.0 m。文献[40]采用6景机载P波段数据结合SKP成像方法对热带雨林区(丘陵)进行林下地形提取，与LiDAR相比，反演精度为2.41 m。
3.4 综合分析
总体上，目前林下地形反演还处于方法测试阶段。实验主要采用在相对平坦地形条件下获取的机载数据。从已有研究成果可以看出，3种技术反演的林下地形精度与地形条件、森林属性及数据属性有关。此外，在林下地形反演精度方面，TomoSAR技术最优，InSAR技术最差，PolInSAR技术折中。但是，TomoSAR技术要求在短时间内获取同一地区的大量影像，且对影像基线空间分布要求严格，致使该技术难以应对大范围的林下地形测绘，但可用于局部区域精细化监测。基于树高扣除法的InSAR技术反演林下地形方法只需要少量SAR数据，适用于大范围林下地形测绘，但是该方法受制于树高产品的精度。从现有遥感技术来看，微波遥感是获取大范围、高精度树高产品的最优候选技术之一，但是树高与林下地形对于微波遥感来说是同源问题。基于相干散射模型的InSAR技术反演林下地形受观测信息不足的限制，已被基于相干散射模型的PolInSAR技术取代。尽管同等分辨率条件下，SAR传感器需要牺牲一半的覆盖幅宽来采集全极化数据，但是PolInSAR无需先验信息、无需大量数据的特点已被德国宇航局Tandem-L计划选为绘制全球林下地形的基本方法。然而，基于PolInSAR技术的林下地形反演精度还有很大提升空间。
4 InSAR林下地形测绘面临的挑战
4.1 传感器平台时空基线设计
目前，在轨运行卫星除了TerraSAR-X/Tandem-X星座在全球局部区域开展了森林参数反演飞行计划外，其他卫星受时空基线、波段的限制，尚不具备开展林下地形测绘的功能[49]。尽管规划卫星Tandem-L及BIOMASS卫星在抗时间去相干、强穿透能力方面填补了现有卫星无法应对林下地形测绘的缺陷，但是，时空基线的优化仍是获取高精度林下地形产品的关键技术之一。森林形态特征、生物物理属性在时间上呈现显著动态性，倘若卫星平台时间基线不能很好考虑森林的动态性，其引起的干涉失相关误差将成为林下地形反演的重要误差源。此外，文献[50]系统讨论了空间基线对森林高度反演的影响及优化方案，其优化原则为权衡干涉信号对冠层体散射的敏感度及高相干性，干涉基线可变范围有限。但是，对于林下地形反演而言，长基线干涉系统对地形起伏具有更高敏感度，利于高精度地形测绘。因此，如何优化面向林下地形测绘的基线几何值得探索。
4.2 轨道误差与大气延迟改正
轨道误差、大气延迟是制约高精度地形测绘的重要因素。尽管Tandem-L采用双站模式获取的干涉对不受大气延迟的影响，但是，当需要选择多基线模式增强测高能力时，还是要处理不同时刻获取的SAR影像，此时，仍需考虑大气延迟误差。文献[19]已详细论述了解决这两类误差策略及所面临的挑战。此外，地形测绘中处理对象往往为小数据集，传统相位定标技术及时空滤波分析难以适用[51-52]。鉴于此，如何利用不同极化方式干涉图之间的自洽关系实现轨道误差与大气延迟改正值得深入研究。
4.3 散射机理解译及建模
目前，森林生物物理参数与微波信号散射机理之间的关联尚未清晰。散射机理解译的精细化程度及准确建模是基于相干散射模型反演策略的基础，同时也是建立层析谱与生物物理参数之间关联的重要理论支撑。因此，如何揭示生物物理参数与散射机理之间的关系，是保证复杂森林覆盖及地形条件下准确获取林下地形的关键。
5 结语
传统遥感技术获取的包含树高信号的DEM已不能满足众多研究领域对林下地形的迫切需求。鉴于此，研究人员根据微波信号在植被区具备穿透监测的特点并结合干涉测高原理，发展了基于InSAR、PolInSAR及TomoSAR技术的林下地形测高理论与方法。但受外部先验信息不可靠、散射机理理论体系不完备、传感器平台参数不合理、误差控制体系不完善等因素限制，致使林下地形测绘还处于理论与方法研究探索阶段。随着传感器平台的不断升级及相应数据处理理论的不断提升，林下地形提取未来将成为地形测绘领域的重点研究方向之一。
-
![]()
图 3 研究区主要用地类型及商圈分布图与街道产业类型分布图
Figure 3. Land Use and Industrial Categories of Sub-districts
![]()
图 5 不同街道特色产业类型就业流入地、流出地流量分布箱线图s
Figure 5. Boxplot of Outflow and Inflow Logarithmic Population in Different Categories Sub-districts
![]()
图 6 不相邻街道与相邻街道的跨街道人口流动对数频率分布直方图
Figure 6. Frequency Distribution Histograms of Logarithmic Flow Population Between Adjacent and Nonadjacent Sub-districts
![]()
图 7 研究区就业人口流动与可达性散点图
Figure 7. Scatter Diagrams of Accessibility and Logarithmic Flow Population Between Sub-districts
![]()
图 8 典型街道就业人口分布图
Figure 8. Employment Distributions of Some Typical Sub-districts in the Study Area
![]()
图 9 就业流动网络社区发现结果空间分布图
Figure 9. Community Structures Based on the Employment Linkages in the Study Area
![]()
图 10 城市内部区域间就业人口流动形成机制图
Figure 10. Formation Mechanism of Intra-urban Employment Flows
表 1 研究区3大文化区就业网络密度
Table 1 Network Density Analysis of Employment Flows Network in Hankou, Hanyang and Wuchang
文化区 汉口(流出) 汉阳(流出) 武昌(流出) 汉口(流入) 0.288 0.050 0.030 汉阳(流入) 0.036 0.404 0.005 武昌(流入) 0.049 0.038 0.348 
下载: 导出CSV
表 2 逻辑回归结果表
Table 2 Logistic Regression Results
变量 系数估计 标准化系数 z-值 P-值 显著性 常数项 -1.169 0 0 -7.577 0 3.54×10-14 *** 流出地居住人口 5.50×10-5 0.961 4 25.916 0 ＜2×10-16 *** 流入地工作机会 3.48 ×10-5 1.024 0 26.245 0 ＜2×10-16 *** 两街道是否直接相邻 2.257 0 0.293 3 13.177 0 ＜2×10-16 *** 两区域公交可达性 -0.046 5 -0.672 5 -8.343 0 ＜2×10-16 *** 两区域驾车可达性 -0.124 1 -0.735 0 -8.782 0 ＜2×10-16 *** 两街道是否跨越文化区 -2.233 0 -0.595 1 -16.917 0 ＜2×10-16 *** 流出街道是否具有工业特色 -0.305 5 -0.052 7 -1.594 0 0.111 0 流出街道是否具有商业特色 -1.135 0 -0.177 2 -5.873 0 4.28×10-9 *** 流出街道是否具有科教特色 -1.785 0 -0.307 7 -11.237 0 ＜2×10-16 *** 流入街道是否具有工业特色 0.661 4 0.114 0 3.521 0 0.000 4 *** 流入街道是否具有商业特色 0.095 5 0.014 9 0.565 0 0.571 8 流入街道是否具有科教特色 0.213 1 0.036 8 1.539 0 0.123 7
下载: 导出CSV
-
[1] 陆锋, 刘康, 陈洁.大数据时代的人类移动性研究[J].地球信息科学学报, 2014, 16(5):665-672 http://d.old.wanfangdata.com.cn/Periodical/dqxxkx201405001 Lu Feng, Liu Kang, Chen Jie. Research on Human Mobility in Big Data Era[J]. Journal of Geo-Information Science, 2014, 16(5):665-672 http://d.old.wanfangdata.com.cn/Periodical/dqxxkx201405001
[2] González M C, Hidalgo C A. Understanding Individual Human Mobility Patterns[J]. Nature, 2008, 453(7196):779-782 doi: 10.1038/nature06958
[3] Kaltenbrunner A, Meza R, Grivolla J, et al. Urban Cycles and Mobility Patterns:Exploring and Predicting Trends in a Bicycle-Based Public Transport System[J]. Pervasive & Mobile Computing, 2010, 6(4):455-466 http://dl.acm.org/citation.cfm?id=1860501
[4] Long Y, Han H, Tu Y, et al. Evaluating the Effectiveness of Urban Growth Boundaries Using Human Mobility and Activity Records[J]. Cities, 2015, 46:76-84 doi: 10.1016/j.cities.2015.05.001
[5] Xie K, Xiong H, Li C. The Correlation Between Human Mobility and Socio-demographic in Megacity[C]. IEEE Smart Cities Conference, Trento, Italy, 2016 https://ieeexplore.ieee.org/document/7580733
[6] 周素红, 杨利军.广州城市居民通勤空间特征研究[J].城市交通, 2005, 3(1):62-67 doi: 10.3969/j.issn.1672-5328.2005.01.015 Zhou Suhong, Yang Lijun. Study on the Spatial Characteristic of Commuting in Guangzhou[J].Urban Transport of China, 2005, 3(1):62-67 doi: 10.3969/j.issn.1672-5328.2005.01.015
[7] 孟斌.北京城市居民职住分离的空间组织特征[J].地理学报, 2009, 64(12):1457-1466 http://d.old.wanfangdata.com.cn/Periodical/dlxb200912006 Meng Bin. The Spatial Organization of the Separation Between Jobs and Residential Locations in Beijing[J].Acta Geographica Sinica, 2009, 64(12):1457-1466 http://d.old.wanfangdata.com.cn/Periodical/dlxb200912006
[8] 刘志林, 张艳, 柴彦威.中国大城市职住分离现象及其特征——以北京市为例[J].城市发展研究, 2009, 16(9):110-117 doi: 10.3969/j.issn.1006-3862.2009.09.020 Liu Zhilin, Zhang Yan, Chai Yanwei. Home-Work Separation in the Context of Institutional and Spatial Transformation in Urban China:Evidence from Beijing Household Survey Data[J].Urban Development Studies, 2009, 16(9):110-117 doi: 10.3969/j.issn.1006-3862.2009.09.020
[9] 刘志林, 王茂军.北京市职住空间错位对居民通勤行为的影响分析——基于就业可达性与通勤时间的讨论[J].地理学报, 2011, 66(4):457-467 http://d.old.wanfangdata.com.cn/Conference/8070703 Liu Zhilin, Wang Maojun. Job Accessibility and Its Impacts on Commuting Time of Urban Residents in Beijing:From a Spatial Mismatch Perspective[J].Acta Geographica Sinica, 2011, 66(4):457-467 http://d.old.wanfangdata.com.cn/Conference/8070703
[10] 周素红, 闫小培.基于居民通勤行为分析的城市空间解读——以广州市典型街区为案例[J].地理学报, 2006, 61(2):179-189 doi: 10.3321/j.issn:0375-5444.2006.02.007 Zhou Suhong, Yan Xiaopei. The Impact of Commuters' Travel Pattern on Urban Structure:A Case Study in Some Typical Communities in Guangzhou[J].Acta Geographica Sinica, 2006, 61(2):179-189 doi: 10.3321/j.issn:0375-5444.2006.02.007
[11] 徐卞融, 吴晓.基于"居住-就业"视角的南京市流动人口职住空间分离[J].城市规划学刊, 2010(5):87-97 doi: 10.3969/j.issn.1000-3363.2010.05.011 Xu Bianrong, Wu Xiao. Quantitative Analysis of the Separation of Residential and Employment of the Floating Population in Nanjing[J].Urban Planning Forum, 2010(5):87-97 doi: 10.3969/j.issn.1000-3363.2010.05.011
[12] Calabrese F, Lorenzo G D, Liu L, et al. Estimating Origin-Destination Flows Using Mobile Phone Location Data[J].IEEE Pervasive Computing, 2011, 10(4):36-44 doi: 10.1109/MPRV.2011.41
[13] Alexander L, Jiang S, Murga M, et al. Origin-Destination Trips by Purpose and Time of Day Inferred from Mobile Phone Data[J]. Transportation Research Part C Emerging Technologies, 2015, 58:240-250 doi: 10.1016/j.trc.2015.02.018
[14] 钮心毅, 丁亮.利用手机数据分析上海市域的职住空间关系——若干结论和讨论[J].上海城市规划, 2015(2):39-43 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=shcsgh201502010 Niu Xinyi, Ding Liang. Analyzing Job-Housing Spatial Relationship in Shanghai Using Mobile Phone Data:Some Conclusions and Discussions[J].Shanghai Urban Planning Review, 2015(2):39-43 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=shcsgh201502010
[15] 高志刚, 唐小勇, 赵必成, 等.基于手机信令的组团城市通勤特征分析——以重庆主城区为例[J].城乡规划:城市地理学术版, 2016(3):72-77 http://www.cqvip.com/QK/71524X/201603/669451808.html Gao Zhigang, Tang Xiaoyong, Zhao Bicheng, et al. Analysis on the Characteristics of Urban Commuter Flow Based on Mobile Phone Signal Data:A Case Study of Chongqing[J]. Urban and Rural Planning, 2016(3):72-77 http://www.cqvip.com/QK/71524X/201603/669451808.html
[16] Phithakkitnukoon S, Horanont T, Lorenzo G D, et al. Activity-Aware Map: Identifying Human Daily Activity Pattern Using Mobile Phone Data[C]. Human Behavior Understanding, 1st International Workshop, HBU 2010, Istanbul, Turkey, 2010 http://www.springerlink.com/content/jj21508881433584
[17] Jiang S, Ferreira J, Gonzalez M C. Activity-Based Human Mobility Patterns Inferred from Mobile Phone Data:A Case Study of Singapore[J].IEEE Transactions on Big Data, 2017, 3(2):208-219 doi: 10.1109/TBDATA.2016.2631141
[18] 赵志远, 尹凌, 方志祥, 等.轨迹数据的时间采样间隔对停留识别和出行网络构建的影响[J].武汉大学学报·信息科学版, 2018, 43(8):1152-1158 http://ch.whu.edu.cn/CN/abstract/abstract6167.shtml Zhao Zhiyuan, Yin Ling, Fang Zhixiang, et al. Impacts of Temporal Sampling Intervals on Stay Detection and Movement Network Construction in Trajectory Data[J].Geomatics and Information Science of Wuhan University, 2018, 43(8):1152-1158 http://ch.whu.edu.cn/CN/abstract/abstract6167.shtml
[19] Gao S, Liu Y, Wang Y, et al. Discovering Spatial Interaction Communities from Mobile Phone Data[J].Transactions in GIS, 2013, 17(3):463-481 http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0230166666/
[20] Louail T, Lenormand M, Cantu Ros O G, et al. From Mobile Phone Data to the Spatial Structure of Cities[J]. Sci Rep, 2014, 4(2973):5276 http://go.nature.com/mObYEV
[21] 杨喜平, 方志祥, 赵志远, 等.顾及手机基站分布的核密度估计城市人群时空停留分布[J].武汉大学学报·信息科学版, 2017, 42(1):49-55 http://ch.whu.edu.cn/CN/abstract/abstract5634.shtml Yang Xiping, Fang Zhixiang, Zhao Zhiyuan, et al. Analyzing Space-Time Variation of Urban Human Stay Using Kernel Density Estimation by Conside-ring Spatial Distribution of Mobile Phone Towers[J].Geomatics and Information Science of Wuhan University, 2017, 42(1):49-55 http://ch.whu.edu.cn/CN/abstract/abstract5634.shtml
[22] Ravenstein E G. The Laws of Migration[J]. Journal of the Statistic Society of London, 1885, 48(2):167-235 doi: 10.2307/2979181
[23] Mathä T, Wintr L. Commuting Flows Across Bordering Regions:A Note[J]. Applied Economics Letters, 2009, 16(7):735-738 doi: 10.1080/13504850701221857
[24] Fitzgerald J, Leblang D, Teets J C. Defying the Law of Gravity:The Political Economy of International Migration[J].World Politics, 2014, 66(3):406-445 doi: 10.1017/S0043887114000112
[25] Beine M, Bertoli S. A Practitioners' Guide to Gra-vity Models of International Migration[J]. World Economy, 2016, 39(4):496-512 doi: 10.1111/twec.2016.39.issue-4
[26] Ramos R, Suriñach J. A Gravity Model of Migration Between the ENC and the EU[J]. Tijdschrift Voor Economische EnSociale Geografie, 2017, 108(1):21-35 doi: 10.1111/tesg.12195
[27] Hatton T J, Williamson J G. What Fundamentals Drive World Migration?[J].Cepr Discussion Papers, 2002, 159(2):139-158 doi: 10.1057/9780230522534_2
[28] Hatton T J, Williamson J G. Demographic and Economic Pressure on Emigration out of Africa[J].Scandinavian Journal of Economics, 2003, 105(3):465-486 doi: 10.1111/sjoe.2003.105.issue-3
[29] 于文丽, 蒲英霞, 陈刚, 等.基于空间自相关的中国省际人口迁移模式与机制分析[J].地理与地理信息科学, 2012, 28(2):44-49 http://d.old.wanfangdata.com.cn/Periodical/dlxygtyj201202011 Yu Wenli, Pu Yingxia, Chen Gang, et al. Spatial Analysis of the Paterns and Mechanism of Inter-provincial Migration Flows in China[J].Geography and Geo-Information Science, 2012, 28(2):44-49 http://d.old.wanfangdata.com.cn/Periodical/dlxygtyj201202011
[30] 巫锡炜, 郭静, 段成荣.地区发展、经济机会、收入回报与省际人口流动[J].南方人口, 2013, 28(6):54-61 doi: 10.3969/j.issn.1004-1613.2013.06.007 Wu Xiwei, Guo Jing, Duan Chengrong. Regional Development, Economic Opportunities, Income Returns and Inter-provincial Migration[J]. South China Population, 2013, 28(6):54-61 doi: 10.3969/j.issn.1004-1613.2013.06.007
[31] 王珏, 陈雯, 袁丰.基于社会网络分析的长三角地区人口迁移及演化[J].地理研究, 2014, 33(2):385-400 http://d.old.wanfangdata.com.cn/Periodical/dlyj201402016 Wang Jue, Chen Wen, Yuan Feng. Human Mobility and Evolution Based on Social Network:An Empirical Analysis of Yangtze River Delta[J]. Geographical Research, 2014, 33(2):385-400 http://d.old.wanfangdata.com.cn/Periodical/dlyj201402016
[32] 陈伟, 刘卫东, 柯文前, 等.基于公路客流的中国城市网络结构与空间组织模式[J].地理学报, 2017, 72(2):224-241 http://d.old.wanfangdata.com.cn/Periodical/dlxb201702004 Chen Wei, Liu Weidong, Ke Wenqian, et al. The Spatial Structures and Organization Patterns of China's City Networks Based on the Highway Passenger Flows[J]. Acta Geographica Sinica, 2017, 72(2):224-241 http://d.old.wanfangdata.com.cn/Periodical/dlxb201702004
[33] Liu X, Derudder B, Wu K. Measuring Polycentric Urban Development in China:An Intercity Transportation Network Perspective[J].Regional Stu-dies, 2016(8):1302-1315 doi: 10.1080/00343404.2015.1004535
[34] 刘军.社会网络分析导论[M].北京:社会科学文献出版社, 2004 Liu Jun. An Introduction to Social Network Analysis[M].Beijing:Social Sciences Academic Press, 2004
[35] Wang J, Mo H, Wang F, et al. Exploring the Network Structure and Nodal Centrality of China's Air Transport Network:A Complex Network Approach[J]. Journal of Transport Geography, 2011, 19(4):712-721 doi: 10.1016/j.jtrangeo.2010.08.012
[36] Bagler G. Analysis of the Airport Network of India as a Complex Weighted Network[J]. Physica A Statistical Mechanics and Its Applications, 2004, 387(12):2972-2980 http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_cond-mat%2f0409773
[37] Ran B. Daily O-D Matrix Estimation Using Cellular Probe Data[C]. Transportation Research Board 89th Annual Meeting, Washington D C, USA, 2010
[38] Newman M E, Girvan M. Finding and Evaluating Community Structure in Networks[J]. Physical Review E, Statistical Nonlinear and Soft Matter Physics, 2004, 69(2):026113 doi: 10.1103/PhysRevE.69.026113
[39] Blondel V D, Guillaume J L, Lambiotte R, et al. Fast Unfolding of Communities in Large Networks[J]. Journal of Statistical Mechanics Theory & Experiment, 2008(10):155-168 doi: 10.1088-1742-5468-2008-10-P10008/
[40] Hassan Z, Keymanesh M, Khabiri M. Locating Stations of Public Transportation Vehicles for Improving Transit Accessibility[J]. Transport, 2007, 22(2):99-104 http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_70ceb7a60892abd66f75b7e3452d2783
[41] Ford A C, Barr S L, Dawson R J, et al. Transport Accessibility Analysis Using GIS:Assessing Sustainable Transport in London[J]. ISPRS International Journal of Geo-Information, 2015, 4(1):124-149 doi: 10.3390/ijgi4010124
[42] Karou S, Hull A. Accessibility Modelling:Predicting the Impact of Planned Transport Infrastructure on Accessibility Patterns in Edinburgh, UK[J]. Journal of Transport Geography, 2014, 35:1-11 doi: 10.1016/j.jtrangeo.2014.01.002
[43] Niedzielski M A, Boschmann E E. Travel Time and Distance as Relative Accessibility in the Journey to Work[J].Annals of the Association of American Geographers, 2014, 104(6):1156-1182 doi: 10.1080/00045608.2014.958398
[44] 钟水映.人口流动与社会经济发展[M].武汉:武汉大学出版社, 2000 Zhong Shuiying. Population Flow and Socioecono-mic Development[M]. Wuhan:Wuhan University Press, 2000
[45] 俞路, 张善余.基于空间统计的人口迁移流分析——以我国三大都市圈为例[J].华东师范大学学报(哲学社会科学版), 2005, 37(5):25-31 doi: 10.3969/j.issn.1000-5579.2005.05.005 Yu Lu, Zhang Shanyu. An Analysis of the Population Resettlement Flow Based on the Integration of Spatial Statistics-A Case Study of Three Megalopoleis in Contemporary China[J].Journal of East China Normal University (Philosophy and Social Sciences), 2005, 37(5):25-31 doi: 10.3969/j.issn.1000-5579.2005.05.005
-
期刊类型引用(16)
1. 胡华参，朱建军，付海强，Lopez-Sanchez Juan Manuel，Gómez Cristina，张涛，刘奎. 基于机器学习联合TanDEM-X InSAR和ICESat-2数据估计大范围林下地形. 遥感学报. 2025(01): 191-202 . 
百度学术
2. 朱珺，彭葳，付海强，欧蔓，雷缮诚，张石平. 观测信息不足条件下TanDEM-X InSAR的大范围林下地形反演. 测绘学报. 2024(05): 959-966 . 百度学术
3. 席梦丹，张永洪. 基于GF-3精细条带数据生产DEM的精度分析. 测绘标准化. 2024(02): 101-106 . 百度学术
4. 潘紫阳，万阿芳，王烽，邹明普，温康锋，邹梅芳. 基于时变多项式的InSAR轨道误差去除和DEM估计. 测绘通报. 2024(08): 145-150+176 . 百度学术
5. 卢丽君，郭玄烨，邱志伟，李俊峰. 复杂地形条件下的PolInSAR森林模拟和特征参数分析. 测绘科学. 2024(07): 112-122 . 百度学术
6. 牛宵，江娜. 基于地理国情监测数据的景观生态评价研究. 测绘标准化. 2024(03): 103-109 . 百度学术
7. 刘志卫，赵蓉，朱建军，付海强，周璀，周亦. 联合TanDEM-X DEM与Sentinel-2多光谱数据的林下地形提取. 遥感学报. 2024(10): 2540-2552 . 百度学术
8. 张永鑫，张王菲，姬永杰，赵含. 星载X波段InSAR数据森林高度估测及反演研究. 武汉大学学报(信息科学版). 2024(12): 2279-2289 . 百度学术
9. 邓兴升，王清阳. 基于梯度提升决策树的植被高度模型研究. 长沙理工大学学报(自然科学版). 2023(01): 65-74 . 百度学术
10. 欧蔓，朱建军，张涛，胡华参. 双基线TanDEM-X InSAR数据的森林高度反演. 测绘工程. 2023(03): 8-13 . 百度学术
11. 周鹏，刘伟，李鹏飞，刘战合，赵青，阎泽恒. 基于广义似然比不规则建筑物层析SAR成像方法. 物联网技术. 2022(06): 84-88+92 . 百度学术
12. 朱建军，付海强，汪长城. 极化干涉SAR地表覆盖层“穿透测绘”技术进展. 测绘学报. 2022(06): 983-995 . 百度学术
13. 田方. 基于激光点云数据的地形地貌测量研究. 电子测量技术. 2021(01): 15-19 . 百度学术
14. 岳国栋，刘扬，郑春雨. 地形复杂区域InSAR技术的DEM融合方法研究. 测绘与空间地理信息. 2021(07): 205-209+214 . 百度学术
15. 熊俊麟，范宣梅，窦向阳，杨莹辉. 藏东南然乌湖流域雅弄冰川流速季节性变化. 武汉大学学报(信息科学版). 2021(10): 1579-1588 . 百度学术
16. 刘雅佳，朱建军，付海强. 双极化干涉SAR林下地形反演. 测绘工程. 2020(05): 20-26 . 百度学术
其他类型引用(19)
计量
- 文章访问数: 2499
- HTML全文浏览量: 259
- PDF下载量: 400
- 被引次数: 35
