基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法

欧阳宁 ,  黄辰钰 ,  林乐平

吉林大学学报(工学版) ›› 2025, Vol. 55 ›› Issue (11) : 3727 -3735.

PDF (2357KB)
吉林大学学报(工学版) ›› 2025, Vol. 55 ›› Issue (11) : 3727 -3735. DOI: 10.13229/j.cnki.jdxbgxb.20240252
计算机科学与技术

基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法

作者信息 +

Hyperspectral image classification based on hybrid spectral enhancement and multi-scale spatial aggregation

Author information +
文章历史 +
PDF (2412K)

摘要

由于高光谱图像存在同物异谱和异物同谱现象,仅依赖光谱信息无法充分表征高光谱图像的特征,因此可引入空间信息以更准确地捕捉物体特征。为此,本文提出一种基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法。该方法设计了混合光谱增强模块,利用小波变换构建光谱的多尺度局部特征,通过Transformer架构生成光谱的全局特征,以增强光谱特征的类内一致性。同时,设计了多尺度空间聚合模块,用于提取空间特征固有的多尺度信息,并建立不同尺度间的交互关系,以生成更具鲁棒性的土地覆盖表示,从而进一步提升分类性能。实验结果表明:本文方法相较于其他先进网络表现出显著的优越性,能有效获取更丰富的光谱信息和空间特征表示。

Abstract

Due to the presence of the same object with different spectrum and the same spectrum with different objects in hyperspectral images, spectral information alone cannot fully reflect the features of hyperspectral images, and spatial information can be introduced to capture the features of objects more accurately. Therefore, a hyperspectral image classification method based on hybrid spectral enhancement and multi-scale spatial aggregation is proposed in this paper. In this method, a hybrid spectral enhancement module is designed, multi-scale local features of the spectrum are constructed using wavelet transform, and global features of the spectrum are generated by Transformer architecture, so as to enhance the intra-class consistency of the spectral features. At the same time, a multi-scale spatial aggregation module is designed to extract the inherent multi-scale information of spatial features and establish the interaction between different scales, so as to generate a more robust land cover representation, thereby further improve the classification performance. The experimental results show that the proposed method is superior to other advanced networks, indicating that the method can effectively obtain more abundant spectral information and spatial feature representation.

Graphical abstract

关键词

高光谱图像分类 / 混合光谱增强模块 / 小波变换 / 多尺度空间聚合模块

Key words

hyperspectral image classification / hybrid spectral enhancement module / wavelet transform / multi-scale spatial aggregation module

引用本文

引用格式 ▾
欧阳宁,黄辰钰,林乐平. 基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法[J]. 吉林大学学报(工学版), 2025, 55(11): 3727-3735 DOI:10.13229/j.cnki.jdxbgxb.20240252

登录浏览全文

4963

注册一个新账户 忘记密码

0 引 言

高光谱图像以数据立方体的形式呈现,是一种包含丰富光谱和空间信息的图像,拥有数百个相邻的光谱波段。与常规的RGB图像相比,高光谱图像有更精细的光谱信息,在分类任务中具有很大的潜力。但光谱波段的高维度、冗余性和标记样本数量的有限性,也给高光谱图像分类带来了挑战。特别是在土地覆盖复杂时,类内光谱一致性较弱,原始数据空间难以对其进行有效表达。

随着深度学习的发展,卷积神经网络(Convolutional neural network, CNN)凭借其对光谱空间特征的局部建模能力,在高光谱图像分类领域受到广泛关注。基于CNN的方法12通过叠加卷积核的方式,对光谱空间特征进行从浅层到深层的特征提取。在多层次空-谱融合网络3中,通过结合不同层次的空间与光谱信息,能有效地获取高光谱图像的精细特征。在空谱分组卷积密集网络4中,采用分组卷积重复利用高分辨率信息,可防止细节信息丢失。Hamida等5提出了一种新的三维深度学习方法,实现了光谱和空间信息的联合处理。因此,基于深度学习的分类网络对高光谱图像的空间和光谱信息有较好的表征能力,能有效提取更具判别性的分类特征,获得更优越的分类结果。

