metradar 0.1.0__py3-none-any.whl

metradar/trans_nc_pgmb.py

@@ -0,0 +1,62 @@
+'''
+测试将雷达拼图输出成pgmb格式
+
+'''
+
+# %%
+from pgmb_io import pgmb_write
+import os
+from read_new_mosaic import decode_mosaic
+import xarray as xr
+import numpy as np
+import gzip
+
+def zip_file(path, new_path):
+    with open(new_path, 'wb') as wf:
+        with open(path, 'rb') as rf:
+            data = rf.read()
+            # 压缩数据
+            data_comp = gzip.compress(data,compresslevel=5)
+        wf.write(data_comp)
+ 
+ 
+filepath = '/Users/wenjianzhu/Downloads/CMADAAS/rdmosaic_bin/202307/'
+for filename in os.listdir(filepath):
+    if not filename.endswith('.nc'):
+        continue
+    # filename = 'ACHN.CREF000.20230421.180009.nc'
+    obstimestr = filename.split('.')[2] + filename.split('.')[3][0:4]
+    outname = filename.split('.')[2] + filename.split('.')[3][0:4] + '_fmi.radar.composite.lowest_FIN_SUOMI1.pgm'
+    outpath = '/Users/wenjianzhu/pysteps_data/radar/fmi/pgm/%s'%filename.split('.')[2] 
+    if not os.path.exists(outpath):
+        os.makedirs(outpath)
+
+    file_name = outpath + os.sep + outname
+
+    org_data = xr.open_dataset(filepath+os.sep+filename)
+    data = org_data.sel(lat=slice(33, 37), lon=slice(109, 117))
+
+    width = data.cref.shape[1]
+    height = data.cref.shape[0]
+    maxval = 255
+    gray = data.cref*2+66
+
+    aa = np.isnan(gray.data)
+    gray.data[aa]=255
+    cref = np.flipud(gray.data.astype(np.uint8))
+    params={}
+    params['obstimestr'] = obstimestr
+    params['left_lon'] = 109
+    params['right_lon'] = 117
+    params['bottom_lat'] = 33
+    params['upper_lat'] = 37
+
+    pgmb_write ( file_name, params, width, height, maxval, cref )
+
+    # gzip file
+
+    # zip_file(file_name, file_name + '.gz')
+
+
+    print('write pgmb file %s successfully!'%(file_name))
+    # break

metradar/trans_new_mosaic_nc.py

