[发明专利]Landsat TM遥感影像数据除云方法及系统

权利要求书

1.一种Landsat TM遥感影像数据除云方法，其特征在于，该方法包括如下步骤：

a.输入Landsat TM遥感影像数据；

b.采用基于面向对象的方法，对上述Landsat TM遥感影像数据进行多尺度分割；

c.对上述多尺度分割后的Landsat TM遥感影像数据，根据厚云的光谱特点，建立运算特征“ThickC”；

d.利用阈值分类法，将满足阈值条件的运算特征“ThickC”的对象由“unclassified”分类至“Thick Cloud”；

e.对上述多尺度分割后的Landsat TM遥感影像数据，根据薄云的光谱特点和分布特点，建立关系特征“ThinC”；

f.在剩余的“unclassified”对象中，利用双阈值分类法将满足条件的对象“unclassified”分类至“Thin Cloud”；

g.将上述分类至“Thick Cloud”和“Thin Cloud”的对象统一分类至“Cloud”，完成数据除云；

所述的步骤c具体包括：

对于每一个对象ω＝{P₁，P₂，P₃，…，P_n}，P为该对象所包含的像元，n为该对象所包含的像元数量，其运算特征“ThickC”的具体计算公式如下：

其中，

其中，ThickC表示该对象的运算特征“ThickC”的数值；m’表示该影像对象层中的对象数量；表示一个像元的B₂亮度值；表示一个像元的B₄亮度值；表示一个像元的B₇亮度值；C表示一个缩放常数，根据用户需求设定；

对影像对象层中的每一个对象{ω₁，ω₂，ω₃，…，ω_n}都分别计算运算特征“ThickC”。

2.如权利要求1所述的方法，其特征在于，所述Landsat TM遥感影像数据包括7个谱段：B₁蓝绿谱段、B₂绿谱段、B₃红谱段、B₄近红外谱段、B₅近短波红外谱段、B₆热红外谱段、B₇近短波红外谱段。

3.如权利要求2所述的方法，其特征在于，所述的步骤b具体包括：

根据输入的Landsat TM遥感影像数据，建立影像对象层；

选择参与多尺度分割计算的各个波段数据的权重；

设置多尺度分割的尺度参数；

设置多尺度分割的形状参数；

设置多尺度分割的紧致度参数；

在影像对象层上计算多尺度分割的结果，将分割后产生的对象赋类为“unclassified”。

4.如权利要求3所述的方法，其特征在于，所述的步骤e具体包括：

“ThinC”表示若以一个对象为中心对象，则边缘距离这个中心对象边缘的最短距离不超过d个像元的其它对象将作为计算该中心对象关系特征“ThinC”值的元素，其描述如下：

其中，d表示预设距离值；d(ω_i，ω_j)表示对象ω_i与对象ω_j之间的距离；表示参与计算关系特征“ThinC”的所有对象的集合；{ω|d(ω_i，ω_j)≤d，i＝1，2，3，…，n，j＝1，2，3，…，n}表示满足判断条件d(ω_i，ω_j)≤d，i＝1，2，3，…，n，j＝1，2，3，…，n的所有对象的集合；

根据运算特征及其统计算法计算关系特征“ThinC”：

其中，ThinC表示该对象的关系特征值；ThickC表示参与计算关系特征“ThinC”的对象的运算特征“ThickC”值；m表示参与计算关系特征“ThinC”的对象的数量。