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

shancx/Plot/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/Plot/radar_nmc_china_map_f.py

@@ -1,121 +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)
-
- 
-# 示例用法
-
-
-def drawCNLonLat(array_dt,Lon,Lat,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
-        ax.scatter(list(Lon), list(Lat), s=1, c='red', marker='o', transform=ccrs.PlateCarree())
-        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.scatter(list(Lon), list(Lat), s=1, c='red', marker='o', transform=ccrs.PlateCarree())
-        # 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)

shancx/Plot/radar_nmc_us_map_f.py

@@ -1,128 +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)
-
-def drawUSLonLat(array_dt,Lon,Lat,ty="CR"):   ###  x_coords2 维度 
-    # 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')
-        ax.scatter(list(Lon), list(Lat), s=1, c='red', marker='o', transform=ccrs.PlateCarree())
-        # ax.scatter(list(x_coords2), list(y_coords2), 'ro', markersize=1)
-        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=1, c='red', marker='o', transform=ccrs.PlateCarree())
-        # 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/Plot/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/Plot/subplots_single_china_map.py

@@ -1,45 +0,0 @@
-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()
-    """
-    本地路径读取 #'/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)
-
-def mkDir(filepath):
-    directory = os.path.dirname(filepath)
-    if not os.path.exists(directory):
-        os.makedirs(directory)
-
-def drawpic(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%S")
-    outpath = f"./aplot_image/{output_name}_{now_str}.png"
-    mkDir(outpath)
-    fig.tight_layout()
-    fig.savefig(outpath, dpi=300)
-    plt.close(fig)

shancx/S/__init__.py

@@ -1,138 +0,0 @@
-#!/usr/bin/python
-# -*- coding: utf-8 -*-
-import os
-def start():
-    print("import successful")
-# constants
- 
-
-__author__ = 'shancx'
- 
-__author_email__ = 'shancx@126.com'
-
-
-
-# @Time : 2023/09/27 下午8:52
-# @Author : shanchangxi
-# @File : util_log.py
-import time
-import logging  
-from logging import handlers
-
-def mkDir(path):
-    if "." in path:
-        os.makedirs(os.path.dirname(path),exist_ok=True)
-    else:
-        os.makedirs(path, exist_ok=True)
-
-loggers = logging.getLogger()
-loggers.setLevel(logging.INFO) 
-log_name =  '/home/scx/logs/project_S.log'
-mkDir(log_name)
-logfile = log_name
-time_rotating_file_handler = handlers.TimedRotatingFileHandler(filename=logfile, when='D', encoding='utf-8')
-time_rotating_file_handler.setLevel(logging.INFO)   
-formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
-time_rotating_file_handler.setFormatter(formatter)
-loggers.addHandler(time_rotating_file_handler)
-
-
-"""      
-import tracback
-try:
-
-except Exception as e:
-   
-     logger.info(traceback.format_exc())
-"""
-
-
-'''
- ##定義一個streamHandler
-# print_handler = logging.StreamHandler()  
-# print_handler.setFormatter(formatter) 
-# loggers.addHandler(print_handler)
-'''
-
-"""
-from main import makeAll,options
-from multiprocessing import Pool
-import datetime
-from config import logger,output
-import time
-import pandas as pd
-import os
-from itertools import product
-import threading
-
-def excuteCommand(cmd):
-    print(cmd)
-    os.system(cmd)
-
-def gpuPro(makeListUTC, isPhase, isDebug, gpu, isOverwrite):
-    productList = product(makeListUTC, [isPhase], [isDebug], [gpu], [isOverwrite])
-
-    with Pool(4) as p:
-        p.map(makeAll, productList)
-
-if __name__ == '__main__':
-    cfg = options()
-    isPhase = cfg.isPhase
-    isDebug = cfg.isDebug
-    sepSec = cfg.sepSec
-    gpu = cfg.gpu
-    pool = cfg.pool
-    isOverwrite = cfg.isOverwrite
-    timeList = pd.date_range(cfg.times[0], cfg.times[-1], freq=f"{sepSec}s")
-    logger.info(f"时间段check {timeList}")
-    gpuNum = 2
-    eachGPU = 4
-
-    makeListUTC = []
-    for UTC in timeList:
-        UTCStr = UTC.strftime("%Y%m%d%H%M")
-        outpath = f"{output}/{UTCStr[:4]}/{UTCStr[:8]}/MSP2_WTX_AIW_QPF_L88_CHN_{UTCStr}_00000-00300-00006.nc"
-        if not os.path.exists(outpath) or isOverwrite:
-            makeListUTC.append(UTC)
-    [print(element) for element in makeListUTC]
-
-    phaseCMD = "--isPhase" if isPhase else ""
-    debugCMD = "--isDebug" if isDebug else ""
-    OverwriteCMD = "--isOverwrite"  
-    gpuCMD = f"--gpu={gpu}"
-    # cmdList = list(map(lambda x:f"python main.py --times={x.strftime('%Y%m%d%H%M')} {phaseCMD} {debugCMD} {OverwriteCMD} {gpuCMD}",makeListUTC))
-    cmdList = list(map(lambda x:f"python main.py --times={x.strftime('%Y%m%d%H%M')} {phaseCMD} {debugCMD} {gpuCMD}",makeListUTC))
-
-    with Pool(pool) as p:
-        p.map(excuteCommand, cmdList)
-"""
-
-
-'''
-# @Time : 2023/09/27 下午8:52
-# @Author : shanchangxi
-# @File : util_log.py
-import time
-import logging  
-from logging import handlers
- 
-logger = logging.getLogger()
-logger.setLevel(logging.INFO) 
-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)   
-formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
-time_rotating_file_handler.setFormatter(formatter)
-logger.addHandler(time_rotating_file_handler)
-print_handler = logging.StreamHandler()   
-print_handler.setFormatter(formatter)   
-logger.addHandler(print_handler)
-
-'''
-
-'''
-###解决方法  pip install torch==2.4.0  torchvision    torchaudio三个同时安装  python 3.12  解决cuda启动不了的问题
-
-Res网络
-'''

shancx/W/__init__.py

@@ -1,132 +0,0 @@
-#!/usr/bin/python
-# -*- coding: utf-8 -*-
-import os
-def start():
-    print("import successful")
-# constants 
-
-__author__ = 'shancx'
- 
-__author_email__ = 'shancx@126.com'
-
-# @Time : 2023/09/27 下午8:52
-# @Author : shanchangxi
-# @File : util_log.py
-import time
-import logging  
-from logging import handlers
-
-def mkDir(path):
-    if "." in path:
-        os.makedirs(os.path.dirname(path),exist_ok=True)
-    else:
-        os.makedirs(path, exist_ok=True)
-
-loggers = logging.getLogger()
-loggers.setLevel(logging.INFO) 
-log_name =  '/home/scx/logs/project_W.log'
-mkDir(log_name)
-logfile = log_name
-time_rotating_file_handler = handlers.TimedRotatingFileHandler(filename=logfile, when='D', encoding='utf-8')
-time_rotating_file_handler.setLevel(logging.INFO)   
-formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
-time_rotating_file_handler.setFormatter(formatter)
-loggers.addHandler(time_rotating_file_handler)
-
-"""      
-import tracback
-try:
-
-except Exception as e:
-   
-     logger.info(traceback.format_exc())
-"""
-'''
- ##定義一個streamHandler
-# print_handler = logging.StreamHandler()  
-# print_handler.setFormatter(formatter) 
-# loggers.addHandler(print_handler)
-'''
-
-"""
-from main import makeAll,options
-from multiprocessing import Pool
-import datetime
-from config import logger,output
-import time
-import pandas as pd
-import os
-from itertools import product
-import threading
-
-def excuteCommand(cmd):
-    print(cmd)
-    os.system(cmd)
-
-def gpuPro(makeListUTC, isPhase, isDebug, gpu, isOverwrite):
-    productList = product(makeListUTC, [isPhase], [isDebug], [gpu], [isOverwrite])
-
-    with Pool(4) as p:
-        p.map(makeAll, productList)
-
-if __name__ == '__main__':
-    cfg = options()
-    isPhase = cfg.isPhase
-    isDebug = cfg.isDebug
-    sepSec = cfg.sepSec
-    gpu = cfg.gpu
-    pool = cfg.pool
-    isOverwrite = cfg.isOverwrite
-    timeList = pd.date_range(cfg.times[0], cfg.times[-1], freq=f"{sepSec}s")
-    logger.info(f"时间段check {timeList}")
-    gpuNum = 2
-    eachGPU = 4
-
-    makeListUTC = []
-    for UTC in timeList:
-        UTCStr = UTC.strftime("%Y%m%d%H%M")
-        outpath = f"{output}/{UTCStr[:4]}/{UTCStr[:8]}/MSP2_WTX_AIW_QPF_L88_CHN_{UTCStr}_00000-00300-00006.nc"
-        if not os.path.exists(outpath) or isOverwrite:
-            makeListUTC.append(UTC)
-    [print(element) for element in makeListUTC]
-
-    phaseCMD = "--isPhase" if isPhase else ""
-    debugCMD = "--isDebug" if isDebug else ""
-    OverwriteCMD = "--isOverwrite"  
-    gpuCMD = f"--gpu={gpu}"
-    # cmdList = list(map(lambda x:f"python main.py --times={x.strftime('%Y%m%d%H%M')} {phaseCMD} {debugCMD} {OverwriteCMD} {gpuCMD}",makeListUTC))
-    cmdList = list(map(lambda x:f"python main.py --times={x.strftime('%Y%m%d%H%M')} {phaseCMD} {debugCMD} {gpuCMD}",makeListUTC))
-
-    with Pool(pool) as p:
-        p.map(excuteCommand, cmdList)
-"""
-
-
-'''
-# @Time : 2023/09/27 下午8:52
-# @Author : shanchangxi
-# @File : util_log.py
-import time
-import logging  
-from logging import handlers
- 
-logger = logging.getLogger()
-logger.setLevel(logging.INFO) 
-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)   
-formatter = logging.Formatter("%(asctime)s - %(filename)s[line:%(lineno)d] - %(levelname)s: %(message)s")
-time_rotating_file_handler.setFormatter(formatter)
-logger.addHandler(time_rotating_file_handler)
-print_handler = logging.StreamHandler()   
-print_handler.setFormatter(formatter)   
-logger.addHandler(print_handler)
-
-'''
-
-'''
-###解决方法  pip install torch==2.4.0  torchvision    torchaudio三个同时安装  python 3.12  解决cuda启动不了的问题
-
-Res网络
-'''