@@ -0,0 +1,309 @@
+'''
+读取新的雷达拼图数据
+朱文剑
+20210830
+
+数据格式说明
+
+//
+/////////////////////////////////////////////////////////////////////
+//  组网定量估测降水系统拼图产品文件规则
+//  1. 文件里命名为:
+//		ACHN_XXX_20121012_090000_nn.bin
+// 2. 文件格式:文件头 + 数据块 
+// 3. 数据从西北，往东南写
+//***************************************
+typedef struct tagMQPE_MosiacProductHeaderStructure
+{
+	// 文件卷标，12字节(0-11), 
+	char  label[4];	    //0-3	文件固定标识：MOC 
+	char  Version[4];	//4-7	文件格式版本代码，如: 1.0，1.1，etc
+	int   FileBytes;	//8-11	包含头信息在内的文件字节数，不超过2M
+
+						//产品描述, 84字节(12-95) 
+	short MosaicID;		  //12-13 拼图产品编号
+	short coordinate;	  //14-15 坐标类型: 2=笛卡儿坐标,3=等经纬网格坐标
+	char  varname[8];     //16-23 产品代码,如: ET,VIL,CR,CAP,OHP,OHPC
+	char  description[64];//24-87 产品描述,如Composite Reflectivity mosaic
+	int   BlockPos; //88-91 产品数据起始位置(字节顺序)
+	int   BlockLen; //92-95  产品数据字节数
+
+					//数据时间, 28字节(96-123) 
+	int TimeZone;   //96-99 数据时钟,0=世界时,28800=北京时
+	short yr;		//100-101	观测时间中的年份
+	short mon;	    //102-103	观测时间中的月份（1－12）
+	short day;		//104-105	观测时间中的日期（1－31）
+	short hr;		//106-107	观测时间中的小时（00－23）
+	short min;	    //108-109	观测时间中的分（00－59）
+	short sec;	    //110-111	观测时间中的秒（00－59）
+	int   ObsSeconds; //112-115	观测时间的seconds 
+	unsigned short ObsDates; //116-117 观测时间中的Julian dates
+	unsigned short GenDates; //118-119 产品处理时间的天数
+	int   GenSeconds;   //120-123 产品处理时间的描述
+
+						//数据区信息, 48字节(124-171) 
+	int  edge_s;  //124-127 数据区的南边界，单位：1/1000度，放大1千倍
+	int  edge_w;  //128-131 数据区的西边界，单位：1/1000度，放大1千倍
+	int  edge_n;  //132-135 数据区的北边界，单位：1/1000度，放大1千倍
+	int  edge_e;  //136-139 数据区的东边界，单位：1/1000度，放大1千倍
+	int  cx;      //140-143 数据区中心坐标，单位：1/1000度，放大1千倍
+	int  cy;      //144-147 数据区中心坐标，单位：1/1000度，放大1千倍
+	int  nX;      //148-151 格点坐标为列数,
+	int  nY;      //152-155 格点坐标为行数,
+	int  dx;      //156-159 格点坐标为列分辨率，单位：1/10000度，放大1万倍
+	int  dy;      //160-163 格点坐标为行分辨率，单位：1/10000度，放大1万倍
+	short height; //164-165 雷达高度
+	short Compress; //166-167 数据压缩标识, 0=无,1=bz2,2=zip,3=lzw
+	int   num_of_radars; //168-171有多少个雷达进行了拼图
+
+	int UnZipBytes; //数据段压缩前的字节数, //172-175
+	short scale;     //176-177
+	short unUsed; //178-179
+	char  RgnID[8];  //180-187
+	char units[8];  //180-187
+	char  reserved[60];
+}MOSAIC_PRODUCTHEADER;
+
+-32768表示空白
+
+
+'''
+
+# %%
+CODE1 = 'B'
+CODE2 = 'H'
+INT1 = 'B'
+INT2 = 'H'
+INT4 = 'I'
+REAL4 = 'f'
+REAL8 = 'd'
+SINT1 = 'b'
+SINT2 = 'h'
+SINT4 = 'i'
+
+MOSAIC_HEADER = (
+    ('label', '4s'),
+    ('Version', '4s'),
+    ('FileBytes', INT4),
+    ('MosaicID', INT2),
+    ('coordinate', INT2),
+    ('varname', '8s'),
+	('description', '64s'),
+    ('BlockPos', INT4),
+    ('BlockLen', INT4),
+    ('TimeZone', INT4),
+    ('year', INT2),
+    ('month', INT2),
+	('day', INT2),
+    ('hour', INT2),
+    ('min', INT2),
+    ('sec', INT2),
+    ('ObsSeconds', INT4),
+    ('ObsDates', INT2),
+	('GenDates', INT2),
+    ('GenSeconds', INT4),
+    ('edge_s', INT4),
+    ('edge_w', INT4),
+    ('edge_n', INT4),
+    ('edge_e', INT4),
+	('cx', INT4),
+    ('cy', INT4),
+    ('nX', INT4),
+    ('nY', INT4),
+    ('dx', INT4),
+    ('dy', INT4),
+	('height', INT2),
+    ('Compress', INT2),
+    ('num_of_radars', INT4),
+    ('UnZipBytes', INT4),
+    ('scale', INT2),
+    ('unUsed', INT2),
+	('RgnID', '8s'),
+	('units', '8s'),
+	('reserved', '60s'),
+
+)
+
+import os
+import numpy as np
+import struct
+import bz2
+import xarray as xr
+from datetime import datetime,timedelta
+
+
+from multiprocessing import cpu_count, Pool,freeze_support
+
+def _structure_size(structure):
+        """ Find the size of a structure in bytes. """
+        return struct.calcsize('<' + ''.join([i[1] for i in structure]))
+
+
+def _unpack_from_buf( buf, pos, structure):
+    """ Unpack a structure from a buffer. """
+    size = _structure_size(structure)
+    return _unpack_structure(buf[pos:pos + size], structure)
+
+
+def _unpack_structure(string, structure):
+    """ Unpack a structure from a string """
+    fmt = '<' + ''.join([i[1] for i in structure])  # little-endian
+    lst = struct.unpack(fmt, string)
+    return dict(zip([i[0] for i in structure], lst))
+
+def decode_mosaic(filepath,filename,minv=10,maxv = 80):
+
+	if not os.path.exists(filepath+os.sep + filename):
+		print(filepath+os.sep+filename + ' not exists!')
+		return None,None
+	try:
+		fin = open(filepath + os.sep + filename,'rb')
+	except:
+		return None,None
+
+	buf = fin.read()
+	fin.close()
+	pos = 0
+	if len(buf)==0:
+		return None
+	dic_gh = _unpack_from_buf(buf, pos, MOSAIC_HEADER)
+	pos +=_structure_size(MOSAIC_HEADER)
+	databuf = bz2.decompress(buf[pos:])
+	data = np.frombuffer(databuf, 'int16')
+	data = np.reshape(data,[dic_gh['nY'],dic_gh['nX']])
+	data=data.astype('float')
+	# data[data == -32768] = np.NAN #  空白
+	# data[data == -1280] = np.NAN # 有覆盖但无数据
+	data[data < -300] = np.NAN
+	data = np.flipud(data)/dic_gh['scale']
+	data[data < int(minv)] = np.NAN
+	data[data > int(maxv)] = np.NAN
+
+	# set longitude and latitude coordinates
+	lat = dic_gh['edge_s']/1000 + np.arange(dic_gh['nY'])*dic_gh['dy']/10000
+	lon = dic_gh['edge_w']/1000 + np.arange(dic_gh['nX'])*dic_gh['dx']/10000
+
+	# set time coordinates
+	time = datetime(dic_gh['year'],dic_gh['month'],dic_gh['day'],dic_gh['hour'],dic_gh['min'],dic_gh['sec'])
+	time = np.array([time], dtype='datetime64[m]')
+	# data = np.expand_dims(data, axis=0)
+	tmpdata = np.reshape(data, (1, data.shape[0], data.shape[1]))
+	# define coordinates
+	lev_coord = ('lev', [dic_gh['height'],], {'long_name':'height', 'units':'m', '_CoordinateAxisType':'Height'})
+	lon_coord = ('lon', lon, {
+	    'long_name':'longitude', 'units':'degrees_east', '_CoordinateAxisType':'Lon'})
+	lat_coord = ('lat', lat, {
+	    'long_name':'latitude', 'units':'degrees_north', '_CoordinateAxisType':'Lat'})
+
+	# create xarray
+	varattrs = {'long_name': 'Refelectivity', 
+	            'short_name': 'ref', 'units': 'dBZ',
+				'maxv':maxv,
+				'minv':minv}
+	data = xr.Dataset({'ref':(['lev','lat', 'lon'], tmpdata, varattrs)},coords={ 'lev':lev_coord,'lat':lat_coord, 'lon':lon_coord})
+
+	# add attributes
+	data.attrs['Conventions'] = "CF-1.6"
+	data.attrs['Origin'] = 'cmadaas'
+	data.attrs['author'] = 'ZWJ'
+	
+	return data
+
+def trans_single(param):
+
+	filepath = param['filepath']
+	filename = param['filename']
+	outpath = param['outpath']
+	pass
+	varflag = 'unknown'
+	if filename.find('CREF') > 0:
+		varflag = 'CREF000'
+	elif filename.find('QREF') > 0:
+		varflag = 'QREF000'
+	elif filename.find('_CAP_') > 0:
+		varflag = 'CAP'
+	if varflag == 'CAP':
+		outname = 'ACHN.' + varflag + '.' + filename.split('_')[4][0:8] +'.' + \
+			       filename.split('_')[4][8:14] + '_' + filename.split('_')[-1].split('.')[0] + '.nc'
+	else:
+		# outname = 'ACHN.' + varflag + '.' + filename.split('_')[4][0:8] +'.' + filename.split('_')[4][8:14] + '.nc'
+		outname = 'ACHN.' + varflag + '.' + filename.split('_')[2] +'.' + filename.split('_')[3][0:6] + '.nc'
+	if os.path.exists(outpath + os.sep + outname):
+		print(outpath + os.sep + outname + ' already exists!')
+		return True
+	try:
+		data=decode_mosaic(filepath,filename,minv=0,maxv=80)
+		if data is None:
+			return False
+		datakey = list(data.keys())
+		if 'cref' not in datakey:
+			#将ref变量rename为cref
+			data=data.rename({datakey[0]:'cref'})
+		if not os.path.exists(outpath):
+			os.makedirs(outpath)
+		data.to_netcdf(outpath + os.sep + outname,encoding={'cref': {'zlib': True}})
+		print(outpath + os.sep + outname + ' done!')
+	except Exception as e:
+		print(e)
+		print(outpath + os.sep + outname + ' make error!')
+
+# 批量转nc文件
+def trans_bin_nc_batch(inpath,outpath,):
+
+	if not os.path.exists(inpath):
+		print(inpath + ' not exists!')
+	if not os.path.exists(outpath):
+		os.makedirs(outpath)
+
+	files = os.listdir(inpath)
+	files=sorted(files)
+	if len(files)==0:
+		print('no valid files')
+		return False
+	
+	params=[]
+	for filename in files:
+		if filename.startswith('.') or filename.startswith('..'):
+			continue
+		if not filename.endswith('.bin'):
+			continue
+		curparam=dict()
+		curparam['filepath'] = inpath
+		curparam['filename'] = filename
+		curparam['outpath'] = outpath# + os.sep + filename.replace('.bin','.nc')#filename.split('_')[9][0:4]
+		params.append(curparam)
+		
+
+	# 多进程处理
+	MAXP = int(cpu_count()*0.6)
+	pools = Pool(MAXP)
+
+	pools.map(trans_single, params)
+	pools.close()
+	pools.join()
+
+if __name__ == '__main__':
+
+
+	sourcepath = r'/Users/wenjianzhu/Downloads/CMADAAS/rdmosaic_bin'
+	outpath = r'/Users/wenjianzhu/Downloads/CMADAAS/mosaicdata_nc'
+	starttime = '20230701'
+	endtime = '20230702'
+	startt = datetime.strptime(starttime,'%Y%m%d')
+	endt = datetime.strptime(endtime,'%Y%m%d')
+	curt = startt
+	# ACHN.CREF000.20221001.143000.nc
+	prepath = ''
+	while curt <= endt:
+		pass
+		curpath = sourcepath + os.sep + curt.strftime('%Y%m')
+		if prepath == curpath:
+			curt  = curt + timedelta(days=1)
+			continue
+		if not os.path.exists(curpath):
+			curt  = curt + timedelta(days=1)
+			continue
+		
+		trans_bin_nc_batch(curpath,outpath)
+		curt  = curt + timedelta(days=1)
+		prepath = curpath

