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

shancx/Path1.py

@@ -0,0 +1,161 @@
+
+    
+import os 
+from pathlib import Path
+def curPath_():
+    current_file_path = os.path.abspath(__file__) 
+    current_folder_path = os.path.dirname(current_file_path)
+    parent = Path(__file__).parent
+    return current_folder_path,parent
+
+import sys
+def curPath():
+    # 获取当前执行文件的绝对路径
+    current_file_path = os.path.abspath(sys.argv[0])
+    current_folder_path = os.path.dirname(current_file_path)
+    return current_folder_path
+
+from pathlib import Path 
+def RootFilePaths(root_path=None):
+    root_dir = Path(f'{root_path}')
+    npy_files = list(root_dir.rglob('*.*'))
+    return npy_files
+""" 
+import gzip
+with gzip.open(file_path=None, 'rt') as file:
+    first_line = file.readline()
+"""
+
+"""
+CSV
+
+import gzip
+import csv
+with gzip.open(file_path, 'rt') as file:
+    reader = csv.reader(file)
+    for row in reader:
+        print(row)
+
+XML
+
+import gzip
+import xml.etree.ElementTree as ET
+with gzip.open(file_path, 'rt') as file:
+    tree = ET.parse(file)
+    root = tree.getroot()
+    for child in root:
+        print(child.tag, child.attrib)
+
+YAML        
+
+import gzip
+import yaml
+with gzip.open(file_path, 'rt') as file:
+    data = yaml.safe_load(file)
+    print(data)
+
+TEXT
+
+import gzip
+with gzip.open(file_path, 'rt') as file:
+    for line in file:
+        print(line.strip())
+
+"""
+ 
+"""   model 2
+import matplotlib.pyplot as plt
+import datetime
+import time
+import os
+import time
+import shutil 
+from dateutil.relativedelta import relativedelta
+import glob
+import argparse
+from multiprocessing import Pool
+from itertools import product
+import pandas as pd  
+import numpy as np
+import copy
+import traceback  
+from shancx.NN import _loggers
+logger = _loggers()
+from shancx import crDir,Mul_sub
+
+base_source_dir = "/mnt/wtx_weather_forecast/CMA_DATA/NAFP/EC/C1D/2023/2023123120/"
+base_source_dir = "/mnt/wtx_weather_forecast/CMA_DATA/NAFP/EC/C1D"
+# 目标文件夹路径
+target_dir = "/mnt/wtx_weather_forecast/scx/CMA_DATA/NAFP/EC/C1D/2023/2023123120/"
+base_target_dir = "/mnt/wtx_weather_forecast/scx/CMA_DATA/NAFP/EC/C1D"
+
+def map_data(conf):
+    UCT = conf[0]
+    # UCT = CST + relativedelta(hours=-8)
+    UCTstr = UCT.strftime("%Y%m%d%H%M") 
+    # 文件名
+    file_name = "ECMFC1D_WIV_100"
+    source_dir = f"{base_source_dir}/{UCTstr[:4]}/{UCTstr[:10]}/*"
+    target_dir = f"{base_target_dir}/{UCTstr[:4]}/{UCTstr[:10]}"
+    crDir(target_dir)
+    source_fileL = [i for i in glob.glob(source_dir) if "ECMFC1D_WIV_100"  in  i ]
+    if not source_fileL:
+        return
+    source_file = source_fileL[0]
+    filename = os.path.basename(source_file)
+    target_file = f"{target_dir}/{filename}"
+    if os.path.exists(target_file):
+        print(f"目标文件已存在: {target_file}")
+    else:
+        try:
+            shutil.copy(source_file, target_dir)
+            print(f"文件已成功复制到: {target_file}")
+            logger.info(f"文件已成功复制到: {source_file} ---> {target_file}")
+        except Exception as e:
+            print(f"复制文件时出错: {traceback.format_exc()}")
+            logger.error(f"复制文件时出错: {traceback.format_exc()}")
+def options():
+    parser = argparse.ArgumentParser(description='scx')
+    # parser.add_argument('--times', type=str, default='202406290000,202406300000')
+    parser.add_argument('--times', type=str, default='202407210000,202407220000')
+    parser.add_argument('--pac', type=str, default='100000')
+    # parser.add_argument('--combine', action='store_true', default=False)
+    parser.add_argument('--combine',action='store_true',default=False)
+    parser.add_argument('--isDebug',action='store_true',default=False)
+    parser.add_argument('--isDraw',action='store_true',default=False)
+    parser.add_argument('--freq', type=str, default="1H")
+    parser.add_argument('--tag',type=str, default=datetime.datetime.now().strftime("%Y%m%d%H%M"))
+    config= parser.parse_args()
+    print(config)
+    config.times = config.times.split(",")
+    config.pac = config.pac.split(",")
+    if len(config.times) == 1:
+        config.times = [config.times[0], config.times[0]]
+    config.times = [datetime.datetime.strptime(config.times[0], "%Y%m%d%H%M"),
+                    datetime.datetime.strptime(config.times[1], "%Y%m%d%H%M")]
+    return config
+if __name__ == '__main__':
+    cfg = options()
+    sUTC = cfg.times[0] + relativedelta(hours=-8)
+    eUTC = cfg.times[-1] + relativedelta(hours=-8)
+    sUTCstr = sUTC.strftime("%Y%m%d")
+    cfg = options()
+    # 假设 cfg.times[0] 和 cfg.times[-1] 是 datetime 对象 
+    # 生成每天的时间序列
+    date_range = pd.date_range(start=sUTC.date(), end=eUTC.date(), freq="D")
+    # 生成每天 08:00 和 20:00 的时间点
+    timeList = []
+    for date in date_range:
+        timeList.append(date.replace(hour=8, minute=0, second=0, microsecond=0))
+        timeList.append(date.replace(hour=20, minute=0, second=0, microsecond=0))
+
+    # 过滤超出范围的时间点
+    timeList1 = [t for t in timeList if sUTC <= t <= eUTC]
+    logger.info(timeList1)
+    # for CST in timeList: 
+    try: 
+        Mul_sub(map_data,[timeList1],31)
+    except Exception as e:
+        print(traceback.format_exc())
+"""
+

