2024年4月27日发(作者:)

山东科学

第34卷第6期2021年12月出版

Vol.34No.6Dec.2021

SHANDONGSCIENCE

0.2m分辨率航空数字正射影像制作技术

曲莉莉ꎬ朱丰琪

(山东省国土测绘院ꎬ山东济南250013)

摘要:针对高分辨率数字正射影像制作中的数字高程模型编辑环节自动化程度低、编辑效果差引起的地物变形、移位及

匀光匀色环节中影像整体匀色不一致、接边痕迹明显、去雾气效果差的问题ꎬ从实际应用角度ꎬ分析了高分辨率数字正射

影像图制作的关键技术ꎮ利用滨州市域9600km

0.2m分辨率DMCIII航空影像数据ꎬ基于测区地理地貌特征ꎬ探讨了利

用Inpho、GEOWAYCIPS、PixelFactory等集群式影像处理系统制作1

:

2000数字正射影像图的技术难点及解决方法ꎮ对

高分辨率航空影像的处理提出建设性意见ꎬ有助于提高今后类似项目的质量与效率ꎬ为基础测绘助力国家建设奠定基

础ꎮ

关键词:数字正射影像图ꎻ0.2m分辨率ꎻ数字高程模型ꎻ制作技术ꎻ质量控制

中图分类号:P20   文献标志码:A   文章编号:1002 ̄4026(2021)06 ̄0127 ̄07

开放科学(资源服务)标志码(OSID):

Productiontechnologyofmappinganaerialdigitalorthophoto

witha0.2mresolution

QULi ̄liꎬZHUFeng ̄qi

(ShandongProvincialInstituteofLandSurveyingandMappingꎬJinan250013ꎬChina)

Abstract∶Focusingontheproductionmethodformappingahigh ̄resolutiondigitalorthophotoandtheproblemsthatoften

occurinthisproductionprocessꎬthekeytechnologyofproducinghigh ̄resolutiondigitalorthophotomapsfromthepractical

perspectiveisproposedandthetechnicalchallengesandcorrespondingsolutionsaresummarizedusingthe9600km

0.2mresolutionDMCⅢdigitalaerialimagedataofBinzhouareaandclusterimageprocessingsystems(suchasInphoꎬ

GEOWAYCIPSꎬandPixelFactory)toproduce1

:

2000digitalorthophotomaps.Constructivesuggestionsontheacquisition

andprocessingofhigh ̄resolutionaerialimagesarecarefullysuggested.Theproposalswillhelpinimprovingthequality

andefficiencyoffutureprojectsinthisfieldandwillcontributetothefundamentalsurveyingandmappingꎬwhichwillbenefit

theeconomicdevelopmentofChina.

Keywords∶digitalorthophotomapꎻ0.2mresolutionꎻdigitalelevationmodelꎻproductiontechnologyꎻqualitycontrol

收稿日期:2021 ̄01 ̄07

作者简介:曲莉莉(1963—)ꎬ女ꎬ高级工程师ꎬ研究方向为测绘产品质量验收ꎮTel:134****5352ꎬE ̄mail:2692006272@qq.com

通信作者ꎬ朱丰琪(1982—)ꎬ男ꎬ高级工程师ꎬ研究方向为航空摄影测量ꎮTel:158****2842ꎬE ̄mail:46594649@qq.com

Copyright©博看网 . All Rights Reserved.

  随着«测绘地理信息事业“十三五”规划»

[1]

的实施ꎬ高分辨率数字正射影像图(digitalorthophotomapꎬ

DOM)在基础测绘、国土调查、国情监测、城市规划、地籍测量等项目中发挥了重要作用ꎬ对高分辨率的研究

(digitalelevationmodelꎬDEM)编辑和DOMꎬ并提出了多种处理方式方法ꎮ渠甲源

[6]

曾就高分辨率航空正射

影像技术设计难点探讨了利用MicroStation软件分层编辑DEM数据ꎬ但编辑效率及效果不佳ꎬ自动化程度较

低ꎬ需使用大量人工编辑量ꎮ目前也无较为成熟的针对大面积、多时相、不同天气引起的航片色彩差异进行

整体均衡匀色的处理方法ꎮ随着高科技测绘技术的不断发展ꎬ传统处理模式已逐渐被先进的自动化的集群

式影像处理系统所替代ꎬ与传统处理模式相比ꎬ自动化集群式影像处理系统已有数倍的效率提升ꎬ但DEM编

辑及DOM匀色这两个环节的技术难点是目前数字正射影像快速生产急需解决的问题ꎬ将极大地影响自动

化处理系统的影像生产效率及产品质量ꎮ

PixelFactory集群式影像处理系统ꎬ采用自动化处理和人工干预相结合的方式生产DOMꎬ研究DEM编辑及

匀色工序中的关键技术ꎬ提升DEM自动化处理效率ꎬ降低建筑物变形、线状地物错位等错误ꎻ消除因天气、航

摄分区多、航摄时间跨度长、地表植被覆盖季节性变化、水体反光等原因造成的云雾遮盖、色调不一致、色彩

不真实的问题ꎬ以期突破影响影像自动化处理系统效率及质量的瓶颈ꎮ该项目对整个滨州辖区进行正射影

像处理ꎬ影像数据量及覆盖面积大ꎬ且辖区内有平原和丘陵地貌ꎬ北部临海ꎬ地物信息丰富ꎬ影像数据特征典

型ꎬ所采用的处理系统GEOWAYCIPS和PixelFactory为目前国内外最为先进的集群式影像处理系统ꎮ

针对以上技术难点ꎬ本文将以滨州市9600km

:

2000DOM制作项目为例ꎬ探讨利用GEOWAYCIPS、

显得十分必要ꎮ研究表明

[2 ̄5]

ꎬ在航空正射影像生产过程中ꎬ影响正射影像质量的重要因素是数字高程模型

1 主要技术指标及仪器设备

1.1 主要技术指标

带ꎬ中央经线117°ꎮ

的地面分辨率要求为0.2mꎻ彩色影像灰阶不应低于24bitꎬ灰度直方图应基本呈正态分布ꎻ采用50cm

×

50cm

正方形分幅ꎻDOM平面位置中误差平地、丘陵地不应大于1.2mꎬ山地、高山地不大于1.6mꎮ

无明显失真ꎮ

DOM成果应无明显拼接痕迹并保证建筑物等实体的影像完整ꎬ影像色彩应接近真实自然ꎬ纹理应清晰、

DOM成果由分幅DOM数据、元数据文件组成ꎮDOM数据格式为不压缩的TIFF格式ꎻ1

:

2000比例尺DOM

平面坐标系采用2000国家大地坐标系(CGCS2000)ꎬ投影方式为高斯

克吕格投影ꎬ分带方式为3°分

1.2 仪器设备

本文使用的主要仪器设备及软件系统主要有PC ̄6飞机、DMCIII航摄仪、Inpho、GEOWAYCIPS、Pixel

表1 主要设备与软件

Table1 Listofequipmentandsoftware

序号

名称

DMCIII航摄仪

GEOWAYCIPS

PixelFactory

Photoshop

Inpho

PC ̄6

用途

航摄平台

空中三角测量

航摄仪

Factory等ꎮ详见表1ꎮ

DEM、DOM数据生产、匀色

调色、DOM图面修复

匀色、镶嵌、分幅

2 技术路线及工艺流程

采用PC ̄6飞机搭载DMCIII数字航摄仪进行航空摄影ꎬ获取滨州摄区9600km

0.2m分辨率航空遥感

Copyright©博看网 . All Rights Reserved.

影像ꎻ利用GEOWAYCIPS、PixelFactory等集群式影像处理系统ꎬ经空三加密、DEM修测、正射校正、镶嵌、匀

光匀色等工序制作1

:

2000DOMꎮ工艺流程见图1ꎮ

图1 航空图像采集和处理流程图

Fig.1 Flowchartoftheacquisitionandprocessingofaerialimages

3 DEM优化研究

试验证明DEM的质量是决定数字正射影像质量的重要因素之一ꎬDEM优化质量的好坏直接影响DOM

精度ꎮ因此ꎬ在生产中对DEM编辑工序必须投入较大精力ꎬ在保证DOM的精度质量的前提下提高DEM编

辑效率ꎮ

3.1 DEM难点分析

滨州市地势南高北低ꎬ大致上由西南向东北倾斜ꎬ渐次过渡到大海ꎮ小清河以南的邹平南部属丘陵区ꎬ

地势高峻ꎬ其余均为山前倾斜平原ꎬ地势平缓ꎬ海拔高程一般在8~800mꎮ小清河以北为黄河冲积平原ꎬ海拔