基于CNN的方法虽然具有强大的局部特征捕捉能力,但其对全局信息的关注仍显不足。Dosovitskiy等6提出视觉Transformer(Vision transformer, ViT)网络,将Transformer的特性引入自然图像处理中,此后越来越多的学者开始将ViT应用于高光谱图像分类领域。Hong等7提出SpectralFormer网络,从分组的相邻波段中学习光谱信息作为ViT的输入,并构建跨层编码输出的相互连接,深度挖掘和表示光谱特征的序列属性;Sun等8设计了Transformer框架,对浅层的空-谱特征进行特征变换,并输入Transformer构建特征表示。总的来说,ViT具有全局感受野,可获得像素间的远距离依赖关系,能为所有像素分类提供更精准的注意力权重,从而避免遗漏重要信息。

为充分发挥CNN的局部接受域优势和Transformer的全局特征提取优势,本文将CNN与ViT相结合,提出了基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法(Hyperspectral image classification based on hybrid spectral enhancement and multi-scale spatial aggregation, HEMA)。根据高光谱图像数据特点,本文借鉴小波变换思想和ViT网络特性,构建了混合光谱增强模块(Hybrid spectral enhancement module, HSEM),通过配置模块参数,在保留空间信息的前提下有效增强了光谱的特征信息。同时,在此基础上构建了多尺度空间聚合模块(Multi-scale spatial aggregation module, MSAM),通过依次相加的聚合构造函数,增强了不同尺度特征映射间的交互性,提升了网络的特征表达能力。实验结果表明,相较于一些前沿的深度网络方法,本文方法可有效提升高光谱图像分类网络的性能,进而得到较高的分类精度。

1 相关研究

原始高光谱图像数据光谱维度较高,当数据非常复杂时,相对较差的类内光谱一致性对分类结果有很大影响。因此,降低光谱变化带来的影响是高光谱图像分类的关键问题之一。正交变换方法,如主成分分析(Principal component analysis, PCA),可通过正交变换识别光谱波段间的相关性,进而提取其内在特征9。虽然正交变换方法在提取主成分方面具有一定的优势,但在捕捉非线性特征方面相对较弱,特别是在处理高度变化的光谱信息时灵活性不足。小波变换(Discrete wavelet transform, DWT)10是一种多层次分解信号的方法,能揭示信号在不同时间尺度上的频率成分。小波变换凭借其在频率域分布规律处理和多尺度分析方面的优势脱颖而出,尤其在处理具有复杂光谱结构的场景时表现更为出色。Tang等11使用小波变换对光谱数据进行处理,将其作为特征提取网络的输入,以减少异常光谱数据对分类的影响。Ahmad等12利用小波变换进行可逆下采样以保持数据的完整性,并设计基于小波的CNN提取多尺度空间光谱特征。小波变换的多尺度分解特性可有效捕捉光谱信息的详细结构和关键特征。因此,本文引入小波变换对光谱数据进行多尺度分解,将光谱分解成各种近似分量和细节分量,获取分离频谱中更细节的特征。当类的频谱相似时,添加该小波特征可增强类间区分能力。

为了获取更多感受野,现如今很多具有多尺度光谱空间特征结构的CNN被应用于高光谱图像分类13-15。Ghaderizadeh等16通过构建多尺度残差模块获取更多的局部和全局特征;Wang等17利用多尺度块考虑多范围上下文信息。受多尺度分析在提升模型适应性方面的启发,本文构建了多尺度空间聚合模块。与简单叠加感受野不同,该模块通过对多个感受野进行逐特征融合,充分整合了多尺度信息。这使得模型能更全面、准确地捕捉空间特征中的多样化信息,进而显著提升对复杂图像的识别和分类能力。

本文将混合光谱增强模块和多尺度空间聚合模块相结合,提出了基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法,有效提升了高光谱图像的分类性能和效率。

2 本文方法

2.1 混合光谱增强与多尺度空间聚合网络分类框架

本文所提HEMA的框架如图1所示。首先,该网络以高光谱图像像素n的图像邻域块作为输入,将其转换为二维向量PnR(w×w)×1×B,再使用混合光谱增强模块对光谱信息进行信息增强和全局光谱特征提取,得到混合增强特征SnR(w×w)×1×B,旨在降低类内光谱变化对分类的影响。由于将三维高光谱数据直接展平为二维矩阵会忽略高光谱图像原本的结构特性,为从光谱特征中恢复鲁棒性更强的空间特征,本文对二维特征状态下的Sn执行Reshape操作,将Sn重塑为原始空间尺寸的三维特征图InRw×w×B。其次,通过多尺度空间聚合模块获取包含多尺度空间信息的特征FnRw×w×b,提高空间特征的表达能力。最后,使用多层感知机MLP对Fn进行分类,并采用交叉熵构建损失函数,用于度量预测标签与真实标签之间的差异。

2.2 混合光谱增强模块

本文设计的混合光谱增强模块结构如图2所示。光谱特征提取是构建本文网络的重要环节,会直接影响后续的分类结果。由于高光谱图像数据频带数量多,冗余性较强,若直接将光谱波段输入Transformer中,并不能充分利用其自注意力机制的优势。因此,本文将原始光谱信息与小波变换处理后的光谱信息相结合,共同输入Transformer中,目的是在保留原始光谱特征的前提下,降低类内光谱变化对分类的影响,获取包含更多细节的光谱特征Xn。本文采用经典的Haar小波,如图3所示。具体来说,首先使用低通滤波器fL=(1/2,1/2)和高通滤波器fH=(1/2,-1/2)对输入频谱带I沿水平方向进行滤波和下采样,生成子带系数矩阵ILIH。其次采用同样的低通和高通滤波器,沿ILIH的垂直方向进行过滤和下采样,生成4个亚带矩阵ILLILHIHLIHH。最后,对ILL继续分解,以获得更细节的频谱特征,最终取得的光谱序列分别表示为ILLJ,[IHLk,ILHk,IHHk](k=1,2,,J),本实验中J=3。由于小波变换具有很好的时频局域性,可将光谱的不同底层特征变换成对应的小波系数,有利于提取更丰富的特征。

之后,引入位置编码EspeposR1×B,将其与Xn相加,得到新的位置嵌入ZnR(w×w)×2×B,再将Zn输入Transformer编码器。其中,自注意力机制的公式为:

Q=Zn·WQ=(Xn+Espepos)WQ
K=Zn·WK=(Xn+Espepos)WK
V=Zn·WV=(Xn+Espepos)WV
Attention=softmaxQKTdV=
softmaxZnWQWKTZnTdV

式中:QKV分别为查询向量、键向量和值向量。

Transformer编码器的表达式为:

Anl=Tnl-1+MHALN(Tnl-1)
Tnl=Anl+MLPLN(Anl)

式中:Tnl为空间位置n的第l个Transformer编码器输出;MHA和MLP分别为多头注意和多层感知器;LN为归一化层;当l=1时,Tnl-1=Zn

为进一步聚合不同Transformer编码器的尺度特征,获得更有效的特征表达,本文参考了SpectralFormer中的跨层自适应融合(Context-aware fusion, CAF)模块。该模块的输出结果为:

T¯nl=w¨TnlTnl-2

式中:w¨R1×2为自适应融合的可学习网络参数。

2.3 多尺度空间聚合模块

本文设计的多尺度空间聚合模块结构如图4所示。首先,通过一维卷积对特征图In进行处理,提取浅层空间特征并减少特征映射的通道数,生成新的空间模块ĬRb×w×w×C。其次,为进一步增强光谱特征信息,对Ĭ执行通道洗牌操作,包括重塑、转置和再次重塑。完成一系列迭代洗牌后,原本的通道顺序被打乱,通道间建立了新的关联。这一过程最终生成新的特征向量Sshuffle,它在保留原始信息的同时,增强了通道间的交互和融合。

较大的卷积核能捕捉更丰富的特征信息,但可能导致物体边界处的预测误差增大。因此,在设计卷积核尺寸时需进行权衡,既要确保能有效提取空间光谱特征,又要保证有足够的感受野以捕捉更广泛的特征信息。为拓展感知范围,本文首先采用不同尺度的深度可分离卷积,提取包含4个接收域的空间特征S1S2S3S4Rb×w¯×w¯×c,在提高效率的同时更全面地捕捉不同尺度的空间信息。其次,引入了依次相加的聚合构造函数,通过元素加法操作改善了不同尺度间的交互作用。再次,通过通道拼接策略生成了多尺度融合的空间特征图,计算多尺度空间信息 fRb×w¯×w¯×4c

f=cat(S1,S1+S2,S2+S3,S3+S4)

最后,通过全局平均池化处理和相应的维度处理,得到最终的特征信息y=Rb×4c,增强模型的鲁棒性。

2.4 数据集

在网络模型训练中,采用以下两个高光谱数据集验证本文方法的有效性:

Indian pines (IN)数据集18是由AVIRIS传感器在印第安纳州西北部采集的影像,波长范围为400~2 500 nm,去除20个噪声严重的波段后共有200个波段,图像尺寸为145×145像素,空间分辨率约为20 m。该场景涉及16种类别,其中包括许多天然植被。

WHU-hi-longkou (LK)数据集1920是由Headwall Nano-Hyperspec图像传感器在中国湖北省龙口镇采集的影像,波长范围为400~1 000 nm,共有270个波段,图像尺寸为550×400像素,空间分辨率约为0.109 m。该场景涉及9种类别,其中包括玉米、棉花、芝麻、阔叶大豆、窄叶大豆和水稻6种作物。

这两个数据集都具有较长的波段,光谱波段间存在较高的冗余性,往往会使图像的类内一致性变得较弱,这为研究和解决类内一致性较弱的问题提供了试验场景。本文以此为基础,设计了HEMA方法以应对高光谱图像分类中类内一致性较弱的问题。

2.5 损失函数

在网络模型的训练过程中,采用交叉熵作为损失函数,用于度量预测标签向量与实际标签向量之间的差异,以评估本文方法的泛化性能。交叉熵的计算公式为:

Loss=-1Nyklog(y¯k)+(1-yk)log(1-y¯k)

式中:N为类别数;yk为真实样本标签中属于第k类的概率值;y¯k为预测样本标签中属于第k类的概率值。

3 实验仿真

3.1 实验配置

实验使用少量样本进行训练,在IN数据集中,除第1、7、9类只选择15个样本作为训练集外,其他类别都随机选择50个样本作为训练集;在LK数据集中,所有类别都随机选择50个样本作为训练集。设置批量大小为64,学习率为0.000 5,对所有数据集进行10次训练,并采用总体准确度(Overall accuracy, OA)、平均准确度(Average accuracy, AA)、Kappa系数及单类分类精度作为模型性能评估指标,以获得对结果更准确的估计。实验平台配置为RTX2080Ti显卡及Windows操作系统,并采用PyTorch深度学习框架实现。

3.2 消融实验

为验证各模块的有效性,本文在IN数据集上通过添加不同模块进行消融实验。不同模块组成下的OA、AA及Kappa值如表1所示,可以看出,SpectralFormer的分类效果值最低;将原SpectralFormer的输入部分替换为引入小波变换的光谱增强信息后,OA提升了2.02百分点。这证实了小波变换的引入能更有效地利用每个频谱中包含的丰富信息,提高分类准确性。多尺度空间聚合模块的引入,与基线相比OA提高了7.26百分点,说明相较于仅进行频谱特征提取的SpectralFormer,增加多尺度空间聚合模块有助于提高模型性能。在同时包含小波变换和多尺度空间聚合模块的实验中,各评估指标较SpectralFormer都提高了10百分点以上,说明了两个模块的必要性。综上所述,各模块均对网络性能有积极影响。

4 实验结果与分析

4.1 结果对比分析

将HEMA方法与最新的基于深度学习和Transformer的高光谱图像分类方法进行对比,对比算法包括ABLSTM21、CNN3D5、ViT6、SpectralFormer7、SSFTT8。对比算法均使用相同规模的训练数据,并从原始数据中选取大小为7×7×B的立方体作为输入,其中B为光谱维度。

HEMA与对比方法的分类结果如表2表3所示,其中加粗数字为该行最大值,即性能最优。从表中可以看出,HEMA的分类精度相较于ABLSTM和CNN3D明显提升,OA、AA及Kappa值在IN数据集中至少提高了16.91百分点、21.1百分点和19.15,在LK数据集中至少提高了5.28百分点、13.37百分点和6.75。此外,ABLSTM和CNN3D的性能整体低于基于Transformer的方法,表明Transformer在获取全局特征方面具有优势;改进后的ViT方法能更好地提取全局和局部空-谱特征,获取更高级的语义信息。相较于ViT、SpectralFormer和SSFTT等基于Transformer的方法,HEMA的OA、AA和Kappa值均有提升,在IN数据集中至少提高了11.03百分点、8.28百分点和2.5,在LK数据集中至少提高了1.37百分点、5.16百分点和1.77。此外,HEMA在大多数类别上的分类精度都优于对比算法。这表明HEMA在原始光谱信息基础上引入小波变换,可从多尺度光谱信息中获取噪声点更小的光谱信息,并增强类内光谱一致性。同时,多尺度空间聚合模块对空间信息的多尺度特征提取,促进了不同尺度空间信息的交互,从而获得更细节的特征信息,提升了网络的分类精度。

4.2 分类效果图对比分析

HEMA方法与对比方法在IN数据集和LK数据集上的高光谱图像分类效果如图5图6所示。可以观察到,HEMA方法生成的分类效果图最为清晰,这意味着其正确分类的样本点数量最多。如图5(f)相较于图5(a)~(e),HEMA将Soybean-mintill类别错分为Corn-mintill类别的像素点更少。相应地,如表2表3所示,HEMA方法在两个数据集的各个类别上的分类精度几乎都达到了最佳,并且其AA也达到了最优。由此可见,HEMA能更有效地捕捉光谱的局部和全局特征,以及空间的多尺度交互信息,最终得到的空-谱特征使分类效果图能更准确地表示真实的分类。

4.3 训练和测试时间对比分析

HEMA和对比方法的平均训练时间和平均测试时间如表4所示。可以看出,HEMA的训练时间比ABLSTM和SpectralFormer更短,这是因为HEMA在进行空间特征提取时采用了深度可分离卷积,减轻了网络的负担。CNN3D、ViT和SSFTT的训练时间则比HEMA更短,这是因为对比网络的规模较小,本文采用的多尺度空间聚合模块增加了网络的复杂性,导致训练时间和测试时间并非最短。但HEMA通过增加一定的训练时间获取了更高的精确度,在效率和准确性之间实现了很好的平衡,具有一定的应用优势。

5 结束语

为综合利用高光谱图像中的光谱信息和空间信息,本文提出了一种基于混合光谱增强与多尺度空间聚合的高光谱图像分类方法。其中,混合光谱增强模块能有效增强光谱信息,并提取全局光谱特征,从而更全面地表征复杂的光谱特性;多尺度空间聚合模块能对不同尺度的空间信息进行交互融合,使不同尺度的特征相互引导,从而有效提高空间特征的表达能力。实验结果表明,本文提出的方法能更有效地提升高光谱图像分类任务的准确性和鲁棒性。未来研究将在保持分类精度的前提下,考虑如何降低网络模型的复杂度并提高计算效率。

参考文献

[1]

Yang K, Sun H, Zou C B, et al. Cross-attention spectral-spatial network for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 1-14.

[2]

