气象卫星光学遥感器二维点扩散函数处理方法技术

本发明专利技术涉及一种气象卫星光学遥感器二维点扩散函数(PSF)的处理方法，属于气象技术领域。在新定义的极坐标系内，利用在两个相互垂直方向上的己知一维PSF，按照沿极角和极轴两个方向上的嵌套一维插值，得到光学遥感器二维PSF分布，在确保已知方向上PSF信息不丢失的前提下，获得在空间上分布连续、光滑的PSF估计结果。经复原处理后，可以获得图像质量更为清晰的观测目标信息，从而可进一步提高了气象应用的准确性。

【技术实现步骤摘要】

本专利技术涉及气象卫星光学遥感器二维点扩散函数的处理方法，属于气象

技术介绍
在当前全球气候变化、灾害性极端天气事件频发的大背景下，如何充分发挥星载观测，尤其是气象卫星装载的光学遥感器的高品质观测数据，是实现短临天气预报、中短期气候预测等气象防灾减灾的重要观测保障之一，并日益发挥着巨大的社会和经济效益。气象卫星是天基气象观测手段的核心，是气象要素综合观测系统的重要组成部分。对于气象卫星而言，需要重点探测的气象要素包括在特定时间和空间尺度上，大气温度、湿度垂直分布，云、水体(含海洋、湖泊等)、地表等目标在不同谱段内辐射信息及其相对变化特征等。对于光学遥感器而言，一方面，在不同的探测波段上，遥感器接收到的目标辐射可能来自于目标自身的发射辐射或目标反射的太阳辐射，或二者兼而有之，探测结果将作为辨识、分析目标特性的重要依据；另一方面，观测目标在空间上总是有一定的大小，且自身特性也将随时间有一定的变化，这就决定了遥感器需要在合适的时间和空间尺度上对目标进行观测，即设计合适的时间分辨率和空间分辨率。例如，对于以监视日常天气变化为主的静止轨道气象卫星而言，其时间分辨率通常在IOtl-1O1分钟量级；对于以监测大尺度/全球气候与环境变化为主的极地轨道气象卫星而言，其时间分辨率通常在数小时到天量级。同时，空间分辨率的选取通常与目标的空间变化特性有关。例如，对于大气温度和湿度等气象要素，通常在IOtl-1O1公里范围内可认为是相对均匀的，因此，以探测大气温度和湿度信息为主的遥感器，其空分辨率选择在10公里左右较为常见。对于云、地表等目标的探测，由于目标细节变化丰富，综合考虑应用与技术实现约束等条件，空间分辨率常选择在IO2-1O3米量级。遥感器的空间分辨率，表征了仪器对目标在空间尺度上变化特性的刻画能力。简单的来说，空间分辨率越高(即可辨识的目标空间距离越小)，以成像型遥感器为例，其观测图像的细节就越丰富。然而，实际经验表明，即使两个遥感器具有相同的空间分辨率，当观测相同目标时，获得的观测数据/图像质量(如清晰度等)也可能存在差异，这其中的原因主要是不同遥感器的空间响应特性(或称为点扩散函数，Point Spread Function, PSF)不一致。当前，随着国内外众多气象卫星的成功发射以及所装载的遥感器稳定运行，综合应用不同遥感器数据的迫切性越发显得突出。这其中，如何有效抑制不同遥感器间PSF带来的观测差异，是多源卫星观测数据应用面临的重要挑战之一。为了实现这一目标，首先需要在轨准确的评价或测量PSF的空间分布。以气象卫星为例，其主流光学遥感器仍然采用单元或多元探测器。为了实现对大范围地球目标的观测，一般都采用扫描方式来实现。虽然遥感器的PSF有明显的二维分布，但由于米用了扫描成像方式，不同方向(如扫描方向及其垂直方向)上的空间米样和时间间隔存在较大差异，因此，认为PSF在二维方向上各向同性是不够准确的。实际应用中，通常分别计算两个方向，一般取扫描方向和垂直于扫描方向(即卫星飞行方向或观测区域改变方向)，分别计算在两个方向上的一维PSF。同时，在假定两个方向上PSF正交条件下，利用简单的矢量相乘，获得二维PSF分布特性。应用表明，采用经典的维纳滤波方法，若利用以上获得的二维PSF来改善观测数据质量，将在目标变化剧烈的地方出现复原结果异常，这其中的主要原因是由上述方法所获得的PSF在二维空间分布上存在显著的不连续点。为此，本专利技术将提出新的二维PSF构造方法，解决由简单矢量相乘来构造二维PSF引起的空间分布不连续问题，为气象卫星观测定量化应用奠定基础。
技术实现思路
本专利技术的目的是提出一种气象卫星光学遥感器二维点扩散函数(PSF)的处理方法。在新定义的极坐标系内，利用在两个相互垂直方向上的已知一维PSF，按照沿极角和极轴两个方向上的嵌套一维插值，得到光学遥感器二维PSF分布，在确保已知方向上PSF信息不丢失的前提下，获得在空间上分布连续、光滑的PSF估计结果。本专利技术提出的气象卫星光学遥感器二维PSF处理方法，包括以下步骤(I)将两个已知的相互垂直方向上的一维PSF所在方向，分别记为水平⑶方向和垂直方向(V),将两个一维PSF分布的中心重合并记为原点,建立直角坐标系H0V。同时，以直角坐标系HOV的原点为新建立的极坐标系的原点，+H方向为极轴变化的正方向，+H到+V沿逆时针转动的方向为极角变化的正方向，建立极坐标系(r，Θ)。所建直角坐标系HOV和极坐标系Cr, Θ )如图2所示；(2)将两个己知方向上的一维PSF，在直角坐标系内记为PSFid (h)和PSF 1D( V)，则二者在上述定义的极坐标系下可表示为本文档来自技高网...

