| 主权项 |
一种基于时空联合调制干涉成像光谱仪数据坏像元的校正方法，其特征在于：所述校正方法包括以下几个步骤：步骤一：获得时空联合调制型干涉成像光谱仪中坏像元的位置；步骤二：通过时空联合调制型干涉成像光谱仪，获取包含感兴趣目标地物点的图像，即推扫图像序列；采集图像的数量根据香农采样定理确定，采集图像的数量为N，对于感兴趣目标地物点，每一幅图像对应一个光程差，N幅图像对应N个不同的光程差；步骤三：通过数据处理，生成一个数据立方体；数据立方体的三个维度分别为两个空间维和一个光谱维，两个空间维的数据构成一幅幅的推扫图像，光谱维的任一组数据都表示某一个地物点的对应不同光程差的完整的干涉数据；假设单帧推扫图像大小为N1*N2，目标地物点在第一帧上的位置为(x，y)，则目标地物点在第K帧上的位置为(x+K，y)；当满足香农采样定理时，数据立方体的大小为(N1+N)*N2*N；步骤四：判断目标地物点的完整的干涉强度数据是否存在坏点，如果存在，确定坏点位置；目标地物点在不同帧图像中，出现在不同坐标位置，理想情况下，位置连接为一条沿着推扫方向的直线段；结合坏像元的位置数据，判断目标地物点的完整的干涉强度数据是否经过坏像元点，并确定其在完整的干涉强度数据中的位置；步骤五：获取感兴趣目标地物所在区域的数据立方体数据，并对坏像元位置进行插值校正；假设感兴趣目标地物点的完整干涉强度曲线上有N个坏像元位置，取三个坏像元A、B和C，A、B、C三个坏点位置到零光程差位置的距离分别为x<sub>1</sub>、x<sub>2</sub>和x<sub>3</sub>，设定阈值T<sub>1</sub>和T<sub>2</sub>，且0<x<sub>1</sub><T<sub>1</sub><x<sub>2</sub><T<sub>2</sub><x<sub>3</sub>，阈值的确定为：由于数据立方体中光谱维的每一列数据都一一对应于地面上一个地物点的干涉强度曲线，已知目标地物点位置，从数据立方体中提取目标地物的干涉强度曲线，设定T<sub>1</sub>为次级大所在波峰后的第二个零点位置到零光程差位置之间的距离，T<sub>2</sub>为11级大波峰后的第二个零点位置到零光程差位置之间的距离；插值校正的方法具体为：当坏像元为A点时，在空间维对坏像元的强度值进行校正，在交轨方向和沿轨方向对数据进行三次立方卷积插值，分别获得坏像元位置在不同方向的强度估计值；设交轨方向权值为p<sub>1</sub>，沿轨方向权值为p<sub>2</sub>，p<sub>1</sub>>p<sub>2</sub>，p<sub>1</sub>+p<sub>2</sub>＝1；当坏像元为C点时，此时坏像元位置远离零光程差，则在光谱维对坏像元的强度值进行校正，提取目标地物点完整的干涉强度数据，对坏像元位置进行三次立方卷积插值，获得坏像元位置的强度估计值；当坏像元为B点时，则使用光谱维和空间维加权平均方法进行插值，且对于不同的x<sub>2</sub>位置，光谱维的权重因子q<sub>1</sub>和空间维的权重因子q<sub>2</sub>的取值不同，q<sub>1</sub>+q<sub>2</sub>＝1，且满足距离为T<sub>1</sub>时，q<sub>1</sub>＝0,在距离为T<sub>2</sub>时，q<sub>2</sub>＝0；步骤六：将感兴趣目标地物点的插值后的完整干涉强度数据，进行傅里叶变换复原得到该感兴趣目标地物点的光谱。 |