metradar/trans_polor2grid_func.py

@@ -0,0 +1,203 @@
+'''
+本程序将极坐标格式的ppi数据转换为格点数据，并保存为nc格式
+
+'''
+
+# %%
+import numpy as np
+import xarray as xr
+from decode_fmt_pyart import read_cnrad_fmt
+import os
+import time
+import shutil
+
+def do_ppi(params):
+
+    # sweep_slice,radar,outname
+    sweep_slice = params['sweep_slice']
+    target_reso = params['out_reso']
+    # print(sweep_slice)
+    radar = params['radar']
+
+    ori_reso = radar.range['meters_between_gates']
+    # target_reso = 500
+    reso_list = list(np.arange(0,11,1)*ori_reso)
+    # 从reso_list中找到最接近target_reso的值
+    out_scale = np.argmin(np.abs(np.array(reso_list)-target_reso))
+    # out_scale = int(target_reso/ori_reso)
+    target_reso = ori_reso * out_scale
+
+    # sweep_slice = slice(0,384,None)
+    sweep_azimuths = radar.azimuth['data'][sweep_slice]
+    elv = np.mean(radar.elevation['data'][sweep_slice])
+
+    # 输出的距离库只保留一半
+    ngates = int(len(radar.gate_x['data'][sweep_slice][0])/2) #//2*2 
+    grid_x = np.arange(-1*ngates,ngates+1,out_scale)
+    grid_y = np.arange(-1*ngates,ngates+1,out_scale)
+    total_gates = len(grid_x)
+
+    aa = np.meshgrid(grid_x,grid_y)
+    azi_grid = np.arctan2(aa[0],aa[1])*180/np.pi
+    azi_grid[azi_grid<0]+=360
+    
+    azi_reso = 360/len(sweep_azimuths)
+
+    # 求方位角索引
+    new_azi = azi_grid.flatten()
+    t = new_azi-sweep_azimuths[0]
+    t[t<0]+=360
+    ray_number = np.round(t/azi_reso,0).astype(int)
+    ray_number[ray_number==len(sweep_azimuths)]=0
+    ray_number = np.reshape(ray_number,(total_gates,total_gates))
+
+
+    # 求距离索引
+    dis_grid = np.sqrt(aa[0]**2 + aa[1]**2)
+    dis_grid = np.round(dis_grid.flatten(),0).astype(int)
+    dis_grid = np.reshape(dis_grid,(total_gates,total_gates))
+
+  
+  
+    sweepdata = radar.fields['reflectivity']['data'][sweep_slice]
+    # 对数据进行截断，在径向方向上
+    sweepdata = sweepdata[:,0:ngates]
+    # if np.max(sweepdata.data.flatten()) < -30:
+    #     # 表示该层无有效数据
+    #     return None
+    # data_grid = np.zeros((total_gates,total_gates)) - 9999
+    data_grid = np.zeros((total_gates,total_gates))
+    new_data = data_grid.flatten()
+    new_spdata = sweepdata.flatten()
+
+    pos_out = [i+j*total_gates for i in range(total_gates) for j in range(total_gates) if dis_grid[i,j] < ngates]
+    d_out = [rn*ngates+dg for rn,dg in zip(ray_number.flatten(),dis_grid.flatten()) if dg < ngates]
+   
+    new_data[pos_out]=new_spdata[d_out]
+    new_data[new_data<-10]=np.nan
+    new_data[new_data>=95]=np.nan
+    new_data = new_data *2 + 64
+    new_data[np.isnan(new_data)]=255
+    new_data = new_data.astype(np.uint8)
+    data_grid = np.reshape(new_data,(total_gates,total_gates))
+    outdic={'data':data_grid.transpose(),'elv':elv,'tol_gates':total_gates,
+            'out_reso':target_reso,'datatype':'uint8','scale_factor':2,'add_offset':64,
+            'missing_value':255}
+    print('ele: ' + str(elv) + ' done!')
+
+    del data_grid,new_data,new_spdata,pos_out,d_out,sweepdata
+    return outdic
+
+def do_radar(param):
+    pass
+    try:
+        filename = param['filename']
+        outname = param['outname']
+        out_reso = param['out_reso']
+        if 'bcopy' in param.keys():
+            bcopy = param['bcopy']
+        else:
+            bcopy = 0
+        if 'outpath2' in param.keys():
+            outpath2 = param['outpath2']
+        else:
+            outpath2 = ''
+            print('outpath2 not set!')
+        print('start trans ' + filename)
+
+        # 读取雷达数据
+        radar = read_cnrad_fmt(filename)
+
+        if 'reflectivity' in radar.fields.keys():
+            print('ref in data!')
+
+
+        start  = time.time()
+
+        all_data = []
+        elvs = []
+        
+        nnf=0
+        pre_ele=-999
+        for sweep_slice in radar.iter_slice():
+            # if nnf ==2:
+            #     break
+            params={'sweep_slice':sweep_slice,'radar':radar,'out_reso':out_reso,}
+            cur_ele = np.mean(radar.elevation['data'][sweep_slice])
+            if abs(cur_ele - pre_ele) < 0.2:
+                # print('该仰角已处理过！')
+                # 实际上是跳过了最低仰角的两层径向速度
+                continue
+
+            outdic = do_ppi(params)
+            if outdic is None:
+                continue
+            all_data.append(outdic['data'])
+            elvs.append(outdic['elv'])
+            pre_ele = cur_ele
+            nnf +=1
+        end = time.time()
+       
+
+        #构建xarray数据
+        total_gates = outdic['tol_gates']
+        if total_gates % 2 == 0:
+            out_grid = np.arange(int(-total_gates/2),int(total_gates/2))
+        else:
+            out_grid = np.arange(int(-(total_gates-1)/2),int((total_gates-1)/2)+1)
+        
+        if len(all_data) == 0:
+            print('no data in ' + filename)
+
+            return
+        all_data = np.array(all_data)
+        all_data = all_data.astype(np.uint8)
+        data = xr.DataArray(np.array(all_data),coords=[elvs,out_grid,out_grid],dims=['z','x','y'],name='ref')
+        data.attrs['units'] = 'dBZ'
+        data.attrs['standard_name'] = 'equivalent_reflectivity_factor'
+        data.attrs['long_name'] = 'equivalent_reflectivity_factor'
+        data.attrs['radar_lat'] = radar.latitude['data'][0]
+        data.attrs['radar_lon'] = radar.longitude['data'][0]
+        data.attrs['radar_altitude'] = radar.altitude['data'][0]
+        data.attrs['grid_num'] = total_gates
+        data.attrs['grid_reso'] = outdic['out_reso']
+        data.attrs['elevation'] = elvs
+        data.attrs['obs_range'] = int(outdic['out_reso'] * (total_gates-1)/2)
+        data.attrs['distance_unit'] = 'meter'
+        data.attrs['missing_value'] = 255
+        data.attrs['scale_factor'] = 2
+        data.attrs['add_offset'] = 64
+        data.attrs['datatype'] = 'uint8'
+        data.attrs['decode_method'] = 'dbz = (data - add_offset) / scale_factor'
+
+        if not os.path.exists(os.path.dirname(outname)):
+            os.makedirs(os.path.dirname(outname))
+
+        data.to_netcdf(outname,encoding={'ref': {'zlib': True,'dtype':'uint8'},})
+        print(outname + ' saved!')
+      
+        # copy to outpath2
+        if bcopy == 1:
+            if not os.path.exists(outpath2):
+                os.makedirs(outpath2)
+            shutil.copy(outname,outpath2 + os.sep + os.path.basename(outname))
+            print(outpath2 + os.sep + os.path.basename(outname) + ' copied!')
+         
+
+    except Exception as e:
+        print(e)
+    
+if __name__ == "__main__":
+    rootpath = 'testdata/FMT/'
+    filename = 'Z_RADR_I_Z9200_20240831130000_O_DOR_SAD_CAP_FMT.bin.bz2'
+
+    outname = '/Users/wenjianzhu/Downloads/polargrid/' + filename + '_new.nc'
+    params={'filename':rootpath + filename,'outname':outname,'out_reso':500,'bcopy':0}
+
+    
+    st = time.time()
+    do_radar(params)
+    et = time.time()
+    print(et-st)
+
+    

metradar-0.1.0.dist-info/METADATA

@@ -0,0 +1,12 @@
+Metadata-Version: 2.2
+Name: metradar
+Version: 0.1.0
+Summary: radar data processing
+Author-email: Wenjian Zhu <kevin2075@163.com>
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+Requires-Dist: numpy>=1.18
+Requires-Dist: requests>=2.24.0

metradar-0.1.0.dist-info/RECORD

@@ -0,0 +1,41 @@
+metradar/__init__.py,sha256=K7a1zThYS0yYHgWTDHbJ5jcIHZVtm5_RYwkzcKEWi0g,183
+metradar/cnrad_level2.py,sha256=vy4d-CCAigJGZrm-Dr3_NsmzSs1Tg_ufQvVSZSH57ic,44916
+metradar/comm_func.py,sha256=SByR9y5y9qdEcoDooSbEAN0Mxcw8QnPFN78n-s-INbg,4170
+metradar/construct_aws_refvpr_mainprog.py,sha256=CJgEZRt9Bcqj-KkiUnKfdTnz2CMG47rkJljnT6-4bWg,22205
+metradar/construct_aws_refvpr_mainprog_cams.py,sha256=9xCa5_PuqZxPEW-A-657D2-sQ7ZFRQYmlT387mWYWWY,12404
+metradar/construct_aws_refvpr_mainprog_datan3d.py,sha256=1jP4Kg86wfVX3Xal5MBG9T8IgQ5NqKedW9WsvIz7mtE,15857
+metradar/construct_aws_refvpr_mainprog_swan.py,sha256=uP1RqMX5191tRDh9bOKS_JM86Td1M66dr-RkpTepHy8,12175
+metradar/decode_fmt_pyart.py,sha256=Fjl6tSr53AUnir_HeDEOD0zQ8kuA5kPzWo-77zaWsiY,7542
+metradar/decode_pup_rose.py,sha256=d5ajuJuM6X8kCiNGUBOrCeub6j2K9Y-Qgvjq4H-Qemc,82826
+metradar/draw_mosaic_new.py,sha256=1PyeNGa0Rpj_u-veXorer_bUeP7ogLMCpNbEcooKhus,16193
+metradar/draw_radar_aws_jilin_new.py,sha256=57RIy1u-1ENv_iJrSOxVLMI5Fcf4Erg5bWben4bDxpE,7904
+metradar/draw_radar_comp_func.py,sha256=Hc7TtZMLobW4VCPYPlYyq5MarVYBYDqgBhVJrYMgcHI,63969
+metradar/exceptions.py,sha256=C5jsADZix3E6Ah1sUPyBtgrGygy69lFc0QD5wO07a00,1300
+metradar/geo_transforms_pyart.py,sha256=7yOJOVLUaBnUQiesgvKZuESSNOx1uEYDilcRQAw5avM,21941
+metradar/get_cross_section_from_pyart.py,sha256=x4fSAUa8AzhuXsn48lvduA1ELlL-L7wZQZq8r4FniyU,11662
+metradar/get_tlogp_from_sharppy.py,sha256=oYGFOHjkSbeRYS4yFXk9APfTsuqN3CxUVuv5xR6TUAs,2832
+metradar/grid.py,sha256=2ARummd819SCjWi87toZgOZ0A8obcIUFfAmRV-11EU0,11165
+metradar/grid_data.py,sha256=9e4-x_82ZwcwtjRM6BWUGWQXfg1PhuAH3F59OAd8wxo,1878
+metradar/main_pydda.py,sha256=i32-YF0ebm4im7UkQcUXKcz0TJvRscunzdxSsz56gqU,28659
+metradar/make_gif.py,sha256=rDNqhZXKAyQB_B2adU-ven07QgOTXwdY6eRFVsjTrdI,563
+metradar/make_mosaic_mp_archive.py,sha256=AXiTSSOL1ik0aLy-u4B_LjO4-sPVinhBqkTtPcMVlYI,19556
+metradar/mosaic_merge.py,sha256=eTTBJ5w2LchMBhjMHp1g97PGVcSzadoSKQwgtqjE4-M,1628
+metradar/mosaic_quickdraw.py,sha256=k0WI5hUKifCNSlbMrS8oJSBuM5XjWu5QDOkXoR50wFI,11910
+metradar/nowcast_by_pysteps.py,sha256=OcDY5Gb0pd1EM4pPD-bhIpoln3bH3aZ74J4DTCa4INk,7250
+metradar/oa_couhua.py,sha256=Gs85Lf1s5zuVdUQ3cuVUVF_h5KAQHOoNsQU4nvyUdHs,4414
+metradar/oa_dig_func.py,sha256=k3myqBg8jM901vrpKVZWgEL3RvZ06yQLJPV7O3l0Lng,34342
+metradar/parse_pal.py,sha256=dtF0hwXNoRFVHQygRs8lUePHPVRfoCvwvdAl1gRPjwg,5072
+metradar/pgmb_io.py,sha256=4D7beTSazhKTAzGwsQHaTxSn6PMqJpF4oOmMDIMuGio,5459
+metradar/prepare_for_radar_draw.py,sha256=59PRnUZ0olef272z8yYl8cC8H7lRQTRQc2VUhaTQXY8,6512
+metradar/read_new_mosaic.py,sha256=zE5Se0One0ev0mKn-_YOw3TfNf-ATPVGA8wVZg3_zbo,1205
+metradar/read_new_mosaic_func.py,sha256=dw_-b3H5n49G47PHbEtAxboECIgt28MLyZc0leLcAIs,7259
+metradar/retrieve_cmadaas.py,sha256=apwhisjHKfN6keahMgn73cVcA9jmEFGPOEmd6iX-xXc,137463
+metradar/retrieve_micaps_server.py,sha256=cb-a_OKM1zCSw0tZ_Gk7RofE16RiRX6qq3oTUQt7lMc,85607
+metradar/rose_structer.py,sha256=3EAZ18OXg0FXR4U1DNEvwU4gxP8mSRGkGUKjAW85IPs,25293
+metradar/trans_nc_pgmb.py,sha256=EfjaswKix4f8G-T1O6UKJ8lX9F_8z_IILucX_R2xdSk,1738
+metradar/trans_new_mosaic_nc.py,sha256=dZPr625hV7msFOgUVPwFgjTrTgQOEsIywMABzansPYA,9705
+metradar/trans_polor2grid_func.py,sha256=Ad6Ii_Ew3dwbRmN7RrRpDD-wU33xm7jbZAJbIVM3EsI,6864
+metradar-0.1.0.dist-info/METADATA,sha256=Tm0STHZgVcIUibDsffJQQN8IjkKwKOpHDW9lVIO3NPo,398
+metradar-0.1.0.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+metradar-0.1.0.dist-info/top_level.txt,sha256=vwxtTxux8rWm6wgfnwf0YWYrWjN16TZQCOy4h_OVwoE,9
+metradar-0.1.0.dist-info/RECORD,,

metradar-0.1.0.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (75.8.0)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

metradar-0.1.0.dist-info/top_level.txt

@@ -0,0 +1 @@
+metradar