CN102129669A - 一种航空遥感影像的最小二乘区域网匀色方法 - Google Patents

一种航空遥感影像的最小二乘区域网匀色方法 Download PDF

Info

Publication number
CN102129669A
CN102129669A CN 201110044025 CN201110044025A CN102129669A CN 102129669 A CN102129669 A CN 102129669A CN 201110044025 CN201110044025 CN 201110044025 CN 201110044025 A CN201110044025 A CN 201110044025A CN 102129669 A CN102129669 A CN 102129669A
Authority
CN
China
Prior art keywords
image
images
pixel
orthorectified
area
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.)
Granted
Application number
CN 201110044025
Other languages
English (en)
Other versions
CN102129669B (zh
Inventor
孙明伟
张祖勋
张剑清
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN 201110044025 priority Critical patent/CN102129669B/zh
Publication of CN102129669A publication Critical patent/CN102129669A/zh
Application granted granted Critical
Publication of CN102129669B publication Critical patent/CN102129669B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

本发明公开了一种航空遥感影像的最小二乘区域网匀色方法，该方法借鉴光束法区域网空中三角测量的思想，首先统计单像正射纠正影像的相邻影像叠加区域的像素色调信息，再构建每幅影像像素色调信息补偿参数的方程组，利用最小二乘平差方法，解算每幅影像的像素色调信息补偿参数，最后根据解算得到的色调补偿参数调整每幅影像的色调信息。本发明方法能够对大测区以及测区内影像内容差异较大的影像进行色调差异消除处理，在消除色差的过程中还能够方便地指定测区内部任意位置、任意数量的影像作为色调控制影像来控制整个测区的基本色调。

Description

一种航空遥感影像的最小二乘区域网匀色方法
技术领域
本发明属于遥感影像处理领域，尤其涉及一种航空遥感影像的匀色方法。
背景技术
在测绘生产中，数字正射影像（DOM）制作的一个主要步骤是影像镶嵌，影像镶嵌就是将完成单像正射纠正的相邻影像镶嵌成标准图幅的正射影像成果，或者将完成单像正射纠正的相邻影像镶嵌成整个测区的正射影像然后输出为标准图幅。在影像镶嵌过程中不仅要处理好相邻影像间几何接边问题，还要保证镶嵌后的影像整体色调一致、色彩过渡均匀。影像的颜色处理效果直接关系到最终正射影像成果的质量，研究并解决影像镶嵌过程中颜色问题的自动处理，对于提高测绘生产效率，保证正射影像产品质量具有重要意义。
一般来说，影响影像颜色的因素主要有两个：一是单幅影像内部的亮度分布不均匀；二是相邻影像之间存在颜色差异。已有大量学者对单幅影像内部亮度分布不均匀问题的处理（也称匀光处理）进行了研究。对相邻影像之间颜色差异问题的处理（也称匀色处理），较早提出的有基于Wallis滤波与色调基准影像相结合的方法[1]，该方法处理影像内容比较一致的测区时效果较为理想，但对较大测区、大比例尺的航空摄影影像或低空摄影影像而言，测区内影像之间的内容往往差异较大（如图1），此时，采用本方法进行匀色处理则无法取得理想效果。潘俊[2]等提出根据影像间的邻接关系、采用最短路径算法确定影像间的传递路径、并利用线性关系对相邻影像进行匀色处理，由于航空影像内容的复杂性，传递路径上的影像不一定严格满足线性关系，所以在测区较大的情况下容易出现传递路径上的误差累积，从而影响影像整体的匀色处理效果。
文中文中涉及的参考文献如下：
[1] 李治江. 彩色影像色调重建的理论与实践[D]. 武汉: 武汉大学:2005.4.
[2] 潘俊. 自动化的航空影像色彩一致性处理及接缝线网络生成方法研究[D]. 武汉：武汉大学，2008.11.
发明内容                     
针对现有技术存在的不足，本发明的目的是提供一种航空遥感影像的最小二乘区域网匀色方法，该方法能有效改善大测区、大比例尺的航空遥感影像的自动匀色效果。
为实现上述目的，本发明采用如下技术方案：
一种航空遥感影像的最小二乘区域网匀色方法，包括以下步骤：
步骤1，确认测区内所有待镶嵌的单像正射纠正影像在镶嵌过程中的接边关系，并获得相应的接边影像对；
步骤2，对步骤1获得的接边影像对，按照地理坐标对每对接边影像对中的两张单像正射纠正影像进行叠加，叠加后排除叠加区域内单像正射纠正影像的背景像素和纹理差异较大的像素，然后，分别统计接边影像对的叠加区域内未被排除的像素数和包括像素灰度均值和像素灰度方差的像素色调信息；
步骤3，计算测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数，该步骤进一步包括以下子步骤：
3-1、利用步骤2统计得到的像素色调信息构建矩阵方程，设测区内所有待镶嵌的单像正射纠正影像数为p，获得相应的接边影像对的数量为q，则：                                                
Figure 724775DEST_PATH_IMAGE001
，
其中，
，
Figure 276979DEST_PATH_IMAGE003
、
Figure 2011100440258100002DEST_PATH_IMAGE004
、
Figure 322165DEST_PATH_IMAGE005
、
Figure 2011100440258100002DEST_PATH_IMAGE006
、…
Figure 934412DEST_PATH_IMAGE007
分别为测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数；
为
Figure 143676DEST_PATH_IMAGE009
矩阵，且
Figure 2011100440258100002DEST_PATH_IMAGE010
；
Figure 109358DEST_PATH_IMAGE011
，为第i对接边影像对中的第一张影像x在叠加区域内的像素色调信息，
Figure 621767DEST_PATH_IMAGE013
为第i对接边影像对中的第二张影像y在叠加区域内的像素色调信息；
3-2、将p张单像正射纠正影像中色调较好的t张影像设置为色调控制影像，并构建矩阵方程：
Figure 2011100440258100002DEST_PATH_IMAGE014
，
其中，
Figure 936073DEST_PATH_IMAGE002
，、
Figure 588957DEST_PATH_IMAGE004
、
Figure 687363DEST_PATH_IMAGE005
、
Figure 47938DEST_PATH_IMAGE006
、…
Figure 762953DEST_PATH_IMAGE007
分别为
Figure 929492DEST_PATH_IMAGE015
张单像正射纠正影像的像素色调信息补偿参数；
[i,j]=
Figure 757770DEST_PATH_IMAGE017
，i
Figure 2011100440258100002DEST_PATH_IMAGE018
[0, t-1]，j [0, p-1]，i表示在t张色调控制影像中的第i张影像，当
Figure 285365DEST_PATH_IMAGE019
时，
Figure 458858DEST_PATH_IMAGE017
＝1，其它情况下，
Figure 797435DEST_PATH_IMAGE017
＝0；
Figure DEST_PATH_IMAGE020
为t行零常数的列向量，
Figure 952342DEST_PATH_IMAGE021
；
3-3、令，
Figure 173107DEST_PATH_IMAGE023
，构建矩阵方程
Figure DEST_PATH_IMAGE024
，则
Figure 681449DEST_PATH_IMAGE015
张单像正射纠正影像的像素色调信息补偿参数
Figure 546637DEST_PATH_IMAGE025
；
3-4、构建权矩阵
Figure DEST_PATH_IMAGE026
修正步骤3-3得到的像素色调信息补偿参数，则p张单像正射纠正影像的色调补偿参数
Figure 954003DEST_PATH_IMAGE027
，
其中，
Figure DEST_PATH_IMAGE028
，
Figure 334169DEST_PATH_IMAGE029
，
Figure DEST_PATH_IMAGE030
为叠加区域内的有效像素数，为阶单位矩阵；
步骤4，根据步骤3所得的像素色调补偿参数，调整每张单向正射纠正影像的像素色调信息。
上述步骤1具体如下：
基于测区内所有待镶嵌的单像正射纠正影像的摄站坐标构建Voronoi图，利用Voronoi图的邻接关系来确认镶嵌过程中所有待镶嵌的单像正射纠正影像之间的接边关系。
上述步骤2中排除叠加区域内纹理差异较大的像素采用如下的相关系数法：
利用预设尺寸的窗口对叠加区域进行逐个窗口内相关系数
Figure 84005DEST_PATH_IMAGE033
的计算，如果窗口位置的相关系数值小于预先设定的阈值，则排除掉两幅影像上该窗口位置内的像素，相关系数的计算公式如下：
Figure DEST_PATH_IMAGE034
其中，
Figure 589122DEST_PATH_IMAGE035
为窗口尺寸，
Figure DEST_PATH_IMAGE036
、
Figure 376949DEST_PATH_IMAGE037
分别为叠加区域内的两张影像在窗口内处的像素灰度值。
上述步骤4中是采用Wallis滤波算子调整镶嵌过程中所有待镶嵌的单像正射纠正影像的像素色调信息。
与现有技术相比，本发明具有以下优点和有益效果：
本发明方法能够对大测区以及测区内影像内容差异较大的影像进行色调差异消除处理，在消除色差的过程中还能够方便地指定测区内部任意位置、任意数量的影像作为色调控制影像来控制整个测区的基本色调。图7为利用本发明方法和常规方法对航空遥感影像进行匀色处理后的正射镶嵌结果对比，从图中可以看出，与常规匀色方法相比，采用本发明方法匀色处理后的镶嵌结果消除了影像间的色调差异。
附图说明     
图1为165张航空遥感影像的正射镶嵌图；
图2为二维平面上一系列输入点的Voronoi图；
图3为原始航空遥感影像与其对应的单像正射纠正影像， (a)为原始航空遥感影像，(b)为图(a)的单像正射纠正影像；
图4为利用相关系数法对单像正射纠正影像进行纹理差异较大的像素排除结果，(a)为左影像正射纠正后的局部区域，(b)为右影像正射纠正后的局部区域，(c)为左右正射影像叠加图，(d)为纹理差异较大的像素的排除结果；
图5为利用本发明方法选用不同的色调控制影像进行匀色处理的结果，(a)为不同位置的色调控制方案，(b)为选用黑框位置的色调控制影像的匀色处理结果，(c)为选用白框位置的色调控制影像的匀色处理结果；
图6为利用本发明方法对253张原始洞窟壁画影像匀色处理前后的拼接效果，(a)为253张原始洞窟壁画影像的正射镶嵌效果，(b)为利用本发明方法进行匀色处理后的正射镶嵌效果；
图7为利用本发明方法和常规方法处理后的正射镶嵌结果对比，(a)为原始影像的正射镶嵌效果，(b)为利用常规方法匀色处理后的正射镶嵌效果，(c)为采用本发明方法匀色处理后的正射镶嵌效果。
具体实施方式
本发明采用基于最小二乘的区域网匀色方法解决航空影像间的颜色差异问题，该方法借鉴光束法区域网空中三角测量的思想，首先统计单像正射纠正影像的相邻影像叠加区域的像素色调信息，再构建每幅影像像素色调信息补偿参数的方程组，利用最小二乘平差方法，解算每幅影像的像素色调信息补偿参数，最后根据解算得到的色调补偿参数调整每幅影像的色调信息。
为了更好地理解本发明的技术方案，下面将结合附图对本发明做进一步的详细说明，本发明的具体步骤如下：
步骤1，
对测区内所有待镶嵌的p张单像正射纠正影像，构建基于它们摄站坐标的Voronoi图，利用Voronoi图的邻接关系确定正射影像镶嵌过程中单像正射纠正影像两两之间是否存在接边，并获得q对接边影像对。
Voronoi图是计算几何学中的一个基础数据结构，它的输入是分布于二维平面上的一系列点，输出是对应于每个点的一系列多边形区域，输出的多边形区域应满足多边形区域内部的任意位置的点到与多边形区域所对应的输入点的几何距离最短，具体如图2所示。在正射影像镶嵌过程中，往往需要选取靠近摄影镜头中心位置的影像拼接到大图象中，所以可以利用Voronoi图的邻接关系确定哪些影像存在接边关系。
步骤2，
本步骤包括以下两个子步骤：
2-1、对步骤1获得的q对接边影像对，按照地理坐标对每对接边影像对中的两张单像正射纠正影像进行叠加，叠加后排除叠加区域内单像正射纠正影像的背景像素和纹理差异较大的像素，
叠加区域内单像正射纠正影像的背景像素，即原始航空影像在正射纠正过程中边缘不能覆盖的区域中的像素（如图3（b）所示的黑色区域）。叠加区域内单像正射纠正影像的背景像素的排除可采用本领域公知的技术手段来实现，本实施例中排除叠加区域内单像正射纠正影像的背景像素采用如下方法：对叠加区域内的两张影像，若其中一张影像的像素灰度值与背景像素灰度值相同，则排除该像素，同时排除另一张影像中的对应像素。
单像正射纠正影像上纹理差异较大的叠加区域，比如由于投影差造成的房屋屋顶位置的不同（如图4（a）、（b）所示）。对存在纹理差异的重叠区域，可利用相关系数法进行鉴别排除，本实施例中采用的具体方法如下：将重叠区域按照一定的窗口尺寸（比如5×5）进行逐个窗口内相关系数
Figure 541696DEST_PATH_IMAGE033
的计算，如果两幅影像上该窗口位置的相关系数
Figure 22355DEST_PATH_IMAGE033
小于一定的阈值（比如0.7），则将两幅影像上该窗口位置内的像素排除。相关系数
Figure 111534DEST_PATH_IMAGE033
的计算公式为：
   （1）
其中，
Figure 990814DEST_PATH_IMAGE035
为窗口尺寸，
Figure 829326DEST_PATH_IMAGE036
、分别为叠加区域内的两张影像在窗口内
Figure 66590DEST_PATH_IMAGE038
处的像素灰度值。
2-2、统计q对接边影像对的叠加区域内未被排除的像素数t r （即有效像素数）和像素色调信息，所述的像素色调信息包括像素灰度均值和像素灰度方差。
对第r对接边影像对中的第一张影像在叠加区域范围内的有效像素，利用下列公式分别计算其像素灰度均值
Figure 943279DEST_PATH_IMAGE039
和像素灰度方差
Figure DEST_PATH_IMAGE040
：
Figure 828058DEST_PATH_IMAGE041
              （2）
Figure DEST_PATH_IMAGE042
       （3）
其中，
t r 为第r对接边影像对在叠加区域内的有效像素数；
k为图像的色阶数，对每个通道为8Bit的图像，k=255；
h i 为叠加区域内第一张影像中像素灰度值为i且未被排除的像素个数，i
Figure 32775DEST_PATH_IMAGE018
[0, k]，h i 可由公式表示，w、h分别为第r对单像正射纠正影像对的叠加区域的像素宽和高；
Figure DEST_PATH_IMAGE044
，
Figure 762538DEST_PATH_IMAGE045
，
Figure DEST_PATH_IMAGE046
；g 0(x, y)表示第一张图像重叠区域内(x, y)坐标处的像素灰度值。
采用相同的方法计算得到第r对接边影像对的叠加区域内第二张影像的像素灰度均值和像素灰度方差
Figure DEST_PATH_IMAGE048
。重复本步骤，逐一统计所有接边影像对在叠加区域内像素灰度均值与像素灰度方差
Figure DEST_PATH_IMAGE050
，其中，
Figure 152303DEST_PATH_IMAGE029
，。
步骤3，计算p张单像正射纠正影像的像素色调信息补偿参数，影像的像素色调信息包括像素灰度均值和像素灰度方差，则像素色调信息补偿参数也包括像素灰度均值补偿参数和像素灰度方差补偿参数，下面将以像素灰度均值补偿参数的计算为例来进行具体说明：
3-1、利用影像叠加区域的像素灰度均值信息构建矩阵方程
Figure 902270DEST_PATH_IMAGE001
，矩阵方程
Figure 815999DEST_PATH_IMAGE001
的构建具体如下：
首先，假设第r对接边影像对是由第x、y两张影像构成，根据影像镶嵌拼接过程中叠加区域应该具有相同色调的基本要求，针对第x、y张影像的叠加区域的像素灰度均值，得到如下方程：
Figure DEST_PATH_IMAGE052
       （4）
其中，
Figure 681668DEST_PATH_IMAGE053
，
Figure DEST_PATH_IMAGE054
，且
Figure 449773DEST_PATH_IMAGE055
；
Figure 784939DEST_PATH_IMAGE039
、
Figure 700812DEST_PATH_IMAGE047
分别为第x张影像和第y张影像在叠加区域内的像素灰度均值；
Figure DEST_PATH_IMAGE056
、分别表示第x、y张影像的像素灰度均值的补偿参数。
将上述方程（4）写为：              
Figure DEST_PATH_IMAGE058
   （5）
方程（5）等价为：   （6）
其中，
Figure DEST_PATH_IMAGE060
是由
Figure 642726DEST_PATH_IMAGE061
组成的一维系数矩阵，
Figure DEST_PATH_IMAGE062
；当i＝x时，＝1；当i＝y时，＝－1；其它情况下，＝0；
然后，若按矩阵方程（6）写出q对接边影像对中所有单像正射纠正影像的矩阵方程，并将所得矩阵方程联立，即可构建矩阵方程：
          （7）
上述矩阵方程（7）可写为
Figure 116826DEST_PATH_IMAGE001
形式，，，
Figure DEST_PATH_IMAGE066
，
其中，
Figure 134646DEST_PATH_IMAGE065
，
Figure 779254DEST_PATH_IMAGE067
为第i张单像正射纠正影像的像素灰度均值的补偿参数，
Figure 704485DEST_PATH_IMAGE062
；
Figure 866476DEST_PATH_IMAGE064
，
Figure 541956DEST_PATH_IMAGE008
可看作由
Figure DEST_PATH_IMAGE068
构成的矩阵，i
Figure 623045DEST_PATH_IMAGE018
[0, q-1]，j
Figure 363468DEST_PATH_IMAGE018
[0, p-1]，
Figure 125887DEST_PATH_IMAGE068
中的i表示q对接边影像对中的第i对影像对，j表示p张单像正射纠正影像中的第j张影像。假设，第i对影像对由第x、y张影像构成，x
Figure 494420DEST_PATH_IMAGE018
[0, p-1]，y [0, p-1]，则当j＝x时，＝1；当j＝y时，
Figure 540240DEST_PATH_IMAGE068
＝－1，其它情况下，
Figure 373067DEST_PATH_IMAGE068
＝0；
Figure 999220DEST_PATH_IMAGE066
，为第i对接边影像对中的第一张影像在叠加区域内的像素灰度均值，为第i对接边影像对中的第二张影像在叠加区域内的像素灰度均值。
3-2、将p张单像正射纠正影像中色调较好的t张影像设为色调控制影像，并构建矩阵方程：
Figure 759070DEST_PATH_IMAGE014
；
将测区中一些色调较好的影像设为色调控制影像，使整个测区的最终色调与控制影像基本一致，色调控制影像的像素灰度均值补偿参数应该较小或者为0。假定第j张影像为色调控制影像（
Figure 711983DEST_PATH_IMAGE071
），设定其像素灰度均值补偿参数为0，可用如下方程表示：
Figure DEST_PATH_IMAGE072
     （8）
其中，
Figure 836933DEST_PATH_IMAGE073
为由组成的一维系数矩阵，
Figure 632720DEST_PATH_IMAGE062
，当i＝j时，
Figure 602950DEST_PATH_IMAGE074
＝1；其它情况下，
Figure 941528DEST_PATH_IMAGE074
＝0。
假设选定t张影像设为色调控制影像，则可列矩阵方程：
Figure 706221DEST_PATH_IMAGE075
      （9）
上述矩阵方程（9）可写为
Figure 740036DEST_PATH_IMAGE014
形式，
Figure 511028DEST_PATH_IMAGE065
，
Figure DEST_PATH_IMAGE076
、
Figure 438532DEST_PATH_IMAGE077
…
Figure DEST_PATH_IMAGE078
分别为
Figure 905286DEST_PATH_IMAGE015
张单像正射纠正影像的像素灰度均值补偿参数；
Figure 675664DEST_PATH_IMAGE079
，可看作由（i
Figure 832025DEST_PATH_IMAGE018
[0, t-1]，j
Figure 141784DEST_PATH_IMAGE018
[0, p-1]）构成的矩阵，
Figure 150191DEST_PATH_IMAGE017
中的i表示在t张色调控制影像中的第i张影像（其影像序号为z，有z
Figure 534424DEST_PATH_IMAGE018
[0, p-1]）。当j＝z时，
Figure 374204DEST_PATH_IMAGE017
＝1；其它情况下，
Figure 651601DEST_PATH_IMAGE017
＝0；
Figure 271938DEST_PATH_IMAGE021
。
3-3、将步骤3-1和步骤3-2得到的矩阵方程联立，得到
Figure 129036DEST_PATH_IMAGE081
（10），令，
Figure 458572DEST_PATH_IMAGE023
，则方程（10）可简化成
Figure 503888DEST_PATH_IMAGE024
，则
Figure 226994DEST_PATH_IMAGE015
张单像正射纠正影像的像素灰度均值补偿参数；
3-4、引入步骤2-2中统计的叠加区域的有效像素个数t r ，构建权矩阵P：
Figure DEST_PATH_IMAGE082
其中：
Figure 598249DEST_PATH_IMAGE028
，
Figure 865283DEST_PATH_IMAGE029
，
Figure 123570DEST_PATH_IMAGE030
为叠加区域内的有效像素数，
Figure 589186DEST_PATH_IMAGE031
为
Figure 848129DEST_PATH_IMAGE032
阶单位矩阵。
采用权矩阵P对步骤3-3得到像素灰度均值补偿参数
Figure 664776DEST_PATH_IMAGE025
进行修正，即影像的最终像素灰度均值补偿参数为
Figure 854317DEST_PATH_IMAGE027
，即可得到每张影像的像素灰度均值补偿参数m i 。
同理， 将矩阵
Figure 377703DEST_PATH_IMAGE065
中的
Figure 338705DEST_PATH_IMAGE067
改成影像的像素方差补偿参数v i 。
Figure 580331DEST_PATH_IMAGE062
；将矩阵
Figure 448930DEST_PATH_IMAGE066
中的
Figure 623559DEST_PATH_IMAGE083
改成
Figure DEST_PATH_IMAGE084
，
Figure 161988DEST_PATH_IMAGE085
，其中，
Figure DEST_PATH_IMAGE086
为第i对接边影像对中的第一张影像在叠加区域内的像素方差，
Figure 956156DEST_PATH_IMAGE087
为第i对接边影像对中的第二张影像在叠加区域内的像素方差。然后，采用步骤3-1～3-4解算出每张影像的像素灰度方差补偿参数v i 。
步骤4，
根据步骤3所得的影像像素灰度均值补偿参数和像素灰度方差补偿参数，并利用Wallis等滤波算子调整p张单向正射纠正影像的像素色调信息。
如果单像正射纠正影像为彩色图像，则首先将单像正射纠正影像按R、G、B三个波段分解为灰度图像，再采用上述步骤分别对相同波段的图像进行匀色处理，并将匀色处理完毕后的灰度图像按照R、G、B波段重新组合成彩色图像。
本发明的方法的匀色效果可见图5所示。图7为利用本发明方法和常规方法匀色处理后的正射镶嵌结果对比，从图中可以看出，跟常规匀色方法相比，采用本发明方法匀色处理后的镶嵌结果消除了影像间的色调差异。

Claims (4)

1.一种航空遥感影像的最小二乘区域网匀色方法，其特征在于，依次包括以下步骤：
步骤1，确认测区内所有待镶嵌的单像正射纠正影像在镶嵌过程中的接边关系，并获得相应的接边影像对；
步骤2，对步骤1获得的接边影像对，按照地理坐标对每对接边影像对中的两张单像正射纠正影像进行叠加，叠加后排除叠加区域内单像正射纠正影像的背景像素和纹理差异较大的像素，然后，分别统计接边影像对的叠加区域内未被排除的像素数和包括像素灰度均值和像素灰度方差的像素色调信息；
步骤3，计算测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数，该步骤进一步包括以下子步骤：
3-1、利用步骤2统计得到的像素色调信息构建矩阵方程，设测区内所有待镶嵌的单像正射纠正影像数为p，获得相应的接边影像对的数量为q，则：                                                
Figure 2011100440258100001DEST_PATH_IMAGE001
，
其中，
，
Figure 2011100440258100001DEST_PATH_IMAGE003
、
Figure 852852DEST_PATH_IMAGE004
、
Figure 2011100440258100001DEST_PATH_IMAGE005
、
Figure 429326DEST_PATH_IMAGE006
、…
Figure DEST_PATH_IMAGE007
分别为测区内所有待镶嵌的单像正射纠正影像的像素色调信息补偿参数；
Figure 762219DEST_PATH_IMAGE008
为
Figure DEST_PATH_IMAGE009
矩阵，且
Figure 698295DEST_PATH_IMAGE010
；
Figure DEST_PATH_IMAGE011
，
Figure 848654DEST_PATH_IMAGE012
为第i对接边影像对中的第一张影像x在叠加区域内的像素色调信息，
Figure DEST_PATH_IMAGE013
为第i对接边影像对中的第二张影像y在叠加区域内的像素色调信息；
3-2、将p张单像正射纠正影像中色调较好的t张影像设置为色调控制影像，并构建矩阵方程：，
其中，
Figure 298275DEST_PATH_IMAGE002
，
Figure 967154DEST_PATH_IMAGE003
、
Figure 288414DEST_PATH_IMAGE004
、
Figure 839481DEST_PATH_IMAGE005
、
Figure 576493DEST_PATH_IMAGE006
、…
Figure 37561DEST_PATH_IMAGE007
分别为
Figure DEST_PATH_IMAGE015
张单像正射纠正影像的像素色调信息补偿参数；
Figure 798231DEST_PATH_IMAGE016
[i,j]=
Figure DEST_PATH_IMAGE017
，i
Figure 164491DEST_PATH_IMAGE018
[0, t-1]，j
Figure 439614DEST_PATH_IMAGE018
[0, p-1]，i表示在t张色调控制影像中的第i张影像，当
Figure DEST_PATH_IMAGE019
时，＝1，其它情况下，
Figure 401940DEST_PATH_IMAGE017
＝0；
为t行零常数的列向量，；
3-3、令
Figure 334310DEST_PATH_IMAGE022
，
Figure DEST_PATH_IMAGE023
，构建矩阵方程
Figure 769970DEST_PATH_IMAGE024
，则
Figure 866584DEST_PATH_IMAGE015
张单像正射纠正影像的像素色调信息补偿参数；
3-4、构建权矩阵
Figure 613960DEST_PATH_IMAGE026
修正步骤3-3得到的像素色调信息补偿参数，则p张单像正射纠正影像的色调补偿参数
Figure DEST_PATH_IMAGE027
，
其中，
Figure 824361DEST_PATH_IMAGE028
，，
Figure 160534DEST_PATH_IMAGE030
为叠加区域内的有效像素数，为
Figure 430978DEST_PATH_IMAGE032
阶单位矩阵；
步骤4，根据步骤3所得的像素色调补偿参数，调整每张单向正射纠正影像的像素色调信息。
2.根据权利要求1所述的航空遥感影像的最小二乘区域网匀色方法，其特征在于：
所述步骤1具体如下：
基于测区内所有待镶嵌的单像正射纠正影像的摄站坐标构建Voronoi图，利用Voronoi图的邻接关系来确认镶嵌过程中所有待镶嵌的单像正射纠正影像之间的接边关系。
3.根据权利要求1或2所述的航空遥感影像的最小二乘区域网匀色方法，其特征在于：
所述步骤2中排除叠加区域内纹理差异较大的像素采用如下的相关系数法：
利用预设尺寸的窗口对叠加区域进行逐个窗口内相关系数
Figure DEST_PATH_IMAGE033
的计算，如果窗口位置的相关系数值小于预先设定的阈值，则排除掉两幅影像上该窗口位置内的像素，相关系数的计算公式如下：
Figure 23950DEST_PATH_IMAGE034
    
其中，
Figure DEST_PATH_IMAGE035
为窗口尺寸，
Figure 827346DEST_PATH_IMAGE036
、
Figure DEST_PATH_IMAGE037
分别为叠加区域内的两张影像在窗口内
Figure 737533DEST_PATH_IMAGE038
处的像素灰度值。
4.根据权利要求1或2所述的航空遥感影像的最小二乘区域网匀色方法，其特征在于：
所述步骤4中是采用Wallis滤波算子调整镶嵌过程中所有待镶嵌的单像正射纠正影像的像素色调信息。
CN 201110044025 2011-02-24 2011-02-24 一种航空遥感影像的最小二乘区域网匀色方法 Expired - Fee Related CN102129669B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110044025 CN102129669B (zh) 2011-02-24 2011-02-24 一种航空遥感影像的最小二乘区域网匀色方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110044025 CN102129669B (zh) 2011-02-24 2011-02-24 一种航空遥感影像的最小二乘区域网匀色方法

Publications (2)

Publication Number Publication Date
CN102129669A true CN102129669A (zh) 2011-07-20
CN102129669B CN102129669B (zh) 2012-07-11

Family

ID=44267745

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110044025 Expired - Fee Related CN102129669B (zh) 2011-02-24 2011-02-24 一种航空遥感影像的最小二乘区域网匀色方法

Country Status (1)

Country Link
CN (1) CN102129669B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103955897A (zh) * 2014-04-22 2014-07-30 中国资源卫星应用中心 一种基于最邻近统计的ccd图像色差消除方法
CN109272465A (zh) * 2018-09-15 2019-01-25 武汉智觉空间信息技术有限公司 一种航空影像色彩一致性处理算法
CN110176053A (zh) * 2019-05-24 2019-08-27 武汉大势智慧科技有限公司 一种大规模实景三维整体匀色方法
CN111402167A (zh) * 2020-03-19 2020-07-10 重庆市勘测院 一种影像色彩归一方法
CN111447426A (zh) * 2020-05-13 2020-07-24 中测新图(北京)遥感技术有限责任公司 一种影像色彩校正方法以及装置
CN111754590A (zh) * 2020-05-14 2020-10-09 北京吉威空间信息股份有限公司 基于全球色彩特征库遥感影像自动匀色的方法
CN112884676A (zh) * 2021-03-18 2021-06-01 国家海洋信息中心 一种基于分空间渐进控制的大范围航空遥感影像调色方法
CN113012276A (zh) * 2021-01-27 2021-06-22 中国科学院空天信息创新研究院 基于辐射度的地表高分辨率光谱信息遥感反演方法
CN117593388A (zh) * 2023-12-04 2024-02-23 北京数慧时空信息技术有限公司 遥感影像匀色方法、装置、电子设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1496705A2 (fr) * 2003-07-11 2005-01-12 Alcatel Compression contextuelle d'images numérique
CN101216555A (zh) * 2007-12-27 2008-07-09 武汉大学 Rpc模型参数提取方法和几何纠正方法
CN101251926A (zh) * 2008-03-20 2008-08-27 北京航空航天大学 一种基于局部轮廓协方差矩阵的遥感图像配准方法
CN101876701A (zh) * 2010-07-02 2010-11-03 中国测绘科学研究院 一种侧视雷达遥感影像定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1496705A2 (fr) * 2003-07-11 2005-01-12 Alcatel Compression contextuelle d'images numérique
CN101216555A (zh) * 2007-12-27 2008-07-09 武汉大学 Rpc模型参数提取方法和几何纠正方法
CN101251926A (zh) * 2008-03-20 2008-08-27 北京航空航天大学 一种基于局部轮廓协方差矩阵的遥感图像配准方法
CN101876701A (zh) * 2010-07-02 2010-11-03 中国测绘科学研究院 一种侧视雷达遥感影像定位方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103955897A (zh) * 2014-04-22 2014-07-30 中国资源卫星应用中心 一种基于最邻近统计的ccd图像色差消除方法
CN109272465A (zh) * 2018-09-15 2019-01-25 武汉智觉空间信息技术有限公司 一种航空影像色彩一致性处理算法
CN110176053A (zh) * 2019-05-24 2019-08-27 武汉大势智慧科技有限公司 一种大规模实景三维整体匀色方法
CN111402167B (zh) * 2020-03-19 2023-04-18 重庆市勘测院 一种影像色彩归一方法
CN111402167A (zh) * 2020-03-19 2020-07-10 重庆市勘测院 一种影像色彩归一方法
CN111447426A (zh) * 2020-05-13 2020-07-24 中测新图(北京)遥感技术有限责任公司 一种影像色彩校正方法以及装置
CN111754590A (zh) * 2020-05-14 2020-10-09 北京吉威空间信息股份有限公司 基于全球色彩特征库遥感影像自动匀色的方法
CN111754590B (zh) * 2020-05-14 2024-04-02 北京吉威空间信息股份有限公司 基于全球色彩特征库遥感影像自动匀色的方法
CN113012276A (zh) * 2021-01-27 2021-06-22 中国科学院空天信息创新研究院 基于辐射度的地表高分辨率光谱信息遥感反演方法
CN113012276B (zh) * 2021-01-27 2021-09-24 中国科学院空天信息创新研究院 基于辐射度的地表高分辨率光谱信息遥感反演方法
CN112884676A (zh) * 2021-03-18 2021-06-01 国家海洋信息中心 一种基于分空间渐进控制的大范围航空遥感影像调色方法
CN117593388A (zh) * 2023-12-04 2024-02-23 北京数慧时空信息技术有限公司 遥感影像匀色方法、装置、电子设备及介质
CN117593388B (zh) * 2023-12-04 2024-08-20 北京数慧时空信息技术有限公司 遥感影像匀色方法、装置、电子设备及介质

Also Published As

Publication number Publication date
CN102129669B (zh) 2012-07-11

Similar Documents

Publication Publication Date Title
CN102129669B (zh) 一种航空遥感影像的最小二乘区域网匀色方法
US11783446B2 (en) Large-field-angle image real-time stitching method based on calibration
KR101742120B1 (ko) 영상 처리 장치 및 방법
CN101276465B (zh) 广角图像自动拼接方法
EP3111647B1 (en) Image stitching and automatic-color correction
CN109859105B (zh) 一种无参数影像自然拼接方法
US8374428B2 (en) Color balancing for partially overlapping images
US20080143744A1 (en) Gradient-domain compositing
CN106060493A (zh) 多源投影无缝边缘拼接方法及系统
CN104794683B (zh) 基于围绕渐变拼缝区域平面扫描的视频拼接方法
US8577139B2 (en) Method of orthoimage color correction using multiple aerial images
CN104574501A (zh) 一种针对复杂三维场景的高质量纹理映射方法
CN108550113A (zh) 图像扫描输出方法、装置、计算机设备和存储介质
CN102509294B (zh) 一种基于单幅图像的全局深度估计方法
US20230125649A1 (en) Image inpainting method and electronic device
CN115239820A (zh) 一种分体式飞行车辆鸟瞰图实时拼接及车位检测方法
CN105005964A (zh) 基于视频序列影像的地理场景全景图快速生成方法
CN107909544A (zh) 一种图像校正方法及系统
CN108629742B (zh) 真正射影像阴影检测与补偿方法、装置及存储介质
CN106780326A (zh) 一种提高全景图像清晰度的融合方法
CN107403448B (zh) 代价函数生成方法和代价函数生成装置
CN102905077B (zh) 一种图像暗角亮度调节方法及装置
CN102802020B (zh) 监测双目立体视频视差信息的方法和装置
CN106412461A (zh) 视频拼接方法
CN116109681A (zh) 图像融合方法、装置、电子设备及可读存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120711