shancx 1.8.92__py3-none-any.whl → 1.9.33.218__py3-none-any.whl

shancx/makenetCDFN.py

@@ -1,42 +0,0 @@
-import xarray as xr
-import numpy as np
-def convert_longitude(nc_input_path, nc_output_path):
-    ds = xr.open_dataset(nc_input_path)
-    lon = 'lon'
-    if lon not in ds.coords:
-        raise ValueError("输入文件中没有'lon'坐标")
-    ds[lon] = xr.where(ds[lon] > 180, ds[lon] - 360, ds[lon])
-    ds = ds.sortby(lon)
-    ds = ds.sortby('lat', ascending=False)
-    ds.to_netcdf(nc_output_path)
-    print(f"成功将数据转换并保存到 {nc_output_path}")
-convert_longitude('CMORPH2_0.25deg-30min_202410160100.RT.nc', 'CMORPH2_0.25deg-30min_202410160100.RT_N.nc')
-
-import netCDF4 as nc
-import numpy as np
-with nc.Dataset(path) as dataNC:
-    hourlyPrecipRateGC = dataNC["hourlyPrecipRateGC"][:] 
-    latArr = dataNC["Latitude"][:]  
-    lonArr = dataNC["Longitude"][:] 
-    latArr_flipped = latArr[::-1]
-    hourlyPrecipRateGC_flipped = hourlyPrecipRateGC[::-1, :]
-
-with nc.Dataset(path) as dataNC:
-    ref = dataNC["var"][:][::-1]  # 读取数据并翻转第一个维度（通常是纬度）
-    latArr = dataNC["lat"][:][::-1]  # 翻转纬度数组
-    lonArr = dataNC["lon"][:]  # 
-
-
-import netCDF4 as nc
-with nc.Dataset(path) as dataNC:
-    # 翻转数据的纬度维度
-    ref = dataNC["var"][:][::-1]  # 假设第一个维度是纬度    
-    # 翻转纬度数组
-    latArr = dataNC["lat"][:][::-1]    
-    # 经度转换：从 [0, 359] 转换为 [-179, 179]
-    lonArr = dataNC["lon"][:]
-    lonArr = ((lonArr + 180) % 360) - 180  # 将 [0, 359] 转换为 [-180, 180)    
-    # 排序经度数组，并相应调整数据
-    lon_order = lonArr.argsort()  # 获取排序索引
-    lonArr = lonArr[lon_order]    # 按照索引重新排序经度数组
-    ref = ref[:, lon_order]       # 重新排序数据的经度维度

shancx/mkIMGSCX.py

@@ -1,92 +0,0 @@
-
-import numpy as np
-from matplotlib.colors import ListedColormap
-import matplotlib.pyplot as plt
-import datetime
-cmp_hjnwtx={}
-
-newcolorsNMC = np.array([
-    [68,157,237, 255],
-    [98,230,234, 255],
-    [104,249,82, 255],
-    [0,215,46, 255],
-    [0,143,27, 255],
-    [254,254,63, 255],
-    [231,192,48, 255],
-    [255,154,41, 255],
-    [255,19,27, 255],
-    [215,14,21, 255],
-    [193,11,18, 255],
-    [255,28,236, 255],
-    [152,15,177, 255],
-    [175,145,237, 255]])/255
-cmp_hjnwtx["radar_nmc"] = ListedColormap(newcolorsNMC)
-
-newcolorsNMC = np.array([
-    [68,157,237, 255],
-    [98,230,234, 255],
-    [104,249,82, 255],
-    [0,215,46, 255],
-    [0,143,27, 255],
-    [254,254,63, 255],
-    [231,192,48, 255],
-    [255,154,41, 255],
-    [255,19,27, 255],
-    [215,14,21, 255],
-    [193,11,18, 255],
-    [255,28,236, 255],
-    [152,15,177, 255],
-    [175,145,237, 255]])/255
-cmp_hjnwtx["radar_moc"] = ListedColormap(newcolorsNMC)
-
-
-newcolorsPRE = np.array([
-    [128, 255, 255, 255],
-    [35, 182, 254, 255],
-    [0, 120, 180, 255],
-    [0, 82, 202, 255],
-    [0, 16, 220, 255],
-    [150, 2, 244, 255],
-    [110, 0, 182, 255],
-    [77, 0, 130, 255]])/255
-cmp_hjnwtx["pre_tqw"] = ListedColormap(newcolorsPRE)
-
-newcolorsWS = np.array([
-    [75, 140, 244, 255],
-    [0, 89, 235, 255],
-    [36, 173, 0, 255],
-    [18, 129, 1, 255],
-    [3, 64, 4, 255],
-    [218, 183, 5, 255],
-    [179, 125, 1, 255],
-    [155, 70, 16, 255],
-    [253, 3, 127, 255],
-    [255, 0, 55, 255],
-    [233, 0, 3, 255]])/255
-cmp_hjnwtx["ws_nmic"] = ListedColormap(newcolorsWS)
-
-import os
-def mkDir(path):
-    if "." in path:
-        os.makedirs(os.path.dirname(path),exist_ok=True)
-    else:
-        os.makedirs(path, exist_ok=True)
-
-def Radar_Nmc(array_dt,temp = "850"):
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M")
-    if len(array_dt.shape)==3:
-        for i , img_ch_nel in enumerate(array_dt): 
-            plt.imshow(img_ch_nel,vmin=0,vmax=100,cmap=cmp_hjnwtx["radar_nmc"])
-            plt.colorbar()
-            outpath = f"./radar_nmc/{temp}_{now_str}.png"
-            mkDir(outpath)
-            plt.savefig(outpath)
-            plt.close()
-    if len(array_dt.shape)==2:
-            plt.imshow(array_dt,vmin=0,vmax=100,cmap=cmp_hjnwtx["radar_nmc"])
-            plt.colorbar()
-            outpath = f"./radar_nmc/{temp}_{now_str}.png"
-            mkDir(outpath)
-            plt.savefig(outpath)
-            plt.close() 
-            

shancx/netCDF.py

@@ -1,130 +0,0 @@
-
-
-import netCDF4 as nc
-import numpy as np
-with nc.Dataset(path) as dataNC:
-    # 获取 'time' 变量
-    name_data = dataNC.variables[list(dataNC.variables)[0]][:] 
-    print(list(dataNC.variables))   
-    # 打印时间数据
-    name = list(dataNC.variables)[0]
-    print(f"{name}数据:")
-    print(name_data)
-
-with nc.Dataset(path) as dataNC:
-    ref = dataNC["var"][:][::-1]*100
-    latArr = dataNC["lat"][:][::-1]
-    lonArr = dataNC["lon"][:]
-
-with nc.Dataset(path) as dataNC:
-    ref = dataNC["var"][:] 
-    latArr = dataNC["lat"][:] 
-    lonArr = dataNC["lon"][:]
-
-def readnetCDF4(path):
-    with nc.Dataset(path) as dataNC:
-        # 获取 'time' 变量
-        name_data = dataNC.variables[list(dataNC.variables)[0]][:] 
-        print(list(dataNC.variables))   
-        # 打印时间数据
-        name = list(dataNC.variables)[0]
-        print(f"{name}数据:")
-        print(name_data)
-def readnetCDF4all(path):
-    with nc.Dataset(path) as dataNC:
-        # 获取 'time' 变量
-        print("开始读取呢台CDF4数据")
-        print(list(dataNC.variables))  
-        for i in list(dataNC.variables):
-            name_data = dataNC.variables[i][:]              
-            # 打印时间数据
-            name = i
-            print(f"{name}数据:")
-            print(name_data)
-            
-import pygrib
-def grib(path):    
-    grbs = pygrib.open(path) #"/mnt/wtx_weather_forecast/gfs_110/20231030/18/gfs.t18z.pgrb2.0p25.f002"
-    for grb in grbs:
-        print(grb)
- 
-import json
-import gzip
-import pandas as pd
-# Assuming the file is a gzipped JSON
-file_path = '/mnt/wtx_weather_forecast/scx/GeoEnvData/rawData/ZW_1.gz'
-with gzip.open(file_path, 'rt') as file:  # 'rt' mode to read as text
-    data = json.load(file)
-
-# Normalize JSON data into a DataFrame
-df = pd.json_normalize(data)
-
-"""
-def read_grib(fn, vars, var_name,high,stepRange=30):
-    # logger.info(f"{fn}, {vars}, {var_name},{high},{stepRange}")
-    fid = pygrib.open(fn)
- 
-    data = {}
-    for ivar, ivar_name in zip(vars, var_name): 
-        grbs = fid.select(shortName = ivar,level=high,stepRange=str(stepRange)) 
-        d = []
-        for grb in grbs:
-            d.append(grb.values)
-        d = np.squeeze(np.array(d))
-        data.update({ivar_name : d})
-    lat, lon = fid[1].latlons()
-    latArr = lat[:, 0]
-    lonArr = lon[0]
-    d_clip = clip(data[ivar_name], env, latArr[0], lonArr[0], 0.25)    
-    # d = zoom(d_clip, [4201/169,6201/249], order=1)[:-1, :-1]
-    data.update({'lon':lon})
-    data.update({'lat': lat})
-    return d_clip
-"""
-
-"""
-
-import netCDF4 as nc
-
-# 设定NetCDF文件的路径
-path = 'your_file.nc'  # 请替换为您的实际文件路径
-
-# 打开NetCDF文件
-with nc.Dataset(path) as dataNC:
-    # 获取所有变量的名称
-    variables = dataNC.variables.keys()    
-    # 输出所有变量
-    print("变量列表:")
-    for var in variables:
-        print(var)
-
-with nc.Dataset(path) as dataNC:
-    for var_name in dataNC.variables:
-        var = dataNC.variables[var_name]
-        print(f"变量名称: {var_name}")
-        print(f"数据类型: {var.dtype}")
-        print(f"维度: {var.dimensions}")
-        print(f"属性: {var.__dict__}")
-        print()  # 打印空行以
-
-# 设定NetCDF文件的路径
-path = 'your_file.nc'  # 请替换为您的实际文件路径
-
-with nc.Dataset(path) as dataNC:
-    time_data = dataNC.variables['time'][:]    
-    # 打印时间数据
-    print("时间数据:")
-    print(time_data)
-    # 如果需要转换为可读的日期时间格式
-    time_units = dataNC.variables['time'].units
-    print(f"时间单位: {time_units}")
-
-    # 如果单位是 "days since YYYY-MM-DD" 形式，进行转换
-    if "days since" in time_units:
-        base_time = nc.num2date(time_data, units=time_units)
-        print("转换后的时间数据:")
-        for t in base_time:
-            print(t)
-
-"""
-

shancx/radar_nmc_china_map_compare1.py

@@ -1,50 +0,0 @@
-import matplotlib.pyplot as plt
-import numpy as np
-import datetime
-from hjnwtx.colormap import cmp_hjnwtx  # 假设这是您的自定义颜色映射库
-import cartopy.crs as ccrs
-import cartopy.feature as cfeature
-import cartopy.io.shapereader as shpreader
-from shancx import crDir
-def drawpic_com(array_dt,array_dt1,name="temp"):
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/{name}_{now_str}.png"
-    crDir(outpath)
-
-    # 创建绘图和设置坐标系
-    fig = plt.figure(figsize=(10, 8))
-    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
-    
-    # 设置图像显示的范围
-    ax.set_extent([73, 135, 18, 54], ccrs.PlateCarree())  # 根据需要调整
-
-    # 添加中国地图的边界和特征，包括省份轮廓
-    add_china_map(ax)
-    
-    # 添加数据层
-    if len(array_dt.shape) == 3:
-        for i, img_ch_nel in enumerate(array_dt):
-            ax.imshow(img_ch_nel, vmin=50, vmax=500, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=[73, 134.99, 12.21, 54.2])
-            plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-            plt.savefig(f"{outpath}_layer_{i}.png")
-            plt.clf()  # 清除图形以绘制下一个通道图像
-    elif len(array_dt.shape) == 2:
-        ax.imshow(array_dt, vmin=0, vmax=100, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=[73, 134.99, 12.21, 54.2], alpha=1)
-        ax.imshow(array_dt1, vmin=0, vmax=100, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=[73, 134.99, 12.21, 54.2], alpha=0.3)
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    plt.close(fig)
-
-def add_china_map(ax):
-    # 在地图上添加地形特征
-    ax.add_feature(cfeature.COASTLINE, edgecolor='gray')
-    ax.add_feature(cfeature.BORDERS, linestyle=':', edgecolor='gray')
-    ax.add_feature(cfeature.LAKES, alpha=0.8)
-    # 添加省份轮廓
-    provinces = shpreader.natural_earth(resolution='10m', category='cultural', name='admin_1_states_provinces')
-    provinces_features = shpreader.Reader(provinces).geometries()
-    ax.add_geometries(provinces_features, ccrs.PlateCarree(), facecolor='none', edgecolor='gray', linestyle=':', linewidth=0.5, alpha=0.8)
-
- 
-
-# 示例用法

shancx/radar_nmc_china_map_f.py

@@ -1,125 +0,0 @@
-import matplotlib.pyplot as plt
-import numpy as np
-import datetime
-from hjnwtx.colormap import cmp_hjnwtx  # 假设这是您的自定义颜色映射库
-import cartopy.crs as ccrs
-import cartopy.feature as cfeature
-import cartopy.io.shapereader as shpreader
-import os
-from shancx import crDir
-
-def drawCN(array_dt,ty="CR"):
-    env = [73, 135, 18, 54]   #US  [-119, -64, 22, 50]
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/CN_{now_str}.png"
-    crDir(outpath)
-
-    # 创建绘图和设置坐标系
-    fig = plt.figure(figsize=(10, 8))
-    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
-    
-    # 设置图像显示的范围
-    ax.set_extent(env, ccrs.PlateCarree())  # 根据需要调整    [-119, -64, 22, 50]
-
-    # 添加中国地图的边界和特征，包括省份轮廓
-    add_china_map(ax)
-    
-    # 添加数据层
-    if len(array_dt.shape) == 3:
-        for i, img_ch_nel in enumerate(array_dt):
-            ax.imshow(img_ch_nel, vmin=50, vmax=500, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-            plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-            plt.savefig(f"{outpath}_layer_{i}.png")
-            plt.clf()  # 清除图形以绘制下一个通道图像
-    elif len(array_dt.shape) == 2 and ty=="pre":
-        # ax.imshow(array_dt, vmin=0, vmax=100, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=[-119, -64, 22, 50])   #pre_tqw
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    else :
-        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-        # ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=[-119, -64, 22, 50])   #pre_tqw
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    plt.close(fig)
-
-def add_china_map(ax):
-    # 在地图上添加地形特征
-    ax.add_feature(cfeature.COASTLINE, edgecolor='gray')
-    ax.add_feature(cfeature.BORDERS, linestyle=':', edgecolor='gray')
-    ax.add_feature(cfeature.LAKES, alpha=0.8)
-    # 添加省份轮廓
-    if os.path.exists('/home/scx/ne_10m_admin_1_states_provinces.shp'):
-        provinces = '/home/scx/ne_10m_admin_1_states_provinces.shp'
-    else :
-        provinces = shpreader.natural_earth(resolution='10m', category='cultural', name='admin_1_states_provinces')   ###自动下载 
-    provinces_features = shpreader.Reader(provinces).geometries()
-    """
-    本地路径读取 #'/home/scx/.local/share/cartopy/shapefiles/natural_earth/cultural/ne_10m_admin_1_states_provinces.shp'   #ne_10m_admin_1_states_provinces.shx 
-    # shapefile_path = '/home/scx/.local/share/cartopy/shapefiles/natural_earth/cultural/admin_1_states_provinces.shp'
-    # provinces_features = shpreader.Reader(shapefile_path).geometries()
-    """
-    ax.add_geometries(provinces_features, ccrs.PlateCarree(), facecolor='none', edgecolor='gray', linestyle=':', linewidth=0.5, alpha=0.8)
-
- 
-# 示例用法
-
-
-
-import matplotlib.pyplot as plt
-import numpy as np
-import datetime
-from hjnwtx.colormap import cmp_hjnwtx  # 假设这是您的自定义颜色映射库
-import cartopy.crs as ccrs
-import cartopy.feature as cfeature
-import cartopy.io.shapereader as shpreader
-from shancx import crDir
-import os
-
-def drawCNXY(array_dt,XY = [3,2]):  #x_coords2,y_coords2
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/{now_str}.png"
-    crDir(outpath)
-
-    # 创建绘图和设置坐标系
-    fig = plt.figure(figsize=(10, 8))
-    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
-    
-    # 设置图像显示的范围
-    ax.set_extent([73, 135, 18, 54], ccrs.PlateCarree())  # 根据需要调整
-
-    # 添加中国地图的边界和特征，包括省份轮廓
-    add_china_map(ax)    
-    # 添加数据层
-    if len(array_dt.shape) == 3:
-        for i, img_ch_nel in enumerate(array_dt):
-            ax.imshow(img_ch_nel, vmin=50, vmax=500, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=[73, 134.99, 12.21, 54.2])
-            plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-            plt.savefig(f"{outpath}_layer_{i}.png")
-            plt.clf()  # 清除图形以绘制下一个通道图像
-    elif len(array_dt.shape) == 2:
-        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=[73, 134.99, 12.21, 54.2])
-        x_coords2 = XY[0]
-        y_coords2 = XY[1]
-        plt.plot(list(x_coords2), list(y_coords2), 'ro', markersize=1)
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    plt.close(fig)
-
-def add_china_map(ax):
-    # 在地图上添加地形特征
-    ax.add_feature(cfeature.COASTLINE, edgecolor='gray')
-    ax.add_feature(cfeature.BORDERS, linestyle=':', edgecolor='gray')
-    ax.add_feature(cfeature.LAKES, alpha=0.8)
-    # 添加省份轮廓
-    if os.path.exists('/home/scx/ne_10m_admin_1_states_provinces.shp'):
-        provinces = '/home/scx/ne_10m_admin_1_states_provinces.shp'
-    else :
-        provinces = shpreader.natural_earth(resolution='10m', category='cultural', name='admin_1_states_provinces')   ###自动下载 
-    provinces_features = shpreader.Reader(provinces).geometries()
-    """
-    本地路径读取 #'/home/scx/.local/share/cartopy/shapefiles/natural_earth/cultural/ne_10m_admin_1_states_provinces.shp'   #ne_10m_admin_1_states_provinces.shx 
-    # shapefile_path = '/home/scx/.local/share/cartopy/shapefiles/natural_earth/cultural/admin_1_states_provinces.shp'
-    # provinces_features = shpreader.Reader(shapefile_path).geometries()
-    """
-    ax.add_geometries(provinces_features, ccrs.PlateCarree(), facecolor='none', edgecolor='gray', linestyle=':', linewidth=0.5, alpha=0.8)

shancx/radar_nmc_us_map_f.py

@@ -1,67 +0,0 @@
-import matplotlib.pyplot as plt
-import numpy as np
-import datetime
-from hjnwtx.colormap import cmp_hjnwtx  # Assuming this is your custom colormap library
-import cartopy.crs as ccrs
-import cartopy.feature as cfeature
-import cartopy.io.shapereader as shpreader
-import os
-from shancx import crDir
-
-def drawUS(array_dt,ty="CR"):
-    # env = [-119, -64, 22, 50]
-    env = [-132.0, -47.0, 0, 57.0]
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/US_{now_str}.png"
-    crDir(outpath)
-
-    # Create figure and set the coordinate system
-    fig = plt.figure(figsize=(10, 8))
-    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
-
-    # Set the extent for the United States
-    ax.set_extent(env, ccrs.PlateCarree())  # Adjust as needed
-
-    # Add the US map boundaries and features
-    add_us_map(ax)
-    
-    # Add data layers
-    if len(array_dt.shape) == 3:
-        for i, img_ch_nel in enumerate(array_dt):
-            ax.imshow(img_ch_nel, vmin=50, vmax=500, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-            plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-            plt.savefig(f"{outpath}_layer_{i}.png")
-            plt.clf()  # Clear the figure to draw the next channel image
-    elif len(array_dt.shape) == 2 and ty =="pre":
-        # ax.imshow(array_dt, vmin=0, vmax=100, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)   #pre_tqw
-        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   #pre_tqw
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    else :
-        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-        # ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=[-119, -64, 22, 50])   #pre_tqw
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-
-    plt.close(fig)
-
-def add_us_map(ax):
-    # Add geographic features for the US
-    ax.add_feature(cfeature.COASTLINE, edgecolor='gray')
-    ax.add_feature(cfeature.BORDERS, linestyle=':', edgecolor='gray')
-    ax.add_feature(cfeature.LAKES, alpha=0.8)
-    
-    # Adding state boundaries
-    if os.path.exists('/home/scx/ne_10m_admin_1_states_provinces.shp'):
-        states = '/home/scx/ne_10m_admin_1_states_provinces.shp'
-    else:
-        states = shpreader.natural_earth(resolution='10m', category='cultural', name='admin_1_states_provinces')  # Automatically download 
-
-    states_features = shpreader.Reader(states).geometries()
-    
-    ax.add_geometries(states_features, ccrs.PlateCarree(), facecolor='none', edgecolor='gray', linestyle=':', linewidth=0.5, alpha=0.8)
-
-# Example usage
-# Assuming array_dt is your data array, pass it to drawUS
-# array_dt = np.random.rand(3, 100, 100)  # Example random data; replace with your actual data
-# drawUS(array_dt)

shancx/subplots_compare_devlop.py

@@ -1,36 +0,0 @@
-
-import glob
-import numpy as np
-from hjnwtx.colormap import cmp_hjnwtx
-PathList = glob.glob("/mnt/wtx_weather_forecast/scx/testDBscan/2024_p/20240312/*")
-path = PathList[0]
-print(path) 
-
-import numpy as np 
-import matplotlib.pyplot as plt
-x = np.load(path)
-# for i in range(35):
-    # plt.imshow(x[i])
-    # plt.savefig(f"./image/a{str(i)}.png")
-    # plt.show()
-
-
-from hjnwtx.colormap import cmp_hjnwtx
-import matplotlib.pyplot as plt
-import datetime
-def drawpic_com(base_up,base_down,shape_len, name="temp"): 
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    Count = shape_len
-    data_all = np.concatenate([base_up,base_down], axis=0)
-    fig, axs = plt.subplots(2, Count, figsize=(10*shape_len, 10))
-    # data_all = data_all * 70
-    for i in range(2):
-        for j in range(Count):
-            index = i * Count + j
-            axs[i, j].imshow(data_all[index, :, :],vmax=70,vmin=0,cmap=cmp_hjnwtx["radar_nmc"])
-            axs[i, j].axis('off')
-    plt.tight_layout()
-    plt.savefig(name +now_str + '.png')
-    plt.close() 
-#drawpic(x,x,35)
-print("done")

shancx/tool.py

@@ -1,18 +0,0 @@
-
-import netCDF4 as nc
-import numpy as np
-def getPoint(pre, df, lat0, lon0, resolution, decimal=1):
-    latIdx = ((lat0 - df["Lat"]) / resolution + 0.5).astype(np.int64)
-    lonIdx = ((df["Lon"] - lon0) / resolution + 0.5).astype(np.int64)
-    return pre[...,latIdx, lonIdx].round(decimals=decimal)
-def Get_Lat_Lon_QPF(path,Lon_data,Lat_data):
-    with nc.Dataset(path) as dataNC:
-        latArr = dataNC["lat"][:]
-        lonArr = dataNC["lon"][:]
-        if "AIW_QPF" in  path:
-            pre = dataNC[list(dataNC.variables.keys())[3]][:]    
-        elif "AIW_REF" in path:
-            pre = dataNC[list(dataNC.variables.keys())[4]][:]   
-    data = getPoint(pre , {"Lon":Lon_data,"Lat":Lat_data} , latArr[0], lonArr[0], 0.01)
-    data = getPoint(pre , {"Lon":Lon_data,"Lat":Lat_data} , latArr[0], lonArr[0], 0.01)
-    return data