import nmc_verification.nmc_vf_base.io as io import nmc_verification.nmc_vf_base.function.gxy_sxy as gxy_sxy import os def get_tests_path(path=__file__): path = os.path.dirname(path) if path[-5:] == 'tests': return path elif len(path) == 3: return else: path = get_tests_path(path) return path if __name__ == '__main__': path = get_tests_path() grib_data = io.read_griddata.read_from_nc(path + '/data/BT18010120.012.nc') sta_data = io.read_stadata.read_station( path + '/data/station_table.txt', columns=['id', 'lat', 'lon', 'alt']) in_sta = gxy_sxy.interpolation_linear(grib_data, sta_data) plot.frequency_histogram_muti_model(in_sta, [in_sta, in_sta, in_sta], clevs=[0, 0.1, 0.3, 0.6, 0.9], save_path='aa.png')
if path[-5:] == 'tests': return path elif len(path) == 3: return else: path = get_tests_path(path) return path if __name__ == '__main__': # path = get_tests_path() # # grib_data = io.read_griddata.read_from_nc(path + '/data/BT18010120.012.nc') # # sta_data = io.read_stadata.read_station(path + '/data/station_table.txt', columns=['id', 'lat', 'lon', 'alt']) # # in_sta = gxy_sxy.interpolation_linear(grib_data, sta_data) # # plot.frequency_histogram_muti_model(in_sta, [in_sta, in_sta, in_sta], clevs=[0, 0.1, 0.3, 0.6, 0.9], # save_path='aa.png') path = get_tests_path() grib_data = io.read_griddata.read_from_nc(path + '\data\BT19060108.024.nc') ob_data = io.read_stadata.read_from_micaps3(path + '\data\BT19060208.000') # print(sta_data) fo_data = gxy_sxy.interpolation_linear(grib_data, ob_data) ob_data['time'] = fo_data['time'] plot.frequency_histogram_muti_model(ob_data, [fo_data], clevs=[0, 0.1, 0.3, 0.6, 0.9], save_path='aa.png')