CN115546311B - 一种基于场景信息的投影仪标定方法 - Google Patents

一种基于场景信息的投影仪标定方法 Download PDF

Info

Publication number
CN115546311B
CN115546311B CN202211188222.1A CN202211188222A CN115546311B CN 115546311 B CN115546311 B CN 115546311B CN 202211188222 A CN202211188222 A CN 202211188222A CN 115546311 B CN115546311 B CN 115546311B
Authority
CN
China
Prior art keywords
projector
coordinate system
point
image
focal length
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.)
Active
Application number
CN202211188222.1A
Other languages
English (en)
Other versions
CN115546311A (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.)
Communication University of China
Original Assignee
Communication University of China
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 Communication University of China filed Critical Communication University of China
Priority to CN202211188222.1A priority Critical patent/CN115546311B/zh
Publication of CN115546311A publication Critical patent/CN115546311A/zh
Application granted granted Critical
Publication of CN115546311B publication Critical patent/CN115546311B/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
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Transforming Electric Information Into Light Information (AREA)

Abstract

本发明公开一种基于场景信息的投影仪标定方法,包括:计算投影仪镜头的实际焦距,基于所述实际焦距,得到投影仪的内参矩阵;计算投影仪的外参旋转矩阵;基于外参旋转矩阵,得到投影仪的外参平移矢量,完成标定。本发明不需要借助照相机实现投影仪标定,也不需要借助标定物,方法简单有效,可用于大场景文旅演出、展览展示时投影仪标定,为远距离、非平面幕投影几何校正提供技术支撑。

Description

