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

shancx/E/__init__.py

@@ -1,152 +1,178 @@
-#!/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_E.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网络
-'''
-
-'''
-#!/bin/bash
-subject="Daily System Report"
-to_email="shanhe12@163.com"
-temp_file="/home/scx/logs/111.log"  # 替换为日志文件的实际路径
-search_terms=("error" "ERROR")  # 可以添加多个搜索词
-for term in "${search_terms[@]}"; do
-    if grep -q "$term" "$temp_file"; then
-        grep "$term" "$temp_file" | mail -s "$subject - $term" "$to_email"
-    else
-        echo "No matches found for '$term'." >&2 
-    fi
-done
-
-'''
+#!/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 =  './project_E.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 logging
+def setup_logger(logger_name=loggers, root="./", phase="project", level=logging.INFO, screen=False):
+    '''set up logger'''
+    l = logging.getLogger(logger_name)
+    formatter = logging.Formatter(
+        '%(asctime)s.%(msecs)03d - %(levelname)s: %(message)s', datefmt='%y-%m-%d %H:%M:%S')
+    log_file = os.path.join(root, '{}.log'.format(phase))
+    fh = logging.FileHandler(log_file, mode='w')
+    fh.setFormatter(formatter)
+    l.setLevel(level)
+    l.addHandler(fh)
+    if screen:
+        sh = logging.StreamHandler()
+        sh.setFormatter(formatter)
+        l.addHandler(sh)
+    return l
+
+def greetJU(**kwargs):
+    if 'name' in kwargs:
+        print(f"Hello, {kwargs['name']}!")
+    else:
+        print("Hello, World!")
+# greet(name='Alice')  # 输出：Hello, Alice!
+# greet()  # 
+
+"""      
+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网络
+'''
+
+'''
+#!/bin/bash
+subject="Daily System Report"
+to_email="shanhe12@163.com"
+temp_file="/home/scx/logs/111.log"  # 替换为日志文件的实际路径
+search_terms=("error" "ERROR")  # 可以添加多个搜索词
+for term in "${search_terms[@]}"; do
+    if grep -q "$term" "$temp_file"; then
+        grep "$term" "$temp_file" | mail -s "$subject - $term" "$to_email"
+    else
+        echo "No matches found for '$term'." >&2 
+    fi
+done
+
+'''
+
+"""
+manager = UserManager("./key/user_data.json")
+manager.add_user("001", '[f"{key}","https://qyapi.weixin.qq.com/cgi-bin/webhook/send?key="]')    
+user_info = manager.get_user("user002")
+"""

shancx/Fillmiss/imgidwJU.py

@@ -0,0 +1,46 @@
+import numpy as np
+from scipy.spatial import cKDTree
+import matplotlib.pyplot as plt
+def fill_missing_idw(image, power=2, max_neighbors=8):
+    filled_image = image.copy()
+    known_mask = ~np.isnan(image)
+    missing_mask = np.isnan(image)
+    known_points = np.column_stack(np.where(known_mask))
+    missing_points = np.column_stack(np.where(missing_mask))
+    known_values = image[known_mask]
+    if len(missing_points) == 0:
+        print("There are no missing points to fill.")
+        return filled_image
+    tree = cKDTree(known_points)
+    distances, indexes = tree.query(missing_points, k=max_neighbors)
+    with np.errstate(divide='ignore'): 
+        weights = 1 / distances**power
+    weights[~np.isfinite(weights)] = 0  
+    weight_sums = np.sum(weights, axis=1)
+    weight_sums[weight_sums == 0] = 1  
+    interpolated_values = np.sum(weights * known_values[indexes], axis=1) / weight_sums
+    for i, (x, y) in enumerate(missing_points):
+        filled_image[x, y] = interpolated_values[i]
+    return filled_image
+def test_fill_missing_idw():
+    np.random.seed(42)
+    image = np.random.rand(20, 20)
+    image[5:10, 5:10] = np.nan  
+    plt.figure(figsize=(12, 6))
+    plt.subplot(1, 2, 1)
+    plt.imshow(image, cmap='viridis', interpolation='none')
+    plt.title("before")
+    plt.colorbar()
+    plt.savefig("before.png")
+    filled_image = fill_missing_idw(image, power=2, max_neighbors=8)
+    plt.subplot(1, 2, 2)
+    plt.imshow(filled_image, cmap='viridis', interpolation='none')
+    plt.title("uper")
+    plt.colorbar()
+
+    plt.tight_layout()
+    plt.savefig("uper.png")
+    plt.show()
+
+if __name__ == "__main__":
+    test_fill_missing_idw()

shancx/Fillmiss/imgidwLatLonJU.py

@@ -0,0 +1,82 @@
+import numpy as np
+from netCDF4 import Dataset
+from scipy.spatial import cKDTree
+import os
+
+def fill_missing_idw(variable, lat, lon, power=2, max_neighbors=8):
+    filled_variable = variable.copy()
+    ny, nx = variable.shape
+    if isinstance(variable, np.ma.MaskedArray):
+        known_mask = ~variable.mask
+        missing_mask = variable.mask
+        known_values = variable.data[known_mask]
+    else:
+        known_mask = ~np.isnan(variable)
+        missing_mask = np.isnan(variable)
+        known_values = variable[known_mask]
+    known_points = np.column_stack((lon[known_mask], lat[known_mask]))
+    missing_points = np.column_stack((lon[missing_mask], lat[missing_mask]))
+    if len(known_points) == 0:
+        print("There are no known points available for interpolation.")
+        return filled_variable
+    if len(missing_points) == 0:
+        print("There are no missing points to fill.")
+        return filled_variable
+    tree = cKDTree(known_points)
+    distances, indexes = tree.query(missing_points, k=max_neighbors)
+    if distances.size == 0 or indexes.size == 0:
+        print("There are no valid neighbor points for interpolation.")
+        return filled_variable
+    if max_neighbors == 1:
+        distances = distances[:, np.newaxis]
+        indexes = indexes[:, np.newaxis]
+    with np.errstate(divide='ignore'): 
+        weights = 1 / distances**power
+    weights[~np.isfinite(weights)] = 0  
+    weight_sums = np.sum(weights, axis=1)
+    weight_sums[weight_sums == 0] = 1  
+    interpolated_values = np.sum(weights * known_values[indexes], axis=1) / weight_sums
+    exact_matches = distances[:, 0] == 0
+    if np.any(exact_matches):
+        interpolated_values[exact_matches] = known_values[indexes[exact_matches, 0]]
+    filled_variable[missing_mask] = interpolated_values
+    n_filled = np.sum(~missing_mask)
+    print(f"尝试填补了 {len(missing_points)} 个缺失点，成功填补了 {n_filled} 个。")    
+    return filled_variable
+def process_nc_file(input_path, output_path):
+    with Dataset(input_path, 'r') as src:
+        with Dataset(output_path, 'w', format=src.file_format) as dst:
+            dst.setncatts({attr: src.getncattr(attr) for attr in src.ncattrs()})
+            for dim_name, dim in src.dimensions.items():
+                dst.createDimension(dim_name, len(dim) if not dim.isunlimited() else None)
+            for var_name, var in src.variables.items():
+                fill_value = getattr(var, '_FillValue', None)
+                dst_var = dst.createVariable(var_name, var.datatype, var.dimensions, fill_value=fill_value)
+                dst_var.setncatts({attr: var.getncattr(attr) for attr in var.ncattrs()})
+                data = var[:]
+                if var_name in ['CR']:
+                    if data.ndim == 3: 
+                        filled_data = np.ma.array(data) 
+                        lat = src.variables['lat'][:]
+                        lon = src.variables['lon'][:]
+                        for t in range(data.shape[0]):
+                            print(f"Processing time index {t} for variable {var_name}")
+                            var_data = data[t, :, :]
+                            if not isinstance(var_data, np.ma.MaskedArray):
+                                var_data = np.ma.masked_where(np.isnan(var_data), var_data)
+                            filled_var = fill_missing_idw(var_data, lat, lon)
+                            filled_data[t, :, :] = filled_var
+                        fill_val = getattr(var, '_FillValue', np.nan)
+                        dst_var[:] = filled_data.filled(fill_val)
+                    else:
+                        raise ValueError(f"Variable {var_name} 维度不符合预期，预期为 3D 但实际为 {data.ndim}D。")
+                else:
+                    dst_var[:] = data
+    print(f"The data after filling has been saved to {output_path}")
+if __name__ == "__main__":
+    input_nc = "/mnt/wtx_weather_forecast/scx/WTX_DATA/RADA/MQPF_1204_diffu12/2024/20240908/MSP2_WTX_AIW_REF_L88_CHN_202409080448_00000-00300-00006.nc"
+    output_nc = "1aaa_filled.nc"
+    if not os.path.exists(input_nc):
+        print(f"输入文件 {input_nc} 不存在。")
+    else:
+        process_nc_file(input_nc, output_nc)

shancx/Gpu/__init__.py

@@ -0,0 +1,55 @@
+#!/usr/bin/python
+# -*- coding: utf-8 -*-
+# @Time : 2024/10/17 上午午10:40
+# @Author : shancx
+# @File : __init__.py
+# @email : shanhe12@163.com
+import torch
+def CheckGpuPlus(num=1):
+    if torch.cuda.is_available():
+        print(f"CUDA is available. Number of GPUs available: {torch.cuda.device_count()}")
+        
+        for i in range(torch.cuda.device_count()):
+            print(f"GPU {i}: {torch.cuda.get_device_name(i)}")
+        device = torch.device(f'cuda:{num}' if torch.cuda.is_available() else 'cpu')   
+        return device
+    else:
+        print("CUDA is not available. Using CPU.")
+        return None
+
+#pd.concat(filter(None, results))
+#valid_results = [df for df in results if isinstance(df, pd.DataFrame) and not df.empty]
+
+import os
+def visDevices(device_ids):
+    if isinstance(device_ids, int):
+        device_ids = str(device_ids)
+    elif isinstance(device_ids, list):
+        device_ids = ",".join(map(str, device_ids))
+    os.environ["CUDA_VISIBLE_DEVICES"] = device_ids
+    if torch.cuda.is_available():
+        print(f"Visible GPUs: {device_ids}")
+        print(f"Current visible GPUs: {os.environ.get('CUDA_VISIBLE_DEVICES')}")
+        for device_id in device_ids.split(","):
+            if not torch.cuda.device(int(device_id)):
+                print(f"Warning: GPU {device_id} is not available.")
+    else:
+        print("No GPU available. Using CPU.")
+
+import torch
+import torch.nn as nn
+
+def multiGpu(model, gpu_ids):
+    # 检查是否有可用的 GPU
+    if not torch.cuda.is_available():
+        print("CUDA is not available. Using CPU.")
+        device = torch.device("cpu")
+        return model.to(device), device
+    device = torch.device(f"cuda:{gpu_ids[0]}")
+    if len(gpu_ids) > 1:
+        print(f"Using {len(gpu_ids)} GPUs: {gpu_ids}")
+        model = nn.DataParallel(model, device_ids=gpu_ids)
+    else:
+        print(f"Using GPU: {gpu_ids[0]}")
+    model = model.to(device)
+    return model, device

shancx/H9/__init__.py

@@ -0,0 +1,126 @@
+from shancx.H9.ahisearchtable import ahisearchtable
+from shancx.H9.ahi_read_hsd import ahi_read_hsd
+from tqdm import tqdm
+import os
+import numpy as np 
+import datetime 
+import time
+import warnings
+warnings.filterwarnings("ignore", category=UserWarning)
+class AHIScene(ahi_read_hsd, ahisearchtable) :
+    def __init__(self, subpoint=140.7, resolution=0.02):
+        super().__init__(subpoint=subpoint, resolution=resolution)
+        self.Tempfile = []
+    def hsdBlock(self, srcHSDfiles, tmppath, fillvalue=65535) :
+        ''' 对H8、H9的HSD文件进行解析、拼接成NOM '''
+        # HS_H09_20230115_0400_B01_FLDK_R10_S0110.DAT.bz2
+        BandID, BlockIDMin, BlockIDMax, SegmentTotal = self.setHSDInfo(srcHSDfiles)
+        outdata = None
+        BlockIDs = []
+        with tqdm(total=len(srcHSDfiles), iterable='iterable',
+                  desc = '正在进行第%i波段块合成' %(BandID), mininterval=1) as pbar:
+            for hsdname in srcHSDfiles :
+                if not os.path.isfile(hsdname):
+                    print('文件不存在【%s】' %(hsdname))
+                    pbar.update(1)
+                    continue
+
+                # 获取文件名信息
+                nameinfo = self.getHSDNameInfo(hsdname)
+                if nameinfo is None :
+                    pbar.update(1)
+                    continue
+                SegmentNum = nameinfo['SegmemtID']
+
+                # print('正在解压bz2文件【%s】' %(hsdname))
+                self._unzipped = self.unzip_file(hsdname, tmppath)
+                if self._unzipped:
+                    self.is_zipped = True
+                    filename = self._unzipped
+
+                    self.Tempfile.append(filename)
+                else:
+                    filename = hsdname
+
+                if filename.endswith('.bz2') :
+                    print('解压bz2文件失败【%s】' %(filename))
+                    pbar.update(1)
+                    continue
+
+                # 根据块号对数据进行拼接
+                data = self.readhsd(filename, SegmentNum)
+                if data is None :
+                    pbar.update(1)
+                    continue
+
+                if outdata is None :
+                    line, pixel = data.shape
+                    outdata = np.full(shape=(line*SegmentTotal, pixel),
+                                      fill_value=fillvalue, dtype=np.uint16)
+
+                data[np.isnan(data)] = fillvalue/100.0
+                outdata[(SegmentNum-BlockIDMin)*line:(SegmentNum-BlockIDMin+1)*line, :] \
+                    = np.array(data*100.0, dtype=np.uint16)
+                BlockIDs.append(SegmentNum)
+                pbar.update(1)
+        pbar.close()
+        self.__del__()
+        return outdata 
+  
+    def setHSDInfo(self, filelist):
+
+        BandID = None
+        BlockIDs = []
+        for filename in filelist :
+            nameinfo = self.getHSDNameInfo(filename)
+            if nameinfo is None :
+                continue
+
+            if BandID is None :
+                BandID = nameinfo['BandID']
+            elif BandID != nameinfo['BandID'] :
+                raise Exception('输入的文件列表中有多个波段的块数据文件【%s】' %(filename))
+            BlockIDs.append(nameinfo['SegmemtID'])
+
+        BlockIDMin = np.nanmin(BlockIDs)
+        BlockIDMax = np.nanmax(BlockIDs)
+
+        SegmentTotal = int(BlockIDMax-BlockIDMin+1)
+
+        return BandID, BlockIDMin, BlockIDMax, SegmentTotal
+
+    def getHSDNameInfo(self, filename):
+
+        basename = os.path.basename(filename)
+        basename = basename.split('.')[0]
+        if len(basename) != 39 :
+            print('非标准文件名，需要输入文件名【HS_H09_YYYYMMDD_HHMM_BXX_FLDK_R20_S0810】')
+            return None
+
+        nameinfo = {}
+        namelist = basename.split('_')
+
+        nameinfo['SatID'] = namelist[1]
+        nameinfo['StartTime'] = datetime.datetime.strptime('%s %s' %(namelist[2], namelist[3]), '%Y%m%d %H%M')
+        nameinfo['BandID'] = int(namelist[4][1:])   # 2-digit band number (varies from "01" to "16");
+        nameinfo['ObsType'] = namelist[5]
+        nameinfo['Resolution'] = float(namelist[6][1:])/10.0/100    # spatial resolution ("05": 0.5km, "10": 1.0km, "20": 2.0km);
+        nameinfo['SegmemtID'] = int(namelist[7][1:3])
+        nameinfo['SegmemtTotal'] = int(namelist[7][3:5])    # total number of segments (fixed to "10")
+
+        return nameinfo
+
+    def __del__(self):
+        # pass
+        for filename in self.Tempfile :
+            if os.path.isfile(filename) :
+                try:
+                    os.remove(filename)
+                except BaseException as e :
+                    time.sleep(1)
+                    try:
+                        fp = open(filename, 'r')
+                        fp.close()
+                        os.remove(filename)
+                    except BaseException as e :
+                        pass