基于多波段遥感影像的边缘异常像素裁剪处理方法及系统技术方案

本申请公开了一种基于多波段遥感影像的边缘异常像素裁剪处理方法及系统。首先获取遥感影像文件，并将遥感影像文件所对应的栅格数据转化为三维数组进行存储；将三维数组中各波段影像求和并进行二值化处理得到二值化影像矩阵C；然后根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，得到去除边缘异常部分的二值化影像矩阵F；最后将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像。本申请解决卫星传感器在拍摄过程中产生的边缘抖动引起的异常值导致影像有效范围不精确的问题。题。题。

【技术实现步骤摘要】
基于多波段遥感影像的边缘异常像素裁剪处理方法及系统


[0001]本专利技术涉及图像处理
，特别涉及一种基于多波段遥感影像的边缘异常像素裁剪处理方法及系统。

技术介绍

[0002]多波段遥感影像是指通过遥感技术获取的具有多个波长段（通常超过3个）的图像。这些波长段在电磁光谱中不同，可以包括可见光、红外线、紫外线、雷达等。每个波段都对应着不同的能量和反射率，从而呈现出不同的地物特征。
[0003]在数据处理过程中，现在所有处理过程针对图像内部异常数据进行处理，常见的处理手段有图像平滑，图像增强，邻域填充等方法。然而现有处理过程中会出现拍摄过程所产生的边缘抖动引起的异常值导致影像有效范围不精确的问题。

技术实现思路

[0004]基于此，本申请实施例提供了一种基于多波段遥感影像的边缘异常像素裁剪处理方法及系统，能够解决现有处理过程中存在的问题。
[0005]第一方面，提供了一种基于多波段遥感影像的边缘异常像素裁剪处理方法，该方法包括：
[0006]获取遥感影像文件，并将所述遥感影像文件所对应的栅格数据转化为三维数组进行存储；其中，三维数组包含了遥感影像的行、列和波段信息；
[0007]将三维数组中各波段影像求和得到影像矩阵，并将所述影像矩阵进行二值化处理得到二值化影像矩阵C；其中，具体将矩阵中每个位置所有波段上的像元值相加得到每个位置的一个总像元值；
[0008]根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，得到去除边缘异常部分的二值化影像矩阵F；
>[0009]将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像。
[0010]可选地，所述获取遥感影像文件，并将所述遥感影像文件所对应的栅格数据转化为三维数组进行存储，包括：
[0011]设置遥感影像文件背景值填充为0，有效值范围为1~255。
[0012]可选地，根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，所述方法包括：
[0013]构造一个2*2的全1矩阵E；
[0014]将全1矩阵E在二值化影像矩阵C边缘按照先行i，后列j，最后波段k的顺序进行平移；
[0015]基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值。
[0016]可选地，基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值，包括：
[0017]全1矩阵E依顺序在二值化影像矩阵C上移动当其覆盖二值化影像矩阵C的区域为[1,1;1,1]时，二值化影像矩阵C输出图像的对应位置赋值为1，其余情况赋值为0。
[0018]可选地，将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像，包括：
[0019]二值化影像矩阵F与影像矩阵如果相同位置为真，则保留影像矩阵原值，否则赋值为0；
[0020]依次计算所有波段，得到边缘异常像素裁剪后的目标影像。
[0021]第二方面，提供了一种基于多波段遥感影像的边缘异常像素裁剪处理系统，该系统包括：
[0022]获取模块，用于获取遥感影像文件，并将所述遥感影像文件所对应的栅格数据转化为三维数组进行存储；其中，三维数组包含了遥感影像的行、列和波段信息；
[0023]二值化模块，用于将三维数组中各波段影像求和得到影像矩阵，并将所述影像矩阵进行二值化处理得到二值化影像矩阵C；其中，具体将矩阵中每个位置所有波段上的像元值相加得到每个位置的一个总像元值；
[0024]边缘处理模块，用于根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，得到去除边缘异常部分的二值化影像矩阵F；
[0025]赋值模块，用于将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像。
[0026]可选地，所述获取模块包括：
[0027]设置遥感影像文件背景值填充为0，有效值范围为1~255。
[0028]可选地，所述边缘处理模块包括：
[0029]构造一个2*2的全1矩阵E；
[0030]将全1矩阵E在二值化影像矩阵C边缘按照先行i，后列j，最后波段k的顺序进行平移；
[0031]基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值。
[0032]可选地，所述基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值，包括：
[0033]全1矩阵E依顺序在二值化影像矩阵C上移动当其覆盖二值化影像矩阵C的区域为[1,1;1,1]时，二值化影像矩阵C输出图像的对应位置赋值为1，其余情况赋值为0。
[0034]可选地，所述赋值模块包括：
[0035]二值化影像矩阵F与影像矩阵如果相同位置为真，则保留影像矩阵原值，否则赋值为0；
[0036]依次计算所有波段，得到边缘异常像素裁剪后的目标影像。
[0037]本申请实施例提供的技术方案中首先获取遥感影像文件，并将遥感影像文件所对应的栅格数据转化为三维数组进行存储；将三维数组中各波段影像求和并进行二值化处理得到二值化影像矩阵C；根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，得到去
除边缘异常部分的二值化影像矩阵F；将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像。可以看出，本申请在常规影像生产过程中需要对影像正射成果进行精确的范围确定提供基础，解决卫星传感器在拍摄过程中产生的边缘抖动引起的异常值导致影像有效范围不精确的问题。
附图说明
[0038]为了更清楚地说明本专利技术的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引申获得其它的实施附图。
[0039]图1为本申请实施例提供的一种基于多波段遥感影像的边缘异常像素裁剪处理流程图；
[0040]图2为本申请实施例中进行边缘处理的示意图；
[0041]图3为本申请实施例中进行条件与运算赋值的示意图。
具体实施方式
[0042]为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
[0043]在本专利技术的描述中，术语“包括”、“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包括了一系列步骤或单元的过程、方法、系统、产品或设备不必限于已明确列出的那些步骤或单元，而是还可包含虽然并未明确列出的但对于这些过程、方法、产品或设备固有的其它步骤或单元，或者基于本专利技术构思进一步的优化方案所增加的步骤或单元。
...

