Methods of ASTER remote sensing data used in extracting mineralization and alteration information in the Hutouya mining area of the Qimantag metallogenic belt, Qinghai Province
-
摘要: 以ASTER遥感影像为数据源，对青海祁漫塔格成矿带虎头崖矿区，应用成像光谱法和主成分分析法对矿区矿化蚀变异常信息进行提取分析。成像光谱方法利用最小噪音分量（MNF）变换、像元纯度指数（PPI），n维可视化（n-Dimensional Visualization）端元识别，并进行混合谐调匹配滤波（MTMF），最终得到矿区的矿物信息分布图。而主成分分析法则选择ASTER1、2、3、4波段提取铁染异常信息，1、3、4、5波段提取矽卡岩化蚀变异常信息，1、3、4、（5+6）/2波段提取蒙脱石、伊利石与绢云母等矿物的蚀变异常信息，1、3、4、8波段提取方解石、黑云母、绿泥石等矿物的蚀变异常信息。利用野外实地踏勘确定的矿床、矿（化）点，在ArcGIS10.1软件中对遥感提取蚀变异常信息进行验证显示，二者吻合效果良好，表明提取结果与数据处理方法较可靠，对于扩大找矿规模和发现矿化富集地段具有重要的指导意义。Abstract: Using ASTER remote sensing images as the data source, the extraction of mineralization and alteration anomaly in the Hutouya mining area of the Qimantag metallogenic belt in Qinghai Province is studied through two methods. One method is using the minimum noise component (MNF) transform, pixel purity index (PPI), n-Dimensional Visualization to do the end element recognition and mixture tuned matched filtering (MTMF), in order to obtain the mineral distribution map of the mining area. The other is principal component analysis, through which ferric contamination anomaly is extracted in ASTER1, 2, 3, 4 band, skarnization alteration anomaly in 1, 3, 4, 5 band, montmorillonite, illite and sericite in 3, 4, (5+6)/2 band and calcite, biomica, chlorite and other minerals in 1, 3, 4, 8 band. The results of remote sensing extraction of alteration anomaly verified by ArcGIS10.1 software are in good consistency with the mineral deposits and mineralization points determined by the field survey, indicating that the extraction results and the data processing method are reliable, which is of great guiding significance for the expansion of the prospecting scale and the discovery of mineralization enrichment areas.
-
虎头崖多金属矿区位于柴达木盆地西南缘，东昆仑祁漫塔格山北坡，隶属于东昆仑祁漫塔格成矿带。矿区自然环境差，海拔高，地质构造复杂，是该区比较典型的矽卡岩型多金属矿床(李智明等，2007；张爱奎等，2010；高晓峰等，2010；胡杏花等，2011；伍跃中等，2011；马圣钞等，2013；刘渭等，2014；易立文等，2019；张大明等，2020)，矿区侵入岩分布较广，岩体与地层的接触带附近围岩蚀变普遍强烈，其蚀变类型主要有绿帘石化、绿泥石化、黄铁矿化、透闪石化、阳起石化、石榴石化、透辉石化、绢云母化和大理岩化等(张爱奎等，2010；高晓锋等，2010；胡杏花等，2011；刘渭等，2014；杨兴科等，2016)。一些学者已经广泛开展了采用ASTER数据进行遥感蚀变异常信息提取的研究(张玉君等，2006；毛晓长等，2005；耿新霞等，2008；杨长保等，2009；张玉君和姚佛军，2009；高万里等，2010；燕守勋等，2011；姚佛军等，2012；陈菁和周萍，2015；周迪等，2018)，然而对于单个矿区虎头崖矿区遥感找矿方法和手段报道较少。因此现以ASTER遥感影像为数据源，应用成像光谱法和主成分分析法蚀变遥感异常提取技术方法对青海祁漫塔格成矿带虎头崖矿区进行蚀变异常信息提取，并对部分异常点进行了实地查证，以期为扩大找矿规模和发现矿化富集地段提供重要的指导意义。
1. 研究区成矿地质特征
虎头崖多金属矿床位于青海省祁漫塔格地区，大地构造位置上处于东昆仑西段，北昆仑(祁漫塔格)岩浆弧带(李荣社等，2008), 是该区比较典型的矽卡岩型多金属矿床，主要有Fe、Cu、Pb、Zn、Wu和Sn等矿种(张晓飞等2012；刘渭等，2014)。
矿区地层属柴达木地层区的柴达木南缘分区，主要为蓟县系狼牙山组(Jxl)，奥陶系—志留系滩涧山群(OST)，下石炭统大干沟组(C1dg)，上石炭统缔敖苏组(C2d)，中生界上三叠统鄂拉山组(T3e)以及山麓地带及河床等广泛发育的第四系(Q)松散堆积物。矿区内褶皱、断裂等构造非常发育，构造形迹主要表现为东—西向断裂构造(刘渭等，2014)。矿区侵入岩主要出露于矿区中部，以中酸性侵入岩类为主，侵入于狼牙山组(Jxl)、滩涧山群(OST)、大干沟组(C1dg)及缔敖苏组(C2d)等地层中(李侃等，2015)(图 1)。狼牙山组上岩段(碳酸盐段Jxlb)、滩间山群火山岩(OST)、下石炭统大干沟组(C1dg)和上石炭统缔敖苏组(C2d)，大理岩、含铁石英砂岩、灰岩、硅质岩岩性段是主要的含矿岩性段(苏松，2011)。矿区内典型矿床类型为矽卡岩型，褶皱、断裂和岩体与地层接触带是矿体赋存的有利部位。蚀变类型主要有透闪石化、透辉石化、阳起石化、绿泥石化、绿帘石化、石榴子石化、高岭土化、硅化、绢云母化等(杨兴科等，2016)。
图 1 虎头崖矿区地质简图(刘渭等，2014)1—第四系；2—蓟县系狼牙山组；3—奥陶志留系滩间山群；4—下石炭统大干沟组；5—上石炭统缔敖苏组；6—上三叠统鄂拉山组；7—花岗闪长岩；8—二长花岗岩；9—石英斑岩；10—似斑状二长花岗岩；11—闪长岩；12—闪长玢岩脉；13—推断断层；14—逆冲断层；15—褶皱；16—矿带；17—青海祁漫塔格成矿带及虎头崖矿区Figure 1. Simplified geological map of the Hutouya mining area(Liu et al, 2014)2. 遥感数据介绍及预处理
2.1 ASTER数据介绍
ASTER数据有14个波段：可见光、近红外和短波红外到热红外，光谱分辨率较高，特别是在短波红外有6个波段、热红外有5个波段，它们分别对粘土矿物、碳酸盐和硅酸盐类矿物具有较好的识别能力(杨日红等，2012)。
ASTER数据主要有标准HDF格式数据和dat格式地面站数据两种。文中所用ASTER地面站数据采集日期为2003年9月23日，数据级别为L1A，为ASTL1A 0309230444060310050159B，包括了5个文件：prbr0180.brs、prbr0181.jpg、prbr0182. jpg、0183.jpg、prdat018.dat。其中prbr0181.jpg、prbr0182.jpg、prbr0183.jpg 3个文件为多光谱影像的示意图，prbr0180.brs为数据头文件，prdat018.dat为数据源文件。
2.2 数据预处理
2.2.1 串扰校正与辐射定标
因ASTER短波红外通道传感器的相互干扰，B5对B4和B9具有串扰现象，使得B4和B9反射率值升高，因此须对数据进行串扰校正(杨日红等，2012；高慧等，2013)。同时ASTER数据在VNIR波段(1—3)空间分辨率为15 m，在SWIR波段(4—9)空间分辨率为30 m，需对4—9波段进行重采样将分辨率变为15 m，再利用ASTER Radiance模块实现1—9波段图像的辐射定标(刘知和周萍，2017)。
2.2.2 大气校正
将上述生成的1—9波段Radiance图像转换生成BIL格式，之后进行FLAASH大气校正。具体过程包括，一是生成的Radiance图像头文件会丢失中心波长半极值宽度(FWHM)，需要打开头文件依次添加1—9波段头文件相关信息；二是对于ASTER数据FLAASH大气校正来说，其所需的量纲与辐射半径定标后的数据量纲不同，这里辐射亮度值的比例因子应设置为10，完成量纲转换(陈建明等，2009)。三是通过对Radiance图像属性读取功能，设置全部的参数，从而完成FLAASH大气校正。对校正后的影像裁剪，得到稍大于研究区范围的预处理数据。
3. 成像光谱方法提取矿化信息
光谱变换技术对于减少冗余信息十分必要，光谱数据可以通过最小噪音分量(MNF)变换进行降维和压缩。选取预处理后ASTER的6个短波红外波段进行矿物信息提取，基于最小噪音分量MNF变换，像元纯度指数(PPI，Pixel Purity Index)进行波谱端元选择，采用混合谐调匹配滤波(MTMF，mixture tuned matched filtering)等方法进行矿物填图(Pour and Hashim, 2012)。
3.1 最小噪音分量(MNF)变换
最小噪音分量(MNF)变换是一种光谱数据减维技术，可以分离数据中的噪声，减小进一步处理所需的运算量(高慧等，2013)。通过主组分变换。将数据分为两类：一类为噪声图像，另一类为特征图像，与大特征值有关。通过观察最终特征值和MNF图像(特征图像)来确定数据的固有维数，即只利用其中的相关部分将数据中的噪声分离，改善光谱处理的结果。根据预处理后ASTER的6个短波红外波段完成MNF变换，得到特征值(表 1)和较大特征值MNF波段1—3灰度图像(图 2—图 4)。可见MNF波段特征值越大，图像所含有用的信息就越多。
表 1 各波段平均值和特征值Table 1. Average value and eigenvalue of each band波段名称 平均值 特征值 Band 1 21.173028 28.297260 Band 2 7.034879 9.762023 Band 3 6.628618 4.759365 Band 4 6.013828 2.901527 Band 5 4.625094 2.367740 Band 6 3.821772 1.830600 
图 2 MNF变换后波段1灰度图像Figure 2. Grayscale image of the band 1 after the MNF transformation
图 3 MNF变换后波段2灰度图像Figure 3. Grayscale image of the band 2 after the MNF transformation
图 4 MNF变换后波段3灰度图像Figure 4. Grayscale image of the band 3 after the MNF transformation3.2 像元纯度指数(PPI)
像元纯度指数(PPI)是在多光谱或成像光谱图像中寻找“光谱上最纯的”或极端像元的一种工具，它用于从混合像元中提取较纯像元，从而减少确定端元所要分析的像元数目，容易分离和识别端元。基于上述MNF变换后的图像，忽略特征值低的MNF波段，进一步处理高特征值波段(高慧等，2013)。通过多次反复投影N维散点图到随机单元矢量，记录每次投影的极值像元，最后注记每个像元被标定为极值的总次数，最终产生PPI图像(图 5)。
3.3 n维可视化(n-Dimensional Visualization)进行端元识别
n维可视化(n-Dimensional Visualization)，光谱可视为n维散点图中的一个点，n是波段数。对给定的像元，在n维空间中，每个波段对应像元的光谱反射率组成各个点位的坐标，利用点位在n维空间中的位置可以估计分析光谱端元数和它们的纯光谱特征。混合的端元落在纯端元之间，处在纯端元勾画的多面体内，这种混合光谱的凸面几何特征，可用于确定端元光谱数并估算它们的光谱特征，经过PPI选择出的潜在端元光谱输入n维散点图中进行反复多次旋转可识别出纯端元(Pour and Hashim, 2012)。
这里采用了4维可视化散点图，并勾画较好的端元(图 6)。根据前面的分析，较好的端元通常会出现在n维散点图的顶点和拐角处，当一系列的端元点被确定后，就可以将其输入到图像中的感兴趣区(ROI)，从图像中提取每个感兴趣区平均反射率光谱曲线(图 7)作为成像光谱矿物填图的端元。
图 7 n-D Mean绘图窗口获取的平均波谱曲线Figure 7. Average spectral curve obtained by the n-D Mean drawing window3.4 混合谐调匹配滤波(MTMF)
混合谐调匹配滤波(MTMF)是线性混合理论与匹配滤波技术相结合的一种方法，融合了线性光谱分解和匹配滤波的技术优势，能探测微细矿物成分。利用最终确定的端元光谱(它使用端元的MNF光谱而非反射率光谱)进行混合谐调匹配滤波(MTMF)的结果值为0—1范围的灰度图像，可以对矿物的相对丰度进行估计(亮度值高的像素代表较高的矿物含量)(Pour and Hashim, 2012)。最终得到矿区的矿物信息分布图(图 8)。
图 8 虎头崖矿区的矿化信息提取图粉红色—蒙脱石；红色—绿泥石；黄色—方解石Figure 8. Extraction map of mineralization information in the Hutouya mining area4. ASTER主成分分析方法提取矿化蚀变信息
基于研究区围岩蚀变发育情况和主要的蚀变类型(苏松，2011；杨兴科等，2016)，从USGS标准波谱数据库中，有针对性地选择了透闪石、方解石、绿泥石、伊利石、蒙脱石；铁的蚀变矿物代表黄铁矿等的连续光谱，将其重采样到研究区ASTER数据光谱，从而得到研究区主要的蚀变矿物的光谱(图 9)。根据研究区主要蚀变矿物的波谱特征(图 9)显示，伊利石、蒙脱石等在6波段形成一个明显的吸收谷，该吸收谷与Al-OH在2.2 μm处的吸收谷一致；透闪石矿物在5波段有一个微小的吸收谷；绿泥石在3波段有一个微小的吸收谷，方解石也在8波段有一个明显的吸收谷；黄铁矿为代表蚀变矿物在1波段处有吸收谷，在4波段中也存在一个明显的吸收谷。
图 9 研究区ASTER数据主要蚀变矿物的光谱Figure 9. Spectra of major altered minerals in ASTER data of the study area通过对蚀变矿物及所对应的离子或离子团特征分析，基于ASTER数据采用ASTER1、ASTER3、ASTER4、ASTER (N)的波段组合进行主成分分析(耿新霞等，2008；陈晔等，2014；辜平阳等，2016)，分别选择1、2、3、4波段提取铁染信息；1、3、4、5波段提取矽卡岩化蚀变信息；1、3、4、(5+6)/2波段提取蒙脱石、伊利石等Al-OH离子团蚀变矿物的异常信息；利用1、3、4、8波段提取方解石、黑云母、绿泥石等Mg-OH离子团蚀变矿物的异常信息。
4.1 含羟基矿化信息提取
利用ASTER1、ASTER3、ASTER4、ASTER(5+6)/2的波段组合进行主成分分析(张玉君等，2006)，提取包括蒙脱石、伊利石与绢云母等Al-OH离子团蚀变矿物的遥感异常信息。根据Al-OH离子团蚀变矿物主分量的判断标准及对参数值PC1-PC4(表 2)分析可知，ASTER4贡献系数符号为负，ASTER3与ASTER(5+6)/2贡献系数与ASTER4符号相反，PC4为含有蚀变矿物信息的主分向量，聚集了蚀变遥感异常信息。而表 2所示的PC4中ASTER4特征向量的符号为负，故需作PC4×(-1)的处理，使得ASTER4的特征向量符号为正。对虎头崖矿区PC4图像进行低值切割，且进行了异常滤波，即以2倍σ(标准离差)作为主分量输出的动态范围，获得虎头崖地区蚀变遥感异常信息(图 10)。
表 2 虎头崖矿区ASTER1、3、4、(5+6)/2主分量变换统计值Table 2. Statistical values of ASTER1, 3, 4, (5 + 6)/2 principal component transformation in the Hutouya mining areaASTER1 ASTER3 ASTER4 ASTER(5+6)/2 PC1 0.78201 0.61025 0.12076 0.03833 PC2 -0.62199 0.75182 0.21406 0.04560 PC3 -0.03748 0.24894 -0.91790 -0.30674 PC4 -0.01380 0.01967 -0.31155 0.94993 
图 10 虎头崖矿区ASTER遥感图像信息提取总异常图深紫色—Al-OH离子团矿化蚀变信息；黄色—Mg-OH离子团矿化蚀变信息；红色—铁染矿化蚀变信息；蓝色—矽卡岩化蚀变矿物信息Figure 10. Total anomalies of the ASTER remote sensing image extraction in the Hutouya mining area采用ASTER1、ASTER3、ASTER4、ASTER8波段参与主成分分析，提取Mg-OH离子团蚀变矿物的遥感异常信息。根据Mg-OH蚀变矿物主分量的判断标准和信息提取主成分的物理参数可知，ASTER4贡献系数符号为负，ASTER3与ASTER8贡献系数与ASTER4符号相反，PC4为含有蚀变矿物信息的主分向量，聚集了蚀变矿物信息。而表 3所示的PC4中ASTER4特征向量的符号为负，故需作PC4×(-1)的处理，使得ASTER4的特征向量符号为正。对虎头崖矿区第四主分量图像进行低值切割，且进行了异常滤波来获得Mg-OH离子团蚀变矿物的异常信息图像。即以2倍σ(标准离差)作为主分量输出的动态范围，获得虎头崖矿区蚀变遥感异常信息(图 10)。
表 3 虎头崖矿区ASTER1、3、4、8主分量变换统计值Table 3. Statistical values of ASTER1, 3, 4, 8 principal component transformation in the Hutouya mining areaASTER1 ASTER3 ASTER4 ASTER8 PC1 0.78225 0.61043 0.12078 0.02951 PC2 -0.62173 0.75383 0.21140 0.02240 PC3 -0.03460 0.24181 -0.93799 -0.24599 PC4 -0.01824 0.02538 -0.24679 0.96856. 4.2 铁染矿化蚀变矿物的蚀变遥感异常信息提取
采用ASTER1、2、3、4波段做主成分分析，提取与含Fe3+矿物相关的铁染异常主分量，主成分变换的统计结果如表 4所示，在PC4图像中包含了较多的铁染矿化蚀变信息，并且表现为亮色调。
表 4 虎头崖矿区ASTER1、2、3、4主分量变换统计值Table 4. Statistical values of ASTER1, 2, 3, 4 principal component transformation in the Hutouya mining areaASTER1 ASTER2 ASTER3 ASTER4 PC1 0.61068 0.62611 0.47567 0.09382 PC2 0.53342 0.12541 -0.80518 -0.22677 PC3 0.45937 -0.56826 0.02367 0.68227 PC4 0.36264 -0.51898 0.35338 -0.68868 对虎头崖矿区第四主分量图像，进行高值切割获得铁染异常信息图像。以2倍σ(标准离差)作为主分量输出的动态范围，获得虎头崖地区蚀变遥感异常主分量(图 10)。
4.3 矽卡岩化蚀变矿物的蚀变遥感异常信息提取
采用图像ASTER1、3、4、5波段做主成分分析，提取与含矽卡岩化矿物相关的蚀变遥感异常信息，主成分变换的统计结果如表 5所示，在PC4图像中包含了一些的矽卡岩化信息，并且表现为亮色调。
表 5 虎头崖矿区ASTER1、3、4、5主分量变换统计值Table 5. Statistical values of ASTER1, 3, 4, 5 principal component transformation in the Hutouya mining areaASTER1 ASTER3 ASTER4 ASTER5 PC1 0.78211 0.61033 0.12077 0.03498 PC2 -0.62190 0.75288 0.21263 0.03464 PC3 -0.03584 0.24549 -0.92974 -0.27209 PC4 -0.01620 0.02015 -0.27530 0.96101 对虎头崖矿区第四主分量图像，进行高值切割获得矽卡岩化异常信息图像。以2倍σ(标准离差)作为主分量输出的动态范围，获得虎头崖矿区蚀变遥感异常主分量(图 10)。
4.4 异常的中值滤波处理
蚀变提取后图像上会出现比较均匀的某一类图斑上分布孤立异常点的现象，可以去掉图中过于零星孤立的异常像元，或将孤立像元归并到包围它们或与它们相邻的较连续分布的那些异常中去(张玉君等，2006)。这里需要采用滤波的方法，对于每一类异常，给定一个最小连片像元数，去除小于此数的所有连片像元和孤立像元，或者整合到周围的较大的连片异常中。这里采用了3×3像元的窗口对所得异常进行中值滤波，得到总异常图(图 10)。
5. 讨论
为了验证提取的矿化蚀变遥感异常的准确性，特进行野外实地踏勘验证(图 11)。矽卡岩化蚀变验证点多出露透辉石、石榴石、绿泥石、透闪石，多见于褶皱、断裂和岩体与地层接触带中，外带碳酸盐岩裂隙中常见透辉石化、绿帘石化、磁铁矿化、绿泥石化。矽卡岩矿化强烈，一般伴有铜蓝、孔雀石化、磁铁矿化、褐铁矿化；受断层控制明显的矽卡岩化较弱，提取出来的矽卡岩化蚀变信息不明显，大部分还是多见于断裂、褶皱和岩体与地层接触带中。
野外实地踏勘与已有资料表明，多见围岩蚀变有硅化、钠长石化、绿帘石化、绿泥石化、黄铁矿化、角岩化、透闪石化、阳起石化、石榴石化、透辉石化、绢云母化和大理岩化等。Al-OH离子团蚀变矿物蚀变异常多见于绢云母、高岭土、黄铁绢英岩、石英斑岩脉等。蚀变区域里的矿点常有黄铁矿、磁铁矿体产出；Mg-OH离子团蚀变矿物蚀变异常多表现为绿泥石化、绿帘石化、白云母化和碳酸盐化，多见于大理岩和灰岩。
通过野外实地踏勘获得研究区已知矿床、矿(化)点坐标如表 6所示，结合地质图资料(刘渭等，2014)，在ArcGIS10.1软件将已知矿点与蚀变异常信息图进行叠加分析，可以看出提取的蚀变异常与野外验证结果吻合度较高(图 12)。
表 6 虎头崖矿区矿化点坐标(因涉密故总分用*号代替)Table 6. Coordinates of the mineralization points of the Hutouya mining area (The * sign is used instead due to confidentiality reasons)点号 X Y 矿化类型 HT126 E91°38′*″ N37°5′*″ 褐铁矿化、黄铁矿化 HT145 E91°38′*″ N37°4′*″ 铜矿化、磁铁矿化 HT137 E91°37′*″ N37°6′*″ 黄铜矿化 HT25 E91°36′*″ N37°5′*″ 铜铅锌矿化、锡矿化、钨锑矿化 HT27 E91°36′*″ N37°5′*″ 黄铁矿化、黄铜矿化、方铅矿化、铅锌矿化 HT48 E91°36′*″ N37°4′*″ 褐铁矿化、黄铁矿化、磁铁矿化 HT23 E91°37′*″ N37°5′*″ 黄铜矿化、孔雀石化 DY06 E91°36′*″ N37°5′*″ 黄铁矿化、褐铁矿化 DT108 E91°36′*″ N37°4′*″ 矽卡岩化、孔雀石化 HT163 E91°36′*″ N37°4′*″ 闪锌矿化、黄铜矿化 HT189 E91°34′*″ N37°4′*″ 石英斑岩、黄铜矿化、黄铁矿化 HT193 E91°34′*″ N37°4′*″ 黄铁矿化、黄铜矿化 
图 12 虎头崖矿区ASTER遥感图像信息提取总异常与野外验证点叠加图粉红色—成像光谱方法提出的异常信息；深紫色—Al-OH离子团矿化蚀变信息；黄色—Mg-OH离子团矿化蚀变信息；红色—铁染矿化蚀变信息；蓝色—矽卡岩化蚀变矿物信息；黄点—野外验证点坐标Figure 12. Superposition diagram of total anomalies of ASTER remote sensing image extraction and field verification points in the Hutouya mining area6. 结论
(1) 应用成像光谱方法初步得到矿区的矿物分布图，采用ASTER1、ASTER3、ASTER4、ASTER (N)的波段组合进行主成分分析，提取铁染异常信息，矽卡岩化蚀变异常信息，蒙脱石、伊利石与绢云母等矿物的蚀变异常信息，方解石、黑云母、绿泥石等矿物的蚀变异常信息。
(2) 通过野外实地踏勘获得已知矿床、矿(化)点坐标，基于ArcGIS10.1软件将其与基于ASTER数据采用成像光谱法和主成分分析法所提取的蚀变异常信息叠加对比分析，野外验证效果良好。表明提取结果与数据处理方法技术较可靠，对扩大找矿规模和发现富集地段具有重要的指导意义。
(3) 成像光谱方法提取的蚀变异常信息也有不确定性，同时主成分分析法提取出来的矽卡岩化蚀变异常信息不明显，可能是因为USGS标准波谱数据库的矿物波谱没有达到对该区矿物更精细的识别要求，因此应加强野外样品的波谱更加精细地测试工作。
致谢: 长安大学晁会霞博士在论文野外露头照片的矿化类型确认方面给予了悉心指导，同时两位审稿专家的建设性意见对文章的最终定稿帮助很大，在此一并表示感谢。 -
图 1 虎头崖矿区地质简图(刘渭等，2014)
1—第四系；2—蓟县系狼牙山组；3—奥陶志留系滩间山群；4—下石炭统大干沟组；5—上石炭统缔敖苏组；6—上三叠统鄂拉山组；7—花岗闪长岩；8—二长花岗岩；9—石英斑岩；10—似斑状二长花岗岩；11—闪长岩；12—闪长玢岩脉；13—推断断层；14—逆冲断层；15—褶皱；16—矿带；17—青海祁漫塔格成矿带及虎头崖矿区
Figure 1. Simplified geological map of the Hutouya mining area(Liu et al, 2014)
图 2 MNF变换后波段1灰度图像
Figure 2. Grayscale image of the band 1 after the MNF transformation
图 3 MNF变换后波段2灰度图像
Figure 3. Grayscale image of the band 2 after the MNF transformation
图 4 MNF变换后波段3灰度图像
Figure 4. Grayscale image of the band 3 after the MNF transformation
图 7 n-D Mean绘图窗口获取的平均波谱曲线
Figure 7. Average spectral curve obtained by the n-D Mean drawing window
图 8 虎头崖矿区的矿化信息提取图
粉红色—蒙脱石；红色—绿泥石；黄色—方解石
Figure 8. Extraction map of mineralization information in the Hutouya mining area
图 9 研究区ASTER数据主要蚀变矿物的光谱
Figure 9. Spectra of major altered minerals in ASTER data of the study area
图 10 虎头崖矿区ASTER遥感图像信息提取总异常图
深紫色—Al-OH离子团矿化蚀变信息；黄色—Mg-OH离子团矿化蚀变信息；红色—铁染矿化蚀变信息；蓝色—矽卡岩化蚀变矿物信息
Figure 10. Total anomalies of the ASTER remote sensing image extraction in the Hutouya mining area
图 12 虎头崖矿区ASTER遥感图像信息提取总异常与野外验证点叠加图
粉红色—成像光谱方法提出的异常信息；深紫色—Al-OH离子团矿化蚀变信息；黄色—Mg-OH离子团矿化蚀变信息；红色—铁染矿化蚀变信息；蓝色—矽卡岩化蚀变矿物信息；黄点—野外验证点坐标
Figure 12. Superposition diagram of total anomalies of ASTER remote sensing image extraction and field verification points in the Hutouya mining area
表 1 各波段平均值和特征值
Table 1. Average value and eigenvalue of each band
波段名称 平均值 特征值 Band 1 21.173028 28.297260 Band 2 7.034879 9.762023 Band 3 6.628618 4.759365 Band 4 6.013828 2.901527 Band 5 4.625094 2.367740 Band 6 3.821772 1.830600 
下载: 导出CSV
表 2 虎头崖矿区ASTER1、3、4、(5+6)/2主分量变换统计值
Table 2. Statistical values of ASTER1, 3, 4, (5 + 6)/2 principal component transformation in the Hutouya mining area
ASTER1 ASTER3 ASTER4 ASTER(5+6)/2 PC1 0.78201 0.61025 0.12076 0.03833 PC2 -0.62199 0.75182 0.21406 0.04560 PC3 -0.03748 0.24894 -0.91790 -0.30674 PC4 -0.01380 0.01967 -0.31155 0.94993
下载: 导出CSV
表 3 虎头崖矿区ASTER1、3、4、8主分量变换统计值
Table 3. Statistical values of ASTER1, 3, 4, 8 principal component transformation in the Hutouya mining area
ASTER1 ASTER3 ASTER4 ASTER8 PC1 0.78225 0.61043 0.12078 0.02951 PC2 -0.62173 0.75383 0.21140 0.02240 PC3 -0.03460 0.24181 -0.93799 -0.24599 PC4 -0.01824 0.02538 -0.24679 0.96856.
下载: 导出CSV
表 4 虎头崖矿区ASTER1、2、3、4主分量变换统计值
Table 4. Statistical values of ASTER1, 2, 3, 4 principal component transformation in the Hutouya mining area
ASTER1 ASTER2 ASTER3 ASTER4 PC1 0.61068 0.62611 0.47567 0.09382 PC2 0.53342 0.12541 -0.80518 -0.22677 PC3 0.45937 -0.56826 0.02367 0.68227 PC4 0.36264 -0.51898 0.35338 -0.68868
下载: 导出CSV
表 5 虎头崖矿区ASTER1、3、4、5主分量变换统计值
Table 5. Statistical values of ASTER1, 3, 4, 5 principal component transformation in the Hutouya mining area
ASTER1 ASTER3 ASTER4 ASTER5 PC1 0.78211 0.61033 0.12077 0.03498 PC2 -0.62190 0.75288 0.21263 0.03464 PC3 -0.03584 0.24549 -0.92974 -0.27209 PC4 -0.01620 0.02015 -0.27530 0.96101
下载: 导出CSV
表 6 虎头崖矿区矿化点坐标(因涉密故总分用*号代替)
Table 6. Coordinates of the mineralization points of the Hutouya mining area (The * sign is used instead due to confidentiality reasons)
点号 X Y 矿化类型 HT126 E91°38′*″ N37°5′*″ 褐铁矿化、黄铁矿化 HT145 E91°38′*″ N37°4′*″ 铜矿化、磁铁矿化 HT137 E91°37′*″ N37°6′*″ 黄铜矿化 HT25 E91°36′*″ N37°5′*″ 铜铅锌矿化、锡矿化、钨锑矿化 HT27 E91°36′*″ N37°5′*″ 黄铁矿化、黄铜矿化、方铅矿化、铅锌矿化 HT48 E91°36′*″ N37°4′*″ 褐铁矿化、黄铁矿化、磁铁矿化 HT23 E91°37′*″ N37°5′*″ 黄铜矿化、孔雀石化 DY06 E91°36′*″ N37°5′*″ 黄铁矿化、褐铁矿化 DT108 E91°36′*″ N37°4′*″ 矽卡岩化、孔雀石化 HT163 E91°36′*″ N37°4′*″ 闪锌矿化、黄铜矿化 HT189 E91°34′*″ N37°4′*″ 石英斑岩、黄铜矿化、黄铁矿化 HT193 E91°34′*″ N37°4′*″ 黄铁矿化、黄铜矿化
下载: 导出CSV
-
CHEN J, ZHOU P, 2015. Extraction of linear structure and alteration information using ASTER data and prospecting prognosis for Tudun, Hami[J]. Journal of Geomechanics, 21(2):209-217. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlxxb201502011 CHEN J M, SUN W D, YAN B K, et al., 2009. The application and research of the anomaly extraction process base on the ASTER multi-spectral remote sensing in Tianhu iron ore mine[J]. Xinjiang Geology, 27(4):368-372. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xjdz200904014 CHEN Y, HE Z W, DENG H, et al., 2014. Identification and extraction of mineralized alteration information by ASTER remote sensing data[J]. Journal of Guilin University of Technology, 34(1):51-57. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=glgxy201401008 GAO H, ZHANG J L, OUYANG Y, et al., 2013. Lithological classification based on ASTER data by minimum noise fraction transform[J]. Journal of Guilin University of Technology, 33(2):259-265. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=glgxy201302010 GAO W L, ZHANG X J, WANG Z G, et al., 2010. Extraction of lithologic information from the east Kulun orogenic belt using ASTER remote sensing image[J]. Journal of Geomechanics, 16(1):59-69. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlxxb201001008 GAO X F, XIAO P X, XIE C R, et al., 2010. Discussion on tectonic-magmatic process and mineralization in the Qimantage Region[J]. Northwestern Geology, 43(4):119-123. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdz201004014 GENG X X, YANG J M, ZHANG Y J, et al., 2008. The application of ASTER remote sensing data for extraction of alteration anomalies information in shallow overburden area——a case study of the Baoguto porphyry copper deposit intrusion in western Junggar, Xinjiang[J]. Geological Review, 54(2):184-191. (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZLP200802006.htm GU P Y, CHEN R M, CHA X F, et al., 2016. Exploration and practice of 1:50000 geological mapping techniques for alpine-gorge area:a case study in Beishan Area of Wushi, Xinjiang[J]. Journal of Geomechanics, 22(4):837-855. (in Chinese with English abstract) http://www.en.cnki.com.cn/Article_en/CJFDTotal-DZLX201604004.htm HU X H, ZHU G C, LIU H, et al., 2011. Characteristics and mineralization of the Hutouya polymetallic deposit in the Qimantage Metallogenetic Belt[J]. Geology and Prospecting, 47(2):216-221. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dzykt201102011 LI K, GAO Y B, QIAN B, et al., 2015. Geochronology, geochemical characteristics and Hf isotopic compositions of granite in the Hutouya deposit, Qimantag, East Kunlune[J]. Geology in China, 42(3):630-645. (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTotal-DIZI201503017.htm LI R S, JI W H, YANG Y C, et al., 2008. Geology of Kunlun Mountains and adjacent areas[M]. Beijing:Geological Publishing House:1-400. (in Chinese) LI Z M, XUE C J, WANG X H, et al., 2007. Features of regional mineralization and analysis of the exploration development in the eastern Kunlun Mountains[J]. Geological Review, 53(5):708-718. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/OA000004811 LIU W, YANG X K, WANG S L, et al., 2014. Characteristics of ore-controlling structures of the Hutouya orefield in the Qimantage metallogenic belt, Qinghai province[J]. Geology in China, 41(1):222-234. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdizhi201401018 LIU Z, ZHOU P, 2017. Information extraction of alteration of Huaniushan region, Gansu province based on ASTER data[J]. Journal of North University of China (Natural Science Edition), 38(2):124-132, 139. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hbgxyxb201702006 MA S C, FENG C Y, ZHANG D J, et al., 2013. Alteration and mineralization zoning of Hutouya polymetallic deposit in Qimantag area, Qinghai province[J]. Mineral Deposits, 32(1):109-121. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kcdz201301007 MAO X C, LIU W C, DU J G, et al., 2005. Comparison between ETM+ and ASTER data for extraction of alteration information:a case study of Fenghuangshan Orefield, Tongling, Anhui province[J]. Geoscience, 19(2):309-314. (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTOTAL-XDDZ20050200L.htm POUR A B, HASHIM M, 2012. The application of ASTER remote sensing data to porphyry copper and epithermal gold deposits[J]. Ore Geology Reviews, 44:1-9. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=781743811e925019f3090bf4d69edba5 SU S, 2011. Study of mineralization law and analysis of economic value of Hutouya lead-zinc mine in the Qimatage Mine Margine of East Kulun[D]. Beijing: China University of Geosciences (Beijing): 1-127. (in Chinese with English abstract) https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CDFD&dbname=CDFDLAST2015&filename=1011077457.nh&v=MjE2MTJxSkViUElSOGVYMUx1eFlTN0RoMVQzcVRyV00xRnJDVVI3cWZZT2RuRnlIa1ZyckpWRjI2SDdPL0dkWEo= WU Y Z, QIAO G B, CHEN D H, 2011. A preliminary study on relationship between tectonic magmatism and mineralization in Qimantage Area, eastern Kunlun Mountains[J]. Geotectonica et Metallogenia, 35(2):232-241. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ddgzyckx201102008 YAN S X, WU X B, ZHOU C X, et al., 2011. Remote sensing and spectral geology and their applications to mineral exploration[J]. Advances in Earth Science, 26(1):13-29. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dqkxjz201101002 YANG C B, JIANG Q G, ZHU Q, et al., 2009. Altered minerals discriminating using ASTER Data at low vegetation cover area[J]. Geology and Exploration, 45(6):761-766. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dzykt200906019 YANG R H, LI Z Z, CHEN X F, 2012. Information extraction of typical alteration mineral assemblage in porphyry copper using ASTER satellite data, Arequipa province of South Peru[J]. Journal of Geo-information Science, 14(3):411-418. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dqxxkx201203019 YANG X K, HAN K, FAN Y, et al., 2016. Tectonic background and magma-thermodynamic structural features of typical orefields in Qimantage Ore Belt, Eastern Kunlun[J]. Geotectonica et Metallogenia, 40(2):201-212. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ddgzyckx201602002 YAO F J, ZHANG Y J, YANG J M, et al., 2012. Application of ASTER remote sensing data to extraction of alterration zoning information from Dexing porphyry copper deposit[J]. Mineral Deposits, 31(4):881-890. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kcdz201204016 YI L W, LU A H, GU X P, et al., 2019. Evaluation of the methods for quantitative determination of Fe2+/Fe3+ ratios of magnetite from Qimantag metallogenic belts[J]. Journal of Geomechanics, 25(S1):1-4. (in Chinese with English abstract) https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFDLAST2019&filename=DZLX2019S1001&v=MDUwNjg3RGgxVDNxVHJXTTFGckNVUjdxZllPZG5GeUhrVnJ2QklUZkhkckc0SDlpdnJvOUZaWVI4ZVgxTHV4WVM= ZHANG A K, MO X X, LI Y P, et al., 2010. New progress and significance in the Qimantage metallogenic belt prospecting, western Qinghai, China[J]. Geological Bulletin of China, 29(7):1062-1074. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgqydz201007013 ZHANG D M, ZHANG A K, QU G J, et al., 2020. Metallogenic and prospecting model of Kaerqueka Iron-copper polymetallic deposit in western segment of eastern Kunlun[J]. Northwestern Geology, 53(1):91-106. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdz202001009 ZHANG X F, LI Z M, CHEN B, et al., 2012. The contribution of the Tanjianshan group stratum to Skarnization in Qimantage region, Qinghai province[J]. Northwestern Geology, 45(1):40-47. (in Chinese with English abstract) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdz201201006 ZHANG Y J, YANG J M, YAO F J, 2006. The extraction of OHA from different types of deposits by using ASTER data[J]. Mineral Deposits, 25(S1):507-510. (in Chinese with English abstract) ZHANG Y J, YAO F J, 2009. Application study of multi-spectral ASTER data for determination of ETM remote sensing anomaly property:Taking Wulonggou region of eastern Kunlun mountain range as example[J]. Acta Petrologica Sinica, 25(4):963-970. (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTOTAL-YSXB200904020.htm ZHOU D, NI Z Y, YANG Z Y, 2018. Object-oriented remote sensing information extraction method for rocky desertification in Karst area-a case study of Dafang county, Guizhou[J]. Journal of Geomechanics, 24(2):263-273. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dzlxxb201802014 陈建明, 孙卫东, 闫柏琨, 等, 2009. Aster多光谱遥感异常提取在新疆天湖铁矿中的应用[J].新疆地质, 27(4):368-372. http://d.old.wanfangdata.com.cn/Periodical/xjdz200904014 陈菁, 周萍, 2015.基于ASTER的哈密土墩地区构造及蚀变信息提取与找矿预测[J].地质力学学报, 21(2):209-217. https://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20150211&journal_id=dzlxxb 陈晔, 何政伟, 邓辉, 等, 2014.利用ASTER影像识别和提取矿化蚀变信息[J].桂林理工大学学报, 34(1):51-57. http://d.old.wanfangdata.com.cn/Periodical/glgxy201401008 高慧, 张建龙, 欧阳渊, 等, 2013.基于最小噪声分量变换的ASTER遥感数据岩性分类[J].桂林理工大学学报, 33(2):259-265. http://d.old.wanfangdata.com.cn/Periodical/glgxy201302010 高万里, 张绪教, 王志刚, 等, 2010.基于ASTER遥感图像的东昆仑造山带岩性信息提取研究[J].地质力学学报, 16(1):59-69. https://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20100108&journal_id=dzlxxb 高晓峰, 校培喜, 谢从瑞, 等, 2010.祁漫塔格地区构造-岩浆作用与成矿[J].西北地质, 43(4):119-123. http://d.old.wanfangdata.com.cn/Periodical/ddgzyckx201102008 耿新霞, 杨建民, 张玉君, 等, 2008. ASTER数据在浅覆盖区蚀变遥感异常信息提取中的应用:以新疆西准噶尔包古图斑岩铜矿岩体为例[J].地质论评, 54(2):184-191. http://d.old.wanfangdata.com.cn/Periodical/dzlp200802005 辜平阳, 陈瑞明, 查显峰, 等, 2016.高山峡谷区1:50000地质填图技术方法探索与实践:以新疆乌什北山为例[J].地质力学学报, 22(4):837-855. https://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20160404&journal_id=dzlxxb 胡杏花, 朱谷昌, 刘欢, 等, 2011.祁漫塔格矿带虎头崖多金属矿矿床特征和成矿作用分析[J].地质与勘探, 47(2):216-221. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzykt201102011 李侃, 高永宝, 钱兵, 等, 2015.东昆仑祁漫塔格虎头崖铅锌多金属矿区花岗岩年代学、地球化学及Hf同位素特征[J].中国地质, 42(3):630-645. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201503017 李荣社, 计文化, 杨永成, 等, 2008.昆仑山及邻区地质[M].北京:地质出版社:1-400. 李智明, 薛春纪, 王晓虎, 等, 2007.东昆仑区域成矿特征及有关找矿突破问题分析[J].地质评论, 53(5):708-718. http://d.old.wanfangdata.com.cn/Periodical/dzlp200705017 刘渭, 杨兴科, 王守良, 等, 2014.青海省祁漫塔格矿带虎头崖矿田构造控矿特征[J].中国地质, 41(1):222-234. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201401018 刘知, 周萍, 2017.基于ASTER数据的甘肃花牛山地区蚀变信息提取[J].中北大学学报(自然科学版), 38(2):124-132, 139. http://d.old.wanfangdata.com.cn/Periodical/hbgxyxb201702006 马圣钞, 丰成友, 张道俊, 等, 2013.青海虎头崖矽卡岩型多金属矿床蚀变矿化分带特征研究[J].矿床地质, 32(1):109-121. http://d.old.wanfangdata.com.cn/Periodical/kcdz201301007 毛晓长, 刘文灿, 杜建国, 等, 2005. ETM+和ASTER数据在遥感矿化蚀变信息提取应用中的比较-以安徽铜陵凤凰山矿田为例[J].现代地质, 19(2):309-314. http://d.old.wanfangdata.com.cn/Periodical/xddz200502022 苏松, 2011.东昆仑祁漫塔格成矿带虎头崖铅锌矿矿化规律和经济评价[D].北京: 中国地质大学(北京): 1-127. http://cdmd.cnki.com.cn/Article/CDMD-11415-1011077457.htm 伍跃中, 乔耿彪, 陈登辉, 2011.东昆仑祁漫塔格地区构造岩浆作用与成矿关系初步探讨[J].大地构造与成矿学, 35(2):232-241. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ddgzyckx201102008 燕守勋, 武晓波, 周朝宪, 等, 2011.遥感和光谱地质进展及其对矿产勘查的实践应用[J].地球科学进展, 26(1):13-29. http://d.old.wanfangdata.com.cn/Periodical/mtjs201111059 杨长保, 姜琦刚, 朱群, 等, 2009.浅覆盖区基于ASTER数据的蚀变矿物识别方法[J].地质与勘探, 45(6):761-766. http://d.old.wanfangdata.com.cn/Periodical/dzykt200906019 杨日红, 李志忠, 陈秀法, 2012. ASTER数据的斑岩铜矿典型蚀变矿物组合信息提取方法:以秘鲁南部阿雷基帕省斑岩铜矿区为例[J].地球信息科学学报, 14(3):411-418. http://d.old.wanfangdata.com.cn/Periodical/dqxxkx201203019 杨兴科, 韩珂, 范阅, 等, 2016.东昆仑祁漫塔格矿带典型矿田构造背景与岩浆-热力构造特征[J].大地构造与成矿学, 40(2):201-212. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ddgzyckx201602002 姚佛军, 张玉君, 杨建民, 等, 2012.利用ASTER提取德兴斑岩铜矿遥感蚀变分带信息[J].矿床地质, 31(4):881-890. http://d.old.wanfangdata.com.cn/Periodical/kcdz201204016 易立文, 鲁安怀, 谷湘平, 等, 2019.祁漫塔格成矿带磁铁矿Fe2+/Fe3+比值测定方法对比研究[J].地质力学学报, 25(S1):1-4. https://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=2019S101&journal_id=dzlxxb 张爱奎, 莫宣学, 李云平, 等, 2010.青海西部祁漫塔格成矿带找矿新进展及其意义[J].地质通报, 29(7):1062-1074. http://d.old.wanfangdata.com.cn/Periodical/zgqydz201007013 张大明, 张爱奎, 屈光菊, 等, 2020.东昆仑西段卡而却卡铁铜多金属矿床成矿模式及找矿模型[J].西北地质, 53(1):91-106. http://d.old.wanfangdata.com.cn/Periodical/xbdz202001009 张晓飞, 李智明, 陈博, 等, 2012.东昆仑祁漫塔格地区滩间山群矽卡岩化成矿作用[J].西北地质, 45(1):40-47. http://d.old.wanfangdata.com.cn/Periodical/xbdz201201006 张玉君, 杨建民, 姚佛军, 2006.用ASTER数据进行不同类型矿床蚀变异常提取研究[J].矿床地质, 25(S1):507-510. http://d.old.wanfangdata.com.cn/Conference/6278224 张玉君, 姚佛军, 2009.应用多光谱ASTER数据对ETM遥感异常的定性判别研究:以东昆仑五龙沟为例[J].岩石学报, 25(4):963-970. http://d.old.wanfangdata.com.cn/Periodical/ysxb98200904020 周迪, 倪忠云, 杨振宇, 2018.面向对象的喀斯特地区石漠化遥感信息提取研究:以贵州省大方地区为例[J].地质力学学报, 24(2):263-273. https://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20180213&journal_id=dzlxxb -
下载:
下载:
计量
- 文章访问数: 877
- HTML全文浏览量: 208
- PDF下载量: 62
- 被引次数: 0
