【技术实现步骤摘要】
本申请涉及多光谱原始影像边缘增强技术，特别涉及一种影像增强方法及装置。
技术介绍
1、在遥感影像处理中，基于调制传递函数(modulation transfer function,mtf)相关的影像增强技术被广泛用于卫星原始影像的校正中。由于卫星传感器于在轨运行中性能衰减，导致实验室测得传感器的mtf无法满足长时序、高精度的要求，需要对其在轨测量。通常，mtf的在轨测量需要获取具有固定角度特征(如0°、45°或90°)的地面标志物，但在遥感影像中，往往难以找到具有精确固定角度特征的标志物。
技术实现思路
1、本申请提供一种影像增强方法及装置。
2、第一方面，提供一种影像增强方法。该方法包括：
3、获取原始影像中存在的第一直线段；
4、根据所述第一直线段的长度，从所述第一直线段中获得第二直线段，并根据所述第二直线段于所述原始影像中的位置，从所述原始影像中确定第一影像块；
5、对所述第一影像块旋转和/或平移，获得第二影像块，并根据所述第二影像块的像素坐标系，获得所述第一影像块中各个像素各自对应在所述像素坐标系中的参考坐标；
6、根据所述第一影像块中各个像素的参考坐标，以及所述第一影像块中各个像素的灰度值，获得所述第二影像块中各个像素的灰度值；其中，所述第二影像块中各个像素的灰度值为亚像素级灰度值；
7、根据所述第二影像块的灰度均值和纹理均匀性值，获得地面标志物，并计算所述地面标志物所对应的mtf函数；
9、在一种可能的实施方式中，所述根据所述第二影像块的灰度均值和纹理均匀性值，获得地面标志物，包括：从所述第二影像块中，获得符合预设区域灰度差异阈值范围以及预设纹理均匀性阈值范围的地面标志物；其中，所述纹理均匀性值由第二影像块的idm值表示。
10、在一种可能的实施方式中，所述计算所述地面标志物所对应的mtf函数，包括：
11、对所述地面标志物中的任意第i行像素序列，根据所述第i行像素序列中各个像素的灰度值，计算所述第i行像素序列中第二边缘点的位置；
12、根据所述第i行像素序列中第二边缘点的位置，获得所述第i行像素序列的第二灰度分布函数；
13、对所述地面标志物中各行像素序列的第二灰度分布函数进行累加求平均，获得所述地面标志物的第三灰度分布函数；
14、根据所述地面标志物的第三灰度分布函数，计算获得所述地面标志物的mtf函数。
15、在一种可能的实施方式中，所述对所述地面标志物中的任意第i行像素序列，根据第i行像素序列中各个像素的灰度值，计算所述第i行像素序列中第二边缘点的位置，包括：
16、对所述地面标志物中任意第i行像素序列，根据所述第i行像素序列中各个像素的灰度值，获得所述第i行像素序列中的第一边缘点；
17、根据所述第i行像素序列中各个像素的灰度值，计算所述第i行像素序列的灰度阈值；
18、对所述第一边缘点两侧的若干数量的像素的灰度值拟合，获得所述第i行像素序列的第一灰度分布函数；
19、将所述灰度阈值代入所述第一灰度分布函数中，获得所述第i行像素序列的第二边缘点的亚像素级位置。
20、在一种可能的实施方式中，所述根据所述第i行像素序列中各个像素的灰度值，计算所述第i行像素序列的灰度阈值，包括：根据所述第i行像素序列中各个像素的灰度值，通过大津法，计算所述第i行像素序列的灰度阈值。
21、在一种可能的实施方式中，所述对所述地面标志物中各行像素序列的第二灰度分布函数进行累加求平均，获得所述地面标志物的第三灰度分布函数；包括：
22、对所述地面标志物的任意第i行像素序列的灰度阈值进行检测，获得所述地面标志物的第二像素序列；
23、对所述地面标志物的各行第二像素序列的第二灰度分布函数进行累加求平均，获得所述地面标志物的第三灰度分布函数。
24、在一种可能的实施方式中，对所述地面标志物中各行像素序列的第二灰度分布函数进行累加求平均，获得所述地面标志物的第三灰度分布函数；包括：
25、根据随机抽样一致性函数，对所述地面标志物的各行像素序列的第二边缘点检测，获得所述地面标志物的第三像素序列；
26、根据所述地面标志物的各行第三像素序列的第二灰度分布函数，获得所述地面标志物的第三灰度分布函数。
27、在一种可能的实施方式中，根据所述地面标志物所对应的mtf函数，计算获得二维mtf矩阵，包括：根据所述地面标志物所对应的mtf函数，通过二维高斯函数模型进行模拟，构建二维mtf矩阵。
28、第二方面，提供一种影像增强装置，其特征在于，包括：
29、线段获取单元，用于获取原始影像中存在的第一直线段；
30、影像确定单元，根据所述第一直线段的长度，从所述第一直线段中获得第二直线段，并根据所述第二直线段于所述原始影像中的位置，从所述原始影像中确定第一影像块；
31、坐标计算单元，用于对所述第一影像块旋转和/或平移，获得第二影像块，并根据所述第二影像块的像素坐标系，获得所述第一影像块中各个像素各自对应在所述像素坐标系中的参考坐标；
32、灰度计算单元，用于根据所述第一影像块中各个像素的参考坐标，以及所述第一影像块中各个像素的灰度值，获得所述第二影像块中各个像素的灰度值；其中，所述第二影像块中各个像素的灰度值为亚像素级灰度值；
33、函数获取单元，用于根据所述第二影像块的灰度均值和纹理均匀性值，获得地面标志物，并计算所述地面标志物所对应的mtf函数；
34、影像增强单元，用于根据所述mtf函数获得二维mtf矩阵，并使用所述二维mtf矩阵对原始影像增强。
35、第三方面，提供一种计算设备，包括存储器和处理器，所述存储器中存储有可执行代码，所述处理器执行所述可执行代码时，实现第一方面中任一项所述的方法。
36、本申请提供的技术方案中，通过根据原始影像中存在的任意角度第一直线段的长度筛选，获取符合长度范围的第二直线段，并于原始影像中确认第二直线段所对应的影像块；对第二直线段所对应的影像块进行平移和/或旋转操作；并根据影像块的灰度值与纹理均匀度，从影像块中获得符合能够计算mtf函数特征的地面标志物，从而解决了只能根据固定角度特征的地面标志物计算mtf函数的问题。
本文档来自技高网...【技术保护点】
1.一种影像增强方法，其特征在于，包括：
2.根据权利要求1所述的方法，其中，所述根据所述第二影像块的灰度均值和纹理均匀性值，获得地面标志物，包括：从所述第二影像块中，获得符合预设区域灰度差异阈值范围以及预设纹理均匀性阈值范围的地面标志物；其中，所述纹理均匀性值由第二影像块的IDM值表示。
3.根据权利要求1所述的方法，其中，所述计算所述地面标志物所对应的MTF函数，包括：
4.根据权利要求3所述的方法，其中，所述对所述地面标志物中的任意第i行像素序列，根据第i行像素序列中各个像素的灰度值，计算所述第i行像素序列中第二边缘点的位置，包括：
5.根据权利要求4所述的方法，其中，所述根据所述第i行像素序列中各个像素的灰度值，计算所述第i行像素序列的灰度阈值，包括：根据所述第i行像素序列中各个像素的灰度值，通过大津法，计算所述第i行像素序列的灰度阈值。
6.根据权利要求3所述的方法，其中，所述对所述地面标志物中各行像素序列的第二灰度分布函数进行累加求平均，获得所述地面标志物的第三灰度分布函数；包括：
7.根据权利要
8.根据权利要求1-7中任一项所述的方法，其中，根据所述地面标志物所对应的MTF函数，计算获得二维MTF矩阵，包括：根据所述地面标志物所对应的MTF函数，通过二维高斯函数模型进行模拟，构建二维MTF矩阵。
9.一种影像增强装置装置，其特征在于，包括：
10.一种计算设备，包括存储器和处理器，所述存储器中存储有可执行代码，所述处理器执行所述可执行代码时，实现权利要求1-8中任一项所述的方法。
...【技术特征摘要】
1.一种影像增强方法，其特征在于，包括：
2.根据权利要求1所述的方法，其中，所述根据所述第二影像块的灰度均值和纹理均匀性值，获得地面标志物，包括：从所述第二影像块中，获得符合预设区域灰度差异阈值范围以及预设纹理均匀性阈值范围的地面标志物；其中，所述纹理均匀性值由第二影像块的idm值表示。
3.根据权利要求1所述的方法，其中，所述计算所述地面标志物所对应的mtf函数，包括：
4.根据权利要求3所述的方法，其中，所述对所述地面标志物中的任意第i行像素序列，根据第i行像素序列中各个像素的灰度值，计算所述第i行像素序列中第二边缘点的位置，包括：
5.根据权利要求4所述的方法，其中，所述根据所述第i行像素序列中各个像素的灰度值，计算所述第i行像素序列的灰度阈值，包括：根据所述第i行像素序列中各个像素的灰度值，通过大津法，计算所述...
【专利技术属性】
技术研发人员：刘永坤，龙腾飞，焦伟利，杜意鸿，王桂周，何国金，
申请(专利权)人：中国科学院空天信息创新研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。