[发明专利]一种区域无云影像合成方法

权利要求书

1.一种区域无云影像合成方法，其特征在于，包括以下步骤：

步骤1，确定每幅待合成影像I_i及其对应的云检测结果影像I'_i，其中i＝1,2,...,n，所述云检测结果影像I'_i为二值影像；

步骤2，根据每幅待合成影像I_i的左上角和右下角的地理坐标，确定最终合成结果的覆盖范围，并创建一个同样大小的模板影像以及合成影像，同时将像元值初始化为0；

步骤3，根据合成影像左上角第一个像素的地理位置，在待合成影像中选出在该地理位置具有唯一灰度值覆盖的影像，以影像编号和灰度值分别填充到模板影像、合成影像对应位置；

步骤4，对于下一个待填充的位置点，根据地理位置在待合成影像中寻找对该位置有灰度值覆盖的影像序列，并根据对应的云检测结果影像将影像序列分为有云覆盖和无云覆盖两类；从无云覆盖的影像中挑选在待填充点处周边像素最相似的影像，将影像编号填充到模板影像，将灰度值填充到合成影像；

步骤5，重复步骤是4，直至合成影像上每一个像元都被填充，得到合成的目标区域的无云影像和记录了合成影像每个像元来源的模板影像；

步骤6，设定羽化宽度，根据模板影像记录的信息对合成后的无云影像，以羽化宽度矩形范围内各影像编号出现的次数为权重对每一个像素做羽化；

所述步骤6的实现方式如下，

设拼接缝羽化的宽度为w，从行列号为(w/2，w/2)开始到(Rows-w/2,Cols-w/2)结束，对每一个像素做羽化处理，其中Rows、Cols为模板影像或者合成影像的行数和列数；

步骤6.1，以待羽化点(i，j)为中心，以羽化宽度为边长确定矩形区域，统计出现的像元来源的影像编号i₁、i₂、…、i_q及其对应的像素个数n_i1、n_i2、….、n_iq；

步骤6.2，以步骤6.1中的影像编号确定影像序列，根据在待羽化点是否有数据覆盖以及是否为云层剔除部分影像，保留在该点无云覆盖的有效影像Ij₁、Ij₂、...、Ij_r，其中j₁、j₂、...、j_r为影像编号i₁、i₂、…、i_q的子集；

步骤6.3，根据有效影像各自在羽化矩形区域内出现的像素个数n_j1、n_j2、….、n_jr，设置权重加权得到羽化后的像素值，加权羽化计算方式为：

其中，H为创建的合成影像在(i，j)处羽化后的灰度值，I_j1为影像Ij₁在对应地理位置的灰度值；

步骤6.4，对合成影像每一个波段都重复步骤6.3的方式处理；

步骤7，重复步骤6，对每一个像素做拼接缝羽化，得到经过拼接缝消除的目标区域的无云影像。