基于夜光遥感影像的“一带一路”沿线国家城市发展时空格局分析

李德仁, 余涵若, 李熙

李德仁, 余涵若, 李熙. 基于夜光遥感影像的“一带一路”沿线国家城市发展时空格局分析[J]. 武汉大学学报 ( 信息科学版), 2017, 42(6): 711-720. DOI: 10.13203/j.whugis20170100
引用本文: 李德仁, 余涵若, 李熙. 基于夜光遥感影像的“一带一路”沿线国家城市发展时空格局分析[J]. 武汉大学学报 ( 信息科学版), 2017, 42(6): 711-720. DOI: 10.13203/j.whugis20170100
LI Deren, YU Hanruo, LI Xi. The Spatial-Temporal Pattern Analysis of City Development in Countries along the Belt and Road Initiative Based on Nighttime Light Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(6): 711-720. DOI: 10.13203/j.whugis20170100
Citation: LI Deren, YU Hanruo, LI Xi. The Spatial-Temporal Pattern Analysis of City Development in Countries along the Belt and Road Initiative Based on Nighttime Light Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(6): 711-720. DOI: 10.13203/j.whugis20170100

基于夜光遥感影像的“一带一路”沿线国家城市发展时空格局分析

基金项目: 

湖北省自然科学基金 2014CFB726

间数据挖掘与信息共享教育部重点实验室(福州大学)开放基金 2016LSDMIS03

015测绘地理信息公益性行业科研专项基 201512026

中央高校基本科研业务费专项资金 2042016kf0162

详细信息
    作者简介:

    李德仁, 博士, 教授, 中国科学院院士, 中国工程院院士, 欧亚科学院院士, 主要研究方向为摄影测量与遥感。drli@whu.edu.cn

    通讯作者:

    李熙, 博士, 副教授。 E-mail:lixi@whu.edu.cn

  • 中图分类号: P231

The Spatial-Temporal Pattern Analysis of City Development in Countries along the Belt and Road Initiative Based on Nighttime Light Data

Funds: 

The National Natural Science Foundation of Hubei Province 2014CFB726

Open Research Fund Program of Key Laboratory of Spatial Data Mining and Information Sharing of the Ministry of Education (Fuzhou University) 2016LSDMIS03

Special Scientific Research Fund of Public Welfare Profession on Surveying, Mapping and GeoInformation 201512026

Fundamental Research Funds for the Central Universities 2042016kf0162

