CN105931210A - 一种高分辨率图像重建方法 - Google Patents

一种高分辨率图像重建方法 Download PDF

Info

Publication number
CN105931210A
CN105931210A CN201610234555.1A CN201610234555A CN105931210A CN 105931210 A CN105931210 A CN 105931210A CN 201610234555 A CN201610234555 A CN 201610234555A CN 105931210 A CN105931210 A CN 105931210A
Authority
CN
China
Prior art keywords
pixel
resolution
image
low
resolution image
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
CN201610234555.1A
Other languages
English (en)
Other versions
CN105931210B (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.)
Luoyang Institute of Electro Optical Equipment AVIC
Original Assignee
Luoyang Institute of Electro Optical Equipment AVIC
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 Luoyang Institute of Electro Optical Equipment AVIC filed Critical Luoyang Institute of Electro Optical Equipment AVIC
Priority to CN201610234555.1A priority Critical patent/CN105931210B/zh
Publication of CN105931210A publication Critical patent/CN105931210A/zh
Application granted granted Critical
Publication of CN105931210B publication Critical patent/CN105931210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20028Bilateral filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Abstract

本发明涉及一种高分辨率图像重建方法,在同一场景采集两帧以上的低分辨率图像,从采集的图像中选取一帧图像作为基准帧图像;计算其余帧的低分辨率图像的各像素点相对于基准帧图像的偏移量,根据该偏移量将其余帧的低分辨率图像以基准帧图像为基准进行配准;设定分辨率放大倍数,将低分辨率图像映射到待复原的高分辨率图像栅格上;对映射到各个栅格上的像素点的像素值通过融合计算得到各个栅格顶点的像素值,从而得到该场景的高分辨率图像。本发明在保证高分辨率图像细节丰富的同时,算法简单、速度快,容易实现。

Description

一种高分辨率图像重建方法
技术领域
本发明涉及一种高分辨率图像重建方法,属于数字图像处理的技术领域。
背景技术
高分辨率图像在图像处理的技术领域中具有十分重要的作用,广泛应用于在卫星影像,交通监控或者目标识中。但是,由于受到各种拍摄条件的限制,如:拍摄设备工艺、周边环境气候等,往往获取的图像分辨率并不理想,这就需要对其进行后期处理,使其转换成满足要求的高分辨率图像。对于图像超分辨率重建领域中,通常要将多帧同一场景的低分辨率图像通过配准方法映射到高分辨率图像栅格中,实现将低分辨率图像转换成高分辨率图像。但是现有技术方法虽然可以获取高分辨率图像,但是,普遍方法都比较复杂,计算效率低。
发明内容
本发明的目的在于克服现有技术的不足,提出了一种高分辨率图像重建方法,该方法可以快速获取细节丰富的高分辨率图像。
本发明是通过如下方案予以实现的:
一种高分辨率图像重建方法,步骤如下:
步骤1,在同一场景采集两帧以上的低分辨率图像,从采集的图像中选取一帧图像作为基准帧图像;计算其余帧的低分辨率图像的各像素点相对于基准帧图像的偏移量,根据该偏移量将余帧的低分辨率图像以基准帧图像为基准进行配准;
步骤2,设定分辨率放大倍数,以基准帧图像各像素点乘以所述分辨率放大倍数为目标位置,将基准帧图像的像素点映射到待复原的高分辨率图像栅格上;将根据偏移量配准后的其余帧的低分辨率图像的像素点乘以所述分辨率放大倍数为目标位置,将配准后的其余帧的低分辨率图像的像素点映射到待复原的高分辨率图像栅格上;
步骤3,对步骤2中映射到各个栅格上的像素点的像素值通过融合计算得到各个栅格顶点的像素值,从而得到该场景的高分辨率图像。
进一步的,步骤1中,其余帧的低分辨率图像的各像素点相对于基准帧图像的偏移量包括:位置偏移量和角度偏移量;位置偏移量又包括:水平偏移值和垂直偏移值。
进一步的,采集的低分辨率图像数量的为分辨率放大倍数的2倍到4倍。
进一步的,步骤2中用(xb,yb)表示基准帧图像的像素点;f表示设定的分辨率放大倍数;则基准帧图像的像素点映射到待复原的高分辨率图像栅格上后得到的像素点为(f·xb,f·yb);
用(xa,ya)表示其余的低分辨率图像的像素点;用xd和yd分别表示该其余低分辨率图像相对于基准帧图像在X轴和Y轴的位置偏移量;其余帧的低分辨率图像的像素点映射到待复原的高分辨率图像栅格上后得到的像素点为(f·xa+f·xd,f·ya+f·yd)。
进一步的,步骤3中利用双边滤波器对映射到各个栅格上的像素点的像素值进行融合,表达式如下:
g ( i , j ) = Σ k , l f ( k , l ) w ( i , j , k , l ) Σ k , l w ( i , j , k , l )
其中,g(i,j)表示通过融合后得到的高分辨率图像的像素值;f(k,l)表示映射到各个栅格上的像素点的像素值;w(i,j,k,l)表示对应f(k,l)的权重系数。
本发明和现有技术相比的有益效果是:
本发明提出了一种高分辨率图像重建方法,对同一场景采集多帧低分辨率图像,选取其中一帧图像作为基准帧图像,其余图像以基准帧图像为基准进行配准。然后,对所有低分辨率图像进行放大并映射到高分辨率图像栅格上,再通过数据融合的方法得到高分辨率图像。通过采集多幅图像,获取丰富的图像信息,利用本发明方法实现对高分辨率图像的重建;本发明在保证高分辨率图像细节丰富的同时,算法简单、速度快,容易实现。
附图说明
图1是本发明实施例中基准图像映射到高分辨率图像栅格示意图;
图2是本发明实施例中所有低分辨率图像映射到高分辨率图像栅格示意图;
图3是本发明实施例中一种高分辨率图像重建方法的流程示意图;
图4是本发明实施例中典型应用场景示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步详细的说明。
一种高分辨率图像重建方法,该方法的主要思路为采集多帧同一场景的低分辨率图像,将所述的低分辨率图像通过配准方法映射到高分辨图像栅格中,之后再通过信息融合的方式将采集的低分辨率图像像素转换为高分辨率图像像素,进而得到高分辨率图像。下面分别就采集低分辨率图像、低分辨率图像配准、低分辨率图像映射到高分辨率图像栅格和图像融合进行依次说明,具体步骤如下:
1.采集低分辨率图像:通过摄像头采集同一场景多帧的低分辨率图像,在采集的过程中使摄像头具有微位移,这样获取得到的图像之间信息可互补,所提供的互补信息越多,获得的高分辨率图像的有效信息就越多。采集图像越多,提供的信息就越丰富,但由于受到分辨率放大倍数的限制,当图像帧数达到一定数量后,信息冗余度就会升高。因此,本实施例中设置图像帧数是分辨率放大倍数的2倍或者4倍。本实施例中采集的低分辨率图像帧数仅为优选的实施方式,作为其他实施方式,也可以根据实际情况选取其他数量的低分辨率图像。
2.低分辨率图像配准:从采集的低分辨率图像中选取一幅图像作为基准图像,本实施例中选取的为采集的第一帧图像;然后,对其余帧图像采用基于光流的方法,计算每一帧图像相对于基准图像的偏移量,本发明的偏移量包括位置偏移量和角度偏移量,位置偏移量又包括水平偏移值和垂直偏移值,实现对其余每副低分辨率图像的配准,所述的配准方法需要达到亚像素级的配准精度。
用(xb,yb)表示基准图像的像素点,用(xa,ya)表示其余帧的低分辨率图像的像素点。本实施例中利用位置偏移量对其余帧的低分辨率图像的像素点进行配准,则其余帧的低分辨率图像的像素点配准后得到的像素点为(xa+xd,ya+yd)。xd和yd分别表示该其余低分辨率图像帧相对于基准帧图像在X轴和Y轴的偏移量。
3.低分辨率图像映射到高分辨率图像栅格:
设置分辨率的放大倍数f(本实施例中设置放大倍数为2),将低分辨率图像中的像素映射到高分辨率图像栅格中,具体方法如下:
如图1所示,以基准帧图像为基础,将基准帧图像中的像素点(xb,yb)乘以f倍,即而得到图像点为(f·xb,f·yb);然后,将得到的图像点(f·xb,f·yb)赋到该栅格上,从而该基准帧图像转换成一幅高分辨率栅格图像。
如图2所示,将配准后的其余低分辨率图像进行上述同样转换并映射到对应高分辨率栅格上,得到映射后的图像点为(f·xa+f·xd,f·ya+f·yd)。
4.图像数据融合:采用保边滤波器的方法对上述映射到高分辨率栅格顶点像素值进行融合,从而实现将上述低分辨率图像像素转换成高分辨率图像。本实施例中采用的保边滤波器为双边滤波器,即在高分辨率图像栅格上的每个像素点上使用双边滤波器计算其像素值,表达式如下:
g ( i , j ) = Σ k , l f ( k , l ) w ( i , j , k , l ) Σ k , l w ( i , j , k , l ) - - - ( 1 )
其中,g(i,j)表示通过融合后得到的高分辨率图像的像素值;f(k,l)表示映射到各个栅格上的像素点的像素值;w(i,j,k,l)表示对应f(k,l)的权重系数。
其中,w(i,j,k,l)值的大小取决于定义域核d(i,j,k,l)和值域核r(i,j,k,l)的乘积,定义域核d(i,j,k,l)和值域核r(i,j,k,l)的表达式分别如下;
d ( i , j , k , l ) = exp ( - ( i - k ) 2 + ( j - l ) 2 2 σ d 2 ) - - - ( 2 )
r ( i , j , k , l ) = exp ( - | | f ( i , j ) - f ( k , l ) | | 2 2 σ r 2 ) - - - ( 3 )
则计算w(i,j,k,l)的表达式为:
ω ( i , j , k , l ) = exp ( - ( i - k ) 2 + ( j - l ) 2 2 σ d 2 - | | f ( i , j ) - f ( k , l ) | | 2 2 σ r 2 ) - - - ( 4 )
其中,上式(4)中邻域大小可根据实际情况取值,通常选取的典型值为3、5和7。
本实施例中还可以对低分辨率图像配准方法上选用基于光流的方法,在对低分辨率图像点进行融合时利用双边滤波器;上述方法仅为优选的实施方式而已,并不是对本发明的限制,作为其他实施方式还可以选择基于直方图投影或基于特征的方法进行配准,对低分辨率图像的融合方式也可以利用非局部均值滤波器(NLM)等方法。