高程一般在1~20mꎬ总体上地势低平ꎬ呈现出滨州地区整体地貌较平坦但局部地势高峻的地貌特征ꎮ因此ꎬ

滨州项目存在以下难点问题ꎬ影响DEM的精度ꎮ

DSM密集匹配时高层建筑物会存在匹配漏洞ꎬ导致此类建筑物在点云形态与建筑物形态匹配度不吻合ꎬ在

对建筑物滤波时会存在漏滤现象ꎮ

滨州市周边存在大量的工业园区ꎬ部分工业园区厂房较大ꎬ厂房内露天设施外型复杂ꎮ对于人工构筑物

较大及造型较为复杂等情况ꎬ点云滤波参数不宜设置过大ꎬ防止滤波时地貌特征损失过多ꎮ

征采用同一过滤参数效果不佳ꎮ

(2)滨州市整体地势较为平坦ꎬ基本以平原为主ꎬ但最南边邹平市存在丘陵地貌ꎮ对于此类多种地貌特

滨州市整体以平原为主ꎬ航摄时间为2020年4月—6月ꎬ测区内存在大面积的小麦种植区ꎮ测区北部

沿海存在大量盐田、海水养殖场等ꎬ大片纹理类似地物在DSM匹配时会存在大量的错误匹配点ꎬ这些错误点

在滤波时影响滤波效果ꎮ

(1)滨州城区内高层建筑物多ꎬ航摄影像在航向重叠约60%左右ꎬ旁向重叠约30%左右ꎮ此类重叠度在

Copyright©博看网 . All Rights Reserved.

其余非地面点需人工过滤ꎮ

3.2 DEM生产流程

(3)集群系统在密集匹配时可能会产生噪声高点或低点ꎬ而自动滤波只能有效过滤约80%的非地面点ꎬ

对原始航空影像数据进行空三加密ꎬ然后将空三加密成果导入到GEOWAYCIPS、PixelFactory集群式

影像处理系统中ꎬ在集群中采用逐像素密集匹配同名点的方式进行数字地表模型(digitalsurfacemodelꎬ

DSM)提取ꎬ然后经DSM滤波、DSM矢量化、DEM编辑等工序ꎬ制作数字高程模型DEMꎮ

DSM进行植被及建筑物滤波ꎬ去除植被及房屋等处的点云ꎮ

后续编辑ꎮ

(1)DSM匹配及滤波:利用空三成果ꎬ生成立体模型ꎬ对立体模型进行核线影像密集匹配生成DSMꎬ对

(2)DSM矢量化:经过滤波的DSM进行构TIN重采样ꎬ生成DEM过程数据ꎬ再生成单片正射影像ꎬ用于

(3)DEM编辑:对单片正射影像进行变形检查ꎬ对相应位置的DEM数据进行编辑处理ꎮ除去滤波过程

未过滤掉的高程异常值ꎮ

3.3 试验结果分析

在滤波前根据地貌地物特征ꎬ将不同地貌及地物区域分开ꎮ如将平地与丘陵地区分开ꎬ各自调整过滤参

数进行过滤ꎮ滤波参数将最大建筑物长度、最大地表高差、迭代距离、是否进行山脊填补设为关键参数值ꎮ

其中平地过滤参数中的建筑物长度可适当放宽ꎬ默认参数为60mꎬ可根据测区实际建筑长度及宽度调整至

120mꎬ丘陵地区建筑物长度可适当调小至40mꎬ即该测区如果平地有建筑物长度超过120mꎬ丘陵地区有超

过40m的建筑ꎬ将不会进行过滤ꎻ最大地表高差默认参数为1.2mꎬ平地调整至0.6mꎬ丘陵调整至1.4mꎬ地

表高差是判断地形起伏程度的参数ꎬ平地变化小ꎬ丘陵地变化大ꎻ迭代距离默认参数为1.8mꎬ平地可调整至

0.8mꎬ丘陵调整至1.2mꎬ该参数主要是配合地表高差参数ꎬ形成滤波窗口尺寸ꎻ滤波时有时会将比较平坦的

山脊过滤掉ꎮ平地无山脊ꎬ不进行山脊填补ꎮ将滨州市邹平区单独分开ꎬ采用丘陵地区的滤波参数进行过

滤ꎬ其余区域采用平地滤波参数过滤ꎮ同理将城区有大量高层建筑的区域与周边乡镇及农村区域分开滤波ꎮ

滤波参数过滤效果成功率约在80%左右ꎬ滤波后仍有较多的高程异常值需要处理ꎮ针对这种情况ꎬ为

减小人工干预工作量ꎬ将平地区域与丘陵区域分开ꎮ平地经滤波后ꎬ采用比较大的平滑窗口ꎬ对整体地形进

行平滑处理ꎬ可获得较好的平滑结果ꎬ同时能避免一些地物轻微变形ꎬ见图2ꎮ丘陵地区采用较小的平滑窗

口ꎬ以免平滑损失过多的地形信息ꎬ从而影响后续正射影像纠正精度ꎮ针对不同地貌特征采用合理的平滑优

化ꎬ可极大程度减少异常高程值引起的地物变形ꎬ节省大量人工编辑工作量ꎮ

图2 建筑物变形修正

Fig.2 Buildingdeformationcorrection

全自动化处理时不可避免会存留部分粗差ꎬ采用人机交互的方式对DEM过程数据进行编辑ꎬ解决建筑

物扭曲变形、地形地貌失真等问题ꎮ设置平地参数处理后的地区有可能造成山头缺失ꎬ需要人工添加匹配

点ꎻ大面积水域可手工沿水域的边缘选择编辑区域ꎬ输入水面高程ꎬ用平面拟合算法进行平滑处理ꎻ将异常错

误的点云手动人工删除ꎬ手动添加缺失的特征点ꎬ使高程网准确贴合地面ꎮ

Copyright©博看网 . All Rights Reserved.

4 DOM匀色优化研究

4.1 DOM难点问题分析

滨州测区航摄分区6个ꎬ原始航摄数据量大小为3.4TBꎮ航摄时间跨度为2020年4月中旬至6月底ꎮ

测区南北狭长ꎬ整体为平原地貌ꎬ且北临渤海湾ꎮ滨州地区4月至5月份以晴朗天气为主ꎬ6月份天气逐渐

向多云多雾变化ꎮ以上情况会导致下列问题ꎮ

镶嵌时不同架次接边处相同地物、植被存在不同的色彩及纹理特征ꎬ影响影像接边质量及效果ꎮ

的曝光现象ꎬ影响水体图面效果ꎮ

4.2 DOM制作流程

DOMꎮ

(1)因航摄时间跨度较长ꎬ地表植被覆盖生长情况差异明显ꎬ造成各架次影像色调严重不一致ꎬ在影像

(2)测区北部临渤海海域ꎬ存在大面积盐池及海水养殖场ꎬ水面整体反光严重ꎬ局部存在因强反光产生

(3)部分架次影像存在轻度雾气ꎬ整体视觉通透性差ꎬ影像噪点多ꎬ地物纹理表现力差ꎮ

基于制作完成的DEM成果和空三加密成果ꎬ进行单片正射纠正、匀光匀色、镶嵌、图面修复等工序ꎬ制作

(1)单片纠正:利用DEM成果和空三加密成果对滨州原始航摄影像进行逐片微分正射校正ꎮ

正后的单片正射影像采用匀色模板进行匀光匀色处理ꎬ使所有单片影像颜色趋于一致ꎮ

(2)匀光匀色:调整滨州测区较为合适的匀色模板ꎬ对因跨时相、跨架次、天气等引起的颜色差异基于纠

(3)影像镶嵌:对匀色后的正射影像自动化镶嵌处理ꎬ对镶嵌线切割各类线状地物、建筑物造成明显错

位等问题进行人工镶嵌线编辑ꎬ使影像各要素完整、拼接痕迹不明显、纹理清晰、城区高层建筑物投影方向一

致ꎬ避免建(构)筑物相互压盖、交错、反向等现象ꎮ

各类图面问题进行修复ꎮ

4.3 试验结果分析

(4)图面修复:对标准分幅影像成果进行人工检查ꎬ对因DEM编辑不当造成的影像变形问题以及其他

首先利用GEOWAYCIPS对测区所有单片先进行第一次匀色处理ꎬ使每张影像的色彩色调趋于一致ꎮ

然后将匀色后的单片导入到PixelFactory中进行镶嵌和二次匀色ꎬ提高了影像成果的质量和显示效果ꎬ方法

如下ꎮ

据选出几张典型影像ꎬ微调匀色模板测试匀色效果直至影像地物纹理清晰、层次分明、色彩均衡一致、直方图

基本呈正态分布为止ꎬ然后用调试好的匀色模板对所有轻度雾气影像进行匀色ꎻ大面积水域反光或有曝光现

象的影像降低匀色模板亮度后进行匀色ꎮ

羽化参数ꎬ过渡镶嵌效果ꎬ此步骤可进一步优化基于第一步匀色后的影像匀色效果ꎻ对于大面积盐池和水域

反光ꎬ尽量将镶嵌线从盐池或养殖场的堤岸处过渡ꎬ避免镶嵌线切割水面ꎬ造成阴阳水面ꎮ

像处理软件对其调整色彩色调ꎬ然后将其重新导入至PixelFactory中ꎬ将其颜色以地理映射的方式应用于全

局ꎬ以小带大的方式获得整个测区最佳匀色效果ꎬ见图3ꎮ

(3)镶嵌线编辑完毕后将整个滨州测区导出一个低分辨率、数据量小的镶嵌影像模板ꎬ在Photoshop图

(2)在PixelFactory中编辑镶嵌线ꎬ使其沿块状或线状地物的边沿ꎬ避免切割同一地物ꎬ同时设置合理的

(1)合理选择匀色模板ꎮ时相或影像质量一致或接近的为一组ꎬ采用同一匀色模板ꎻ存在轻度雾气的数

Copyright©博看网 . All Rights Reserved.

图3 测区匀色前和匀色后

  镶嵌线沿街道、公路、河流线状地物进行编辑ꎬ尽量避免镶嵌线从线状地物中间穿过分割建筑物、山脊

等ꎬ使影像镶嵌后色彩自然合理、地物完整ꎮ

Fig.3 Beforeandafterimageuniformityoftestingzone

4 结论

本文采用较先进的技术方案ꎬ利用集群式影像处理系统ꎬ针对正射影像生产流程中对精度及效果影响较

大的DEM编辑及影像匀光匀色关键性技术进行研究ꎬ得出了如下结论:

DEM质量的同时最大化自动化处理ꎬ减少了大量人工干预环节ꎬ节省人力成本和时间ꎮ处理后的DEM成果

平面精度及高程精度满足后续DOM生产要求ꎬ即无明显高程异常值ꎬ地形表述与实际地形一致ꎬ为后续生

产提供了质量较优的DEM数据ꎮ

地物合理编辑镶嵌线及羽化效果ꎬ采用地理映射的方式对整测区进行二次匀色ꎮ匀色后的影像地物细节清

晰ꎬ反差适中ꎬ层次分明ꎬ色彩均衡ꎬ镶嵌线两侧和相邻影像色调应基本一致ꎮ

通过以上试验方法制作完成DOMꎬ平面精度及色彩纹理质量均达到了预期的0.2m分辨率的标准要

求ꎮ试验方法合理有效ꎬ基于目前市面最先进的集群式影像处理系统ꎬ进一步提升生产了DEM编辑自动化

程度、减少了人工干预ꎻ匀光匀色消除因原始影像引起的质量问题ꎬ形成一套快速规模化数字正射影像处理

方案ꎬ能有效指导后续此类项目的生产工作ꎮ

(2)DOM处理中的影像匀光匀色环节按照影像质量合理使用不同的匀色模板进行初次匀色ꎬ根据测区

(1)DEM编辑处理中采用按地貌类别分区滤波、合理使用平滑优化参数、人机交互编辑等方法ꎬ在保证

参考文献:

[1]发展改革委ꎬ测绘地信局.关于印发«测绘地理信息事业“十三五”规划»的通知[EB/OL].[2020–08–31].http://www.

gov.cn/gongbao/content/2017/content_5194895.htm.

DevelopmentandReformCommissionꎬBureauofSurveyingꎬMappingandInformationTechnology.Noticeonprintingand

distributingthe13thfiveyearplanforsurveyingꎬmappingandgeographicinformationindustry[EB/OL].[2020–08–31].

[2]刘学杰.像控布设方案对无人机航测精度影响的测试[J].地理信息世界ꎬ2016ꎬ23(5):109 ̄112.DOI:10.3969/j.issn.

1672 ̄1586.2016.05.026.