【技术保护点】
一种气象卫星光学遥感器二维点扩散函数处理方法，其特征在于该方法包括如下步骤：?(1)将两个已知的相互垂直方向上的一维PSF所在方向，分别记为水平(H)方向和垂直方向(V)，将两个一维PSF分布的中心重合并记为原点，建立直角坐标系HOV；同时，以直角坐标系HOV的原点为新建立的极坐标系的原点，+H方向为极轴变化的正方向，+H到+V沿逆时针转动的方向为极角变化的正方向，建立极坐标系(r，θ)，所建直角坐标系HOV和极坐标系(r，θ)如图2所示。?(2)将两个已知方向上的一维PSF，在直角坐标系内记为PSF1D(h)和PSF1D(v)，则二者在上述定义的极坐标系下可表示为：?式中，h和v分别表示一维PSF值所在位置到原点的距离。?(3)在极坐标系(r，θ)内分别选取相同矢径，r可取1，2，3……，极角θ分别取0，π/2，π，3π/2，根据需要选取一维插值函数(Interpolation?Function，IPF)，得到r取定值(1，2，3……)时任意极角θ下的PSF2D值，并有?(4)对于直角坐标系内任一待估计的二维PSF值对应的位置(h0，v0)，计算其在极坐标系(r，θ)下的坐标(r0，θ0)，并有?利用步骤(3)得到的PSF2D(r，θ)|r∈{1，2，3，......}结果，可获得任一位置(h0，v0)的PSF2D(h0，v0)值，并有?PSF2D(h0，v0)＝IPF(PSF2D(r，θ0)|r∈{1，2，3，......}，r0)?(5)将步骤(4)计算得到的PSF2D(h0，v0)记录在二维PSF分布表中；?(6)重复步骤(4)、(5)，得到光学遥感器二维PSF分布值。?dest_path_FSB00001021531400011.tif,dest_path_FSB00001021531400012.tif,dest_path_FSB00001021531400013.tif,dest_path_FSB00001021531400014.tif,dest_path_FSB00001021531400015.tif...

【技术特征摘要】
1.ー种气象卫星光学遥感器ニ维点扩散函数处理方法，其特征在于该方法包括如下步骤 (1)将两个已知的相互垂直方向上的ー维PSF所在方向，分别记为水平(H)方向和垂直方向(V),将两个ー维PSF分布的中心重合并记为原点,建立直角坐标系HOV ;同时，以直角坐标系HOV的原点为新建立的极坐标系的原点，+H方向为极轴变化的正方向，+H到+V沿逆时针转动的方向为极角变化的正方向，建立极坐标系(r，0)，所建直角坐标系HOV和极坐标系Cr，0 )如图2所示。(2)将...

【专利技术属性】
技术研发人员：郭强，
申请(专利权)人：国家卫星气象中心，
类型：发明
国别省市：