遥感原理学习补充材料.docVIP

一、星历表法解算卫星坐标 有了轨道6参数，可解算某一瞬间的卫星坐标，还须测定卫星在该瞬间的精确时间，解算卫星位置。 （1）卫星在地心直角坐标系中的坐标 （2）卫星在大地地心直角坐标系中的坐标 （3）卫星在地理坐标系中的坐标 真近点角v和偏近点角E的关系如下图： （1）真近点角V，AQ与AS夹角 （2）假想一颗虚拟卫星S///以恒角速度在轨道上运动，S///的真近点角称为平近点角 （3）偏近点角E 首先以轨道椭圆O为圆心，分别以椭圆的长半径a、短半径b作两个圆。然后通过卫星S作椭圆短轴的平行线与a为半径的大圆交于S/, 再通过S作长轴的平行线与小圆交于S//。 根据椭圆特性，S/, S//, O三点位于同一条直线上。E即OQ旋转至OS/ 的角度。平近点角M和真近点角V是通过偏近点角E联系的。 计算卫星在轨道坐标系中的位置 首先建立一个轨道坐标系，该坐标系的坐标原点位于地心，位于轨道平面上， 轴和轨道平面的法线矢量重合。轨道坐标系是一个右手坐标系。计算步骤如下： 用下式计算平近点角 为卫星过近地点的时刻；为卫星的平均角速度，用下式计算： 为轨道椭圆的长半径，（注：引力常数，此为地球质量） 解开普勒方程,计算偏近点角 解算时采用角度制， （离心率） 代入开普勒方程反复迭代，直至＜时为止。 （当偏心率很小时，迭代法的收敛速度很快） 计算卫星至地心的距离 计算真近点角 计算卫星在轨道坐标系中的坐标 或跳过3、4直接计算： 轨道坐标和大地坐标的换算 将上式化算到大地坐标系中去，一是用地心空间直角坐标系（）来表示点的位置，二是用经纬度和大地高（）来表示点的位置，只要确定椭球体的参数和定位，（）和（）之间就可以换算。 轨道坐标系只需经三次旋转即可和大地坐标系（）重合。首先绕轴反时针旋转一个角，使旋至（指向升交点）。再绕反时针旋转角，这样与重合。最后绕反时针旋转一个（）,这两个坐标系就重合了。角（P25，图2-1中角）是轴与春分点X方向的夹角，即为格林尼治恒星时角。于是有 其中： 地心坐标系与地理经纬度坐标系间转换 或写为： 地理坐标与地图坐标间的转换（略） 【例题】已知卫星的轨道根数如下，计算卫星在时的位置。 在时刻时角，。 1. 计算平近点角M 将GM代入求卫星的平均角速度 弧度/秒 弧度 化为度计算： 2. 解开普勒方程求偏近点角E 解一：由于偏心率e较小所以采用迭代法进行，计算时采用角度制，所以先求出再依次迭代： 解二： 对开普勒方程两侧求微分： （分母以弧度计算） 取M作为初值代入方程： 再求E2 : 然后求E3 : 结果是相同的。其它再往下做。 3. 计算时卫星到地心的距离 4. 计算真近点角 计算卫星在轨道坐标系中的坐标 亦可直接用偏近点角E来计算卫星的坐标 最后再将求得的卫星位置从轨道坐标转换到大地坐标系中去： 就求得了卫星在9：00：00.0000时的大地坐标。计算大地坐标时未考虑极移改正。 参看教材59页图4-2，其它计算（略）。 二、图像的频谱表示与傅里叶变换 将图像从空间域变入频率域是采用傅里叶变换，反之，则采用傅里叶逆变换。当M与N不相等时采用以下变换对： 傅里叶变换： 傅里叶逆变换： 当M=N时，即为方阵，此时傅里叶变换对如下式： 【例题】对下列图像数据，求二维傅里叶变换。 1 2 3 4 4 1 2 3 3 4 1 2 2 3 4 1 解：1.基本公式： 复数表达： 欧拉公式： 若： 则： 傅氏变换函数为周期函数：如： 2.由于是方阵，故使用