• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 13
    • 下载费用:30 金币  

    重庆时时彩是否合法: 一种资源三号卫星三线阵影像和多光谱影像的拼接方法.pdf

    关 键 词:
    一种 资源 三号 卫星 三线 影像 光谱 拼接 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410095059.3

    申请日:

    2014.03.14

    公开号:

    CN103914808A

    公开日:

    2014.07.09

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G06T 3/40申请日:20140314|||公开
    IPC分类号: G06T3/40 主分类号: G06T3/40
    申请人: 国家测绘地理信息局卫星测绘应用中心
    发明人: 张过; 唐新明; 蒋永华; 高小明; 祝小勇; 胡芬; 方辰; 邓明军
    地址: 101300 北京市顺义区机场东路国家地理信息科技产业园
    优先权:
    专利代理机构: 北京纪凯知识产权代理有限公司 11245 代理人: 徐宁
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410095059.3

    授权公告号:

    ||||||

    法律状态公告日:

    2017.02.01|||2014.08.06|||2014.07.09

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明涉及一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其包括以下步骤:建立真实CCD线阵的严密成像几何模型;根据真实CCD线阵的严密成像几何模型,建立虚拟CCD线阵的严密成像几何模型;建立拼接影像的有理多项式模型;完成资源三号卫星三线阵影像和多光谱影像的拼接。利用基于虚拟CCD线阵的内视场的拼接方案,在焦平面上构造出无畸变的CCD阵列,对原始多个CCD成像影像按照线中心投影成像方式进行重成像,以实现多个CCD成像影像的线中心投影无缝拼接。本发明可以广泛应用于资源三号卫星影像处理中。

    权利要求书

    权利要求书
    1.  一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其具体包括以下步骤:
    1)建立真实CCD线阵的严密成像几何模型;
    2)根据真实CCD线阵的严密成像几何模型,建立虚拟CCD线阵的严密成像几何模型;
    3)建立拼接影像的有理多项式模型;
    4)完成资源三号卫星三线阵影像和多光谱影像的拼接。

    2.  如权利要求1所述的一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其特征在于:所述步骤1)中,采用以下步骤建立真实CCD线阵的严密成像几何模型:
    ①利用全色或多光谱相机的像主点(x0,y0),像主距f,像元(x,y)大小,确定像元(x,y)在全色或多光谱相机坐标系下的坐标为:
    x-x0y-y0-f;]]>
    ②根据像元(x,y)的行y,采用线性内插方法内插得到成像时间t,通过成像时间t,分别采用拉格朗日内插和线性内插方法得到某成像时刻的轨道XsYsZsWGS84]]>和姿态Q;
    ③根据由步骤②得到的轨道XsYsZsWGS84]]>和姿态Q,建立像点坐标与地面点坐标的对应关系,即真实CCD线阵的严密成像几何模型为:
    XYZWGS84=XsYsZsWGS84+mRbody2WGS84Rcamera2bodyx-x0y-y0-f;]]>
    式中,XYZWGS84]]>是影像拍摄地面点的地面坐标;m为成像比例尺;XsYsZsWGS84]]>是某成像时刻的轨道;Rbody2WGS84是由卫星下传的定姿设备测量获得的描述卫星成像姿态的数据;Rcamera2body是载荷在卫星本体坐标系下的安装矩阵,描述载荷与卫星本体的相对关系;
    根据真实CCD线阵的严密成像几何模型以及影像的长度和宽度,计算得到真实CCD线阵对应的地面影像范围。

    3.  如权利要求1所述的一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其特征在于:所述步骤2)中,采用以下步骤建立虚拟CCD线阵的严密成像几何模型:
    ①确定虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标为:
    x=(xmax-xmin)/wvirtual·s+xminy=ax+b-f,]]>
    式中,wvirtual为虚拟CCD影像的宽度,a和b均为虚拟CCD线阵直线方程的系数,f为像主距,xmax=max(xi),xmin=min(xi),i≤n;
    ②建立t时刻卫星姿态模型为:

    式中,ωn,κn为卫星不同时刻姿态拟合的多项式系数,n为采用的多项式模型阶数,n=0,1,…;
    ③将步骤①确定的虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标以及步骤②建立的t时刻卫星姿态模型均代入真实CCD线阵的严密成像几何模型中,得到虚拟CCD线阵的严密成像几何模型。

    4.  如权利要求1或2或3所述的一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其特征在于:所述步骤3)中,建立拼接影像的有理多项式模型的过程包括:
    ①根据据虚拟CCD线阵的严密成像几何模型,分别计算虚拟CCD影像的上下左右四个角点的地理坐标(lon,lat,h),即某区域的地理范围,从全球1km格网SRTM数据上读取该区域地理范围的影像,计算该区域的高程范围;
    ②在虚拟CCD的影像空间内,根据虚拟影像大小划分虚拟控制格网,根据虚拟CCD线阵的严密成像几何模型计算格网上各个网格点的地理坐标,组成控制点(X,Y,P,L,H),并求解控制点(X,Y,P,L,H)对应的RFM模型:
    Y=NL(P,L,H)DL(P,L,H)X=NS(P,L,H)DS(P,L,H),]]>
    式中,X,Y为影像坐标,P,L,H分别为地面坐标纬度、经度和高程;NL(P,L,H)为:
    NL(P,L,H)=a1+a2L+a3P+a4H+a5LP+a6LH+a7PH+a8L2+a9P2
              +a10H2+a11PLH+a12L3+a13LP2+a14LH2+a15L2P+a16P3+a17PH2,
              +a18L2H+a19P2H+a20H3
    将NL(P,L,H)中的参数ai换成bi,ci和di,分别得到DL(P,L,H),Ns(P,L,H)和Ds(P,L,H),i=1,2,…,20。

    5.  如权利要求1或2或3所述的一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其特征在于:所述步骤4)中,采用以下步骤完成资源三号卫星三线阵影像和多光谱影像的拼接:
    ①确定虚拟CCD影像的高度和宽度,其中,将原始影像的高度设定为虚拟CCD影像的高度,虚拟CCD影像的宽度为wvirtual;
    ②对控制点(X,Y,P,L,H)对应的RFM模型进行逆变换,并计算得到拼接后影像上任一点O对应高程h的地面点坐标(lat,lon,h);
    ③利用反算模型,将地面点坐标(lat,lon,h)投影到真实CCD影像上,得到地面点坐标(lat,lon,h)对应的像点坐标(x',y');
    ④采用距离加权方法,计算得到像点坐标(x',y')对应的灰度值g(x',y')为:
    g(x',y')=(1-dx)·(1-dy)·g1+dx·(1-dy)·g2+dy·(1-dx)·g3+dx·dy·g4;
    ⑤对已知高度和宽度的虚拟CCD影像上的所有像素,重复步骤②~步骤④,完成资源三号卫星三线阵影像和多光谱影像的拼接。

    6.  如权利要求5所述的一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其特征在于:所述步骤③中,计算得到地面点坐标(lat,lon,h)对应的像点坐标(x',y'),其具体包括以下步骤:
    Ⅰ、根据真实CCD线阵的严密成像几何模型,计算得到原始影像上下左右四个角点的地理坐标(lat,lon),根据像元(x,y)与地理坐标(lat,lon)之间的关系式:
    x=f0+f1lat+f2lony=g0+g1lat+g2lon,]]>
    计算仿射变换参数f0、f1、f2、g0、g1、g2;
    Ⅱ、将计算得到的地面点坐标(lat,lon,h)代入像元(x,y)与地理坐标(lat,lon)之间的关系式,得到地面点坐标(lat,lon,h)对应的影像点的平面坐标(xp,yp);
    Ⅲ、将影像点的平面坐标(xp,yp)代入真实CCD线阵的严密成像几何模型,计算得到影像点坐标(xp,yp,h)的地面点坐标(latp,lonp);
    Ⅳ、重复步骤Ⅰ~步骤Ⅲ,直到满足迭代条件:

    地面点坐标(latp,lonp)所对应的影像点的平面坐标(xp,yp)即为地面点对应的像点坐标(x',y')。

    说明书

    说明书一种资源三号卫星三线阵影像和多光谱影像的拼接方法
    技术领域
    本发明涉及一种影像拼接方法,特别是关于一种资源三号卫星三线阵影像和多光谱影像的拼接方法。
    背景技术
    在高分辨光学线阵推扫式卫星相机的设计中,为获取较大的幅宽,将多个CCD(Charge-coupled Device,电荷耦合元件)线阵在焦面上交错安装,通过分光等模式形成近似直线的一条连续CCD阵列。但是由于多个CCD线阵的安装精度和TDI-CCD(Time Delayed and IntegrationCharge-coupled Device时间延迟和积分电荷耦合元件)沿着轨道方向宽度等因素的限制,多个CCD线阵在焦面不可能形成一条严格意义上的CCD线阵。在摄影测量数据生产中,均使用卫星影像供应商提供的按照成像几何特性和辐射特性拼接后的影像,例如IKONOS伊科诺斯、QuickBird快鸟和WorldView-2等。而影像的拼接精度会直接影响后续摄影测量生产的精度,如在利用IKONOS立体像对匹配生产DEM(Digital Elevation Model,数字高程模型)时,拼接精度不高会导致DEM精度有损失。多CCD线阵的亚像素级拼接精度是保证DEM生产精度的前提。如何对交错安装的多个CCD线阵影像进行高精度的拼接成为了高分辨率遥感卫星应用中亟待解决的问题。
    多光谱相机一般采用多谱合一的TDI-CCD,那么同一时刻各个谱段对应不同的摄影点。即各个波段对应同一地物的成像时间不一致,而卫星在飞行过程中存在高频颤振,那么各个波段对应同一地物的成像姿态也不一致,若是波段直接相加,那么各个波段的影像无法保证完全叠合。而且,多光谱影像上也存在着多个CCD影像拼接的问题,如何对多光谱数据进行拼接,同时保证其无缝是有效利用多光谱影像的关键。
    而当前多CCD影像拼接的方法主要可以分为两类:基于像方的影像拼接和基于物方的影像拼接?;谙穹降挠跋衿唇尤衔跋裰浠旧现淮嬖谄苹蛘叩徒谆?,从相邻CCD影像的重叠区域中获取一些同名点,再利用同名点之间的位置关系,建立起平移甚至仿射变换等低阶畸变模型,从而建立起影像与影像间的对应关系,完成拼接。但是,该拼接方案必须基于一个假设前提,那就是卫星飞行平台稳定,姿态不存在较大跳变。然而,事实上卫星在飞行过程中,为了保持三轴稳定,姿轨控制系统会不断调整卫星的姿态,由于这些因素的存在,使得卫星姿态中存在高频抖动,从而导致不位于同一直线的CCD之间会出现成像时间不一致的现象,该现象称为卫星姿态的高频噪声,该高频噪声是无法通过简单的低阶模型消除的。而且,基于像方的影像拼 接严重依赖于同名点匹配的精度,一旦拼接处纹理单一,同名点匹配可靠性不足时,就会出现拼接错误。除此之外,更为严重的问题是,该拼接方法生产的影像丢失了严密成像几何模型。如果直接利用通用传感器模型,如RFM(Rational Function Model,有理多项式模型)来替代原始模型的话,又会存在模型替代精度不高的问题,如ALOS(日本的对地观测卫星)的1B产品。与之不同,基于物方的影像拼接方法采用基于各CCD线阵的严密成像几何模型,根据各CCD的物方定位建立相邻CCD线阵关系,完成拼接;因此,采用该方案建立的相邻CCD位置关系更为严密,且其拼接精度不依赖于同名点匹配精度。
    发明内容
    针对上述问题,本发明的目的是提供一种不依赖于同名点匹配精度的资源三号卫星三线阵影像和多光谱影像的拼接方法,该方法是一种物方拼接方法,能够解决常用像方拼接方法的不足,并通过虚拟CCD的严密成像几何模型,建立拼接影像的严密成像几何模型,在完成影像拼接的同时,实现通用几何模型的建立,使得拼接影像能够进一步应用于摄影测量等处理。
    为实现上述目的,本发明采取以下技术方案:一种资源三号卫星三线阵影像和多光谱影像的拼接方法,其具体包括以下步骤:1)建立真实CCD线阵的严密成像几何模型;2)根据真实CCD线阵的严密成像几何模型,建立虚拟CCD线阵的严密成像几何模型;3)建立拼接影像的有理多项式模型;4)完成资源三号卫星三线阵影像和多光谱影像的拼接。
    所述步骤1)中,采用以下步骤建立真实CCD线阵的严密成像几何模型:①利用全色或多光谱相机的像主点(x0,y0),像主距f,像元(x,y)大小,确定像元(x,y)在全色或多光谱相机坐标系下的坐标为:
    x-x0y-y0-f;]]>
    ②根据像元(x,y)的行y,采用线性内插方法内插得到成像时间t,通过成像时间t,分别采用拉格朗日内插和线性内插方法得到某成像时刻的轨道XsYsZsWGS84]]>和姿态Q;
    ③根据由步骤②得到的轨道XsYsZsWGS84]]>和姿态Q,建立像点坐标与地面点坐标的对应关 系,即真实CCD线阵的严密成像几何模型为:
    XYZWGS84=XsYsZsWGS84+mRbody2WGS84Rcamera2bodyx-x0y-y0-f;]]>
    式中,XYZWGS84]]>是影像拍摄地面点的地面坐标;m为成像比例尺;XsYsZsWGS84]]>是某成像时刻的轨道;Rbody2WGS84是由卫星下传的定姿设备测量获得的描述卫星成像姿态的数据;Rcamera2body是载荷在卫星本体坐标系下的安装矩阵,描述载荷与卫星本体的相对关系;根据真实CCD线阵的严密成像几何模型以及影像的长度和宽度,计算得到真实CCD线阵对应的地面影像范围。
    所述步骤2)中,采用以下步骤建立虚拟CCD线阵的严密成像几何模型:①确定虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标为:
    x=(xmax-xmin)/wvirtual·s+xminy=ax+b-f,]]>
    式中,wvirtual为虚拟CCD影像的宽度,a和b均为虚拟CCD线阵直线方程的系数,f为像主距,xmax=max(xi),xmin=min(xi),i≤n;②建立t时刻卫星姿态模型为:

    式中,ωn,κn为卫星不同时刻姿态拟合的多项式系数,n为采用的多项式模型阶数,n=0,1,…;③将步骤①确定的虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标以及步骤②建立的t时刻卫星姿态模型均代入真实CCD线阵的严密成像几何模型中,得到虚拟CCD线阵的严密成像几何模型。
    所述步骤3)中,建立拼接影像的有理多项式模型的过程包括:①根据据虚拟CCD线阵的严密成像几何模型,分别计算虚拟CCD影像的上下左右四个角点的地理坐标(lon,lat,h),即某区域的地理范围,从全球1km格网SRTM数据上读取该区域地理范围的影像,计算该区域的高程范围;②在虚拟CCD的影像空间内,根据虚拟影像大小划分虚拟控制格网,根据虚拟CCD线阵的严密成像几何模型计算格网上各个网格点的地理坐标,组成控制点(X,Y,P,L,H),并求解控制点(X,Y,P,L,H)对应的RFM模型:
    Y=NL(P,L,H)DL(P,L,H)X=NS(P,L,H)DS(P,L,H),]]>
    式中,X,Y为影像坐标,P,L,H分别为地面坐标纬度、经度和高程;NL(P,L,H)为:
    NL(P,L,H)=a1+a2L+a3P+a4H+a5LP+a6LH+a7PH+a8L2+a9P2
              +a10H2+a11PLH+a12L3+a13LP2+a14LH2+a15L2P+a16P3+a17PH2,
              +a18L2H+a19P2H+a20H3
    将NL(P,L,H)中的参数ai换成bi,ci和di,分别得到DL(P,L,H),Ns(P,L,H)和Ds(P,L,H),i=1,2,…,20。
    所述步骤4)中,采用以下步骤完成资源三号卫星三线阵影像和多光谱影像的拼接:①确定虚拟CCD影像的高度和宽度,其中,将原始影像的高度设定为虚拟CCD影像的高度,虚拟CCD影像的宽度为wvirtual;②对控制点(X,Y,P,L,H)对应的RFM模型进行逆变换,并计算得到拼接后影像上任一点O对应高程h的地面点坐标(lat,lon,h);③利用反算模型,将地面点坐标(lat,lon,h)投影到真实CCD影像上,得到地面点坐标(lat,lon,h)对应的像点坐标(x',y');④采用距离加权方法,计算得到像点坐标(x',y')对应的灰度值g(x',y')为:
    g(x',y')=(1-dx)·(1-dy)·g1+dx·(1-dy)·g2+dy·(1-dx)·g3+dx·dy·g4;
    ⑤对已知高度和宽度的虚拟CCD影像上的所有像素,重复步骤②~步骤④,完成资源三号卫星三线阵影像和多光谱影像的拼接。
    所述步骤③中,计算得到地面点坐标(lat,lon,h)对应的像点坐标(x',y'),其具体包括以下步骤:Ⅰ、根据真实CCD线阵的严密成像几何模型,计算得到原始影像上下左右四个角点的地理坐标(lat,lon),根据像元(x,y)与地理坐标(lat,lon)之间的关系式:
    x=f0+f1lat+f2lony=g0+g1lat+g2lon,]]>
    计算仿射变换参数f0、f1、f2、g0、g1、g2;Ⅱ、将计算得到的地面点坐标(lat,lon,h)代入像元(x,y)与地理坐标(lat,lon)之间的关系式,得到地面点坐标(lat,lon,h)对应的影像点的平面坐标(xp,yp);Ⅲ、将影像点的平面坐标(xp,yp)代入真实CCD线阵的严密成像几何模型,计算得到影像点坐标(xp,yp,h)的地面点坐标(latp,lonp);Ⅳ、重复步骤Ⅰ~步骤Ⅲ,直到满足迭代条件:

    地面点坐标(latp,lonp)所对应的影像点的平面坐标(xp,yp)即为地面点对应的像点坐标(x',y')。
    本发明由于采取以上技术方案,其具有以下优点:1、本发明是基于物方的拼接方案,其拼接精度不依赖同名点匹配精度。2、采用本发明生成的拼接影像,能够同时提供真实CCD线阵和虚拟CCD线阵的严密成像几何模型。3、本发明通过建立理想直线的虚拟CCD线阵完成影像拼接,所生成的拼接影像符合理想的线中心投影,为无畸变影像,便于后续几何处理?;谝陨嫌诺?,本发明可以广泛应用于资源三号卫星影像处理中。
    附图说明
    图1是本发明的流程图
    图2是灰度双线性内插示意图
    具体实施方式
    下面结合附图和实施例对本发明进行详细的描述。
    如图1所示,本发明资源三号卫星三线阵影像和多光谱影像的拼接方法,具体包括以下步骤:
    1)建立真实CCD线阵的严密成像几何模型;
    ①利用全色或多光谱相机的像主点(x0,y0),像主距f,像元(x,y)大小,确定像元(x,y)在全色或多光谱相机坐标系下的坐标为:
    x-x0y-y0-f.]]>
    ②根据像元(x,y)的行y,采用线性内插方法内插得到成像时间t。通过成像时间t,分别采用拉格朗日内插和线性内插方法得到某成像时刻的轨道XsYsZsWGS84]]>和姿态Q。
    ③根据由步骤②得到的轨道XsYsZsWGS84]]>和姿态Q等参数,建立像点坐标与地面点坐标的对应关系,即真实CCD线阵的严密成像几何模型为:
    XYZWGS84=XsYsZsWGS84+mRbody2WGS84Rcamera2bodyx-x0y-y0-f---(1)]]>
    式(1)中,XYZWGS84]]>是影像拍摄地面点的地面坐标;m为成像比例尺;XsYsZsWGS84]]>是某成像时刻的轨道;Rbody2WGS84是由卫星下传的定姿设备测量获得的描述卫星成像姿态的数据;Rcamera2body是载荷在卫星本体坐标系下的安装矩阵,描述载荷与卫星本体的相对关系;
    根据式(1)以及原始影像的长度和宽度,计算得到真实CCD线阵对应的地面影像范围。
    2)根据真实CCD线阵的严密成像几何模型,建立虚拟CCD线阵的严密成像几何模型;
    ①确定虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标,其具体包括以下步骤:
    利用步骤1)中确定的所有像元(x,y)在全色或多光谱相机坐标系下的坐标(x-x0,y-y0),通过最小二乘法求解虚拟CCD线阵直线方程的解,拟合出最佳的虚拟CCD线阵。
    Ⅰ、将步骤1)确定的所有像元(x,y)在全色或多光谱相机坐标系下的坐标(x-x0,y-y0)作为观测值,虚拟CCD线阵直线方程y=ax+b的系数a和b作为未知数,建立间接平差方程:
    V=AX-L    (2)
    式中,A=x11x21......xn-11xn1,X=ab,L=y1y2...yn-1yn,]]>则X=(ATA)-1ATL。
    假设共有n对观测值(xi,yi),i≤n,根据式(2),采用最小二乘法,求解得到虚拟CCD线阵直线方程的系数a和b。
    Ⅱ、根据原始影像的宽度wreal和真实CCD影像间的重叠像素woverlap,确定虚拟CCD影像的宽度wvirtual为:
    wvirtual=wreal-woverlap    (3)
    Ⅲ、对于虚拟CCD线阵上任一像素s,由于各真实CCD分布在与卫星运动方向垂直的方向上,结合步骤Ⅰ计算得到的系数a和b,得到任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标y'为:y'=ax'+b。结合虚拟CCD影像的宽度wvirtual,得到 任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标为:
    x=(xmax-xmin)/wvirtual·s+xminy=ax+b-f---(4)]]>
    式中,xmax=max(xi),xmin=min(xi)。
    ②建立t时刻卫星姿态模型,其具体包括以下步骤:
    Ⅰ、为解决积分时间跳变的问题,将成像时刻设置为等间隔,建立虚拟CCD成像时间模型为:
    t=t0+Δt·(y-y0)    (5)
    式中,y0为影像起始行,t0是影像起始行的成像时间,Δt为平均积分时间,t为第y行影像的成像时刻。
    Ⅱ、为解决卫星姿态跳变的问题,在重成像时,利用卫星不同时刻姿态拟合的多项式系数ωn,κn,(其中,n=0,1,…)建立t时刻卫星姿态模型为:

    ③将步骤①确定的虚拟CCD线阵上任一像素s所对应的像元在全色或多光谱相机坐标系下的坐标以及步骤②建立的t时刻卫星姿态模型均代入式(1)中,得到虚拟CCD线阵的严密成像几何模型。
    3)建立拼接影像的有理多项式模型;
    ①根据据虚拟CCD线阵的严密成像几何模型,分别计算虚拟CCD影像的上下左右四个角点的地理坐标(lon,lat,h),即某区域的地理范围,从全球1km格网SRTM(Shuttle Radar Topography Mission)数据上读取该区域地理范围的影像,计算该区域的高程范围。
    ②在虚拟CCD的影像空间内,根据虚拟影像大小划分虚拟控制格网,例如划分成200×200×15,根据虚拟CCD线阵的严密成像几何模型计算格网上各个网格点的地理坐标,组成形如(X,Y,P,L,H)的控制点,求解控制点(X,Y,P,L,H)对应的RFM模型:
    Y=NL(P,L,H)DL(P,L,H)X=NS(P,L,H)DS(P,L,H)---(7)]]>
    式中,X,Y为影像坐标,P,L,H分别为地面坐标纬度、经度和高程;NL(P,L,H)为:
    NL(P,L,H)=a1+a2L+a3P+a4H+a5LP+a6LH+a7PH+a8L2+a9P2
               +a10H2+a11PLH+a12L3+a13LP2+a14LH2+a15L2P+a16P3+a17PH2    (8)
               +a18L2H+a19P2H+a20H3
    分别将式(8)中待求的RFM模型的参数ai(i=1,2,…,20)换成bi,ci和di得到DL(P,L,H),Ns(P,L,H)和Ds(P,L,H)。
    4)完成资源三号卫星三线阵影像和多光谱影像的拼接;
    ①确定虚拟CCD影像的高度和宽度,其中,将原始影像的高度设定为虚拟CCD影像的高度,虚拟CCD影像的宽度为wvirtual。
    ②对式(7)进行逆变换并计算得到拼接后影像上任一点O对应高程h的地面点坐标(lat,lon,h)。
    ③利用反算模型,将地面点坐标(lat,lon,h)投影到真实CCD影像上,得到地面点坐标(lat,lon,h)对应的像点坐标(x',y'),其具体包括以下步骤:
    Ⅰ、根据式(1)计算得到原始影像上下左右四个角点的地理坐标(lat,lon),根据像元(x,y)与地理坐标(lat,lon)之间的关系式:
    x=f0+f1lat+f2lony=g0+g1lat+g2lon---(9)]]>
    计算仿射变换参数f0、f1、f2、g0、g1、g2;
    Ⅱ、将步骤②计算得到的地面点坐标(lat,lon,h)代入式(9),得到地面点坐标(lat,lon,h)对应的影像点的平面坐标(xp,yp);
    Ⅲ、将影像点的平面坐标(xp,yp)代入式(1),计算得到影像点坐标(xp,yp,h)的地面点坐标(latp,lonp);
    Ⅳ、重复步骤Ⅰ~步骤Ⅲ,直到满足以下迭代条件时,地面点坐标(latp,lonp)所对应的影像点的平面坐标(xp,yp)即为地面点对应的像点坐标(x',y');

    ④如图2所示,采用距离加权方法,计算得到像点坐标(x',y')对应的灰度值g(x',y')为:
    g(x',y')=(1-dx)·(1-dy)·g1+dx·(1-dy)·g2+dy·(1-dx)·g3+dx·dy·g4    (11)
    ⑤对已知高度和宽度的虚拟CCD影像上的所有像素,重复步骤②~步骤④,完成 资源三号卫星三线阵影像和多光谱影像的拼接。
    上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和方法步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的?;し段е?。

    关于本文
    本文标题:一种资源三号卫星三线阵影像和多光谱影像的拼接方法.pdf
    链接地址://www.4mum.com.cn/p-6130853.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 广西快乐10分 十分钟彩票官网 双色球首页 河南十一选五基本走势图百度彩票 篮球竞猜让分胜负 杭州麻将胡牌组合 安徽25选5走势图 云南十一选五开奖查询 牛牛碰在线观看大香 雷锋彩票安卓 精准单双中特 福利彩票投注站申请 快乐扑克奖金 娱乐场游戏 捕鱼大富翁赢话费官网 陕西快乐十分单式开奖