More Information
    Author Bio:

    LI Deren, PhD, professor, academician of Chinese Academy of Science, Chinese Academy of Engineering, and Eurasian Academy of Sciences, specialize in photogrammetry and remote sensing. drli@whu.edu.cn

    Corresponding author:

    LI Xi, associate PhD, professor. E-mail:lixi@whu.edu.cn

  • 摘要: 基于1993~2012年长时间序列夜光遥感数据,利用空间统计、空间标准差椭圆和位序-规模分布方法研究“一带一路”(the Belt and Road Initiative,B & R)沿线国家的城市体系空间结构和时空演化趋势。结果表明:“一带一路”沿线各国城市夜间灯光增长最快的是进行经济改革的发展中国家和战后重建的国家,城市夜间灯光锐减主要发生在社会经济动荡地区。“一带一路”区域内城市夜间灯光的规模分布总体呈现空间扩张趋势,夜间灯光重心向东南亚方向移动。“一带一路”区域内城市夜光规模前2 000位的城市夜间灯光符合位序-规模准则,高位序城市夜间灯光规模不够突出而中小城市发育相对较好,集中趋势大于分散趋势。
    Abstract: Based on the DMSP/OLS nighttime light data for the years 1993-2012 and spatial analysis methods including standard deviational ellipse and rank-size distribution, this paper systematically analyzes the spatial structure and spatiotemporal dynamics of the urban system in countries along B & R (The Belt and Road Initiative). We found that nighttime light increased in most countries along B & R. These fast growing countries are undergoing economic reforms and post-war reconstruction, while nighttime light reduction occurs in areas of social and economic unrest. The trend of size distribution of nighttime light in B&R is continuous spatial expansion, and the center of the nighttime light is moving to southeast Asia. The nighttime light distribution in the top 2000 urban places in B&R follows the rank-size distribution, thus urban land distribution is more concentrated than the past. The high rank cities are fairly well developed, but the development of the small cities is lagging behind. The general distribution trend toward, concentration is stronger than decentralization in B&R.
  • 低轨卫星在遥感、通信、导航等一系列科学研究和工程应用中发挥着不可或缺的作用。低轨卫星精密位置信息是完成上述任务的前提和基础。自GPS成功应用于TOPEX卫星精密定轨后[1-2],利用GPS进行低轨卫星定轨已广泛应用于重力场恢复和气候试验卫星(gravity recovery and climate experiment, GRACE)[3]、挑战性微型卫星载荷(challenging mini-satellite payload, CHAMP)[4-5]、海洋[6]、风云[7-8]等卫星。利用GPS精密轨道和钟差产品,低轨卫星浮点解定轨精度可达1~3 cm[9-10]

    固定载波相位模糊度[11-14]能够进一步提升定轨精度。为实现单接收机模糊度固定,国内外学者先后提出了小数周偏差法[15-16]、整数相位钟法[17]和解耦钟差法[18]。文献[19]利用美国喷气推进实验室(Jet Propulsion Laboratory, JPL)提供的相位偏差产品固定模糊度,实现Jason-2卫星径向定轨精度优于1 cm;文献[20]研究了固定模糊度的GRACE卫星运动学定轨,相对于浮点解,三轴精度提升19%~50%;文献[11]利用法国国家空间研究中心(Centre National d’Etudes Spatiales, CNES)发布的整数钟产品[21]对哨兵(Sentinel)-3A卫星模糊度进行固定,轨道重叠弧段和卫星激光测距(satellite laser ranging,SLR)检核精度均有明显提升;文献[14]实现了基于单接收机模糊度固定的GRACE后续星(GRACE follow‑on, GRACE-FO)简化动力学定轨,相对于浮点解,定轨精度提升18%;文献[22]利用CNES整数钟产品实现GRACE-FO卫星简化动力学和运动学定轨中模糊度参数的固定,相对于浮点解,固定解轨道重叠弧段精度提升20%。

    上述研究多基于星载GPS观测数据,文献[7-8]研究了星载GPS/北斗卫星导航系统(BeiDou navigation satellite system, BDS)组合定轨,但尚未实现模糊度固定。Sentinel-6A卫星于2020年11月21日发射,旨在提供连续的海面测高数据,其轨道高度为1 336 km,轨道倾角为66°[23-24]。为满足精密定轨需求,Sentinel-6A卫星搭载了多普勒无线电定位定轨系统(Doppler orbitography and radio positioning integrated by sate‑llite, DORIS)、GPS/Galileo双模接收机和SLR反射棱镜[13]。星载GPS/Galileo数据为研究基于多模全球导航卫星系统(global navigation satellite system, GNSS)和单接收机模糊度固定的低轨卫星精密定轨提供了契机。文献[13]分别利用欧洲定轨中心(Center for Orbit Determination in Europe, CODE)和CNES产品开展多种模式的Sentinel-6A精密定轨,结果表明,固定模糊度后,单GPS与单Galileo简化动力学定轨结果与GPS/Galileo双系统组合定轨一致性优于1.0 cm;文献[25]利用星载GNSS数据开展基于广播星历的实时定轨研究,得益于Galileo较优的空间信号精度,Sentinel-6A三维实时定轨精度优于10 cm。

    由于不依赖力学模型,运动学精密定轨信息可直接用于地球重力场反演,具有重要的科学意义和研究价值[20],同时运动学定轨算法简单、计算资源需求较小,适用于大规模星座自主运行管理。本文利用Sentinel-6A卫星实测数据开展精密定轨,在分析星载GPS、Galileo观测数据质量的基础上,分别利用CODE[26]、CNES[21]、德国地学研究中心[27](German Research Centre for Geosciences, GFZ)和中国武汉大学[28](Wuhan University, WHU)发布的观测值偏差(observation specific biases, OSB)及对应的轨道和钟差产品实现单接收机模糊度固定,研究多模GNSS以及模糊度固定对低轨卫星运动学定轨精度的影响。

    GNSS伪距和载波相位观测值可表示为[29]

    Pr,js=ρrs+c(dtr-dts)+Ir,js+Trs+br,j-    bjs+εLr,js=ρrs+c(dtr-dts)-Ir,js+Trs+λj(Nr,js+    Br,j-Bjs)+ξ (1)

    式中,Pr,jsLr,js分别表示伪距和载波相位观测值;rs分别表示星载接收机和卫星;j表示信号频率;ρrs为接收机和卫星间的几何距离;c为光速,dtrdts分别为接收机和卫星钟误差;Ir,js为电离层延迟;Trs表示对流层延迟,对于低轨观测值,该项为零;br,jbjs分别为接收机和卫星端伪距硬件偏差;Br,jBjs为载波相位硬件偏差,以周表示;λj为频率j对应的载波波长;Nr,js为载波相位模糊度;εξ分别为伪距和载波观测值噪声。

    构建无电离层组合观测值,以伪距为钟差基准,此时伪距硬件偏差被钟差吸收,式(1)可写为:

    Pr,IFs=ρrs+c(dt¯r-dt¯s)Lr,IFs=ρrs+c(dt¯r-dt¯s)+λ1N¯r,IFsdt¯r=dtr+br,IF/cdt¯s=dts+bIFs/cN¯r,IFs=Nr,IFs+(Br,IF-br,IF/λ1)-(BIFs-bIFs/λ1) (2)

    式中,Pr,IFsLr,IFs为无电离层组合伪距和载波相位观测值;d¯trd¯ts分别为吸收了伪距硬件偏差的接收机和卫星钟差;N¯r,IFs为利用伪距和载波相位观测值定位定轨时解算得到的浮点解模糊度;Nr,IFs为无电离层组合模糊度;br,IFbIFs分别为接收机和卫星端无电离层组合伪距硬件偏差;Br,IFBIFs为无电离层组合载波相位偏差。

    无电离层组合模糊度N¯r,IFs可分解为宽巷(wide-lane,WL)和窄巷(narrow-lane,NL)模糊度[29]

    N¯r,IFs=f1f2f12-f22Nr,WLs+f1f1+f2(Nr,NLs+    dr,NL-dNLs)Nr,WLs=Nr,1s-Nr,2sNr,NLs=Nr,1sdr,NL=f1+f2f1(Br,IF-br,IF/λ1)dNLs=f1+f2f1(BIFs-bIFs/λ1) (3)

    式中,f1f2分别为载波L1L2的频率;Nr,WLsNr,NLs分别为整数宽巷和窄巷模糊度;Nr,1sNr,2s分别为载波L1L2的模糊度;dr,NLdNLs分别为接收机和卫星端伪距和载波硬件偏差的窄巷组合。整数宽巷模糊度Nr,WLs可利用Melbourne-Wübbena组合[30-31]解算:

    (f1Lr,1s-f2Lr,2sf1-f2-f1Pr,1s+f2Pr,2sf1+f2)/λWL=Nr,WLs+dr,WL-dWLs (4)

    式中,dr,WLdWLs分别为接收机和卫星端宽巷硬件偏差;λWL为宽巷波长。分别对接收机和卫星端宽巷偏差进行改正,可实现宽巷模糊度固定。宽巷模糊度固定后,代入式(3)中可得:

    f2f1-f2Nr,WLs-f1+f2f1N¯r,IFs=Nr,NLs+dr,NL-dNLs (5)

    对接收机和卫星端窄巷偏差进行改正可以实现窄巷模糊度固定。实际处理时,卫星端宽巷和窄巷偏差可利用地面站网解算并向用户提供,用户端偏差可采用星间单差消除。固定星间单差宽巷和窄巷模糊度后,可根据式(6)恢复无电离层组合模糊度:

    ΔNr,IFs1,s2=f1f2f12-f22ΔNr,WLs1,s2+f1f1+f2(ΔNr,NLs1,s2-ΔdNLs1,s2) (6)

    式中,s1s2为导航卫星;ΔNr,IFs1,s2为单差无电离层组合模糊度;ΔNr,WLs1,s2ΔNr,NLs1,s2分别为已固定的单差宽巷和窄巷模糊度;ΔdNLs1,s2为卫星端窄巷偏差的星间单差值。之后,利用恢复的单差无电离层组合模糊度作为虚拟观测值对非差模糊度进行约束,得到固定解轨道。

    Sentinel-6A卫星搭载了澳大利亚RUAG公司开发的PODRIX GPS/Galileo双模接收机[13],支持接收L1/E1、L2和L5/E5a信号,根据不同类型的GPS卫星输出不同的伪距观测值,具体观测数据类型及对应的接收机独立交换格式版本3(receiver independent exchange formation version 3, Rinex 3)标识见表1。由于缺乏Galileo E14和E18卫星的历书数据,接收机未跟踪该两颗卫星的信号[13]

    表  1  Sentinel-6A接收机输出的GPS和Galileo观测值类型
    Table  1.  GPS and Galileo Observation Types Supported by Sentinel-6A Receiver
    GNSS卫星伪距载波
    GPS IIRC1C、C1W、C2WL1C、L2W
    GPS IIR-M、IIF、ⅢC1C、C2LL1C、L2L
    GalileoC1C、C5QL1C、L5Q
    下载: 导出CSV 
    | 显示表格

    利用2021年年积日第175—250天的GPS C1W、C2W、L1C、L2W和Galileo C1C、C5Q、L1C、L5Q观测数据进行运动学定轨。利用CODE的OSB产品将GPS IIR-M、IIF、Ⅲ卫星的C1C、C2L、L2L数据转换成C1W、C2W、L2W观测值。Sentinel-6A星载GNSS天线参考点、相位中心偏差(phase center offset, PCO)、卫星质心和SLR反射器在星固系下的坐标如表2所示[32]

    表  2  星载GNSS天线PCO、卫星质心和SLR反射器在星固系下的坐标/mm
    Table  2.  Positions of Onboard GNSS Antenna, SLR Retroreflector and Center of Mass in Satellite Reference Frame/mm
    设备名称XYZ
    GNSS天线参考点+2 474.8+0.1-1 080.3
    GNSS天线GPS PCO00+75.0
    GNSS天线Galileo PCO00+93.0
    SLR 反射棱镜+1 624.8-400.6+664.8
    卫星质心+1 533.0-7.0+37.0
    下载: 导出CSV 
    | 显示表格

    伪距观测质量将影响定轨精度和模糊度固定效果[33]。利用Anubis软件[34]对GPS和Galileo

    观测值伪距多径进行分析,结果如图1所示。可以看出,GPS C1W和C2W伪距多径均值分别为12.7和15.3 cm,Galileo C1C和C5Q伪距多径均值分别为7.4和3.9 cm,主要是由于Galileo E1和E5信号采用二进制偏置载波调制,相对于GPS采用的二进制相移键控调制,具有更强的抗干扰和多路径抑制能力[35]。Galileo较高伪距精度有利于后续模糊度的固定。

    图  1  GPS C1W/C2W和Galileo C1C/C5Q伪距多径
    Figure  1.  Code Multipath of GPS C1W/C2W and Galileo C1C/C5Q Observations

    利用Sentinel-6A星载GPS/Galileo数据开展3种模式的运动学定轨,包括GPS单系统、Galileo单系统以及GPS/Galileo双系统组合定轨,其中每种定轨模式根据模糊度是否固定又分为浮点解和固定解。GPS和Galileo卫星轨道和钟差分别采用CODE[26]、CNES[21]、GFZ[27]以及WHU[28]分析中心发布的多模GNSS实验(multi-GNSS experiment, MGEX)产品,同时利用OSB产品对载波相位偏差进行改正,以实现模糊度的固定。定轨时首先基于无电离层组合观测值得到浮点解轨道,之后根据§1介绍的方法依次固定宽巷和窄巷模糊度并恢复无电离层组合模糊度,最终得到固定解轨道。Sentinel-6A卫星运动学定轨策略见表3。由于目前尚无机构提供Sentinel-6A精

    表  3  Sentinel-6A卫星运动学定轨策略
    Table  3.  Sentinel-6A Precise Orbit Determination Strategy Based on Kinematic Method
    模型参数说明
    使用数据和产品观测数据非差无电离层组合
    定轨弧长/h24
    采样间隔/s10
    截止高度角/(°)3
    GNSS轨道和钟差CODE, CNES, GFZ, WHU MGEX产品
    GNSS相位偏差CODE, CNES, GFZ, WHU OSB产品
    GNSS天线igs14.atx
    接收机天线改正
    相位缠绕改正
    相对论改正IERS 2010
    参数估计卫星位置随机游走,过程噪声5 m/s,每历元估计
    接收机钟差随机游走,过程噪声30 m/s,每历元估计
    模糊度常数,每跟踪弧段估计1个
    下载: 导出CSV 
    | 显示表格

    密轨道产品,为评估不同模式运动学定轨精度,利用CODE产品开展GPS/Galileo组合简化动力学定轨并固定模糊度参数。

    可视卫星数和位置精度衰减因子(position dilution of precision, PDOP)与运动学定轨精度密切相关,首先分析GPS/Galileo双系统组合对可视卫星数和PDOP值的影响。对75 d的观测数据进行统计,GPS平均可视卫星数为7.5颗,Galileo为6.4颗,主要是因为Galileo系统卫星数较少,且Sentinel-6A接收机未跟踪E14和E18卫星[13]图2为运动学定轨中GPS单系统和GPS/Galileo组合系统的PDOP值分布,GPS单系统平均PDOP为2.2,Galileo单系统为2.5,双系统组合PDOP为1.5,显著提升了定轨的几何构型。

    图  2  PDOP值全球分布
    Figure  2.  Global Distribution of PDOP

    观测值残差能够反映所使用的模型和测量值之间的一致性。图3为使用CODE产品时GPS和Galileo运动学固定解定轨中伪距和载波相位残差均方根(root mean square,RMS)的全球分布图。图3(a)和图3(c)为伪距残差RMS,其中GPS为0.85 m,而Galileo为0.35 m,明显优于GPS,与图1中伪距多径分析结果一致。图3(b)和图3(d)为载波残差RMS值,GPS为5.5 mm,Galileo为4.8 mm,较小的观测值残差显示出Galileo较优的信号性能。

    图  3  运动学固定解定轨GPS和Galileo伪距和载波残差RMS分布
    Figure  3.  Code and Carrier Phase Residual RMS of GPS and Galileo Derived from Kinematic Orbit Determination with Ambiguity Resolution

    模糊度固定率为定轨弧段内固定为整数的模糊度数量与总的模糊度个数的比值,可以反映算法模型和定轨结果的正确性及精度。利用取整法[36]对浮点模糊度进行固定,宽巷、窄巷模糊度固定阈值均为0.2周。图4(a)、4(b)分别给出了使用CODE产品Sentinel-6A卫星宽巷和窄巷模糊度的固定率。红色、绿色和蓝色点分别表示GPS单系统(G_kin)、Galileo单系统(E_kin)和GPS/Galileo双系统组合(GE_kin)运动学定轨模糊度固定结果。GPS和Galileo宽巷模糊度固定率分别为99.6%和99.9%,双系统组合固定率为99.7%;GPS和Galileo窄巷固定率分别为93.5%和95.2%,双系统组合可进一步提升窄巷固定率至96.7%,原因在于双系统组合提升了运动学定轨精度以及浮点模糊度的估计精度。图4中同时给出了GPS/Galileo双系统组合动力学参考轨道(GE_dyn)的模糊度固定结果,以粉色点表示,宽巷和窄巷模糊度固定率分别为99.8%和97.9%。

    图  4  使用CODE产品不同定轨模式的宽巷和窄巷模糊度固定率
    Figure  4.  Wide Lane and Narrow Lane Ambiguity Fixing Rates of Different Orbit Determination Solutions Using CODE Product

    图5(a)和图5(b)为使用各分析中心产品时GPS单系统运动学定轨中宽巷和窄巷模糊度固定率,红色、绿色、蓝色和粉色点分别对应CODE、CNES、GFZ和WHU分析中心产品结果,

    图  5  基于各分析中心产品的GPS运动学定轨宽巷和窄巷模糊度固定率
    Figure  5.  GPS Wide Lane and Narrow Lane Ambiguity Fixing Rates for Kinematic Orbit Determination Using Products of Different Analysis Centers

    GFZ从2021年年积日第190天开始提供OSB产品。各分析中心GPS宽巷模糊度固定率均优于98%(CODE产品为99.6%,GFZ为99.5%,CNES为99.6,WHU为98.5%),窄巷固定率优于92%(CODE产品为93.5%,GFZ为94.0%,CNES为93.0%,WHU为92.2%),体现出较好的产品性能。WHU产品宽巷固定率相较于其他分析中心低1%左右。图6(a)为Galileo单系统的宽巷模糊度固定率,其中基于CODE和CNES产品的固定率为99.9%,GFZ产品为99.7%。基于WHU产品的Galileo宽巷模糊度固定率偏低,为94%,并最终影响了窄巷模糊度的固定。图6(b)为Galileo窄巷模糊度固定率,CODE和CNES为95.3%,GFZ为95.2%,WHU为89.5%,由于宽巷模糊度固定直接采用Melbourne-Wübbena组合观测值,尚未涉及GNSS卫星轨道和钟差,因此原因可能是伪距OSB精度影响了WHU产品宽巷模糊度的固定。

    图  6  基于各分析中心产品的Galileo运动学定轨宽巷和窄巷模糊度固定率
    Figure  6.  Galileo Wide Lane and Narrow Lane Ambiguity Fixing Rates for Kinematic Orbit Determination Using Products of Different Analysis Centers

    为进一步验证WHU精密轨道、钟差以及窄巷模糊度固定性能,图7分别给出了使用CODE和WHU产品时Galileo的模糊度总数(图7(a))、固定的宽巷模糊度数(图7(b))以及固定的窄巷与宽巷模糊度个数的比值(图7(c))。可以看出,使用CODE和WHU产品时Galileo模糊度总数基本一致,使用WHU产品固定的宽巷模糊度明

    图  7  使用CODE和WHU产品时Galileo运动学定轨模糊度总数、固定的宽巷模糊度数和固定的窄巷与宽巷模糊度数的比值
    Figure  7.  Galileo Total Ambiguity, Fixed Wide-Lane Ambiguity and Ratio of Fixed Narrow-Lane to Wide-Lane Ambiguity Using CODE and WHU Products

    显少于CODE产品,固定的窄巷与宽巷模糊度个数的比值为95%,说明WHU产品模糊度固定性能主要受伪距OSB影响,窄巷模糊度固定性能与CODE无明显差异。

    以GPS/Galileo双系统组合动力学固定解轨道为参考,对不同模式的运动学定轨结果进行评估。图8为2021年年积日第175天使用CODE产品运动学浮点解轨道与参考轨道在切向、法向和径向的差异,其中红色、蓝色和绿色点分别表示GPS单系统(G)、Galileo单系统(E)和GPS/Galileo组合(GE)定轨结果。可以看出,相对于单系统,双系统浮点解定轨精度提升较为明显,Galileo单系统在14时和20时处可视卫星数只有4颗,导致定轨差异较大,双系统组合增加了可视卫星数,定轨结果较为平稳。

    图  8  2021年年积日第175天运动学浮点解轨道与参考轨道的差异
    Figure  8.  Orbit Differences Between Kinematic Orbits with Float Ambiguity and Reference Orbits on Day of Year 175 in 2021

    图9为2021年年积日第175—250天运动学浮点解轨道与参考轨道在切向(图9(a))、法向(图9(b))和径向(图9(c))的差异,三轴差异统计值在表4中列出。与参考轨道相比,运动学浮点解定轨结果系统性偏差较小,为1~4 mm,表明动力学定轨中的力学参数较好地描述了卫星的受力情况。GPS单系统(G)浮点解轨道三维精度为38.4 mm,Galileo单系统(E)为57.1 mm,双系统组合(GE)明显提升了浮点解定轨精度,三维精度优于30 mm,相对于GPS提升超过20%。在分析时段内双系统浮点解轨道未出现明显离散点,提升了运动学定轨的可用性。

    图  9  使用CODE产品时运动学浮点解轨道与参考轨道的差异
    Figure  9.  Orbit Differences Between Kinematic Orbits with Float Ambiguity and Reference Orbits Using CODE Product
    表  4  使用CODE产品时运动学定轨结果与参考轨道差异统计值/mm
    Table  4.  Statistics of Orbit Differences Between Kinematic Orbits and Reference Orbits Using CODE Product/mm
    定轨模式切向法向径向三维
    GPS浮点解-3.1±22.5-0.1±18.9-1.1±24.538.4
    Galileo浮点解-2.7±30.5-3.7±31.53.1±36.257.1
    GPS/Galileo浮点解-2.9±16.2-2.2±15.11.5±17.728.6
    GPS固定解0.4±9.5-1.0±10.4-1.5±20.324.8
    Galileo固定解-0.8±9.2-0.7±12.01.7±23.828.3
    GPS/Galileo固定解-0.3±6.0-0.9±7.81.2±13.116.5
    下载: 导出CSV 
    | 显示表格

    为分析模糊度固定对运动学定轨精度的影响,图10给出了使用CODE产品时运动学固定解轨道与参考轨道的差异,三轴差异统计值也在表4中列出。与图9相比,固定模糊度后,定轨精度明显提升,浮点解轨道法向上的波动得到显著消除,同时与参考轨道切向和法向的系统性偏差在1 mm之内,径向小于2 mm。GPS(G)和Galileo(E)单系统固定解轨道三维精度分别为24.8和28.3 mm,相对于浮点解结果分别提升了35%和50%。双系统组合(GE)进一步提升定轨精度至16.5 mm,相对于GPS单系统提升超30%。

    图  10  使用CODE产品时运动学固定解轨道与参考轨道的差异
    Figure  10.  Orbit Differences Between Kinematic Orbits with Fixed Ambiguity and Reference Orbits Using CODE Product

    分别利用各分析中心产品进行Sentinel-6A卫星运动学定轨并固定模糊度,包括GPS单系

    统、Galileo单系统以及GPS/Galileo双系统组合模式。图11为各分析中心不同模式运动学固定解轨道与参考轨道的三维差异。GPS/Galileo双系统组合明显提升了定轨精度,与参考轨道的三维差异优于20 mm,相对于GPS单系统平均提升30%。Galileo单系统定轨精度偏低,主要原因是Galileo可用卫星较少,导致DOP值偏大,基于WHU产品的Galileo固定解轨道由于模糊度固定率较低,导致与参考轨道差异较大。

    图  11  基于各分析中心产品的运动学固定解轨道与参考轨道的三维差异
    Figure  11.  Orbit Differences Between Kinematic Orbits with Fixed Ambiguity and Reference Orbits Using Products of Different Analysis Centers

    Sentinel-6A卫星搭载了SLR反射器,SLR测距值可作为独立手段对卫星定轨结果进行评估。Sentinel-6A卫星激光反射棱镜阵列由一个中心棱镜和其周围的8个棱镜组成,对于波长532 nm的激光可产生约49 mm的系统性误差[13],文中计算SLR残差时对该系统性偏差进行改正。选取11个质量较好的SLR测站(Arequipa、Beijing、 Graz、Greenbelt、Haleakala、Hartebeesthoek、Herstmonceux、Mount Stromlo、Papeete、Wettzell、Yarragadee)对运动学定轨结果进行评估。

    图12为使用CODE产品时GPS单系统(G_kin)、Galileo单系统(E_kin)以及GPS/Galileo双系统组合运动学(GE_kin)浮点解轨道的SLR残差。GPS和Galileo单系统浮点解轨道SLR残差的RMS分别为19和23 mm,双系统组合则优于15 mm,相对于GPS单系统提升了20%。图13为使用各分析中心产品时Sentinel-6A运动学固定解轨道SLR残差的RMS。与浮点解结果相比,固定模糊度可明显提升定轨精度。GPS单系统RMS为13~15 mm,GPS/Galileo双系统组合RMS为12~14 mm,相对单系统提升超10%。与前述分析一致,使用WHU产品的Galileo单系统定轨精度稍差,除此之外,各分析中心结果之间的差异小于1.5 mm。

    图  12  使用CODE产品时运动学浮点解轨道SLR残差
    Figure  12.  SLR Residuals of Kinematic Orbits with Float Ambiguity Using CODE Product
    图  13  基于各分析中心产品的固定解轨道SLR残差RMS
    Figure  13.  SLR Residuals RMS of Different Orbits with Fixed Ambiguity Using Products of Different Analysis Centers

    Sentinel-6A卫星搭载的GPS/Galileo双模接收机为研究基于多模GNSS的低轨卫星精密定轨提供了数据支撑,使用2021年年积日第175—250天在轨实测数据,开展GPS单系统、Galileo单系统以及GPS/Galileo组合运动学定轨。同时,利用CODE、CNES、GFZ和WHU发布的GNSS精密轨道、钟差和OSB产品,实现单接收机模糊度固定。以动力学固定解轨道为参考,GPS/Galileo双系统组合运动学浮点解和固定解轨道三维精度分别优于30 和20 mm,相对于GPS单系统结果分别提升了20%和30%。基于CODE、CNES和GFZ分析中心产品的GPS模糊度固定率优于93%,Galileo优于95%,利用WHU产品的模糊度固定率则相对较低,具体原因需进一步分析。GPS单系统运动学固定解轨道SLR残差RMS为13~15 mm,双系统组合结果为12~14 mm,提升超过10%。相对于GPS单系统,GPS/Galileo组合能够显著提升低轨卫星运动学定轨的精度和可用性,对大规模低轨星座自主定轨和载荷研制具有借鉴意义。

  • 图  1   “一带一路”沿线国家示意图(资料来源中国新闻网)

    Figure  1.   "B&R" Framework

    图  2   “一带一路”沿线国家区域夜光遥感影像(WGS1984大地坐标系)

    Figure  2.   Nighttime Light Image in "B&R" (WGS1984)

    图  3   1993年~2012年夜光总量前20位国家

    Figure  3.   Nighttime Light for Top 20 Countries in 1993 and 2012

    图  4   1993年~2012年夜光增长率前20位与后20位国家

    Figure  4.   Growth Rates for Top and Last 20 Countries from 1993 to 2012

    图  5   1993~2012年夜光增量各国百分比

    Figure  5.   Percentage of Countries from 1993 to 2012

    图  6   1993~2012年夜光增量各地区百分比

    Figure  6.   Percentage of Subareas from 1993 to 2012

    图  7   “一带一路”夜间灯光规模标准差椭圆及其重心变动(古德陆地投影)

    Figure  7.   Specific Ellipse of Nighttime Light Distribution in "B&R"(Projection:Goode Homolosine (Land))

    图  8   分区域夜光规模标准差椭圆及其重心变动(古德陆地投影)

    Figure  8.   Specific Ellipse of Nighttime Light Distribution in Subareas(Projection:Goode Homolosine (Land))

    图  9   “一带一路”区域1993与2012年夜光城市体系位序-规模双对数图

    Figure  9.   Cities' Rank-size Plots of Nighttime Light in "B&R" for 1993 and 2012

    图  10   “一带一路”前2000区域位序-规模|q|值时间序列图

    Figure  10.   |q| of Rank-size Double Logarithm in Top 2000 Areas of "B&R" from 1993 to 2012

    表  1   “一带一路”沿线国家范围

    Table  1   Countries in "B&R" Area

    主要地区国家
    蒙中俄蒙古、俄罗斯、中国
    中亚5国哈萨克斯坦、吉尔吉斯斯坦、塔吉克斯坦、乌兹别克斯坦、土库曼斯坦
    南亚8国印度、巴基斯坦、孟加拉国、阿富汗、尼泊尔、不丹、斯里兰卡、马尔代夫
    东南亚11国越南、老挝、柬埔寨、泰国、马来西亚、新加坡、印度尼西亚、文莱、菲律宾、缅甸、东帝汶
    西亚非洲20国土耳其、伊朗、叙利亚、伊拉克、阿联酋、沙特阿拉伯、卡塔尔、巴林、科威特、黎巴嫩、阿曼、也门、约旦、以色列、巴勒斯坦、亚美尼亚、格鲁吉亚、阿塞拜疆、埃及、肯尼亚
    欧洲20国波兰、捷克、斯洛伐克、匈牙利、斯洛文尼亚、克罗地亚、罗马尼亚、保加利亚、塞尔维亚、黑山、马其顿、波黑、阿尔巴尼亚、爱沙尼亚、立陶宛、拉脱维亚、乌克兰、白俄罗斯、摩尔多瓦、希腊
    下载: 导出CSV

    表  2   子区域标准差椭圆变化

    Table  2   Specific Ellipse of Nighttime Light Distribution in Subareas

    主要地区年份范围/km2重心坐标空间增长
    率/%
    空间
    变化
    重心移动
    距离/km
    重心主要
    移动方向
    中亚1993
    2012
    1 936 862
    2 061 522
    (67°57′E,44°00′N)
    (68°31′E,44°27′N)
    6扩张70西南
    南亚1993
    2012
    2 511 894
    2 429 953
    (76°57′E,23°00′N)
    (77°06′E,22°14′N)
    -3收缩95
    东南亚1993
    2012
    4 977 804
    4 497 249
    (105°29′E,5°14′N)
    (105°01′E,6°41′N)
    -10收缩178西北
    西亚埃及1993
    2012
    3 502 919
    3 572 633
    (43°03′E,30°57′N)
    (43°27′E,30°28′N)
    2扩张69东南
    欧洲1993
    2012
    1 422 474
    1 238 823
    (22°44′E,48°17′N)
    (22°36′E,48°36′N)
    -13收缩167西
    中国1993
    2012
    3 253 791
    3 487 545
    (115°14′E,34°09′N)
    (114°04′E,33°32′N)
    7扩张86西南
    下载: 导出CSV
  • [1] 明浩. "一带一路"与"人类命运共同体[J].中央民族大学学报(哲学社会科学版), 2015, 42(6): 23-30 http://www.cnki.com.cn/Article/CJFDTOTAL-ZYMD201506003.htm

    Ming Hao. "One Belt and One Road" and "Community of Human Destiny"[J].Journal of University of China Minzu (Philosophy and Social Sciences Edition), 2015, 42(6): 23-30 http://www.cnki.com.cn/Article/CJFDTOTAL-ZYMD201506003.htm

    [2] 匡文慧, 陈利军, 刘纪远, 等.亚洲人造地表覆盖遥感精细化分类与分布特征分析[J].中国科学:地球科学, 2016, 46(09): 1 162-1 179 http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201609002.htm

    Kuang Wenhui, Chen Lijun, Liu Jiyuan, et al. Remote Sensing-based Artificial Surface Cover Classification in Asia and Spatial Pattern Analysis[J].Science China, 2016, 46(09): 1 162-1 179 http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201609002.htm

    [3] 成晨, 傅文学, 胡召玲, 等.基于遥感技术的近30年中亚地区主要湖泊变化[J].国土资源遥感, 2015, 27(1): 146-152 http://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201501024.htm

    Cheng Chen, Fu Wenxue, Hu Zhaoling, et al. Changes of Major Lakes in Central Asia over the Past 30 Years Revealed by Remote Sensing Technology[J]. Remote Sensing for Land & Resources, 2015, 27(1): 146-152 http://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201501024.htm

    [4] 李斌, 燕琴, 张丽, 等.长江中游洪涝灾害特征的MODIS时序监测与分析[J].武汉大学学报·信息科学版, 2013, 38(7): 789-792 http://ch.whu.edu.cn/CN/abstract/abstract2691.shtml

    Li Bin, Yan Qin, Zhang Li, et al. Flood Monitoring and Analysis over the Middle Reaches of Yangtze River Basin with MODIS Time-Series Imagery[J].Geomatics and Information Science of Wuhan University, 2013, 38(7): 789-792 http://ch.whu.edu.cn/CN/abstract/abstract2691.shtml

    [5] 李德仁, 李熙.论夜光遥感数据挖掘[J].测绘学报, 2015, 44(6): 591-601 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201506002.htm

    Li Deren, Li Xi. An Overview on Data Mining of Nighttime Light Remote Sensing[J].Acta Geodaetica et Cartographica Sinica, 2015, 44(6): 591-601 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201506002.htm

    [6]

    Li Xi, Ge Linlin, Chen Xiaoling. Detecting Zimbabwe's Decadal Economic Decline Using Nighttime Light Imagery[J].Remote Sensing, 2013, 5(9): 4 551-4 570 doi: 10.3390/rs5094551

    [7]

    Xie Yanhua, Weng Qihao. Updating Urban Extents with Nighttime Light Imagery by Using an Object-Based Thresholding Method[J].Remote Sensing of Environment, 2016, 187: 1-13 doi: 10.1016/j.rse.2016.10.002

    [8]

    Marcantonio M, Pareeth S, Rocchini D, et al. The Integration of Artificial Night-Time Lights in Landscape Ecology: A Remote Sensing Approach[J].Ecological Complexity, 2015, 22: 109-120 doi: 10.1016/j.ecocom.2015.02.008

    [9]

    He Chunyang, Liu Zhifeng, Tian J I E, et al. Urban Expansion Dynamics and Natural Habitat Loss in China: A Multiscale Landscape Perspective[J].Global Change Biology, 2014, 9(20): 2 886-2 902 http://www.researchgate.net/profile/Zhifeng_Liu4/publication/260916582_Urban_expansion_dynamics_and_natural_habitat_loss_in_China_a_multiscale_landscape_perspective/links/0c960534a3d886657b000000.pdf?disableCoverPage=true

    [10]

    Gao Bin, Huang Qingxu, He Chunyang, et al. Dynamics of Urbanization Levels in China From 1992 to 2012: Perspective From DMSP/OLS Nighttime Light Data[J].Remote Sensing, 2015, 7(2): 1 721-1 735 doi: 10.3390/rs70201721

    [11]

    Shi Kaifang, Chen Yun, Yu Bailang, et al. Modeling Spatiotemporal CO2 (Carbon Dioxide) Emission Dynamics in China From DMSP-OLS Nighttime Stable Light Data Using Panel Data Analysis[J]. Applied Energy, 2016, 168: 523-533 doi: 10.1016/j.apenergy.2015.11.055

    [12]

    Nghiem S V, Balk D, Rodriguez E, et al. Observations of Urban and Suburban Environments with Global Satellite Scatterometer Data[J].ISPRS Journal of Photogrammetry and Remote Sensing, 2009, 64(4): 367-380 doi: 10.1016/j.isprsjprs.2009.01.004

    [13]

    Yu Bailang, Shi Kaifang, Hu Yingjie, et al. Poverty Evaluation Using Npp-VIIRS Nighttime Light Composite Data at the County Level in China[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(3): 1-13 doi: 10.1109/JSTARS.2015.2414853

    [14]

    Small C, Elvidge C D. Night On Earth: Mapping Decadal Changes of Anthropogenic Night Light in Asia[J].International Journal of Applied Earth Observation and Geoinformation, 2013, 22: 40-52 doi: 10.1016/j.jag.2012.02.009

    [15]

    Zhou Naijun, Hubacek K, Roberts M. Analysis of Spatial Patterns of Urban Growth Across South Asia Using DMSP-OLS Nighttime Lights Data[J].Applied Geography, 2015, 63: 292-303 doi: 10.1016/j.apgeog.2015.06.016

    [16] 李德仁, 李熙.夜光遥感技术在评估经济社会发展中的应用——兼论其对"一带一路"建设质量的保障[J].宏观质量研究, 2015, 3(4): 1-8 http://www.cnki.com.cn/Article/CJFDTOTAL-HGLY201504002.htm

    Li Deren, Li Xi. Applications of Nighttime Light Remote Sensing Evaluating of Socioeconomic Development[J].Journal of Macro-Quality Research, 2015, 3(4): 1-8 http://www.cnki.com.cn/Article/CJFDTOTAL-HGLY201504002.htm

    [17]

    Ioannides Y M, Overman H G. Zipf's Law for Cities: An Empirical Examination[J].Regional Science and Urban Economics, 2003, 33(2): 127-137 doi: 10.1016/S0166-0462(02)00006-6

    [18] 顾朝林.中国城镇体系等级规模分布模型及其结构预测[J].经济地理, 1990, 10(03): 54-56 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL199003011.htm

    Gu Chaolin. The Scale Distribution Model and Structure Prediction of Urban System in China[J].Economic Geography, 1990, 10(03): 54-56 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL199003011.htm

    [19]

    Gan Li, Li Dong, Song Shunfeng. Is the ZIPF Law Spurious in Explaining City-Size Distributions[J].Economics Letters, 2006, 92: 256-262 doi: 10.1016/j.econlet.2006.03.004

    [20]

    Anderson G, Ge Ying. The Size Distribution of Chinese Cities[J].Regional Science and Urban Economics, 2005, 35(6): 756-776 doi: 10.1016/j.regsciurbeco.2005.01.003

    [21]

    Gangopadhyay K, Basu B. City Size Distributions for India and China[J].Physica a-Statistical Mechanics and Its Applications, 2009, 388: 2 682-2 688 doi: 10.1016/j.physa.2009.03.019

    [22] 张超, 王春杨, 吕永强, 等.长江经济带城市体系空间结构—基于夜间灯光数据的研究[J].城市发展研究, 2015, 22(03): 19-27 http://www.cnki.com.cn/Article/CJFDTOTAL-CSFY201503004.htm

    Zhang Chao, Wang Chunyang, Lu Yongqiang, et al. Research on City System Spatial Structure of the Yangtze River Economic Belt: Based on DMSP/DMSP Night Time Light Data[J].Urban Development Studies, 2015, 22(03): 19-27 http://www.cnki.com.cn/Article/CJFDTOTAL-CSFY201503004.htm

    [23] 张虹鸥, 叶玉瑶, 陈绍愿.珠江三角洲城市群城市规模分布变化及其空间特征[J].经济地理, 2006, 26(5): 806-809 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL200605018.htm

    Zhang Hong'ou, Ye Yuyao, Chen Shaoyuan. The Changes of Population Size Hierarchy of the Urban Group of Pearl River Delta and Its Spatial Character Since 20 Years Ago[J].Economic Geography, 2006, 26(5): 806-809 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL200605018.htm

    [24] 杨洋, 李雅静, 何春阳, 等.环渤海地区三大城市群城市规模分布动态比较——基于1992-2012年夜间灯光数据的分析和透视[J].经济地理, 2016, 36(4): 59-69 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL201604009.htm

    Yang Yang, LI Yajing, HE Chunyang, et al. A Comparative Study on Spatiotemporal Dynamics of City Size Distribution Among Three Urban Agglomerations in Bohai Rim: Based on the Analysis and Perspective of Nighttime Light Data from 1992 to 2012[J].Economic Geography, 2016, 36(4): 59-69 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL201604009.htm

    [25]

    Small C, Elvidge C D, Balk D, et al. Spatial Scaling of Stable Night Lights[J].Remote Sensing of Environment, 2011, 115(2): 269-280 doi: 10.1016/j.rse.2010.08.021

    [26]

    Small C, Pozzi F, Elvidge C. Spatial Analysis of Global Urban Extent from DMSP-OLS Night Lights[J].Remote Sensing of Environment, 2005, 96(3/4): 277-291 http://www.researchgate.net/publication/222409862_Spatial_analysis_of_global_urban_extent_from_DMSP-OLS_night_light

    [27]

    Jiang Bin, Yin Junjun, Liu Qingling. Zipf'S Law for all the Natural Cities Around the World[J].Geographical Information Science, 2015, 3(29): 498-522 https://www.researchgate.net/publication/260167281_Zipf%27s_Law_for_All_the_Natural_Cities_around_the_World

    [28] 邹嘉龄, 刘春腊, 尹国庆, 等.中国与"一带一路"沿线国家贸易格局及其经济贡献[J].地理科学进展, 2015, 34(5): 598-605 http://youxian.cnki.com.cn/yxdetail.aspx?filename=ZLXC2017030800L&dbname=CAPJ2015

    Zou Jialing, Liu chula, Yin Guoqing. Spatial Patterns and Economic Effects of China's Trade with Countries along the Belt and Road[J].Progress in Geography, 2015, 34(5): 598-605 http://youxian.cnki.com.cn/yxdetail.aspx?filename=ZLXC2017030800L&dbname=CAPJ2015

    [29]

    Elvidge C D, Ziskin D, Baugh K E, et al. A Fifteen Year Record of Global Natural Gas Flaring Derived from Satellite Data[J].Energies, 2009, 2(3): 595-622 doi: 10.3390/en20300595

    [30]

    Zhang Qingling, Pandey B, Seto K C. A Robust Method to Generate a Consistent Time Series from DMSP/OLS Nighttime Light Data[J].IEEE Transactions On Geoscience and Remote Sensing, 2016, 54(10): 5 821-5 831 doi: 10.1109/TGRS.2016.2572724

    [31]

    Elvidge C, Zhizhin M, Baugh K, et al. Methods for Global Survey of Natural Gas Flaring from Visible Infrared Imaging Radiometer Suite Data[J].Energies, 2016, 9(1): 14

    [32] 陈军, 陈晋, 廖安平, 等.全球30 m地表覆盖遥感制图的总体技术[J].测绘学报, 2014, 43(6): 551-557 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201406002.htm

    Chen Jun, Chen Jin, Liao Anping, et al. Concepts and Key Techniques for 30 m Global Land Cover Mapping[J].Acta Geodaetica et Cartographica Sinic, 2014, 43(6): 551-557 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201406002.htm

    [33] 赵璐, 赵作权.基于特征椭圆的中国经济空间分异研究[J].地理科学, 2014, 34(08): 979-986 http://www.cnki.com.cn/Article/CJFDTOTAL-DLKX201408011.htm

    Zhao Lu, Zhao Zuoquan. Projecting the Spatial Variation of Economic Based on the Specific Ellipses in China[J].Scientia Geographica Sinica, 2014, 34(08): 979-986 http://www.cnki.com.cn/Article/CJFDTOTAL-DLKX201408011.htm

    [34] 赵璐, 赵作权, 王伟.中国东部沿海地区经济空间格局变化[J].经济地理, 2014, 34(02): 14-18 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL201402003.htm

    Zhao Lu, Zhao Zuoquan, Wang Wei. The Spatial Pattern of Economy in Coastal Area of China[J].Economic Geography, 2014, 34(02): 14-18 http://www.cnki.com.cn/Article/CJFDTOTAL-JJDL201402003.htm

    [35] 杨开忠, 陈良文.中国区域城市体系演化实证研究[J].城市问题, 2008(03): 6-12 http://www.cnki.com.cn/Article/CJFDTOTAL-CSWT200803002.htm

    Yang Kaizhong, Chen Liangwen. An Empirical Study on the Evolution of Regional Urban System in China[J].Urban Problems, 2008(03): 6-12 http://www.cnki.com.cn/Article/CJFDTOTAL-CSWT200803002.htm

    [36] 舒松, 余柏蒗, 吴健平, 等.基于夜间灯光数据的城市建成区提取方法评价与应用[J].遥感技术与应用, 2011, 26(02): 169-176 http://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201102007.htm

    Shu Song, Yu Bolang, Wu Jianping, et al. Methods for Deriving Urban Built-up Area Using Night-light Data: Assessment and Application[J].Remoting Sensing Technology and Application, 2011, 26(02): 169-176 http://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201102007.htm

    [37] 吴健生, 赫胜彬, 彭建, 等.基于DMSP/OLS数据的城市发展空间特征研究[J].地理与地理信息科学, 2014, 30(2): 20-25 http://youxian.cnki.com.cn/yxdetail.aspx?filename=HDJJ201705003&dbname=CJFDPREP

    Wu Jiansheng, He Shengbin, Peng Jian, et al. Research on Spatial Characteristics of Urban Development Based on DMSP/OLS Data[J].Geography and GeoInformation Science, 2014, 30(2): 20-25 http://youxian.cnki.com.cn/yxdetail.aspx?filename=HDJJ201705003&dbname=CJFDPREP

    [38]

    Jiang Bin, Jia Tao. Zipf's Law for All the Natural Cities in the United States: A Geospatial Perspective[J].International Journal of Geographical Information Science, 2011, 25(8): 1 269-1 281 doi: 10.1080/13658816.2010.510801

    [39]

    Small C, Sousa D. Humans on Earth: Global Extents of Anthropogenic Land Cover from Remote Sensing[J].Anthropocene, 2016, 14: 1-33 doi: 10.1016/j.ancene.2016.04.003

图(10)  /  表(2)
计量
  • 文章访问数:  5479
  • HTML全文浏览量:  718
  • PDF下载量:  1645
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-04-10
  • 网络出版日期:  2023-07-26
  • 发布日期:  2017-06-04

目录

/

返回文章
返回