Claims (5)

1.一种高分辨率图像重建方法,其特征在于,步骤如下:
步骤1,在同一场景采集两帧以上的低分辨率图像,从采集的图像中选取一帧图像作为基准帧图像;计算其余帧的低分辨率图像的各像素点相对于基准帧图像的偏移量,根据该偏移量将余帧的低分辨率图像以基准帧图像为基准进行配准;
步骤2,设定分辨率放大倍数,以基准帧图像各像素点乘以所述分辨率放大倍数为目标位置,将基准帧图像的像素点映射到待复原的高分辨率图像栅格上;将根据偏移量配准后的其余帧的低分辨率图像的像素点乘以所述分辨率放大倍数为目标位置,将配准后的其余帧的低分辨率图像的像素点映射到待复原的高分辨率图像栅格上;
步骤3,对步骤2中映射到各个栅格上的像素点的像素值通过融合计算得到各个栅格顶点的像素值,从而得到该场景的高分辨率图像。
2.根据权利要求1所述的一种高分辨率图像重建方法,其特征在于,步骤1中,其余帧的低分辨率图像的各像素点相对于基准帧图像的偏移量包括:位置偏移量和角度偏移量;位置偏移量又包括:水平偏移值和垂直偏移值。
3.根据权利要求1所述的一种高分辨率图像重建方法,其特征在于,采集的低分辨率图像数量的为分辨率放大倍数的2倍到4倍。
4.根据权利要求2所述的一种高分辨率图像重建方法,其特征在于,步骤2中用(xb,yb)表示基准帧图像的像素点;f表示设定的分辨率放大倍数;则基准帧图像的像素点映射到待复原的高分辨率图像栅格上后得到的像素点为(f·xb,f·yb);
用(xa,ya)表示其余的低分辨率图像的像素点;用xd和yd分别表示该其余低分辨率图像相对于基准帧图像在X轴和Y轴的位置偏移量;其余帧的低分辨率图像的像素点映射到待复原的高分辨率图像栅格上后得到的像素点为(f·xa+f·xd,f·ya+f·yd)。
5.根据权利要求4所述的一种高分辨率图像重建方法,其特征在于,步骤3中利用双边滤波器对映射到各个栅格上的像素点的像素值进行融合,表达式如下:
g ( i , j ) = Σ k , l f ( k , l ) w ( i , j , k , l ) Σ k , l w ( i , j , k , l )
其中,g(i,j)表示通过融合后得到的高分辨率图像的像素值;f(k,l)表示映射到各个栅格上的像素点的像素值;w(i,j,k,l)表示对应f(k,l)的权重系数。
CN201610234555.1A 2016-04-15 2016-04-15 一种高分辨率图像重建方法 Active CN105931210B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610234555.1A CN105931210B (zh) 2016-04-15 2016-04-15 一种高分辨率图像重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610234555.1A CN105931210B (zh) 2016-04-15 2016-04-15 一种高分辨率图像重建方法

