- 孙伟;刘经洲;
针对组合导航数据融合中的粒子滤波(PF)存在粒子退化现象，以及容积卡尔曼滤波(CKF)不能准确估计非高斯噪声下的状态等问题，该文提出容积粒子滤波(CPF)方案并将其应用于捷联惯性导航系统/全球定位系统(SINS/GPS)松组合数据处理。通过将容积卡尔曼滤波作为粒子滤波中的重要性密度函数对SINS/GPS松组合系统中的导航参数进行估计，基于容积滤波与粒子滤波对非线性非高斯的优势，得到该条件下的SINS/GPS松组合导航后验信息。搭建仿真平台，分别开展基于CPF、CKF及扩展粒子滤波(EPF)的组合系统估计，结果表明，容积粒子滤波可进一步提高SINS/GPS松组合系统的姿态、速度及位置参数估计精度。在对粒子滤波做出更好优化的同时也为SINS/GPS松组合导航的质量控制策略提供了一种新思路。
2022年09期 v.47;No.291 1-7页 [查看摘要][在线阅读][下载 607K] [下载次数：199 ] |[网刊下载次数：0 ] |[引用频次：3 ] |[阅读次数：0 ] - 祝会忠;李骏鹏;李军;
针对智能手机定位的研究大都集中于单频定位性能与算法，而缺乏对多系统多个频段定位实验的分析的问题，该文对智能手机GNSS多系统多个频段的观测数据进行质量分析，根据应用程序接口(API)信息，检核后的多普勒积分法与电离层残差法对异常数据进行处理，同时根据不同系统不同频段的观测数据质量，实时修改滤波观测值噪声方差以提高定位精度。开发了一款基于安卓平台的实时定位应用程序，通过实验对比分析了GPS、BDS、Galileo、GLONASS各个频段不同组合的动态定位性能。实验结果表明，GNSS多系统组合的收敛速度与定位精度最佳，固定解平面与高程方向精度最高可达0.04、0.07 m,浮点解精度收敛之后平面方向精度为0.15 m、高程方向精度为0.16 m。多系统第一频率(L1/B1/E1/G1)组合后定位结果较单系统平面方向精度提升71.3%,高程方向精度提升55.5%,在多系统第一频率基础上增加第二频率(L5/E5a)后，平面方向精度提升13%,高程方向精度提升32%。
2022年09期 v.47;No.291 8-19页 [查看摘要][在线阅读][下载 865K] [下载次数：381 ] |[网刊下载次数：0 ] |[引用频次：10 ] |[阅读次数：1 ] - 白羽;柳景斌;杨晨云鹂;姚尧;
针对室内活动场景下，额外布设定位基站成本较高，应用范围有限，基于手机的室内高精度位置服务需求难以满足的问题。考虑到Wi-Fi在全球网络中的广泛部署及智能手机内置多种运动传感器，本文实现了基于Wi-Fi RTT测距信息和手机内置多传感器的实时室内定位方法，并在室内进行静态和动态定位实验。实验结果表明：该算法无论是定位精度还是定位准度都优于仅依赖于Wi-Fi RTT测距信息的Wi-Fi RTT LS算法和Wi-Fi RTT UKF算法，平均定位误差为1.16 m,最大定位误差为2.14 m。
2022年09期 v.47;No.291 20-31页 [查看摘要][在线阅读][下载 698K] [下载次数：234 ] |[网刊下载次数：0 ] |[引用频次：4 ] |[阅读次数：0 ] - 李玉星;秘金钟;徐彦田;李建彬;
为了削弱中长基线差分后大气延迟残留误差对实时动态模糊度固定和定位精度的影响，该文提出了一种估计大气延迟参数的卡尔曼滤波算法，采用相对对流层天顶延迟和相对电离层天顶延迟作为残留大气延迟误差参数，与坐标向量和模糊度向量参数一起进行卡尔曼滤波估计。通过41 km参考站数据对算法进行验证。结果表明，BDS-3(B1C/B2a)与GPS精度相当，平面精度优于1.5 cm,高程精度优于2.7 cm; BDS-3(B1C/B2a)/GPS组合平面精度优于1.1 cm,高程精度优于2.7 cm,相比BDS-3(B1C/B2a)单系统平面精度提升了23.8%,估计大气延迟参数的中长基线卡尔曼滤波方法可有效提高模糊度固定成功率，提高中长基线条件下的定位精度。
2022年09期 v.47;No.291 32-37页 [查看摘要][在线阅读][下载 477K] [下载次数：250 ] |[网刊下载次数：0 ] |[引用频次：3 ] |[阅读次数：0 ] - 马志伟;纪兵;钟玉龙;
针对径向基函数建模框架下参考重力场模型对局部大地水准面精化的影响问题，通过移去不同截断阶次的参考场，分别构建大地水准面模型并进行了精度比较。结果显示，截断至420阶的XGM2016模型，更适合作为科罗拉多地区的参考重力场模型。其所得大地水准面，相对于GPS/水准和区域平均校准模型的误差RMS分别为2.2 cm和2.4 cm,比其他截断阶次下的大地水准面精度提高0.5～2 cm。另外，与其他机构的大地水准面模型相比，其精度也是较优的。该文所采用的确定参考场截断阶次的方法，可为局部地区高精度大地水准面模型的构建提供一定参考。
2022年09期 v.47;No.291 38-45页 [查看摘要][在线阅读][下载 409K] [下载次数：179 ] |[网刊下载次数：0 ] |[引用频次：0 ] |[阅读次数：0 ] - 张捍卫;杨永勤;张华;李晓玲;
缔合勒让德函数是地球科学领域中的主要数学工具之一。针对超高阶球谐重力场模型的构建与快速计算，以及跨阶次递推公式的适用性问题，该文利用勒让德函数的原理性公式，给出了3种类型的跨阶次递推公式。当次为0和1时，第1种采用了标准向前按列递推公式，第2和第3种则是采用了列式递推公式。其中，第一种是学界经常采用的公式。研究表明，在要求相对误差为1.0E-11情况下，3种类型的跨阶次递推公式都可以至少递推至15 000阶。在两个极点处，第2种和第3种算法的适用性优于第1种算法，但是第1种算法的计算精度优于后两种算法。
2022年09期 v.47;No.291 46-51+66页 [查看摘要][在线阅读][下载 701K] [下载次数：211 ] |[网刊下载次数：0 ] |[引用频次：0 ] |[阅读次数：2 ] - 岳彩亚;杨强;刘慧娟;
为了研究汶川地震后川滇地区地壳形变的时空分布带来的影响，该文以GPS速度场为约束，通过构建平移-旋转-应变模型计算得到区域内2008-06—2010-04、2010-05—2013-04和2013-05—2016-06这3期应变率场，并结合断裂带分布对不同期地壳应变场进行对比分析。结果表明：川滇地区地壳水平运动整体以继承性运动为主，强震的发生并未改变区域内整体运动趋势与作用力方向，且应变率场和面膨胀率的高值区域主要分布在断裂带间的交叉点或周边，变化较大区域主要位于几条深大断裂带及周边。汶川地震后区域内速度场随时间变化更加复杂且量级有所增加，且该地震对整个川滇地区地壳活动具有重要影响，而玉树地震和雅安地震主要对北半部分地壳起到一定调节作用，雅安地震后该区域应变累计增强，周边仍然有发生强震的危险性。
2022年09期 v.47;No.291 52-59页 [查看摘要][在线阅读][下载 607K] [下载次数：219 ] |[网刊下载次数：0 ] |[引用频次：0 ] |[阅读次数：0 ] - 晏慧能;戴吾蛟;温亚鑫;
针对GNSS位移时序数据中常存在由各种因素引起的异常值，使得难以准确提取GNSS站点的变形信号，而现有的GNSS位移时序处理方法主要考虑异常值的探测和变形特征参数的估计，但是对异常值进行修复的研究较少这一问题，该文提出了一种基于t检验和先验变形特征函数的GNSS位移时序处理方法，同时实现异常值的自动识别和修复以及变形特征参数的估计。利用该方法对旧金山湾区和南加州的GNSS位移时序数据进行了处理。结果表明，该方法在两个研究区域均能准确判断并有效修复GNSS位移时序中的异常值，提高GNSS位移时序对异常干扰信息的抵抗能力。将计算得到的E、N、U方向的变形速率与MIDAS算法进行对比，发现大部分GNSS站点的平面速率差异和高程速率差异都很小。对比未采用时序处理方法估计的变形速率值，发现该方法处理后的速率估值更接近MIDAS算法的结果，并且速率估值的方差更小，说明该方法能有效提高变形速率估计的精度和可靠性。
2022年09期 v.47;No.291 60-66页 [查看摘要][在线阅读][下载 426K] [下载次数：445 ] |[网刊下载次数：0 ] |[引用频次：5 ] |[阅读次数：0 ] - 周慧君;
针对目前航船轨迹预测方法中对影响船体运动及船员驾驶决策的海洋地理环境缺乏考虑的问题，提出了一种顾及胁迫反应机制的船舶轨迹预测方法，该方法不仅考虑到了航船轨迹本身的运动特性，更通过构建船舶航行中所处的实时海洋环境，将直接影响船体运动的海洋环境特征与影响船员决策的海洋环境特征应用到航船轨迹的预测中。结果表明：考虑航船所处海洋地理环境对船体运动状态的影响与船员面对环境威胁所采取的应对决策可以有效提升船舶轨迹预测的准确度，且本文提出的顾及胁迫反应机制的船舶轨迹预测方法较现有常见方法具有更好的预测能力。
2022年09期 v.47;No.291 67-75+94页 [查看摘要][在线阅读][下载 377K] [下载次数：183 ] |[网刊下载次数：0 ] |[引用频次：3 ] |[阅读次数：0 ] - 吴庆;刘长建;张熙;冯绪;王超;
针对PPP/INS组合导航初始对准阶段中静基座对准方法选择问题，研究了常见的解析法、间接双矢量、间接多矢量和Kalman精对准4种对准方法的原理及特点，并采用不同级别IMU的仿真和实测数据进行了对比分析。结果表明：相同IMU级别的静态数据下，间接双矢量和间接多矢量对准精度一致，Kalman精对准可以达到更高精度；间接多矢量比间接双矢量对准稍快，Kalman精对准的速度最慢；IMU级别越高，对准精度越高，无论何种对准方法和IMU级别，横滚角和俯仰角对准精度和速度一致，航向角对准精度最低，速度最慢；间接多矢量对准方法满足对准精度要求，对准速度最快并且对基座具有一定的角运动和线运动抗干扰能力，适于组合导航在实际中的应用。
2022年09期 v.47;No.291 76-85+111页 [查看摘要][在线阅读][下载 880K] [下载次数：187 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ] - 刘笃学;贺凯飞;丁志鹏;杨金权;
针对基于BDS-SMC的远海高精度动态定位中存在的数据中断或通讯延迟现象，从稳定性与精度两方面对GPS、GLONASS、Galileo与BDS-3的钟差改正数进行质量评估，并分析其数据结构特征。结果表明，在中断区间较短且IOD不变的情况下，由于钟差改正数的数量级保持在纳秒级，钟差改正数的预报误差低于钟差数据，故钟差改正数预报可有效应对该现象。针对BDS-SMC通讯成本高昂的问题，本文提出利用钟差改正数预报进行数据压缩，仅利用BDS-SMC传输模型参数可显著降低数据量。以原始数据序列与一次差分序列为训练样本，选取结构简单、参数较少的3种模型进行预报效果的评估。结果表明，由于预报区间较短，差分序列预报模型的误差累积效应不明显，预报精度高于原始序列预报模型。
2022年09期 v.47;No.291 86-94页 [查看摘要][在线阅读][下载 373K] [下载次数：52 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ]
- 田枭;林敦灵;刘洋;韩森;董梦;李航昊;
针对永定河流域(北京段)地表沉降及其与地下水位变化、生态补水之间的关系问题，利用小基线集合成孔径雷达干涉测量技术对覆盖研究区域的32景升轨雷达影像进行处理，提取了地表形变时空分布特征，分析了地表形变与地下水位变化的关系，并探讨了生态补水对区域地表沉降的影响。结果表明，流域上部具有一定的沉降区域，主要沉降速率为5～20 mm/a。地表形变与地下水位变化之间存在显著相关性：流域上部地下水位整体呈下降趋势，区域地表沉降相对严重；流域中下部地下水位整体呈上涨趋势，区域地表沉降得到遏制。生态补水对区域地表沉降有显著影响，但具体类型与区域地形地质、河道特征相关。
2022年09期 v.47;No.291 95-101+137页 [查看摘要][在线阅读][下载 563K] [下载次数：265 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ] - 张龙;黄国满;王章朋;杨书成;
针对近年来人们对海岸带的开发愈加频繁，海岸线的活动变化愈加剧烈的问题，该文基于合成孔径雷达在监测多云多雾的海岸带区域所具有得天独厚的优势，梳理了国内外SAR影像水边界线提取方法的研究进展及各种方法的提取结果的优劣比较，介绍了海岸线的基本概念与分类、SAR影像海岸线监测波段的选择、从SAR影像提取的瞬时海岸线到真实海岸线的校正过程及对未来SAR影像海岸线提取的研究方向进行了展望。结果表明，联合极化特征与神经网络的提取方法，综合考虑了多极化更丰富的目标信息和神经网络分类的准确性，为目前最优的提取方法。
2022年09期 v.47;No.291 102-111页 [查看摘要][在线阅读][下载 664K] [下载次数：737 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：3 ] - 向淇文;郭金城;汪娟;刘东烈;岳发政;
针对传统地质灾害隐患排查手段无法进行全天时、全天候、大范围识别监测的问题，该文以ALOS-2升轨和Sentinel-1A升降轨雷达影像为数据源，利用时序InSAR技术开展了贵州省六盘水和黔西南区域地质灾害隐患识别监测，同时基于SAR卫星入射观测方向和地形地貌之间的关系，对研究区域侧视成像几何可视性进行判别分析。最后综合利用3种数据获取的形变结果以及结合光学影像、地形等数据进行了形变解译，共识别出108个地质灾害疑似隐患点，并通过外业核查验证了地质灾害隐患识别的准确性。实验结果表明，时序InSAR技术能够有效进行地质灾害隐患识别监测，可为地质灾害防治提供科学依据。
2022年09期 v.47;No.291 112-119页 [查看摘要][在线阅读][下载 729K] [下载次数：744 ] |[网刊下载次数：0 ] |[引用频次：3 ] |[阅读次数：0 ] - 张立亭;孔文学;罗亦泳;邓先金;夏文生;
针对现有的遥感影像建筑物提取方法存在着效率低下、精度不高等问题，该文利用轻量型分割网络LinkNet框架构建出新的建筑物提取全卷积网络。设计三层卷积模块替换LinkNet中的残差层作为新的编码块，有效减少网络参数，加快了网络训练速度；融合增强感受野模块聚合多尺度上下文信息，有利于图像特征细节的恢复，从而提高网络分割精度；综合上述两点构建出基于深度学习的高性能建筑物自动提取网络。在相应建筑物数据集上进行实验结果表明，本文构建的全卷积网络比已有的建筑物提取网络SE-Unet综合预测精度更高，取得82.80%的均交并比和95.99%的召回率，同时，在提取建筑物的完整度、边界分割精度等方面具有较好的效果。
2022年09期 v.47;No.291 120-127+145页 [查看摘要][在线阅读][下载 750K] [下载次数：350 ] |[网刊下载次数：0 ] |[引用频次：4 ] |[阅读次数：0 ] - 潘春天;邱志伟;卢丽君;牛原;吴振宇;
针对高级陆地观测卫星2号提供的多极化数据源，不同极化方式影像的永久散射体点密度存在较大差别，而寻求极化最优的散射机制为永久散射体的识别和提取提供了一种新思路。该文以珠穆朗玛峰地区24景双极化影像为例，提出了改进后的模拟退火算法，退火过程采用不同的退火速率节约了退火总时间，通过不同的初始温度和稳定迭代次数实现双极化最优，然后利用干涉点目标分析技术反演地表时序形变，最后证实了双极化最优与HH极化形变结果的一致性。实验结果表明双极化最优结果实现了永久散射体数量最大化，能在大面积地表形变监测与灾害治理领域提供广泛的应用价值。
2022年09期 v.47;No.291 128-137页 [查看摘要][在线阅读][下载 966K] [下载次数：171 ] |[网刊下载次数：0 ] |[引用频次：0 ] |[阅读次数：0 ] - 何云海;国巧真;乔悦;吴欢欢;朱丽;
针对高空间分辨率卫星影像空间分辨率高、波段较少这一特点，该文以高分六号多光谱影像为数据源，通过分析各类地物的光谱特征曲线，建立了一种新指数。该指数通过增强植被和裸土信息，抑制不透水面信息来扩大不透水面和植被、裸土的差异，实现对不透水面的分离提取。验证阶段将该指数用于天津市河西区和河北省易县中心城区，并将提取结果与垂直不透水面指数(PISI)、比值居民地指数(RRI)的提取结果进行对比。结果表明：该方法在河西区和易县的提取效果都要优于PISI和RRI,并且两块实验区提取结果的总体精度分别为94.0%和91.7%,具有一定的稳定性。
2022年09期 v.47;No.291 138-145页 [查看摘要][在线阅读][下载 729K] [下载次数：267 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ] - 李梅;沈麒凯;陈启浩;刘修国;
针对合成孔径雷达(SAR)图像不同类型建筑物的区分问题，该文提出了一种基于U-Net的极化SAR图像建筑物分类方法。该方法将极化SAR数据的Pauli分解参数、规范化圆极化相关系数和G~0统计纹理参数作为U-Net的输入，建立建筑物分类U-Net模型，同时考虑建筑物的高度和单体面积的情况下，将建筑物分为高层、中层、低层小面积、低层厂房类大面积建筑物4类。对武汉市城区GF-3极化SAR数据的各类建筑物分类精度均在80%以上，最高达94.2%。该方法与仅使用单类别特征的U-Net网络方法以及卷积神经网络方法相比，分类结果更完整、建筑边界更准确，也更适合于复杂中心城区的不同类型建筑物的分类。
2022年09期 v.47;No.291 146-153+162页 [查看摘要][在线阅读][下载 915K] [下载次数：257 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ] - 李世伟;余洁;朱琳;韩丽华;陈蜜;
针对InSAR-BM3D滤波在干涉图条纹密集、噪声严重情况下，存在的对干涉条纹保护较差、易丢失细节信息的问题，该文利用圆周期中值滤波对条纹频率估计方法进行优化以提高条纹频率估计精度，并利用优化的条纹频率方法对InSAR-BM3D滤波方法进行改进。为验证本文改进滤波方法的有效性，分别进行模拟数据与真实数据的滤波实验，并将改进滤波方法与未改进InSAR-BM3D滤波以及常用的Lee滤波、Goldstein滤波、小波硬阈值滤波、小波软阈值滤波进行滤波结果的定性和定量比较。结果表明，改进滤波方法在保持InSAR-BM3D滤波高效抑制噪声优势的同时，提高了对干涉条纹边缘和细节的保护能力，能获得更好的滤波效果。
2022年09期 v.47;No.291 154-162页 [查看摘要][在线阅读][下载 1508K] [下载次数：224 ] |[网刊下载次数：0 ] |[引用频次：1 ] |[阅读次数：0 ] - 任超;刘桃林;
针对遥感影像进行水体提取存在普适性不高和难以取得最佳阈值问题，该文提出一种布谷鸟搜索算法与OTSU结合组成大津布谷鸟自适应阈值找寻方法(Cuckoo Search OTSU,CS-OTSU)。首先利用大津法进行多级阈值分割，然后在布谷鸟搜索算法的基础上增加一种新的步长搜索机制代替OTSU穷举式多级阈值搜索，并利用取得最佳阈值对遥感影像进行水体提取。通过实验表明，CS-OTSU可快速对不同环境进行水体提取，在精度上优于其他阈值方法，总体精度均大于97%,Kappa系数高于0.93。
2022年09期 v.47;No.291 163-169+180页 [查看摘要][在线阅读][下载 419K] [下载次数：393 ] |[网刊下载次数：0 ] |[引用频次：3 ] |[阅读次数：0 ] - 朱洪波;张在岩;秦育罗;宋伟东;张晋赫;
针对实际采集场景下路面影像中病害受背景纹理噪声影响程度大、病害边缘模糊导致分割不准确的问题，该文提出了一种基于Res＿UNet和全连接条件随机场的路面病害像素级精准检测方法：(1)对路面影像进行灰度化、中值滤波和自适应直方图均衡化等预处理；(2)根据辽宁省多年份实测路面影像制作大规模、多场景、像素级路面病害数据集，然后融合注意力机制及Dense Crf优化Res＿UNet网络结构完成模型训练；并引入损失函数dice loss增强了该方法对细小病害提取的能力；(3)将深度卷积神经网络分割后的路面病害特征图导入全连接条件随机场，对预测的路面病害结果进行轮廓优化，其检测结果为获取路面裂缝宽度，进而评估路面病害等级奠定了基础。该文选用2 000张辽宁省农村公路实测路面影像，并以人工判读作为标准，分别从准确率、召回率和精确率3个方面验证本文方法、分水岭算法和Res＿UNet模型在实际工作环境下的农村公路路面病害分割性能。结果表明，方法的准确率为91.3%,召回率为87.8%,精确率为87.5%,路面病害轮廓提取更加精细，能够适应于复杂路面条件下病害高鲁棒分割。
2022年09期 v.47;No.291 170-180页 [查看摘要][在线阅读][下载 697K] [下载次数：249 ] |[网刊下载次数：0 ] |[引用频次：2 ] |[阅读次数：1 ] - 赵莉;危小建;程朋根;谢名睿;王慧敏;
为了探究景观破碎化的驱动机制，该文运用景观格局指数、主成分分析、地理探测器和多尺度地理加权回归模型(MGWR)等方法，定量分析了大南昌都市圈景观破碎化分布特征，探究了驱动因子与景观破碎度间的空间相互作用。结果表明：研究区整体的景观破碎度在研究期内呈增加趋势，各县市和格网尺度上景观破碎度的空间分布差异较为显著；地理探测器结果显示，各因素均不同程度影响着研究区景观破碎化程度，且各驱动因素间的交互作用均有所增强；与地理加权回归模型(GWR)相比，MGWR模型的拟合效果显著提高，同时，由于空间位置的不同，影响研究区景观破碎化程度变化的主导因素也有所差异。研究结果可为区域土地利用合理规划和生态环境保护提供理论依据和决策参考。
2022年09期 v.47;No.291 181-191+223页 [查看摘要][在线阅读][下载 895K] [下载次数：1041 ] |[网刊下载次数：0 ] |[引用频次：8 ] |[阅读次数：1 ] - 任建设;王枭轩;张向军;禹小伟;余海坤;
针对变化检测方法对差异遥感影像去除斑点噪声和地物边缘检测效果不佳的问题，该文提出一种基于显著性卷积神经网络(SCWNN)的变化检测方法。该方法采用上下文感知显著性检测方法，提高遥感影像中差异较大区域显著性值，减小差异较小区域显著性值，有效优化了地物边缘，同时降低斑点噪声；然后采用卷积神经网络方法，克服了由于样本不足对变化检测精度的影响。通过超高分辨率影像和合成孔径雷达两组双时相遥感影像数据进行实验，结果表明方法的有效性和鲁棒性。
2022年09期 v.47;No.291 192-199+234页 [查看摘要][在线阅读][下载 628K] [下载次数：437 ] |[网刊下载次数：0 ] |[引用频次：4 ] |[阅读次数：0 ] - 侯赛男;党国锋;
针对城市热岛研究中MODIS地表温度数据云覆盖像元值问题，该文构建SVR、RF和XGBoost 3种机器学习模型以实现对MODIS地表温度产品缺失像元值的最优模拟，以最优模型获得郑州市2009—2018年无缝地表温度数据；基于距离加权不透水面聚集密度，提取出郑州市逐年主要建成区范围；结合地表温度数据、建成区范围及郊区范围动态分析郑州市热岛面积及热岛强度变化。结果表明：XGBoost模型回归拟合系数为0.95,均方根误差为0.06,填补MODIS地表温度数据像元值空缺效果优于SVR和RF模型；研究期间内，郑州市每年主要建成区面积增长速率均超过8%;郑州市这10 a间地表温度的标准差呈平稳下降的趋势，中温区面积呈上升趋势，热岛强度呈下降趋势，表明了郑州市“高温化”的缓解。
2022年09期 v.47;No.291 200-207页 [查看摘要][在线阅读][下载 414K] [下载次数：755 ] |[网刊下载次数：0 ] |[引用频次：8 ] |[阅读次数：1 ]