【技术保护点】

【技术特征摘要】
1.一种基于多波段遥感影像的边缘异常像素裁剪处理方法，其特征在于，所述方法包括：获取遥感影像文件，并将所述遥感影像文件所对应的栅格数据转化为三维数组进行存储；其中，三维数组包含了遥感影像的行、列和波段信息；将三维数组中各波段影像求和得到影像矩阵，并将所述影像矩阵进行二值化处理得到二值化影像矩阵C；其中，具体将矩阵中每个位置所有波段上的像元值相加得到每个位置的一个总像元值；根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，得到去除边缘异常部分的二值化影像矩阵F；将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像。2.根据权利要求1所述的边缘异常像素裁剪处理方法，其特征在于，所述获取遥感影像文件，并将所述遥感影像文件所对应的栅格数据转化为三维数组进行存储，包括：设置遥感影像文件背景值填充为0，有效值范围为1~255。3.根据权利要求1所述的边缘异常像素裁剪处理方法，其特征在于，根据预先构造的全1矩阵对二值化影像矩阵C进行边缘处理，所述方法包括：构造一个2*2的全1矩阵E；将全1矩阵E在二值化影像矩阵C边缘按照先行i，后列j，最后波段k的顺序进行平移；基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值。4.根据权利要求3所述的边缘异常像素裁剪处理方法，其特征在于，基于全1矩阵E和二值化影像矩阵C对应元素位置的数值对二值化影像矩阵C边缘进行重新赋值，包括：全1矩阵E依顺序在二值化影像矩阵C上移动当其覆盖二值化影像矩阵C的区域为[1,1;1,1]时，二值化影像矩阵C输出图像的对应位置赋值为1，其余情况赋值为0。5.根据权利要求1所述的边缘异常像素裁剪处理方法，其特征在于，将去除边缘异常部分的二值化影像矩阵F与影像矩阵进行条件与运算赋值得到边缘异常像素裁剪后的目标影像，包括：二值化影像矩阵F与影像矩阵如果相同位置为真，则保留影像矩阵原值，否则赋值为0...

【专利技术属性】
技术研发人员：高占平，徐华，
申请(专利权)人：北京爱特拉斯信息科技有限公司，
类型：发明
国别省市：