metradar 0.1.6__py3-none-any.whl → 0.1.8.2__py3-none-any.whl

metradar/{decode_pup_rose.py → io/decode_pup_rose.py}

@@ -1,12 +1,5 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-'''
-@File    :   decode_pup_rose.py
-@Time    :   2023/04/27 23:12:54
-@Author  :   Wenjian Zhu 
-@Version :   1.0
-@Email   :   kevin2075@163.com
-'''
+
+# _*_ coding: utf-8 _*_
 
 '''
 解析pup_rose文件
@@ -16,18 +9,18 @@
 # %%
 import os
 import numpy as np
-import geo_transforms_pyart as geotrans
+import metradar.util.geo_transforms_pyart as geotrans
 from datetime import datetime,timedelta
-# from rose_structer import _unpack_from_buf,_structure_size,SIT_HEADER_BLOCK,GENERIC_HEADER,\
-#                 SITE_CONFIG,TASK_CONFIG,SCAN_CONFIG,PRODUCT_HEADER_BLOCK,PRODUCT_ID_NAME_MAP,\
-#                 PRODUCT_DEPENDENT_PARAMETER,STORM_MOTION_BLOCK,STORM_FST_HIS_NUM,STORM_FST_HIS_BLOCK, \
-#                 STORM_PROPERTY,STORM_TRACK_PARAM,STORM_COMPONENT
+from .rose_structer import _unpack_from_buf,_structure_size
+from .rose_structer import *
+
 import xarray as xr
-import rose_structer as rs
+from pyart.core import Radar
+from pyart.config import FileMetadata
+from pyart.io.common import make_time_unit_str
 
 class READ_ROSE(object):
 
-
     def __init__(self,):
         self.stiinfo = None
         self.mesoinfo = None
@@ -149,81 +142,81 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
         # pprint(prod_type)
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
 
         # 获取产品数据块
 
         #获取STI头信息
-        dic_sti_header = rs._unpack_from_buf(buf, pos, rs.SIT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.SIT_HEADER_BLOCK)
+        dic_sti_header = _unpack_from_buf(buf, pos, SIT_HEADER_BLOCK)
+        pos = pos + _structure_size(SIT_HEADER_BLOCK)
 
         # 风暴追踪信息块
         # 风暴移动信息
         storm_motion_block=[]
         for nn in np.arange(dic_sti_header['storm_number']):
-            storm_motion_block.append(rs._unpack_from_buf(buf, pos, rs.STORM_MOTION_BLOCK))
-            pos = pos + rs._structure_size(rs.STORM_MOTION_BLOCK)
+            storm_motion_block.append(_unpack_from_buf(buf, pos, STORM_MOTION_BLOCK))
+            pos = pos + _structure_size(STORM_MOTION_BLOCK)
 
         # 风暴预报信息
         all_storm_fst_block=[]
         for nn in np.arange(dic_sti_header['storm_number']):
-            dic_storm_fst_num = rs._unpack_from_buf(buf, pos, rs.STORM_FST_HIS_NUM)
-            pos = pos + rs._structure_size(rs.STORM_FST_HIS_NUM)
+            dic_storm_fst_num = _unpack_from_buf(buf, pos, STORM_FST_HIS_NUM)
+            pos = pos + _structure_size(STORM_FST_HIS_NUM)
 
             storm_fst_block=[]
             for nn in np.arange(dic_storm_fst_num['position_number']):
-                storm_fst_block.append(rs._unpack_from_buf(buf, pos, rs.STORM_FST_HIS_BLOCK))
-                pos = pos + rs._structure_size(rs.STORM_FST_HIS_BLOCK)
+                storm_fst_block.append(_unpack_from_buf(buf, pos, STORM_FST_HIS_BLOCK))
+                pos = pos + _structure_size(STORM_FST_HIS_BLOCK)
             all_storm_fst_block.append(storm_fst_block)
 
         # 风暴历史信息
         all_storm_his_block=[]
         for nn in np.arange(dic_sti_header['storm_number']):
-            dic_storm_his_num = rs._unpack_from_buf(buf, pos, rs.STORM_FST_HIS_NUM)
-            pos = pos + rs._structure_size(rs.STORM_FST_HIS_NUM)
+            dic_storm_his_num = _unpack_from_buf(buf, pos, STORM_FST_HIS_NUM)
+            pos = pos + _structure_size(STORM_FST_HIS_NUM)
 
             storm_his_block=[]
             for nn in np.arange(dic_storm_his_num['position_number']):
-                storm_his_block.append(rs._unpack_from_buf(buf, pos, rs.STORM_FST_HIS_BLOCK))
-                pos = pos + rs._structure_size(rs.STORM_FST_HIS_BLOCK)
+                storm_his_block.append(_unpack_from_buf(buf, pos, STORM_FST_HIS_BLOCK))
+                pos = pos + _structure_size(STORM_FST_HIS_BLOCK)
             all_storm_his_block.append(storm_his_block)
 
         # 风暴属性表块数据
         # 风暴属性
         all_storm_prop=[]
         for nn in np.arange(dic_sti_header['storm_number']):
-            dic_storm_prop = rs._unpack_from_buf(buf, pos, rs.STORM_PROPERTY)
-            pos = pos + rs._structure_size(rs.STORM_PROPERTY)
+            dic_storm_prop = _unpack_from_buf(buf, pos, STORM_PROPERTY)
+            pos = pos + _structure_size(STORM_PROPERTY)
 
             # 将风暴数字ID转换为字符串ID
             dic_storm_prop['id_char'] = self.get_id_char(dic_storm_prop['id'])
@@ -237,13 +230,13 @@ class READ_ROSE(object):
         # 风暴构成表
         all_storm_comp=[]
         for nn in np.arange(dic_sti_header['storm_number']):
-            dic_storm_comp = rs._unpack_from_buf(buf, pos, rs.STORM_COMPONENT)
-            pos = pos + rs._structure_size(rs.STORM_COMPONENT)
+            dic_storm_comp = _unpack_from_buf(buf, pos, STORM_COMPONENT)
+            pos = pos + _structure_size(STORM_COMPONENT)
             all_storm_comp.append(dic_storm_comp)
 
         # 风暴追踪适配数据
-        dic_storm_track_param = rs._unpack_from_buf(buf, pos, rs.STORM_TRACK_PARAM)
-        pos = pos + rs._structure_size(rs.STORM_TRACK_PARAM)
+        dic_storm_track_param = _unpack_from_buf(buf, pos, STORM_TRACK_PARAM)
+        pos = pos + _structure_size(STORM_TRACK_PARAM)
 
 
         allresult={}
@@ -331,49 +324,49 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
         # pprint(prod_type)
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
 
         # 获取产品数据块
 
         #获取SS头信息
-        dic_ss_header = rs._unpack_from_buf(buf, pos, rs.SS_HEAD_BLOCK)
-        pos = pos + rs._structure_size(rs.SS_HEAD_BLOCK)
+        dic_ss_header = _unpack_from_buf(buf, pos, SS_HEAD_BLOCK)
+        pos = pos + _structure_size(SS_HEAD_BLOCK)
 
         # SS结构信息
         ss_tab=[]
         for nn in np.arange(dic_ss_header['storm_number']):
-            tmp_tab = rs._unpack_from_buf(buf, pos, rs.SS_TAB)
-            pos = pos + rs._structure_size(rs.SS_TAB)
+            tmp_tab = _unpack_from_buf(buf, pos, SS_TAB)
+            pos = pos + _structure_size(SS_TAB)
             # 将风暴数字ID转换为字符串ID
             tmp_tab['id_char'] = self.get_id_char(tmp_tab['storm_id'])
             
@@ -384,29 +377,29 @@ class READ_ROSE(object):
         cell_trend=[]
         for nn in np.arange(dic_ss_header['storm_number']):
             cur_cellinfo={}
-            cell_info=rs._unpack_from_buf(buf, pos, rs.CELL_TREND)
-            pos = pos + rs._structure_size(rs.CELL_TREND)
+            cell_info=_unpack_from_buf(buf, pos, CELL_TREND)
+            pos = pos + _structure_size(CELL_TREND)
             cur_cellinfo['head_info'] = cell_info
             # 历史体扫信息
             tmpcell_his=[]
             for im in np.arange(cell_info['his_vol_num']):
-                tmpcell_his.append(rs._unpack_from_buf(buf, pos, rs.HIS_VOL))
-                pos = pos + rs._structure_size(rs.HIS_VOL)
+                tmpcell_his.append(_unpack_from_buf(buf, pos, HIS_VOL))
+                pos = pos + _structure_size(HIS_VOL)
             cur_cellinfo['cell_info'] = tmpcell_his
 
             cell_trend.append(cur_cellinfo)
 
         # 风暴段适配数据
-        seg_adapt = rs._unpack_from_buf(buf, pos, rs.SEG_ADAPT)
-        pos = pos + rs._structure_size(rs.SEG_ADAPT)
+        seg_adapt = _unpack_from_buf(buf, pos, SEG_ADAPT)
+        pos = pos + _structure_size(SEG_ADAPT)
 
         # 风暴质心适配数据
-        centroid_adapt = rs._unpack_from_buf(buf, pos, rs.CENTROIDS_ADAPT)
-        pos = pos + rs._structure_size(rs.CENTROIDS_ADAPT)
+        centroid_adapt = _unpack_from_buf(buf, pos, CENTROIDS_ADAPT)
+        pos = pos + _structure_size(CENTROIDS_ADAPT)
 
         # 风暴追踪适配数据
-        storm_track_adapt = rs._unpack_from_buf(buf, pos, rs.STORM_TRACK_PARAM)
-        pos = pos + rs._structure_size(rs.STORM_TRACK_PARAM)
+        storm_track_adapt = _unpack_from_buf(buf, pos, STORM_TRACK_PARAM)
+        pos = pos + _structure_size(STORM_TRACK_PARAM)
         
 
         # 将所有的方位角和距离转换成经纬度
@@ -436,59 +429,58 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
         # pprint(prod_type)
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 获取产品数据块
 
         #获取中气旋头信息
-        dic_meso_header = rs._unpack_from_buf(buf, pos, rs.MESO_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.MESO_HEADER_BLOCK)
+        dic_meso_header = _unpack_from_buf(buf, pos, MESO_HEADER_BLOCK)
+        pos = pos + _structure_size(MESO_HEADER_BLOCK)
 
         # 中气旋表块
         meso_tab=[]
         for nn in np.arange(dic_meso_header['meso_number']):
-            meso_tab.append(rs._unpack_from_buf(buf, pos, rs.MESO_TABLE))
-            pos = pos + rs._structure_size(rs.MESO_TABLE)
+            meso_tab.append(_unpack_from_buf(buf, pos, MESO_TABLE))
+            pos = pos + _structure_size(MESO_TABLE)
 
         # 中气旋特征表
         meso_feature_tab=[]
         for nn in np.arange(dic_meso_header['feature_number']):
-            meso_feature_tab.append(rs._unpack_from_buf(buf, pos, rs.MESO_FEATURE_TAB))
-            pos = pos + rs._structure_size(rs.MESO_FEATURE_TAB)
+            meso_feature_tab.append(_unpack_from_buf(buf, pos, MESO_FEATURE_TAB))
+            pos = pos + _structure_size(MESO_FEATURE_TAB)
 
         # 中气旋适配数据
-        dic_meso_adapt_param = rs._unpack_from_buf(buf, pos, rs.MESO_ADAPTATION_DATA)
-        pos = pos + rs._structure_size(rs.MESO_ADAPTATION_DATA)
-
+        dic_meso_adapt_param = _unpack_from_buf(buf, pos, MESO_ADAPTATION_DATA)
+        pos = pos + _structure_size(MESO_ADAPTATION_DATA)
 
         
         for nn in np.arange(dic_meso_header['meso_number']):
@@ -534,52 +526,52 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
         # pprint(prod_type)
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 获取产品数据块
 
         #获取TVS头信息
-        dic_tvs_header = rs._unpack_from_buf(buf, pos, rs.TVS_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.TVS_HEADER_BLOCK)
+        dic_tvs_header = _unpack_from_buf(buf, pos, TVS_HEADER_BLOCK)
+        pos = pos + _structure_size(TVS_HEADER_BLOCK)
 
         # TVS表块
         tvs_tab=[]
         for nn in np.arange(dic_tvs_header['tvs_number']):
-            tvs_tab.append(rs._unpack_from_buf(buf, pos, rs.TVS_TAB))
-            pos = pos + rs._structure_size(rs.TVS_TAB)
+            tvs_tab.append(_unpack_from_buf(buf, pos, TVS_TAB))
+            pos = pos + _structure_size(TVS_TAB)
 
         # TVS适配数据
-        dic_tvs_adapt_param = rs._unpack_from_buf(buf, pos, rs.TVS_ADAPTATION_DATA)
-        pos = pos + rs._structure_size(rs.TVS_ADAPTATION_DATA)
+        dic_tvs_adapt_param = _unpack_from_buf(buf, pos, TVS_ADAPTATION_DATA)
+        pos = pos + _structure_size(TVS_ADAPTATION_DATA)
 
 
         # 将所有的方位角和距离转换成经纬度
@@ -616,50 +608,50 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
             return None
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_radial_header = rs._unpack_from_buf(buf, pos, rs.RADIAL_HEADER)
-        pos = pos + rs._structure_size(rs.RADIAL_HEADER)
+        dic_radial_header = _unpack_from_buf(buf, pos, RADIAL_HEADER)
+        pos = pos + _structure_size(RADIAL_HEADER)
 
-        print('数据类型: %d-%s'%(dic_radial_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
+        print('数据类型: %d-%s'%(dic_radial_header['data_type'],PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取径向数据RADIAL_DATA
         data=[]
         data_azi=[]
         binnum = None
         for nn in np.arange(dic_radial_header['radial_number']):
-            radial_data_info = rs._unpack_from_buf(buf, pos, rs.RADIAL_DATA)
-            pos = pos + rs._structure_size(rs.RADIAL_DATA)
+            radial_data_info = _unpack_from_buf(buf, pos, RADIAL_DATA)
+            pos = pos + _structure_size(RADIAL_DATA)
             binnum = radial_data_info['num_bins']
             rdata = buf[pos:pos+dic_radial_header['bin_length']*radial_data_info['num_bins']]
             data.extend(rdata)
@@ -677,11 +669,11 @@ class READ_ROSE(object):
         print(nx,ny)
         # 将data转成numpy数组
         data = data.reshape((nx,ny))
-
+        
+        
         # decode data
         # ppiarray = (ppiarray - dic_radial_header['offset']) / dic_radial_header['scale'] 
 
-        
         # 将方位角坐标转换成经纬度坐标
         sweep_azimuths =np.array(data_azi)
         ele = dic_prod_param['ele']
@@ -795,42 +787,42 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
             return None
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         print('产品类型: %d-%s'%(dic_prod_header['product_type'],prod_type))
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_grid_header = rs._unpack_from_buf(buf, pos, rs.GRID_HEADER)
-        pos = pos + rs._structure_size(rs.GRID_HEADER)
-        print('数据类型: %d-%s'%(dic_grid_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
+        dic_grid_header = _unpack_from_buf(buf, pos, GRID_HEADER)
+        pos = pos + _structure_size(GRID_HEADER)
+        print('数据类型: %d-%s'%(dic_grid_header['data_type'],PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取栅格数据GRID_DATA
         
@@ -907,41 +899,41 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         print('产品类型: %d-%s'%(dic_prod_header['product_type'],prod_type))
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_grid_header = rs._unpack_from_buf(buf, pos, rs.GRID_HEADER)
-        pos = pos + rs._structure_size(rs.GRID_HEADER)
-        print('数据类型: %d-%s'%(dic_grid_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
+        dic_grid_header = _unpack_from_buf(buf, pos, GRID_HEADER)
+        pos = pos + _structure_size(GRID_HEADER)
+        print('数据类型: %d-%s'%(dic_grid_header['data_type'],PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取栅格数据GRID_DATA
         
@@ -1018,41 +1010,41 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         print('产品类型: %d-%s'%(dic_prod_header['product_type'],prod_type))
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_grid_header = rs._unpack_from_buf(buf, pos, rs.GRID_HEADER)
-        pos = pos + rs._structure_size(rs.GRID_HEADER)
-        print('数据类型: %d-%s'%(dic_grid_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
+        dic_grid_header = _unpack_from_buf(buf, pos, GRID_HEADER)
+        pos = pos + _structure_size(GRID_HEADER)
+        print('数据类型: %d-%s'%(dic_grid_header['data_type'],PRODUCT_DATA_TYPE[dic_grid_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取栅格数据GRID_DATA
         
@@ -1135,49 +1127,49 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         print('产品类型: %d-%s'%(dic_prod_header['product_type'],prod_type))
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_radial_header = rs._unpack_from_buf(buf, pos, rs.RADIAL_HEADER)
-        pos = pos + rs._structure_size(rs.RADIAL_HEADER)
-        print('数据类型: %d-%s'%(dic_radial_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
+        dic_radial_header = _unpack_from_buf(buf, pos, RADIAL_HEADER)
+        pos = pos + _structure_size(RADIAL_HEADER)
+        print('数据类型: %d-%s'%(dic_radial_header['data_type'],PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取径向数据RADIAL_DATA
         data=[]
         data_azi=[]
         binnum = None
         for nn in np.arange(dic_radial_header['radial_number']):
-            radial_data_info = rs._unpack_from_buf(buf, pos, rs.RADIAL_DATA)
-            pos = pos + rs._structure_size(rs.RADIAL_DATA)
+            radial_data_info = _unpack_from_buf(buf, pos, RADIAL_DATA)
+            pos = pos + _structure_size(RADIAL_DATA)
             binnum = radial_data_info['num_bins']
             # rdata = buf[pos:pos+dic_radial_header['bin_length']*radial_data_info['num_bins']]
             rdata = np.frombuffer(buf[pos:pos+dic_radial_header['bin_length']*radial_data_info['num_bins']], '<u2')
@@ -1290,50 +1282,50 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
             return None
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         print('产品类型: %d-%s'%(dic_prod_header['product_type'],prod_type))
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 读取产品径向头信息
-        dic_radial_header = rs._unpack_from_buf(buf, pos, rs.RADIAL_HEADER)
-        pos = pos + rs._structure_size(rs.RADIAL_HEADER)
-        print('数据类型: %d-%s'%(dic_radial_header['data_type'],rs.PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
+        dic_radial_header = _unpack_from_buf(buf, pos, RADIAL_HEADER)
+        pos = pos + _structure_size(RADIAL_HEADER)
+        print('数据类型: %d-%s'%(dic_radial_header['data_type'],PRODUCT_DATA_TYPE[dic_radial_header['data_type']]))
         # Value=(Code-Offset)/Scale
         # 读取径向数据RADIAL_DATA
         data=[]
         data_azi=[]
         binnum = None
         for nn in np.arange(dic_radial_header['radial_number']):
-            radial_data_info = rs._unpack_from_buf(buf, pos, rs.RADIAL_DATA)
-            pos = pos + rs._structure_size(rs.RADIAL_DATA)
+            radial_data_info = _unpack_from_buf(buf, pos, RADIAL_DATA)
+            pos = pos + _structure_size(RADIAL_DATA)
             binnum = radial_data_info['num_bins']
             # rdata = buf[pos:pos+dic_radial_header['bin_length']*radial_data_info['num_bins']]
             rdata = np.frombuffer(buf[pos:pos+dic_radial_header['bin_length']*radial_data_info['num_bins']], '<u2')
@@ -1450,51 +1442,51 @@ class READ_ROSE(object):
         pos = 0
 
         # 获取通用信息头
-        dic_gh = rs._unpack_from_buf(buf, pos, rs.GENERIC_HEADER)
+        dic_gh = _unpack_from_buf(buf, pos, GENERIC_HEADER)
         # pprint(dic_gh)
         if dic_gh['magic_number'] != 0x4D545352:
             print('源数据格式错误！')
-        pos = pos + rs._structure_size(rs.GENERIC_HEADER)
+        pos = pos + _structure_size(GENERIC_HEADER)
 
         # 获取站点信息
-        dic_scfg = rs._unpack_from_buf(buf, pos, rs.SITE_CONFIG)
-        pos = pos + rs._structure_size(rs.SITE_CONFIG)
+        dic_scfg = _unpack_from_buf(buf, pos, SITE_CONFIG)
+        pos = pos + _structure_size(SITE_CONFIG)
 
         # 获取任务信息
-        dic_tcfg = rs._unpack_from_buf(buf, pos, rs.TASK_CONFIG)
-        pos = pos + rs._structure_size(rs.TASK_CONFIG)
+        dic_tcfg = _unpack_from_buf(buf, pos, TASK_CONFIG)
+        pos = pos + _structure_size(TASK_CONFIG)
 
         # 获取扫描信息
         cutinfo = []
         for im in np.arange(dic_tcfg['cut_number']):
-            dic_cutcfg = rs._unpack_from_buf(buf, pos, rs.SCAN_CONFIG)
+            dic_cutcfg = _unpack_from_buf(buf, pos, SCAN_CONFIG)
             cutinfo.append(dic_cutcfg)
-            pos = pos + rs._structure_size(rs.SCAN_CONFIG)
+            pos = pos + _structure_size(SCAN_CONFIG)
 
         # 获取产品头信息
-        dic_prod_header = rs._unpack_from_buf(buf, pos, rs.PRODUCT_HEADER_BLOCK)
-        pos = pos + rs._structure_size(rs.PRODUCT_HEADER_BLOCK)
+        dic_prod_header = _unpack_from_buf(buf, pos, PRODUCT_HEADER_BLOCK)
+        pos = pos + _structure_size(PRODUCT_HEADER_BLOCK)
 
         # 获取产品参数信息
-        prod_type = rs.PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
+        prod_type = PRODUCT_ID_NAME_MAP[dic_prod_header['product_type']]
         # pprint(prod_type)
-        dic_prod_param = rs._unpack_from_buf(buf, pos, rs.PRODUCT_DEPENDENT_PARAMETER[prod_type])
+        dic_prod_param = _unpack_from_buf(buf, pos, PRODUCT_DEPENDENT_PARAMETER[prod_type])
         pos = pos + 64 # 产品参数长度固定为64个
 
         # 获取产品数据块
         # 冰雹个数
-        hail_number = rs._unpack_from_buf(buf, pos, rs.HAIL_NUMBER)
-        pos = pos + rs._structure_size(rs.HAIL_NUMBER)
+        hail_number = _unpack_from_buf(buf, pos, HAIL_NUMBER)
+        pos = pos + _structure_size(HAIL_NUMBER)
 
         # TVS表块
         hail_tab=[]
         for nn in np.arange(hail_number['hail_number']):
-            hail_tab.append(rs._unpack_from_buf(buf, pos, rs.HAIL_TAB))
-            pos = pos + rs._structure_size(rs.HAIL_TAB)
+            hail_tab.append(_unpack_from_buf(buf, pos, HAIL_TAB))
+            pos = pos + _structure_size(HAIL_TAB)
 
         # TVS适配数据
-        dic_hail_adapt_param = rs._unpack_from_buf(buf, pos, rs.HAIL_ADAPTATION_DATA)
-        pos = pos + rs._structure_size(rs.HAIL_ADAPTATION_DATA)
+        dic_hail_adapt_param = _unpack_from_buf(buf, pos, HAIL_ADAPTATION_DATA)
+        pos = pos + _structure_size(HAIL_ADAPTATION_DATA)
 
         # 将所有的方位角和距离转换成经纬度
         for nn in np.arange(hail_number['hail_number']):
@@ -1514,440 +1506,10 @@ class READ_ROSE(object):
         self.hailinfo = allresult
         return allresult
     
-    # 将风暴追踪产品转换为gr2格式
-    def trans_sti_gr2(self,outpath_gr2=None,outname_gr2=None,m_pret=None,time_type='UTC',bhistory=False,bdebug=False):
-        
-        if self.stiinfo is None:
-            print('mesoinfo is None, please run read_sti first!')
-            return False
-
-        # try:
-        # duration 代表该时次产品的可视时间段
-        if not os.path.exists(outpath_gr2):
-            os.makedirs(outpath_gr2)
-
-
-        if outpath_gr2 is None or outname_gr2 is None:
-            print('outpath_gr2 and outname_gr2 params should be set')
-            return False
-
-        if not os.path.exists(outpath_gr2):
-            try:
-                os.makedirs(outpath_gr2)
-            except:
-                print(outpath_gr2 + ' not exist, and fail to create this path!')
-                return False
-
-        of = open(outpath_gr2 + os.sep + outname_gr2,'wt',encoding='GBK')
-
-        # 单个文件就不用信息头了
-        if bdebug and not bhistory:
-            of.write('Title: ROSE2风暴追踪信息\n')
-            of.write('Refreshseconds: 30\n')
-            of.write('Threshold: 500\n')
-            of.write('Font: 1, 20, 1, "Courier New"\n')
-            # of.write('IconFile:  1, 18, 18, 9, 9, 2016_144dpi.png\n')
-            of.write('IconFile:  1, 24, 24, 12, 12, 2016_192dpi.png\n')
-
-        # 这里直接用m_pret 可以保证不会出现重叠的记录
-        if not bdebug:
-            if not m_pret is None and bhistory:
-                pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                endtime = m_pret + timedelta(minutes=1*int(6/2))
-            elif not bhistory:
-                pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                endtime = m_pret + timedelta(minutes=30)
-                pass
-            
-            
-            timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-
-            # 实时模式下，改时间一定要放到timerange后面
-            if time_type == "BJT":
-                pre_time = pre_time + timedelta(hours=8)
-                endtime = endtime + timedelta(hours=8)
-
-            if bhistory:
-                timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-                
-            of.write(timerangestr + '\n')
-
-        for dd in self.stiinfo['track']:
-            if len(dd)  > 1:
-                of.write('Line: 2,0\n')     
-                for nn in dd:
-                    of.write('%.4f,%.4f\n'%(nn[0],nn[1])) 
-                of.write('End:\n\n')
-
-        for dd in self.stiinfo['marker_past']:
-            of.write('Icon: %.4f,%.4f,0,1,88\n'%(dd[0],dd[1]))
-        for dd in self.stiinfo['marker_fst']:
-            of.write('Icon: %.4f,%.4f,0,1,91\n'%(dd[0],dd[1]))
-
-        for nn in np.arange(len(self.stiinfo['all_storm_prop'])):
-
-            curstorm = self.stiinfo['all_storm_prop'][nn]
-            clon = curstorm['lon']
-            clat = curstorm['lat']
-
-            lbstr=r'数据来源: ROSE风暴识别产品\n'
-            lbstr += r'风暴编号: %s\n'%curstorm['id_char']
-            lbstr += r'方位角/距离: %.1f度/%.1f公里\n'%(curstorm['azi'],curstorm['range']/1000.0)
-            lbstr += r'风暴移动方向 : %s\n'%(self.get_wind_dir_name((curstorm['mv_dir']+180)%360))
-            lbstr += r'风暴移动速度 : %.1f 公里/小时\n'%(curstorm['mv_spd']*3.6)
-            
-            lbstr += r'风暴低高: %.1f 公里\n'%(curstorm['base_height']/1000)
-            lbstr += r'风暴顶高: %.1f 公里\n'%(curstorm['top_height']/1000)
-            lbstr += r'风暴体VIL: %.1f 千克/平方米\n'%(curstorm['vil'])
-            lbstr += r'最大反射率: %d dBZ\n'%(curstorm['max_ref'])
-            lbstr += r'最大反射率高度 : %.1f 公里\n'%(curstorm['max_ref_height']/1000)        
-            
-            of.write('Icon: %.4f,%.4f,0,1,40,"%s"\n'%(clat,clon,lbstr))
-
-            
-        of.write('End:\n\n')
-
-        of.close()
-        print(outpath_gr2 + os.sep + outname_gr2 + ' saved!')
-        return outpath_gr2 + os.sep + outname_gr2
-        # except:
-        #     return False
-    
-
-    # 将中气旋的信息转换为gr2格式
-    def trans_mda_gr2(self,outpath_gr2=None,outname_gr2=None,m_pret=None,time_type='UTC',bhistory=False,bdebug=False):
-        if self.mesoinfo is None:
-            print('mesoinfo is None, please run read_mda first!')
-            return False
-        try:
-            # duration 代表该时次产品的可视时间段
-            if not os.path.exists(outpath_gr2):
-                os.makedirs(outpath_gr2)
-
-
-            if outpath_gr2 is None or outname_gr2 is None:
-                print('outpath_gr2 and outname_gr2 params should be set')
-                return False
-
-            if not os.path.exists(outpath_gr2):
-                try:
-                    os.makedirs(outpath_gr2)
-                except:
-                    print(outpath_gr2 + ' not exist, and fail to create this path!')
-                    return False
-
-            of = open(outpath_gr2 + os.sep + outname_gr2,'wt',encoding='GBK')
-
-            # 单个文件就不用信息头了
-            if bdebug and not bhistory:
-                of.write('Title: ROSE2中气旋识别产品\n')
-                of.write('Refreshseconds: 30\n')
-                of.write('Threshold: 500\n')
-                of.write('Font: 1, 20, 1, "Courier New"\n')
-                # of.write('IconFile:  1, 47, 47, 23, 23, mda_icons_144.png\n')
-                of.write('IconFile:  1, 63, 63, 31, 31, mda_icons_192.png\n')
-
-            # 这里直接用m_pret 可以保证不会出现重叠的记录
-            if not bdebug:
-                if not m_pret is None and bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=1*int(6/2))
-                elif not bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=30)
-                    pass
-                
-                
-                timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-
-                # 实时模式下，改时间一定要放到timerange后面
-                if time_type == "BJT":
-                    pre_time = pre_time + timedelta(hours=8)
-                    endtime = endtime + timedelta(hours=8)
-
-                if bhistory:
-                    timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-                    
-                of.write(timerangestr + '\n')
-
-            for nn in range(len(self.mesoinfo['meso'])):
-                curitem = self.mesoinfo['meso'][nn]
-                clon = curitem['lon']
-                clat = curitem['lat']
-                lbstr=''
-                icon_index=0
-                baseindex=0
-                if curitem['base'] > 1500: # 米
-                    baseindex = 0
-                else:
-                    baseindex = 1*4
-                '''
-                关于风切变的单位换算问题
-                https://www.theweatherprediction.com/habyhints/193/
-                The wind speed is 20 knots at 850 mb while the wind speed is 30 knots at 700 mb. The distance between 850 and 700 mb is 1,600 meters. What is the value of wind shear?
-
-                The change in wind speed is 30 - 20 = 10 knots = 5 m/s
-                The change is distance is 1,600 m
-
-                The change in wind speed over distance = 5/1,600 s^-1 = 0.003125 s^-1
-
-                '''
-                shear_ms = curitem['ave_shear'] * (curitem['azi_diam'] + curitem['radial_diam'])/2/1000.0/2
-                mrank = self.get_mda_rank(curitem['range']/1000,shear_ms)
-
-                if mrank is None:
-                    continue
-                else:
-                    icon_index = baseindex + mrank
-
-                mlevel =''
-                if mrank==1:
-                    mlevel = '弱切变'
-                elif mrank == 2:
-                    mlevel = '弱中气旋'
-                elif mrank == 3:
-                    mlevel = '中等强度中气旋'
-                elif mrank == 4:
-                    mlevel = '强中气旋'
-
-                mtype='未知'
-                if curitem['feature_type'] == 1:
-                    mtype = '中气旋'
-                elif curitem['feature_type'] == 2:
-                    mtype = '3D相关切变'
-                elif curitem['feature_type'] == 3:
-                    mtype = '非相关切变'
-
-                lbstr=''
-                lbstr += r'数据来源: ROSE2中气旋识别产品\n'
-                lbstr += r'所属风暴编号: %s\n'%curitem['storm_id_char']
-                lbstr += r'中气旋编号: %s\n'%curitem['feature_id_char']
-                lbstr += r'类型: %s\n'%mtype
-                lbstr += r'中气旋所在方位/距离: %.1f度 / %.1f 公里\n'%(curitem['azi'],curitem['range']/1000.0)
-                lbstr += r'强度等级: %s\n'%mlevel
-                lbstr += r'中气旋底高: %.1f 公里\n'%(curitem['base']/1000.0)
-                lbstr += r'中气旋顶高: %.1f 公里\n'%(curitem['top']/1000.0)
-                lbstr += r'中气旋沿径向直径: %.1f 公里\n'%(curitem['radial_diam']/1000.0)
-                lbstr += r'中气旋沿方位直径: %.1f 公里\n'%(curitem['azi_diam']/1000.0)
-                lbstr += r'平均切变值: %.1f (E-3/S)\n'%(curitem['ave_shear'])
-
-                if shear_ms > 25:
-                    lbstr += r'===============提醒============\n换算后的切变值超过了25米/秒\n可能是很强的中气旋或者是错误数据\n可能是存在速度模糊对PUP算法的影响。\n'
-
-                of.write('Icon: %.4f,%.4f,0,1,%d,"%s"\n'%(clat,clon,icon_index, lbstr))
-                pass
-            of.write('End:\n\n')
-
-            of.close()
-            print(outpath_gr2 + os.sep + outname_gr2 + ' saved!')
-            return outpath_gr2 + os.sep + outname_gr2
-        except:
-            return False
-        
-
-    # 将tvs的信息转换为gr2格式
-    def trans_tvs_gr2(self,outpath_gr2=None,outname_gr2=None,m_pret=None,time_type='UTC',bhistory=False,bdebug=False):
-        
-        if self.tvsinfo is None:
-            print('tvsinfo is None, please run read_tvs  first')
-            return False
-        
-        try:
-            # duration 代表该时次产品的可视时间段
-            if not os.path.exists(outpath_gr2):
-                os.makedirs(outpath_gr2)
-
-
-            if outpath_gr2 is None or outname_gr2 is None:
-                print('outpath_gr2 and outname_gr2 params should be set')
-                return False
-
-            if not os.path.exists(outpath_gr2):
-                try:
-                    os.makedirs(outpath_gr2)
-                except:
-                    print(outpath_gr2 + ' not exist, and fail to create this path!')
-                    return False
-
-            of = open(outpath_gr2 + os.sep + outname_gr2,'wt',encoding='GBK')
-
-            # 单个文件就不用信息头了
-            if bdebug and not bhistory:
-                of.write('Title: ROSE2龙卷涡旋特征识别产品\n')
-                of.write('Refreshseconds: 30\n')
-                of.write('Threshold: 500\n')
-                of.write('Font: 1, 20, 1, "Courier New"\n')
-                of.write('IconFile:  1, 51, 51, 26, 26, lsr_icons_192.png\n')
-
-            # 这里直接用m_pret 可以保证不会出现重叠的记录
-            if not bdebug:
-                if not m_pret is None and bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=1*int(6/2))
-                elif not bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=30)
-                    pass
-                
-                
-                timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-
-                # 实时模式下，改时间一定要放到timerange后面
-                if time_type == "BJT":
-                    pre_time = pre_time + timedelta(hours=8)
-                    endtime = endtime + timedelta(hours=8)
-
-                if bhistory:
-                    timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-                    
-                of.write(timerangestr + '\n')
-
-            for nn in range(len(self.tvsinfo['tvs'])):
-                curitem = self.tvsinfo['tvs'][nn]
-                clon = curitem['lon']
-                clat = curitem['lat']
-                curtype = '未知'
-                if curitem['type'] == 1:
-                    curtype = 'TVS'
-                elif curitem['type'] == 2:
-                    curtype = 'ETVS'
-                icon_index=1
-                
-                lbstr=''
-                lbstr += r'数据来源: ROSE2龙卷涡旋识别产品\n'
-                lbstr += r'所属风暴编号: %s\n'%curitem['storm_id_char']
-                lbstr += r'类型: %s\n'%curtype
-                lbstr += r'TVS所在方位/距离: %.1f度 / %.1f 公里\n'%(curitem['azi'],curitem['range']/1000.0)
-
-                lbstr += r'TVS底高: %.1f 公里\n'%(curitem['base']/1000.0)
-                lbstr += r'TVS顶高: %.1f 公里\n'%(curitem['top']/1000.0)
-                lbstr += r'TVS厚度: %.1f 公里\n'%(curitem['depth']/1000.0)
-                lbstr += r'低层速度差LLDV: %.1f 米/秒\n'%(curitem['lldv']) # *1.609344/3.6
-                lbstr += r'最大速度差MXDV: %.1f 米/秒\n'%(curitem['mxdv']) # *1.609344/3.6
-                lbstr += r'最大速度差所在高度: %.1f 公里\n'%(curitem['hmdv']/1000.0)
-
-                lbstr += r'最大切变值MASHR: %.1f (E-3/S)\n'%(curitem['max_shear'])
-                lbstr += r'最大切变值所在高度: %.1f 公里\n'%(curitem['h_max_shear']/1000.0)
-            
-                if curtype=='ETVS':
-                    icon_index=20
-                    of.write('Icon: %.4f,%.4f,0,1,%d,"%s"\n'%(clat,clon,icon_index, lbstr))
-                else:
-                    icon_index=11
-                    of.write('Icon: %.4f,%.4f,0,1,%d,"%s"\n'%(clat,clon,icon_index, lbstr))
-
-            of.write('End:\n\n')
-
-            of.close()
-            print(outpath_gr2 + os.sep + outname_gr2 + ' saved!')
-            return outpath_gr2 + os.sep + outname_gr2
-        except:
-            return False
-        
-
-    # 将hail的信息转换为gr2格式
-    def trans_hda_gr2(self,outpath_gr2=None,outname_gr2=None,m_pret=None,time_type='UTC',bhistory=False,bdebug=False):
-
-        if self.hailinfo is None:
-            print('hailinfo is None, please run read_hda first!')
-            return False
-
-        try:
-            # duration 代表该时次产品的可视时间段
-            if not os.path.exists(outpath_gr2):
-                os.makedirs(outpath_gr2)
-
-
-            if outpath_gr2 is None or outname_gr2 is None:
-                print('outpath_gr2 and outname_gr2 params should be set')
-                return False
-
-            if not os.path.exists(outpath_gr2):
-                try:
-                    os.makedirs(outpath_gr2)
-                except:
-                    print(outpath_gr2 + ' not exist, and fail to create this path!')
-                    return False
-
-            of = open(outpath_gr2 + os.sep + outname_gr2,'wt',encoding='GBK')
-
-            # 单个文件就不用信息头了
-            if bdebug and not bhistory:
-                of.write('Title: ROSE2冰雹指数产品\n')
-                of.write('Refreshseconds: 30\n')
-                of.write('Threshold: 500\n')
-                of.write('Font: 1, 20, 1, "Courier New"\n')
-                # of.write('IconFile:  1, 37, 37, 18, 18, hail_icons_144.png\n')
-                of.write('IconFile:  1, 51, 51, 26, 26, hail_icons_192.png\n')
-
-            # 这里直接用m_pret 可以保证不会出现重叠的记录
-            if not bdebug:
-                if not m_pret is None and bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=1*int(6/2))
-                elif not bhistory:
-                    pre_time = m_pret + timedelta(minutes=-1*int(6/2))
-                    endtime = m_pret + timedelta(minutes=30)
-                    pass
-                
-                
-                timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-
-                # 实时模式下，改时间一定要放到timerange后面
-                if time_type == "BJT":
-                    pre_time = pre_time + timedelta(hours=8)
-                    endtime = endtime + timedelta(hours=8)
-
-                if bhistory:
-                    timerangestr = 'TimeRange: ' + pre_time.strftime('%Y-%m-%dT%H:%M:%S') + ' ' + endtime.strftime('%Y-%m-%dT%H:%M:%S')
-                    
-                of.write(timerangestr + '\n')
-
-            for nn in range(len(self.hailinfo['hail'])):
-                curitem = self.hailinfo['hail'][nn]
-                clon = curitem['lon']
-                clat = curitem['lat']
-                
-                lbstr=''
-                lbstr=r'数据来源: ROSE2冰雹指数产品\n'
-                lbstr += r'冰雹编号: %s\n'%curitem['hail_id_char']
-                lbstr += r'大冰雹概率: %d%%\n'%curitem['psh']
-                lbstr += r'冰雹概率: %d%%\n'%curitem['ph']
-                lbstr += r'最大冰雹尺寸估计: %.1f 厘米 或 %.1f 英寸\n'%(curitem['hsize'],curitem['hsize']/2.54,)
-                
-                icon_index=0
-                hsize_inch = curitem['hsize']/2.54
-                
-                if hsize_inch < 1:
-                    icon_index = 1
-                elif hsize_inch < 2:
-                    icon_index = 5
-                elif hsize_inch < 3:
-                    icon_index = 9
-                elif hsize_inch < 4:
-                    icon_index = 13
-                elif hsize_inch >= 4:
-                    icon_index = 17
-                else:
-                    icon_index = 1
-                if curitem['ph'] >= 30:
-                    of.write('Icon: %.4f,%.4f,0,1,%d,"%s"\n'%(clat,clon,icon_index, lbstr))
-
-            of.write('End:\n\n')
-
-            of.close()
-            print(outpath_gr2 + os.sep + outname_gr2 + ' saved!')
-            return outpath_gr2 + os.sep + outname_gr2
-        except:
-            return False
         
 if __name__ == "__main__":
-    rootpath = '/Users/wenjianzhu/Library/CloudStorage/OneDrive-个人/PythonCode/MyWork/nmc_met_radar/testdata/rose_product'
+    rootpath = 'testdata/rose_product'
     staid = 'Z9852'
-    # filepath = 'D:\Downloads\CMADAAS\pup\Z9852\STI'
-    outpath_gr2 = 'output/gr2data/sti' + os.sep + staid
-    if not os.path.exists(outpath_gr2):
-        os.makedirs(outpath_gr2)
 
     decoder = READ_ROSE()
 
@@ -1961,33 +1523,28 @@ if __name__ == "__main__":
     
 
     # sti
-    filepath_sti = rootpath + os.sep + staid + os.sep + 'STI'
-    filename_sti = 'Z_RADR_I_Z9852_20230421164954_P_DOR_CDD_STI_NUL_300_NUL_FMT.bin'
-    allresult_sti = decoder.read_sti(filepath_sti, filename_sti)
-    outname_sti = 'pls_sti.txt'
-    decoder.trans_sti_gr2(outpath_gr2,outname_sti,time_type='BJT',bdebug = True)
+    # filepath_sti = rootpath + os.sep + staid + os.sep + 'STI'
+    # filename_sti = 'Z_RADR_I_Z9852_20230421164954_P_DOR_CDD_STI_NUL_300_NUL_FMT.bin'
+    # allresult_sti = decoder.read_sti(filepath_sti, filename_sti)
     
     
-    # meso
-    filepath_meso = rootpath + os.sep + staid + os.sep + 'M'
-    filename_meso = 'Z_RADR_I_Z9852_20230421135042_P_DOR_CDD_M_NUL_200_NUL_FMT.bin'
-    allresult_meso = decoder.read_mda(filepath_meso, filename_meso)
-    outname_meso = 'pls_meso.txt'
-    decoder.trans_mda_gr2(outpath_gr2,outname_meso,time_type='BJT',bdebug = True)
-
-    # tvs
-    filepath_tvs = rootpath + os.sep + staid + os.sep + 'TVS'
-    filename_tvs = 'Z_RADR_I_Z9852_20230421201451_P_DOR_CDD_TVS_NUL_200_NUL_FMT.bin'
-    allresult_tvs = decoder.read_tvs(filepath_tvs, filename_tvs)
-    outname_tvs = 'pls_tvs.txt'
-    decoder.trans_tvs_gr2(outpath_gr2,outname_tvs,time_type='BJT',bdebug = True)
-
-
-    # hail
-    filepath_hail = rootpath + os.sep + staid + os.sep + 'HI'
-    filename_hail = 'Z_RADR_I_Z9852_20230421181642_P_DOR_CDD_HI_NUL_200_NUL_FMT.bin'
-    allresult_hail = decoder.read_hda(filepath_hail, filename_hail)
-    outname_hail = 'pls_hail.txt'
-    decoder.trans_hda_gr2(outpath_gr2,outname_hail,time_type='BJT',bdebug = True)
+    # # meso
+    # filepath_meso = rootpath + os.sep + staid + os.sep + 'M'
+    # filename_meso = 'Z_RADR_I_Z9852_20230421135042_P_DOR_CDD_M_NUL_200_NUL_FMT.bin'
+    # allresult_meso = decoder.read_mda(filepath_meso, filename_meso)
+   
+
+    # # tvs
+    # filepath_tvs = rootpath + os.sep + staid + os.sep + 'TVS'
+    # filename_tvs = 'Z_RADR_I_Z9852_20230421201451_P_DOR_CDD_TVS_NUL_200_NUL_FMT.bin'
+    # allresult_tvs = decoder.read_tvs(filepath_tvs, filename_tvs)
+   
+
+
+    # # hail
+    # filepath_hail = rootpath + os.sep + staid + os.sep + 'HI'
+    # filename_hail = 'Z_RADR_I_Z9852_20230421181642_P_DOR_CDD_HI_NUL_200_NUL_FMT.bin'
+    # allresult_hail = decoder.read_hda(filepath_hail, filename_hail)
+   
 
 # %%