卫星产品数据解析

stream2011

  本帖卫星数据产品主要分静止气象卫星图像产品，极轨气象卫星图像产品，格点场，离散场，其中极轨和离散场本文中没考虑。如下图为格点场数据解析，左边为解析结果，右边为meteoinfo软件显示结果：file:///F:\文档\378628923\Image\C2C\0YXG0B79RIXI90S4FZ${U22.png
下图为静止气象卫星图像产品，左边为解析结果，右边为meteoinfo软件显示结果：


部分源码如下，详见py文件：

def get_stat_data(self):  # 按照像素读取
    with open(self.filepath, "rb") as f:
        #  ----------------------  第一级文件头  -------------------------
        print(f.__sizeof__())
        self.ini_stat_first(f)
        #  ---------------------- 静止气象卫星 第二级文件头  -------------------------
        if self.product_kind == 1:  # 静止气象卫星图像产品
            self.staic_stat_second(f)
            self.get_static_data(f)

        elif self.product_kind == 2:  # 极轨气象卫星图象产品
            self.polar_stat_second(f)
        elif self.product_kind == 3:  # 格点场
            self.grid_stat_second(f)
            self.get_grid_data(f)

        elif self.product_kind == 4:  # 离散场
            self.dispersed_stat_second(f)
if __name__ == '__main__':
    filepath = r'sat96/FY2C_TBB_IR1_OTG_20080605_0600.AWX'
    filepath = r'sat96/FY2C_TBB_IR1_OTG_20080605_0700.AWX'

    filepath = r'sat96/FY2C_ANI_IR1_R01_20080605_1400.AWX'

    path = r'D:\PycharmProjects\study\qixiangtai_2020\satellite_data_resolving\sat96'
    filelist = fu.get_filelist(path, '.AWX')
    for f in filelist:
        cs = CStat2004(f)
        if cs.product_kind == 1:  # 静止型
            print('静止型  lat:   ' + str(cs.South_lat) + '  ——   ' + str(cs.North_lat))
            print('静止型  lon:   ' + str(cs.West_lon) + '  ——  ' + str(cs.East_lon))
        elif cs.product_kind == 3:  # 格点型
            print('格点型  lat:   ' + str(cs.rightdown_lat) + '  ——   ' + str(cs.leftup_lat))
            print('格点型  lon:   ' + str(cs.leftup_lon) + '  ——  ' + str(cs.rightdown_lon))

        plt.imshow(cs.data)
        # plt.grid()
        # plt.show()
        plt.draw()
        plt.pause(5)  # 间隔的秒数：6s


来自群组: 农业气象，服务三农