一种基于场景信息的投影仪标定方法
技术领域
本发明涉及异形面动态投影技术领域,特别涉及一种基于场景信息的投影仪标定方法。
背景技术
投影仪内参、外参标定是非平面幕投影时投影仪正确投影的基础。由于投影仪不具备主动获取环境信息的能力,因此常用的投影仪标定方法需要借助照投影仪来实现,也就是把投影仪和照投影仪连接成一个刚性的系统,先标定照投影仪,再借助照投影仪标定投影仪。采用该方法标定投影仪,一方面会把投影仪标定时的误差传递到投影仪标定,另外,投影仪标定通常需要借助标定物来实现。大尺度场景下,标定物的尺寸会变得很大,标定物制作和使用困难。
发明内容
为解决上述现有技术中所存在的问题,本发明提供一种基于场景信息的投影仪标定方法,不需要借助照相机实现投影仪标定,也不需要借助标定物,方法简单有效,可用于大场景文旅演出、展览展示时投影仪标定,为远距离、非平面幕投影几何校正提供技术支撑。
为了实现上述技术目的,本发明提供了一种基于场景信息的投影仪标定方法,包括:
计算投影仪镜头的实际焦距,基于所述实际焦距,得到投影仪的内参矩阵;
计算所述投影仪的外参旋转矩阵;
基于所述外参旋转矩阵,得到所述投影仪的外参平移矢量,完成标定。
可选地,所述实际焦距的获取过程为:
设置所述投影仪的光心和投影仪图像中心;
将世界坐标系中x方向的平行线投影到投影仪图像的平面后,相交于第一消隐点;将所述世界坐标系中y方向的平行线投影到所述投影仪图像的平面后,相交于第二消隐点;
连接所述第一消隐点和所述第二消隐点,得到消隐点直线;
所述投影仪图像中心投影到所述消隐点直线上,得到投影点;
基于所述光心、所述投影仪图像中心和所述投影点,得到所述实际焦距。
可选地,所述实际焦距的计算公式为:
其中:
式中,f为实际焦距;Op为投影仪的光心;Oi为投影仪图像中心;Vi为投影仪图像中心投影到消隐点直线上的投影点的坐标;V1为第一消隐点的坐标;V2为第二消隐点的坐标。
可选地,所述内参矩阵为:
式中,K为内参矩阵;fx=f/dx、fy=f/dy分别为投影仪在投影仪图像水平方向和垂直方向的有效焦距,f为投影仪镜头实际焦距,dx和dy为成像平面上一个像素在水平方向和垂直方向所对应的实际物理尺寸;s为倾斜因子,可作为常数0处理;u0和v0为投影仪镜头主光轴与成像平面交点的图像坐标。
可选地,所述外参旋转矩阵的获取过程为:
构建一个新坐标系,所述新坐标系与世界坐标系方向相同,所述新坐标系和投影仪坐标系之间的旋转矩阵与所述世界坐标系和投影仪坐标系之间的旋转矩阵相同。
可选地,所述外参旋转矩阵为:
式中,v1i、v1j为第一消隐点的坐标值;v2i、v2j为第二消隐点的坐标值;z′cx、z′cy和z′cz分别为新坐标系的坐标轴在x方向、y方向和z方向的旋转矢量。
可选地,所述外参平移矢量的获取过程为:
选取场景中的长度已知线段,将所述长度已知线段的一个端点作为世界坐标系原点;
已知所述外参旋转矩阵,将所述长度已知线段变换到投影仪坐标系,再变换到投影仪图像坐标系,最后映射到像素坐标系,得到第一像素点和第二像素点;
基于所述第一像素点、所述第二像素点、投影仪图像水平方向分辨率和投影仪图像垂直方向分辨率,得到所述投影仪图像坐标系下的点与所述像素坐标系下的点的变换关系;
基于所述变换关系,平移所述长度已知线段,得到平移后线段;
计算所述世界坐标系原点到所述投影仪坐标系的原点间距离;
基于所述原点间距离和所述平移后线段,得到所述外参平移矢量。
可选地,所述外参平移矢量为:
式中,D为原点间距离;Op为投影仪光心;P1′为平移后线段的一端点的坐标。
本发明具有如下技术效果:
本发明不需要借助照相机实现投影仪标定,也不需要借助标定物,方法简单有效,可用于大场景文旅演出、展览展示时投影仪标定,为远距离、非平面幕投影几何校正提供技术支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例基于场景信息的投影仪标定方法的流程框图;
图2为本发明实施例的投影仪成像位置关系示意图;
图3为本发明实施例场景中长度已知线段的投影图;
图4为本发明实施例的实验场地图;
图5为本发明实施例中投影仪在物理世界投射出三个方向的平行线的示意图;
图6为本发明实施例投影仪图像示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明公开一种基于场景信息的投影仪标定方法,包括:
投影仪投影是相机成像的逆过程,都服从小孔成像原理。
如图2所示,设投影仪光心为Op,投影仪图像中心为Oi。世界坐标系中x方向的平行线投影到投影仪图像平面后相交于第一消隐点V1=(v1i,v1j),世界坐标系中y方向的平行线投影到投影仪图像平面后相交于第二消隐点V2=(v2i,v2j)。Oi在直线V1V2上的投影为Vi
投影仪内参矩阵为:
其中,fx=f/dx、fy=f/dy分别为投影仪在投影仪图像水平方向和垂直方向的有效焦距(以像素为单位),f为投影仪镜头实际焦距,dx和dy为成像平面上一个像素在水平方向和垂直方向所对应的实际物理尺寸;s为倾斜因子,用于描述图像传感器水平方向和垂直方向的垂直程度,可作为常数0处理;u0和v0为投影仪镜头主光轴与成像平面交点的图像坐标。
投影仪像素水平方向物理尺寸和垂直方法物理尺寸通常是一样的,即dx=dy,所以fx=fy=f。光心通常被认为和投影仪图像中心重叠,即u0、v0分别为投影仪图像水平方向分辨率/2、投影仪图像垂直方向分辨率/2。因此,要求得投影仪内参矩阵K,只需要求出投影仪镜头实际焦距f即可。
投影仪镜头实际焦距f的计算过程包括:
光心Op和投影仪图像中心Oi在光轴上,则实际焦距f为:
其中:
投影仪的外参矩阵包括旋转矩阵R和平移矢量t,旋转矩阵R的获取过程包括:
被投射到物理世界的同一方向的平行线对应的投影仪图像中的线段,在投影仪图像平面相交于同一个消隐点。消隐点V1和V2在世界坐标系的两个正交轴的方向上,世界坐标系的原点为OW。构造一个新坐标系,该新坐标系原点和投影仪坐标系原点重合,该新坐标系原点为Op,新坐标系的坐标轴为矢量X′c,Y′c,Z′c
Z′c=X′c×Y′c
构造的新坐标系与世界坐标系具有相同的方向,因此,新坐标系和投影仪坐标系之间的旋转矩阵,与世界坐标系和投影仪坐标系之间的旋转矩阵相同,可得旋转矩阵R为:
其中,z′cx、z′cy和z′cz分别为新坐标系的坐标轴Z′c在x方向、y方向和z方向的旋转矢量;v1i、v1j分别为投影仪图像中第一消隐点的横坐标和纵坐标;v2i、v2j分别为投影仪图像中第二消隐点的横坐标和纵坐标。
平移矢量t的获取过程包括:
平移矢量t为沿着右手笛卡尔坐标系的x,y,z轴分别的平移距离。不失一般性,世界坐标系的原点可以选择场景中的任意一点。如图3所示,已知场景中一条线段P1P2,p1P2长度已知。不妨选择P1作为世界坐标系的原点,这个特殊的点称为世界坐标原点Wd,则:
p1=[0,0,0]T
其中,分别表示P2在世界坐标系下的坐标。
由于旋转矩阵R已知,我们可以将线段P1P2转换到投影仪坐标系下,利用如下公式将世界坐标转换为投影仪坐标:
其中,分别表示P1点和P2点在投影仪坐标系下的坐标。
线段P1P2由世界坐标系变换到投影仪坐标系,再到投影仪图像坐标系,最后映射到像素坐标系,得到第一像素点和第二像素点/>投影仪图像坐标系下的点与像素坐标系下的点的变换关系为:
因此,能够平移投影仪图像平面上的线段,平移后的线段为P′1P′2,其中:
如图3所示,可以得到两个相似三角形△PpP1P2和△OpP′1Q。根据相似三角形的性质,则:
因此,投影仪坐标系的原点到世界坐标系的原点(即特殊点Wd)的距离D为:
得到平移矢量t为:
投影仪不具备主动获取物理世界信息的能力,因此,本发明采用投影仪主动投射出三组线,这三组线投射到物理世界的平行线上,投射出的这三组线在物理世界符合曼哈顿世界假设,即每组线互相平行,3组线相互垂直,符合世界坐标x、y、z方向。具体实验时可以选物理世界墙角三边作为x、y、z方向。实验场地如图4所示,图中包含3个方向的平行线,3个方向相互垂直。图5中投影仪在物理世界即实验场地中投射出了3个方向的平行线。图6为投影出图5中的3个方向的平行线时对应的投影仪图像。
进一步地,我们记画线时对应现实世界中平行线的几条线段为同一组。对图6中的投影仪图像进行处理,延长图6中的线段的长度,得到同一组线段的交点,即得到3个方向的消隐点,选择其中的两个消隐点作为实验用消隐点V1和V2
在进行上述所述工作时,同时要测量并记录其中几条手动标记线段所对应的物理世界的线段的长度,并且要选择其中一条线段的起始端点(在准备阶段做特殊标记,也就是投影仪标定中提到的世界特殊点Wd)作为世界坐标原点。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (5)

1.一种基于场景信息的投影仪标定方法,其特征在于,包括以下步骤:
计算投影仪镜头的实际焦距,基于所述实际焦距,得到投影仪的内参矩阵;
计算所述投影仪的外参旋转矩阵;
基于所述外参旋转矩阵,得到所述投影仪的外参平移矢量,完成标定;
所述实际焦距的获取过程为:
设置所述投影仪的光心和投影仪图像中心;
将世界坐标系中x方向的平行线投影到投影仪图像的平面后,相交于第一消隐点;将所述世界坐标系中y方向的平行线投影到所述投影仪图像的平面后,相交于第二消隐点;
连接所述第一消隐点和所述第二消隐点,得到消隐点直线;
所述投影仪图像中心投影到所述消隐点直线上,得到投影点;
基于所述光心、所述投影仪图像中心和所述投影点,得到所述实际焦距;
所述外参旋转矩阵为:
式中,v1i、v1j为第一消隐点的坐标值;v2i、v2j为第二消隐点的坐标值;z'cx、z'cy和z'cz分别为新坐标系的坐标轴在x方向、y方向和z方向的旋转矢量;
所述外参平移矢量的获取过程为:
选取场景中的长度已知线段,将所述长度已知线段的一个端点作为世界坐标系原点;
已知所述外参旋转矩阵,将所述长度已知线段变换到投影仪坐标系,再变换到投影仪图像坐标系,最后映射到像素坐标系,得到第一像素点和第二像素点;
基于所述第一像素点、所述第二像素点、投影仪图像水平方向分辨率和投影仪图像垂直方向分辨率,得到所述投影仪图像坐标系下的点与所述像素坐标系下的点的变换关系;
基于所述变换关系,平移所述长度已知线段,得到平移后线段;
计算所述世界坐标系原点到所述投影仪坐标系的原点间距离;
基于所述原点间距离和所述平移后线段,得到所述外参平移矢量。
2.根据权利要求1所述的基于场景信息的投影仪标定方法,其特征在于,所述实际焦距的计算公式为:
其中:
式中,f为实际焦距;Op为投影仪的光心;Oi为投影仪图像中心;Vi为投影仪图像中心投影到消隐点直线上的投影点的坐标;V1为第一消隐点的坐标;V2为第二消隐点的坐标。
3.根据权利要求1所述的基于场景信息的投影仪标定方法,其特征在于,所述内参矩阵为:
式中,K为内参矩阵;fx=f/dx、fy=f/dy分别为投影仪在投影仪图像水平方向和垂直方向的有效焦距,f为投影仪镜头实际焦距,dx和dy为成像平面上一个像素在水平方向和垂直方向所对应的实际物理尺寸;s为倾斜因子,可作为常数0处理;u0和v0为投影仪镜头主光轴与成像平面交点的图像坐标。
4.根据权利要求1所述的基于场景信息的投影仪标定方法,其特征在于,所述外参旋转矩阵的获取过程为:
构建一个新坐标系,所述新坐标系与世界坐标系方向相同,所述新坐标系和投影仪坐标系之间的旋转矩阵与所述世界坐标系和投影仪坐标系之间的旋转矩阵相同。
5.根据权利要求1所述的基于场景信息的投影仪标定方法,所述外参平移矢量为:
式中,D为原点间距离;Op为投影仪光心;P1'为平移后线段的一端点的坐标。
CN202211188222.1A 2022-09-28 2022-09-28 一种基于场景信息的投影仪标定方法 Active CN115546311B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211188222.1A CN115546311B (zh) 2022-09-28 2022-09-28 一种基于场景信息的投影仪标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211188222.1A CN115546311B (zh) 2022-09-28 2022-09-28 一种基于场景信息的投影仪标定方法

Publications (2)

Publication Number Publication Date
CN115546311A CN115546311A (zh) 2022-12-30
CN115546311B true CN115546311B (zh) 2023-07-25

Family

ID=84730460

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211188222.1A Active CN115546311B (zh) 2022-09-28 2022-09-28 一种基于场景信息的投影仪标定方法

Country Status (1)

Country Link
CN (1) CN115546311B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112598747A (zh) * 2020-10-15 2021-04-02 武汉易维晟医疗科技有限公司 一种单目相机与投影仪的联合标定方法
CN113160339A (zh) * 2021-05-19 2021-07-23 中国科学院自动化研究所苏州研究院 一种基于沙姆定律的投影仪标定方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4501239B2 (ja) * 2000-07-13 2010-07-14 ソニー株式会社 カメラ・キャリブレーション装置及び方法、並びに、記憶媒体
JP4340923B2 (ja) * 2007-02-23 2009-10-07 セイコーエプソン株式会社 プロジェクタ、プログラムおよび情報記憶媒体
CN101894366B (zh) * 2009-05-21 2014-01-29 北京中星微电子有限公司 一种获取标定参数的方法、装置及一种视频监控系统
CN108230399B (zh) * 2017-12-22 2019-11-08 清华大学 一种基于结构光技术的投影仪标定方法
CN111540004B (zh) * 2020-04-16 2023-07-14 北京清微智能科技有限公司 单相机极线校正方法及装置
CN113298886B (zh) * 2021-07-27 2021-10-08 光量信息科技(宁波)有限公司 一种投影仪的标定方法
CN113686262B (zh) * 2021-08-13 2022-10-11 桂林电子科技大学 一种线结构光扫描仪标定方法、装置及存储介质
CN113920206B (zh) * 2021-12-15 2022-03-08 深圳市艾视铂智能技术有限公司 透视移轴相机的标定方法
CN114842090A (zh) * 2022-04-21 2022-08-02 天津大学 基于精密角度基准的视觉惯性标定系统及其标定方法
CN114792345B (zh) * 2022-06-27 2022-09-27 杭州蓝芯科技有限公司 一种基于单目结构光系统的标定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112598747A (zh) * 2020-10-15 2021-04-02 武汉易维晟医疗科技有限公司 一种单目相机与投影仪的联合标定方法
CN113160339A (zh) * 2021-05-19 2021-07-23 中国科学院自动化研究所苏州研究院 一种基于沙姆定律的投影仪标定方法

Also Published As

Publication number Publication date
CN115546311A (zh) 2022-12-30

Similar Documents

Publication Publication Date Title
CN107767422B (zh) 一种鱼眼镜头的校正方法、装置及便携式终端
CN107749979B (zh) 一种投影机左右梯形校正方法
US5878174A (en) Method for lens distortion correction of photographic images for texture mapping
JP4715334B2 (ja) 車両用画像生成装置および方法
US5898438A (en) Texture mapping of photographic images to CAD surfaces
US20100259624A1 (en) Method and apparatus for calibrating video camera
JP4010754B2 (ja) 画像処理装置と画像処理方法及びコンピュータ読み取り可能な記録媒体
KR20160116075A (ko) 카메라로부터 획득한 영상에 대한 자동보정기능을 구비한 영상처리장치 및 그 방법
CN114727081B (zh) 投影仪投影校正方法、装置及投影仪
CN111540022B (zh) 一种基于虚拟相机的图像一致化方法
WO2021208486A1 (zh) 一种相机坐标变换方法、终端以及存储介质
JPH11118425A (ja) キャリブレーション方法及び装置、並びにキャリブレーション用データ生成方法
JP2009042162A (ja) キャリブレーション装置及びその方法
CN112734860A (zh) 一种基于弧形幕先验信息的逐像素映射投影几何校正方法
CN110555813B (zh) 一种无人机遥感影像的快速几何校正方法及系统
WO2019232793A1 (zh) 双摄像头标定方法、电子设备、计算机可读存储介质
CN110766763A (zh) 一种基于消失点运算的相机标定法
CN113160333B (zh) 参数优化相机标定方法
CN111696141B (zh) 一种三维全景扫描采集方法、设备及存储设备
CN115546311B (zh) 一种基于场景信息的投影仪标定方法
CN109685800B (zh) 一种校准模板和快速标定相机外部参数的方法
CN116743973A (zh) 无感投影图像自动校正方法
CN113962853B (zh) 一种旋转线阵扫描图像位姿自动化精密解算方法
CN113989378A (zh) 一种双目线阵相机位置调节以及标定的方法
JP2006338167A (ja) 画像データ作成方法

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