metradar 0.1.0__py3-none-any.whl

metradar/draw_mosaic_new.py

@@ -0,0 +1,421 @@
+'''
+绘制雷达拼图
+'''
+
+
+# %%
+import os
+import matplotlib.pyplot as plt 
+import cartopy.crs as ccrs
+from parse_pal import parse
+import xarray as xr
+import numpy as np
+from matplotlib.transforms import offset_copy
+import matplotlib.patheffects as path_effects
+from matplotlib.font_manager import FontProperties
+import json
+import matplotlib as mpl
+import matplotlib
+matplotlib.use('pdf')
+# matplotlib.use('MacOSX')
+from cartopy.io.shapereader import Reader
+import pandas as pd
+from multiprocessing import cpu_count, Pool,freeze_support
+
+mappath  = 'resources/国省市县审图号(GS (2019) 3082号'
+
+g_fontfile='resources/fonts/msyhbd.ttc'
+g_font_max = FontProperties(fname=g_fontfile, size=8)
+g_font_mid = FontProperties(fname=g_fontfile, size=6)
+g_font_min = FontProperties(fname=g_fontfile, size=4)
+
+mpl.rcParams['font.size'] = 7
+mpl.rcParams['font.weight'] = 'bold'
+mpl.rcParams['font.family'] = 'Times'
+
+
+def add_gisinfo(ax,font=None):
+    city_file = 'resources/stations/2018年国家级地面气象站-2423.xlsx'
+    # city_file = '../common/resources/stations/cma_county_station_info.dat'# cma_city_station_info cma_national_station_info
+    
+    if not os.path.exists(city_file):
+        print(city_file +' not exist!')
+        return False
+
+    # cities = pd.read_csv(city_file,  delimiter=r"\s+") # cma_city_station_info
+    # cities = pd.read_csv(city_file,  delimiter=",") # cma_national_station_info
+    # cities = pd.read_table(city_file,sep='\s+',header=None,names=['stanum','lat','lon','name','city','province']) # cma_county_station_info
+    data = pd.read_excel(city_file, sheet_name = 0,header=0,skiprows=1)
+
+    names=[]
+    lat=[]
+    lon=[]
+    for i in range(data.shape[0]):
+        lat_deg = int(data['纬度'][i]/10000)
+        lat_min = int(str(data['纬度'][i])[-4:-2])
+        lat_sec = int(str(data['纬度'][i])[-2::])
+        tmp_lat = lat_deg + lat_min/60 + lat_sec/3600
+        
+        lon_deg = int(data['经度'][i]/10000)
+        lon_min = int(str(data['经度'][i])[-4:-2])
+        lon_sec = int(str(data['经度'][i])[-2::])
+        tmp_lon = lon_deg + lon_min/60 + lon_sec/3600 
+        
+        lat.append(tmp_lat)
+        lon.append(tmp_lon)
+        names.append(data['站名'][i])
+
+
+    cities = pd.DataFrame(data={'name':names, 'lat':lat, 'lon':lon})
+    # draw station information
+    font = g_font_min
+    geodetic_transform = ccrs.Geodetic()._as_mpl_transform(ax)
+    for _, row in cities.iterrows():
+        text_transform = offset_copy(geodetic_transform, units='dots', x=-5)
+        ax.plot(row['lon'], row['lat'], marker='o', color=[0.2,0.2,0.2], linestyle='None',
+                markersize=2, alpha=0.8, transform=ccrs.PlateCarree())
+        ax.text(row['lon']-0.01, row['lat'], row['name'], clip_on=True, # city_name
+                verticalalignment='center', horizontalalignment='right',
+                transform=text_transform, fontproperties=font, color=[0.2,0.2,0.2],
+                path_effects=[path_effects.Stroke(linewidth=1, foreground=[1,1,1]),
+                path_effects.Normal()])
+                
+
+def draw_gisinfo(ax,font=None):
+
+    city_file = 'resources/stations/cma_city_station_info.dat'# cma_city_station_info cma_national_station_info
+    
+    if not os.path.exists(city_file):
+        print(city_file +' not exist!')
+        return False
+
+    cities = pd.read_csv(city_file,  delimiter=r"\s+") # cma_city_station_info
+    # cities = pd.read_csv(city_file,  delimiter=",") # cma_national_station_info
+    # cities = pd.read_table(city_file,sep='\s+',header=None,names=['stanum','lat','lon','name','city','province']) # cma_county_station_info
+    # extract subregion
+
+    # targetprovince = '吉林省'
+    # flag1 = cities['province']==targetprovince
+    # cityname = np.unique(cities[flag1]['city'])
+    # cityname = np.unique(cities['city'])
+    cityname = np.unique(cities['city_name'])
+    names=[]
+    lat=[]
+    lon=[]
+
+    curdata = cities
+    for name in cityname:
+        flag2 =  curdata['city_name'] == name
+        tmp_lat = np.mean(curdata[flag2]['lat'])
+        tmp_lon = np.mean(curdata[flag2]['lon'])
+        # if name.find('藏族') >0: name = name.strip(name[name.find('藏族')::])
+        # if name.find('蒙古') >0: name = name.strip(name[name.find('蒙古')::])
+        names.append(name)
+        lon.append(tmp_lon)
+        lat.append(tmp_lat)
+
+    cities = pd.DataFrame(data={'name':names, 'lat':lat, 'lon':lon})
+    # draw station information
+    font = g_font_min
+    geodetic_transform = ccrs.Geodetic()._as_mpl_transform(ax)
+    for _, row in cities.iterrows():
+        text_transform = offset_copy(geodetic_transform, units='dots', x=-5)
+        ax.plot(row['lon'], row['lat'], marker='o', color=[0.2,0.2,0.2], linestyle='None',
+                markersize=2, alpha=0.8, transform=ccrs.PlateCarree())
+        ax.text(row['lon']-0.05, row['lat'], row['name'], clip_on=True, # city_name
+                verticalalignment='center', horizontalalignment='right',
+                transform=text_transform, fontproperties=font, color=[0.2,0.2,0.2],
+                path_effects=[path_effects.Stroke(linewidth=1, foreground=[1,1,1]),
+                path_effects.Normal()])
+                
+    
+def add_china_map_2cartopy(ax, facecolor='none',transform=None,name='province',
+                           edgecolor='k', lw=0.8, **kwargs):
+    """
+    Draw china boundary on cartopy map.
+
+    :param ax: matplotlib axes instance.
+    :param name: map name.
+    :param facecolor: fill color, default is none.
+    :param edgecolor: edge color.
+    :param lw: line width.
+    :return: None
+    """
+
+    # map name
+    names = {'nation': "BOUL_G", 'province': "BOUL_S",
+            'city': "BOUL_D", 'county': "BOUL_X",}
+
+    if not name in {'nation','province','city','county'}:
+        print(name + ' should be nation or province or city or county !')
+        return False
+
+
+    # add 省界
+    if transform is None:
+        transform = ccrs.PlateCarree()
+
+    if name == 'county':
+        # 添加县边界
+        shpfile = mappath + os.sep + names[name] + ".shp"
+        ax.add_geometries(
+            Reader(shpfile).geometries(), transform,
+            path_effects=[path_effects.Stroke(linewidth=lw, foreground=[1,1,1]),path_effects.Normal()],
+            facecolor=facecolor, edgecolor=edgecolor, lw=lw-0.5, **kwargs)
+
+    if name == 'city':
+        # 添加市边界
+        shpfile = mappath + os.sep + names[name] + ".shp"
+        ax.add_geometries(
+            Reader(shpfile).geometries(), transform,
+            path_effects=[path_effects.Stroke(linewidth=lw+0.2, foreground=[1,1,1]),path_effects.Normal()],
+            facecolor=facecolor, edgecolor=edgecolor, lw=lw-0.2, **kwargs)
+        
+    # 添加省边界
+    if name == 'province':
+        shpfile = mappath + os.sep + names[name] + ".shp"
+        ax.add_geometries(
+            Reader(shpfile).geometries(), transform,
+            path_effects=[path_effects.Stroke(linewidth=lw+0.5, foreground=[1,1,1]),path_effects.Normal()],
+            facecolor=facecolor, edgecolor=edgecolor, lw=lw+0.2, **kwargs)
+
+def add_private_shp(ax,facecolor='none',transform=None,name='province',edgecolor='k', lw=0.8, **kwargs):
+    mappath  = '/Users/wenjianzhu/Library/CloudStorage/OneDrive-个人/推广应用/GR2/吉林/地理信息/归档'
+
+    shpfile = mappath + os.sep + "jilin.shp"
+
+    ax.add_geometries(
+            Reader(shpfile).geometries(), transform,
+            path_effects=[path_effects.Stroke(linewidth=lw+0.2, foreground=[1,1,1]),path_effects.Normal()],
+            facecolor=facecolor, edgecolor=edgecolor, lw=lw-0.2, **kwargs)
+
+def cm_precip():
+    """
+    Standardized colormaps from National Weather Service
+    https://github.com/blaylockbk/pyBKB_v2/blob/master/BB_cmap/NWS_standard_cmap.py
+
+    Range of values:
+        metric: 0 to 762 millimeters
+        english: 0 to 30 inches
+    """
+    # The amount of precipitation in inches
+    # a = [0,.01,.1,.25,.5,1,1.5,2,3,4,6,8,10,15,20,30]
+    a = [0,0.01,0.5,2.5,5,10,15,20,30,40,60,80,100,150,200,250]
+
+    # Normalize the bin between 0 and 1 (uneven bins are important here)
+    norm = [(float(i)-min(a))/(max(a)-min(a)) for i in a]
+
+    # Color tuple for every bin
+    C = np.array(# [255,255,255],
+        [ [199,233,192], [161,217,155], [116,196,118], [49,163,83], 
+        [0,109,44], [255,250,138], [255,204,79], [254,141,60], [252,78,42], 
+        [214,26,28], [173,0,38], [112,0,38], [59,0,48], [76,0,115], [255,219,255]])
+
+    # Create a tuple for every color indicating the normalized 
+    # position on the colormap and the assigned color.
+    # COLORS = []
+    # for i, n in enumerate(norm):
+    #     COLORS.append((n, np.array(C[i])/255.))
+
+    # Create the colormap
+    # cmap = mpl.colors.LinearSegmentedColormap.from_list("precipitation", COLORS)
+    cmap, norm = mpl.colors.from_levels_and_colors(np.array(a), np.array(C)/255.0, extend='neither')
+    return cmap, norm
+
+
+def draw_mosaic(data,lat,lon,slat=None,nlat=None,wlon=None,elon=None,outpath=None,outname=None,tstr=None,prefix_title='雷达组合反射率拼图 ',units='dBZ', 
+                subtitle='',titlecolor='r',thred=None,cb_ratio=0.618,dpi=400,add_title=False):
+    
+    if units == 'dBZ':
+        var_cmap,norm = parse('resources/gr2_colors/BR_WDTB_Bright.pal')#BR_WDTB_Bright.pal  default_BR_PUP2.pal BR_AVL_BroadcastNegatives.pal
+    elif units == 'mm':
+        var_cmap,norm = cm_precip()
+    elif units == '%':
+        var_cmap,norm = parse('resources/gr2_colors/default_CC_ROSE.pal')
+    elif units == 'dB':
+        var_cmap,norm = parse('resources/gr2_colors/default_ZDR_ROSE.pal')
+    elif units == 'deg/km':
+        var_cmap,norm = parse('resources/gr2_colors/default_KDP_ROSE.pal')
+    elif units == 'cat':
+        var_cmap,norm = parse('resources/gr2_colors/default_HCAS.pal')
+    else:
+        pass
+        
+    # newdata.data.shape[1],newdata.data.shape[0]
+    fig= plt.figure(figsize=[6,4],dpi=dpi,clear=True)
+    map_panel_axes = [0.1, 0.1, .8, .8]
+    # plt.subplots_adjust(top=1,bottom=0,left=0,right=1,hspace=0,wspace=0)
+    # projs=ccrs.LambertConformal(central_longitude=(wlon + elon)/2,central_latitude=(slat + nlat)/2,standard_parallels=(20,50))
+    projs = ccrs.PlateCarree()
+    transforms = ccrs.PlateCarree()
+    
+    ax = plt.axes(map_panel_axes, projection=projs)
+
+    extent = [wlon, elon, slat, nlat]
+    ax.set_extent(extent, crs=ccrs.PlateCarree()) 
+
+    # add terrain
+    # try:
+    #     request = TDT_Hillshade()  # 地形阴影
+    #     ax.add_image(request, 10)  # level=10 缩放等
+    # except:
+    #     pass
+    if not thred is None:
+        newdata = np.ma.masked_array(data, data < thred)
+    else:
+        newdata = data
+    grid_lat,grid_lon=np.meshgrid(lat,lon)
+    pcm = ax.pcolormesh(grid_lon,grid_lat,newdata.T,cmap=var_cmap,norm=norm,transform=transforms)
+
+    # add_private_shp(ax)
+
+    # add_china_map_2cartopy(ax,name='county')
+    # add_china_map_2cartopy(ax,name='city')
+    add_china_map_2cartopy(ax,name='province',lw=0.5)
+
+    draw_gisinfo(ax)
+
+
+    # 添加地图网格
+    # ax.minorticks_on()
+    # ax.tick_params(direction='out', length=6, width=2, colors='r',grid_color='r', grid_alpha=0.5)
+    if elon - wlon > 20:
+        lon_step = 5
+    elif elon - wlon > 5:
+        lon_step = 1
+    elif elon - wlon > 1:
+        lon_step = 0.5
+    else:
+        lon_step = 0.2
+
+    if nlat - slat > 20:
+        lat_step = 5
+    elif nlat - slat > 5:
+        lat_step = 1
+    elif nlat - slat > 1:
+        lat_step = 0.5
+    else:
+        lat_step = 0.2
+
+    lb=ax.gridlines(draw_labels=True,x_inline=False, y_inline=False,xlocs=np.arange(int(wlon)-1,int(elon)+1,lon_step), ylocs=np.arange(int(slat)-1,int(nlat)+1,lat_step),
+                    linewidth=0.1, color='k', alpha=0.8, linestyle='--')
+    # point_lat = 30.72
+    # point_lon = 108.57
+    # ax.plot(point_lon,point_lat,marker='o',color='r',markersize=4,transform=ccrs.PlateCarree())
+    lb.top_labels = None
+    lb.right_labels = None
+    lb.rotate_labels = False
+    # 'Composite Reflectivity '
+    titlestr1 = prefix_title + tstr[0:8] + 'T' + tstr[8:12] + ' UTC ' + subtitle
+
+    if add_title:
+        ax.set_title(titlestr1,fontproperties=g_font_mid,loc='left',verticalalignment='top',color=titlecolor)
+    # add colorbar
+    ax_pos = ax.get_position().bounds
+    cbar_axes = [ax_pos[0]+ ax_pos[2]+0.005,(ax_pos[1]+ax_pos[1]+ax_pos[3])/2-ax_pos[3]*cb_ratio/2,0.02,ax_pos[3]*cb_ratio]
+    cb_ax = fig.add_axes(cbar_axes)
+    if units == 'cat':
+        plt.colorbar(mappable=pcm,cax=cb_ax,extend='both',ticks=np.arange(0.5,11.5,1)) # extend='both',ticks=np.arange(-5,80,5)
+    elif units == 'dBZ':
+        plt.colorbar(mappable=pcm,cax=cb_ax,extend='both',ticks=np.arange(0,80,5))
+    else:
+        plt.colorbar(mappable=pcm,cax=cb_ax,extend='both')
+    cb_ax.text(1.8, 1.035, units, transform=cb_ax.transAxes,va="center", ha="center", font=g_font_min ) # 0.5,1.02
+
+    if units == 'cat':
+        cb_ax.set_yticklabels(['空','地物','晴空','干雪','湿雪','冰晶','霰','大雨滴','小到中雨','大雨','冰雹',],font=g_font_min)
+
+    
+    # DPI=1
+    # fig,ax = plt.subplots(figsize=[newdata.data.shape[1],newdata.data.shape[0]],dpi=DPI,clear=True)
+    # plt.subplots_adjust(top=1,bottom=0,left=0,right=1,hspace=0,wspace=0)
+
+    # data = np.ma.masked_array(newdata.data, newdata.data < 10)
+    # plt.pcolormesh(data,cmap=ref_cmap,norm=norm)
+    # plt.axis("off")
+    if os.path.exists(outpath):
+        if outname.endswith('.eps'):
+            plt.savefig(outpath + os.sep  + outname, transparent=False, dpi=dpi,papertype='a4')
+        else:
+            fig.savefig(outpath + os.sep  + outname, bbox_inches="tight",transparent=False, dpi=dpi,)
+        # plt.close()
+        print(outpath + os.sep + outname + ' saved!')
+        # plt.tight_layout()
+        # plt.show()
+
+def draw_single(param):
+    filename = param['filename']
+    filepath = param['filepath']
+    outpath = param['outpath']
+    slat = param['slat']
+    nlat = param['nlat']
+    wlon = param['wlon']
+    elon = param['elon']
+    dpi = param['dpi']
+    thred = param['thred']
+    
+
+    outname = filename.replace('.nc','.png')
+    tstr = filename.split('.')[2] + filename.split('.')[3]
+    data = xr.open_dataset(filepath + os.sep + filename)
+    try:
+        newdata = data.cref.sel(latitude=slice(slat,nlat),longitude=slice(wlon,elon))
+        draw_mosaic(newdata,newdata.latitude.data,newdata.longitude.data,slat,nlat,wlon,elon,outpath,outname,tstr,dpi=dpi,thred=thred)
+    except:
+        newdata = data.cref.sel(lat=slice(slat,nlat),lon=slice(wlon,elon))
+        draw_mosaic(newdata,newdata.lat.data,newdata.lon.data,slat,nlat,wlon,elon,outpath,outname,tstr,dpi=dpi,thred=thred)
+    
+    
+    
+
+if __name__ == "__main__":      
+    freeze_support()
+    filepath = '/Users/wenjianzhu/Downloads/CMADAAS/rdmosaic_bin/202307'
+ 
+    outpath = '/Users/wenjianzhu/Downloads/20230704万州滑坡/pic'
+    
+    if not os.path.exists(outpath):
+        os.makedirs(outpath)
+
+    # curparam=dict()
+    # curparam['filepath'] = filepath
+    # curparam['filename'] = filename
+    # curparam['outpath'] = outpath
+    # curparam['slat'] = 33
+    # curparam['nlat'] = 37
+    # curparam['wlon'] = 109
+    # curparam['elon'] =117
+    # curparam['dpi'] = 600
+    # curparam['thred'] = 10
+
+    # draw_single(curparam)
+
+    
+    params=[]
+    files = os.listdir(filepath)
+    files = sorted(files)
+    for filename in files:
+        if filename.startswith('.') or filename.startswith('..'):
+            continue
+        # if filename.split('.')[3][2:4] != '00' and filename.split('.')[3][2:4] != '12' and filename.split('.')[3][2:4] != '30' and filename.split('.')[3][2:4] != '48':
+        #     continue
+        curparam=dict()
+        curparam['filepath'] = filepath
+        curparam['filename'] = filename
+        curparam['outpath'] = outpath
+        curparam['slat'] = 28
+        curparam['nlat'] = 32.5
+        curparam['wlon'] = 105.5
+        curparam['elon'] =110.5
+        curparam['dpi'] = 300
+        curparam['thred'] = 15
+
+        params.append(curparam)
+
+    MAXP = int(cpu_count()*0.8)
+    pools = Pool(MAXP)
+
+    pools.map(draw_single, params)
+    pools.close()
+    pools.join()
+

metradar/draw_radar_aws_jilin_new.py

@@ -0,0 +1,206 @@
+'''
+将雷达和自动站叠加绘图
+朱文剑
+'''
+
+# %%
+import pandas as pd
+import warnings
+warnings.filterwarnings('ignore')
+
+from draw_radar_comp_func import DRAW_RADAR_OTHER,ini_params
+from multiprocessing import cpu_count, Pool,freeze_support
+
+
+def draw_all(params):
+    # sourcepath,outpath,filename,start_lat,end_lat,start_lon,end_lon
+    # 河南北部区域
+    sourcepath = params['sourcepath']
+    outpath = params['outpath']
+    radfilename = params['radfilename']
+    slat = params['slat']
+    nlat = params['nlat']
+    wlon = params['wlon']
+    elon = params['elon']
+
+    timestr = radfilename[5:13] + '.' +  radfilename[14:20]
+
+    fontname='resources/fonts/msyhbd.ttc'
+
+    # 添加中文地名
+
+    filename = 'resources/中文地理信息原始文件/吉林省所有站点信息.xls'
+  
+    data_guojia = pd.read_excel(filename, sheet_name = '国家站55',header=0,skiprows=0,index_col=None,)
+    for jj in range(data_guojia.shape[0]):
+        curidx = data_guojia['站名'][jj].find('国家')
+        data_guojia['站名'][jj] = data_guojia['站名'][jj].replace( data_guojia['站名'][jj][curidx:],'')
+        data_guojia['站名'][jj] = data_guojia['站名'][jj].replace( data_guojia['组织机构'][jj],'')
+
+
+    data_guojiatianqi = pd.read_excel(filename, sheet_name = '国家天气站333',header=0,skiprows=0,index_col=None,)
+                        
+    for jj in range(data_guojiatianqi.shape[0]):
+        data_guojiatianqi['站名'][jj] = data_guojiatianqi['站名'][jj].replace('国家气象观测站','')
+        data_guojiatianqi['站名'][jj] = data_guojiatianqi['站名'][jj].replace( data_guojiatianqi['组织机构'][jj],'')
+        if data_guojiatianqi['站名'][jj].find('尔罗斯')>=0:
+            data_guojiatianqi['站名'][jj] = data_guojiatianqi['站名'][jj].replace('尔罗斯','')
+
+    data_quyuzhan = pd.read_excel(filename, sheet_name = '区域站1048',header=0,skiprows=0,index_col=None,)
+    for jj in range(data_quyuzhan.shape[0]):
+        data_quyuzhan['站名'][jj] = data_quyuzhan['站名'][jj].replace('气象观测站','')
+        
+        data_quyuzhan['站名'][jj] = data_quyuzhan['站名'][jj].replace(data_quyuzhan['组织机构'][jj],'')
+        if data_quyuzhan['站名'][jj].find('尔罗斯')>=0:
+            data_quyuzhan['站名'][jj] = data_quyuzhan['站名'][jj].replace('尔罗斯','')
+        
+
+
+    lat=[]
+    lon=[]
+    staname=[]
+    for ng in range(len(data_guojia['经度'])):
+        lon.append(data_guojia['经度'][ng])
+        lat.append(data_guojia['纬度'][ng])
+        staname.append(data_guojia['站名'][ng])
+    
+    for ng in range(len(data_guojiatianqi['经度'])):
+        lon.append(data_guojiatianqi['经度'][ng])
+        lat.append(data_guojiatianqi['纬度'][ng])
+        staname.append(data_guojiatianqi['站名'][ng])
+    
+    # for ng in range(len(data_quyuzhan['经度'])):
+    #     lon.append(data_quyuzhan['经度'][ng])
+    #     lat.append(data_quyuzhan['纬度'][ng])
+    #     staname.append(data_quyuzhan['站名'][ng])
+
+
+
+    # radfilepath='/Users/wenjianzhu/Downloads/ZZHN'
+    radfilepath = sourcepath
+    params['radarfile_path'] = radfilepath
+    params['radarfile_name'] = radfilename
+    params['pic_path'] = outpath
+    params['timestr'] = timestr
+    params['slat'] = slat
+    params['nlat'] = nlat
+    params['wlon'] = wlon
+    params['elon'] = elon
+    params['fontfile'] = fontname
+    params['gis_lats'] = lat
+    params['gis_lons'] = lon
+    params['gis_name'] = staname
+    params['breplace'] = True
+    params['bdraw_crs'] = True
+    params['ref_colorfile'] = '../common/gr2_colors/default_BR_PUP2.pal'
+    params['vel_colorfile'] = '../common/gr2_colors/default_BV_PUP2.pal'
+    params['figsize_width'] = 4
+    params['fontsize_gis'] = 5
+    params['fontsize_colorbar'] = 5
+    params['fontsize_title'] = 6
+    params['mapcolor'] = [0/255,0/255,0/255]
+    params['dpi'] = 800
+    params['pic_format'] = 'jpg'
+    params['bdraw_title_ppi'] = False
+
+    # params={'radarfile_path':radfilepath,
+    #         'radarfile_name':radfilename,
+    #         'mosaicfile_path':'',
+    #         'mosaicfile_name':'',
+    #         'pic_path':outpath,
+    #         'timestr':timestr,
+    #         'aws_min_file_path':'',
+    #         'aws_min_file_name':'',
+    #         'aws_hour_file_path':'',
+    #         'aws_hour_file_name':'',
+    #         'gis_name':staname,
+    #         'gis_lats':lat,
+    #         'gis_lons':lon,
+    #         'slat':slat,
+    #         'nlat':nlat,
+    #         'wlon':wlon,
+    #         'elon':elon,
+    #         'ref_colorfile':'../common/gr2_colors/default_BR_PUP2.pal',
+    #         'vel_colorfile':'../common/gr2_colors/default_BV_PUP2.pal',
+    #         'fontfile':fontname,
+    #         'dpi':800,
+    #         'pic_format':'png',
+    #         'figsize_width':4,
+    #         'fontsize_gis':5,
+    #         'fontsize_colorbar':5,
+    #         'fontsize_title':6,
+    #         'mapcolor':[0/255,0/255,0/255],
+    #         'breplace':True, #如果图片文件已存在，是否重新绘制
+    #         'bdraw_crs':False
+    # }
+
+    _draw_radar_other = DRAW_RADAR_OTHER(params)
+
+    _draw_radar_other.read_vol_data()
+
+    # _draw_radar_other.draw_ref_alone(subdir='回波强度',tilt=0,thred=-5)
+    _draw_radar_other.draw_ref_alone(subdir='回波强度',tilt=1,thred=-5) 
+
+
+    # _draw_radar_other.draw_vel_alone(subdir='径向速度',tilt=0)
+    # _draw_radar_other.draw_vel_alone(subdir='径向速度',tilt=1)
+    # _draw_radar_other.draw_vel_alone(subdir='径向速度',tilt=2)
+
+    # _draw_radar_other.draw_vel_pre()
+    # _draw_radar_other.draw_vel_wind_barb()
+    # _draw_radar_other.draw_vel_wind_quiver()
+    # _draw_radar_other.draw_ref_pre()
+    # _draw_radar_other.draw_ref_pre_wind_barb()
+
+    # _draw_radar_other.get_cref_from_radar([_draw_radar_other.g_rad_lat,_draw_radar_other.g_rad_lon])
+    # _draw_radar_other.get_cref_from_mosaicfile()
+    # _draw_radar_other.get_cref_from_radar([35.6,114.0])
+    # _draw_radar_other.draw_cref_pre()
+    # _draw_radar_other.draw_cref_wind_barb()
+    # _draw_radar_other.draw_cref_wind_quiver()
+    # _draw_radar_other.draw_cref_pre_wind_barb()
+
+
+# %%
+import os
+if __name__ == '__main__':
+    pass
+    freeze_support()
+    # Pool不支持跨CPU的虚拟服务器，会出现页面不足的错误提示
+
+    paramfilepath = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/绘图参数/回波强度'
+    # paramfilepath = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/绘图参数/径向速度'
+    # drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20120612.csv',encoding='gb18030')
+    # drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20120701.csv',encoding='gb18030')
+    # drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20150608.csv',encoding='gb18030')
+    # drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20170905.csv',encoding='gb18030')
+    drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20210909.csv',encoding='gb18030')
+    # drawinfo = pd.read_csv(paramfilepath + os.sep + 'radardrawlist_20190602.csv',encoding='gb18030')
+
+    params = []
+    nums=drawinfo.shape[0]
+    # nums=1
+    for nn in range(nums):
+        pass
+        curparam=ini_params()
+        curparam['sourcepath'] = drawinfo['sourcepath'].iloc[nn]
+        curparam['outpath'] = drawinfo['outpath'].iloc[nn]
+        curparam['radfilename'] = drawinfo['filename'].iloc[nn]
+        curparam['slat'] = drawinfo['start_lat'].iloc[nn]
+        curparam['nlat'] = drawinfo['end_lat'].iloc[nn]
+        curparam['wlon'] = drawinfo['start_lon'].iloc[nn]
+        curparam['elon'] = drawinfo['end_lon'].iloc[nn]
+        # aws_min_delta_t_file_path
+        params.append(curparam)
+
+    # MAXP = int(cpu_count()*0.5)
+    MAXP=1#nums
+    pools = Pool(MAXP)
+
+    pools.map(draw_all, params)
+    pools.close()
+    pools.join()
+
+
+
+