metradar 0.1.3__py3-none-any.whl → 0.1.5__py3-none-any.whl

metradar/draw_radar_aws_jilin_new.py

@@ -1,206 +0,0 @@
-'''
-将雷达和自动站叠加绘图
-朱文剑
-'''
-
-# %%
-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()
-
-
-
-    

metradar/nowcast_by_pysteps.py

@@ -1,219 +0,0 @@
-#!/bin/env python
-"""
-Extrapolation nowcast
-=====================
-
-This tutorial shows how to compute and plot an extrapolation nowcast using 
-Finnish radar data.
-
-"""
-
-from datetime import datetime
-import matplotlib.pyplot as plt
-import numpy as np
-from pprint import pprint
-from pysteps import io, motion, nowcasts, rcparams, verification
-from pysteps.utils import conversion, transformation
-from pysteps.visualization import plot_precip_field, quiver
-import os
-import xarray as xr
-from draw_mosaic_new import draw_mosaic
-from datetime import datetime, timedelta
-import matplotlib
-matplotlib.use('MacOSX')
-###############################################################################
-# Read the radar input images
-# ---------------------------
-#
-# First, we will import the sequence of radar composites.
-# You need the pysteps-data archive downloaded and the pystepsrc file
-# configured with the data_source paths pointing to data folders.
-
-# Selected case
-# date = datetime.strptime("201609281600", "%Y%m%d%H%M")
-date = datetime.strptime("202304211800", "%Y%m%d%H%M")
-data_source = rcparams.data_sources["fmi"]
-n_leadtimes = 12
-
-###############################################################################
-# Load the data from the archive
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-
-root_path = data_source["root_path"]
-path_fmt = data_source["path_fmt"]
-fn_pattern = data_source["fn_pattern"]
-fn_ext = data_source["fn_ext"]
-importer_name = data_source["importer"]
-importer_kwargs = data_source["importer_kwargs"]
-timestep = data_source["timestep"]
-
-# Find the input files from the archive
-fns = io.archive.find_by_date(
-    date, root_path, path_fmt, fn_pattern, fn_ext, timestep, num_prev_files=2
-)
-
-# Read the radar composites
-importer = io.get_method(importer_name, "importer")
-Z, _, metadata = io.read_timeseries(fns, importer, **importer_kwargs)
-
-# write Z to ds
-# set longitude and latitude coordinates
-slat = 33
-nlat = 37
-wlon = 109
-elon = 117
-lat = np.linspace(slat,nlat,Z.shape[1])
-lon = np.linspace(wlon,elon,Z.shape[2])
-
-# set time coordinates
-curtime = metadata['timestamps'][-1]
-curtime = np.array([curtime], dtype='datetime64[m]')
-# data = np.expand_dims(data, axis=0)
-
-# define coordinates
-time_coord = ('time', curtime)
-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': 'Composite Refelectivity', 
-            'short_name': 'cref', 'units': 'dBZ',
-            'maxv':80,
-            'minv':0}
-data = xr.Dataset({'cref':(['lat', 'lon'], np.flipud(Z[-1]), varattrs)},
-    coords={ 'lat':lat_coord, 'lon':lon_coord})
-
-outpath = 'pic'
-
-tstr = metadata['timestamps'][-1].strftime('%Y%m%d%H%M')
-outname = '%s_obs.png'%tstr
-dpi = 600
-thred=10
-# draw_mosaic(data.cref,data.lat.data,data.lon.data,slat,nlat,wlon,elon,outpath,outname,tstr,subtitle='实况',titlecolor='k',dpi=dpi,thred=thred)
-# print(outpath + os.sep + outname + ' done!')
-# Convert to rain rate
-# R, metadata = conversion.to_rainrate(Z, metadata)
-
-# # Plot the rainfall field
-# fig1 = plt.figure(figsize=(6, 6))
-# plot_precip_field(Z[-1, :, :], geodata=metadata)
-# # plt.show()
-
-# Store the last frame for plotting it later later
-# R_ = R[-1, :, :].copy()
-
-# Log-transform the data to unit of dBR, set the threshold to 0.1 mm/h,
-# set the fill value to -15 dBR
-# R, metadata = transformation.dB_transform(R, metadata, threshold=0.1, zerovalue=-15.0)
-
-# Nicely print the metadata
-# pprint(metadata)
-
-###############################################################################
-# Compute the nowcast
-# -------------------
-#
-# The extrapolation nowcast is based on the estimation of the motion field,
-# which is here performed using a local tracking approach (Lucas-Kanade).
-# The most recent radar rainfall field is then simply advected along this motion
-# field in oder to produce an extrapolation forecast.
-
-# Estimate the motion field with Lucas-Kanade
-# st = time.time()
-oflow_method = motion.get_method("LK")
-V = oflow_method(Z[-3:, :, :])
-
-# Extrapolate the last radar observation
-extrapolate = nowcasts.get_method("extrapolation")
-# R[~np.isfinite(R)] = metadata["zerovalue"]
-Z_f = extrapolate(Z[-1, :, :], V, n_leadtimes)
-for nn in range(len(Z_f)):
-    fsttime = metadata['timestamps'][-1] + timedelta(minutes=10*(nn+1))
-    tstr = fsttime.strftime('%Y%m%d%H%M')
-    # set time coordinates
-    fsttime = np.array([fsttime], dtype='datetime64[m]')
-    # data = np.expand_dims(data, axis=0)
-
-    # define coordinates
-    time_coord = ('time', fsttime)
-    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': 'Composite Refelectivity', 
-                'short_name': 'cref', 'units': 'dBZ',
-                'maxv':80,
-                'minv':0}
-    data = xr.Dataset({'cref':(['lat', 'lon'], np.flipud(Z_f[nn]), varattrs)},
-        coords={ 'lat':lat_coord, 'lon':lon_coord})
-
-    outpath = 'pic'
-
-    
-    outname = '%s_fst.png'%tstr
-    dpi = 600
-    thred=10
-    draw_mosaic(data.cref,data.lat.data,data.lon.data,slat,nlat,wlon,elon,outpath,outname,tstr,subtitle='预报',titlecolor='r',dpi=dpi,thred=thred,add_title=1,prefix_title='雷达组合反射率拼图')    
-    # print(outpath + os.sep + outname + ' done!')
-
-pass
-# Back-transform to rain rate
-# R_f = transformation.dB_transform(R_f, threshold=-10.0, inverse=True)[0]
-
-# et = time.time()
-# print("Execution time(s): ", et - st)
-# Plot the motion field
-# fig2 = plt.figure(figsize=(6, 6))
-# plot_precip_field(Z_f[-1,:,:], geodata=metadata)
-# quiver(V, geodata=metadata, step=50)
-# plt.show()
-
-###############################################################################
-# Verify with FSS
-# ---------------
-#
-# The fractions skill score (FSS) provides an intuitive assessment of the
-# dependency of skill on spatial scale and intensity, which makes it an ideal
-# skill score for high-resolution precipitation forecasts.
-
-# Find observations in the data archive
-fns = io.archive.find_by_date(
-    date,
-    root_path,
-    path_fmt,
-    fn_pattern,
-    fn_ext,
-    timestep,
-    num_prev_files=0,
-    num_next_files=n_leadtimes,
-)
-# Read the radar composites
-Z_o, _, metadata_o = io.read_timeseries(fns, importer, **importer_kwargs)
-# R_o, metadata_o = conversion.to_rainrate(R_o, metadata_o, 223.0, 1.53)
-
-# Compute fractions skill score (FSS) for all lead times, a set of scales and 1 mm/h
-fss = verification.get_method("FSS")
-scales = [2, 4, 8, 16, 32, 64, 128]
-thr = 5.0
-score = []
-for i in range(n_leadtimes):
-    score_ = []
-    for scale in scales:
-        score_.append(fss(Z_f[i, :, :], Z_o[i + 1, :, :], thr, scale))
-    score.append(score_)
-
-# plt.figure()
-# fig3 = plt.figure(figsize=(6, 6))
-# x = np.arange(1, n_leadtimes + 1) * timestep
-# plt.plot(x, score)
-# plt.legend(scales, title="Scale [km]")
-# plt.xlabel("Lead time [min]")
-# plt.ylabel("FSS ( > 5 dBZ ) ")
-# plt.title("Fractions skill score")
-# plt.show()
-
-# sphinx_gallery_thumbnail_number = 3