Publications (2)

Publication Number Publication Date
CN105931210A true CN105931210A (zh) 2016-09-07
CN105931210B CN105931210B (zh) 2019-10-15

Family

ID=56839131

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610234555.1A Active CN105931210B (zh) 2016-04-15 2016-04-15 一种高分辨率图像重建方法

Country Status (1)

Country Link
CN (1) CN105931210B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106686310A (zh) * 2017-01-13 2017-05-17 浪潮(苏州)金融技术服务有限公司 一种动态提升摄像设备实际效果的方法及系统
CN108205804A (zh) * 2016-12-16 2018-06-26 阿里巴巴集团控股有限公司 图像处理方法、装置及电子设备
CN108596855A (zh) * 2018-04-28 2018-09-28 国信优易数据有限公司 一种视频图像质量增强方法、装置以及视频画质增强方法
CN110827224A (zh) * 2019-11-12 2020-02-21 南京工程学院 一种高分辨率高温锻件形状检测方法
CN113284077A (zh) * 2020-02-19 2021-08-20 华为技术有限公司 图像处理方法、装置、通信设备及可读存储介质
WO2023279920A1 (zh) * 2021-07-05 2023-01-12 腾讯科技(深圳)有限公司 基于显微镜的超分辨率方法、装置、设备及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050013509A1 (en) * 2003-07-16 2005-01-20 Ramin Samadani High resolution image reconstruction
US20100086227A1 (en) * 2008-10-04 2010-04-08 Microsoft Corporation Image super-resolution using gradient profile prior
CN101980289A (zh) * 2010-10-25 2011-02-23 上海大学 基于频域配准和凸集投影的多帧图像超分辨率重建方法
CN102651127A (zh) * 2012-04-01 2012-08-29 深圳市万兴软件有限公司 一种超分辨率重建的图像处理方法及系统
CN104933678A (zh) * 2015-06-30 2015-09-23 西安理工大学 一种基于像素强度的图像超分辨率重建方法
CN105279746A (zh) * 2014-05-30 2016-01-27 西安电子科技大学 一种基于双边滤波的多曝光图像融合方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050013509A1 (en) * 2003-07-16 2005-01-20 Ramin Samadani High resolution image reconstruction
US20100086227A1 (en) * 2008-10-04 2010-04-08 Microsoft Corporation Image super-resolution using gradient profile prior
CN101980289A (zh) * 2010-10-25 2011-02-23 上海大学 基于频域配准和凸集投影的多帧图像超分辨率重建方法
CN102651127A (zh) * 2012-04-01 2012-08-29 深圳市万兴软件有限公司 一种超分辨率重建的图像处理方法及系统
CN105279746A (zh) * 2014-05-30 2016-01-27 西安电子科技大学 一种基于双边滤波的多曝光图像融合方法
CN104933678A (zh) * 2015-06-30 2015-09-23 西安理工大学 一种基于像素强度的图像超分辨率重建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李展 等: "多分辨率图像序列的超分辨率重建", 《自动化学报》 *
王静 等: "多帧图像超分辨率重建算法", 《天津职业技术师范大学学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108205804A (zh) * 2016-12-16 2018-06-26 阿里巴巴集团控股有限公司 图像处理方法、装置及电子设备
CN108205804B (zh) * 2016-12-16 2022-05-31 斑马智行网络(香港)有限公司 图像处理方法、装置及电子设备
CN106686310A (zh) * 2017-01-13 2017-05-17 浪潮(苏州)金融技术服务有限公司 一种动态提升摄像设备实际效果的方法及系统
CN108596855A (zh) * 2018-04-28 2018-09-28 国信优易数据有限公司 一种视频图像质量增强方法、装置以及视频画质增强方法
CN110827224A (zh) * 2019-11-12 2020-02-21 南京工程学院 一种高分辨率高温锻件形状检测方法
CN113284077A (zh) * 2020-02-19 2021-08-20 华为技术有限公司 图像处理方法、装置、通信设备及可读存储介质
WO2021164329A1 (zh) * 2020-02-19 2021-08-26 华为技术有限公司 图像处理方法、装置、通信设备及可读存储介质
WO2023279920A1 (zh) * 2021-07-05 2023-01-12 腾讯科技(深圳)有限公司 基于显微镜的超分辨率方法、装置、设备及介质

