CN103632360A - 无人机航摄影像的拼接方法 - Google Patents
无人机航摄影像的拼接方法 Download PDFInfo
- Publication number
- CN103632360A CN103632360A CN201310692990.5A CN201310692990A CN103632360A CN 103632360 A CN103632360 A CN 103632360A CN 201310692990 A CN201310692990 A CN 201310692990A CN 103632360 A CN103632360 A CN 103632360A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- coordinate
- reference images
- registration
- images
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及航测领域,特别涉及无人机航摄影像的拼接方法。本发明提供的无人机航摄影像的拼接方法,其得到基准影像的像空间坐标系与物空间坐标系的夹角数据;根据夹角数据形成的旋转矩阵,将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,进而得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,并根据其拼接基准影像与待配准影像,由于待配准影像及基准影像所在的像坐标系均为与物坐标系平行的平面坐标系,因此避免了拼接的影像产生扭曲变形。
Description
技术领域
本发明涉及航测领域,具体而言,涉及无人机航摄影像的拼接方法。
背景技术
近年来,我国自然灾害频发,采用常规的航摄手段周期长费用高,无法及时有效地满足应急测绘、小面积高分辨率地理信息数据更新的要求,而轻小型无人机的航摄***作为传统航空摄影测量手段的有力补充,通过无人机搭载的视频摄像头获取航摄影像,既可以实时获取的视频数据实时查看,获取灾区实时情况,又可以截取视频流数据对其进行拼接,以便快速获取灾区影像,为了解灾情、决策指挥、抢险救灾、灾后恢复重建等工作提供了重要的应急服务保障。具有灵活机动、现势性强、高效快速、作业成本低以及适用范围广等特点,在小区域和飞行困难地区高分辨率影像快速获取方面更具有明显优势。
相关技术中无人机航摄影像的拼接方法,包括:获取无人机上设有的视频摄像头的航摄视频流数据;按照时间间隔、航向重叠度及旁向重叠度截取视频流数据的关键帧影像;将关键帧影像分别设定为基准影像及待配准影像;匹配基准影像的特征点与待配准影像的特征点,获取基准影像与待配准影像中具有重叠度影像的同名点;通过粗差剔除算法剔除同名点中的误配特征点,得到基准影像与待配准影像剔除误配特征点后的同名点;设定基准影像的像坐标系,以及待配准影像的像坐标系;根据剔除误配特征点后的同名点,得到待配准影像的像坐标系转换到基准影像的像坐标系的单应矩阵;通过单应矩阵,得到待配准影像在基准影像的像坐标系中的坐标;获取基准影像的像点在其坐标系中的坐标;根据基准影像的像点在其坐标系中的坐标,以及待配准影像在基准影像的像坐标系中的坐标,拼接待配准影像与基准影像。
但在执行上述无人机航摄影像的拼接方法时,发明人发现至少存在以下问题:由于待配准影像与基准影像所在的像坐标系不同,通过待配准影像到基准影像的单应矩阵,将一次修正的待配准影像,转换至一次修正的基准影像的像坐标系中进行拼接时,拼接的影像会扭曲变形。
发明内容
本发明的目的在于提供无人机航摄影像的拼接方法,以解决上述的问题。
在本发明的实施例中提供了无人机航摄影像的拼接方法,包括:
获取拍摄基准影像的瞬间及多个拍摄待配准影像的瞬间,基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标;
获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标;
根据所述基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标,及所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的夹角数据;
根据所述夹角数据形成旋转矩阵;
根据所述旋转矩阵,将所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像。
进一步地,所述获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,包括:
获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像坐标系中的坐标;
根据拍摄焦距,以及基准影像的像主点在基准影像的像坐标系中的坐标,得到基准影像的像主点在基准影像的像空间坐标系的坐标;
根据拍摄焦距,以及多个待配准影像的像主点分别在基准影像的像坐标系中的坐标,得到多个待配准影像的像主点分别在基准影像的像空间坐标系的坐标。
进一步地,根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,及其对应的GPS天线相位中心在地心坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的旋转矩阵,包括:
预先设定基准影像的像空间坐标系与物空间坐标系的初始夹角数据均为0°,根据所述初始夹角数据形成旋转矩阵初始值:
由
预先设定比例系数的初始值其中H为航摄飞行时的航高,f为摄像机焦距;
根据基准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs1 Ys1 Zs1),以及每一个待配准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs2 Ys2 Zs2),得到每一个平移矩阵的初始值
根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(X Y Z),所述旋转矩阵初始值R0,比例系数的初始值λ0,平移矩阵的初始值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第一次计算概略值比例系数第一次计算概略值dλ1以及平移矩阵第一次计算概略值
根据所述平移矩阵的初始值 及所述平移矩阵第一次计算概略值 得到平移矩阵的第一次迭代值
根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(XY Z),所述旋转矩阵的第一次迭代值R1,比例系数的第一次迭代值λ1,平移矩阵的第一次迭代值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第二次计算概略值、比例系数第二次计算概略值dλ2以及平移矩阵第二次计算概略值
其中θ为预设的阈值;
进一步地,根据所述旋转矩阵,将所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,包括:
根据所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标 以及旋转矩阵R,得到所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标
。
进一步地,根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标包括:
根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标 得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标
进一步地,根据所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像包括:
获取所述基准影像以及所述多个待配准影像中,每两张相邻的影像分别在与物坐标系平行的平面坐标系的同名点的坐标;
根据所述基准影像以及所述多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数;
以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,将所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为所述多个待配准影像的像点,分别在基准坐标系的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在基准坐标系的坐标,拼接基准影像与待配准影像。
进一步地,根据所述基准影像以及所述多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,包括:
对于每一个基准影像,根据所述基准影像以及待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标(xi 1′ yi 1′)及(xi 2′ yi 2′),得到待配准影像与基准影像之间的平移参数(△x,△y):
以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,以及所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,将所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为所述多个待配准影像的像点,分别在基准坐标系的坐标,包括:
对于每一个待配准影像,根据待配准影像与基准影像之间的平移参数(△x,△y)以及待配准影像所有像点分别在与物坐标系平行的平面坐标系的坐标(x’,y’),得到待配准影像所有像点在基准坐标系下的坐标(x”,y”):(x”,y”)=(△x,△y)+(x’,y’)。
本发明实施例提供的无人机航摄影像的拼接方法,其得到基准影像的像空间坐标系与物空间坐标系的夹角数据;根据夹角数据形成的旋转矩阵,将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,进而得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,并根据其拼接基准影像与待配准影像,由于待配准影像及基准影像所在的像坐标系均为与物坐标系平行的平面坐标系,因此避免了拼接的影像产生扭曲变形。
附图说明
图1示出了本发明的无人机航摄影像的拼接方法的一种实施例的示意图。
具体实施方式
本发明中的名词解释如下:
基准影像:影像匹配时,作为匹配基准的参考影像。
待配准影像:影像匹配时,需要将影像坐标纠正到以参考影像为基准的坐标系下的待匹配影像。
特征点:特征点是反映地物类型或区域地理分布特征的点。在地图上具有准确的地理位置和明确的地理属性及含义。大致包含:①独立地物点;②线型要素或面状要素边界线的拐点或折点;③各种线状要素及面状要素边线之间的交叉点等。通过这些点,在地理信息***中建立起不同地理类型间的拓扑关系。
同名点:即同名像点,同一目标点在不同像片上的构像点。
下面通过具体的实施例子并结合附图对本发明做进一步的详细描述。
无人机航摄影像的拼接方法,如图1所示,包括:
S101:获取拍摄基准影像的瞬间及多个拍摄待配准影像的瞬间,基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标;
S102:获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标;
S103:根据基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标,及基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的夹角数据;
S104:根据夹角数据形成旋转矩阵;
S105:根据旋转矩阵,将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标;
S106:根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标;
S107:根据基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像。
其中,获取拍摄基准影像的瞬间及多个拍摄待配准影像的瞬间,基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标,是指:获取拍摄基准影像的瞬间,基准影像对应的GPS天线相位中心在地心坐标系的坐标,以及对于多个待配准影像中的每一个待配准影像,获取拍摄待配准影像的瞬间,待配准影像对应的GPS天线相位中心在地心坐标系的坐标。
本发明实施例提供的无人机航摄影像的拼接方法,其根据基准影像的像空间坐标系与物空间坐标系的夹角数据形成的旋转矩阵,将将所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,从而得到所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,再根据所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像,从而使得待配准影像与基准影像所在的像坐标系相同,避免了拼接的影像产生扭曲变形。
进一步地,在上述的无人机航摄影像的拼接方法中,获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,包括:获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像坐标系中的坐标;根据拍摄焦距,以及基准影像的像主点在基准影像的像坐标系中的坐标,得到基准影像的像主点在基准影像的像空间坐标系的坐标;根据拍摄焦距,以及多个待配准影像的像主点分别在基准影像的像坐标系中的坐标,得到多个待配准影像的像主点分别在基准影像的像空间坐标系的坐标。
通过明确获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标的方法,提高获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标的效率,进而提高拼接影像的效率。
进一步地,在上述的无人机航摄影像的拼接方法中,根据基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,及其对应的GPS天线相位中心在地心坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的旋转矩阵,包括:预先设定基准影像的像空间坐标系与物空间坐标系的初始夹角数据均为0°,根据初始夹角数据形成旋转矩阵初始值:由 其中为得到的基准影像的像空间坐标系与物空间坐标系的夹角数据;预先设定比例系数的初始值其中H为航摄飞行时的航高,f为摄像机焦距;根据基准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs1 Ys1 Zs1),以及每一个待配准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs2 Ys2 Zs2),得到每一个平移矩阵的初始值 根据基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(X Y Z),旋转矩阵初始值R0,比例系数的初始值λ0,平移矩阵的初始值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第一次计算概略值比例系数第一次计算概略值dλ1以及平移矩阵第一次计算概略值 其中, 第一次迭代时:根据基准影像的像空间坐标系与物空间坐标系的夹角数据第一次计算概略值形成旋转矩阵的第一次迭代值R1:根据比例系数的初始值及比例系数第一次计算概略值dλ1,得到比例系数的第一次迭代值λ1:根据平移矩阵的初始值 及平移矩阵第一次计算概略值 得到平移矩阵的第一次迭代值 根据基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(X Y Z),旋转矩阵的第一次迭代值R1,比例系数的第一次迭代值λ1,平移矩阵的第一次迭代值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第二次计算概略值比例系数第二次计算概略值dλ2以及平移矩阵第二次计算概略值 :
通过明确根据基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,及其对应的GPS天线相位中心在地心坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的旋转矩阵的方法,提高得到基准影像的像空间坐标系与物空间坐标系的旋转矩阵的效率,进而提高拼接影像的效率。
进一步地,在上述的无人机航摄影像的拼接方法中,根据旋转矩阵,将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,包括:根据基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标 以及旋转矩阵R,得到基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标
其中,由于航摄中航摄的焦距相同,因此f=f’。
通过明确根据旋转矩阵,将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标的过程,提高了将基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标的效率,进而提高拼接影像的效率。
进一步地,在上述的无人机航摄影像的拼接方法中,根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标包括:根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标 得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标
其中, 是由物空间坐标系与无坐标系平行的平面坐标系的定义可以得到。
通过明确根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标的方法,提高根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标的效率,进而提高拼接影像的效率。
进一步地,在上述的无人机航摄影像的拼接方法中,根据基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像包括:获取基准影像以及多个待配准影像中,每两张相邻的影像分别在与物坐标系平行的平面坐标系的同名点的坐标;根据基准影像以及多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数;以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数,将多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为多个待配准影像的像点,分别在基准坐标系的坐标;根据基准影像的像点以及多个待配准影像的像点,分别在基准坐标系的坐标,拼接基准影像与待配准影像。
通过明确根据基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像的具体拼接方法,将多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,通过待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数,转换为多个待配准影像的像点,分别在基准坐标系的坐标,进而根据基准影像的像点以及多个待配准影像的像点,分别在基准坐标系的坐标,拼接基准影像与待配准影像,提高了拼接基准影像与待配准影像的效率。
进一步地,在上述的无人机航摄影像的拼接方法中,根据基准影像以及多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数,包括:对于每一个基准影像,根据基准影像以及待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标(xi 1′ yi 1′)及(xi 2′ yi 2′),得到待配准影像与基准影像之间的平移参数(△x,△y): 以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,以及待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数,将多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为多个待配准影像的像点,分别在基准坐标系的坐标,包括:对于每一个待配准影像,根据待配准影像与基准影像之间的平移参数(△x,△y)以及待配准影像所有像点分别在与物坐标系平行的平面坐标系的坐标(x’,y’),得到待配准影像所有像点在基准坐标系下的坐标(x”,y”):(x”,y”)=(△x,△y)+(x’,y’)。
通过明确获取平移参数的方法,以及通过平移参数,将多个待配准影像的像点转换至与基准影像的像点同一坐标系的方法,便于本领域技术人员实现根据基准影像以及多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到待配准影像与基准影像在与物坐标系平行的平面坐标系的平移参数,并通过平移参数,得到多个待配准影像的像点转换至与基准影像的像点同一坐标系的坐标,进而提高拼接基准影像与待配准影像的效率。
显然,本领域的技术人员应该明白,上述的本发明的各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.无人机航摄影像的拼接方法,其特征在于,包括:
获取拍摄基准影像的瞬间及多个拍摄待配准影像的瞬间,基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标;
获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标;
根据所述基准影像及多个待配准影像分别对应的GPS天线相位中心在地心坐标系的坐标,及所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的夹角数据;
根据所述夹角数据形成旋转矩阵;
根据所述旋转矩阵,将所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像。
2.根据权利要求1所述的无人机航摄影像的拼接方法,其特征在于,所述获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,包括:
获取基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像坐标系中的坐标;
根据拍摄焦距,以及基准影像的像主点在基准影像的像坐标系中的坐标,得到基准影像的像主点在基准影像的像空间坐标系的坐标;
根据拍摄焦距,以及多个待配准影像的像主点分别在基准影像的像坐标系中的坐标,得到多个待配准影像的像主点分别在基准影像的像空间坐标系的坐标。
3.根据权利要求2所述的无人机航摄影像的拼接方法,其特征在于,根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标,及其对应的GPS天线相位中心在地心坐标系的坐标,得到基准影像的像空间坐标系与物空间坐标系的旋转矩阵,包括:
预先设定基准影像的像空间坐标系与物空间坐标系的初始夹角数据均为0°,根据所述初始夹角数据形成旋转矩阵初始值:
由
预先设定比例系数的初始值其中H为航摄飞行时的航高,f为摄像机焦距;
根据基准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs1 Ys1 Zs1),以及每一个待配准影像对应的GPS天线相位中心在地心坐标系的坐标(Xs2 Ys2 Zs2),得到每一个平移矩阵的初始值
根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(X Y Z),所述旋转矩阵初始值R0,比例系数的初始值λ0,平移矩阵的初始值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第一次计算概略值比例系数第一次计算概略值dλ1以及平移矩阵第一次计算概略值
第一次迭代时:根据基准影像的像空间坐标系与物空间坐标系的夹角数据第一次计算概略值形成旋转矩阵的第一次迭代值R1:
根据所述平移矩阵的初始值 及所述平移矩阵第一次计算概略值 得到平移矩阵的第一次迭代值
根据所述基准影像的像主点以及多个待配准影像的像主点,分别在基准影像的像空间坐标系的坐标(x y -f),以及基准影像的像主点及多个待配准影像的像主点分别对应的GPS天线相位中心在地心坐标系的坐标(X Y Z),所述旋转矩阵的第一次迭代值R1,比例系数的第一次迭代值λ1,平移矩阵的第一次迭代值 得到基准影像的像空间坐标系与物空间坐标系的夹角数据第二次计算概略值比例系数第二次计算概略值dλ2以及平移矩阵第二次计算概略值
其中θ为预设的阈值;
为基准影像的像空间坐标系与物空间坐标系的夹角数据;R为基准影像的像空间坐标系与物空间坐标系的旋转矩阵。
4.根据权利要求3所述的无人机航摄影像的拼接方法,其特征在于,根据所述旋转矩阵,将所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标,转换为所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,包括:
根据所述基准影像的像点以及多个待配准影像的像点,分别在基准影像的像空间坐标系中的坐标 以及旋转矩阵R,得到所述基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标
。
5.根据权利要求4所述的无人机航摄影像的拼接方法,其特征在于,根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标,得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标包括:
根据基准影像的像点以及多个待配准影像的像点,分别在物空间坐标系中的坐标 得到基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标
6.根据权利要求1所述的无人机航摄影像的拼接方法,其特征在于,根据所述基准影像的像点以及多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,拼接基准影像与待配准影像包括:
获取所述基准影像以及所述多个待配准影像中,每两张相邻的影像分别在与物坐标系平行的平面坐标系的同名点的坐标;
根据所述基准影像以及所述多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数;
以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,将所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为所述多个待配准影像的像点,分别在基准坐标系的坐标;
根据所述基准影像的像点以及多个待配准影像的像点,分别在基准坐标系的坐标,拼接基准影像与待配准影像。
7.根据权利要求6所述的无人机航摄影像的拼接方法,其特征在于,根据所述基准影像以及所述多个待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标,得到所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,包括:
对于每一个基准影像,根据所述基准影像以及待配准影像,分别在与物坐标系平行的平面坐标系的同名点的坐标(xi 1′yi 1′)及(xi 2′yi 2′),得到待配准影像与基准影像之间的平移参数(△x,△y):
以基准影像在与物坐标系平行的平面坐标系的坐标为基准建立基准坐标系,根据所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,以及所述待配准影像与所述基准影像在与物坐标系平行的平面坐标系的平移参数,将所述多个待配准影像的像点,分别在与物坐标系平行的平面坐标系的坐标,转换为所述多个待配准影像的像点,分别在基准坐标系的坐标,包括:
对于每一个待配准影像,根据待配准影像与基准影像之间的平移参数(△x,△y)以及待配准影像所有像点分别在与物坐标系平行的平面坐标系的坐标(x’,y’),得到待配准影像所有像点在基准坐标系下的坐标(x”,y”):(x”,y”)=(△x,△y)+(x’,y’)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310692990.5A CN103632360B (zh) | 2013-12-17 | 2013-12-17 | 无人机航摄影像的拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310692990.5A CN103632360B (zh) | 2013-12-17 | 2013-12-17 | 无人机航摄影像的拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103632360A true CN103632360A (zh) | 2014-03-12 |
CN103632360B CN103632360B (zh) | 2016-08-17 |
Family
ID=50213377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310692990.5A Active CN103632360B (zh) | 2013-12-17 | 2013-12-17 | 无人机航摄影像的拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103632360B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105847750A (zh) * | 2016-04-13 | 2016-08-10 | 中测新图(北京)遥感技术有限责任公司 | 基于地理编码的无人机视频影像实时显示的方法及装置 |
CN106897962A (zh) * | 2016-12-26 | 2017-06-27 | 中国科学院长春光学精密机械与物理研究所 | 大视场空间对地观测图像的圆盘投影和拼接方法 |
CN109492525A (zh) * | 2018-09-26 | 2019-03-19 | 中睿通信规划设计有限公司 | 一种测量基站天线工程参数的方法 |
CN113179673A (zh) * | 2018-12-27 | 2021-07-27 | 世界性创新电子合作企业 | 应用多摄像头移动路径追踪技术的影像监视装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20010075948A (ko) * | 2000-01-21 | 2001-08-11 | 오길록 | 영상 특징점 매칭을 이용한 변환 행렬 추정 장치 및 그방법과 그를 이용한 모자이크 영상 생성 장치 및 그 방법 |
-
2013
- 2013-12-17 CN CN201310692990.5A patent/CN103632360B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20010075948A (ko) * | 2000-01-21 | 2001-08-11 | 오길록 | 영상 특징점 매칭을 이용한 변환 행렬 추정 장치 및 그방법과 그를 이용한 모자이크 영상 생성 장치 및 그 방법 |
Non-Patent Citations (1)
Title |
---|
王涛: "线阵CCD传感器实验场几何定标的理论与方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105847750A (zh) * | 2016-04-13 | 2016-08-10 | 中测新图(北京)遥感技术有限责任公司 | 基于地理编码的无人机视频影像实时显示的方法及装置 |
CN105847750B (zh) * | 2016-04-13 | 2019-01-15 | 中测新图(北京)遥感技术有限责任公司 | 基于地理编码的无人机视频影像实时显示的方法及装置 |
CN106897962A (zh) * | 2016-12-26 | 2017-06-27 | 中国科学院长春光学精密机械与物理研究所 | 大视场空间对地观测图像的圆盘投影和拼接方法 |
CN106897962B (zh) * | 2016-12-26 | 2019-06-25 | 中国科学院长春光学精密机械与物理研究所 | 大视场空间对地观测图像的圆盘投影和拼接方法 |
CN109492525A (zh) * | 2018-09-26 | 2019-03-19 | 中睿通信规划设计有限公司 | 一种测量基站天线工程参数的方法 |
CN109492525B (zh) * | 2018-09-26 | 2021-10-01 | 中通服中睿科技有限公司 | 一种测量基站天线工程参数的方法 |
CN113179673A (zh) * | 2018-12-27 | 2021-07-27 | 世界性创新电子合作企业 | 应用多摄像头移动路径追踪技术的影像监视装置 |
CN113179673B (zh) * | 2018-12-27 | 2023-09-19 | 世界性创新电子合作企业 | 应用多摄像头移动路径追踪技术的影像监视装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103632360B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yahyanejad et al. | A fast and mobile system for registration of low-altitude visual and thermal aerial images using multiple small-scale UAVs | |
CN103198524B (zh) | 一种大规模室外场景三维重建方法 | |
CN107808362A (zh) | 一种基于无人机pos信息与图像surf特征结合的图像拼接方法 | |
JP3776787B2 (ja) | 三次元データベース生成システム | |
AU2005271385B2 (en) | Method of preparing a composite image with non-uniform resolution | |
US20170293216A1 (en) | Aerial panoramic oblique photography apparatus | |
CN107527328B (zh) | 一种兼顾精度与速度的无人机影像几何处理方法 | |
JP5134784B2 (ja) | 空中写真測量方法 | |
CN112819903A (zh) | 基于l型标定板的相机和激光雷达联合标定的方法 | |
KR102195051B1 (ko) | 드론의 영상 정보를 이용한 공간 정보 생성 시스템 및 방법과, 이를 위한 컴퓨터 프로그램 | |
KR20120105452A (ko) | 복수의 검출기 어레이를 구비한 다해상도 디지털 대형 카메라 | |
TW202214487A (zh) | 協同3d測繪及表面套合技術 | |
CN103632360A (zh) | 无人机航摄影像的拼接方法 | |
CN109325913A (zh) | 无人机图像拼接方法及装置 | |
He et al. | Automated relative orientation of UAV-based imagery in the presence of prior information for the flight trajectory | |
CN105761257A (zh) | 交叉航带上无人机影像匹配中粗差的剔除方法和装置 | |
Sadeq | Accuracy assessment using different UAV image overlaps | |
JP5304465B2 (ja) | モザイク画像生成方法、装置及びプログラム | |
Moussa et al. | A fast approach for stitching of aerial images | |
Liu et al. | A new approach to fast mosaic UAV images | |
JP3808833B2 (ja) | 空中写真測量方法 | |
JP6155091B2 (ja) | モザイク画像生成装置及び生成方法並びにモザイク画像生成プログラム | |
JP2000276045A (ja) | 撮影画像を用いた地図作成方法、装置および撮影画像の歪み補正方法 | |
CN108195359B (zh) | 空间数据的采集方法及*** | |
KR102426282B1 (ko) | 도로선형 기반의 3차원 공간객체 위치보정 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |