CN104881850B - 一种多时相遥感影像的厚云自动去除方法 - Google Patents
一种多时相遥感影像的厚云自动去除方法 Download PDFInfo
- Publication number
- CN104881850B CN104881850B CN201510274174.1A CN201510274174A CN104881850B CN 104881850 B CN104881850 B CN 104881850B CN 201510274174 A CN201510274174 A CN 201510274174A CN 104881850 B CN104881850 B CN 104881850B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- msub
- spissatus
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种多时相遥感影像的厚云自动去除方法,由以下步骤组成:步骤1:采集遥感卫星厚云影像和T幅与其同区域的多时相影像;步骤2:厚云区域检测,得到厚云区域指示模板 ;步骤3:自动选择参考影像:自动选择一幅多时相影像作为参考影像;步骤4:采用泊松方程修复方法去除厚云,得到初步去云结果;步骤5:将所述初步去云结果和参考影像纳入变分模型,再次去除厚云,得到最终的去云结果。本发明能通过厚云影像和多幅多时相影像在梯度值之间的均方根误差()确定参考影像。它无需人工交互,自动地去除厚云及其阴影,它兼顾了原始影像的像素亮度与参考影像的梯度信息,对像素值有较好保真性。
Description
技术领域
本发明涉及一种遥感影像的厚云自动去除方法,尤其是一种多时相遥感影像的厚云自动去除方法,属于遥感影像预处理技术领域。
背景技术
遥感影像中,云覆盖是造成遥感数据缺乏的重要因素之一。大量的遥感影像由于云覆盖的干扰,降低了感兴趣信息的清晰度,从而降低了利用率。有效地减少或去除云的影响,是提高遥感数据利用率的一个重要途径,也是遥感影像预处理中的一个重要问题。遥感影像厚云去除可以恢复不完整的影像,增加遥感影像数据来源、降低数据成本,为促进遥感影像军民两用提供技术支持,以获取较大的经济效益。
由于厚云区域缺乏可利用的信息,直接去除比较困难。目前国内外去除厚云的方法主要有单幅影像的厚云去除方法和多时相影像的厚云去除方法。单幅影像的厚云去除方法主要运用影像修复和合成技术。它虽然能够产生合理的视觉效果,但是无法保证影像像素值的真实性。多时相影像的厚云去除方法是利用同一区域不同时间的影像,根据影像间的时间和空间相关性来去除厚云。但是现有多时相影像厚云去除方法缺乏有效稳定的选择参考影像的标准。同时,这些方法需要很多人工干预,难以同时去除厚云和云的阴影。所以这些方法无法得到实际应用。
发明内容
本发明的目的在于针对已有技术存在的缺陷,提供一种能够有效的选择参考影像,自动去除厚云及其阴影的多时相遥感影像的厚云自动去除方法。
为实现上述目的,本发明采用下述技术方案:
一种多时相遥感影像的厚云自动去除方法,由以下步骤组成:
步骤1:采集遥感卫星厚云影像Iraw(i,j,k),1≤k≤K和T幅与其同区域的多时相影像Iraw,t(i,j,k),1≤i≤M,1≤j≤N,1≤t≤T,1≤k≤K,;其中i和j分别为图像中像素的行坐标和列坐标,k为图像的波段编号,M、N和K分别代表厚云影像的行数、列数和波段数。
步骤2:厚云区域检测,得到厚云区域指示模板mask(i,j);
步骤3:自动选择参考影像:自动选择一幅多时相影像作为参考影像Iref(i,j,k);
步骤4:采用泊松方程修复方法去除厚云,得到初步去云结果I0(i,j,k);
步骤5:将所述初步去云结果I0(i,j,k)和参考影像Iref(i,j,k)纳入变分模型,再次去除厚云,得到最终的去云结果I'0(i,j,k)。
上述步骤2由以下具体步骤组成:
步骤2-1:将厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)与阈值TC比较,计算厚云区域初始指示模板mask0(i,j):
式中mask0(i,j)为1表示对应像素属于厚云区域,mask0(i,j)为0表示对应像素属于无云区域;
步骤2-2:对所述厚云区域初始指示模板mask0(i,j)进行形态学处理,去除孤立的斑点并填补细小的空洞,得到厚云区域指示模板mask(i,j)。
上述步骤2-1中阈值TC的计算方法为:
TC=m+λσ (2)
其中,m为所述厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素均值:
σ为厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素标准方差:
λ为权重系数:
上述步骤3由以下具体步骤组成:
步骤3-1:计算厚云影像的各波段的图像Iraw(i,j,k),1≤k≤K在x方向的梯度值gx,0(x,y,k)和y方向的梯度值gy,0(x,y,k):
步骤3-2:计算所述各多时相影像Iraw,t(i,j,k),1≤t≤T在x方向的梯度值gx,t(x,y,k)和y方
向的梯度值gy,t(x,y,k):
步骤3-3:计算所述厚云影像Iraw(i,j,k),1≤k≤K与各多时相影像Iraw,t(i,j,k),1≤t≤T在无
云区域对应位置梯度值之间的均方根误差RMSEt:
步骤3-4:选择RMSEt最小的多时相影像Iram,t(i,j,k),作为所述厚云影像Iraw(i,j,k)的去云参考影像Iref(i,j,k)。
上述步骤4中所述初步去云结果I0(i,j,k)中无云区域像素值与所述厚云影像Iraw(i,j,k)的无云区域像素值相等,其厚云区域图像的像素值由下式计算得到:
其中,Ω代表厚云区域内部;代表厚云区域边界,表示梯度算子,I1表示I0在区域Ω的像素。
上述步骤4中厚云区域图像的像素值可以根据以下递推公式计算得到:
其中
式中pi和pj分别代表初步去云结果I0(i,j,k)中第i个和第j个厚云区域像素;N(pi)代表所述像素pi的四邻域像素集;
式中I0(qi)表示初步去云结果I0(i,j,k)中像素qi的像素值,Iref(pi)和Iref(qi)分别表示所述去云参考影像Iref(i,j,k)中像素pi和qi的像素值,像素qi为N(pi)与Ω交集内的第i个像素;C为A的Choleskey分解,A=CCT,其中C为一个下三角矩阵,其下三角部分具有和A完全相同的结构。
上述步骤5中的变分模型代价函数为:
其中β为用于权衡两个约束项贡献度参数,β≥0。
上述步骤5中的变分模型代价函数最小值的求解方法为采用下述迭代公式计算:
其中Δt是迭代步长,α是正则化参数,αn的计算公式为:
式中var_n为设定的标准参数。
本发明的有益效果在于:
1、本发明能通过厚云影像和多幅多时相影像在梯度值之间的均方根误差(RMSE)确定参考影像。它无需人工交互,自动地去除厚云及其阴影
2、本发明兼顾了原始影像的像素亮度与参考影像的梯度信息,对像素值有较好保真性。
附图说明
图1为本发明的流程图;
图2为本发明步骤2的流程图;
图3为本发明步骤3的流程图。
图4为本发明的遥感卫星采集的厚云影像;
图5为本发明的遥感卫星采集的参考影像;
图6为本发明检测出的厚云区域的指示模板;
图7为本发明的初步去云结果;
图8为本发明的最终去云结果。
具体实施方式
实施例1:
如图1所示,一种多时相遥感影像的厚云自动去除方法,由以下步骤组成:
步骤1:采集遥感卫星厚云影像Iraw(i,j,k),1≤k≤K和T幅与其同区域的多时相影像Iraw,t(i,j,k),1≤i≤M,1≤j≤N,1≤t≤T,1≤k≤K,;其中i和j分别为图像中像素的行坐标和列坐标,k为图像的波段编号,M、N和K分别代表厚云影像的行数、列数和波段数;;在本实施例中M、N和K分别为400、400和7;
步骤2:厚云区域检测,得到厚云区域指示模板mask(i,j);
步骤3:自动选择参考影像:自动选择一幅多时相影像作为参考影像Iref(i,j,k);
步骤4:采用泊松方程修复方法去除厚云,得到初步去云结果I0(i,j,k);
步骤5:将所述初步去云结果I0(i,j,k)和参考影像Iref(i,j,k)纳入变分模型,再次去除厚云,得到最终的去云结果I'0(i,j,k);
如图2所示,上述步骤2由以下具体步骤组成:
步骤2-1:将厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)与阈值TC比较,计算厚云区域初始指示模板mask0(i,j):
式中mask0(i,j)为1表示对应像素属于厚云区域,mask0(i,j)为0表示对应像素属于无云区域;
步骤2-2:对所述厚云区域初始指示模板mask0(i,j)进行形态学处理,去除孤立的斑点并填补细小的空洞,得到厚云区域指示模板mask(i,j)。
上述步骤2-1中阈值TC的计算方法为:
TC=m+λσ (2)
其中,m为所述厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素均值:
σ为厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素标准方差:
λ为权重系数:
如图3所示,上述步骤3由以下具体步骤组成:
步骤3-1:计算厚云影像的各波段的图像Iraw(i,j,k),1≤k≤K在x方向的梯度值gx,0(x,y,k)和y方向的梯度值gy,0(x,y,k):
步骤3-2:计算所述各多时相影像Iraw,t(i,j,k),1≤t≤T在x方向的梯度值gx,t(x,y,k)和y方
向的梯度值gy,t(x,y,k):
步骤3-3:计算所述各厚云影像Iraw(i,j,k),1≤k≤K与各多时相影像Iraw,t(i,j,k),1≤t≤T在
无云区域对应位置梯度值之间的均方根误差RMSEt:
步骤3-4:选择RMSEt最小的多时相影像Iram,t(i,j,k),作为所述厚云影像Iram(i,j,k)的去云参考影像Iref(i,j,k)。
上述步骤4由以下具体步骤组成:
步骤4-1:所述初步去云结果I0(i,j,k)中无云区域像素值与所述厚云影像Iraw(i,j,k)的无云区域像素值相等,其厚云区域图像的像素值由下式计算得到:
其中,Ω代表厚云区域内部;代表厚云区域边界,表示梯度算子,I1表示I0在区域Ω的像素。
上述步骤4-1中厚云区域图像的像素值可以根据以下递推公式计算得到:
其中
式中pi和pj分别代表初步去云结果I0(i,j,k)中第i个和第j个厚云区域像素;N(pi)代表所述像素pi的四邻域像素集;
式中I0(qi)表示初步去云结果I0(i,j,k)中像素qi的像素值,Iref(pi)Iref(qi)分别表示所述去云参考影像Iref(i,j,k)中像素pi和qi的像素值,像素qi为N(pi)与Ω交集内的第i个像素;C为A的Choleskey分解,A=CCT,其中C为一个下三角矩阵,其下三角部分具有和A完全相同的结构。
上述步骤5中的变分模型代价函数为:
其中β为用于权衡两个约束项贡献度参数,β≥0。
上述步骤5中的变分模型代价函数最小值的求解方法为采用下述迭代公式计算:
其中Δt是迭代步长,α是正则化参数,αn的计算公式为:
式中var_n为设定的标准参数。在本实施例中,var_n为519.84。
Claims (7)
1.一种多时相遥感影像的厚云自动去除方法,其特征在于:由以下步骤组成:
步骤1:采集遥感卫星厚云影像Iraw(i,j,k),1≤k≤K和T幅与其同区域的多时相影像Iraw,t(i,j,k),1≤i≤M,1≤j≤N,1≤t≤T,1≤k≤K,其中i和j分别为图像中像素的行坐标和列坐标,k为图像的波段编号,M、N和K分别代表厚云影像的行数、列数和波段数;
步骤2:厚云区域检测,得到厚云区域指示模板mask(i,j);
步骤3:自动选择参考影像:自动选择一幅多时相影像作为参考影像Iref(i,j,k),具体步骤为:
步骤3-1:计算厚云影像的各波段的图像Iraw(i,j,k),1≤k≤K在x方向的梯度值gx,0(x,y,k)和y方向的梯度值gy,0(x,y,k):
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<mi>K</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3-2:计算各所述多时相影像Iraw,t(i,j,k),1≤t≤T在x方向的梯度值gx,t(x,y,k)和y方向的梯度值gy,t(x,y,k):
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3-3:计算各所述厚云影像Iraw(i,j,k),1≤k≤K与各多时相影像Iraw,t(i,j,k),1≤t≤T在无云区域对应位置梯度值之间的均方根误差RMSEt:
<mrow>
<msub>
<mi>RMSE</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>M</mi>
<mi>N</mi>
<mi>K</mi>
</mrow>
</mfrac>
</msqrt>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>s</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3-4:选择RMSEt最小的多时相影像Iram,t(i,j,k),作为所述厚云影像Iraw(i,j,k)的去云参考影像Iref(i,j,k);
步骤4:采用泊松方程修复方法去除厚云,得到初步去云结果I0(i,j,k);
步骤5:将所述初步去云结果I0(i,j,k)和参考影像Iref(i,j,k)纳入变分模型,再次去除厚云,得到最终的去云结果I'0(i,j,k)。
2.根据权利要求1所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤2由以下具体步骤组成:
步骤2-1:将厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)与阈值TC比较,计算厚云区域初始指示模板mask0(i,j):
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>s</mi>
<mi>k</mi>
<mn>0</mn>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>N</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中mask0(i,j)为1表示对应像素属于厚云区域,mask0(i,j)为0表示对应像素属于无云区域;
步骤2-2:对所述厚云区域初始指示模板mask0(i,j)进行形态学处理,去除孤立的斑点并填补细小的空洞,得到厚云区域指示模板mask(i,j)。
3.根据权利要求2所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤2-1中阈值TC的计算方法为:
TC=m+λσ (2)
其中,m为厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素均值:
<mrow>
<mi>m</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
σ为厚云影像Iraw(i,j,k)的蓝光波段图像I(i,j)的像素标准方差:
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
λ为权重系数:
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>(</mo>
<mi>I</mi>
<mo>)</mo>
<mo>-</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mi>I</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
4.根据权利要求1所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤4中所述初步去云结果I0(i,j,k)中无云区域像素值与所述厚云影像Iraw(i,j,k)的无云区域像素值相等,其厚云区域图像的像素值由下式计算得到:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>arg</mi>
<mi> </mi>
<mi>min</mi>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<munder>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</munder>
<mo>|</mo>
<mo>&dtri;</mo>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>&dtri;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mo>&part;</mo>
<mi>&Omega;</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mo>&part;</mo>
<mi>&Omega;</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>s</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Ω代表厚云区域内部;代表厚云区域边界,表示梯度算子,I1表示I0在区域Ω的像素。
5.根据权利要求4所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤4中厚云区域图像的像素值根据以下递推公式计算得到:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mo><</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>></mo>
</mrow>
<mrow>
<mo><</mo>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>></mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mo><</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>></mo>
</mrow>
<mrow>
<mo><</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>></mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<mi>B</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mi>B</mi>
<mo>-</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
式中pi和pj分别代表初步去云结果I0(i,j,k)中第i个和第j个厚云区域像素;N(pi)代表所述像素pi的四邻域像素集;
<mrow>
<mi>B</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&cap;</mo>
<mo>&part;</mo>
<mi>&Omega;</mi>
</mrow>
</munder>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式中I0(qi)表示初步去云结果I0(i,j,k)中像素qi的像素值,Iref(pi)和Iref(qi)分别表示所述去云参考影像Iref(i,j,k)中像素pi和qi的像素值,像素qi为N(pi)与Ω交集内的第i个像素;C为A的Choleskey分解,A=CCT,其中C为一个下三角矩阵,其下三角部分具有和A完全相同的结构。
6.根据权利要求1所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤5中的变分模型代价函数为:
<mrow>
<mi>J</mi>
<mo>&lsqb;</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<munder>
<mrow>
<mo>&Integral;</mo>
<mo>&Integral;</mo>
</mrow>
<mi>&Gamma;</mi>
</munder>
<mfrac>
<mi>&beta;</mi>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mo>&dtri;</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mo>&dtri;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>x</mi>
<mi>d</mi>
<mi>y</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
其中β为用于权衡两个约束项贡献度参数,β≥0。
7.根据权利要求6所述的多时相遥感影像的厚云自动去除方法,其特征在于:所述步骤5中的变分模型代价函数最小值的求解方法为采用下述迭代公式计算:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mn>0</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mn>0</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mn>0</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>-</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&alpha;</mi>
<mi>n</mi>
</msup>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Δt是迭代步长,表示梯度算子,α是正则化参数,αn的计算公式为:
<mrow>
<msup>
<mi>&alpha;</mi>
<mi>n</mi>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&dtri;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>|</mo>
<mo>&dtri;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>I</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>var</mi>
<mo>_</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
式中var_n为设定的标准参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510274174.1A CN104881850B (zh) | 2015-05-26 | 2015-05-26 | 一种多时相遥感影像的厚云自动去除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510274174.1A CN104881850B (zh) | 2015-05-26 | 2015-05-26 | 一种多时相遥感影像的厚云自动去除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104881850A CN104881850A (zh) | 2015-09-02 |
CN104881850B true CN104881850B (zh) | 2017-12-22 |
Family
ID=53949335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510274174.1A Expired - Fee Related CN104881850B (zh) | 2015-05-26 | 2015-05-26 | 一种多时相遥感影像的厚云自动去除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104881850B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780353B (zh) * | 2016-11-10 | 2019-10-22 | 哈尔滨工业大学 | 一种基于时相光谱角度量的多时相云遮挡数据恢复方法 |
CN107230195B (zh) * | 2017-07-12 | 2020-09-18 | 中国科学院遥感与数字地球研究所 | 一种影像处理方法和装置 |
CN108765329B (zh) * | 2018-05-21 | 2020-12-04 | 北京师范大学 | 一种遥感影像的厚云去除方法及系统 |
CN110335208B (zh) * | 2019-06-10 | 2022-06-07 | 武汉大学 | 一种基于逐步校正的高分辨率遥感影像厚云去除方法 |
CN114511786B (zh) * | 2022-04-20 | 2022-07-19 | 中国石油大学(华东) | 融合多时相信息和分通道密集卷积的遥感图像去云方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246545B (zh) * | 2008-02-22 | 2010-07-07 | 华南师范大学 | 一种光学遥感图像去云的泊松方法 |
CN103020939A (zh) * | 2012-12-18 | 2013-04-03 | 武汉大学 | 利用多时相数据去除光学遥感影像大面积厚云的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9977978B2 (en) * | 2011-11-14 | 2018-05-22 | San Diego State University Research Foundation | Image station matching, preprocessing, spatial registration and change detection with multi-temporal remotely-sensed imagery |
-
2015
- 2015-05-26 CN CN201510274174.1A patent/CN104881850B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246545B (zh) * | 2008-02-22 | 2010-07-07 | 华南师范大学 | 一种光学遥感图像去云的泊松方法 |
CN103020939A (zh) * | 2012-12-18 | 2013-04-03 | 武汉大学 | 利用多时相数据去除光学遥感影像大面积厚云的方法 |
Non-Patent Citations (2)
Title |
---|
A Variational Gradient-based Fusion Method for Visible and SWIR Imagery;Huifang Li et al;《Photogrammetric Engineering & Remote Sensing》;20120930;第79卷(第9期);第2页左栏倒数第1-2段,第2页右栏第1-4段,第4页左栏倒数第1-2段,第9页右栏第2段,图5,图6 * |
全色遥感影像上厚云去除的若干关键技术研究;董保根 等;《信息工程大学学报》;20110831;第12卷(第4期);第2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104881850A (zh) | 2015-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104881850B (zh) | 一种多时相遥感影像的厚云自动去除方法 | |
CN104809698A (zh) | 一种基于改进三边滤波的Kinect深度图像修复方法 | |
CN105574827B (zh) | 一种图像去雾的方法、装置 | |
CN106204447A (zh) | 基于总变差分和卷积神经网络的超分辨率重建方法 | |
JP2010187376A5 (zh) | ||
CN103985098B (zh) | 一种证件图像的高光去除方法及系统 | |
CN103699900B (zh) | 卫星影像中建筑物水平矢量轮廓自动批量提取方法 | |
CN103971338B (zh) | 一种基于显著图的可变块图像修复方法 | |
CN104299263B (zh) | 一种基于单幅图像建模云场景的方法 | |
CN105303616B (zh) | 基于单张照片的浮雕建模方法 | |
CN104715451B (zh) | 一种基于颜色及透明度一致优化的图像无缝融合方法 | |
CN104899845A (zh) | 一种基于lαβ空间场景迁移的多曝光图像融合方法 | |
CN107808386A (zh) | 一种基于图像语义分割的海天线检测方法 | |
CN104143186B (zh) | 一种slic超像素分块优化方法 | |
CN105243670A (zh) | 一种稀疏和低秩联合表达的视频前景对象精准提取方法 | |
CN103954970B (zh) | 一种地形要素采集方法 | |
CN102129669B (zh) | 一种航空遥感影像的最小二乘区域网匀色方法 | |
CN103761734A (zh) | 一种时域一致性保持的双目立体视频场景融合方法 | |
CN108197283A (zh) | 利用河道实测水位进行洪水二、三维动态展示的方法 | |
CN110335208B (zh) | 一种基于逐步校正的高分辨率遥感影像厚云去除方法 | |
CN108416803A (zh) | 一种基于深度神经网络的多信息融合的场景深度恢复方法 | |
CN107610092B (zh) | 基于视频流的路面裂缝动态检测方法 | |
CN103578085A (zh) | 基于可变块的图像空洞区域修补方法 | |
CN101984465B (zh) | 一种图像无缝复制的方法 | |
CN105931210A (zh) | 一种高分辨率图像重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171222 Termination date: 20200526 |
|
CF01 | Termination of patent right due to non-payment of annual fee |