CN115439365A - 高塔相机图像几何校正方法和装置 - Google Patents
高塔相机图像几何校正方法和装置 Download PDFInfo
- Publication number
- CN115439365A CN115439365A CN202211087790.2A CN202211087790A CN115439365A CN 115439365 A CN115439365 A CN 115439365A CN 202211087790 A CN202211087790 A CN 202211087790A CN 115439365 A CN115439365 A CN 115439365A
- Authority
- CN
- China
- Prior art keywords
- image
- panoramic
- rectangular image
- resampling
- row
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012952 Resampling Methods 0.000 claims abstract description 72
- 238000011426 transformation method Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims description 17
- 230000003287 optical effect Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 abstract description 14
- 238000003384 imaging method Methods 0.000 abstract description 12
- 230000000007 visual effect Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008569 process Effects 0.000 abstract description 6
- 241001465754 Metazoa Species 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种高塔相机图像几何校正方法和装置,属于遥感技术领域。本发明首先获取高塔相机360度拍摄得到的全景矩形图像,然后计算全景矩形图像每一列像素的重采样比例系数,并将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,然后将重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。本发明在将全景矩形图像变换为圆环形真实图像照片的过程中,通过重新计算全景矩形图像每一列像素的重采样比例系数并进行重采样,解决了高塔相机成像尺度随视野距离增长呈非线性变化的问题,提高了高塔相机图像几何校正的准确性,高塔相机成像数据在人类活动监控、植被生长监测、大型动物跟踪监测等方面的应用研究奠定了基础。
Description
技术领域
本发明涉及遥感技术领域,特别是指一种高塔相机图像几何校正方法和装置。
背景技术
在通信铁塔等高塔上安装光学相机等传感器用于高塔周边区域人类活动监控、植被生长监测等工作,具有时效快、分辨率高、成本低等特点,是生态环境监管的另一种技术手段。高塔光学相机通过扫描拍摄的方式,获得视野范围内360度的周边地物照片的全景矩形照片。通过几何校正将高塔光学相机采集的原始矩形全景图像转换为空间地理位置和尺度准确的圆环状图像,是高塔光学相机开展进一步业务应用的前提。
目前,国内外专门针对高塔相机的几何校正方法几乎空白。传统的图像处理方法通常是将矩形图片通过坐标变换的方式转换为圆环形状,这种方式只考虑了常规图片形状上的变换,照片重采样尺度不变。但是高塔相机的成像尺度随视野距离增长呈非线性变化,因此传统重采样尺度不变的坐标变换方式不能满足高塔相机照片尺度校正的需求。
发明内容
为解决上述技术问题,本发明提供一种高塔相机图像几何校正方法和装置,通过重新计算全景矩形图像每一列像素的重采样比例系数并进行重采样,解决了高塔相机成像尺度随视野距离增长呈非线性变化的问题,提高了高塔相机图像几何校正的准确性。
本发明提供技术方案如下:
一种高塔相机图像几何校正方法,所述方法包括:
S1:获取高塔相机360度拍摄得到的全景矩形图像;其中,所述全景矩形图像的行表示所述高塔相机360度拍摄的圆环形区域的径向,所述全景矩形图像的列表示所述高塔相机360度拍摄的圆环形区域的周向;
Rring=Rout-Rin
S3:计算所述全景矩形图像每一列像素的重采样比例系数Scalei;
其中,Scalei为所述全景矩形图像第i列像素的重采样比例系数,i=1,2,…,OriImgcol,OriImgcol为所述全景矩形图像的列数,Pixceli为所述全景矩形图像第i列像素的采样尺度;
S4:将所述全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像;
S5:将所述重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。
进一步的,所述S5包括:
S51:计算所述重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Pring和圆环形区域的外圆半径Rout的像素数Pout;
S52:设定所述圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
S53:将所述重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到所述圆环形图像对应的坐标点(Newx,Newy);
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol;
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow;
ScaleImgrow和ScaleImgcol分别为所述重采样全景矩形图像的行数和列数;
S54:将所述重采样全景矩形图像的坐标点(row,col)的像素值赋值到所述圆环形图像对应的坐标点(Newx,Newy),得到所述圆环形图像。
进一步的,所述S4包括:
将所述全景矩形图像的行保持不变,按照每一列的重采样比例系数逐列进行重采样,并依次追加在前一列后面,直至所有列重采样完毕为止,得到所述重采样全景矩形图像。
进一步的,所述S1包括:
获取高塔相机360度拍摄的图像,并去除无效区域,得到所述全景矩形图像;其中,所述无效区域包括开始拍摄和结束拍摄时的重复拍摄区域以及天际线周边及以外区域。
进一步的,所述方法还包括:
S6:将所述圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
进一步的,所述S6包括:
S61:以所述遥感影像为基准,将所述圆环形图像的圆心的经纬度坐标设定为所述高塔相机中心点在地面的投影点的经纬度坐标;
S62:在所述圆环形图像和所述遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对;
S63:通过几何校正,进行圆环形图像投影变换和地理配准,使得所述圆环形图像和所述遥感影像的几何校正控制点对对准,得到空间位置准确的高塔相机图像。
一种高塔相机图像几何校正装置,所述装置包括:
全景矩形图像获取模块,用于获取高塔相机360度拍摄得到的全景矩形图像;其中,所述全景矩形图像的行表示所述高塔相机360度拍摄的圆环形区域的径向,所述全景矩形图像的列表示所述高塔相机360度拍摄的圆环形区域的周向;
Rring=Rout-Rin
第二计算模块,用于计算所述全景矩形图像每一列像素的重采样比例系数Scalei;
其中,Scalei为所述全景矩形图像第i列像素的重采样比例系数,i=1,2,…,OriImgcol,OriImgcol为所述全景矩形图像的列数,Pixceli为所述全景矩形图像第i列像素的采样尺度;
重采样模块,用于将所述全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像;
图像转换模块,用于将所述重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。
进一步的,所述图像转换模块包括:
第一计算单元,用于计算所述重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Rring和圆环形区域的外圆半径Rout的像素数Pout;
中心点设定单元,用于设定所述圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
坐标变换单元,用于将所述重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到所述圆环形图像对应的坐标点(Newx,Newy);
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol;
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow;
ScaleImgrow和ScaleImgcol分别为所述重采样全景矩形图像的行数和列数;
赋值单元,用于将所述重采样全景矩形图像的坐标点(row,col)的像素值赋值到所述圆环形图像对应的坐标点(Newx,Newy),得到所述圆环形图像。
进一步的,所述重采样模块用于:
将所述全景矩形图像的行保持不变,按照每一列的重采样比例系数逐列进行重采样,并依次追加在前一列后面,直至所有列重采样完毕为止,得到所述重采样全景矩形图像。
进一步的,所述全景矩形图像获取模块用于:
获取高塔相机360度拍摄的图像,并去除无效区域,得到所述全景矩形图像;其中,所述无效区域包括开始拍摄和结束拍摄时的重复拍摄区域以及天际线周边及以外区域。
进一步的,所述装置还包括:
地理配准模块,用于将所述圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
进一步的,所述地理配准模块包括:
中心点匹配单元,用于以所述遥感影像为基准,将所述圆环形图像的圆心的经纬度坐标设定为所述高塔相机中心点在地面的投影点的经纬度坐标;
几何校正控制点对设置单元,用于在所述圆环形图像和所述遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对;
几何校正单元,用于通过几何校正,进行圆环形图像投影变换和地理配准,使得所述圆环形图像和所述遥感影像的几何校正控制点对对准,得到空间位置准确的高塔相机图像。
本发明具有以下有益效果:
本发明首先获取高塔相机360度拍摄得到的全景矩形图像,然后计算全景矩形图像每一列像素的重采样比例系数,并将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,然后将重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。本发明在将全景矩形图像变换为圆环形真实图像照片的过程中,通过重新计算全景矩形图像每一列像素的重采样比例系数并进行重采样,解决了高塔相机成像尺度随视野距离增长呈非线性变化的问题,提高了高塔相机图像几何校正的准确性,高塔相机成像数据在人类活动监控、植被生长监测、大型动物跟踪监测等方面的应用研究奠定了基础。
附图说明
图1为本发明的高塔相机图像几何校正方法的流程图;
图2为全景矩形图像的示例图;
图3为图2所示的全景矩形图像几何校正得到的圆环形图像的示例图;
图4为高塔相机成像的简示图;
图5为本发明的高塔相机图像几何校正装置的示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
实施例1:
本发明实施例提供一种高塔相机图像几何校正方法,如图1所示,该方法包括:
S1:获取高塔相机360度拍摄得到的全景矩形图像;其中,该全景矩形图像的行表示高塔相机360度拍摄的圆环形区域的径向,全景矩形图像的列表示高塔相机360度拍摄的圆环形区域的周向。
本步骤用于获取待研究区域的全景矩形图像,全景矩形图像的获取方式不受限制,其中的一个实现方式为:
高塔相机360度拍摄图像后,根据待研究区域裁切掉无效区域,得到待研究区域标准化的全景矩形图像,记为OriImg,全景矩形图像的一个示例如图2所示。所述的无效区域包括相机开始拍摄和结束拍摄时的重复拍摄区域、天际线周边及以外区域、其他无效区域等。高塔相机360度拍摄的全景矩形图像对应的圆环形区域如图3所示。
本发明可以适用于装载在高塔上的各种相机的图像,例如全色相机、多光谱相机、高光谱相机、热红外相机及其他光学相机的图像。
假设高塔周边的观测区域为平坦地区,相机在高塔上装载,高塔相机的光轴相对于水平面的俯视角度为α,高塔相机的视场角为高塔相机距离地面的垂直高度为h。由于高塔观测区域相对较小,因此可以认为观测成像面为水平平面。
Rring=Rout-Rin
S3:计算全景矩形图像每一列像素的重采样比例系数Scalei。
本步骤中,需要首先获取全景矩形图像的行数OriImgrow和列数OriImgcol。全景矩形图像的行数OriImgrow也即为全景矩形图像对应的圆环形区域的内圆的像素个数,全景矩形图像的列数OriImgcol也即为全景矩形图像对应的圆环形区域的圆环环宽的像素数。
再计算全景矩形图像第i列像素的采样尺度Pixceli,Pixceli也即为全景矩形图像对应的圆环形区域的相邻圈像素的采样尺度。
其中,i=1,2,…,OriImgcol。
最后计算全景矩形图像每一列像素的重采样比例系数Scalei。
其中,Scalei为全景矩形图像第i列像素的重采样比例系数,Scalei也即为全景矩形图像对应的圆环形区域的第i圈像素的重采样比例系数。
S4:将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像。
本步骤的一个具体实现方式为:将全景矩形图像的行保持不变,按照每一列的重采样比例系数逐列进行重采样,并依次追加在前一列后面,直至所有列重采样完毕为止,得到重采样全景矩形图像ScaleImg,重采样全景矩形图像ScaleImg即为尺度校准后的图像。
S5:将重采样全景矩形图像ScaleImg通过坐标变换的方法转换为圆环形图像CircleImg。
图2所示的全景矩形图像经过本发明的高塔相机图像几何校正方法处理后,得到图3所示的圆环形图像。
本发明首先获取高塔相机360度拍摄得到的全景矩形图像,然后计算全景矩形图像每一列像素的重采样比例系数,并将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,然后将重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。本发明在将全景矩形图像变换为圆环形真实图像照片的过程中,通过重新计算全景矩形图像每一列像素的重采样比例系数并进行重采样,解决了高塔相机成像尺度随视野距离增长呈非线性变化的问题,提高了高塔相机图像几何校正的准确性,高塔相机成像数据在人类活动监控、植被生长监测、大型动物跟踪监测等方面的应用研究奠定了基础。
作为本发明实施例的一种改进,前述的S5包括:
S51:计算所述重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Rring和圆环形区域的外圆半径Rout的像素数Pout。
S52:设定圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
S53:以中心点坐标(Centerx,Centery)为中心,将重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到圆环形图像对应的坐标点(Newx,Newy);
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol;
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow;
ScaleImgrow和ScaleImgcol分别为重采样全景矩形图像的行数和列数;
S54:将重采样全景矩形图像的坐标点(row,col)的像素值赋值到圆环形图像对应的坐标点(Newx,Newy),得到圆环形图像CircleImg。
本步骤中,利用重采样全景矩形图像的坐标点(row,col)的像素值ScaleImg(row,col),按照(row,col)与(Newx,Newy)的坐标对应关系,对圆环形图像的像素点按照如下公式进行逐像素赋值,得到尺度校准后的圆环形图像。
CircleImg(Newx,Newy)=ScaleImg(row,col)
本发明以圆环环宽内圆像素空间分辨率(即重采样全景矩形图像最内侧一列像素空间分辨率)为基准,通过极坐标与直角坐标转换,解决了由于高塔相机旋转拍摄引起的图像几何畸变问题,同时最大限度的保留了原始图像的采集信息。
作为本发明实施例的另一种改进,本发明的方法还包括:
S6:将圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
本步骤用于与高分辨率基准遥感影像进行地理配准,得到空间位置准确的几何校正后的高塔相机图像,其中的一个实现方式包括:
S61:以遥感影像为基准,将圆环形图像CircleImg的圆心的经纬度坐标设定为高塔相机中心点在地面的投影点的经纬度坐标。
S62:在圆环形图像和遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对。
高塔相机拍摄的全景矩形图像以及得到的圆环形图像上每个像素点对应着地面上的不同地物,有具体的地理位置;遥感影像上每个像元也对应着不同地物,有具体的地理位置;因此在圆环形图像和遥感影像上可以选取同名地物作为几何校正的控制点对。
S63:通过几何校正,进行圆环形图像投影变换和地理配准,使得圆环形图像和遥感影像的几何校正控制点对对准,得到空间位置准确的几何校正后的高塔相机图像。
下面以一个具体的试验示例对本发明进行详细阐述:
1、高塔相机360度拍摄的行数为22560行,列数为960列,去除无效区域后,得到有效行数为22380行,有效列数为960列的全景矩形图像,如图2所示。
2、装载在高塔上的高塔相机距离地面的垂直高度为65米,高塔相机镜头轴向相对于水平面的俯视角度α为23°,高塔相机镜头的视场角为21.6°。计算得到圆环形区域的外圆观测角为77.8°,内圆观测角为56.2°,外圆半径Rout为300.64米,内圆半径Rin为97.10米,圆环环宽Rring为203.54米。
4、计算重采样后的圆环形区域的圆环环宽Rring的像素数Rring为2467,圆环形区域的外圆半径Rout的像素数Pout为3644。
5、将全景矩形图像按照每一列像素的重采样比例系数重采样,获得重采样全景矩形图像。
6、将重采样全景矩形图像通过坐标转换,转换为圆环形图像。
7、以正射校正后的高分辨率卫星遥感影像为基准,对研究区域尺度校准后的圆环形图像进行地理配准和投影变换,获得几何校正后的高塔相机图像,如图3所示。
实施例2:
本发明实施例提供一种高塔相机图像几何校正装置,如图5所示,该装置包括:
全景矩形图像获取模块1,用于获取高塔相机360度拍摄得到的全景矩形图像;其中,全景矩形图像的行表示高塔相机360度拍摄的圆环形区域的径向,全景矩形图像的列表示高塔相机360度拍摄的圆环形区域的周向。
Rring=Rout-Rin
第二计算模块3,用于计算全景矩形图像每一列像素的重采样比例系数Scalei。
其中,Scalei为全景矩形图像第i列像素的重采样比例系数,i=1,2,…,OriImgcol,OriImgcol为全景矩形图像的列数,Pixceli为全景矩形图像第i列像素的采样尺度。
重采样模块4,用于将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像。
图像转换模块5,用于将重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。
本发明首先获取高塔相机360度拍摄得到的全景矩形图像,然后计算全景矩形图像每一列像素的重采样比例系数,并将全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,然后将重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。本发明在将全景矩形图像变换为圆环形真实图像照片的过程中,通过重新计算全景矩形图像每一列像素的重采样比例系数并进行重采样,解决了高塔相机成像尺度随视野距离增长呈非线性变化的问题,提高了高塔相机图像几何校正的准确性,高塔相机成像数据在人类活动监控、植被生长监测、大型动物跟踪监测等方面的应用研究奠定了基础。
作为本发明实施例的一种改进,前述的图像转换模块包括:
第一计算单元,用于计算重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Pring和圆环形区域的外圆半径Rout的像素数Pout。
中心点设定单元,用于设定圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
坐标变换单元,用于将重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到圆环形图像对应的坐标点(Newx,Newy)。
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol。
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow。
ScaleImgrow和ScaleImgcol分别为重采样全景矩形图像的行数和列数。
赋值单元,用于将重采样全景矩形图像的坐标点(row,col)的像素值赋值到圆环形图像对应的坐标点(Newx,Newy),得到圆环形图像。
前述的重采样模块进一步用于:
将全景矩形图像的行保持不变,按照每一列的重采样比例系数逐列进行重采样,并依次追加在前一列后面,直至所有列重采样完毕为止,得到重采样全景矩形图像。
进一步的,前述的全景矩形图像获取模块进一步用于:
获取高塔相机360度拍摄的图像,并去除无效区域,得到全景矩形图像;其中,无效区域包括开始拍摄和结束拍摄时的重复拍摄区域以及天际线周边及以外区域。
作为本发明实施例的一种改进,本发明的装置还包括:
地理配准模块,用于将圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
具体的,地理配准模块包括:
中心点匹配单元,用于以遥感影像为基准,将圆环形图像的圆心的经纬度坐标设定为高塔相机中心点在地面的投影点的经纬度坐标。
几何校正控制点对设置单元,用于在圆环形图像和遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对。
几何校正单元,用于通过几何校正,进行圆环形图像投影变换和地理配准,使得圆环形图像和遥感影像的几何校正控制点对对准,得到空间位置准确的高塔相机图像。
本发明实施例所提供的装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例1中相应内容。所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,前述描述的装置和单元的具体工作过程,均可以参考上述方法实施例中的对应过程,在此不再赘述。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明。本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种高塔相机图像几何校正方法,其特征在于,所述方法包括:
S1:获取高塔相机360度拍摄得到的全景矩形图像;其中,所述全景矩形图像的行表示所述高塔相机360度拍摄的圆环形区域的径向,所述全景矩形图像的列表示所述高塔相机360度拍摄的圆环形区域的周向;
Rring=Rout-Rin
S3:计算所述全景矩形图像每一列像素的重采样比例系数Scalei;
其中,Scalei为所述全景矩形图像第i列像素的重采样比例系数,i=1,2,…,OriImgcol,OriImgcol为所述全景矩形图像的列数,Pixceli为所述全景矩形图像第i列像素的采样尺度;
S4:将所述全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像;
S5:将所述重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。
2.根据权利要求1所述的高塔相机图像几何校正方法,其特征在于,所述S5包括:
S51:计算所述重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Pring和圆环形区域的外圆半径Rout的像素数Pout;
S52:设定所述圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
S53:将所述重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到所述圆环形图像对应的坐标点(Newx,Newy);
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol;
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow;
ScaleImgrow和ScaleImgcol分别为所述重采样全景矩形图像的行数和列数;
S54:将所述重采样全景矩形图像的坐标点(row,col)的像素值赋值到所述圆环形图像对应的坐标点(Newx,Newy),得到所述圆环形图像。
3.根据权利要求2所述的高塔相机图像几何校正方法,其特征在于,所述S4包括:
将所述全景矩形图像的行保持不变,按照每一列的重采样比例系数逐列进行重采样,并依次追加在前一列后面,直至所有列重采样完毕为止,得到所述重采样全景矩形图像。
4.根据权利要求3所述的高塔相机图像几何校正方法,其特征在于,所述S1包括:
获取高塔相机360度拍摄的图像,并去除无效区域,得到所述全景矩形图像;其中,所述无效区域包括开始拍摄和结束拍摄时的重复拍摄区域以及天际线周边及以外区域。
5.根据权利要求1-4任一所述的高塔相机图像几何校正方法,其特征在于,所述方法还包括:
S6:将所述圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
6.根据权利要求5所述的高塔相机图像几何校正方法,其特征在于,所述S6包括:
S61:以所述遥感影像为基准,将所述圆环形图像的圆心的经纬度坐标设定为所述高塔相机中心点在地面的投影点的经纬度坐标;
S62:在所述圆环形图像和所述遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对;
S63:通过几何校正,进行圆环形图像投影变换和地理配准,使得所述圆环形图像和所述遥感影像的几何校正控制点对对准,得到空间位置准确的高塔相机图像。
7.一种高塔相机图像几何校正装置,其特征在于,所述装置包括:
全景矩形图像获取模块,用于获取高塔相机360度拍摄得到的全景矩形图像;其中,所述全景矩形图像的行表示所述高塔相机360度拍摄的圆环形区域的径向,所述全景矩形图像的列表示所述高塔相机360度拍摄的圆环形区域的周向;
Rring=Rout-Rin
第二计算模块,用于计算所述全景矩形图像每一列像素的重采样比例系数Scalei;
其中,Scalei为所述全景矩形图像第i列像素的重采样比例系数,i=1,2,…,OriImgcol,OriImgcol为所述全景矩形图像的列数,Piexceli为所述全景矩形图像第i列像素的采样尺度;
重采样模块,用于将所述全景矩形图像的每一列分别按照该列的重采样比例系数进行重采样,得到重采样全景矩形图像;
图像转换模块,用于将所述重采样全景矩形图像通过坐标变换的方法转换为圆环形图像。
8.根据权利要求7所述的高塔相机图像几何校正装置,其特征在于,所述图像转换模块包括:
第一计算单元,用于计算所述重采样全景矩形图像对应的圆环形区域的圆环环宽Rring的像素数Pring和圆环形区域的外圆半径Rout的像素数Pout;
中心点设定单元,用于设定所述圆环形图像的中心点坐标(Centerx,Centery),其中:
Centerx≥Pout,Centery≥Pout
坐标变换单元,用于将所述重采样全景矩形图像的坐标点(row,col)进行逐像素坐标变换,得到所述圆环形图像对应的坐标点(Newx,Newy);
Newx=Centerx+C*sinθ-0.5
Newy=Centery-C*cosθ-0.5
其中,C=Pout-col,col=1,2,…,ScaleImgcol;
θ=(row+1)×2π/ScaleImgrow,row=1,2,…,ScaleImgrow;
ScaleImgrow和ScaleImgcol分别为所述重采样全景矩形图像的行数和列数;
赋值单元,用于将所述重采样全景矩形图像的坐标点(row,col)的像素值赋值到所述圆环形图像对应的坐标点(Newx,Newy),得到所述圆环形图像。
9.根据权利要求7或8所述的高塔相机图像几何校正装置,其特征在于,所述装置还包括:
地理配准模块,用于将所述圆环形图像与遥感影像进行地理配准,得到空间位置准确的高塔相机图像。
10.根据权利要求9所述的高塔相机图像几何校正装置,其特征在于,所述地理配准模块包括:
中心点匹配单元,用于以所述遥感影像为基准,将所述圆环形图像的圆心的经纬度坐标设定为所述高塔相机中心点在地面的投影点的经纬度坐标;
几何校正控制点对设置单元,用于在所述圆环形图像和所述遥感影像上均匀选取特征明显的同名地物点,作为几何校正控制点对;
几何校正单元,用于通过几何校正,进行圆环形图像投影变换和地理配准,使得所述圆环形图像和所述遥感影像的几何校正控制点对对准,得到空间位置准确的高塔相机图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211087790.2A CN115439365B (zh) | 2022-09-07 | 2022-09-07 | 高塔相机图像几何校正方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211087790.2A CN115439365B (zh) | 2022-09-07 | 2022-09-07 | 高塔相机图像几何校正方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115439365A true CN115439365A (zh) | 2022-12-06 |
CN115439365B CN115439365B (zh) | 2023-02-17 |
Family
ID=84246540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211087790.2A Active CN115439365B (zh) | 2022-09-07 | 2022-09-07 | 高塔相机图像几何校正方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115439365B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009069997A2 (en) * | 2007-11-27 | 2009-06-04 | Mimos Berhad | Method for image geometric transformation |
CN102905086A (zh) * | 2012-10-17 | 2013-01-30 | 中国人民解放军第四军医大学 | 一种用于全景拍摄的ccd 感光器 |
CN102982516A (zh) * | 2012-10-25 | 2013-03-20 | 西安理工大学 | 一种基于半球环形全景镜头实现全景图像的方法 |
CN104055534A (zh) * | 2013-03-21 | 2014-09-24 | 重庆海扶医疗科技股份有限公司 | 超声图像/超声视频传输方法和装置 |
JP2017053766A (ja) * | 2015-09-10 | 2017-03-16 | 株式会社豊田中央研究所 | 表面撮像装置、表面検査装置、及び表面撮像方法 |
WO2020137174A1 (ja) * | 2018-12-28 | 2020-07-02 | 株式会社Jvcケンウッド | プロジェクタシステム |
CN112330537A (zh) * | 2020-06-19 | 2021-02-05 | 中国科学院空天信息创新研究院 | 一种应急救援活动中无人机航拍影像快速拼接的方法 |
-
2022
- 2022-09-07 CN CN202211087790.2A patent/CN115439365B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009069997A2 (en) * | 2007-11-27 | 2009-06-04 | Mimos Berhad | Method for image geometric transformation |
CN102905086A (zh) * | 2012-10-17 | 2013-01-30 | 中国人民解放军第四军医大学 | 一种用于全景拍摄的ccd 感光器 |
CN102982516A (zh) * | 2012-10-25 | 2013-03-20 | 西安理工大学 | 一种基于半球环形全景镜头实现全景图像的方法 |
CN104055534A (zh) * | 2013-03-21 | 2014-09-24 | 重庆海扶医疗科技股份有限公司 | 超声图像/超声视频传输方法和装置 |
JP2017053766A (ja) * | 2015-09-10 | 2017-03-16 | 株式会社豊田中央研究所 | 表面撮像装置、表面検査装置、及び表面撮像方法 |
WO2020137174A1 (ja) * | 2018-12-28 | 2020-07-02 | 株式会社Jvcケンウッド | プロジェクタシステム |
CN112330537A (zh) * | 2020-06-19 | 2021-02-05 | 中国科学院空天信息创新研究院 | 一种应急救援活动中无人机航拍影像快速拼接的方法 |
Non-Patent Citations (5)
Title |
---|
LISA GOTTESFELD BROWN: "A survey of image registration techniques", 《ACM COMPUTING SURVEYS》 * |
LUCAS-ZHANG: "图像极坐标变换与反变换(Halcon + OpenCV)", 《CSDN》 * |
张淑敏: "全景环形透镜成像系统的研究与应用", 《西安电子科技大学》 * |
曹玲玲: "遥感图像快速预处理方法的研究与实现", 《中北大学》 * |
王权海等: "GDI+图像处理程序中的性能提高技术", 《物探化探计算技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115439365B (zh) | 2023-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110197466B (zh) | 一种广角鱼眼图像矫正方法 | |
CN108020826B (zh) | 多线激光雷达与多路相机混合标定方法 | |
CN106895851B (zh) | 一种光学遥感卫星多ccd多相机统一处理的传感器校正方法 | |
CN106123798B (zh) | 一种数字摄影激光扫描装置 | |
CN111486864B (zh) | 基于立体正八边结构的多源传感器联合标定方法 | |
CN113538595B (zh) | 利用激光测高数据辅助提升遥感立体影像几何精度的方法 | |
CN109242779A (zh) | 一种相机成像模型的构建方法、装置及汽车自动驾驶系统 | |
CN113393499B (zh) | 一种高分七号卫星全色影像和多光谱影像自动配准方法 | |
CN109341720A (zh) | 一种基于恒星轨迹的遥感相机几何标定方法 | |
CN104361563B (zh) | 基于gps的高光谱遥感图像几何精校正方法 | |
CN113793270A (zh) | 一种基于无人机姿态信息的航拍图像几何校正方法 | |
CN108734685A (zh) | 一种无人机载高光谱线阵列遥感影像的拼接方法 | |
CN108830811A (zh) | 一种飞行参数与相机内参相结合的航空影像实时校正方法 | |
CN111340888A (zh) | 一种无需白图像的光场相机检校方法及系统 | |
CN116310127A (zh) | 基于环形激光三角测量的管道内壁三维重建方法及系统 | |
CN111272196A (zh) | 特定拍摄条件下的在轨外方位元素自检校方法及系统 | |
CN111508028A (zh) | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 | |
CN111486868A (zh) | 基于地物特征的光电望远镜免方位标扩展标校法 | |
CN115439365B (zh) | 高塔相机图像几何校正方法和装置 | |
CN108917722A (zh) | 植被覆盖度计算方法及装置 | |
CN103776426A (zh) | 一种旋转平台农田成像几何校正方法 | |
CN110068313A (zh) | 一种基于投影变换的数字天顶仪定向方法 | |
CN105571598A (zh) | 一种卫星激光高度计足印相机姿态的测定方法 | |
CN107705272A (zh) | 一种空间影像的高精度几何校正方法 | |
CN111275773A (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 |