CN107305681A - 一种用于光学场景图像的快速扩充方法 - Google Patents
一种用于光学场景图像的快速扩充方法 Download PDFInfo
- Publication number
- CN107305681A CN107305681A CN201610258381.2A CN201610258381A CN107305681A CN 107305681 A CN107305681 A CN 107305681A CN 201610258381 A CN201610258381 A CN 201610258381A CN 107305681 A CN107305681 A CN 107305681A
- Authority
- CN
- China
- Prior art keywords
- image
- expansion
- point
- clockwise
- expand
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003287 optical effect Effects 0.000 title claims abstract description 20
- 230000003190 augmentative effect Effects 0.000 claims abstract description 21
- 230000000007 visual effect Effects 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 6
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000013041 optical simulation Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Classifications
-
- G06T3/18—
Abstract
本发明公开了一种用于光学场景图像的快速扩充方法,该方法的步骤包括:获取实测图像,并确定该图像的扩充图像尺寸和扩充算子S1、计算图像中单个像素点的数值,并从内到外进行逐层扩充S2、分别获取顺时针和逆时针的扩充图像S3和对获取的顺时针扩充图像和逆时针扩充图像取平均值,获得最终扩充的图像S4。本方案解决了仿真试验中实测数据与光学成像场景生成器视场角不匹配时的应用问题,能够在不影响探测系统工作的前提下按照场景生成器视场角对实测数据进行快速扩充,保障了光学仿真试验的顺利进行,具有广阔的应用前景。
Description
技术领域
本发明涉及一种图像处理方法,特别是涉及一种用于光学场景图像的快速扩充方法。
背景技术
光学场景仿真系统需要加载实测数据对光学探测系统的性能进行测试。由于实测数据采集设备和光学成像场景生成器的视场角可能不同,有时实测数据并不能完全匹配光学成像场景生成器。当实测数据视场角过大时可以直接进行裁剪;但当实测数据视场角过小不能完全覆盖光学成像场景生成器时,需要在不影响现有数据的前提下对实测图像数据进行扩充。简单的单色背景填充方法可能会引起光学探测系统的误操作,因此需要结合实测数据信息进行外延式的扩充,使得实测数据的图像能够平滑的扩充至光学场景生成器的整个视场,同时保证新扩充的图像数据不会对光学探测系统产生影响。
发明内容
本发明要解决的技术问题是提供一种用于光学场景图像的快速扩充方法,解决实测数据与光学成像场景生成器视场角不匹配的问题,在不影响探测系统工作的前提下按照场景生成器视场角对实测数据进行有效扩充。
为解决上述技术问题,本发明采用下述技术方案:
一种用于光学场景图像的快速扩充方法,该方法的步骤包括:
S1、获取实测图像,并确定该图像的扩充图像尺寸和扩充算子;
S2、计算图像中单个像素点的数值,并从内到外进行逐层扩充;
S3、分别获取顺时针和逆时针的扩充图像;
S4、对获取的顺时针扩充图像和逆时针扩充图像取平均值,获得最终扩充的图像。
优选地,所述步骤S1包括:
S11、获取尺寸为m0×n0的实测图像;
S12、根据视场角与像素的转换关系确定扩充后图像的最大尺寸为m1×n1,该扩充后图像的图像中心区域为实测图像m0×n0的图像数据,其余部分为扩充数据;
S13、根据扩充图像的尺寸和像素,选用归一化扩充算子对图像进行扩充;
S14、为扩充图像准备临时图像数据Ia,Ib,Ic,临时图像的大小均为(m0+2d+8)×(n0+2d+8),其中,d为扩充层数,其值为不小于(m1-m0)/2和(n1-n0)/2的最小整数,临时图像数据Ia,Ib,Ic的中心m0×n0区域数值设置为I0,其余区域均为0.
优选地,该方法采用9阶服从高斯分布的归一化扩充算子。
优选地,所述步骤S2包括:
S21、将扩充算子中心移动到待计算的某个点,使算子作用到以该点为中心9乘9大小的区域上,将对应位置的数值相乘并累加可得:
S22、对上述公式进行补偿,可得:
S23、根据步骤S21和S22的计算公式,确定计算的某个点的数值:
优选地,所述步骤S3中所述顺时针扩充图像包括:
S31、以实测图像占据临时图像Ia的中心区域,以中心区域I0的左上角点在Ia上对应的坐标为起点,顺时针逐点计算数值,完成单层扩充;
S32、由内层至外层依次扩充,直至扩充至d层位置,获得顺时针扩充后的临时图像Ia。
优选地,所述步骤S3中所述逆时针扩充图像包括:
S33、以实测图像占据临时图像Ib的中心区域,以中心区域I0的左上角点在Ib上对应的坐标为起点,逆时针逐点计算数值,完成单层扩充;
S34、由内层至外层依次扩充,直至扩充至d层位置,获得逆时针扩充后的临时图像Ib。
优选地,所述步骤S4包括:
S41、将顺时针扩充后的临时图像Ia和逆时针扩充后的临时图像Ib取平均值:Iab=(Ia+Ib)/2;
S42、将扩充算子逐点作用到Iab上,可得:
S43、将Ic中心m0×n0区域数值设置为I0,并按照m1×n1进行裁剪,去掉边缘部分,最终得到扩充后的图像。
本发明的有益效果如下:
本发明所述技术方案,解决了仿真试验中实测数据与光学成像场景生成器视场角不匹配时的应用问题,能够在不影响探测系统工作的前提下按照场景生成器视场角对实测数据进行快速扩充,保障了光学仿真试验的顺利进行,具有广阔的应用前景。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出本发明所述扩充方法的方法流程图;
图2示出本发明所述单点数值计算的示意图;
图3示出本发明所述顺时针图像扩充的示意图。
附图标号
1、扩充算子,2、未扩充的图像点,3、已扩充/原有的图像点,4、待扩充。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
本发明公开了一种用于光学场景图像的快速扩充方法,该方法的步骤包括:
S1、获取实测图像,并确定该图像的扩充图像尺寸和扩充算子;所述步骤S1包括:
S11、获取尺寸为m0×n0的实测图像;
S12、根据视场角与像素的转换关系确定扩充后图像的最大尺寸为m1×n1,该扩充后图像的图像中心区域为实测图像m0×n0的图像数据,其余部分为扩充数据;
S13、根据扩充图像的尺寸和像素,选用9阶服从高斯分布的归一化扩充算子对图像进行扩充;
S14、为扩充图像准备临时图像数据Ia,Ib,Ic,临时图像的大小均为(m0+2d+8)×(n0+2d+8),其中,d为扩充层数,其值为不小于(m1-m0)/2和(n1-n0)/2的最小整数,临时图像数据Ia,Ib,Ic的中心m0×n0区域数值设置为I0,其余区域均为0。
S2、计算图像中单个像素点的数值,并从内到外进行逐层扩充;所述步骤S2包括:
S21、将扩充算子中心移动到待计算的某个点,使算子作用到以该点为中心9乘9大小的区域上,将对应位置的数值相乘并累加可得:
S22、对上述公式进行补偿,可得:
S23、根据步骤S21和S22的计算公式,确定计算的某个点的数值:
S3、分别获取顺时针和逆时针的扩充图像;所述步骤S3中所述顺时针扩充图像包括:
S31、以实测图像占据临时图像Ia的中心区域,以中心区域I0的左上角点在Ia上对应的坐标为起点,顺时针逐点计算数值,完成单层扩充;
S32、由内层至外层依次扩充,直至扩充至d层位置,获得顺时针扩充后的临时图像Ia。
所述步骤S3中所述逆时针扩充图像包括:
S33、以实测图像占据临时图像Ib的中心区域,以中心区域I0的左上角点在Ib上对应的坐标为起点,逆时针逐点计算数值,完成单层扩充;
S34、由内层至外层依次扩充,直至扩充至d层位置,获得逆时针扩充后的临时图像Ib。
S4、对获取的顺时针扩充图像和逆时针扩充图像取平均值,获得最终扩充的图像;所述步骤S4包括:
S41、将顺时针扩充后的临时图像Ia和逆时针扩充后的临时图像Ib取平均值:Iab=(Ia+Ib)/2;
S42、将扩充算子逐点作用到Iab上,可得:
S43、将Ic中心m0×n0区域数值设置为I0,并按照m1×n1进行裁剪,去掉边缘部分,最终得到扩充后的图像。
下面通过一组实施例对本发明做进一步说明:
本实例提供一种光学场景图像快速扩充方法,该方法的具体步骤为:
第一步 确定图像尺寸和扩充算子
假设实测图像尺寸为m0×n0,根据视场角与像素的转换关系确定扩充后的图像尺寸为m1×n1,其中图像中心m0×n0区域对应图像为实测数据,其余部分为扩充数据。设实测图像为I0,表示实测图像中第(i,j)个点的数据值。
选取9阶服从高斯分布的归一化扩充算子:
Pi,j表示算子中第(i,j)个点的数据值。
设定临时图像数据Ia,Ib,Ic,大小均为(m0+2d+8)×(n0+2d+8)。其中d为不小于(m1-m0)/2和(n1-n0)/2的最小整数,d的选取是为保证数据处理过程中图像向四个方向扩充的层数一致,d即为一共需要扩充的层数。将Ia,Ib,Ic中心m0×n0区域数值设置为I0,其余区域均为0。
第二步 计算单点数值
本发明采用从内到外的逐层扩充方法,在扩充过程中,对于临时图像Ia中待计算的某个点其周围8个点的数值不可能同时为0。将扩充算子中心移到该点,使算子作用到以该点为中心9乘9大小的区域上,将对应位置的数值相乘并累加可得:
由于周围9×9的区域中至少有一半是未扩充的零值,扩充算子作用到零值上便不再具有归一化的特性,因此需要对上述公式进行补偿,可得:
其中表示对取符号(零值为0,其他为1),P表示归一化扩充算子。最终得到临时图像Ia中单点(x,y)的数值:
如图2所示,单点的数值计算主要利用了以该点为中心点的模板区域内非零点的信息(图中非零点为模板内中心点左下方的40个点)。
第三步 获取顺时针扩充图像
在第一步获取的图像Ia中,原图像I0的信息占据了中间m0×n0的区域,I0左上角点在Ia上对应的坐标为从该点上方的数据点开始,按照“右-下-左-上”的顺时针顺序,采用第二步的方法进行四个方向的逐点数值计算,完成单层扩充。同理,按照上述方法由内到外依次完成d层扩从,最终得到图像Ia。顺时针扩充的整个扩充方式如图3所示。
第四步 获取逆时针扩充图像
在第一步获取的图像Ib中,原图像I0的信息占据了中间m0×n0的区域,I0左上角点在Ib上对应的坐标为从该点上方的数据点开始,按照“左-下-右-上”的逆时针顺序,采用第二步的方法进行四个方向的逐点数值计算,完成单层扩充。同理,按照上述方法由内到外依次完成d层扩从,最终得到图像Ib,按照逆时针的顺序逐层扩充图像,得到Ib。
第五步 获取最终扩充图像
将第三步顺时针扩充的图像Ia和第四步逆时针扩充的图像Ib取平均可得:
Iab=(Ia+Ib)/2
将扩充算子逐点作用到Iab上,可得:
将Ic中心m0×n0区域数值设置为I0,并按照m1×n1进行裁剪,去掉边缘部分,最终得到扩充后的图像I1。
上述步骤实现了光学场景图像的扩充。由于每个数据点的计算只用到了周边9×9区域内的数据点,并没有利用更多更深层次的图像信息,总体计算量较少。因此本发明所采用的方法是一种快速的图像信息外延式扩充方法。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (8)
1.一种用于光学场景图像的快速扩充方法,其特征在于,该方法的步骤包括:
S1、获取实测图像,并确定该图像的扩充图像尺寸和扩充算子;
S2、计算图像中单个像素点的数值,并从内到外进行逐层扩充;
S3、分别获取顺时针和逆时针的扩充图像;
S4、对获取的顺时针扩充图像和逆时针扩充图像取平均值,获得最终扩充的图像。
2.根据权利要求1所述的快速扩充方法,其特征在于,所述步骤S1包括:
S11、获取尺寸为m0×n0的实测图像,设为I0,其中表示实测图像中第(i,j)个点的数据值;
S12、根据视场角与像素的转换关系确定扩充后图像的最大尺寸为m1×n1,该扩充后图像的图像中心区域为实测图像m0×n0的图像数据,其余部分为扩充数据;
S13、根据扩充图像的尺寸和像素,选用归一化扩充算子对图像进行扩充;
S14、为扩充图像准备临时图像数据Ia,Ib,Ic,临时图像的大小均为(m0+2d+8)×(n0+2d+8),其中,d为扩充层数,临时图像数据Ia,Ib,Ic的中心m0×n0区域数值设置为I0,其余区域均为0。
3.根据权利要求2所述的快速扩充方法,其特征在于,该方法采用9阶服从高斯分布的归一化扩充算子。
4.根据权利要求2所述的快速扩充方法,其特征在于,所述临时图像的大小不小于(m1-m0)/2和(n1-n0)/2的最小整数。
5.根据权利要求1所述的快速扩充方法,其特征在于,所述步骤S2包括:
S21、将扩充算子中心移动到待计算的某个点,使算子作用到以该点为中心9乘9大小的区域上,将对应位置的数值相乘并累加可得:
S22、对上述公式进行补偿,可得:
S23、根据步骤S21和S22的计算公式,确定计算的某个点的数值:
6.根据权利要求2所述的快速扩充方法,其特征在于,所述步骤S3中所述顺时针扩充图像包括:
S31、以实测图像占据临时图像Ia的中心区域,中心区域I0的左上角点在Ia上对应的坐标为现以为起点,顺时针逐点计算数值,完成单层扩充;
S32、由内层至外层依次扩充,直至扩充至d层位置,获得顺时针扩充后的临时图像Ia。
7.根据权利要求6所述的快速扩充方法,其特征在于,所述步骤S3中所述逆时针扩充图像包括:
S33、以实测图像占据临时图像Ib的中心区域,中心区域I0的左上角点在Ib上对应的坐标为现以为起点,逆时针逐点计算数值,完成单层扩充;
S34、由内层至外层依次扩充,直至扩充至d层位置,获得逆时针扩充后的临时图像Ib。
8.根据权利要求7所述的快速扩充方法,其特征在于,所述步骤S4包括:
S41、将顺时针扩充后的临时图像Ia和逆时针扩充后的临时图像Ib取平均值:Iab=(Ia+Ib)/2;
S42、将扩充算子逐点作用到Iab上,可得:
S43、将Ic中心m0×n0区域数值设置为I0,并按照m1×n1进行裁剪,去掉边缘部分,最终得到扩充后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610258381.2A CN107305681B (zh) | 2016-04-22 | 2016-04-22 | 一种用于光学场景图像的快速扩充方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610258381.2A CN107305681B (zh) | 2016-04-22 | 2016-04-22 | 一种用于光学场景图像的快速扩充方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107305681A true CN107305681A (zh) | 2017-10-31 |
CN107305681B CN107305681B (zh) | 2020-09-18 |
Family
ID=60150537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610258381.2A Active CN107305681B (zh) | 2016-04-22 | 2016-04-22 | 一种用于光学场景图像的快速扩充方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107305681B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101496410A (zh) * | 2006-05-29 | 2009-07-29 | Opt株式会社 | 利用广角镜头的摄像数据的压缩方法、压缩装置、广角摄像装置以及监视系统 |
CN101556694A (zh) * | 2009-05-15 | 2009-10-14 | 广东工业大学 | 一种旋转图像的匹配方法 |
US20100054620A1 (en) * | 2008-08-27 | 2010-03-04 | Seiko Epson Corporation | Image processing apparatus, image processing method, and image processing program |
-
2016
- 2016-04-22 CN CN201610258381.2A patent/CN107305681B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101496410A (zh) * | 2006-05-29 | 2009-07-29 | Opt株式会社 | 利用广角镜头的摄像数据的压缩方法、压缩装置、广角摄像装置以及监视系统 |
US20100054620A1 (en) * | 2008-08-27 | 2010-03-04 | Seiko Epson Corporation | Image processing apparatus, image processing method, and image processing program |
CN101556694A (zh) * | 2009-05-15 | 2009-10-14 | 广东工业大学 | 一种旋转图像的匹配方法 |
Non-Patent Citations (1)
Title |
---|
杜常青 等: "二维纹理合成算法优化研究", 《计算机科学》 * |
Also Published As
Publication number | Publication date |
---|---|
CN107305681B (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020108358A1 (zh) | 图像修复方法、装置、计算机设备和存储介质 | |
ES2438141T3 (es) | Aparato para restaurar imágenes y método del mismo | |
CN100458846C (zh) | 一种图像修复方法 | |
CN104463786B (zh) | 一种移动机器人图像拼接方法及装置 | |
CN106709947A (zh) | 一种基于rgbd相机的三维人体快速建模系统 | |
CN104376535A (zh) | 一种基于样本的快速图像修复方法 | |
CN101777178A (zh) | 一种图像修复方法 | |
CN106023230B (zh) | 一种适合变形图像的稠密匹配方法 | |
CN107369204A (zh) | 一种基于深度学习从单幅照片恢复出场景基本三维结构的方法 | |
CN105303616A (zh) | 基于单张照片的浮雕建模方法 | |
CN105160686B (zh) | 一种基于改进sift算子的低空多视角遥感影像匹配方法 | |
CN106780370A (zh) | 一种图像去抖动设备及其方法 | |
CN103413339B (zh) | 十亿像素高动态范围图像重建与显示的方法 | |
CN111881985B (zh) | 立体匹配方法、装置、终端和存储介质 | |
CN106023089A (zh) | 一种基于块匹配的图像修复方法 | |
KR101593316B1 (ko) | 스테레오 카메라를 이용한 3차원 모델 재구성 방법 및 장치 | |
CN102521786B (zh) | 一种基于颜色检测和快速匹配法的照片水印去除方法 | |
CN108010082A (zh) | 一种几何匹配的方法 | |
CN106709883A (zh) | 基于联合双边滤波和尖锐特征骨架提取的点云去噪方法 | |
CN107301674A (zh) | 一种基于窗口插值的三维重建纹理恢复方法 | |
CN106846271A (zh) | 一种去除身份证照片中网纹的方法 | |
CN107945151A (zh) | 一种基于相似性变换的重定位图像质量评价方法 | |
CN113744142B (zh) | 图像修复方法、电子设备及存储介质 | |
CN107358581A (zh) | 快速图像修复方法 | |
CN104836953A (zh) | 多投影仪屏幕特征点自动拍照及去噪识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |