CN103398701B - 一种基于物方投影面的星载非共线tdi ccd影像拼接方法 - Google Patents
一种基于物方投影面的星载非共线tdi ccd影像拼接方法 Download PDFInfo
- Publication number
- CN103398701B CN103398701B CN201310326379.0A CN201310326379A CN103398701B CN 103398701 B CN103398701 B CN 103398701B CN 201310326379 A CN201310326379 A CN 201310326379A CN 103398701 B CN103398701 B CN 103398701B
- Authority
- CN
- China
- Prior art keywords
- image
- colinear
- tdi ccd
- object space
- projecting plane
- 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
Links
Landscapes
- Image Processing (AREA)
- Solid State Image Pick-Up Elements (AREA)
- Studio Devices (AREA)
Abstract
一种基于物方投影面的星载非共线时间延时积分电荷耦合元件(TDI CCD)影像拼接方法,包括以下步骤:(a).利用非共线TDI CCD相机获取非共线TDI CCD影像,建立所述非共线TDI CCD影像的传感器几何模型;(b).定义物方投影面;(c).根据所述的传感器几何模型和所定义的物方投影面建立拼接影像与所述非共线TDI CCD影像之间的像点坐标映射关系;(d).根据所述拼接影像与所述非共线TDI CCD影像之间的像点坐标映射关系对非共线TDI CCD影像进行重采样,生成重采样的拼接影像;(e).构建所述重采样的拼接影像的有理函数模型。本发明理论严密,具有良好的适应性,满足无缝拼接的要求。
Description
技术领域
本发明涉及测绘科学与技术领域,尤其是涉及一种基于物方投影面的星载非共线TDICCD影像拼接方法。
背景技术
时间延时积分电荷耦合元件(TDI CCD)能在不降低成像空间分辨率的同时,有效地提高空间相机的成像灵敏度和信噪比,因而正逐步取代常规CCD器件成为高分辨率光学推扫式相机的主流传感器。由于单片TDI CCD的像元个数有限,在空间相机的设计上,为了保证获取较高空间分辨率和一定幅宽的影像,需要将多片TDI CCD连接排列以获得较大的成像视场。然而,由于TDI CCD在物理结构上是一个小面阵,加上受器件外壳包装等物理因素的限制,多片TDI CCD在焦面上很难直接按照一条直线进行物理排列,而是通常采用非共线TDICCD焦面视场拼接或光学拼接的成像***设计方式。目前,IKONOS、QuickBird等高分辨率商业卫星以及我国资源三号、资源一号02B/02C、遥感二号等卫星上的光学相机均采用了非共线TDI CCD成像***设计。
非共线TDI CCD相机所获取的非共线TDI CCD影像(即原始影像)是按照每片TDICCD成像单独记录的,受传感器焦面摆放位置关系、偏流角控制、卫星姿态角、地面高程以及行积分时间等因素的影响,无法直接形成一幅完整的扫描景影像。
针对非共线TDI CCD影像拼接的现有技术,公开号为CN101799293A的中国专利申请公开了一种基于分段仿射变换的星载三片非共线TDI CCD影像拼接方法,公开号为CN101827223A的中国专利申请公开了一种基于行频归一化的非共线TDI CCD成像数据内视场拼接方法。这两种非共线TDI CCD影像拼接方法都是基于片间匹配的像面平移拼接。
然而,这些传统的基于片间匹配的像面平移拼接方法不仅在理论上缺乏严密性,而且也没能充分考虑非共线TDI CCD相机的成像特点,因而在积分时间跳变、地形起伏剧烈和大侧摆等复杂成像条件下并不适用。多片TDI CCD影像的高精度和高可靠拼接是保证成像数据后期处理和应用质量的基础,因此,围绕该类相机成像数据的严格拼接理论和方法展开研究,具有重要的意义和应用价值。
发明内容
针对现有技术存在的不足,本发明提出一种基于物方投影面的星载非共线TDI CCD影像拼接方法。
本发明提出的技术方案如下:
一种基于物方投影面的星载非共线时间延时积分电荷耦合元件(TDI CCD)影像拼接方法,包括以下步骤:
(a).利用非共线TDI CCD相机获取非共线TDI CCD影像,建立所述非共线TDI CCD影像的传感器几何模型;
(b).定义物方投影面;
(c).根据所述的传感器几何模型和所定义的物方投影面建立拼接影像与所述非共线TDICCD影像之间的像点坐标映射关系;
(d).根据所述的像点坐标映射关系对非共线TDI CCD影像进行重采样,生成重采样的拼接影像;
(e)构建所述重采样的拼接影像的有理函数模型。
本发明基于物方空间地物分布的连续性,沿卫星推扫成像地面目标轨迹方向,对物方投影面上的成像覆盖区域划分等间隔格网,格网大小与非共线TDI CCD影像的地面分辨率相一致,首先,将拼接影像的像点与物方投影面格网点一一对应,然后,基于非共线TDI CCD影像的传感器几何模型便可建立起拼接影像与非共线TDI CCD影像的像点坐标映射关系,最后,采用间接法影像纠正的方式对非共线TDI CCD影像进行重采样生成拼接影像。本发明理论严密,具有良好的适应性。拼接影像的生成过程等效于沿卫星推扫成像的地面目标轨迹方向对非共线TDI CCD影像进行了一次准***几何校正。当非共线TDI CCD片间几何定位精度一致性时,拼接影像上每个像点灰度信息具有唯一可确定性,满足无缝拼接的要求。
附图说明
图1为本发明的基于物方投影面的星载非共线TDI CCD影像拼接方法的原理示意图;
图2为本发明的基于物方投影面的星载非共线TDI CCD影像拼接方法的流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步的说明。
本发明利用非共线TDI CCD的片间摄影几何约束关系,基于物方空间地物分布的连续性,建立拼接影像与非共线TDI CCD影像(即原始影像)的像点坐标换算关系,进而实现对非共线TDI CCD影像的无缝拼接处理。其技术内涵主要体现在两个方面:(1)以某种形式建立拼接影像像点与物方空间的坐标映射关系;(2)基于非共线TDI CCD的成像几何模型建立物方空间与非共线TDI CCD影像像点的坐标映射关系。于是,本发明提出了一种基于物方投影面的星载非共线TDI CCD影像拼接方法。如图1所示,拼接影像上的像点与所建立的物方投影面等间隔格网点一一对应;基于非共线TDI CCD影像的传感器几何模型,可反投影计算物方投影面任一投影点所对应的分片TDI CCD影像像点坐标,这样便建立了拼接影像与非共线TDI CCD影像的像点坐标对应关系。
如附图2所示,本方法的实施可以分为以下几个步骤:
(a).利用非共线TDI CCD相机获取非共线TDI CCD影像,建立所述非共线TDI CCD影像的传感器几何模型:
对于星载线阵推扫式成像数据,其严格成像几何模型由公式(1-1)表达:
其中,α,β分别为成像探元的光线相对于相机焦平面横轴和纵轴方向的指向角;m为摄影比例尺因子; 为地面点的WGS84地心直角坐标; 为成像时刻卫星在WGS84地心直角坐标系下的坐标;RBS为是相机坐标系与卫星本体坐标系之间的旋转矩阵;RFB为卫星本体坐标系到轨道坐标系之间的姿态矩阵;RGF为轨道坐标系与地球惯性坐标系CIS之间的坐标转换矩阵;RT为像点成像时刻地球惯性坐标系CIS到地球固定坐标系CTS的坐标转换矩阵。
由相机的外方位元素和每片TDI CCD的探元指向角,基于公式(1-1)建立每片TDICCD影像的严格成像几何模型,所有分片TDI CCD影像的严格成像几何模型共同构成非共线TDI CCD影像(即原始影像)的传感器几何模型。
(b).定义物方投影面:
将物方投影面定义在成像区域所对应的通用横轴墨卡托(UTM)投影坐标系下,投影面的高程取WGS84大地坐标系下的平均椭球高H。
(c).根据所述的传感器几何模型和所定义的物方投影面建立拼接影像与所述非共线TDI CCD影像之间的像点坐标映射关系:
(c1)基于非共线TDI CCD影像的传感器几何模型和物方投影面的高程值H,分别计算非共线TDI CCD影像第一列的首、尾像素在物方投影面上投影点坐标,将两个投影点的连线方向定义为非共线TDI CCD推扫成像时地面目标轨迹在物方投影面上的近似前进方向,即拼接影像的每一列在物方投影面上的排列方向;
(c2)基于非共线TDI CCD影像的传感器几何模型和物方投影面的高程值H,计算非共线TDI CCD影像在物方投影面上的投影覆盖区域AOP(Area ofProiection);
(c3)确定拼接影像在物方投影面上的覆盖区域AOI(Area of Interest):取AOP在拼接影像列排列方向上的最小外廓矩形;
(c4)令拼接影像上每个像素所对应的物方投影面采样间隔大小相同,且等于非共线TDI CCD影像的地面采样间隔GSD(Ground Sampling Distance),进而确定拼接影像的行数和列数;
(c5)建立拼接影像像点与物方投影面投影点之间的二维坐标转换关系,如公式(1-2):
式中,
G为物方投影面采样间隔;
是根据拼接影像在物方投影面上的列排列方向确定的二维坐标旋转矩阵;
(c6)基于非共线TDI CCD影像的传感器几何模型,确定物方投影面投影点 对应的非共线TDI CCD影像像点的像素坐标 至此便建立了拼接影像与非共线TDI CCD影像的像点坐标映射关系;
(d).根据所述的像点坐标映射关系对非共线TDI CCD影像进行重采样,生成重采样的拼接影像:
根据拼接影像与非共线TDI CCD影像之间的像点坐标映射关系,采用间接法影像纠
正的方式对非共线TDI CCD影像进行重采样,生成重采样的拼接影像;
(e).构建所述重采样的拼接影像的有理函数模型:
(e1)由所述重采样的拼接影像上的任意一点影像坐标 根据所述公式(1-2),得到其对应的物方投影面投影点 坐标;
(e2)基于非共线TDI CCD影像的传感器几何模型,确定物方投影面投影点 对应的非共线TDI CCD影像像点的像素坐标
(e3)通过非共线TDI CCD影像上的影像坐标 加上高程值H,利用非共线TDI CCD影像的传感器几何模型正变换公式生成对应的大地坐标
(e4)根据步骤(e1)、(e2)、(e3)所建立的流程得到拼接影像严密传感器几何模型的正算公式 其中,T表示由拼接影像像点坐标正算到大地经纬度坐标的转换关系;根据该正算公式构建像方均匀分布的控制点立体格网,获得所述重采样的拼接影像的有理函数模型。
本发明基于物方空间地物分布的连续性,沿卫星推扫成像地面目标轨迹方向,对物方投影面上的成像覆盖区域划分等间隔格网,格网大小与非共线TDI CCD影像(即原始影像)的地面分辨率相一致,首先,将拼接影像的像点与物方投影面格网点一一对应,然后,基于非共线TDI CCD影像的传感器几何模型便可建立起拼接影像与非共线TDI CCD影像的像点坐标映射关系,最后,采用间接法影像纠正的方式对非共线TDI CCD影像进行重采样生成拼接影像。本发明理论严密,具有良好的适应性。拼接影像的生成过程等效于沿卫星推扫成像的地面目标轨迹方向对非共线TDI CCD影像进行了一次准***几何校正。当非共线TDI CCD片间几何定位精度一致性时,拼接影像上每个像点灰度信息具有唯一可确定性,满足无缝拼接的要求。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案而非对其限制;尽管参照较佳实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解,依然可对本发明的具体实施方式进行修改或者对部分技术特征进行等同替换而不脱离本发明技术方案的精神,其均应当涵盖本发明请求保护的技术方案范围当中。
Claims (3)
1.一种基于物方投影面的星载非共线时间延时积分电荷耦合元件(TDI CCD)影像拼接方法,包括以下步骤:
(a).利用非共线TDI CCD相机获取非共线TDI CCD影像,建立所述非共线TDI CCD影像的传感器几何模型,其中,所述步骤(a)包括:对于星载线阵推扫式成像数据,其严格成像几何模型由公式(1-1)表达,
其中,α,β分别为成像探元的光线相对于相机焦平面横轴和纵轴方向的指向角;m为摄影比例尺因子;为地面点的WGS84地心直角坐标;为成像时刻卫星在WGS84地心直角坐标系下的坐标;RBS为相机坐标系与卫星本体坐标系之间的旋转矩阵;RFB为卫星本体坐标系到轨道坐标系之间的姿态矩阵;RGF为轨道坐标系与地球惯性坐标系CIS之间的坐标转换矩阵;RT为像点成像时刻地球惯性坐标系CIS到地球固定坐标系CTS的坐标转换矩阵
由相机的外方位元素和每片TDI CCD的探元指向角,基于公式(1-1)建立每片TDI CCD影像的严格成像几何模型,所有分片TDI CCD影像的严格成像几何模型共同构成非共线TDI CCD影像的传感器几何模型;
(b).定义物方投影面;
(c).根据所述的传感器几何模型和所定义的物方投影面建立拼接影像与所述非共线TDI CCD影像之间的像点坐标映射关系,其中,所述步骤(c)包括:
(c1)基于非共线TDI CCD影像的传感器几何模型和物方投影面的高程值H,分别计算非共线TDI CCD影像第一列的首、尾像素在物方投影面上投影点坐标,将两个投影点的连线方向定义为非共线TDI CCD推扫成像时地面目标轨迹在物方投影面上的近似前进方向,即拼接影像的每一列在物方投影面上的排列方向;
(c2)基于非共线TDI CCD影像的传感器几何模型和物方投影面的高程值H,计算非共线TDI CCD影像在物方投影面上的投影覆盖区域(AOP);
(c3)确定拼接影像在物方投影面上的覆盖区域(A01):取AOP在拼接影像列排列方向上的最小外廓矩形;
(c4)令拼接影像上每个像素所对应的物方投影面采样间隔大小相同,且等于非共线TDI CCD 影像的地面采样间隔(GSD),进而确定拼接影像的行数和列数;
(c5)建立拼接影像像点与物方投影面投影点之间的二维坐标转换关系,如公式(1-2):
式中,
为拼接影像像点的像素坐标;
为拼接影像像点在物方投影面上投影点的平面坐标;
G为物方投影面采样间隔;
是根据拼接影像在物方投影面上的列排列方向确定的二维坐标旋转矩阵;为坐标平移量;
(c6)基于非共线TDI CCD影像的传感器几何模型,确定物方投影面投影点对应的非共线TDI CCD影像像点的像素坐标至此便建立了拼接影像与非共线TDI CCD影像的像点坐标映射关系;
(d).根据所述的像点坐标映射关系对非共线TDI CCD影像进行重采样,生成重采样的拼接影像;
(e).构建所述重采样的拼接影像的有理函数模型。
2.如权利要求1所述的基于物方投影面的星载非共线时间延时积分电荷耦合元件(TDI CCD)影像拼接方法,所述步骤(b)包括:将物方投影面定义在成像区域所对应的通用横轴墨卡托(UTM)投影坐标系下,投影面的高程值取WGS84大地坐标系下的平均椭球高H。
3.如权利要求1所述的基于物方投影面的星载非共线时间延时积分电荷耦合元件(TDI CCD)影像拼接方法,所述步骤(d)包括:根据拼接影像与非共线TDI CCD影像之间的像点坐标映射关系,采用间接法影像纠正的方式对非共线TDI CCD影像进行重采样,生成重采样的拼接影像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310326379.0A CN103398701B (zh) | 2013-07-31 | 2013-07-31 | 一种基于物方投影面的星载非共线tdi ccd影像拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310326379.0A CN103398701B (zh) | 2013-07-31 | 2013-07-31 | 一种基于物方投影面的星载非共线tdi ccd影像拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103398701A CN103398701A (zh) | 2013-11-20 |
CN103398701B true CN103398701B (zh) | 2015-07-08 |
Family
ID=49562373
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310326379.0A Active CN103398701B (zh) | 2013-07-31 | 2013-07-31 | 一种基于物方投影面的星载非共线tdi ccd影像拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103398701B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679673B (zh) * | 2013-11-22 | 2016-06-08 | 中国资源卫星应用中心 | 一种宽视场线阵ccd影像几何畸变模拟方法 |
CN103697864B (zh) * | 2013-12-27 | 2015-11-04 | 武汉大学 | 一种基于大虚拟相机的窄视场双相机影像拼接方法 |
CN103778610B (zh) * | 2014-01-24 | 2016-09-14 | 武汉大学 | 一种星载线阵传感器垂轨摆扫影像的几何预处理方法 |
WO2018117985A1 (en) * | 2016-12-21 | 2018-06-28 | Ananbanchachai Nenin | Projector system for spherical platform |
CN110111274B (zh) * | 2019-04-28 | 2020-06-19 | 张过 | 一种星载推扫式光学传感器外方位元素定标方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6782334B1 (en) * | 2003-04-01 | 2004-08-24 | Lockheed Martin Corporation | Method and system for calibration of time delay integration imaging devices |
CN101799293A (zh) * | 2010-03-05 | 2010-08-11 | 武汉大学 | 基于分段仿射变换的星载三片非共线tdiccd影像拼接方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0576005A (ja) * | 1991-09-13 | 1993-03-26 | Nikon Corp | 画像入力装置 |
-
2013
- 2013-07-31 CN CN201310326379.0A patent/CN103398701B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6782334B1 (en) * | 2003-04-01 | 2004-08-24 | Lockheed Martin Corporation | Method and system for calibration of time delay integration imaging devices |
CN101799293A (zh) * | 2010-03-05 | 2010-08-11 | 武汉大学 | 基于分段仿射变换的星载三片非共线tdiccd影像拼接方法 |
Non-Patent Citations (1)
Title |
---|
An Algorithm for Mosaicking Non一collinear TDI CCD Chip Images based on Reference Plane in Object-space;Fen Hu et.al;《第六届国际数字地球会议论文集》;国际数字地球学会第六届国际数字地球会议论文集;20100426;正文第1页第1节第4行、第3页第10-11行、第3页第25-32行、第7页第3-5行 * |
Also Published As
Publication number | Publication date |
---|---|
CN103398701A (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110310248B (zh) | 一种无人机遥感影像实时拼接方法及*** | |
CN103398701B (zh) | 一种基于物方投影面的星载非共线tdi ccd影像拼接方法 | |
CN103914808B (zh) | 一种资源三号卫星三线阵影像和多光谱影像的拼接方法 | |
US9451236B2 (en) | Apparatus for synthesizing three-dimensional images to visualize surroundings of vehicle and method thereof | |
CN101604018B (zh) | 高分辨率遥感影像数据处理方法及其*** | |
CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
CN109115186B (zh) | 一种针对车载移动测量***的360°可量测全景影像生成方法 | |
CN107808362A (zh) | 一种基于无人机pos信息与图像surf特征结合的图像拼接方法 | |
CN104463778B (zh) | 一种全景图生成方法 | |
WO2020062434A1 (zh) | 一种相机外部参数的静态标定方法 | |
CN103900539B (zh) | 一种空中立方体全景成像目标定位方法 | |
CN102667619B (zh) | 具有多个光学***和检测器阵列的大格式数字摄影机 | |
CN105716590A (zh) | 确定大地测量装置的位置和定向偏移的方法及测量装置 | |
CN105678783A (zh) | 折反射全景相机与激光雷达数据融合标定方法 | |
CN105513018A (zh) | 星载摆扫成像几何校正方法和装置 | |
CN105096252B (zh) | 一种带状全方位街景影像图的制作方法 | |
CN102628942A (zh) | 一种雷达影像双视向信息补偿方法 | |
CN108830811A (zh) | 一种飞行参数与相机内参相结合的航空影像实时校正方法 | |
CN103226840A (zh) | 全景影像拼接及量测***及方法 | |
CN115439531A (zh) | 一种获取目标对象的目标空间位置信息的方法与设备 | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
CN104655106A (zh) | 基于gps rtk和全景影像的自主定位定向测图方法 | |
Liu et al. | A new approach to fast mosaic UAV images | |
JP2000276045A (ja) | 撮影画像を用いた地図作成方法、装置および撮影画像の歪み補正方法 | |
CN103778610A (zh) | 一种星载线阵传感器垂轨摆扫影像的几何预处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100830 No. 28 Lianhuachi West Road, Haidian District, Beijing Patentee after: Ministry of Natural Resources Land Satellite Remote Sensing Application Center Address before: 100830 No. 28 Lianhuachi West Road, Haidian District, Beijing Patentee before: Satellite Surveying and Mapping Application Center, NASG |
|
CP01 | Change in the name or title of a patent holder |