CN106780389A - 一种基于坐标转换的鱼眼图像校正方法及装置 - Google Patents
一种基于坐标转换的鱼眼图像校正方法及装置 Download PDFInfo
- Publication number
- CN106780389A CN106780389A CN201611206384.8A CN201611206384A CN106780389A CN 106780389 A CN106780389 A CN 106780389A CN 201611206384 A CN201611206384 A CN 201611206384A CN 106780389 A CN106780389 A CN 106780389A
- Authority
- CN
- China
- Prior art keywords
- fish eye
- eye images
- coordinate
- correction chart
- points
- 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
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 34
- 241000251468 Actinopterygii Species 0.000 claims abstract description 313
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- 230000003287 optical effect Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于坐标转换的鱼眼图像校正方法及装置,从鱼眼图像中任意选取一个像素点P'作为对鱼眼图像进行校正的中心点,通过鱼眼摄像机的成像模型计算出P′在鱼眼半球的入射点P的球面坐标,将平行于以P点为切点与鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标,将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。本发明的方法和装置对高清鱼眼摄相机进行多方位实时校正,便于观察鱼眼图像中各个不同区域的细节。
Description
技术领域
本发明属于视频监控领域,尤其涉及一种基于坐标转换的鱼眼图像校正方法及装置。
背景技术
在很多应用场景中,用户希望在一个视频监控摄像机的监控画面中能够监控到更大的范围,因此超大广角的镜头得到了越来越多的应用,其中,鱼眼镜头拥有能够获得超过180度视场角的图像,实现无死角监控,在视频监控中得到广泛应用。鱼眼摄像机虽然能获取超大视场角的图像,但是鱼眼摄像机拍摄的图像有着严重的畸变,不符合人眼的观察习惯,因此需要对鱼眼摄像机拍摄的图像进行畸变校正,以便人们对鱼眼摄像机的监控区域进行观察。基于透视原理的畸变校正技术可以将鱼眼摄像机拍摄的图像校正为无畸变的图像,目前基于透视原理的畸变校正一般根据球面投影模型,通过多项式求解对鱼眼摄像机拍摄的图像进行校正,这种方法计算复杂,计算量大,在实时应用中有较大的困难。
例如,专利CN201510195083公开了一种基于球面透视投影的鱼眼图像校正方法,该方法首先将鱼眼图像映射成球面图像,据此建立鱼眼图像与校正图像的关系,然后对空白的校正图像内的像素点采用向后映射,以确定空白校正图像中的像素点在鱼眼图像中对应的亚像素坐标,对确定的亚像素坐标进行线性插值,以得到亚像素值,再将亚像素值填充到空白校正图像内,遍历空白校正图像内的每个像素点,直至空白校正图像内每个像素点均赋予对应的亚像素值,得到鱼眼图像的校正图像。
上述方法虽然能对鱼眼图像进行校正,但是该方法首先要将鱼眼图像映射成球面图像,而将鱼眼图像映射成球面图像计算量非常大,并且通过该方法对鱼眼图像进行校正时,并不需要用到全部的鱼眼图像的像素,因此该方法的效率较低。
发明内容
本发明的目的是提供一种基于坐标转换的鱼眼图像校正方法及装置,通过改进坐标转换的球面映射方法,能以鱼眼图像中的任意一个点为中心对鱼眼图像进行快速、精准的校正。
为了实现上述目的,本发明技术方案如下:
一种基于坐标转换的鱼眼图像校正方法,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,所述基于坐标转换的鱼眼图像校正方法,包括:
从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
进一步地,所述根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标,包括:
以所述鱼眼半球中心点O为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面,建立的三维坐标系XYZ,计算P′对应在等距投影模型中鱼眼半球上的入射点P的球面坐标P的公式如下:
f=r*2/π
θ=π-arctan{(d-m),(e-n)}
其中,f为鱼眼摄像机的等效焦距,为OP与Z轴的夹角,P0为P点在XOY平面上的投影,θ为OP0与X轴的夹角,r为鱼眼图像的半径,m为P′点在鱼眼图像中的横坐标,n为P′点在鱼眼图像中的纵坐标,d为鱼眼图像中心点O"在鱼眼图像中的横坐标,e为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
进一步地,所述根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,计算出绘制点在鱼眼图像中对应的像素点的坐标,包括:
将XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使P点落在Z′轴上,得到坐标系X′Y′Z′;
根据校正图像中绘制点M的坐标(x,y)、Q点到O点的距离h,通过如下公式计算出校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球的交点M′的球面坐标(f,θ′)中的θ′:
其中,为OM′与Z轴的夹角,θ′为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为从坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
根据M′的球面坐标(f,θ′)中的θ′以及等距投影模型,通过如下公式计算出绘制点M在鱼眼图像上对应的像素点M"的坐标:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标,r为鱼眼图像的半径。
进一步地,所述基于坐标转换的鱼眼图像校正方法,还包括:
在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
本发明还提出了一种基于坐标转换的鱼眼图像校正装置,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,其特征在于,所述基于坐标转换的鱼眼图像校正装置,包括:
球面坐标计算模块,用于从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
虚拟成像面建立模块,用于将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
像素点坐标计算模块,用于根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
校正图像绘制模块,用于将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
进一步地,所述基于坐标转换的鱼眼图像校正装置,还包括:
校正区域调整模块,用于在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
本发明还提出了一种基于坐标转换的鱼眼图像校正装置,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,所述基于坐标转换的鱼眼图像校正装置,包括:
球面坐标计算模块,用于从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
虚拟成像面建立模块,用于将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
像素点坐标计算模块,用于根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
校正图像绘制模块,用于将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
进一步地,所述球面坐标计算模块根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标,包括:
以所述鱼眼半球中心点O为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面,建立的三维坐标系XYZ,计算P′对应在等距投影模型中鱼眼半球上的入射点P的球面坐标P的公式如下:
f=r*2/π
θ=π-arctan{(d-m),(e-n)}
其中,f为鱼眼摄像机的等效焦距,为OP与Z轴的夹角,P0为P点在在XOY平面上的投影,θ为OP0与X轴的夹角,r为鱼眼图像的半径,m为P′点在鱼眼图像中的横坐标,n为P′点在鱼眼图像中的纵坐标,d为鱼眼图像中心点O"在鱼眼图像中的横坐标,e为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
进一步地,所述像素点坐标计算模块根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,计算出绘制点在鱼眼图像中对应的像素点的坐标,执行如下操作:
将XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使P点落在Z′轴上,得到坐标系X′Y′Z′;
根据校正图像中绘制点M的坐标(x,y)、Q点到O点的距离h,通过如下公式计算出校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球的交点M′的球面坐标(f,θ′)中的θ′:
其中,为OM′与Z轴的夹角,θ′为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为从坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
根据M′的球面坐标(f,θ′)中的θ′以及等距投影模型,通过如下公式计算出绘制点M在鱼眼图像上对应的像素点M"的坐标:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标,r为鱼眼图像的半径。
进一步地,所述基于坐标转换的鱼眼图像校正装置,还包括:
校正区域调整模块,用于在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
本发明提出了一种基于坐标转换的鱼眼图像校正方法及装置,通过等距投影模型和坐标转换对鱼眼图像进行快速、精准的校正,本发明只需获取鱼眼图像的图像半径和中心点的位置,而无需获取鱼眼摄像机的焦距等参数,也无需对鱼眼图像的畸变参数进行标定就可以对鱼眼图像内的任意区域进行校正,得到满足透视原理的校正效果,并且与相似的现有技术比较,具有清晰易懂,计算过程简单,计算量小,并行性高的优点。该校正方法可以对高清鱼眼摄相机进行多方位实时校正,便于观察鱼眼图像中各个不同区域的细节。
附图说明
图1为本发明基于坐标转换的鱼眼图像校正方法的流程图;
图2为本实施例鱼眼摄像机等距投影模型的示意图;
图3为本实施例基于坐标转换的球面投射校正方法的示意图;
图4为本实施例控制校正图像在鱼眼图像中对应的校正区域大小的示意图;
图5为本发明基于坐标转换的鱼眼图像校正装置的结构示意图。
具体实施方式
下面结合附图和实施例对本发明技术方案做进一步详细说明,以下实施例不构成对本发明的限定。
本实施例一种基于坐标转换的鱼眼图像校正方法,如图1所示,包括:
步骤S1、从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标。
鱼眼摄像机的成像模型有等距投影、等立体角投影、体视投影和正交投影等,其中最常用的是等距投影。本实施例鱼眼摄像机等距投影模型的示意图,如图2所示,其中,O点为鱼眼摄像机的光学中心,XYZ为以O点为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面建立的三维坐标系,T为Z轴与鱼眼半球的交点,OT=f,f为鱼眼摄像机的等效焦距,UVO”为鱼眼摄像机的成像面,O"为鱼眼摄像机所成圆形图像(即鱼眼图像)的中心点,O"U与XYZ坐标系中的X轴以及鱼眼图像的横坐标轴i平行,O"V与XYZ坐标系中的Y轴以及鱼眼图像的纵坐标轴j平行。
假设空间中任意一点在鱼眼摄像机中成像时,该点与鱼眼摄像机光学中心O的连线交鱼眼半球于P点,P′点为P点在鱼眼摄像机中的成像点,则根据等距投影的原理得到:
∠XOP0=π-∠UO”P'=θ
其中,P0为P点在在XOY平面上的投影,为OP与Z轴的夹角,θ为OP0与X轴的夹角,为弧长。
从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据等距投影模型及光路可逆的原理可以得到P′在鱼眼半球上的入射点P,然后通过等距投影的逆向过程计算出点P′在鱼眼半球的入射点P的球面坐标P
θ=π-∠UO”P'
P′在鱼眼图像中的坐标为(m,n),鱼眼图像的中心点O"在鱼眼图像中的坐标为(d,e),鱼眼图像的半径为r(即r个像素点间的距离),则得到P′在坐标系UVO”中的坐标为:
u=m-d
v=n-e
其中u为P′在坐标系UVO”中的横坐标,v为P′在坐标系UVO”中的纵坐标。
根据勾股定理及三角函数得到:
∠UO”P'=arctan(u,v)
根据等距投影模型得到,当P′位于鱼眼图像的边缘位置,即P′到O"的距离等于鱼眼图像的半径r时,P位于鱼眼半球圆形底面的边缘位置,即等于π/2,则根据公式得到:
f=r*2/π
因此得到P的球面坐标P为:
f=r*2/π
θ=π-arctan{(m-d),(n-e)}
通过上述方法,本实施例根据鱼眼图像中任意一点P′的坐标(m,n)、鱼眼图像中心点O"的坐标(d,e)、以及鱼眼图像的半径r,通过鱼眼摄像机的等距投影模型计算出P′在鱼眼半球上的入射点P的球面坐标P
步骤S2、将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面。
本实施例计算出鱼眼图像中P′在鱼眼半球上的入射点P后,以P为切点,作一个与鱼眼半球相切的切面,将与该切面平行的任意平面作为鱼眼图像的校正图像所在的平面,校正图像所在的平面与射线OP的交点为Q点,然后以Q点为校正图像的中心点建立校正图像的虚拟成像面。
本实施例以校正图像所在的平面为鱼眼半球的切面为例进行说明,此时Q点即为P点,Q点到O点的距离h等于f。
步骤S3、根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标。
本实施例以Q点为校正图像的中心点建立校正图像的虚拟成像面后,如图3所示,将鱼眼摄像机等距投影模型中的XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使Q点落在Z′轴上,得到坐标系X′Y′Z′,Q点在坐标系X′Y′Z′中的三维坐标为(0,0,h)。
校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球相交于M′点,则根据透视投影的原理,M′为该绘制点M在鱼眼摄像机中成像时的入射点,然后根据鱼眼摄像机等距投影模型可以得到M′在鱼眼图像中对应的像素点M",M"即为绘制点M在鱼眼图像中对应的像素点。
首先,本实施例根据绘制点M在坐标系X′Y′Z′中的坐标计算出M′在XYZ坐标系中的坐标。具体地,设绘制点M在校正图像中的坐标为(x,y),则绘制点M在坐标系X′Y′Z′中的坐标为(x,y,h),根据三角函数得到M′在坐标系XYZ中的坐标为:然后根据坐标系X′Y′Z′与坐标系XYZ的旋转关系得到:
其中,为OM'与Z轴的夹角,θ'为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
则通过转换得到:
本实施例根据上述公式可以计算出M'的球面坐标中和θ′的值,因此得到M′的球面坐标(f,θ′),f为鱼眼摄像机的等效焦距。
然后根据等距投影模型得到M'在鱼眼图像中对应的像素点M"在坐标系UVO”中的坐标为:
其中,s为M"在坐标系UVO”中的横坐标,t为M"在坐标系UVO”中的纵坐标。
则得到M"在鱼眼图像中的坐标为:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
通过上述方法,本实施例根据校正图像中绘制点的坐标、鱼眼图像的半径,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标。
容易理解的是,本实施例对于校正图像中绘制点的坐标,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标,采用了几何学方法中一种计算公式。本领域技术人员还可以采用鱼眼摄像机的其他成像模型来进行计算,如采用等立体角投影、体视投影和正交投影等,其计算的公式与等距投影类似,这里不再赘述。
步骤S4、将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
本实施例计算出校正图像中的绘制点在鱼眼图像中对应的像素点的横坐标和纵坐标后,得到该绘制点在鱼眼图像中对应的像素点,将该绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值,在校正图像中绘制出该绘制点。
本实施例通过上述方法得到校正图像中的每个绘制点在鱼眼图像中对应的像素点,并将各个绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值,在校正图像中绘制出所有的绘制点,生成鱼眼图像对应的校正图像。
需要说明的是,本实施例在对鱼眼图像进行校正时,可以通过鼠标操作来实时选取对鱼眼图像进行校正的中心位置。优选地,本实施例通过鼠标在鱼眼图像中进行点击来实时选取对鱼眼图像进行校正的中心位置,点击的方式可以为单击或双击鼠标左键,也可以是单击或双击鼠标右键,本实施例不限定鼠标点击的方式。本实施例也可以通过按住鼠标左键或右键对校正图像进行拖动来调整对鱼眼图像进行校正的中心位置。本实施例还可以在鱼眼图像中选取多个不同的中心点对鱼眼图像的不同区域同时进行校正,从而实现对鱼眼图像的多个不同区域同时进行观察,使监控画面无死角。
本实施例在对鱼眼图像进行实时校正时,通过调整校正图像的虚拟成像面的中心点Q到O点的距离h,可以控制校正图像在鱼眼图像中对应的校正区域的大小,当增大Q点到O点的距离h时,校正图像在鱼眼图像中对应的校正区域缩小,使校正图像中物体所成的图像放大,图像细节更加清楚;当缩短Q点到O点的距离h时,校正图像的在鱼眼图像中对应的校正区域扩大,使校正图像的监控区域扩大。优选地,本实施例通过鼠标滚轮的滚动来调整Q点到O点的距离,鼠标滚轮向下滚动时Q点到O点的距离变小,鼠标滚轮向上滚动时Q点到O点的距离变大。本实施例控制校正图像在鱼眼图像中对应的校正区域大小的示意图,如图4所示,O点为鱼眼摄像机的光学中心,平面0为校正图像的虚拟成像面与鱼眼半球相切时的位置,OA为平面0的法线方向,平面1和平面2与平面0平行。当鼠标滚轮向下滚动时,校正图像的虚拟成像面沿着法线方向OA,从平面0向平面1移动,逐渐扩大校正图像的在鱼眼图像中对应的校正区域;当鼠标滚轮向上滚动时,校正图像的虚拟成像面沿着法线方向OA,从平面0向平面2移动,逐渐缩小校正图像的在鱼眼图像中对应的校正区域。
本实施例还提出了一种基于坐标转换的鱼眼图像校正装置,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,与上述方法对应,如图5所示,包括:
球面坐标计算模块,用于从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
虚拟成像面建立模块,用于将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
像素点坐标计算模块,用于根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
校正图像绘制模块,用于将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
本实施例球面坐标计算模块球面坐标计算模块根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标,包括:
以所述鱼眼半球中心点O为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面,建立的三维坐标系XYZ,计算P′对应在等距投影模型中鱼眼半球上的入射点P的球面坐标的公式如下:
f=r*2/π
θ=π-arctan{(d-m),(e-n)}
其中,f为鱼眼摄像机的等效焦距,为OP与Z轴的夹角,P0为P点在在XOY平面上的投影,θ为OP0与X轴的夹角,r为鱼眼图像的半径,m为P′点在鱼眼图像中的横坐标,n为P′点在鱼眼图像中的纵坐标,d为鱼眼图像中心点O"在鱼眼图像中的横坐标,e为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
本实施例像素点坐标计算模块根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,计算出绘制点在鱼眼图像中对应的像素点的坐标,执行如下操作:
将XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使P点落在Z′轴上,得到坐标系X′Y′Z′;
根据校正图像中绘制点M的坐标(x,y)、Q点到O点的距离h,通过如下公式计算出校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球的交点M′的球面坐标(f,θ′)中的θ′:
其中,为OM′与Z轴的夹角,θ′为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为从坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
根据M′的球面坐标(f,θ′)中的θ′以及等距投影模型,通过如下公式计算出绘制点M在鱼眼图像上对应的像素点M"的坐标:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标,r为鱼眼图像的半径。
本实施例基于坐标转换的鱼眼图像校正装置,还包括:
校正区域调整模块,用于在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
以上实施例仅用以说明本发明的技术方案而非对其进行限制,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (8)
1.一种基于坐标转换的鱼眼图像校正方法,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,其特征在于,所述基于坐标转换的鱼眼图像校正方法,包括:
从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
2.根据权利要求1所述的基于坐标转换的鱼眼图像校正方法,其特征在于,所述成像模型为等距投影,所述根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标,包括:
以所述鱼眼半球中心点O为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面,建立的三维坐标系XYZ,计算P′对应在等距投影模型中鱼眼半球上的入射点P的球面坐标的公式如下:
f=r*2/π
θ=π-arctan{(d-m),(e-n)}
其中,f为鱼眼摄像机的等效焦距,为OP与Z轴的夹角,P0为P点在XOY平面上的投影,θ为OP0与X轴的夹角,r为鱼眼图像的半径,m为P′点在鱼眼图像中的横坐标,n为P′点在鱼眼图像中的纵坐标,d为鱼眼图像中心点O"在鱼眼图像中的横坐标,e为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
3.根据权利要求2所述的基于坐标转换的鱼眼图像校正方法,其特征在于,所述根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,计算出绘制点在鱼眼图像中对应的像素点的坐标,包括:
将XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使P点落在Z′轴上,得到坐标系X′Y′Z′;
根据校正图像中绘制点M的坐标(x,y)、Q点到O点的距离h,通过如下公式计算出校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球的交点M′的球面坐标中的θ′:
其中,为OM′与Z轴的夹角,θ′为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为从坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
根据M′的球面坐标中的θ′以及等距投影模型,通过如下公式计算出绘制点M在鱼眼图像上对应的像素点M"的坐标:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标,r为鱼眼图像的半径。
4.根据权利要求1所述的基于坐标转换的鱼眼图像校正方法,其特征在于,所述基于坐标转换的鱼眼图像校正方法,还包括:
在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
5.一种基于坐标转换的鱼眼图像校正装置,用于对鱼眼摄像机拍摄的鱼眼图像进行校正,其特征在于,所述基于坐标转换的鱼眼图像校正装置,包括:
球面坐标计算模块,用于从鱼眼图像中任意选取一个像素点P′作为对鱼眼图像进行校正的中心点,根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标;
虚拟成像面建立模块,用于将平行于以P点为切点与所述鱼眼半球相切的切面的任意平面作为校正图像所在的平面,以该平面与所述鱼眼半球中心点O到P点的射线OP的交点Q为校正图像的中心点建立校正图像的虚拟成像面;
像素点坐标计算模块,用于根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,逐个计算出各绘制点在鱼眼图像中对应的像素点的坐标;
校正图像绘制模块,用于将计算出的绘制点在鱼眼图像中对应的像素点的像素值作为该绘制点在校正图像中的像素值在校正图像中进行绘制,得到鱼眼图像的校正图像。
6.根据权利要求5所述的基于坐标转换的鱼眼图像校正装置,其特征在于,所述球面坐标计算模块根据P′点在鱼眼图像中的坐标、鱼眼图像的中心点的坐标以及鱼眼图像的半径,通过鱼眼摄像机的成像模型计算出P′对应在成像模型中鱼眼半球上的入射点P的球面坐标,包括:
以所述鱼眼半球中心点O为坐标原点,以鱼眼半球的圆形底面为XOY所在的平面,建立的三维坐标系XYZ,计算P′对应在等距投影模型中鱼眼半球上的入射点P的球面坐标的公式如下:
f=r*2/π
θ=π-arctan{(d-m),(e-n)}
其中,f为鱼眼摄像机的等效焦距,为OP与Z轴的夹角,P0为P点在在XOY平面上的投影,θ为OP0与X轴的夹角,r为鱼眼图像的半径,m为P′点在鱼眼图像中的横坐标,n为P′点在鱼眼图像中的纵坐标,d为鱼眼图像中心点O"在鱼眼图像中的横坐标,e为鱼眼图像中心点O"在鱼眼图像中的纵坐标。
7.根据权利要求6所述的基于坐标转换的鱼眼图像校正装置,其特征在于,所述像素点坐标计算模块根据校正图像中绘制点的坐标、鱼眼图像中心点的坐标、鱼眼图像的半径、以及Q点到O点的距离,计算出绘制点在鱼眼图像中对应的像素点的坐标,执行如下操作:
将XYZ坐标系绕Z轴旋转θ,使X轴旋转到与OP0重合的位置,然后再绕Y轴旋转使P点落在Z′轴上,得到坐标系X′Y′Z′;
根据校正图像中绘制点M的坐标(x,y)、Q点到O点的距离h,通过如下公式计算出校正图像中的绘制点M与鱼眼光学中心O的连线与鱼眼半球的交点M′的球面坐标中的θ′:
其中,为OM′与Z轴的夹角,θ′为M′在XOY上的投影与X轴的夹角,OM的长度为L,R为从坐标系XYZ到坐标系X′Y′Z′的旋转矩阵,
根据M′的球面坐标中的θ′以及等距投影模型,通过如下公式计算出绘制点M在鱼眼图像上对应的像素点M"的坐标:
其中,i为M"在鱼眼图像中的横坐标,j为M"在鱼眼图像中的纵坐标,m为鱼眼图像中心点O"在鱼眼图像中的横坐标,n为鱼眼图像中心点O"在鱼眼图像中的纵坐标,r为鱼眼图像的半径。
8.根据权利要求5所述的基于坐标转换的鱼眼图像校正装置,其特征在于,所述基于坐标转换的鱼眼图像校正装置,还包括:
校正区域调整模块,用于在对鱼眼图像进行实时校正时,如果需要扩大校正图像在鱼眼图像中对应的校正区域,则缩短校正图像的中心点Q到O点的距离h,如果需要缩小校正图像的在鱼眼图像中对应的校正区域,则增大校正图像的中心点Q到O点的距离h。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611206384.8A CN106780389B (zh) | 2016-12-23 | 2016-12-23 | 一种基于坐标转换的鱼眼图像校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611206384.8A CN106780389B (zh) | 2016-12-23 | 2016-12-23 | 一种基于坐标转换的鱼眼图像校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106780389A true CN106780389A (zh) | 2017-05-31 |
CN106780389B CN106780389B (zh) | 2021-06-15 |
Family
ID=58919756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611206384.8A Active CN106780389B (zh) | 2016-12-23 | 2016-12-23 | 一种基于坐标转换的鱼眼图像校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106780389B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107454326A (zh) * | 2017-08-23 | 2017-12-08 | 广州市九安光电技术股份有限公司 | 利用鱼眼镜头进行全景摄像的方法、摄像机以及全景摄像*** |
CN107665483A (zh) * | 2017-09-27 | 2018-02-06 | 天津智慧视通科技有限公司 | 免定标便捷的单目镜头鱼眼图像畸变矫正方法 |
CN107808367A (zh) * | 2017-11-29 | 2018-03-16 | 天津聚飞创新科技有限公司 | 鱼眼图像校正方法、装置、无人机及存储介质 |
CN109308686A (zh) * | 2018-08-16 | 2019-02-05 | 北京市商汤科技开发有限公司 | 一种鱼眼图像处理方法及装置、设备和存储介质 |
CN111047638A (zh) * | 2018-10-12 | 2020-04-21 | 富华科精密工业(深圳)有限公司 | 单目测距装置 |
CN111353945A (zh) * | 2018-12-21 | 2020-06-30 | 杭州海康威视数字技术股份有限公司 | 鱼眼图像校正方法、装置及存储介质 |
CN111426333A (zh) * | 2020-02-25 | 2020-07-17 | 上海航天控制技术研究所 | 基于几何方法的火星导航敏感器图像形心精确修正方法 |
CN111861904A (zh) * | 2020-06-16 | 2020-10-30 | 浙江大华技术股份有限公司 | 等比例鱼眼矫正方法、装置和计算机设备和可读存储介质 |
CN115049535A (zh) * | 2022-05-24 | 2022-09-13 | 张全胜 | 鱼眼镜头有效区域获取与图像精细校正方法 |
WO2024098405A1 (zh) * | 2022-11-11 | 2024-05-16 | 网易(杭州)网络有限公司 | 一种相机标定方法、立体重建方法及图像校正方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101814181A (zh) * | 2010-03-17 | 2010-08-25 | 天津理工大学 | 一种鱼眼图像复原的展开方法 |
CN104778656A (zh) * | 2015-04-22 | 2015-07-15 | 江苏物联网研究发展中心 | 基于球面透视投影的鱼眼图像校正方法 |
CN105354808A (zh) * | 2015-12-02 | 2016-02-24 | 深圳华强数码电影有限公司 | 鱼眼图像的矫正方法 |
WO2016047893A1 (ko) * | 2014-09-24 | 2016-03-31 | 중앙대학교 산학협력단 | 어안렌즈 영상 보정 방법 및 그 장치 |
CN105488766A (zh) * | 2015-11-23 | 2016-04-13 | 上海汽车集团股份有限公司 | 鱼眼镜头图像校正方法及装置 |
CN105488764A (zh) * | 2015-11-23 | 2016-04-13 | 上海汽车集团股份有限公司 | 鱼眼图像校正方法及装置 |
CN106204496A (zh) * | 2016-07-18 | 2016-12-07 | 北京奇虎科技有限公司 | 一种云台摄像头的鱼眼图像处理方法及装置 |
-
2016
- 2016-12-23 CN CN201611206384.8A patent/CN106780389B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101814181A (zh) * | 2010-03-17 | 2010-08-25 | 天津理工大学 | 一种鱼眼图像复原的展开方法 |
WO2016047893A1 (ko) * | 2014-09-24 | 2016-03-31 | 중앙대학교 산학협력단 | 어안렌즈 영상 보정 방법 및 그 장치 |
CN104778656A (zh) * | 2015-04-22 | 2015-07-15 | 江苏物联网研究发展中心 | 基于球面透视投影的鱼眼图像校正方法 |
CN105488766A (zh) * | 2015-11-23 | 2016-04-13 | 上海汽车集团股份有限公司 | 鱼眼镜头图像校正方法及装置 |
CN105488764A (zh) * | 2015-11-23 | 2016-04-13 | 上海汽车集团股份有限公司 | 鱼眼图像校正方法及装置 |
CN105354808A (zh) * | 2015-12-02 | 2016-02-24 | 深圳华强数码电影有限公司 | 鱼眼图像的矫正方法 |
CN106204496A (zh) * | 2016-07-18 | 2016-12-07 | 北京奇虎科技有限公司 | 一种云台摄像头的鱼眼图像处理方法及装置 |
Non-Patent Citations (2)
Title |
---|
SIXIAN CHAN等: "An Improved Method for Fisheye Camera Calibration and Distortion Correction", 《2016 INTERNATIONAL CONFERENCE ON ADVANCED ROBOTICS AND MECHATRONICS (ICARM)》 * |
孟灵: "等距投影的鱼眼图像畸变校正算法设计", 《数字技术与应用》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107454326B (zh) * | 2017-08-23 | 2021-01-26 | 广州市九安物联科技有限公司 | 利用鱼眼镜头进行全景摄像的方法、摄像机以及全景摄像*** |
CN107454326A (zh) * | 2017-08-23 | 2017-12-08 | 广州市九安光电技术股份有限公司 | 利用鱼眼镜头进行全景摄像的方法、摄像机以及全景摄像*** |
CN107665483B (zh) * | 2017-09-27 | 2020-05-05 | 天津智慧视通科技有限公司 | 免定标便捷的单目镜头鱼眼图像畸变矫正方法 |
CN107665483A (zh) * | 2017-09-27 | 2018-02-06 | 天津智慧视通科技有限公司 | 免定标便捷的单目镜头鱼眼图像畸变矫正方法 |
CN107808367A (zh) * | 2017-11-29 | 2018-03-16 | 天津聚飞创新科技有限公司 | 鱼眼图像校正方法、装置、无人机及存储介质 |
CN109308686A (zh) * | 2018-08-16 | 2019-02-05 | 北京市商汤科技开发有限公司 | 一种鱼眼图像处理方法及装置、设备和存储介质 |
CN109308686B (zh) * | 2018-08-16 | 2022-06-24 | 北京市商汤科技开发有限公司 | 一种鱼眼图像处理方法及装置、设备和存储介质 |
CN111047638A (zh) * | 2018-10-12 | 2020-04-21 | 富华科精密工业(深圳)有限公司 | 单目测距装置 |
CN111047638B (zh) * | 2018-10-12 | 2023-08-18 | 深圳富联富桂精密工业有限公司 | 单目测距装置 |
CN111353945A (zh) * | 2018-12-21 | 2020-06-30 | 杭州海康威视数字技术股份有限公司 | 鱼眼图像校正方法、装置及存储介质 |
CN111353945B (zh) * | 2018-12-21 | 2023-10-20 | 杭州海康威视数字技术股份有限公司 | 鱼眼图像校正方法、装置及存储介质 |
CN111426333A (zh) * | 2020-02-25 | 2020-07-17 | 上海航天控制技术研究所 | 基于几何方法的火星导航敏感器图像形心精确修正方法 |
CN111426333B (zh) * | 2020-02-25 | 2022-03-04 | 上海航天控制技术研究所 | 基于几何方法的火星导航敏感器图像形心精确修正方法 |
CN111861904A (zh) * | 2020-06-16 | 2020-10-30 | 浙江大华技术股份有限公司 | 等比例鱼眼矫正方法、装置和计算机设备和可读存储介质 |
CN115049535A (zh) * | 2022-05-24 | 2022-09-13 | 张全胜 | 鱼眼镜头有效区域获取与图像精细校正方法 |
CN115049535B (zh) * | 2022-05-24 | 2024-01-12 | 杭州追形视频科技有限公司 | 鱼眼镜头有效区域获取与图像精细校正方法 |
WO2024098405A1 (zh) * | 2022-11-11 | 2024-05-16 | 网易(杭州)网络有限公司 | 一种相机标定方法、立体重建方法及图像校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106780389B (zh) | 2021-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106780389A (zh) | 一种基于坐标转换的鱼眼图像校正方法及装置 | |
US10733783B2 (en) | Motion smoothing for re-projected frames | |
CN108122191B (zh) | 鱼眼图像拼接成全景图像和全景视频的方法及装置 | |
CN101783011B (zh) | 一种鱼眼镜头的畸变校正方法 | |
CN110197466B (zh) | 一种广角鱼眼图像矫正方法 | |
CN104574350B (zh) | 三维数据撷取方法及其*** | |
CN102984453B (zh) | 利用单摄像机实时生成半球全景视频图像的方法及*** | |
CN104778656B (zh) | 基于球面透视投影的鱼眼图像校正方法 | |
CN106570938A (zh) | 基于opengl 的全景监控方法及*** | |
CN104537616A (zh) | 鱼眼图像畸变的校正方法 | |
JP2018205074A (ja) | 飛行パラメータ測定装置及び飛行パラメータ測定方法 | |
CN101132535A (zh) | 基于转台的多投影大屏拼接方法 | |
CN107527327B (zh) | 一种基于gpu的鱼眼校正方法 | |
CN106920210B (zh) | 一种基于球面展开模型的鱼眼视频实时全景校正方法 | |
CN107145224B (zh) | 基于三维球面泰勒展开的人眼视线跟踪方法和装置 | |
CN108717704B (zh) | 基于鱼眼图像的目标跟踪方法、计算机装置及计算机可读存储介质 | |
CN109696954A (zh) | 视线追踪方法、装置、设备和存储介质 | |
CN108154536A (zh) | 二维平面迭代的相机标定法 | |
CN106600549A (zh) | 鱼眼图像的矫正方法和装置 | |
CN108230242A (zh) | 一种从全景激光点云到视频流的转换方法 | |
CN104807405A (zh) | 一种基于光线角度标定的三维坐标测量方法 | |
WO2017206824A1 (zh) | 一种智能显示视觉图像的位置信息处理方法及*** | |
CN109376595B (zh) | 基于人眼注意力的单目rgb摄像头活体检测方法及*** | |
CN104751455A (zh) | 一种作物图像稠密匹配方法及*** | |
CN104240294A (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 |