微波遥感集中编程实习报告二.pdfVIP

微波遥感编程实习报告 一、 实习目的 本次实习希望通过编程实现基于R-D 方程的SAR 影像几何校正，比较校正前后SAR 影像特征与区别，分析距离成像几何原理，掌握R-D 校正步骤。希望通过这次实习了解 SAR 成像几何原理，熟悉距离式成像过程。掌握距离-多普勒方程进行校正的原理和难 点，认识校正所需参数意义与读取方法。理解并掌握斜视SAR 成像中距离式成像的机 理与方法，了解 SAR 图像特殊几何特征成因与校正必要性；掌握基于R-D 构像方程的 SAR 图像像素坐标与坐标间的定位关系；在此基础上编程实现SAR 图像的几何校正过程。 二、 实习容 1. 数据介绍 SAR 影像的头文件信息中包括该幅影像上5 个成像点的XYZ 坐标，第一点的成 像时间和5 个点间的成像时间间隔，影像的第一行成像时间与成像间隔，影像距离 向的成像时间与时间间隔，SAR 影像4 个角点的经纬度以及投影椭球的长短半轴。 这些数据用于确定任意时刻下成像点的位置矢量、速度矢量和加速度矢量，确定对 应的DEM 围以及确定DEM 上各点对应SAR 影像的ij 坐标和灰度值。这些参数已经 存储在parameter.txt里，可直接读取使用，但注意单位转换问题。 DEM 影像是ers2dem.img，可以利用DEM 提供的H 结合经纬度坐标BL，利用 投影椭球信息，计算出成像点的坐标XYZ ，但注意裁剪SAR 对应到DEM 的围。 SAR 影像灰度信息则是存储在DAT_01.001，注意SAR 影像的数据格式，这里存 储的是复数形式，需要计算强度信息作为灰度值。满足公式：强度²=实部²+虚部²。 数据的预处理包括利用影像5 个成像点的XYZ 坐标和第一点成像时间、成像时 间间隔，利用多项式拟合的方法，确定轨道系数ABC，从而可求出任意时刻的卫星 位置矢量，位置矢量求导可得到速度矢量，再求导可得到加速度矢量，这些都是在 之后的R-D 方程中有所利用。注意为解求矩阵方便，可将第一点成像时间当作0 。 其余的数据预处理包括单位转换以及由于双程传输引起的时间转换。 考虑到DEM 的分辨率是90m，而SAR 影像分辨率是30m，所以还要进行DEM 加密，这里使用的方法是双线性插方法实现DEM 加密。 2. R-D 方程 已知R-D 方程由3 条方程组成，包括距离方程、多普勒频率方程和椭球方程， 如下图。 图1 R-D 方程 利用R-D 方程进行几何校正时，使用的是间接校正法。即利用DEM 中各点的 经纬度坐标BL结合DEM影像上的高度信息H，加上投影椭球参数，计算出各点的 坐标XYZ，再由该点的位置矢量结合当时卫星的位置矢量、速度矢量和加速度矢量 代入距离-多普勒方程，确定成像时刻，从而找 SAR影像上对应的（i，j），取出 此处的灰度值，赋给DEM对应的坐标处。最后再插得到的灰度影像，即可得到纠 正后的SAR影像。 三、 实习步骤 下面将根据上述实习容一步步地编程实习 SAR 影像几何校正。几何过程可按概括 为以下的流程图： 开始 得到几何校正 将SAR影像上对应 卫星头文件参 后的SAR影像 的灰度值赋给该地 数读取 并保存 面点坐标 多项式拟合 结束 迭代解求DEM上各 求解拟合参