Ahmad M, Khan A M, Mazzara M, et al. Hybrid dense network with attention mechanism for hyperspectral image classification[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2022, 15: 3948-3957.

[3]

欧阳宁,李祖锋,林乐平.基于多层次空-谱融合网络的高光谱图像分类[J].吉林大学学报:工学版,2022,52(10):2438-2446.

[4]

Ouyang Ning, Li Zu-feng, Lin Le-ping. Hyperspectral image classification based on a multi-level spatial-spectral fusion network[J]. Journal of Jilin University (Engineering and Technology Edition), 2022,52(10):2438-2446.

[5]

欧阳宁,李祖锋,林乐平.基于空谱分组卷积密集网络的高光谱图像分类[J].计算机工程与设计,2022,43(7):2031-2039.

[6]

Ouyang Ning, Li Zu-feng, Lin Le-ping. Hyperspectral image classification based on a spatial-spectral group convolution dense network[J]. Computer Engineering and Design, 2022,43(7):2031-2039.

[7]

Hamida A B, Benoit A, Lambert P, et al. 3-D deep learning approach for remote sensing image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(8): 4420-4434.

[8]

Dosovitskiy A, Beyer L, Kolesnikov A, et al. An image is worth 16×16 words: Transformers for image recognition at scale[J]. Arxiv Preprint, 2010, 11: 201011929.

[9]

Hong D F, Han Z, Yao J, et al. Spectralformer: rethinking hyperspectral image classification with transformers[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 1-15.

[10]

Sun L, Zhao G R, Zheng Y H, et al. Spectral–spatial feature tokenization Transformer for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 1-14.

[11]

Uddin M P, Mamun M A, Hossain M A. Feature extraction for hyperspectral image classification[C]∥ IEEE Region 10 Humanitarian Technology Conference. Piscataway, NJ: IEEE. 2017: 379-382.

[12]

Zhang D S. Wavelet Transform[J]. Fundamentals of Image Data Mining: Analysis, Features, Classification and Retrieval, 2019: 35-44.

[13]

Tang Y G, Xie X L, Yu Y H. Hyperspectral classification of two-branch joint networks based on gaussian pyramid multiscale and wavelet transform[J]. IEEE Access, 2022, 10: 56876-56887.

[14]

Ahmad M, Ghous U, Usama M, et al. WaveFormer: spectral-spatial wavelet Transformer for hyperspectral image classification[J]. IEEE Geoscience and Remote Sensing Letters, 2024, 21: 1-5.

[15]

Wang D, Du B, Zhang L P, et al. Adaptive spectral-spatial multiscale contextual feature extraction for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 59(3): 2461-2477.

[16]

Gao H M, Yang Y, Li C M, et al. Multiscale residual network with mixed depthwise convolution for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 59(4): 3396-3408.

[17]

Leng Q M, Yang H O, Jiang J J, et al. Adaptive multiscale segmentations for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(8): 5847-5860.

[18]

Ghaderizadeh S, Abbasi M D, Sharifi A, et al. Multiscale dual-branch residual spectral-spatial network with attention for hyperspectral image classification[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2022, 15: 5455-5467.

[19]

Wang X, Tan K, Du P J, et al. A unified multiscale learning framework for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 1-19.

[20]

Baumgardner M F, Biehl L L, Landgrebe D A. 220 band aviris hyperspectral image data set: June 12, 1992 Indian pine test site 3[J]. Purdue University Research Repository, 2015, 10(7): 991.

[21]

Zhong Y F, Hu X, Luo C, et al. WHU-Hi: UAV-borne hyperspectral with high spatial resolution (H2) benchmark datasets and classifier for precise crop identification based on deep convolutional neural network with CRF[J]. Remote Sensing of Environment, 2020, 250: 112012.

[22]

Zhong Y F, Wang X Y, Xu Y, et al. Mini-UAV-borne hyperspectral remote sensing: From observation and processing to applications[J]. IEEE Geoscience and Remote Sensing Magazine, 2018, 6(4): 46-62.

[23]

Mei S H, Li X G, Liu X, et al. Hyperspectral image classification using attention-based bidirectional long short-term memory network[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 1-12.

基金资助

国家自然科学基金项目(62001133)

广西科技基地和人才专项项目(桂科 AD19110060)

广西自然科学基金项目(2017GXNSFBA198212)

广西无线宽带通信与信号处理重点实验室基金项目(GXKL06200114)

AI Summary AI Mindmap
PDF (2357KB)

4

访问

0

被引

详细

导航
相关文章

AI思维导图

/