CN107121096A - 一种基于图像的空间立体角计算方法 - Google Patents
一种基于图像的空间立体角计算方法 Download PDFInfo
- Publication number
- CN107121096A CN107121096A CN201710203506.6A CN201710203506A CN107121096A CN 107121096 A CN107121096 A CN 107121096A CN 201710203506 A CN201710203506 A CN 201710203506A CN 107121096 A CN107121096 A CN 107121096A
- Authority
- CN
- China
- Prior art keywords
- image
- angle
- solid angle
- pixel
- rectangle
- 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
- 238000000205 computational method Methods 0.000 title abstract 2
- 239000007787 solid Substances 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000004576 sand Substances 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 abstract 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明属于光电检测技术领域,具体为一种基于图像的空间立体角计算方法。本发明方法包括,在获取图像的像素与视场角之后,计算出视场光阑与图像之间在像素坐标系中的垂直距离,并计算图像中任意点与视场光阑中心点连线相对于图像法线方向的角度,将图像中每一个像素看作一个微小的矩形,计算矩形的四个角在像素坐标系中的坐标以及与视场光阑的连线投影到视场光阑所在水平面和竖直面后与图像法线之间的夹角,将某一区域中所有像素对应的立体角进行叠加,即计算图像中任意区域在空间中的立体角。本发明方法误差小、计算复杂度低。
Description
技术领域
本发明属于光电检测技术领域,具体涉及一种在在给定图像像素值及视场角的情况下计算图像中任意区域在空间中的立体角的计算方法。
技术背景
在照明领域,许多光源参数的计算都涉及光源在空间中的立体角,而立体角涉及到三维的空间无法被直接测量,在实际测量过程中,需要首先测量出光源的面积和光源与观察点之间的距离,再通过面积除以距离的平方来近似得到光源在空间中的立体角。
第二种计算立体角的方法为利用图像,首先计算光源中心点与视场光阑连线相对于图像法线的角度,再计算光源投影到连线法平面上的面积,再计算像素坐标系下光源中心点相对于视场光阑的距离,最后通过面积除以距离的平方近似得到光源在空间中的立体角。
上述两种测量方法都存在一定的缺陷,一方面,两种算法中,都是将平面面积近似看作球面面积来计算立体角,当光源的尺寸较大时,会存在比较大的误差;另一方面,在第一种方法中,对于形状不规则的光源,光源面积的测定相对复杂,在实际操作过程中相对复杂。针对上述缺陷,本发明的目的在于提供用基于图像的空间立体角计算方法。
发明内容
本发明的目的在于提供一种计算误差较小、计算过程简单的基于图像的空间立体角计算方法。
本发明提出的基于图像的空间立体角计算方法,具体步骤如下:
第一步,根据图像的像素数与最大视场角计算视场光阑在像素坐标系下的与图像的垂直距离,设图像的像素数为m×n,最大水平视场角为θ,计算视场光阑到图像的垂直距离为:
。
第二步,以视场光阑为原点,过原点与图像平行的平面为xy平面,视场光阑到图像垂直方向为z轴方向建立坐标系,对于图像中任意一个像素点 (a,b)(以左下角的像素点为(0,0)),将其看作一个微小的矩形,矩形的四个顶点素坐标系中的坐标分别为:
。
第三步,计算四个顶点与视场光阑的连线在水平面和垂直面的投影与图像法线方向的夹角,记水平方向夹角为α,垂直方向夹角为β;不失一般性,以为例,水平方向夹角α和垂直方向夹角β分别为:
。
第四步,计算由视场光阑观察,微型矩形在空间中的立体角
设由第三步计算出的四个顶点在水平方向和垂直方向上的夹角分别为:
根据立体角的定义,微型矩形在空间中的立体角为其对应的球面面积与球面半径的平方的比值,不失一般性,设球面的半径r 为1,则球面的方程为:
球面面积的积分为:
易知,积分区域在xy平面上的投影是由四条椭圆曲线围成,方程分别为:
以四条边界进行积分,可得对应球面的面积为:
则该微型矩形对应的空间立体角为:。
第五步,计算图像中任意区域中包含的各个像素点对应的空间立体角,将所有像素点的空间立体角进行线性叠加获得该区域对应的立体角Ω:
。
在本计算方法中,曲面的计算是通过微积分的方式直接计算,不需要将曲面近似为平面计算,因此避免了在计算过程中由于平面近似曲面而带来的误差。另一方面,在本算法中,每一个像素对应的空间立体角可提前存储于设备中,在使用过程中只需要做查表和求和的运算,而不需要进行大量的平方和开方运算,减少了使用过程中的计算量。
附图说明
图1为计算流程图。
图2为空间示意图。
图3为局部放大图。
图4为几何参数示意图。
图5为几何参数示意图。
具体实施方式
下面通过实施例进一步说明本发明。
实施例:当图像的像素数为1280×720,最大水平视场角为120°时,可以求出视场光阑与图像的垂直距离L为:
所求区域包含了四个像素点围成的区域,以像素点 (20,30)为例,其近似的微型矩形的四个顶点在像素坐标系中的坐标分别为:
其与视场光阑的连线在水平方向和垂直方向的投影与图像法线方向的夹角分别为:
该微型矩形对应的球面面积:
相应的空间立体角:。
同理可计算区域中其他像素点对应的空间立体角的大小,将所有像素点对应的立体角相加即可得到该区域对应的立体角:。
若将上述区域对应的曲面面积近似为平面面积计算立体角,在相同的像素分辨率下,计算得到的立体角为:。本例中,在曲面计算过程中本算法由于去除近似过程避免了9.6%的误差,而对于更大面积的区域,本算法将减少更大的误差。
Claims (2)
1.一种基于图像的空间立体角计算方法,其特征在于,具体步骤为:
第一步,根据图像的像素数与最大视场角计算视场光阑在像素坐标系下的与图像的垂直距离,设图像的像素数为m×n,最大水平视场角为θ,计算视场光阑到图像的垂直距离为:
第二步,以视场光阑为原点,过原点与图像平行的平面为xy平面,视场光阑到图像垂直方向为z轴方向建立坐标系,对于图像中任意一个像素点(a,b),其中,以左下角的像素点为(0,0),将其看作一个微小的矩形,矩形的四个顶点像素坐标系中的坐标分别为:
第三步,计算矩形的四个顶点与视场光阑的连线在水平面和垂直面的投影与图像法线方向的夹角,记水平方向夹角为α,垂直方向夹角为β;
第四步,计算由视场光阑观察,微型矩形在空间中的立体角
根据立体角的定义,微型矩形在空间中的立体角为其对应的球面面积S与球面半径r的平方的比值,则该微型矩形对应的空间立体角为:Ω=S/r2;
第五步,计算图像中任意区域中包含的各个像素点对应的空间立体角,将所有像素点的空间立体角进行线性叠加获得该区域对应的立体角Ω:。
2.根据权利要求1所述的基于图像的空间立体角的计算方法,其特征在于,微型矩形在空间中的立体角为其对应的球面面积S为:
其中,分别为由第三步计算出的四个顶点在水平方向和垂直方向上的夹角,球面半径r为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710203506.6A CN107121096B (zh) | 2017-03-30 | 2017-03-30 | 一种基于图像的空间立体角计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710203506.6A CN107121096B (zh) | 2017-03-30 | 2017-03-30 | 一种基于图像的空间立体角计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107121096A true CN107121096A (zh) | 2017-09-01 |
CN107121096B CN107121096B (zh) | 2020-04-24 |
Family
ID=59717456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710203506.6A Active CN107121096B (zh) | 2017-03-30 | 2017-03-30 | 一种基于图像的空间立体角计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107121096B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917556A (zh) * | 2018-03-23 | 2018-11-30 | 东莞理工学院 | 一种基于铣床定位的物体弯曲角度的测量方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86101381A (zh) * | 1986-03-06 | 1987-09-16 | 河南省林业科学研究所 | 立体角规及其测树方法 |
JPH02278104A (ja) * | 1989-04-19 | 1990-11-14 | Fuji Electric Co Ltd | 文書画像の傾き角検出方法 |
JPH10123460A (ja) * | 1996-10-18 | 1998-05-15 | Mitsubishi Electric Corp | 立体表示装置およびその画素ユニット |
JP2004302552A (ja) * | 2003-03-28 | 2004-10-28 | Mitsubishi Space Software Kk | 高輝度反射物体画像処理システム及び船舶画像処理システム及び高輝度反射物体画像処理方法及びプログラムを記録したコンピュータ読み取り可能な記録媒体及びプログラム |
CN101614530A (zh) * | 2009-07-29 | 2009-12-30 | 中国人民解放军国防科学技术大学 | 一种新型光束平漂和角漂实时检测方法及装置 |
CN102128595A (zh) * | 2009-12-30 | 2011-07-20 | 云林科技大学 | 雕刻刀具尺寸参数的检测方法及其装置 |
CN102506769A (zh) * | 2011-11-15 | 2012-06-20 | 益海芯电子技术江苏有限公司 | 旋转角度检测方法 |
CN103139580A (zh) * | 2011-11-29 | 2013-06-05 | 长春理工大学 | 一种三维全景空间立体图像生成方法 |
CN103759728A (zh) * | 2014-01-26 | 2014-04-30 | 长春理工大学 | 空间立体角准均匀分布导航星表划分方法 |
CN103940389A (zh) * | 2014-03-04 | 2014-07-23 | 北京林业大学 | 一种利用电子立体角规代替每木检尺的方法 |
JP2015182241A (ja) * | 2014-03-20 | 2015-10-22 | 京セラ株式会社 | 液体吐出ヘッド、およびそれを用いた記録装置 |
-
2017
- 2017-03-30 CN CN201710203506.6A patent/CN107121096B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86101381A (zh) * | 1986-03-06 | 1987-09-16 | 河南省林业科学研究所 | 立体角规及其测树方法 |
JPH02278104A (ja) * | 1989-04-19 | 1990-11-14 | Fuji Electric Co Ltd | 文書画像の傾き角検出方法 |
JPH10123460A (ja) * | 1996-10-18 | 1998-05-15 | Mitsubishi Electric Corp | 立体表示装置およびその画素ユニット |
JP2004302552A (ja) * | 2003-03-28 | 2004-10-28 | Mitsubishi Space Software Kk | 高輝度反射物体画像処理システム及び船舶画像処理システム及び高輝度反射物体画像処理方法及びプログラムを記録したコンピュータ読み取り可能な記録媒体及びプログラム |
CN101614530A (zh) * | 2009-07-29 | 2009-12-30 | 中国人民解放军国防科学技术大学 | 一种新型光束平漂和角漂实时检测方法及装置 |
CN102128595A (zh) * | 2009-12-30 | 2011-07-20 | 云林科技大学 | 雕刻刀具尺寸参数的检测方法及其装置 |
CN102506769A (zh) * | 2011-11-15 | 2012-06-20 | 益海芯电子技术江苏有限公司 | 旋转角度检测方法 |
CN103139580A (zh) * | 2011-11-29 | 2013-06-05 | 长春理工大学 | 一种三维全景空间立体图像生成方法 |
CN103759728A (zh) * | 2014-01-26 | 2014-04-30 | 长春理工大学 | 空间立体角准均匀分布导航星表划分方法 |
CN103940389A (zh) * | 2014-03-04 | 2014-07-23 | 北京林业大学 | 一种利用电子立体角规代替每木检尺的方法 |
JP2015182241A (ja) * | 2014-03-20 | 2015-10-22 | 京セラ株式会社 | 液体吐出ヘッド、およびそれを用いた記録装置 |
Non-Patent Citations (1)
Title |
---|
陈文成等: "白光LED汽车前照灯的光学设计探讨", 《中国照明电器》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917556A (zh) * | 2018-03-23 | 2018-11-30 | 东莞理工学院 | 一种基于铣床定位的物体弯曲角度的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107121096B (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10643302B2 (en) | Method and apparatus for generating three-dimensional panoramic video | |
US7570280B2 (en) | Image providing method and device | |
JP5058831B2 (ja) | 視体積の外側を分散クリッピングするための方法 | |
WO2017113731A1 (zh) | 360度全景显示方法、显示模块及移动终端 | |
WO2013018173A1 (ja) | 画像処理装置、画像処理方法および画像処理プログラム | |
CN110830781B (zh) | 一种基于双目视觉的投影图像自动校正方法及系统 | |
CN103839227B (zh) | 鱼眼图像校正方法和装置 | |
CN106558017B (zh) | 球形显示图像处理方法及系统 | |
CN107633489A (zh) | 鱼眼镜头圆心重提取映射法畸变矫正方法 | |
CN109949367B (zh) | 一种基于圆形投影的可见光成像定位方法 | |
CN106600549A (zh) | 鱼眼图像的矫正方法和装置 | |
CN106709865A (zh) | 一种深度图像合成方法及装置 | |
KR102449438B1 (ko) | 입방체 복원 장치 및 방법 | |
US10075634B2 (en) | Method and system for generating a surround view | |
WO2018102990A1 (en) | System and method for rectifying a wide-angle image | |
KR20160076271A (ko) | 장애물 검출 장치 및 장애물 검출 방법 | |
CN112015170A (zh) | 运动物体检测及智能驾驶控制方法、装置、介质及设备 | |
US11430206B1 (en) | Branch-and-bound search algorithm for 4D camera localization | |
US20140028665A1 (en) | Method, apparatus and computer for generating reflection texture | |
WO2017024626A1 (zh) | 裸眼3d成像方法及系统 | |
US20180213215A1 (en) | Method and device for displaying a three-dimensional scene on display surface having an arbitrary non-planar shape | |
CN107121096B (zh) | 一种基于图像的空间立体角计算方法 | |
CN107442973A (zh) | 基于机器视觉的焊道定位方法及装置 | |
CN111091594A (zh) | 多点云平面融合方法及装置 | |
CN107346530B (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 |