shancx/Plot/GlobMap.py

@@ -1,116 +1,276 @@
-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 Glob(array_dt,cd="CHN",ty="CR"):
-    if cd == "g":
-        env = [-179.617020, 179.632979,-85.098871,85.051128] 
-    elif cd == "US":
-        env = [-132.0, -47.0, 0, 57.0]
-    elif cd == "CHN":
-        env = [73,134.99,12.21,54.2]   
-    else:
-         env = cd
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/{str(cd)}_{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_glob_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=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
-        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)
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    plt.close(fig)
-
-def add_glob_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)    
-    # 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)
- 
-def GlobLonLat(array_dt,Lon,Lat,cd="CHN",ty="CR"):   ###  x_coords2 维度 
-    if cd == "g":
-        env = [-179.617020, 179.632979,-85.098871,85.051128] 
-    elif cd == "US":
-        env = [-132.0, -47.0, 0, 57.0]
-    elif cd == "CHN":
-        env = [73,134.99,12.21,54.2]   #[73, 134.99, 12.21, 54.2] 
-    else:
-         env = cd    
-    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
-    outpath = f"./radar_nmc/{str(cd)}_{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_glob_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.close()  
-    elif len(array_dt.shape) == 2 and ty =="pre":
-        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        ax.scatter(list(Lon), list(Lat), s=0.5, c='red', marker='o', transform=ccrs.PlateCarree())
-        plt.savefig(outpath)
-    else :
-        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
-        ax.scatter(list(Lon), list(Lat), s=0.5, c='red', marker='o', transform=ccrs.PlateCarree())
-        plt.colorbar(ax.images[0], ax=ax, orientation='vertical')
-        plt.savefig(outpath)
-    plt.close(fig)
-
-def add_us_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)    
-    # 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)
+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 Glob(array_dt=None,cd="CHN",ty="CR",colorbarflag="add"):
+    if cd == "g":
+        env = [-179.617020, 179.632979,-85.098871,85.051128] 
+    elif cd == "US":
+        env = [-132.0, -47.0, 0, 57.0]
+    elif cd == "CHN":
+        env = [73,134.99,12.21,54.2]   
+    else:
+         env = cd
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
+    outpath = f"./radar_nmc/{str(cd)}_{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_glob_map(ax)    
+    # Add data layers
+    if  len(array_dt.shape) == 2 and ty =="pre":
+        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
+        if colorbarflag is not None :
+           plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)
+        plt.savefig(outpath)
+    else :
+        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
+        if colorbarflag is not None :
+           plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)
+        plt.savefig(outpath)
+    plt.close(fig)
+
+def add_glob_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)    
+    # 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)
+ 
+def GlobLonLat(array_dt=None,Lon=None,Lat=None,cd="CHN",ty="CR",colorbarflag="add"):   ###  x_coords2 维度 
+    if cd == "g":
+        env = [-179.617020, 179.632979,-85.098871,85.051128] 
+    elif cd == "US":
+        env = [-132.0, -47.0, 0, 57.0]
+    elif cd == "CHN":
+        env = [73,134.99,12.21,54.2]   #[73, 134.99, 12.21, 54.2] 
+    else:
+         env = cd    
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
+    if array_dt is None :
+        array_dt = np.full([4200,6200],np.nan)
+    outpath = f"./radar_nmc/{str(cd)}_{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_glob_map(ax)
+    
+    # Add data layers
+    if  len(array_dt.shape) == 2 and ty =="pre":
+        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
+        if colorbarflag is not None :
+           plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)
+        ax.scatter(list(Lon), list(Lat), s=0.5, c='red', marker='o', transform=ccrs.PlateCarree())
+    else :
+        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
+        ax.scatter(list(Lon), list(Lat), s=0.5, c='red', marker='o', transform=ccrs.PlateCarree())
+        if colorbarflag is not None :
+           plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)        
+    plt.tight_layout()
+    plt.savefig(outpath)
+    plt.close(fig)
+
+def add_us_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)    
+    # 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)
+
+def GlobLonLatPlus(array_dt=None,Lon=None,Lat=None,cd="CHN",ty="CR",name="temp",markcolor="red",outpath1=None):   ###  x_coords2 维度 
+    if cd == "g":
+        env = [-179.617020, 179.632979,-85.098871,85.051128] 
+    elif cd == "US":
+        env = [-132.0, -47.0, 0, 57.0]
+    elif cd == "CHN":
+        env = [73,134.99,12.21,54.2]   #[73, 134.99, 12.21, 54.2] 
+    else:
+         env = cd    
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
+    if array_dt is None :
+        array_dt = np.full([4200,6200],np.nan)
+    outpath = f"./radar_nmc/{str(cd)}_{now_str}_{name}.png"
+    if outpath1 is not None:
+       outpath = outpath1  
+    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_glob_map(ax)
+    
+    # Add data layers
+    if  len(array_dt.shape) == 2 and ty =="pre":
+        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
+        ax.scatter(list(Lon), list(Lat), s=0.5, c=f'{markcolor}', marker='o', transform=ccrs.PlateCarree())
+        plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)
+        plt.savefig(outpath)
+    else :
+        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
+        ax.scatter(list(Lon), list(Lat), s=0.5, c=f'{markcolor}', marker='o', transform=ccrs.PlateCarree())
+        plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.7)
+        plt.savefig(outpath)
+    plt.close(fig)
+
+def add_us_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)    
+    # 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)
+
+import matplotlib.font_manager as fm
+zh_font_path = '/mnt/wtx_weather_forecast/scx/fonts/truetype/simhei.ttf'   
+zh_font = fm.FontProperties(fname=zh_font_path)
+def GlobLonLatPluss(array_dt,Lon,Lat,Lon1,Lat1,cd="CHN",ty="CR",name="temp",markcolor="red",outpath1=None):   ###  x_coords2 维度 
+    if cd == "g":
+        env = [-179.617020, 179.632979,-85.098871,85.051128] 
+    elif cd == "US":
+        env = [-132.0, -47.0, 0, 57.0]
+    elif cd == "CHN":
+        env = [73,134.99,12.21,54.2]   #[73, 134.99, 12.21, 54.2] 
+    else:
+         env = cd    
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M%S")
+    outpath = f"./radar_nmc/{str(cd)}_{now_str}_{name}.png"
+    array_dt = np.where(array_dt>0,array_dt,np.nan)
+    if outpath1 is not None:
+       outpath = outpath1  
+    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_glob_map(ax)    
+    # Add data layers
+    if  len(array_dt.shape) == 2 and ty =="pre":
+        ax.imshow(array_dt, vmin=0, vmax=10, cmap=cmp_hjnwtx["pre_tqw"], transform=ccrs.PlateCarree(), extent=env)   
+        ax.scatter(list(Lon), list(Lat), s=0.5, c=f'{markcolor}', marker='o', transform=ccrs.PlateCarree())
+    else :
+        ax.imshow(array_dt, vmin=0, vmax=72, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=env)
+        sc1 = ax.scatter(list(Lon), list(Lat), s=2, c=f'blue', marker='o', transform=ccrs.PlateCarree(),edgecolor='black', linewidth=0.3,label="维天信相对彩云漏报 ")
+        sc2 = ax.scatter(list(Lon1), list(Lat1), s=2, c=f'red', marker='o', transform=ccrs.PlateCarree(),edgecolor='black', linewidth=0.3,label="共同漏报 ")
+        plt.colorbar(ax.images[0], ax=ax, orientation='vertical',shrink=0.8)
+        plt.title(f"{name}")
+        ax.legend(loc='upper left', fontsize='small', prop=zh_font)  # 您可以根据需要调整字体大小和位置
+    plt.savefig(outpath)
+    plt.close(fig)
+
+def add_us_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)    
+    # 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)
+
+def globPng(d2,ty=None):
+    import numpy as np
+    import matplotlib.pyplot as plt
+    import cartopy.crs as ccrs
+    import cartopy.feature as cfeature
+    from hjnwtx.colormap import cmp_hjnwtx
+    from shancx.Time import UTCStr
+    CSTstr = UTCStr()
+    data = d2[::10,::10]  # shape: (height, width)
+    lat_min, lat_max = -90, 90
+    lon_min, lon_max = -180, 180
+    fig = plt.figure(figsize=(15, 10))
+    ax = fig.add_subplot(111, projection=ccrs.PlateCarree())
+    im = ax.imshow(data, 
+                  extent=[lon_min, lon_max, lat_min, lat_max],
+                  cmap=cmp_hjnwtx["radar_nmc"],
+                  transform=ccrs.PlateCarree())
+    if ty == "pre":
+       im = ax.imshow(data, 
+                    extent=[lon_min, lon_max, lat_min, lat_max],
+                    vmin=0, 
+                    vmax=10,
+                    cmap=cmp_hjnwtx["pre_tqw"],
+                    transform=ccrs.PlateCarree())
+    ax.add_feature(cfeature.COASTLINE, linewidth=0.8, zorder=2)
+    ax.add_feature(cfeature.BORDERS, linestyle=':', linewidth=0.6, zorder=2)
+    ax.add_feature(cfeature.STATES, linestyle=':', linewidth=0.4, alpha=0.5, zorder=2)
+    ax_pos = ax.get_position()    
+    # 设置色标与主图严格对齐
+    cax = fig.add_axes([
+        ax_pos.x1 + 0.01,  # 主图右侧+0.01的间距
+        ax_pos.y0,         # 与主图底部对齐
+        0.02,              # 色标宽度
+        ax_pos.height      # 与主图同高度
+    ])
+    
+    cbar = plt.colorbar(im, cax=cax)
+    cbar.set_label('Echo Intensity (dBZ)', fontsize=12, labelpad=10)
+    cbar.ax.tick_params(labelsize=10)
+    plt.subplots_adjust(
+        left=0.05, 
+        right=0.85,  # 减小right值以容纳色标
+        top=0.95,
+        bottom=0.125
+    )
+    ax.set_title('Global Radar', 
+                fontsize=18, 
+                pad=25, 
+                weight='bold')
+    ax.set_xlabel('Longitude', fontsize=12, labelpad=10)
+    ax.set_ylabel('Latitude', fontsize=12, labelpad=15)
+    plt.savefig(f"global_radar_aligned_{CSTstr}.png", 
+               dpi=300, 
+               bbox_inches='tight', 
+               facecolor='white')
+    plt.close()