CN103793879A - 数位影像的反扭曲处理方法 - Google Patents
数位影像的反扭曲处理方法 Download PDFInfo
- Publication number
- CN103793879A CN103793879A CN201310429402.9A CN201310429402A CN103793879A CN 103793879 A CN103793879 A CN 103793879A CN 201310429402 A CN201310429402 A CN 201310429402A CN 103793879 A CN103793879 A CN 103793879A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- block
- value
- transformation parameter
- triangle
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title abstract 2
- 238000012545 processing Methods 0.000 claims abstract description 84
- 238000006243 chemical reaction Methods 0.000 claims abstract description 30
- 230000009466 transformation Effects 0.000 claims description 79
- 238000003384 imaging method Methods 0.000 claims description 70
- 238000005520 cutting process Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 13
- 239000006185 dispersion Substances 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/04—Context-preserving transformations, e.g. by using an importance map
- G06T3/047—Fisheye or wide-angle transformations
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
一种数位影像的反扭曲处理方法,包括将第一影像视窗切分多个第一区块;根据这些第一区块的数量与位置将第二影像视窗的欲反扭曲区域划分多个第二区块;以第一区块的一组第一顶点坐标值与对应位置的第二区块的一组第二顶点坐标值进行坐标转换处理,得到坐标转换参数;根据坐标转换参数对第一区块中的每一个坐标值进行坐标转换处理,产生对应于第二区块的推估坐标值;根据推估坐标值从第二影像视窗中取得像素内容值,并将像素内容值运算后写入第一区块的坐标值的像素;重复所有这些第一区块与这些坐标值,直至完成第一影像的像素输出。
Description
技术领域
一种数位影像的处理方法,特別有关于一种数位影像的反扭曲处理方法。
背景技术
随着监控设备的普及化,越来越多人会装设监控设备来保障自身的财产与人身安全。为能无死角的监控,因此有厂商提出可全景(panorama)拍摄的摄影机。而全景摄影机的实现方式可以由多组镜头所接合而成,也可以利用以单一个鱼眼镜头(fisheye lens)而实现。全景拍摄的优点在于可以全方位且无死角的方式观看所拍摄的数位影像。
而透过全景拍摄的数位影像必然会有画面扭曲的情況。因此使用者难以从扭曲的画面中观看特定区域。所以为能提供使用者方便观看,需要将数位影像进行反扭曲(de-warp)的处理。习知的反扭曲处理是将数位影像的每一像素
(pixel)依序的进行坐标转换,并且于转换中对非整数坐标进行插补(interpolation)处理。而反扭曲的处理需要用到大量的计算,因此以嵌入式***(embedded system)搭配单芯片(system on a chip,SoC)设计的摄像装置恐难负担反扭曲处理时所需的运算量。
若欲即时的对全景拍摄的数位影像进行反扭曲处理,则需要将摄像装置连接于其他具有高效能的计算机装置,如桌上型电脑(PC)等。再由连接的高效能计算机装置对全景拍摄的数位影像进行反扭曲处理。而这些的繁复的影像处理势必会增加监控***的设置成本,更重要的,会影响***的可扩充性(scalability)。
发明内容
鉴于以上的问题,本发明在于提供一种数位影像的反扭曲处理方法。
本发明所提供的数位影像的反扭曲处理方法包括以下步骤:将第一影像视窗切分为多个第一区块;根据这些第一区块的数量与位置将第二影像视窗的欲反扭曲区域划分多个第二区块;以第一区块的一组第一顶点坐标值与所对应位置的第二区块的一组第二顶点坐标值进行坐标转换处理,用以得到坐标转换参数;根据坐标转换参数对第一区块中的每一个坐标值进行坐标转换处理,并产生对应于第二区块的推估坐标值;根据推估坐标值从第二影像视窗中取得像素内容值(如顏色、亮度等),并将像素内容值运算后(如interpolation)写入第一区块的坐标值的像素;重复所有这些第一区块与这些坐标值,直至完成第一影像的像素输出。
本发明另提出一种数位影像的反扭曲处理方法包括以下步骤:将第一影像视窗切分为多个第一区块;根据这些第一区块的数量与位置将第二影像视窗划分多个第二区块;以第一区块的一组第一顶点坐标值与所对应位置的第二区块的一组第二顶点坐标值进行坐标转换处理,用以得到坐标转换参数;重复每一第一区块与所对应第二区块的坐标转换处理,并将每一坐标转换参数依序记录于对应表中;将全景影像载入至第二影像视窗,并将全景影像划分成与这些第二区块相同数量的多个第三区块;将对应表的这些坐标转换参数导入并对相应位置的第一区块与第三区块进行坐标转换处理,用以产生第一区块内每一像素坐标对应于第三区块的推估坐标值;根据推估坐标值从第二影像视窗中取得像素内容值,并将像素内容值运算后(如interpolation)写入第一影像视窗中的坐标值的像素。
本发明所提出的数位影像的反扭曲处理方法可以应用在以全景拍摄装置拍摄时所产生的已扭曲影像的还原处理。本发明的处理方法可以降低反扭曲时的运算量,藉以加快数位影像在还原时的处理速度。
有关本发明的特征与实施例,兹配合图式作最佳实施例详细说明如下。
附图说明
图1为本发明的架构示意图;
图2为本发明的反扭曲处理流程示意图;
图3A为本发明的第一区块的切分示意图;
图3B为本发明的第二区块的切分示意图;
图4A为本发明的仿射变换的坐标转换处理流程示意图;
图4B为本发明的第一三角形与第二三角形的选择示意图;
图4C为本发明的透视变换的处理流程示意图;
图5为本发明的推估坐标值的处理流程示意图;
图6为本发明的另一实施例的流程示意图。
其中,附图标记:
数位影像处理装置100
影像输入单元110
处理单元120
储存单元130
反扭曲程序131
对应表132
第一三角形318
第三三角形319
第一影像视窗310
第一影像311
第一区块312
第二影像视窗320
第二影像321
第二区块322
第二三角形328
第四三角形329
具体实施方式
本发明除了应用于数位影像处理装置100之中,亦可透过计算机(computer)执行数位影像的反扭曲处理。而本说明书中是以数位影像处理装置100作为说明。本发明的数位影像处理装置100包括影像输入单元110、处理单元120与储存单元130。请参考图1所示,其为本发明的架构示意图。
处理单元120电性连接于影像输入单元110与储存单元130。影像输入单元110用以输入(或拍摄)数位影像。数位影像的种类为静态图像或动态影片。以数位影像处理装置100而言,影像输入单元110可能为感光元件,或其他可以拍摄全景影像的输入介面。本发明的数位影像由全景或鱼眼的方式拍摄,所以所拍摄的影像画面中被摄物与背景会有变形的情況。
为方便说明本发明的实施方式,因此在本发明中是以第一影像视窗与第二影像视窗作为显示第一影像与第二影像。但第一视窗与第二视窗不限于一定要显示给使用者观看,也可以拿来另做运算而不显示给使用者观看。而第一影像视窗与第二影像视窗的显示范围大小可以相同,也可以依照固定比例进行缩放。举例来说,第一影像视窗与第二影像视窗的大小可以为640*480的像素阵列大小,在本发明中是以第一影像视窗与第二影像视窗为相同解析度作为说明。
第一影像视窗用以播放已经过反扭曲处理后的第一影像。第二影像视窗为全景拍摄时所得到的数位影像。所以第二影像的被摄物(或背景)可能会有扭曲变形的情況。
储存单元130用以储存反扭曲程序131与对应表132。对应表132记录第一影像与第二影像间的多组坐标转换参数(其转换过程将于后文详述)。处理单元120根据第一影像与第二影像进行本发明的反扭曲处理。请参考图2所示,其为本发明的反扭曲处理流程示意图。本发明的反扭曲程序131包括以下步骤:
步骤S210:将第一影像视窗切分为多个第一区块;
步骤S220:根据这些第一区块的数量与位置将第二影像视窗的欲反扭曲区域划分多个第二区块;
步骤S230:选择第一区块与相应位置的第二区块;
步骤S240:以第一区块的一组第一顶点坐标值与所对应位置的第二区块的一组第二顶点坐标值进行坐标转换处理,用以得到坐标转换参数;
步骤S250:根据坐标转换参数对第一区块中的每一个坐标值进行坐标转换处理,并产生对应于第二区块的推估坐标值;
步骤S260:根据推估坐标值从第二影像视窗中取得像素内容值,并将像素内容值运算后写入第一区块的坐标值的像素;以及
步骤S270:重复所有这些第一区块与这些坐标值,直至完成第一影像的像素输出。
请参考图3A所示。首先,处理单元120将第一影像311的所在第一影像视窗310切分为多个第一区块312。第一影像视窗310的第一影像311为矩形,因此将第一影像311划分为多个更小面积的矩形。
由于第二影像321为全景影像,所以第二影像321在划分每一个第二区块322时,每一个第二区块322并非都是相同的四边形。此外,每一个第二区块322的四边并不一定是直线。例如:在第二影像321的边缘部分。由于第二影像321的边缘呈现圆弧型,所以切分为第二区块322时,这部分的第二区块322的边缘也是圆弧。而第二区块322的划分根据第一区块312的数量与位置,如图3B所示。
由于第二区块322为第二影像321视窗的一部分,因此第二区块322可以透过第二影像视窗320的坐标***进一步的界定各像素的坐标值。
如前文所述,第一区块312为矩形,所以第一区块312具有四个顶点坐标值。在此将第一区块312的各顶点坐标值定义为一组第一顶点坐标值。同理,对于第二区块322而言,也可以从第二影像视窗320中取得相应的一组第二顶点坐标值。接著,处理单元120依序选出第一区块312与相应的第二区块322,并根据下述方式进行坐标转换的处理。
坐标转换处理的实施例1:仿射变换(affine transform)
处理单元120以第一区块312的一组第一顶点坐标值与所对应位置的第二区块322的一组第二顶点坐标值进行坐标转换处理。首先,请参考图4A所示,其为本发明的仿射变换的坐标转换处理流程示意图。
步骤S411:根据第一区块的一组第一顶点坐标值的其中三个顶点坐标值决定第一三角形;
步骤S412:根据第一三角形的位置从第二区块的一组第二顶点坐标值中决定一第二三角形;以及
步骤S413:根据第一三角形的各顶点坐标值与相应位置的第二三角形的各顶点坐标值计算坐标转换参数。
由于第一区块312为矩形,所以可以将第一区块312划分为两个不相互重叠的三角形。在此可任一选择其中之一个三角形,并将其所选的定义为第一三角形318。处理单元120根据第一区块312从第二影像视窗320中取得相应位置上的第二区块322。处理单元120并根据第二区块322的第二顶点坐标值与第一三角形318的位置决定第二三角形328,请参考图4B所示。换句话说,若选择第一区块312的左上方部分为第一三角形318,则处理单元120也会对应的选择第二区块322的左上方部分为第二三角形328。
处理单元120根据第一三角形318的各顶点坐标值与相应位置的第二三角形328的各顶点坐标值进行仿射变换。请参考下式1所示,其为仿射变换的计算方式:
(Xi,Yi)为第二三角形328的顶点坐标值,(Xo,Yo)为第一三角形318的顶点坐标值,(a1,a2,a3,b1,b2,b3)为该坐标转换参数。
处理单元120可以将第一三角形318与第二三角形328的三个顶点坐标值带入仿射变换的式1中。如此一来,处理单元120可以得到六元一次方程组。处理单元120再根据方程式解出所有的坐标转换参数。处理单元120可以根据这些坐标转换参数进一步计算第一三角形318中的各像素映射(mapping)到第二三角形328的对应像素(此一处理将于后文进一步详述)。
第一区块312除了第一三角形318外,还另有第三三角形319。而第一三角形318与第三三角形319不相互重叠。第二区块322则包含第二三角形328与第四三角形329,第二三角形328对应于第四三角形329。
在本发明的一实施例中,直接将第一三角形318与第二三角形328间的坐标转换参数套用于第三三角形319与第四三角形329间。
在本发明的另一实施例中,处理单元120判断第二三角形328与第四三角形329的外型差异量是否大于外型门槛值。而外型差异量的判断方式可以藉由两三角形的面积差值所决定,但不限此方式。如果第二三角形328与第四三角形329的外型差异量小于一外型门槛值,则可将第四三角形329视为与第二三角形328相同,并将第一三角形318与第二三角形328间的坐标转换参数套用于第三三角形319与第四三角形329间。如果第二三角形328与第四三角形329的外型差异量大于外型门槛值,则再计算第二三角形328与第四三角形329间的坐标转换参数。
坐标转换处理的实施例2:透视变换(perspective transform)
除了上述的实施例外,本发明另提出透视变换处理方法。请参考图4C所示,其为本发明的透视变换的处理流程示意图。本发明的透视变换处理包括以下步骤:
步骤S421:取得第一区块的四个顶点坐标值;
步骤S422:取得相应位置的第二区块的四个顶点坐标值;以及
步骤S423:根据第一区块的四个顶点坐标值与第二区块的四个顶点坐标值,计算坐标转换参数。
处理单元120取得第一区块312的四个顶点坐标值与第二区块322的四个顶点坐标值。处理单元120根据第一区块312与第二区块322的各顶点坐标值进行下式2的计算:
(Xi,Yi)为第二区块322的一顶点坐标值,(Xo,Yo)为第一区块312的一顶点坐标值,(a1,a2,a3,b1,b2,b3,c1,c2)为该坐标转换参数。
在此将第一区块312的各顶点坐标值定义为一组第一顶点坐标值。同理,对于第二区块322而言,也可以从第二影像视窗320中取得相应的一组第二顶点坐标值。将第一顶点坐标值与第二顶点坐标值带入式2后,可以得到八元一次方程组,并从该些方程式中解出坐标转换参数。处理单元120再根据方程式解出所有的坐标转换参数。
在取得坐标转换参数后,处理单元120可以将第一区块312中的每一坐标值进行坐标转换处理。而经过转换后,可以得到第一区块312的每一坐标值是对应于第二影像视窗320中的哪一个位置。在此将在第二影像视窗320中所得到的坐标定义为推估坐标值。
处理单元120根据推估坐标值,并取得该推估坐标值的像素值。而由于数位影像中的像素位置为整数。例如:第二影像视窗320中的像素的坐标值可能为(100,100)或(123,321)。但由于推估坐标值的计算方式,使得计算结果可能并非整数。所以在本发明中若计算出来的坐标值不是整数时,则取最接近的邻近点(像素)的像素内容值,并将这几个像素进行计算,使得最后的计算结果输出成单一的整数。最后才将此一整数的结果写入第一区块312。换句话说,本发明针对推估坐标值在非整数的情況提出以下的处理。请参考图5所示,其为本发明的推估坐标值的处理流程示意图。
步骤S510:判断推估坐标值是否为整数;
步骤S520:若推估坐标值不为整数,根据推估坐标值取得上限坐标与下限坐标;
步骤S530:从第二影像视窗中取得上限坐标与下限坐标的多个像素值;
步骤S540:将所有像素值进行加权均值处理,并将像素内容值写入第一影像视窗的坐标值的像素;以及
步骤S550:若推估坐标值为整数,则以推估坐标值查找第二影像视窗的像素内容值,并将像素内容值写入第一影像视窗的坐标值的像素。
处理单元120得到推估坐标值后,处理单元120会判断推估坐标值是否为整数。当推估坐标值为整数时,则处理单元120则直接调用第二影像视窗320中推估坐标值上的像素值。处理单元120将该点的像素值写入第一区块312的对应位置中。
当推估坐标值不为整数时,处理单元120根据推估坐标值取得上限坐标与下限坐标。而上限坐标的选取方式是将推估坐标值的数值以上限(ceiling)的方式,选取接近且大于该推估坐标值的整数部分。对于下限坐标的选取的方式是将推估坐标值的数值以下限(floor)的方式,则是选取接近且小于该推估坐标值的整数部分。
举例来说,从第一区块312中选择一坐标点(100,100),意即Xo=100,Yo=100。并假设推估坐标值所计算的结果是Xi=101.79、Yi=102.13。因此推估坐标值并非整数,所以无法从该推估坐标值取得像素值。所以处理单元120会对推估坐标值进行上限、下限的处理。而推估坐标值由二维阵列的方式所构成,此一例子中的推估坐标值具有X、Y轴的两组数字。所以推估坐标值经过上限、下限的处理与排列后将会得到(101,102)、(101,103)、(102,102)、(102,103)四种组合。处理单元120会先从第二影像视窗320的四个坐标上取得该点的像素值。接下来,处理单元120会将四个像素内容值进行加权均值运算,并将均值运算所得到像素内容值写到第一区块312的相应位置上。上述的均值运算,在本实施方式中为双线性内插法(bilinear interpolation),但不限定于此法。
处理单元120根据前文所述的方式再对其他第一影像视窗310的其他坐标进行相应的处理,直至完成第一影像视窗310中的所有像素为止。此一实施例可以由处理单元120即时的将所输入的影像进行反扭曲处理。此一实施例主要被应用在第二影像321所欲反扭曲的范围会改变的情況。一般而言,使用者可以对全景摄影机进行电子式PTZ操作:左右转动(pan)、调整往上(或下)或缩放(zoom in/out)。特別是缩放的过程中,第二影像321中的第二区块322的比例也会随之改变。因此可以藉由本实施例即时的演算出对应的第一影像视窗310。
而除了前述的实施例外,本发明另对固定欲反扭曲范围的第二影像321提出其他的变化例。请参考图6所示,其为本发明的另一实施例的流程示意图。
步骤S610:将第一影像视窗切分为多个第一区块;
步骤S620:根据这些第一区块的数量与位置将第二影像视窗划分多个第二区块;
步骤S630:以第一区块的一组第一顶点坐标值与所对应位置的第二区块的一组第二顶点坐标值进行坐标转换处理,用以得到坐标转换参数;
步骤S640:重复每一第一区块与所对应第二区块的坐标转换处理,并将每一坐标转换参数依序记录于对应表中;
步骤S650:将全景影像载入至第二影像视窗,并将全景影像划分成与这些第二区块相同数量的多个第三区块;
步骤S660:将对应表的这些坐标转换参数导入并对相应位置的第一区块与第三区块进行坐标转换处理,用以产生第一区块内每一像素坐标对应于第三区块的推估坐标值;以及
步骤S670:根据推估坐标值从第二影像视窗中取得像素内容值,并将像素内容值运算后写入第一影像视窗中的坐标值的像素。
处理单元120根据前一实施例所述的方式进行第一区块312的划分。接着,处理单元120再计算第一区块312与第二区块322的坐标转换参数。当处理单元120完成第一区块312的坐标转换参数后,处理单元120会写入对应表132中。
在完成对应表132后,处理单元120会将对应表132储存至储存单元130中。接著,影像输入单元110另外拍摄新的影像时(将其定义为全景影像),处理单元120将全景影像载入至第二影像视窗320,并将全景影像划分成与这些第二区块322相同数量的多个第三区块。处理单元120将对应表132的这些坐标转换参数导入并对相应位置的第一区块312与第三区块进行坐标转换处理,用以产生第一区块312内每一像素坐标对应于第三区块的推估坐标值。处理单元120根据推估坐标值从第二影像视窗320中取得像素内容值,并将像素内容值运算后(如bilinear interpolation)写入第一影像视窗310中的坐标值的像素。重复对每一第一区块312的计算,直至所有第一影像视窗310的像素被填满为止。
上述各实施例,倘若第二区块322的变形程度太过严重或涵盖面积过大,会使得所得到的推估坐标值与实际的坐标值误差太大,而影响反扭曲的效果。我们可藉由缩小第一区块312的面积,来减少这种误差。参考图4B所示,在此定义包围第二区块322的最小矩形为第四区块。第四区块的面积,可做为判断第二区块322变形程度与涵盖面积的依据。可设定一面积门槛值,若第四区块的面积大于此门槛值,则将第一区块312切割为至少两块矩形,其中,切割的方式可以从水平轴向进行切割,也可以由垂直轴向进行切割。但切割的原则是子区块需包含第一区块312的两个顶点坐标值。如此持续切割,使第四区块的面积不大于门槛值为止,再重新计算坐标转换参数。
为避免上述切割步骤次数的不确定性造成问题,可定义一切割次数门槛值。若切割次数已大于等于该切割次数门槛值,则停止切割步骤,再重新计算坐标转换参数。
本发明所提出的数位影像的反扭曲处理可以应用在以全景拍摄装置时所产生的已扭曲影像的还原处理。本发明的处理方法可以降低反扭曲时的运算量,藉以加快数位影像在还原时的处理速度。
Claims (14)
1.一种数位影像的反扭曲处理方法,用以将一第一影像输出于一第一影像视窗中,该第一影像根据一第二影像视窗中的一第二影像的部分区域所输出,其特征在于,该反扭曲处理方法步骤包括:
将该第一影像视窗切分为多个第一区块;
根据该些第一区块的数量与位置将该第二影像视窗划分多个第二区块;
以该第一区块的一组第一顶点坐标值与所对应位置的该第二区块的一组第二顶点坐标值进行一坐标转换参数的运算,用以得到一坐标转换参数;
根据该坐标转换参数对该第一区块中的每一坐标值进行一坐标转换处理,并产生对应该第二区块的一推估坐标值;
根据该推估坐标值从该第二影像视窗中取得一像素内容值,并将该像素内容值写入该第一区块的一坐标值的像素;以及
重复所有该些第一区块与该些坐标值,直至完成该第一影像的像素输出。
2.如权利要求1所述的数位影像的反扭曲处理方法,其特征在于,所述划分该第一区块与该第二区块的步骤后更包括:
将包围该第二区块的范围定义为一第四区块,比较该第四区块的面积是否超过一门槛值;
若面积未超过该门槛值,则计算该第一区块与该第二区块的该坐标转换参数;以及
若面积超过该门槛值,则将该第一区块切割成至少两个子区块,且包围所有的该些子区块范围的面积不超过该门槛值。
3.如权利要求2所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:以该些子区块进行该坐标转换参数的运算,并产生相应的该坐标转换参数。
4.如权利要求1所述的数位影像的反扭曲处理方法,其特征在于,所述进行该坐标转换参数的运算的步骤包括:
根据该第一区块的该组第一顶点坐标值决定一第一三角形;
根据该第一三角形的位置从该第二区块的该组第二顶点坐标值中决定一第二三角形;以及
将该第一三角形的各顶点坐标值与相应位置的该第二三角形的各顶点坐标值进行仿射变换算式,用以进行该坐标转换参数的运算;
其中,根据如下仿射变换算式进行该坐标转换参数的运算
(Xi,Yi)为该第二三角形的顶点坐标值,(Xo,Yo)为该第一三角形的顶点坐标值,(a1,a2,a3,b1,b2,b3)为该坐标转换参数。
5.如权利要求4所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:
从该第一区块选择一第三三角形,且该第一三角形不重叠于该第三三角形;以及
根据该第三三角形的各顶点坐标值与相应位置的一第四三角形的各顶点坐标值进行仿射变换计算,用以计算该坐标转换参数;
其中,根据如下仿射变换算式计算该坐标转换参数
(Xi,Yi)为该第四三角形的顶点坐标,(Xo,Yo)为该第三三角形的顶点坐标,(a1,a2,a3,b1,b2,b3)为该坐标转换参数。
6.如权利要求1所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:
根据该第一区块的各顶点坐标值与相应位置的该第二区块的各顶点坐标值计算该坐标转换参数;
其中,根据下式计算该坐标转换参数,
(Xi,Yi)为该第二区块的顶点坐标值,(Xo,Yo)为该第一区块的顶点坐标值,(a1,a2,b1,b2,c1,c2)为该坐标转换参数。
7.如权利要求1所述的数位影像的反扭曲处理方法,其特征在于,所述在取得该像素内容值的步骤包括:
判断该推估坐标值是否为整数;
若该推估坐标值不为整数,根据该推估坐标值取得一上限坐标与一下限坐标;
从该第二影像视窗中取得该上限坐标与该下限坐标的多个像素值;
将该些像素内容值进行加权均值处理,并将该处理后的像素内容值写入该第一影像视窗的该坐标值的像素;以及
若该推估坐标值为整数,则以该推估坐标值查找该第二影像视窗的该像素内容值,并将该像素内容值写入该第一影像视窗的该坐标值的像素。
8.一种数位影像的反扭曲处理方法,用以将一第一影像输出于一第一影像视窗中,该第一影像根据一第二影像视窗中的一第二影像的部分区域所输出,其特征在于,该反扭曲处理方法包括:
将该第一影像视窗切分为多个第一区块;
根据该些第一区块的数量与位置将该第二影像视窗划分多个第二区块;
以该第一区块的一组第一顶点坐标值与所对应位置的该第二区块的一组第二顶点坐标值进行一坐标转换参数的处理,用以得到一坐标转换参数;
重复每一该第一区块的该坐标转换处理,并将每一该坐标转换参数依序记录于一对应表中;
将一全景影像载入至该第二影像视窗,并将该全景影像划分成与该些第二区块相同数量的多个第三区块;
将该对应表的该些坐标转换参数导入并对相应位置的该第一区块与该第三区块进行该坐标转换处理,用以产生该第一区块内每一像素坐标对应于该第三区块的一推估坐标值;以及
根据该推估坐标值从该第二影像视窗中取得一像素内容值,并将该像素内容值写入该第一影像视窗中的该坐标值的像素。
9.如权利要求8所述的数位影像的反扭曲处理方法,其特征在于,所述划分该第一区块与该第二区块的步骤后更包括:
将包围该第二区块的范围定义为一第四区块,比较该第四区块的面积是否超过一门槛值;
若面积未超过该门槛值,则计算该第一区块与该第二区块的该坐标转换参数;以及
若面积超过该门槛值,则将该第一区块切割成至少两个子区块,且包围所有的该些子区块的面积不超过该门槛值。
10.如权利要求9所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:以该些子区块进行该坐标转换参数的处理,并产生相应的该坐标转换参数。
11.如权利要求8所述的数位影像的反扭曲处理方法,其特征在于,所述进行该坐标转换参数的处理的步骤包括:
根据该第一区块的该组第一顶点坐标值决定一第一三角形;
根据该第一三角形的位置从该第二区块的该组第二顶点坐标值中决定一第二三角形;以及
用该第一三角形的各顶点坐标值与相应位置的该第二三角形的各顶点坐标值进行仿射变换算式运算,用以进行该坐标转换参数计算;
其中,根据如下仿射变换算式进行该坐标转换参数计算
(Xi,Yi)为该第二三角形的顶点坐标值,(Xo,Yo)为该第一三角形的顶点坐标值,(a1,a2,a3,b1,b2,b3)为该坐标转换参数。
12.如权利要求8所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:
从该第一区块选择一第三三角形,且该第一三角形不重叠于该第三三角形;以及
根据该第三三角形的各顶点坐标值与相应位置的一第四三角形的各顶点坐标值进行仿射变换算式运算,用以进行该坐标转换参数计算;
其中,根据如下仿射变换算式进行该坐标转换参数计算
(Xi,Yi)为该第四三角形的顶点坐标值,(Xo,Yo)为该第三三角形的顶点坐标值,(a1,a2,a3,b1,b2,b3)为该坐标转换参数。
13.如权利要求8所述的数位影像的反扭曲处理方法,其特征在于,所述产生该坐标转换参数的步骤包括:
根据该第一区块的各顶点坐标值与相应位置的该第二区块的各顶点坐标值进行该坐标转换参数计算;
其中,根据下式进行该坐标转换参数计算,
(Xi,Yi)为该第二区块的顶点坐标值,(Xo,Yo)为该第一区块的顶点坐标值,(a1,a2,b1,b2,c1,c2)为该坐标转换参数。
14.如权利要求8所述的数位影像的反扭曲处理方法,其特征在于,所述在取得该像素内容值的步骤包括:
判断该推估坐标值是否为整数;
若该推估坐标值不为整数,根据该推估坐标值取得一上限坐标与一下限坐标;
从该第二影像视窗中取得该上限坐标与该下限坐标的多个像素值;
将该些像素内容值进行均值处理,并将处理后的像素内容值写入该第一影像视窗的该坐标值的像素;以及
若该推估坐标值为整数,则以该推估坐标值查找该第二影像视窗的该像素内容值,并将该像素内容值写入该第一影像视窗的该坐标值的像素。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW101140421A TWI476730B (zh) | 2012-10-31 | 2012-10-31 | 數位影像的反扭曲處理方法 |
TW101140421 | 2012-10-31 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103793879A true CN103793879A (zh) | 2014-05-14 |
Family
ID=49322236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310429402.9A Pending CN103793879A (zh) | 2012-10-31 | 2013-09-18 | 数位影像的反扭曲处理方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8942513B2 (zh) |
EP (1) | EP2728544B1 (zh) |
CN (1) | CN103793879A (zh) |
TW (1) | TWI476730B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9984436B1 (en) * | 2016-03-04 | 2018-05-29 | Scott Zhihao Chen | Method and system for real-time equirectangular projection |
US9547883B1 (en) | 2016-08-19 | 2017-01-17 | Intelligent Security Systems Corporation | Systems and methods for dewarping images |
US9609197B1 (en) | 2016-08-19 | 2017-03-28 | Intelligent Security Systems Corporation | Systems and methods for dewarping images |
US10186067B2 (en) * | 2016-10-25 | 2019-01-22 | Aspeed Technology Inc. | Method and apparatus for generating panoramic image with rotation, translation and warping process |
TWI599989B (zh) | 2016-11-29 | 2017-09-21 | 財團法人工業技術研究院 | 用於交通工具之影像處理方法及影像系統 |
US11210840B1 (en) | 2020-08-27 | 2021-12-28 | Aspeed Technology Inc. | Transform method for rendering post-rotation panoramic images |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5319744A (en) * | 1991-04-03 | 1994-06-07 | General Electric Company | Polygon fragmentation method of distortion correction in computer image generating systems |
CN1192312A (zh) * | 1995-07-27 | 1998-09-02 | 传感电子公司 | 适用于无活动机构摄像机的图像形成与处理装置及其方法 |
US6380936B1 (en) * | 1999-10-22 | 2002-04-30 | Microsoft Corportion | System and method for inferring projective mappings |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4862388A (en) * | 1986-12-15 | 1989-08-29 | General Electric Company | Dynamic comprehensive distortion correction in a real time imaging system |
IL79822A (en) * | 1985-12-19 | 1990-03-19 | Gen Electric | Method of comprehensive distortion correction for a computer image generation system |
US6819333B1 (en) * | 2000-05-12 | 2004-11-16 | Silicon Graphics, Inc. | System and method for displaying an image using display distortion correction |
WO2003027766A2 (en) * | 2001-09-27 | 2003-04-03 | Eyesee360, Inc. | System and method for panoramic imaging |
NZ525129A (en) * | 2003-10-03 | 2006-09-29 | Bruce Peter Parker | An improved transformation method for creating pre-distorted images to event surfaces of televised events |
US7242407B2 (en) * | 2004-05-28 | 2007-07-10 | Lockheed Martin Corporation | Reprojecting map images using graphical techniques |
US7768545B2 (en) * | 2007-03-06 | 2010-08-03 | Otto Gregory Glatt | Panoramic image management system and method |
TWI364725B (en) * | 2008-05-06 | 2012-05-21 | Primax Electronics Ltd | Video processing method and video processing system |
JP5075182B2 (ja) * | 2009-10-20 | 2012-11-14 | 楽天株式会社 | 画像処理装置、画像処理方法、および、画像処理プログラム |
TWI443604B (zh) * | 2010-11-01 | 2014-07-01 | Mstar Semiconductor Inc | 影像校正方法及影像校正裝置 |
-
2012
- 2012-10-31 TW TW101140421A patent/TWI476730B/zh active
-
2013
- 2013-03-18 US US13/846,213 patent/US8942513B2/en active Active
- 2013-09-18 CN CN201310429402.9A patent/CN103793879A/zh active Pending
- 2013-10-04 EP EP13187348.1A patent/EP2728544B1/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5319744A (en) * | 1991-04-03 | 1994-06-07 | General Electric Company | Polygon fragmentation method of distortion correction in computer image generating systems |
CN1192312A (zh) * | 1995-07-27 | 1998-09-02 | 传感电子公司 | 适用于无活动机构摄像机的图像形成与处理装置及其方法 |
US6380936B1 (en) * | 1999-10-22 | 2002-04-30 | Microsoft Corportion | System and method for inferring projective mappings |
Also Published As
Publication number | Publication date |
---|---|
TWI476730B (zh) | 2015-03-11 |
US20140119653A1 (en) | 2014-05-01 |
US8942513B2 (en) | 2015-01-27 |
EP2728544B1 (en) | 2018-10-31 |
TW201417044A (zh) | 2014-05-01 |
EP2728544A1 (en) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109618090B (zh) | 对通过使用广角透镜捕获的图像的图像失真校正的方法和*** | |
CN106652972B (zh) | 显示屏的处理电路、显示方法及显示器件 | |
CN103793879A (zh) | 数位影像的反扭曲处理方法 | |
US8849013B2 (en) | System and method of image processing | |
TWI441514B (zh) | 魚眼影像校正及減輕透視變形的影像處理方法及相關影像處理裝置 | |
CN109741463B (zh) | 虚拟现实场景的渲染方法、装置及设备 | |
US9959600B2 (en) | Motion image compensation method and device, display device | |
US9369627B2 (en) | Video processing apparatus, video processing system, and video processing method | |
CN113608641B (zh) | 曲面屏显示位置的调整方法、装置、智能音箱及存储介质 | |
CN108564551A (zh) | 鱼眼图像处理方法及鱼眼图像处理装置 | |
CN110740309A (zh) | 一种图像显示方法、装置、电子设备及存储介质 | |
CN108230454B (zh) | 一种全景图片的切图方法、装置及存储介质 | |
CN109545122B (zh) | Vr显示的补偿方法及补偿装置、显示*** | |
CN108875733B (zh) | 一种红外小目标快速提取*** | |
CN112711984B (zh) | 注视点定位方法、装置和电子设备 | |
US11636708B2 (en) | Face detection in spherical images | |
CN111625210B (zh) | 一种大屏控制方法、装置及设备 | |
CN104754201B (zh) | 一种电子设备及信息处理方法 | |
JP6103942B2 (ja) | 画像データ処理装置及び画像データ処理プログラム | |
CN114257733A (zh) | 具有视点偏移的全向图像的图像处理的方法和*** | |
JP2018041201A (ja) | 表示制御プログラム、表示制御方法および情報処理装置 | |
CN102447829B (zh) | 拍摄参数设置方法及*** | |
CN102065253A (zh) | 一种软硬件协同的图像水平缩放方法及装置 | |
CN109741465A (zh) | 图像处理方法及装置、显示装置 | |
JP5539561B2 (ja) | 画像処理装置およびその方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination |