CN113012047B - 动态摄像头坐标映射建立方法、装置及可读存储介质 - Google Patents

动态摄像头坐标映射建立方法、装置及可读存储介质 Download PDF

Info

Publication number
CN113012047B
CN113012047B CN202110327148.6A CN202110327148A CN113012047B CN 113012047 B CN113012047 B CN 113012047B CN 202110327148 A CN202110327148 A CN 202110327148A CN 113012047 B CN113012047 B CN 113012047B
Authority
CN
China
Prior art keywords
camera
longitude
latitude
coordinates
coordinate
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
CN202110327148.6A
Other languages
English (en)
Other versions
CN113012047A (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.)
Guangzhou Fu'an Digital Technology Co ltd
Original Assignee
Guangzhou Fu'an Digital Technology Co ltd
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 Guangzhou Fu'an Digital Technology Co ltd filed Critical Guangzhou Fu'an Digital Technology Co ltd
Priority to CN202110327148.6A priority Critical patent/CN113012047B/zh
Publication of CN113012047A publication Critical patent/CN113012047A/zh
Application granted granted Critical
Publication of CN113012047B publication Critical patent/CN113012047B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038Image mosaicing, e.g. composing plane images from plane sub-images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Studio Devices (AREA)

Abstract

本发明公开了一种动态摄像头坐标映射建立方法、装置及可读存储介质,涉及摄像监控领域。该方案通过利用摄像头的转动角度可以快速得到转动后的坐标映射,仅需一次测定摄像头画面中心位置的经纬度,无需每次都通过经纬仪获取多组标定物的经纬度坐标和摄像头画面坐标,即可完成坐标映射关系的建立,极大的减少了工作量,通过该方法建立的坐标映射,能快速得到摄像头画面中任意待检测目标的经纬度坐标和摄像头画面坐标,能适用于普通交通工具例如船舶的监控过程。

Description

动态摄像头坐标映射建立方法、装置及可读存储介质
技术领域
本发明涉及摄像机监控领域,具体涉及一种动态摄像头坐标映射建立方法、装置及可读存储介质。
背景技术
一般情况下,摄像头在固定时其画面中的各目标点的坐标、经纬度对应关系不变,但是,当摄像头在转换角度后需要对画面进行重新定位,也即重新建立摄像头画面坐标与实际经纬度的映射关系。
现有技术中,公开号为CN111681190A的发明专利申请公开了一种全景视频高精度坐标映射方法,该发明在枪机视频上选取某一点,获取该点的像素点坐标,针对选取的点使用经纬度测试仪进行测绘,得到该点的经纬度坐标,将选取的点的像素点坐标和经纬度信息相关联,形成一一对应的关系。
在实际应用中,动态摄像头每转动一定角度都需重新建立新的坐标映射关系,每次都需要测绘出选取点的经纬度坐标从而达到像素点坐标和经纬度信息相关联,形成对应关系,测绘成本比较高,且方式不灵活。
因此,需要提出一种动态摄像头坐标映射建立方案以适应摄像头的动态过程。
发明内容
为了解决上述技术问题,本发明的第一目的在于提供一种动态摄像头坐标映射建立方法,其能利用摄像头的转动角度可快速得到转动后的坐标映射,通过该方法建立的坐标映射关系能快速得到摄像头画面任意监控目标的经纬度坐标以及摄像头画面坐标;
本发明的第二目的在于提供一种摄像头坐标映射建立装置;
本发明的第三目的在于提供一种可读存储介质。
为实现上述发明目的,本发明采取的技术方案如下:
一种动态摄像头坐标映射建立方法,包括以下步骤:
步骤S1:获取两个必要夹角,两个所述必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;
每转动一次摄像头,执行一次步骤S2~S5;
步骤S2:摄像头转动后,通过摄像头的视场角得到上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,并获得它们对应的摄像头转角;
步骤S3:利用步骤S1中得到的夹角α0、β0以及步骤S2中得到的四个特殊标记点的摄像头转角,通过计算得到四个特殊标记点的经纬度坐标;
步骤S4:利用摄像头参数,计算S2中得到的四个特殊标记点的摄像头画面坐标;
步骤S5:利用S2中得到的四个特殊标记点的经纬度坐标和摄像头画面坐标,计算得到变换矩阵H,由变换矩阵H建立坐标映射。
作为优选,所述步骤S1包括
步骤S1.1:测定O′的经纬度坐标;
步骤S1.2:根据Haversine公式,利用A′和O′的经纬度坐标计算参数d,d为摄像头所在的位置在水平面上的垂直投影位置A′和摄像头画面中心O′的直线水平距离;
步骤S1.3:根据Haversine公式,利用A′和O′的经纬度坐标计算A′和O′的经度水平距离s;
步骤S1.4:根据d、s、h,计算夹角α0、夹角β0,如下:
Figure BDA0002995083840000021
Figure BDA0002995083840000022
其中,h为摄像头距离水平面的高度。
作为优选,所述步骤S1.2为:
Figure BDA0002995083840000023
其中:r是地球半径,单位为m,a为中间变量,a的计算如下:
Figure BDA0002995083840000024
其中:(λ11)为O′的经纬度坐标;(λ00)为A′的经纬度坐标。
作为优选,所述步骤S1.3为:
Figure BDA0002995083840000031
b为中间变量,b的计算如下:
Figure BDA0002995083840000032
作为优选,所述步骤S2中,当摄像头转动后,假设水平方向的转动角度是Δα,垂直方向的转动角度为Δβ,那么转动后摄像头中心线在水平面上的投影和地理真北方向的夹角α为:α=α0+Δα,转动后摄像头中心线与垂直线的夹角β为:β=β0+Δβ。
摄像头视场角的上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点的摄像头画面坐标记为Ii=(xi,yi)(1≤i≤4),它们对应的地理位置为Oi=(λii)(1≤i≤4),λi为经度,ψi为纬度;同时它们对应的摄像头转角Θi=(ωii)(1≤i≤4),ωi为水平转动角度,γi为垂直转动角度。其中Θ1=(α,β+θ/2)、Θ2=(α,β-θ/2)、Θ3=(α-φ/2,β)、Θ4=(α+φ/2,β);φ为水平视场角度,θ为垂直视场角度。
作为优选,所述步骤S3为:
对于步骤S2所得的每个特殊标记点的地理位置Oi(1≤i≤4),计算A′和Oi的直线水平距离di如下:
di=h·tgγi
计算A′和Oi的经度水平距离si如下:
si=di·sinωi
Oi的经纬度坐标(λii)通过以下公式求得:
Figure BDA0002995083840000033
Figure BDA0002995083840000034
其中,M、N、C1和C2均为中间变量,计算如下:
Figure BDA0002995083840000041
Figure BDA0002995083840000042
Figure BDA0002995083840000043
Figure BDA0002995083840000044
作为优选,所述步骤S4为:
对于步骤S2所得的每个特殊标记点的摄像头画面坐标Ii=(xi,yi)(1≤i≤4)计算如下:
x1=X/2,y1=0
x2=X/2,y2=Y
x3=0,y3=Y/2
x4=X,y4=Y/2
其中:X为图像的像素宽度,Y为像素高度,根据摄像头图像分辨率为X×Y,可得X和Y的参数值。
作为优选,所述步骤S5为:
从经纬度坐标和画面坐标组对集合{<Oi,Ii>|1≤i≤4}中,选择三对数据,构成经纬度坐标和摄像头画面坐标对集合Ω={<λkk,xk,yk>|1≤k≤3},其中:(λii)是经纬度坐标,(xi,yi)是画面坐标;
通过逆矩阵计算得到变换矩阵H:
Figure BDA0002995083840000045
基于H,针对画面中任意坐标(x,y)可通过如下变换获得其对应的观察目标的经纬度坐标(λ,ψ)如下:
Figure BDA0002995083840000051
基于同一发明构思,本发明还采取的技术方案如下:
一种动态摄像头坐标映射建立装置,包括:
夹角获取模块,用于获取两个必要夹角,两个所述必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;
转动角度获取模块,用于在每次摄像头转动后获取摄像头画面上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,以及它们对应的摄像头转角;
计算模块,用于通过四个特殊标记点对应的摄像头转角,计算得到它们对应的地理位置的经纬度;
映射模块,用于通过四个特殊标记点的摄像头画面坐标以及他们的经纬度坐标,计算得到变换矩阵H,由变换矩阵H建立摄像头坐标和经纬度坐标的映射。
坐标转换模块,用于基于变换矩阵H计算摄像头画面任意位置对应的经纬度坐标。
基于同一发明构思,本发明还采取的技术方案如下:
一种可读存储介质,其中,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述的动态摄像头坐标映射建立方法的步骤。
相对于现有技术,本发明取得了有益的技术效果:
本发明采用的动态摄像头坐标映射建立方法、装置及可读存储介质,通过利用摄像头的转动角度可以快速得到转动后的坐标映射,仅需一次测定摄像头画面中心位置的经纬度,无需每次都通过经纬仪获取多组标定物的经纬度坐标和摄像头画面坐标,即可完成坐标映射关系的建立,极大的减少了工作量,通过该方法建立的坐标映射,能快速得到摄像头画面中任意待检测目标的经纬度坐标和摄像头画面坐标,能适用于普通交通工具例如船舶的监控过程。
附图说明
图1是本发明实施例1中动态摄像头坐标映射建立方法整体流程图;
图2是本发明实施例1中经纬度坐标和监控摄像头画面坐标的映射关系计算过程示意图一;
图3是本发明实施例1中经纬度坐标和监控摄像头画面坐标的映射关系计算过程示意图二;
图4是本发明实施例1中监控摄像头视场角的示意图一;
图5是本发明实施例1中监控摄像头视场角的示意图二。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例对本发明进行进一步详细说明,但本发明要求保护的范围并不局限于下述具体实施例。
实施例1
参考图1~图5,本实施例公开了一种动态摄像头坐标映射建立方法,用于实现摄像头动态视角中摄像头画面中的监控目标的经纬度坐标与画面坐标的相互对应,该监控目标可以为普通的交通工具,例如船舶、汽车,还可以是人、动物或建筑物;
该动态摄像头坐标映射建立方法包括以下步骤:
步骤S1:获取两个必要夹角,两个必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;
每转动一次摄像头,执行一次步骤S2~S5;
步骤S2:摄像头转动后,通过摄像头的视场角得到上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,并获得它们对应的摄像头转角;
步骤S3:利用步骤S1中得到的夹角α0、β0以及步骤S2中得到的四个特殊标记点的摄像头转角,通过计算得到四个特殊标记点的经纬度坐标;
步骤S4:利用摄像头参数,计算S2中得到的四个特殊标记点的摄像头画面坐标;
步骤S5:利用S2中得到的四个特殊标记点的经纬度坐标和摄像头画面坐标,计算得到变换矩阵H,由变换矩阵H建立坐标映射。
更优地,参考图2~图3,步骤S1包括:
步骤S1.1:测定O′的经纬度坐标;
步骤S1.2:根据Haversine公式,利用A′和O′的经纬度坐标计算参数d,d为摄像头所在的位置在水平面上的垂直投影位置A′和摄像头画面中心O′的直线水平距离,单位为m;
步骤S1.3:根据Haversine公式,利用A′和O′的经纬度坐标计算A′和O′的经度水平距离s,单位为m;
步骤S1.4:根据d、s、h,计算夹角α0、夹角β0,如下:
Figure BDA0002995083840000071
Figure BDA0002995083840000072
其中,h为摄像头距离水平面的高度,单位为m。
更优地,步骤S1.2为:
Figure BDA0002995083840000073
其中:r是地球半径,单位为m,a为中间变量,a的计算如下:
Figure BDA0002995083840000074
其中:(λ11)为O′的经纬度坐标;(λ00)为A′的经纬度坐标。
更优地,步骤S1.3为:
Figure BDA0002995083840000075
b为中间变量,b的计算如下:
Figure BDA0002995083840000076
更优地,参考图4~图5,步骤S2中,当摄像头转动后,假设水平方向的转动角度是Δα,垂直方向的转动角度为Δβ,那么转动后摄像头中心线在水平面上的投影和地理真北方向的夹角α为:α=α0+Δα,转动后摄像头中心线与垂直线的夹角β为:β=β0+Δβ。
摄像头视场角的上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点的摄像头画面坐标记为Ii=(xi,yi)(1≤i≤4),它们对应的地理位置为Oi=(λii)(1≤i≤4),λi为经度,ψi为纬度;同时它们对应的摄像头转角Θi=(ωii)(1≤i≤4),ωi为水平转动角度,γi为垂直转动角度。其中Θ1=(α,β+θ/2)、Θ2=(α,β-θ/2)、Θ3=(α-φ/2,β)、Θ4=(α+φ/2,β);φ为水平视场角度,θ为垂直视场角度。
更优地,步骤S3为:
对于步骤S2所得的每个特殊标记点的地理位置Oi(1≤i≤4),计算A′和Oi的直线水平距离di如下:
di=h·tgγi
计算A′和Oi的经度水平距离si如下:
si=di·sinωi
Oi的经纬度坐标(λii)通过以下公式求得:
Figure BDA0002995083840000081
Figure BDA0002995083840000082
其中,M、N、C1和C2均为中间变量,计算如下:
Figure BDA0002995083840000083
Figure BDA0002995083840000084
Figure BDA0002995083840000085
Figure BDA0002995083840000086
更优地,步骤S4为:
对于步骤S2所得的每个特殊标记点的摄像头画面坐标Ii=(xi,yi)(1≤i≤4)计算如下:
x1=X/2,y1=0
x2=X/2,y2=Y
x3=0,y3=Y/2
x4=X,y4=Y/2
其中:X为图像的像素宽度,Y为像素高度,根据摄像头图像分辨率为X×Y,可得X和Y的参数值。
更优地,步骤S5为:
从经纬度坐标和画面坐标组对集合{<Oi,Ii>|1≤i≤4}中,选择三对数据,构成经纬度坐标和摄像头画面坐标对集合Ω={<λkk,xk,yk>|1≤k≤3},其中:(λii)是经纬度坐标,(xi,yi)是画面坐标;
通过逆矩阵计算得到变换矩阵H:
Figure BDA0002995083840000091
基于H,针对画面中任意坐标(x,y)可通过如下变换获得其对应的观察目标的经纬度坐标(λ,ψ)如下:
Figure BDA0002995083840000092
实施例2
本实施例公开了一种动态摄像头坐标映射建立装置,用于实现实施例1的方法,包括:
夹角获取模块,用于获取两个必要夹角,两个必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;
转动角度获取模块,用于在每次摄像头转动后获取摄像头画面上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,以及它们对应的摄像头转角;
计算模块,用于通过四个特殊标记点对应的摄像头转角,计算得到它们对应的地理位置的经纬度;
映射模块,用于通过四个特殊标记点的摄像头画面坐标以及他们的经纬度坐标,计算得到变换矩阵H,由变换矩阵H建立摄像头坐标和经纬度坐标的映射。
坐标转换模块,用于基于变换矩阵H计算摄像头画面任意位置对应的经纬度坐标。
以上模块可基于现有的微处理器设计,例如计算模块可设置于一微处理器中,微处理器可采用而不限于工控机、单片机或PLC,夹角获取模块、转动角度获取模块、计算模块、映射模块以及坐标转换模块均可基于实施例1公开的动态摄像头坐标映射建立方法进行程序设计,并按照该方法的步骤依次执行。
另外,在夹角获取模块中,可设置为通过人工测定后向微处理器输入参数或者通过自动测量夹角的设备测量后自动向微处理器输出参数。
实施例3
本实施例公开一种可读存储介质,该可读存储介质上存储有计算机程序,计算机程序被处理器执行时实现实施例1的动态摄像头坐标映射建立方法的步骤。
根据上述说明书的揭示和教导,本发明所属领域的技术人员还可以对上述实施方式进行变更和修改。因此,本发明并不局限于上面揭示和描述的具体实施方式,对发明的一些修改和变更也应当落入本发明的权利要求的保护范围内。此外,尽管本说明书中使用了一些特定的术语,但这些术语只是为了方便说明,并不对发明构成任何限制。

Claims (3)

1.一种动态摄像头坐标映射建立方法,其特征在于,包括以下步骤:
步骤S1:获取两个必要夹角,两个所述必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;
每转动一次摄像头,执行一次步骤S2~S5;
步骤S2:摄像头转动后,通过摄像头的视场角得到上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,并获得它们对应的摄像头转角;
步骤S3:利用步骤S1中得到的夹角α0、β0以及步骤S2中得到的四个特殊标记点的摄像头转角,通过计算得到四个特殊标记点的经纬度坐标;
步骤S4:利用摄像头参数,计算S2中得到的四个特殊标记点的摄像头画面坐标;
步骤S5:利用S2中得到的四个特殊标记点的经纬度坐标和摄像头画面坐标,计算得到变换矩阵H,由变换矩阵H建立坐标映射;
所述步骤S1包括
步骤S1.1:测定O′的经纬度坐标;
步骤S1.2:根据Haversine公式,利用A′和O′的经纬度坐标计算参数d,d为摄像头所在的位置在水平面上的垂直投影位置A′和摄像头画面中心O′的直线水平距离;
步骤S1.3:根据Haversine公式,利用A′和O′的经纬度坐标计算A′和O′的经度水平距离s;
步骤S1.4:根据d、s、h,计算夹角α0、夹角β0,如下:
Figure FDA0003318824980000011
Figure FDA0003318824980000012
其中,h为摄像头距离水平面的高度;
所述步骤S1.2为:
Figure FDA0003318824980000013
其中:r是地球半径,单位为m,a为中间变量,a的计算如下:
Figure FDA0003318824980000014
其中:(λ11)为O′的经纬度坐标;(λ00)为A′的经纬度坐标;
所述步骤S1.3为:
Figure FDA0003318824980000021
b为中间变量,b的计算如下:
Figure FDA0003318824980000022
所述步骤S2中,当摄像头转动后,假设水平方向的转动角度是Δα,垂直方向的转动角度为Δβ,那么转动后摄像头中心线在水平面上的投影和地理真北方向的夹角α为:α=α0+Δα,转动后摄像头中心线与垂直线的夹角β为:β=β0+Δβ;
摄像头视场角的上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点的摄像头画面坐标记为Ii=(xi,yi)(1≤i≤4),它们对应的地理位置为Oi=(λii)(1≤i≤4),λi为经度,ψi为纬度;同时它们对应的摄像头转角Θi=(ωii)(1≤i≤4),ωi为水平转动角度,γi为垂直转动角度;其中Θ1=(α,β+θ/2)、Θ2=(α,β-θ/2)、Θ3=(α-φ/2,β)、Θ4=(α+φ/2,β);φ为水平视场角度,θ为垂直视场角度;
所述步骤S3为:
对于步骤S2所得的每个特殊标记点的地理位置Oi(1≤i≤4),计算A′和Oi的直线水平距离di如下:
di=h·tgγi
计算A′和Oi的经度水平距离si如下:
si=di·sinωi
Oi的经纬度坐标(λii)通过以下公式求得:
Figure FDA0003318824980000023
Figure FDA0003318824980000031
其中,M、N、C1和C2均为中间变量,计算如下:
Figure FDA0003318824980000032
Figure FDA0003318824980000033
Figure FDA0003318824980000034
Figure FDA0003318824980000035
所述步骤S4为:
对于步骤S2所得的每个特殊标记点的摄像头画面坐标Ii=(xi,yi)(1≤i≤4)计算如下:
x1=X/2,y1=0
x2=X/2,y2=Y
x3=0,y3=Y/2
x4=X,y4=Y/2
其中:X为图像的像素宽度,Y为像素高度,根据摄像头图像分辨率为X×Y,可得X和Y的参数值;
所述步骤S5为:
从经纬度坐标和画面坐标组对集合{<Oi,Ii>|1≤i≤4}中,选择三对数据,构成经纬度坐标和摄像头画面坐标对集合Ω={<λkk,xk,yk>|1≤k≤3},其中:(λii)是经纬度坐标,(xi,yi)是画面坐标;
通过逆矩阵计算得到变换矩阵H:
Figure FDA0003318824980000036
基于H,针对画面中任意坐标(x,y)可通过如下变换获得其对应的观察目标的经纬度坐标(λ,ψ)如下:
Figure FDA0003318824980000041
2.一种动态摄像头坐标映射建立装置,其特征在于,包括:
夹角获取模块,用于获取两个必要夹角,两个所述必要夹角分别为A′和O′的连线和地理真北方向的夹角α0、A和O′的连线和垂直线的夹角β0,其中,摄像头所在的位置A在水平面的垂直投影位置为A′,摄像头画面中心位置为O′;获取夹角α0和夹角β0的步骤包括:步骤S1.1:测定O′的经纬度坐标;
步骤S1.2:根据Haversine公式,利用A′和O′的经纬度坐标计算参数d,d为摄像头所在的位置在水平面上的垂直投影位置A′和摄像头画面中心O′的直线水平距离;
步骤S1.3:根据Haversine公式,利用A′和O′的经纬度坐标计算A′和O′的经度水平距离s;
步骤S1.4:根据d、s、h,计算夹角α0、夹角β0,如下:
Figure FDA0003318824980000042
Figure FDA0003318824980000043
其中,h为摄像头距离水平面的高度;
所述步骤S1.2为:
Figure FDA0003318824980000044
其中:r是地球半径,单位为m,a为中间变量,a的计算如下:
Figure FDA0003318824980000045
其中:(λ11)为O′的经纬度坐标;(λ00)为A′的经纬度坐标;
所述步骤S1.3为:
Figure FDA0003318824980000051
b为中间变量,b的计算如下:
Figure FDA0003318824980000052
转动角度获取模块,用于在每次摄像头转动后获取摄像头画面上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点,以及它们对应的摄像头转角;当摄像头转动后,假设水平方向的转动角度是Δα,垂直方向的转动角度为Δβ,那么转动后摄像头中心线在水平面上的投影和地理真北方向的夹角α为:α=α0+Δα,转动后摄像头中心线与垂直线的夹角β为:β=β0+Δβ;摄像头视场角的上边界的中心点、下边界的中心点、左边界的中心点和右边界的中心点四个特殊标记点的摄像头画面坐标记为Ii=(xi,yi)(1≤i≤4),它们对应的地理位置为Oi=(λii)(1≤i≤4),λi为经度,ψi为纬度;同时它们对应的摄像头转角Θi=(ωii)(1≤i≤4),ωi为水平转动角度,γi为垂直转动角度;其中Θ1=(α,β+θ/2)、Θ2=(α,β-θ/2)、Θ3=(α-φ/2,β)、Θ4=(α+φ/2,β);φ为水平视场角度,θ为垂直视场角度;
计算模块,用于通过四个特殊标记点对应的摄像头转角,计算得到它们对应的地理位置的经纬度;对于转动角度获取模块所得的每个特殊标记点的地理位置Oi(1≤i≤4),计算A′和Oi的直线水平距离di如下:
di=h·tgγi
计算A′和Oi的经度水平距离si如下:
si=di·sinωi
Oi的经纬度坐标(λii)通过以下公式求得:
Figure FDA0003318824980000053
Figure FDA0003318824980000061
其中,M、N、C1和C2均为中间变量,计算如下:
Figure FDA0003318824980000062
Figure FDA0003318824980000063
Figure FDA0003318824980000064
Figure FDA0003318824980000065
映射模块,用于通过四个特殊标记点的摄像头画面坐标以及他们的经纬度坐标,计算得到变换矩阵H,由变换矩阵H建立摄像头坐标和经纬度坐标的映射;对于转动角度获取模块所得的每个特殊标记点的摄像头画面坐标Ii=(xi,yi)(1≤i≤4)计算如下:
x1=X/2,y1=0
x2=X/2,y2=Y
x3=0,y3=Y/2
x4=X,y4=Y/2
其中:X为图像的像素宽度,Y为像素高度,根据摄像头图像分辨率为X×Y,可得X和Y的参数值;
坐标转换模块,用于基于变换矩阵H计算摄像头画面任意位置对应的经纬度坐标;从经纬度坐标和画面坐标组对集合{<Oi,Ii>|1≤i≤4}中,选择三对数据,构成经纬度坐标和摄像头画面坐标对集合Ω={<λkk,xk,yk>|1≤k≤3},其中:(λii)是经纬度坐标,(xi,yi)是画面坐标;
通过逆矩阵计算得到变换矩阵H:
Figure FDA0003318824980000071
基于H,针对画面中任意坐标(x,y)可通过如下变换获得其对应的观察目标的经纬度坐标(λ,ψ)如下:
Figure FDA0003318824980000072
3.一种可读存储介质,其特征在于,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1中任一项所述的动态摄像头坐标映射建立方法的步骤。
CN202110327148.6A 2021-03-26 2021-03-26 动态摄像头坐标映射建立方法、装置及可读存储介质 Active CN113012047B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110327148.6A CN113012047B (zh) 2021-03-26 2021-03-26 动态摄像头坐标映射建立方法、装置及可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110327148.6A CN113012047B (zh) 2021-03-26 2021-03-26 动态摄像头坐标映射建立方法、装置及可读存储介质

Publications (2)

Publication Number Publication Date
CN113012047A CN113012047A (zh) 2021-06-22
CN113012047B true CN113012047B (zh) 2021-12-14

Family

ID=76407871

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110327148.6A Active CN113012047B (zh) 2021-03-26 2021-03-26 动态摄像头坐标映射建立方法、装置及可读存储介质

Country Status (1)

Country Link
CN (1) CN113012047B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113569647B (zh) * 2021-06-29 2024-02-20 广州赋安数字科技有限公司 一种基于ais的船舶高精度坐标映射方法
CN113645449B (zh) * 2021-08-10 2024-05-24 深圳市玄羽科技有限公司 工业互联网的数据采集方法、系统及计算机可读存储介质
CN113642463B (zh) * 2021-08-13 2023-03-10 广州赋安数字科技有限公司 一种视频监控和遥感图像的天地多视图对齐方法
CN113838126B (zh) * 2021-09-27 2022-05-10 广州赋安数字科技有限公司 一种视频监控和无人机的图像对齐方法
CN115442531A (zh) * 2022-09-05 2022-12-06 深圳奇迹智慧网络有限公司 一种用于确定摄像头方向的方法以及摄像头方向确定装置
CN116993830B (zh) * 2023-08-17 2024-09-27 广州赋安数字科技有限公司 一种动态摄像头坐标映射的自动校准方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108810473A (zh) * 2018-06-15 2018-11-13 高新兴科技集团股份有限公司 一种在移动平台上实现gps映射摄像机画面坐标的方法及系统
WO2020114231A1 (zh) * 2018-12-06 2020-06-11 杭州海康威视数字技术股份有限公司 一种基于gps的目标跟踪系统、方法及球机
CN111275609A (zh) * 2018-12-05 2020-06-12 杭州海康威视数字技术股份有限公司 一种pt坐标与gps坐标的转换方法、装置及球机
CN111681190A (zh) * 2020-06-18 2020-09-18 深圳天海宸光科技有限公司 一种全景视频高精度坐标映射方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289360A (zh) * 2011-08-25 2011-12-21 浙江大学 一种自适应投影颜色补偿方法
CN106156199B (zh) * 2015-04-22 2022-04-08 清华大学 一种视频监控图像存储检索方法
CN105072414B (zh) * 2015-08-19 2019-03-12 浙江宇视科技有限公司 一种目标检测和跟踪方法及系统
CN105205852B (zh) * 2015-10-27 2017-10-17 中国电子科技集团公司第二十八研究所 一种基于多尺度渲染和拟合的三维船舶动态展现方法
CN105825242B (zh) * 2016-05-06 2019-08-27 南京大学 基于混合网格分层聚类的集群通信终端轨迹实时异常检测方法与系统
CN107808362A (zh) * 2017-11-15 2018-03-16 北京工业大学 一种基于无人机pos信息与图像surf特征结合的图像拼接方法
CN108510445A (zh) * 2018-03-30 2018-09-07 长沙全度影像科技有限公司 一种全景图像拼接方法
CN109523471B (zh) * 2018-11-16 2023-11-10 厦门博聪信息技术有限公司 一种地面坐标和广角摄像机画面坐标的转换方法、系统以及装置
CN111815672B (zh) * 2019-04-12 2023-10-24 浙江宇视科技有限公司 动态跟踪控制方法、装置及控制设备
CN112129281B (zh) * 2019-06-25 2022-02-22 南京航空航天大学 一种基于局部邻域地图的高精度图像导航定位方法
CN110823233B (zh) * 2019-11-28 2021-03-02 广东电网有限责任公司 一种基于坐标变换的地图模型构建方法及系统
CN111127318B (zh) * 2019-12-05 2023-04-07 南京莱斯电子设备有限公司 一种机场环境下的全景图像拼接方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108810473A (zh) * 2018-06-15 2018-11-13 高新兴科技集团股份有限公司 一种在移动平台上实现gps映射摄像机画面坐标的方法及系统
CN111275609A (zh) * 2018-12-05 2020-06-12 杭州海康威视数字技术股份有限公司 一种pt坐标与gps坐标的转换方法、装置及球机
WO2020114231A1 (zh) * 2018-12-06 2020-06-11 杭州海康威视数字技术股份有限公司 一种基于gps的目标跟踪系统、方法及球机
CN111681190A (zh) * 2020-06-18 2020-09-18 深圳天海宸光科技有限公司 一种全景视频高精度坐标映射方法

Also Published As

Publication number Publication date
CN113012047A (zh) 2021-06-22

Similar Documents

Publication Publication Date Title
CN113012047B (zh) 动态摄像头坐标映射建立方法、装置及可读存储介质
CN110057295B (zh) 一种免像控的单目视觉平面距离测量方法
CN107507235B (zh) 一种基于rgb-d设备采集的彩色图像和深度图像的配准方法
CN109764858B (zh) 一种基于单目相机的摄影测量方法及系统
CN101813465B (zh) 非接触精密测转角的单目视觉测量方法
CN105118042B (zh) 对位贴合方法和设备
CN102136140B (zh) 一种基于矩形图样的视频图像距离检测方法
CN113569647B (zh) 一种基于ais的船舶高精度坐标映射方法
CN111220129A (zh) 一种带旋转云台的聚焦测量方法及终端
CN111220120B (zh) 一种动平台双目测距自校准方法及装置
CN106898025A (zh) 一种基于八点编码标志的相机位移变换矩阵标定方法
CN112200875B (zh) 非量测相机交叉耦合误差补偿与影像匹配纠正方法及系统
CN116597013B (zh) 一种基于不同经纬度区域的卫星图像几何定标方法
CN104634248A (zh) 一种双目视觉下的转轴标定方法
CN104655106B (zh) 基于gps rtk和全景影像的自主定位定向测图方法
CN110836662B (zh) 基于相对定向和绝对定向算法的边坡位移监测方法
CN114792345B (zh) 一种基于单目结构光系统的标定方法
CN114663463A (zh) 关节活动度的测量方法、系统、设备、电子设备和存储介质
CN102012213A (zh) 单副图像测量前景高度的方法
JP3282516B2 (ja) 三次元画像計測方法
CN115187612A (zh) 一种基于机器视觉的平面面积测量方法、装置及系统
CN111968182A (zh) 一种双目相机非线性模型参数的标定方法
CN109579796B (zh) 一种投影后影像的区域网平差方法
CN111932642B (zh) 一种结构裂纹体积的测算方法、装置、设备及存储介质
CN109035335B (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
CB02 Change of applicant information

Address after: 510000 room 1301 (Location: room 1301-1), No. 68, yueken Road, Wushan street, Tianhe District, Guangzhou City, Guangdong Province (office only)

Applicant after: Guangzhou Fu'an Digital Technology Co.,Ltd.

Address before: 510000 No. 1501, 68 yueken Road, Tianhe District, Guangzhou City, Guangdong Province

Applicant before: GUANGZHOU FUAN ELECTRONIC TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant