CN101281643B - 一种退化函数随空间变化图像的分块复原和拼接方法 - Google Patents
一种退化函数随空间变化图像的分块复原和拼接方法 Download PDFInfo
- Publication number
- CN101281643B CN101281643B CN2008100604932A CN200810060493A CN101281643B CN 101281643 B CN101281643 B CN 101281643B CN 2008100604932 A CN2008100604932 A CN 2008100604932A CN 200810060493 A CN200810060493 A CN 200810060493A CN 101281643 B CN101281643 B CN 101281643B
- Authority
- CN
- China
- Prior art keywords
- image
- piecemeal
- region
- restoration
- partitioning
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000011084 recovery Methods 0.000 claims description 10
- 230000007704 transition Effects 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 abstract 7
- 230000000694 effects Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000000007 visual effect Effects 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种退化函数随空间变化图像的分块复原和拼接方法,首先以等晕区为依据对退化函数随空间变化的图像进行分块并向外延伸,并对每一个图像块进行图像复原;然后对各复原的图像块去除其包含明显振铃波纹的小部分边缘区域,保留仍大于等晕区的图像块;最后将各图像块拼接成完整的图像。该方法能够消除退化函数随空间变化的图像分块复原处理时产生的振铃波纹。本发明在诸如卫星遥感、飞机航拍图像复原,医学图像处理等众多领域可广泛应用。
Description
技术领域
本发明涉及计算机图像处理技术,尤其是涉及计算机图像处理中退化函数随空间变化图像的分块复原和拼接方法,用以处理分块复原图像边缘产生的振铃效应及灰度突变。
背景技术
光学图像的获取过程中,由于受到光学系统像差,加工装校误差,CCD采样量化,拍摄平台振动,大气扰动,噪声等因素影响,拍摄图像会有一定退化,为了尽可能得到清晰准确的原物体信息,人们采用各种图像复原算法对退化的图像进行处理。如果整幅图像的退化情况类似,能近似认为是在一个等晕区内,那么可以用一个退化函数来模拟退化情况,然后依据这个退化函数恢复图像。但如果图像各视场退化情况不同,不能用近似的等晕区来表达整幅图像,那么就无法用统一的退化函数来描述退化情况,使用任何一个单一视场的退化函数处理图像无法得到理想的效果。于是图像分块处理的方法被提出来解决这一问题。整幅图像依据不同的退化情况划分成不同的块,每个图像块内可认为退化情况一致,即为等晕区,在等晕区内可使用传统的图像处理方法恢复图像,然后将各处理好的图像按原来的排列拼接起来,得到整幅图像。这一方法很好的解决了退化函数随空间变化的图像复原问题,然而图像分块造成了边界截断,复原图像块的边界都会出现振铃波纹,拼接后的图像在各等晕区边界处都有生硬的突变。
发明内容
本发明提供了一种退化函数随空间变化图像的分块复原和拼接方法,该方法为一种能够消除边界振铃效应的算法,对于处理退化函数随空间变化图像的复原有着重大的意义。
一种退化函数随空间变化图像的分块复原和拼接方法,包括:
(1)以等晕区为依据对退化函数随空间变化的图像进行分块,为了消除等晕区边界附近的振铃波纹,把边界外移,即图像的分块并不是在等晕区的边界处即截止,而是往外延伸一段距离,这一距离可以根据图像复原的实际情况进行设定,应大于振铃波纹所分布的区域。
(2)对每一个图像块进行图像复原处理,可根据图像的实际情况选用各种传统图像复原方法,例如维纳滤波,最小二乘复原算法,最大熵复原算法,受限制空间自适应算法,共轭梯度法等等。
(3)对各复原的图像块,首先去除其包含明显振铃波纹的小部分边缘区域,此时剩下的图像块仍大于等晕区。
(4)将各图像块拼接成完整的图像,方法为:各图像块彼此重叠的部分,以其到等晕区边界的距离为基础构造加权系数,依此加权系数计算重叠区域每一像素的光强值。
图像拼接的具体实施可用下面的公式表达。
设相邻两图像块分别以X、Y表示,两者拼接后图像以Z表示。以一维边界为例,X、Y的大小分别是M,N,图像边界自等晕区边界往外延伸L,即等晕区边界在X的第M-L行,在Y的第L行。X、Y在经过图像复原处理之后,去除了包含振铃波纹的l行,那么剩下的两图像块的重叠区宽度d=L-l,即此时的等晕区边界在X的第M-d行,在Y的第d行。两侧图像重叠区各过渡元素的权重系数按递推:
这样得到的图像拼接的边界振铃波纹幅度明显减小,而且等晕区过渡变得自然平滑。
需要说明的是,分块的图像可以是矩形,圆形或者其它任何形状,为了简化描述,本发明仅以一维矩形分块为例,该方法可以很容易的推知二维矩形或者其它形状。
本发明适合于退化函数随空间发生变化的图像复原,比如大视场引起的各点像差分布不同,散焦图像中物距不同的各点的散焦程度不同,非线性匀速运动物体各点运动情况不同等等,在诸如卫星遥感、飞机航拍图像复原,医学图像处理等众多领域可广泛应用。
附图说明
图1是本发明与传统分块图像处理拼接情况的一维简化对比示意图。
图2是本发明与传统分块图像处理的边缘灰度对比示意图。
图3(a)是处理前的原图。
图3(b)是上下两个区域经过不同卷积核模糊的图像(上半部分是衍射卷积核,下半部分是高斯卷积核)。
图4(a)是将图3(b)上半部分取出来用其对应(衍射)卷积核进行复原的结果。
图4(b)是将图3(b)下半部分取出来用其对应(高斯)卷积核进行复原的结果;图4(c)是将图4(a)和图4(b)直接拼接起来的结果。
图5(a)是将图3(b)上半部分以及相邻10个像素取出来用其对应(衍射)卷积核进行复原的结果。
图5(b)是将图3(b)下半部分以及相邻10个像素取出来用其对应(衍射)卷积核进行复原的结果。
图5(c)是将图5(a)和图5(b)以本发明所说明的方法拼接起来的结果。
图6是上下两个区域经过不同卷积核模糊的图像(上半部分是衍射卷积核,下半部分是高斯卷积核),模糊程度比图3(b)更大,图像复原需要的迭代次数更多,复原后的图像边缘振铃波纹更明显。
图7是采用单一模型对图6图像进行恢复的效果图。图7(a)是采用衍射卷积核恢复的结果,图7(b)是采用高斯卷积核恢复的结果。
图8(a)-(c)是对图6的传统分块复原与拼接的结果,含义同图4(a)-(c)。
图9(a)-(c)是用本发明所说明的方法对图6的分块复原并拼接的结果,含义同图5(a)-(c)。
具体实施方式
附图1所示为本发明与传统分块图像处理拼接情况的一维简化对比示意图。图中,(a)是传统分块图像处理一维简化示意图,AB表示整幅图像,C将图像分为AC和BC两个等晕区,分别对AC和CB进行图像复原处理之后,在其边缘,即阴影所示部分会出现振铃波纹,拼接复原之后的图像为AC、CB,得到的最后图像AB等晕区边界部分留有振铃波纹,降低了图像质量。(b)是本发明所描述的图像块拼接的示意图,AB、C的定义同前,D、E表示从边界C处往外延伸,使得图像AD、EB的振铃波纹包含在重叠区域内,被振铃破坏的信息可由另一幅图像的重叠区来补充。但在实际情况中振铃存在的区域并不易判断,不论是经过机器还是人工判读,只能确定其分布的大概区域,不过这对我们后面的复原图像块拼接并不会产生太大影响,因为图像块的振铃波纹即使没有全部去除,其重叠区域还要进过渐变加权系数的处理,会进一步削弱其影响。
附图2为本发明与传统分块图像处理的边缘灰度对比示意图。(a)是传统分块图像处理等晕区边界两边的灰度示意图,在经过图像复原等运算之后,图像块边界附近的灰度均值可能不一致,如图所示将在边界处形成阶梯,这不符合图像数据平滑的准则。(b)是本发明所描述的等晕区边界外延拼接的方法,利用两个图像块的重叠区数据,以其到边界的距离构成加权系数,得到一个灰度渐变过渡的边界。
实施例1:
附图3(a)为实验用的原图像,像素数为256*256。
附图3(b)是将原图像分为上下两个部分,皆为128*256像素,上半部分与15*15的衍射退化核进行卷积,下半部分与15*15的高斯退化核卷积形成的模糊图像,像素数仍为256*256。
其中,衍射退化核的解析表达式为:
式中x,y表示像素所在的位置,J1是一阶贝塞尔函数。
高斯退化核的解析表达式为:
同样x,y表示像素所在的位置,σ是表示高斯函数模糊程度的参数。
模糊的过程即为原图像与这两个退化核卷积的过程。
附图4(a)为取出图3(b)模糊图像的上半部分,使用其对应卷积核,即衍射核进行复原的图像结果,像素数为128*256。
附图4(b)为取出图3(b)模糊图像的下半部分,使用其对应卷积核,即高斯核进行复原的图像结果,像素数为128*256。
附图4(c)是将图4(a)和图4(b)根据图1(a)中所示方法直接拼接起来得到的图像,像素数为256*256。
附图5(a)为取出图3(b)模糊图像上半部分并向下延伸10个像素(10个像素是通过观察复原图像块边缘振铃波纹分布区域而确定的),使用其对应卷积核,即衍射核进行复原的图像结果,像素数为138*256。
附图5(b)为取出图3(b)模糊图像下半部分并向上延伸10个像素,使用其对应卷积核,即高斯核进行复原的图像结果,像素数为138*256。
附图5(c)是将图5(a)和图5(b)去除边界处包含振铃波纹的5个像素,剩下的两图像块的大小皆为133*256像素,然后以本发明所描述的方法在重叠区利用加权系数进行拼接的图像结果,像素数为256*256。
加权系数的计算方法如下:
重叠区各过渡元素的权重系数按50%/5=10%递推,那么:
Z124=X124*90%+Y1*10%
Z125=X125*80%+Y2*20%
Z126=X126*70%+Y3*30%
Z128=X128*50%+Y10*50%
Z129=X129*50%+Y11*50%
Z132=X132*20%+Y9*80%
Z133=X133*10%+Y10*90%
可以看到图5(c)拼接后的图像的振铃波纹得到很好的抑制,边界的过渡也十分平滑,图像复原质量很好。
实施例2:
附图6是模糊程度更大的上下两个区域经过不同卷积核模糊的图像(上半部分是衍射卷积核,下半部分是高斯卷积核),由于模糊程度比图3(b)更大,图像复原所需要的迭代次数更多,复原后的图像边缘振铃效应也更明显。
图7是采用单一退化模型恢复的效果,图7(a)是采用衍射卷积核恢复的结果,因为采用单一适合上半部分的卷积核模型,上半部分图像恢复效果较好,而下半部分仍然模糊。图7(b)是采用高斯卷积核恢复的结果。因为采用了下半部分的卷积核模型,下半部分图像恢复效果较好,而上半部分明显失真。
图8(a)-(c)是对图6图像的传统分块拼接结果,含义同图4(a)-(c),可以看到边界的波纹更多。
图9(a)-(c)是用本发明所说明的方法对图6的分块拼接结果,含义同图5(a)-(c),但图像块等晕区往外延伸了20个像素,去除了包含振铃波纹的10个像素后,重叠区各过渡元素的权重系数按50%/10=5%递推:
Z119=X119*95%+Y1*5%
Z120=X120*90%+Y2*10%
Z121=X121*85%+Y3*15%
Z128=X128*50%+Y10*50%
Z129=X129*50%+Y11*50%
Z137=X137*10%+Y19*90%
Z138=X138*5%+Y20*95%
可以看到图9(c)拼接后的图像的振铃波纹得到有效抑制,边界的过渡也十分平滑,图像得到有效恢复。
Claims (2)
1.一种退化函数随空间变化图像的分块复原和拼接方法,首先以等晕区为依据对退化函数随空间变化的图像进行分块,并将分块向外延伸,再对每一个图像块进行图像复原,然后对各复原的图像块去除其包含明显振铃波纹的小部分边缘区域,保留仍大于等晕区的图像块,最后将各图像块拼接成完整的图像,其特征在于:图像块拼接时利用相邻等晕区图像重叠部分边界附近一段区域按该像素到边界的距离为权重系数对两边图像灰度进行加权运算,得到一个渐变的平缓过渡边界。
2.按权利要求1所述的分块复原和拼接方法,其特征在于:所述的将分块向外延伸是指图像的分块在等晕区的边界处往外延伸一段距离,这一距离应大于振铃波纹所分布的区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100604932A CN101281643B (zh) | 2008-04-23 | 2008-04-23 | 一种退化函数随空间变化图像的分块复原和拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100604932A CN101281643B (zh) | 2008-04-23 | 2008-04-23 | 一种退化函数随空间变化图像的分块复原和拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101281643A CN101281643A (zh) | 2008-10-08 |
CN101281643B true CN101281643B (zh) | 2010-06-02 |
Family
ID=40014090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100604932A Expired - Fee Related CN101281643B (zh) | 2008-04-23 | 2008-04-23 | 一种退化函数随空间变化图像的分块复原和拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101281643B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719251A (zh) * | 2016-01-19 | 2016-06-29 | 浙江大学 | 一种用于大像移线性模糊的压缩降质图像复原方法 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101937561A (zh) * | 2010-09-19 | 2011-01-05 | 华中科技大学 | 一种天基大视场成像点源/斑状目标的图像恢复增强方法 |
CN102867318B (zh) * | 2012-09-25 | 2015-11-11 | 哈尔滨工业大学 | 子图像加权叠加的整幅折返运动离散退化图像构造方法 |
CN103747187A (zh) * | 2013-12-31 | 2014-04-23 | 浙江大学 | 光纤陀螺振动探测和编码曝光的遥感图像快速恢复方法 |
CN104217411B (zh) * | 2014-09-02 | 2017-10-31 | 济南大学 | 一种不规则破碎单面图像快速拼接方法 |
CN104897688B (zh) * | 2015-06-11 | 2017-08-25 | 龚厚志 | 一种化纤纸管外观检测过程中图像精确处理方法及系统 |
CN110047031A (zh) * | 2019-03-26 | 2019-07-23 | 深兰科技(上海)有限公司 | 一种像素段拼接的方法和装置 |
CN110399703A (zh) * | 2019-04-02 | 2019-11-01 | 泰州市海陵区一马商务信息咨询有限公司 | 电子设备数据管理系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6108443A (en) * | 1997-08-25 | 2000-08-22 | Fuji Photo Film Co., Ltd. | Method and apparatus for generating a converted image with reduced quality degradation resulting from the conversion |
CN1264482A (zh) * | 1997-07-16 | 2000-08-23 | 三星电子株式会社 | 信号自适应滤波方法,信号自适应滤波器及其存储程序的计算机可读介质 |
CN1447233A (zh) * | 2003-01-23 | 2003-10-08 | 同济大学 | 多媒体数据保护方法 |
CN1945623A (zh) * | 2006-10-26 | 2007-04-11 | 丹东东方测控技术有限公司 | 数字图像分块局部加密的可去除可见水印方法 |
-
2008
- 2008-04-23 CN CN2008100604932A patent/CN101281643B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1264482A (zh) * | 1997-07-16 | 2000-08-23 | 三星电子株式会社 | 信号自适应滤波方法,信号自适应滤波器及其存储程序的计算机可读介质 |
US6108443A (en) * | 1997-08-25 | 2000-08-22 | Fuji Photo Film Co., Ltd. | Method and apparatus for generating a converted image with reduced quality degradation resulting from the conversion |
CN1447233A (zh) * | 2003-01-23 | 2003-10-08 | 同济大学 | 多媒体数据保护方法 |
CN1945623A (zh) * | 2006-10-26 | 2007-04-11 | 丹东东方测控技术有限公司 | 数字图像分块局部加密的可去除可见水印方法 |
Non-Patent Citations (6)
Title |
---|
傅南翔,郭子祺,袁泉."北京一号"多光谱影像复原方法研究及实现.遥感应用.2007,59-62,及插页. |
傅南翔,郭子祺,袁泉."北京一号"多光谱影像复原方法研究及实现.遥感应用.2007,59-62,及插页. * |
葛西旺.一种精确的图像拼接方法.安徽建筑工业学院学报( 自然科学版)15 1.2007,15(1),69-72. |
葛西旺.一种精确的图像拼接方法.安徽建筑工业学院学报( 自然科学版)15 1.2007,15(1),69-72. * |
陈小蔷等.基于 Neumann 边界条件的图像局部复原.信号处理20 4.2004,20(4),399-402. |
陈小蔷等.基于 Neumann 边界条件的图像局部复原.信号处理20 4.2004,20(4),399-402. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719251A (zh) * | 2016-01-19 | 2016-06-29 | 浙江大学 | 一种用于大像移线性模糊的压缩降质图像复原方法 |
CN105719251B (zh) * | 2016-01-19 | 2018-06-19 | 浙江大学 | 一种用于大像移线性模糊的压缩降质图像复原方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101281643A (zh) | 2008-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101281643B (zh) | 一种退化函数随空间变化图像的分块复原和拼接方法 | |
US10339643B2 (en) | Algorithm and device for image processing | |
Cheng et al. | Inpainting for remotely sensed images with a multichannel nonlocal total variation model | |
CN110675336A (zh) | 一种低照度图像增强方法及装置 | |
US6959117B2 (en) | Method and apparatus for deblurring and re-blurring image segments | |
CN106056629A (zh) | 通过运动物体检测和扩展去除鬼影的高动态范围成像方法 | |
CN104252700A (zh) | 一种红外图像的直方图均衡化方法 | |
CN101853490A (zh) | 一种基于人类视觉特性的仿生图像复原方法 | |
CN102142132A (zh) | 基于模块的图像修复方法 | |
CN101937561A (zh) | 一种天基大视场成像点源/斑状目标的图像恢复增强方法 | |
Shacham et al. | Blind restoration of atmospherically degraded images by automatic best step-edge detection | |
FR2978855A1 (fr) | Procede et dispositif de calcul d'une carte de profondeur a partir d'une unique image | |
US20050073702A1 (en) | Robust recursive envelope operators for fast retinex-type processing | |
Teng et al. | Adaptive morphological filtering method for structural fusion restoration of hyperspectral images | |
Tarbell et al. | On the relation between magnetic field structures and granulation | |
CN115880177A (zh) | 聚合上下文和增强细节的全分辨率低照度图像增强方法 | |
Secker | A ring median filter for digital images | |
CN110738624B (zh) | 一种区域自适应的图像去雾系统及方法 | |
CN102183244A (zh) | 一种航空遥感影像的比值匀光方法 | |
CN117058019A (zh) | 一种基于金字塔增强网络的低光照下目标检测方法 | |
Banerjee et al. | Bacterial foraging-fuzzy synergism based image Dehazing | |
Esmaeilzehi et al. | DSegAN: A deep light-weight segmentation-based attention network for image restoration | |
CN116934583A (zh) | 基于深浅特征融合网络的遥感图像超分辨率算法 | |
Tang et al. | De-ghosting method for image stitching | |
Babic et al. | An efficient noise removal and edge preserving convolution filter |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100602 Termination date: 20120423 |