Also Published As

Publication number Publication date
CN105931210B (zh) 2019-10-15

Similar Documents

Publication Publication Date Title
CN105931210A (zh) 一种高分辨率图像重建方法
CN111080724B (zh) 一种红外和可见光的融合方法
CN107959805B (zh) 基于混合相机阵列的光场视频成像系统及视频处理方法
CN109872401B (zh) 一种无人机视频增强现实实现方法
CN204090039U (zh) 一体化大场景全景视频监控装置
CN103606132A (zh) 基于空域和时域联合滤波的多帧数字图像去噪方法
CN107270910A (zh) 单镜头倾斜摄影航飞线路设计方法、系统及航拍摄影方法
CN103413286A (zh) 一种基于学习的高动态范围与高分辨率图像联合重建方法
CN102122387B (zh) 一种鲁棒的超分辨率图像重建方法
CN101605209A (zh) 摄像装置及图像再生装置
CN103546665A (zh) 基于多摄像头的图像清晰度增强方法
CN103955583B (zh) 一种确定夜晚灯光数据提取城市建成区阈值的方法
CN105005964A (zh) 基于视频序列影像的地理场景全景图快速生成方法
CN110662014B (zh) 一种光场相机四维数据大景深三维显示的方法
CN106254780A (zh) 一种双摄像头拍照控制方法、拍照控制装置及终端
CN106548494A (zh) 一种基于场景样本库的影视图像深度提取方法
CN105069748A (zh) 一种基于微小卫星物方扫描技术获取高分辨率图像的方法
TWI444593B (zh) 地面目標定位系統與方法
CN102930518A (zh) 基于改进的稀疏表示的图像超分辨率方法
CN105606123A (zh) 一种低空航空摄影测量自动纠正数字地面高程模型的方法
CN103808265A (zh) 油菜叶片及菌核病斑形态同步测量的方法、装置及系统
CN109035365B (zh) 一种高分辨率影像的镶嵌处理方法
CN104796623A (zh) 基于金字塔块匹配和泛函优化的拼接视频去结构偏差方法
CN104966277B (zh) 一种基于长短曝光图像对的图像复原方法
CN107330856B (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
GR01 Patent grant
GR01 Patent grant