shancx 1.60__py3-none-any.whl → 1.70__py3-none-any.whl

shancx/__init__.py

@@ -52,3 +52,96 @@ def Mul_Ess(map_fun,tuple_list):
     with Pool(31) as p:
         P_data = p.map(map_fun, tuple_list)
     return P_data
+
+
+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/mkIMGSCX.py

@@ -0,0 +1,92 @@
+
+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/radar_nmc_china_map_compare1.py

@@ -0,0 +1,55 @@
+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
+
+def big_value_fun(array_dt,array_dt1,name="temp"):
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M")
+    outpath = f"./radar_nmc/{name}_{now_str}.png"
+    mkDir(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)
+
+def mkDir(path):
+    # 处理目录创建
+    import os
+    directory = os.path.dirname(path)
+    if not os.path.exists(directory):
+        os.makedirs(directory)
+
+# 示例用法

shancx/radar_nmc_china_map_f.py

@@ -0,0 +1,54 @@
+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
+
+def big_value_fun(array_dt):
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M")
+    outpath = f"./radar_nmc/{now_str}.png"
+    mkDir(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])
+        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)
+
+def mkDir(path):
+    # 处理目录创建
+    import os
+    directory = os.path.dirname(path)
+    if not os.path.exists(directory):
+        os.makedirs(directory)
+
+# 示例用法

shancx/subplots_compare_devlop.py

@@ -0,0 +1,36 @@
+
+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(base_up,base_down,shape_len, name="temp"): 
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M")
+    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/subplots_single_china_map.py

@@ -0,0 +1,40 @@
+import matplotlib.pyplot as plt
+import cartopy.crs as ccrs
+import cartopy.feature as cfeature
+import cartopy.io.shapereader as shpreader
+import datetime
+import os
+
+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)
+
+def mkDir(filepath):
+    directory = os.path.dirname(filepath)
+    if not os.path.exists(directory):
+        os.makedirs(directory)
+
+def draw_radar_images(data, num_images, output_name="temp"):
+    # Validate input
+    if not data.shape[0] >= num_images:
+        raise ValueError("Data does not contain enough entries for the number of images specified.")
+
+    # Setup figure dimensions based on the number of images
+    fig, axs = plt.subplots(1, num_images, figsize=(10 * num_images, 10), subplot_kw={'projection': ccrs.PlateCarree()})
+    
+    for i in range(num_images):
+        axs[i].imshow(data[i, :, :], vmax=70, vmin=0, cmap=cmp_hjnwtx["radar_nmc"], transform=ccrs.PlateCarree(), extent=[73, 135, 18, 54])
+        add_china_map(axs[i])
+        axs[i].axis('off')
+    
+    # Save the figure
+    now_str = datetime.datetime.now().strftime("%Y%m%d%H%M")
+    outpath = f"./aplot_image/{output_name}_{now_str}.png"
+    mkDir(outpath)
+    fig.tight_layout()
+    fig.savefig(outpath, dpi=300)
+    plt.close(fig)

{shancx-1.60.dist-info → shancx-1.70.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: shancx
-Version: 1.60
+Version: 1.70
 Summary: A simple timer decorator
 Home-page: https://gitee.com/shancx
 Author: Your Name

shancx-1.70.dist-info/RECORD

@@ -0,0 +1,12 @@
+my_timer_decorator/__init__.py,sha256=ur0N4vwlv-LfYOCph5A8g6H8aFxxc9EhJdG9sxpa-SY,303
+shancx/__init__.py,sha256=2Y16Cyjk7CxtRJllGNrfjypvDuKZubviN1FNODuVUv4,4312
+shancx/mkIMGSCX.py,sha256=RN0LqYcZ3W8ez4jE5F3KBApw9biNgXVXeFuNVeZvQ2U,2520
+shancx/radar_nmc_china_map_compare1.py,sha256=3PhtDm6sAY8EjmpdT8_TXJYjeGbqi02ce8QeqSMp76s,2493
+shancx/radar_nmc_china_map_f.py,sha256=vgxYKygN1SibHVW4wU2zdlRrBm1r4olo_ln0PZr0178,2302
+shancx/subplots_compare_devlop.py,sha256=ZsWEpv-z9ILPBE-mkleSIvLFx2UbRZyGySiWi7PX8oc,1075
+shancx/subplots_single_china_map.py,sha256=NoxmZoc3o7-FIAGYM-PU2jn75EAdq_7k70lJ-PGY44E,1732
+shancx/util_log.py,sha256=bGu-oCU7hozBBZJB6oDGj9XDDsz5GJVe8ZfJtx5pddA,1204
+shancx-1.70.dist-info/METADATA,sha256=STAvvfrd43RZHkzWpC04zpuUDo6QHipNsmmqK6b3TOM,572
+shancx-1.70.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+shancx-1.70.dist-info/top_level.txt,sha256=akfCS1vKWz3pNmEN_yN9ZiGp-60IQY5ET38mRx_i_-4,7
+shancx-1.70.dist-info/RECORD,,

shancx-1.60.dist-info/RECORD

@@ -1,7 +0,0 @@
-my_timer_decorator/__init__.py,sha256=ur0N4vwlv-LfYOCph5A8g6H8aFxxc9EhJdG9sxpa-SY,303
-shancx/__init__.py,sha256=7E61iUEkmy21JSSyOk6kTRGB9yGJo5dJ5cMqHKqCW4U,1790
-shancx/util_log.py,sha256=bGu-oCU7hozBBZJB6oDGj9XDDsz5GJVe8ZfJtx5pddA,1204
-shancx-1.60.dist-info/METADATA,sha256=Ymku04OYpWTcfJcGi68WJEeG9rQZP_5ZobTgcQU6jXA,572
-shancx-1.60.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-shancx-1.60.dist-info/top_level.txt,sha256=akfCS1vKWz3pNmEN_yN9ZiGp-60IQY5ET38mRx_i_-4,7
-shancx-1.60.dist-info/RECORD,,