shancx 1.60__tar.gz → 1.70__tar.gz

{shancx-1.60 → shancx-1.70}/PKG-INFO

@@ -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.60 → shancx-1.70}/setup.py

@@ -4,7 +4,7 @@ from setuptools import setup, find_packages
 
 setup(
     name='shancx',
-    version='1.60',
+    version='1.70',
     packages=find_packages(),
     description='A simple timer decorator',
     long_description=open('README.md').read(),
@@ -19,3 +19,4 @@ setup(
     ],
     python_requires='>=3.6',
 )
+ 

shancx-1.70/shancx/__init__.py

@@ -0,0 +1,147 @@
+#!/usr/bin/python
+# -*- coding: utf-8 -*-
+# @Time : 2023/09/27 下午8:52
+# @Author : shanchangxi
+# @File : util_log.py
+import time
+import logging  # 引入logging模块
+from logging import handlers
+# from util.util_dir import *
+
+# 第一步，创建一个logger
+loggers = logging.getLogger()
+loggers.setLevel(logging.INFO)  # Log等级开关
+
+# 第二步，创建一个handler，用于写入日志文件
+log_name =  'project_tim_tor.log'
+logfile = log_name
+
+time_rotating_file_handler = handlers.TimedRotatingFileHandler(filename=logfile, when='D', encoding='utf-8')
+time_rotating_file_handler.setLevel(logging.INFO)  # 输出到file的log等级的开关
+
+# 第三步，定义handler的输出格式
+formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
+time_rotating_file_handler.setFormatter(formatter)
+
+# 第四步，将handler添加到logger里面
+loggers.addHandler(time_rotating_file_handler)
+
+# 如果需要同時需要在終端上輸出，定義一個streamHandler
+print_handler = logging.StreamHandler()  # 往屏幕上输出
+print_handler.setFormatter(formatter)  # 设置屏幕上显示的格式
+loggers.addHandler(print_handler)
+
+
+def Tim_Tor(func):
+    def wrapper(*args, **kwargs):
+        start_time = time.time()
+        result = func(*args, **kwargs)
+        end_time = time.time()
+        print(f"{func.__name__} took {end_time - start_time:.4f} seconds")
+        loggers.info(f"{func.__name__} took {end_time - start_time:.4f} seconds")
+        return result
+    return wrapper
+
+
+
+
+
+from itertools import product
+from multiprocessing import Pool
+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-1.70/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-1.70/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-1.70/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-1.70/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-1.70/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 → shancx-1.70}/shancx.egg-info/PKG-INFO

@@ -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/shancx.egg-info/SOURCES.txt

@@ -0,0 +1,12 @@
+README.md
+setup.py
+shancx/__init__.py
+shancx/mkIMGSCX.py
+shancx/radar_nmc_china_map_compare1.py
+shancx/radar_nmc_china_map_f.py
+shancx/subplots_compare_devlop.py
+shancx/subplots_single_china_map.py
+shancx.egg-info/PKG-INFO
+shancx.egg-info/SOURCES.txt
+shancx.egg-info/dependency_links.txt
+shancx.egg-info/top_level.txt

shancx-1.60/shancx/__init__.py

@@ -1,54 +0,0 @@
-#!/usr/bin/python
-# -*- coding: utf-8 -*-
-# @Time : 2023/09/27 下午8:52
-# @Author : shanchangxi
-# @File : util_log.py
-import time
-import logging  # 引入logging模块
-from logging import handlers
-# from util.util_dir import *
-
-# 第一步，创建一个logger
-loggers = logging.getLogger()
-loggers.setLevel(logging.INFO)  # Log等级开关
-
-# 第二步，创建一个handler，用于写入日志文件
-log_name =  'project_tim_tor.log'
-logfile = log_name
-
-time_rotating_file_handler = handlers.TimedRotatingFileHandler(filename=logfile, when='D', encoding='utf-8')
-time_rotating_file_handler.setLevel(logging.INFO)  # 输出到file的log等级的开关
-
-# 第三步，定义handler的输出格式
-formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
-time_rotating_file_handler.setFormatter(formatter)
-
-# 第四步，将handler添加到logger里面
-loggers.addHandler(time_rotating_file_handler)
-
-# 如果需要同時需要在終端上輸出，定義一個streamHandler
-print_handler = logging.StreamHandler()  # 往屏幕上输出
-print_handler.setFormatter(formatter)  # 设置屏幕上显示的格式
-loggers.addHandler(print_handler)
-
-
-def Tim_Tor(func):
-    def wrapper(*args, **kwargs):
-        start_time = time.time()
-        result = func(*args, **kwargs)
-        end_time = time.time()
-        print(f"{func.__name__} took {end_time - start_time:.4f} seconds")
-        loggers.info(f"{func.__name__} took {end_time - start_time:.4f} seconds")
-        return result
-    return wrapper
-
-
-
-
-
-from itertools import product
-from multiprocessing import Pool
-def Mul_Ess(map_fun,tuple_list):
-    with Pool(31) as p:
-        P_data = p.map(map_fun, tuple_list)
-    return P_data

shancx-1.60/shancx.egg-info/SOURCES.txt

@@ -1,7 +0,0 @@
-README.md
-setup.py
-shancx/__init__.py
-shancx.egg-info/PKG-INFO
-shancx.egg-info/SOURCES.txt
-shancx.egg-info/dependency_links.txt
-shancx.egg-info/top_level.txt