metradar/prepare_for_radar_draw.py

@@ -1,197 +0,0 @@
-'''
-prepare for radar_draw_aws.py
-
-'''
-
-import os
-import sys
-import pandas as pd
-
-#递归查找
-def show_files(path, all_files,allpaths):
-    # 首先遍历当前目录所有文件及文件夹
-    file_list = os.listdir(path)
-    # 准备循环判断每个元素是否是文件夹还是文件，是文件的话，把名称传入list，是文件夹的话，递归
-    for file in file_list:
-        # 利用os.path.join()方法取得路径全名，并存入cur_path变量，否则每次只能遍历一层目录
-        cur_path = os.path.join(path, file)
-
-        if file.startswith('.') or file.startswith('..'):
-            continue
-
-        # 判断是否是文件夹
-        if os.path.isdir(cur_path):
-            show_files(cur_path, all_files,allpaths)
-        else:
-            all_files.append(file)
-            allpaths.append(path)
-
-    return all_files
-
-def mainfunc(dic_param):
-
-    allfiles = []
-    allpaths = []  
-    show_files(dic_param['filepath'],allfiles,allpaths)
-
-    start_lat=[]
-    end_lat=[]
-    start_lon=[]
-    end_lon=[]
-    alloutpaths=[]
-    
-    for nn in range(len(allfiles)):
-        curout = dic_param['rootout'] + os.sep + allpaths[nn].split(os.sep)[-1]
-        # print(curout)
-        
-        # alloutpaths.append(curout)
-        # start_lat.append('')
-        # end_lat.append('')
-        # start_lon.append('')
-        # end_lon.append('')
-
-        alloutpaths.append(curout)
-        start_lat.append(dic_param['startlat'])
-        end_lat.append(dic_param['endlat'])
-        start_lon.append(dic_param['startlon'])
-        end_lon.append(dic_param['endlon'])
-
-    outdic = {
-        'sourcepath':allpaths,
-        'outpath':alloutpaths,
-        'filename':allfiles,
-        'start_lat':start_lat,
-        'end_lat':end_lat,
-        'start_lon':start_lon,
-        'end_lon':end_lon,
-    }
-
-    outpd = pd.DataFrame(outdic)
-    newpd = outpd.sort_values(by='filename')
-    newpd.to_csv(dic_param['outname'],encoding='gb18030')
-    print(dic_param['outname'] + ' done!')
-
-if __name__ == "__main__":
-    
-    # dic_param={}
-    # dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20190602-辽源'
-    # dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    # dic_param['outname'] = 'radardrawlist_20190602.csv'
-    
-    # dic_param['startlon'] = 124.2
-    # dic_param['endlon'] = 126
-    # dic_param['startlat'] = 43.5
-    # dic_param['endlat'] = 44.4
-
-    # mainfunc(dic_param=dic_param)
-
-    type = '强度'
-
-    if type == '强度':
-        outpath = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/绘图参数/回波强度'
-    elif type == '速度':
-        outpath = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/绘图参数/径向速度'
-
-    if not os.path.exists(outpath):
-        os.makedirs(outpath)
-    
-
-    #==================================================================================
-    # dic_param={}
-    # dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20120612-白城'
-    # dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    # dic_param['outname'] = outpath + os.sep + 'radardrawlist_20120612.csv'
-    
-    # if type == '强度':
-    #     dic_param['startlon'] = 122.37
-    #     dic_param['endlon'] = 122.97
-    #     dic_param['startlat'] = 45.34
-    #     dic_param['endlat'] = 45.89
-    # elif type == '速度':
-    #     dic_param['startlon'] = 122.6
-    #     dic_param['endlon'] = 122.83
-    #     dic_param['startlat'] = 45.58
-    #     dic_param['endlat'] = 45.76
-
-    # mainfunc(dic_param=dic_param)
-
-
-    # dic_param={}
-    # dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20120701-白城'
-    # dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    # dic_param['outname'] = outpath + os.sep + 'radardrawlist_20120701.csv'
-    
-    # if type == '强度':
-    #     dic_param['startlon'] = 123.1
-    #     dic_param['endlon'] = 124.32
-    #     dic_param['startlat'] = 45.11
-    #     dic_param['endlat'] = 45.97
-    # elif type == '速度':
-    #     dic_param['startlon'] = 123.65
-    #     dic_param['endlon'] = 124.32
-    #     dic_param['startlat'] = 45.45
-    #     dic_param['endlat'] = 45.97
-
-    # mainfunc(dic_param=dic_param)
-
-
-    # dic_param={}
-    # dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20150608-白城'
-    # dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    # dic_param['outname'] = outpath + os.sep + 'radardrawlist_20150608.csv'
-    
-    # if type == '强度':
-    #     dic_param['startlon'] = 122.3
-    #     dic_param['endlon'] = 124.3
-    #     dic_param['startlat'] = 44.2
-    #     dic_param['endlat'] = 45.5
-    # elif type == '速度':
-    #     dic_param['startlon'] = 122.7
-    #     dic_param['endlon'] = 123.5
-    #     dic_param['startlat'] = 44.3
-    #     dic_param['endlat'] = 44.95
-
-    # mainfunc(dic_param=dic_param)
-    
-
-    # dic_param={}
-    # dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20170905-松原'
-    # dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    # dic_param['outname'] = outpath + os.sep + 'radardrawlist_20170905.csv'
-    
-    # if type == '强度':
-    #     dic_param['startlon'] = 125.2
-    #     dic_param['endlon'] = 126.2
-    #     dic_param['startlat'] = 44.9
-    #     dic_param['endlat'] = 45.7
-    # elif type == '速度':
-    #     dic_param['startlon'] = 125.4
-    #     dic_param['endlon'] = 126.2
-    #     dic_param['startlat'] = 45
-    #     dic_param['endlat'] = 45.5
-
-    # mainfunc(dic_param=dic_param)
-
-
-    dic_param={}
-    dic_param['filepath'] = '/Users/wenjianzhu/Downloads/雷达数据-王婷婷/ar2v/20210909-长春'
-    dic_param['rootout']='/Users/wenjianzhu/Downloads/雷达数据-王婷婷/pic'
-    dic_param['outname'] = outpath + os.sep + 'radardrawlist_20210909.csv'
-    
-    if type == '强度':
-        # dic_param['startlon'] = 124.5
-        # dic_param['endlon'] = 125.8
-        # dic_param['startlat'] = 43.3
-        # dic_param['endlat'] = 44.2
-
-        dic_param['startlon'] = 124.9
-        dic_param['endlon'] = 125.6
-        dic_param['startlat'] = 43.6
-        dic_param['endlat'] = 44.1
-    elif type == '速度':
-        dic_param['startlon'] = 124.9
-        dic_param['endlon'] = 125.6
-        dic_param['startlat'] = 43.6
-        dic_param['endlat'] = 44.1
-
-    mainfunc(dic_param=dic_param)