LIUXJ.Thetestontheeffectofcontrolpoints′layoutonaerialphotogrammetryaccuracyofUAV[J].GeomaticsWorldꎬ2016ꎬ

http://www.gov.cn/gongbao/content/2017/content_5194895.htm.

23(5):109 ̄112.DOI:10.3969/j.issn.1672 ̄1586.2016.05.026.

Copyright©博看网 . All Rights Reserved.

[3]黄冠伟.航空摄影测量野外像控点布设方案与测量方法[J].工程建设与设计ꎬ2017(10):209 ̄210.DOI:10.13616/j.cnki.

gcjsysj.2017.05.197.

HUANGGW.Layoutschemeandmeasurementmethodofaerialphotographicsurveyingfield[J].Construction&Designfor

[4]廖瑞钢.浅析航空摄影测量像控点的布设与测量[J].绿色环保建材ꎬ2017(7):99.DOI:10.16767/j.cnki.10 ̄1213/tu.2017.

07.084.

LIAORG.Analysisonthelayoutandmeasurementofimagecontrolpointsinaerialphotogrammetry[J].GreenEnvironmental

ProtectionBuildingMaterialsꎬ2017(7):99.DOI:10.16767/j.cnki.10 ̄1213/tu.2017.07.084.

gnss.2012.04.017.

Engineeringꎬ2017(10):209 ̄210.DOI:10.13616/j.cnki.gcjsysj.2017.05.197.

[5]张汉德ꎬ张海涛ꎬ别君ꎬ等.影响正射影像生成质量的因素分析[J].全球定位系统ꎬ2012ꎬ37(4):60 ̄63.DOI:10.13442/j.

ZHANGHDꎬZHANGHTꎬBIEJꎬetal.AnalysisofthefactorsaffectingthequalityofDOM[J].GNSSWorldofChinaꎬ2012ꎬ

[6]渠甲源.高分辨率航空正射影像技术设计难点探讨[J].测绘与勘探ꎬ2019ꎬ1(1):65 ̄66.DOI:10.36012/se.v1i1.343.

2019ꎬ1(1):65 ̄66.DOI:10.36012/se.v1i1.343.

37(4):60 ̄63.DOI:10.13442/j.gnss.2012.04.017.

QUJY.Discussiononthedesigndifficultiesofhighresolutionaeronauticalorthoimagetechnology[J].SurveyingandExplorationꎬ

[7]陈杏.基于INPHO和像素工厂的真正射影像制作与对比分析[D].湘潭:湘潭大学ꎬ2017.

XiangtanUniversityꎬ2017.

1007/tk.2018.09.033.

CHENX.ProductionandcomparativeanalysisofthetrueorthophotomapbasedonINPHOandpixelfactory[D].Xiangtan:

[8]徐小芹ꎬ张田凤ꎬ许明.航空测量中像控点布设问题探究[J].低碳世界ꎬ2018(9):53 ̄54.DOI:10.16844/j.cnki.cn10 ̄

XUXQꎬZHANGTFꎬXUM.ResearchontheLayoutofimagecontrolpointsinaerialsurveying[J].LowCarbonWorldꎬ2018

[9]周军ꎬ龙盈ꎬ王发艳.低空数码航空摄影测量像控点的布设及构网[J].地矿测绘ꎬ2015ꎬ31(2):35 ̄36.DOI:10.16864/j.

cnki.dkch.2015.02.010.

ZHOUJꎬLONGYꎬWANGFY.Layoutofimagecontrolpointsinlowaltitudedigitalaerialphotography[J].Surveyingand

[10]宋焕胜.DOM、DEM在城市规划设计和管理中的应用[J].智能建筑与智慧城市ꎬ2018(6):33 ̄34.DOI:10.13655/j.cnki.

ibci.2018.06.011.

SONGHS.TheapplicationofDOMandDEMinurbanplanninganddesignandmanagemente[J].IntelligentBuilding&Smart

Cityꎬ2018(6):33 ̄34.DOI:10.13655/j.cnki.ibci.2018.06.011.

MappingofGeologyandMineralResourcesꎬ2015ꎬ31(2):35 ̄36.DOI:10.16864/j.cnki.dkch.2015.02.010.

(9):53 ̄54.DOI:10.16844/j.cnki.cn10 ̄1007/tk.2018.09.033.

Copyright©博看网 . All Rights Reserved.