1、遥感地质学资源学院数地所1参考资料参考资料2主要内容主要内容一一 遥感图像处理二二地质解译标志与遥感图像地学分析方法三三遥感图像地貌解译及其应用四四遥感岩性解译与编图五五遥感构造解译与编图3 1 1 遥感图像处理遥感图像处理o遥感地质学研究方法遥感地质学研究方法o遥感物理基础遥感物理基础o遥感图像类型与特性遥感图像类型与特性o遥感图像处理方法遥感图像处理方法4 1.1 1.1 研究方法研究方法o遥感地质学使用的方法，涉及地物波谱测试方遥感地质学使用的方法，涉及地物波谱测试方法、数理统计、模拟试验、模式识别、数字图法、数理统计、模拟试验、模式识别、数字图像处理、像处理、GIS、GPS及地学（包括
2、地质、地理、及地学（包括地质、地理、地貌，地图学等）的有关研究分析方法等。地貌，地图学等）的有关研究分析方法等。5图图 1.1 遥感地质工作流程图遥感地质工作流程图6 1.2 1.2 遥感物理基础遥感物理基础o电磁波谱与电磁辐射电磁波谱与电磁辐射o太阳辐射及大气对辐射的影响太阳辐射及大气对辐射的影响o地球的辐射与地物波谱地球的辐射与地物波谱7一一 电磁波谱与电磁辐射电磁波谱与电磁辐射1 1 电磁波电磁波n电磁波是指在真空或物质中通过传播电磁场的振动电磁波是指在真空或物质中通过传播电磁场的振动而传输电磁能量的波。而传输电磁能量的波。当电磁振荡进入空间时，变化的磁场激发了涡当电磁振荡进入空间时，变
3、化的磁场激发了涡旋电场，变化的电场激发了涡旋磁场，使电磁振荡旋电场，变化的电场激发了涡旋磁场，使电磁振荡在空间传播，形成电磁波，也称电磁辐射在空间传播，形成电磁波，也称电磁辐射.n在真空以光速传播在真空以光速传播n电磁辐射具有波粒二象性电磁辐射具有波粒二象性 8一一 电磁波谱与电磁辐射电磁波谱与电磁辐射2 电磁波谱电磁波谱 按照电磁波在真空中传播的波长或频率递增或递按照电磁波在真空中传播的波长或频率递增或递减排列形成的一个连续谱带（图减排列形成的一个连续谱带（图1.2)1.2)。电磁波的波段从波长短的一侧开始，电磁波的波段从波长短的一侧开始，依次为：依次为：射线射线XX射线射线紫外线紫外线可见
4、光可见光红外线红外线无线电波无线电波。9图图 1.2 1.2 电磁波谱电磁波谱102.2.电磁波谱电磁波谱 目前，遥感所使用的电磁波的波段有：目前，遥感所使用的电磁波的波段有：o紫外线：紫外线：波长范围为波长范围为0.010.010.38m0.38m，太阳光谱中，只，太阳光谱中，只有有0.30.30.38m0.38m波长的光到达地面，对油污染敏感，波长的光到达地面，对油污染敏感，但探测高度在但探测高度在2000 m2000 m以下。以下。v可见光：可见光：波长范围为波长范围为0.380.380.76m0.76m，人眼对可见光有，人眼对可见光有敏锐的感觉，是遥感技术应用中的重要波段。敏锐的感觉，
5、是遥感技术应用中的重要波段。v红外线红外线：波长范围为波长范围为0.760.761000m1000m，根据性质分为近，根据性质分为近红外、中红外、远红外和超远红外。红外、中红外、远红外和超远红外。v微波：微波：波长范围为波长范围为1 1 mmmm1 m1 m，穿透性好，不受云雾的，穿透性好，不受云雾的影响。影响。一一 电磁波谱与电磁辐射电磁波谱与电磁辐射11一一 电磁波谱与电磁辐射电磁波谱与电磁辐射3 电磁辐射的度量电磁辐射的度量辐射能量辐射能量（W W）：电磁辐射的能量，单位：）：电磁辐射的能量，单位：J J辐射通量辐射通量（）：单位时间内通过某一面积的辐射）：单位时间内通过某一面积的辐射能
6、量，能量，=dw/dtdw/dt，单位是，单位是w w辐射通量密度辐射通量密度（E E）：单位时间通过单位面积上的辐）：单位时间通过单位面积上的辐射能量，射能量，E=d/E=d/dsds，单位：，单位：w/m2w/m2，S S为面积为面积辐照度辐照度（I I）：被辐射物体表面单位面积上的辐射通）：被辐射物体表面单位面积上的辐射通量，量，I=E=d/I=E=d/dsds，单位：，单位：w/m2w/m2，S S为面积为面积辐射出射度辐射出射度（M M）：辐射源物体表面单位面积上的）：辐射源物体表面单位面积上的辐射通量，辐射通量，M=E=d/M=E=d/dsds，单位：，单位：w/m2w/m2，S
7、S为面积为面积辐射亮度辐射亮度(L L):):辐射源在某一方向，辐射源在某一方向，单位投影表面，单位投影表面，单位立体角内的辐射通量。单位立体角内的辐射通量。12一一 电磁波谱与电磁辐射电磁波谱与电磁辐射4 绝对黑体及辐射规律绝对黑体及辐射规律对任何波长的电磁辐射全部吸收的物体。对任何波长的电磁辐射全部吸收的物体。黑体吸收系数黑体吸收系数（，T T）=100%=100%黑体反射系数黑体反射系数（，T T）=0%=0%普朗克公式普朗克公式斯忒藩斯忒藩-玻尔兹曼定律玻尔兹曼定律维恩位移定律维恩位移定律自然界中黑体辐射是不存在，一般地物辐射能量总要自然界中黑体辐射是不存在，一般地物辐射能量总要比黑体
8、辐射能量小。如果利用黑体辐射有关公式，则比黑体辐射能量小。如果利用黑体辐射有关公式，则需要增加一个因子，那就是发射率（比辐射率）。需要增加一个因子，那就是发射率（比辐射率）。依据发射率与波长的关系，将地物分为三种类型：黑依据发射率与波长的关系，将地物分为三种类型：黑体、灰体、选择性发射体。体、灰体、选择性发射体。13二二 太阳辐射及大气对辐射的影响太阳辐射及大气对辐射的影响太阳是被动遥感主要的辐射源，又叫太阳光。太阳是被动遥感主要的辐射源，又叫太阳光。o1 1 太阳辐射太阳辐射 太阳的光谱通常指光球产生的光谱，光球发射太阳的光谱通常指光球产生的光谱，光球发射的能量集中在可见光波段（图的能量集中
9、在可见光波段（图1.3)1.3)。图图1.3 1.3 太阳辐照太阳辐照度分布曲线度分布曲线14二二 太阳辐射及大气对辐射的影响太阳辐射及大气对辐射的影响大气吸收减弱大气吸收减弱大气散射减弱大气散射减弱云层反射减弱云层反射减弱大气折射变向大气折射变向地物吸收减弱地物吸收减弱传感器传感器太阳辐射太阳辐射地物地物地物漫反射减弱地物漫反射减弱地物镜面反射减弱地物镜面反射减弱2 2 大气对辐射的影响大气对辐射的影响15二二 太阳辐射及大气对辐射的影响太阳辐射及大气对辐射的影响3 大气窗口大气窗口o电磁波通过大气层时较少被反射，吸收和散射的，透射率较电磁波通过大气层时较少被反射，吸收和散射的，透射率较高的
10、波段称为高的波段称为大气窗口大气窗口。大气窗口大气窗口波段波段透射率透射率/%/%应用举例应用举例紫外可见光紫外可见光近红外近红外0.30.31.3 m1.3 m9090TM1-4TM1-4、SPOTSPOT的的HRVHRV近红外近红外1.51.51.8 m1.8 m8080TM5TM5近近-中红外中红外2.02.03.5 m3.5 m8080TM7TM7中红外中红外3.53.55.5 m5.5 m3.4-4.2m3.4-4.2m为为90%90%4.6-4.9m4.6-4.9m为为 50-60%50-60%NOAANOAA的的AVHRRAVHRR远红外远红外8 814 m14 m6060707
11、0TM6TM6微波微波0.80.82.5cm2.5cm100100RadarsatRadarsat16三 地球的辐射与地物波谱地球的辐射与地物波谱图图1.5 太阳与地表辐射的电磁波谱太阳与地表辐射的电磁波谱17三 地球的辐射与地物波谱地球的辐射与地物波谱1 1 地球的辐射源地球的辐射源在在0.30.32.5m2.5m波段（主要在可见光和近红外波段），地波段（主要在可见光和近红外波段），地表以反射太阳辐射为主，表以反射太阳辐射为主，地球自身的辐射可以忽略。即地球自身的辐射可以忽略。即在该波段范围内，对地观测遥感主要以太阳的短波辐射在该波段范围内，对地观测遥感主要以太阳的短波辐射对地表进行探测和成
12、像。对地表进行探测和成像。在在2.52.56.0m6.0m波段（主要在中红外波段），地表反射太波段（主要在中红外波段），地表反射太阳辐射和地球自身的热辐射均为被动遥感的辐射源。阳辐射和地球自身的热辐射均为被动遥感的辐射源。在在6.0m6.0m以上的热红外波段，以地球自身的热辐射为主以上的热红外波段，以地球自身的热辐射为主，地表反射太阳辐射可以忽略。地表反射太阳辐射可以忽略。18三 地球的辐射与地物波谱地球的辐射与地物波谱o2 地物波谱地物波谱：地物波谱：地物的电磁波响应特性随电磁波长改变地物的电磁波响应特性随电磁波长改变而变化的规律，称为地表物体波谱，简称地物波谱。而变化的规律，称为地表物体波
13、谱，简称地物波谱。地物波谱特性地物波谱特性是电磁辐射与地物相互作用的一种表是电磁辐射与地物相互作用的一种表现。现。地物波谱的作用：地物波谱的作用：不同类型的地物，其电磁波响应不同类型的地物，其电磁波响应的特性不同，因此地物波谱特征是遥感识别地物的的特性不同，因此地物波谱特征是遥感识别地物的基础。基础。不同电磁波段中地物波谱特性不同电磁波段中地物波谱特性地物反射波谱：地物反射波谱：是研究可见光至近红外波段上地物是研究可见光至近红外波段上地物反射率随波长的变化规律。反射率随波长的变化规律。可见光和近红外波段：可见光和近红外波段：主要表现地物主要表现地物反射作用和地物的吸收作用。反射作用和地物的吸收
14、作用。热红外波段：热红外波段：主要表现地物热辐射作主要表现地物热辐射作用。用。微波波段：微波波段：主动遥感利用地物后向散主动遥感利用地物后向散射；被动遥感利用地物微波辐射。射；被动遥感利用地物微波辐射。地物反射波谱曲线除随不同地物地物反射波谱曲线除随不同地物(反反射率射率)不同外，同种地物在不同内部不同外，同种地物在不同内部结构和外部条件下形态表现结构和外部条件下形态表现(反射率反射率)也不同。一般来说，地物反射率随也不同。一般来说，地物反射率随波长变化有规律可循，从而为遥感波长变化有规律可循，从而为遥感影像的判读提供依据。影像的判读提供依据。19植物的光谱曲线叶绿素吸收叶绿素吸收水分吸收水分
15、吸收叶片色素叶片色素细胞结构细胞结构含水量含水量控制叶片控制叶片反射率的反射率的主要因素主要因素主要吸收主要吸收带带20不同植被类型的光谱曲线比较识别植被的识别植被的理想波段理想波段21三种不同类型土壤在干燥环境下的光谱曲线22叶绿素含量不同时水体的光谱曲线23岩石的光谱曲线24 1.3 1.3 遥感图像类型与特性遥感图像类型与特性o遥感平台遥感平台o摄影成像摄影成像o扫描成像扫描成像o微波遥感与成像微波遥感与成像o遥感图像的特征遥感图像的特征25一一 遥感平台遥感平台o遥感平台（遥感平台（platformplatform）是搭载传感器的工具。）是搭载传感器的工具。o根据运载工具的类型划分：根
16、据运载工具的类型划分：航天平台航天平台150km150km以上以上 卫星、卫星、宇宙飞船宇宙飞船航空平台航空平台100m100m至十余至十余 公里飞机、飞船、气球公里飞机、飞船、气球地面平台地面平台0 050m 50m 车、船、塔车、船、塔26二 摄影成像摄影成像o定义：定义：摄影是通过成像设备获取物体影像的技术。摄影是通过成像设备获取物体影像的技术。o传统摄影：传统摄影：依靠光学镜头及放置在焦平面的感光胶片来依靠光学镜头及放置在焦平面的感光胶片来记录物体影像。记录物体影像。o数字摄影：数字摄影：通过放置在焦平面的光敏元件，经光通过放置在焦平面的光敏元件，经光/电转电转换，以数字信号来记录物体
17、的影像。换，以数字信号来记录物体的影像。o依据探测波段依据探测波段近紫外摄影近紫外摄影可见光摄影可见光摄影红外摄影红外摄影多光谱摄影多光谱摄影27三三 扫描成像扫描成像o扫描成像是依靠探测元件和扫描镜对目标地物扫描成像是依靠探测元件和扫描镜对目标地物以瞬时视场为单位进行的逐点、逐行取样，以以瞬时视场为单位进行的逐点、逐行取样，以得到目标地物电磁辐射特性信息，形成一定谱得到目标地物电磁辐射特性信息，形成一定谱段的图像。段的图像。光机扫描仪光机扫描仪固体扫描仪固体扫描仪成像光谱仪成像光谱仪28摄影成像摄影成像扫描成像扫描成像波谱范围波谱范围光谱分辨率光谱分辨率多光谱获取方式多光谱获取方式数据记录方
18、式数据记录方式投影方式投影方式可见光可见光+近红外近红外相对低相对低多个镜头多个镜头胶片、数字胶片、数字中心投影中心投影可见光可见光+近红外近红外+热红外热红外相对高相对高单镜头，分光单镜头，分光数字数字多中心投影多中心投影29四四 微波遥感与成像微波遥感与成像o侧视雷达侧视雷达是在飞机或卫星平台上由传感器向与飞行方向垂直的是在飞机或卫星平台上由传感器向与飞行方向垂直的侧面，发射一个窄的波束，覆盖地面上这一侧面的一个条带，侧面，发射一个窄的波束，覆盖地面上这一侧面的一个条带，然后接收在这一条带上地物的反射波，从而形成一个图像带。然后接收在这一条带上地物的反射波，从而形成一个图像带。随着飞行器前
19、进，不断地发射这种脉冲波束，又不断地接收回随着飞行器前进，不断地发射这种脉冲波束，又不断地接收回波，从而形成一幅一幅的雷达图像。波，从而形成一幅一幅的雷达图像。o雷达成像的基本条件：雷达成像的基本条件：雷达发射的波束照在目标不同部位时，雷达发射的波束照在目标不同部位时，要有时间先后差异，这样从目标反射的回波也同时出现时间差，要有时间先后差异，这样从目标反射的回波也同时出现时间差，才有可能区分目标的不同部位。才有可能区分目标的不同部位。o合成孔径雷达合成孔径雷达与侧视雷达类似，也是在飞机或卫星平台上由传与侧视雷达类似，也是在飞机或卫星平台上由传感器向与飞行方向垂直的侧面发射信号。所不同的是将发射
20、和感器向与飞行方向垂直的侧面发射信号。所不同的是将发射和接收天线分成许多小单元，每一单元发射和接收信号的时刻不接收天线分成许多小单元，每一单元发射和接收信号的时刻不同。由于天线位置不同，记录的回波相位和强度都不同。同。由于天线位置不同，记录的回波相位和强度都不同。目的：目的：提高图象在飞行方向的分辨率。提高图象在飞行方向的分辨率。30地形畸变地形畸变：山上面向雷达的一面在图象上被压缩，这一部分往往表现山上面向雷达的一面在图象上被压缩，这一部分往往表现为较高的亮度；为较高的亮度；坡底的收缩度比坡顶大；山坡的坡度越大，收缩量越大。坡底的收缩度比坡顶大；山坡的坡度越大，收缩量越大。当面向雷达的山坡很
21、陡时，出现山顶比山底更接近雷达的当面向雷达的山坡很陡时，出现山顶比山底更接近雷达的情况，因此，在图象的距离方向，山顶和山底的相对位置情况，因此，在图象的距离方向，山顶和山底的相对位置颠倒；颠倒；收缩度：坡顶的收缩度比坡底大收缩度：坡顶的收缩度比坡底大当后坡坡度较大，雷达波束不能到达后坡坡面时，没有回当后坡坡度较大，雷达波束不能到达后坡坡面时，没有回波信号产生，图象上出现暗区波信号产生，图象上出现暗区31五五 遥感图像的特征遥感图像的特征o将遥感图像归纳为三个方面的特征：将遥感图像归纳为三个方面的特征：n几何特征几何特征n物理特征物理特征n时间特征时间特征o这个三方面特征的表现参数即为：这个三方
22、面特征的表现参数即为：n空间分辨率空间分辨率n光谱分辨率光谱分辨率n辐射分辨率辐射分辨率n时间分辨辨时间分辨辨32 1.4 1.4 遥感图像处理遥感图像处理o遥感技术的目的是为了获得地物的几何属性和物理属性。遥感技术的目的是为了获得地物的几何属性和物理属性。原始的遥感图像并不能提供实现这个目的所需的准确而完原始的遥感图像并不能提供实现这个目的所需的准确而完备的条件。备的条件。o为了实现这个目的，原始遥感影像需要经过图像处理，来为了实现这个目的，原始遥感影像需要经过图像处理，来消除成像过程中的误差，改善图像质量。消除成像过程中的误差，改善图像质量。o 遥感图像处理包括以下几个阶段：遥感图像处理包
23、括以下几个阶段：图像的校正（预处理）图像的校正（预处理）图像的变换图像的变换图像的增强图像的增强图像的分类图像的分类o所采用的手段有：所采用的手段有：光学图像处理光学图像处理和和数字图像处理数字图像处理两种方法。两种方法。331 1 光学图像处理光学图像处理较成熟的遥感图像光学处理方法：较成熟的遥感图像光学处理方法：n假彩色密度分割假彩色密度分割n光学彩色合成光学彩色合成 n反差增强反差增强 n边缘增强边缘增强 n比值增强比值增强 n影像相减影像相减(差值法差值法)n光学变换处理光学变换处理(低、高、带通滤波、方向滤波低、高、带通滤波、方向滤波)为了突出某一密度等级的色调（或相应地为了突出某一
24、密度等级的色调（或相应地物），将图像（或影像）的色调密度分划成物），将图像（或影像）的色调密度分划成若干个等级，并用不同的颜色分别表示这不若干个等级，并用不同的颜色分别表示这不同的密度等级，得到一幅彩色的等密度分割同的密度等级，得到一幅彩色的等密度分割图像。这一技术过程就叫作密度分割处理，图像。这一技术过程就叫作密度分割处理，或简称密度分割。或简称密度分割。适用于轮廓简单、背景单调、目标特性与适用于轮廓简单、背景单调、目标特性与影像密度有明确对应关系的地面景况的分影像密度有明确对应关系的地面景况的分析判读。析判读。342 2 数字图像处理数字图像处理o图像校正图像校正o图像增强图像增强o图像融
25、合图像融合o图像分类图像分类依据数字图像数学特征，构造各依据数字图像数学特征，构造各种数学模型和相应的算法，由计种数学模型和相应的算法，由计算机进行运算算机进行运算(矩阵变换矩阵变换)处理，处理，获得更有利于实际应用的输出图获得更有利于实际应用的输出图像及有关数据和资料。像及有关数据和资料。35图像校正图像校正 图像校正是指从具有畸变的图像中消除畸变的处理过图像校正是指从具有畸变的图像中消除畸变的处理过程。消除辐射量失真的叫辐射量校正，消除几何畸变的程。消除辐射量失真的叫辐射量校正，消除几何畸变的叫几何校正。叫几何校正。p 辐射校正辐射校正由于由于传感器传感器响应特性和响应特性和大气大气的吸收
26、、散射及其的吸收、散射及其它随机因素影响，导致图象模糊失真，造成图它随机因素影响，导致图象模糊失真，造成图象分辨率和对比度相对下降。这些都需要通过象分辨率和对比度相对下降。这些都需要通过辐射校正辐射校正复原复原.包括：包括：系统辐射校正、大气校正系统辐射校正、大气校正p几何校正几何校正校正遥感图像成像过程中所造成的各种几何畸变称为几何校正。校正遥感图像成像过程中所造成的各种几何畸变称为几何校正。公式法、回归公式法、回归分析法、直方分析法、直方图校正法图校正法最邻近法、双线最邻近法、双线性内插法、三次性内插法、三次卷积法卷积法36图像增强图像增强o数字图像增强的目的是提高图像质量和突出所数字图像
27、增强的目的是提高图像质量和突出所需要信息，有利于分析判断。需要信息，有利于分析判断。o数字图像增强处理方法主要有：对比度增强数字图像增强处理方法主要有：对比度增强（灰度增强）、空间滤波、彩色变换、图像运（灰度增强）、空间滤波、彩色变换、图像运算和多光谱变换。算和多光谱变换。37n图像的差值运算有利于目标与背景反差较小的信息提取:如冰雪覆盖区黄土高原区的界线特征海岸带的潮汐线等。n差值运算还常用于研究同一地区不同时相的动态变化。如监测森林火灾发生前后的变化和计算过火面积；监测水灾发生前后的水域变化和计算受灾面积及损失；监测城市在不同年份的扩展情况及计算侵占农田的比例等。n有时为了突出边缘，也用差
28、值法将两幅图像的行、列各移一位，再与原图像相减，也可起到几何增强的作用。n两幅同样行、列数的图像，对应像元的亮度值相除(除数不为0)就是比值运算n比值运算可以检测波段的斜率信息并加以扩展，以突出不同波段间地物光谱的差异，提高对比度。38n比值运算能扩大不同地物之间的微小差异，有利于岩石、土壤等波谱差异不太明显的地物的区分，也可用于植被类型和分布的研究z铁帽与植被在单波段上不易区分，而通过MSS5/4、MSS7/5的二维比值分析，就可以明显区分出来。z突出遥感影像中的植被特征、提取植被类别或估算植被生物量，常用的植被指数（NDVI）算法近红外波段红波段 (近红外-红)(近红外+红)B5/B4B7
29、/B5铁帽植被39n消除或减弱地形等环境因素的影响，如阳坡和阴坡有不同的亮度，但在MSS4/5，TM1/2上，比值却非常接近，因此消除了地形的影响。TM1 TM2 TM1/240n提取与找矿有关的专题信息。T如黏土矿物在2.2um附近存在强吸收，故TM7为低亮度，而TM5仍为高亮度，因此TM5/7常提取与黏土化有关的矿化信息；T0.48um是铁离子电荷强烈吸收的位置，故TM5/1可提取与铁矿物有关的信息。TM3/1 可用于识别褐铁矿；可用于识别褐铁矿；TM5/4 用于识别有无植被覆盖的土壤和岩石；用于识别有无植被覆盖的土壤和岩石；TM2/3 用于识别绿色植被；用于识别绿色植被；TM5/7 用于
30、识别含羟基矿物、水和硫酸盐和碳酸盐；用于识别含羟基矿物、水和硫酸盐和碳酸盐；TM7/4 用于区分不同种类的特征矿物；用于区分不同种类的特征矿物；TM3/4 用于识别植被的常用波段。用于识别植被的常用波段。41图像融合图像融合o遥感图像信息融合遥感图像信息融合(Fusion)(Fusion)是将多源遥感数据在统一是将多源遥感数据在统一的地理坐标系中，采用一定的算法生成一组新的信息或的地理坐标系中，采用一定的算法生成一组新的信息或合成图像的过程。合成图像的过程。o不同的遥感数据具有不同的空间分辨率、波谱分辨率和不同的遥感数据具有不同的空间分辨率、波谱分辨率和时相分辨率，如果能将它们各自的优势综合起
31、来，可以时相分辨率，如果能将它们各自的优势综合起来，可以弥补单一图像上信息的不足，这样不仅扩大了各自信息弥补单一图像上信息的不足，这样不仅扩大了各自信息的应用范围，而且大大提高了遥感影像分析的精度。的应用范围，而且大大提高了遥感影像分析的精度。o数据融合的概念始于数据融合的概念始于7070年代。数据融合的发展在一定程年代。数据融合的发展在一定程度上解决了多种数据源综合分析的问题度上解决了多种数据源综合分析的问题J不同传感器的遥感数据复合不同传感器的遥感数据复合 J遥感影像与非遥感数据融合遥感影像与非遥感数据融合42图像分类图像分类图象分类过程的图象分类过程的总目标是总目标是，将图，将图象中所有的像元象中所有的像元自动地进行土地自动地进行土地覆盖类型或土地覆盖类型或土地覆盖专题等方面覆盖专题等方面的分类。的分类。原理原理分类方法分类方法43
