CN116185616B 一种fy-3d mersi l1b数据自动化再处理方法 （重庆市气象科学研究所(重庆市生态气象和卫星遥感中心、重庆市农业气象中心)）

号一种FY-3DMERSIL1B数据自动化再处理方法本发明针对风云三号D(FY_3D)气象卫星接种面向FY_3DMERSIL1B数据的自动化再处理方法，其采用Python开发平台和GDAL栅格处理技归档1km和250m空间分辨率等经纬网格多通道图2步骤1：通过风云三省级直收站实时接收并预处理为FY_3D步骤61：获取FY_3D_MERSI待监测图像区域的经纬度坐标范围信息，步骤64：将获取的像元太阳和卫星观测角度、计算步骤21：定期扫描FY_3DMERSIL1B轨道时序观测数据目录下所接3sensor其中，C1和C2为整合后的计算系数，h=6.62606876e_34J.s为普朗克常量，c=步骤51：根据读取的不同空间分辨率的像元地理坐标，获min4步骤53：基于所述构建原始图像像元地理坐标到投影网格行列号之间的关系通过5[0001]本发明涉及FY_3DMERSIL1B数据再处理技术方法领域，具体涉及一种FY_3D星掩星探测仪(GNOS)、红外高光谱大气探测仪(HIRAS)、近红外高光谱温室气体监测仪数据的遥感应用研究，实现了基于FY_3DMERSI_II远红外数据火情监测研究(郑伟、陈洁象部门主要采用FY3DMERSI_II数据产品开展本地气象遥感实[0003]因FY_3DMERSIL1B数据没有等经纬投影信息，光谱波段数据主要以整型辐射亮度灰度值(DigitalNumber,DN)存储，非遥感技术专业人员解析及应用该数据仍需对数据数据的再处理加工方式需大量人工干预和设置，没有形成数据处理的自动化处理工程模平台数据传输进行对接，采用Python和GDAL栅格技术对轨道数据进行再处理，实现FY_3D据处理效率和自动化水平，自动生成和归档250m和1km两种空间分辨率多通道栅格格式遥6[0014]步骤21：定期扫描FY_3DMERSIL1B轨道时序观测数据目录下所接收到的FY_3D7C2为整合后的计算系数，h＝6.62606876e_34J.s为普朗克常量，c=[0049]步骤61：获取FY_3D_MERSI待监测图像区域的经纬度坐标范围信息，对FY_3D8归档250m和1km两种空间分辨率多通道波段栅格遥感图像；通过实验表明经本发明方法处[0061]图7为Band20_Band25中长波热红外辐射TEB波段亮温MAE和RMSE；其中图7(a)为[0062]图8为将原处理方法和本发明处理方法进行拟合的结果；其中图8(a)为波段[0066]本发明要处理的数据为风云三省级直收站系统接收和预处理的FY_3DMERSIL1B[0067]风云三省级直收站接收和预处理的FY_3DMERSIL1B轨道过境观测数据主要有四国家卫星气象中心数据服务网上可查询到数据说明，空间分辨率1km定标观测数据集主要包含4个子集，空间分辨率250m数据集含有6个子集以及不同分辨率地理坐标文件。表1为9[0088]将(4)式带入(3)式简化得到通道大气表观反射率计算公式(5)，将大气表观反射备2个250m空间分辨率热红外波段。该热红外波段数据存储为未定标辐射亮度整型量化值[0094]此处的热红外波段辐射亮度数据的单位是mC2为整合后的计算系数，h＝6.62606876e_34J.s为普朗克常量，c=率图像的最大纬度Latmax和最小经度Lonmin和图像像元行高Res_Line和像元列宽Res_[0110]通过(9)式和(10)式建立了原卫星图像地理坐标到新投影网格图像行列号之间映[0112]本发明设置的FY_3DMERSI图像数据裁剪区范围为104.0E～112.0E和26.0N~新组合到一起，构建1km分辨率大小为29*800*800和250m分辨率大小为10*3200*3200的三[0117]FY3DMERSIL1数据自动再处理技术能够与省级直收站接收预处理平台紧密衔块耦合连接和集约化管理，实现各模块间异步条件触发方式运行模式，提高FY_3DMERSIL1B级数据再处理的智能化水平。数据自动再处理技术流程见图2，输入端为FY_3DMERSI[0119]为了进一步验证本发明所提自动处理方法的有效性，选取了接收完整的2020年8[0121]经本发明方法处理前和处理后的平均数据大小情况见表4，处理加工后的多通道TIF图像的数据量整体要小于处理前数据量，1km空间分辨率数据处理后数据量减少了约[0126]FY_3DMERSIL1数据自动再处理耗时时效变化趋势见图3，可以看出本发明处理[0129]Table7themainmodelsaverageruntimebetweenoldmethodandnew[0132]以2020年8月14日13时33分的FY_3DMERSIL1B数据为例较。从原处理方式和本发明技术的band13(0.709μm)两幅图像的空间分布上看(如图4所[0134]通过计算光谱反射太阳光波段RSB和热辐射波段TEB的像元平均绝对误差MAE和均外辐射TEB波段亮温MAE和RMSE(如图7所示)，发现辐射亮温误差整体不大，其MAE不超过以对前述各实施例所记载的技术方案进行修改，或者对其