shancx 1.9.33.169__py3-none-any.whl → 1.9.33.171__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- shancx/NN/__init__.py +27 -1
- {shancx-1.9.33.169.dist-info → shancx-1.9.33.171.dist-info}/METADATA +1 -1
- {shancx-1.9.33.169.dist-info → shancx-1.9.33.171.dist-info}/RECORD +5 -9
- shancx/H9/__init__.py +0 -126
- shancx/H9/ahi_read_hsd.py +0 -877
- shancx/H9/ahisearchtable.py +0 -298
- shancx/H9/geometry.py +0 -2439
- {shancx-1.9.33.169.dist-info → shancx-1.9.33.171.dist-info}/WHEEL +0 -0
- {shancx-1.9.33.169.dist-info → shancx-1.9.33.171.dist-info}/top_level.txt +0 -0
shancx/NN/__init__.py
CHANGED
|
@@ -15,7 +15,7 @@ import logging
|
|
|
15
15
|
from logging.handlers import RotatingFileHandler
|
|
16
16
|
import os
|
|
17
17
|
from shancx import crDir
|
|
18
|
-
def _loggers(logger_name="loggers", root="./logs", phase="project", level=logging.INFO, screen=True, max_bytes=10*1024*1024, backup_count=5, overwrite=False,handlersflag=
|
|
18
|
+
def _loggers(logger_name="loggers", root="./logs", phase="project", level=logging.INFO, screen=True, max_bytes=10*1024*1024, backup_count=5, overwrite=False,handlersflag=False):
|
|
19
19
|
'''set up logger with rotating file handler'''
|
|
20
20
|
l = logging.getLogger(logger_name)
|
|
21
21
|
if handlersflag:
|
|
@@ -64,6 +64,32 @@ if __name__ == "__main__":
|
|
|
64
64
|
logger.info("这是一条INFO日志")
|
|
65
65
|
'''
|
|
66
66
|
|
|
67
|
+
from itertools import product
|
|
68
|
+
from concurrent.futures import ProcessPoolExecutor as PoolExecutor, as_completed
|
|
69
|
+
import sys
|
|
70
|
+
from tqdm import tqdm
|
|
71
|
+
def validate_param_list(param_list):
|
|
72
|
+
if len(param_list) == 0:
|
|
73
|
+
raise ValueError("param_list cannot be empty.")
|
|
74
|
+
for sublist in param_list:
|
|
75
|
+
if len(sublist) == 0:
|
|
76
|
+
raise ValueError("Sub-lists in param_list cannot be empty.")
|
|
77
|
+
def Mul_sub(task, param_list, num=6):
|
|
78
|
+
print(f"Pro num {num}")
|
|
79
|
+
validate_param_list(param_list)
|
|
80
|
+
if len(param_list) == 1:
|
|
81
|
+
product_list = [(x,) for x in param_list[0]]
|
|
82
|
+
else:
|
|
83
|
+
product_list = list(product(*param_list))
|
|
84
|
+
with PoolExecutor(max_workers=num) as executor:
|
|
85
|
+
try:
|
|
86
|
+
futures = [executor.submit(task, item) for item in product_list]
|
|
87
|
+
for future in tqdm(as_completed(futures), total=len(futures), desc="Processing tasks", unit="task"):
|
|
88
|
+
future.result()
|
|
89
|
+
except KeyboardInterrupt:
|
|
90
|
+
sys.exit(1)
|
|
91
|
+
print("All tasks completed")
|
|
92
|
+
|
|
67
93
|
# import os
|
|
68
94
|
# from datetime import datetime
|
|
69
95
|
# from shancx.NN import _loggers
|
|
@@ -56,13 +56,9 @@ shancx/Fillmiss/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
|
|
|
56
56
|
shancx/Fillmiss/imgidwJU.py,sha256=CyP4ZlhPpXPqGVzzcLR7lohVl8bWjEceJLo7yvKgtEM,1708
|
|
57
57
|
shancx/Fillmiss/imgidwLatLonJU.py,sha256=ltYZj9CwN8cVs4Kow64oVN1EiBERgW9FceTUBFkCOAs,4314
|
|
58
58
|
shancx/Gpu/__init__.py,sha256=4Ahq04phTGVlFWN9Vih0WAh-IqFrhtwM5hj4G1IU2Dk,1950
|
|
59
|
-
shancx/H9/__init__.py,sha256=FCarcXfU2tVD2KrCGKNNUuL51zAWZYHV7lrP7gntvaI,4911
|
|
60
|
-
shancx/H9/ahi_read_hsd.py,sha256=bt9oOOARcXijmyGpmHYXj0NKnuTntZjqx0_tu6Vp2vs,33522
|
|
61
|
-
shancx/H9/ahisearchtable.py,sha256=e2kpz-P5npgL4gzNxn8igERJuWWIysvTNLkptr5_Zcc,9579
|
|
62
|
-
shancx/H9/geometry.py,sha256=ZoCfgP07ANqDxyLf8uFN016DPl7QqW9pAV3sr2ekCng,95893
|
|
63
59
|
shancx/Hug/__init__.py,sha256=NvcqXM0CSwXYtTr5yj4zfkag_7PZcI5tCsIHjuLNXTQ,3040
|
|
64
60
|
shancx/Mos/__init__.py,sha256=IDlCLNN6a_mguYmRcAPcevHSDaud5ZhvYNSuASqztcE,1563
|
|
65
|
-
shancx/NN/__init__.py,sha256=
|
|
61
|
+
shancx/NN/__init__.py,sha256=fVSAQnfE_34zBNw5gz_K5jenY5iE30yoB0eEleMCYqI,4235
|
|
66
62
|
shancx/Plot/GlobMap.py,sha256=WRvZOf5bo8IfECk1-sN4g6K0CrcJhSoEJAiLvZAWDwg,13021
|
|
67
63
|
shancx/Plot/Gray2RGB.py,sha256=yDXu1xWVi7OtrS8_ExxK_ev1988kj1ZxoNCWf-kl4Zc,2734
|
|
68
64
|
shancx/Plot/__init__.py,sha256=R4pywzeVDqEWGMkWeGBAVmopgPBThiLv2RBDlGquNGc,31891
|
|
@@ -83,7 +79,7 @@ shancx/Train/multiGpu.py,sha256=D_oZeiSc7VWktpnVDwrFOC1CYZSt9rxOKY5lngE5vFg,820
|
|
|
83
79
|
shancx/Train/prepare.py,sha256=vL_8UOA66oZCBIwCICtihsGibivtNgaVJGulJxfNdn8,6793
|
|
84
80
|
shancx/Train/renet50.py,sha256=wEhYk1X96WE5zuqHqVxWLJa-A5jDNkz4z6edORNufnA,6428
|
|
85
81
|
shancx/tensBoard/__init__.py,sha256=ga2C5YyJITvvQA1ocpxna_KNFnNRJVwkTjLoIglLZUQ,993
|
|
86
|
-
shancx-1.9.33.
|
|
87
|
-
shancx-1.9.33.
|
|
88
|
-
shancx-1.9.33.
|
|
89
|
-
shancx-1.9.33.
|
|
82
|
+
shancx-1.9.33.171.dist-info/METADATA,sha256=_l_QlL7aOUxagFQpVMYER4Misz3TxtUv_T9Ah8QRXTo,644
|
|
83
|
+
shancx-1.9.33.171.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
|
|
84
|
+
shancx-1.9.33.171.dist-info/top_level.txt,sha256=akfCS1vKWz3pNmEN_yN9ZiGp-60IQY5ET38mRx_i_-4,7
|
|
85
|
+
shancx-1.9.33.171.dist-info/RECORD,,
|
shancx/H9/__init__.py
DELETED
|
@@ -1,126 +0,0 @@
|
|
|
1
|
-
from shancx.H9.ahisearchtable import ahisearchtable
|
|
2
|
-
from shancx.H9.ahi_read_hsd import ahi_read_hsd
|
|
3
|
-
from tqdm import tqdm
|
|
4
|
-
import os
|
|
5
|
-
import numpy as np
|
|
6
|
-
import datetime
|
|
7
|
-
import time
|
|
8
|
-
import warnings
|
|
9
|
-
warnings.filterwarnings("ignore", category=UserWarning)
|
|
10
|
-
class AHIScene(ahi_read_hsd, ahisearchtable) :
|
|
11
|
-
def __init__(self, subpoint=140.7, resolution=0.02):
|
|
12
|
-
super().__init__(subpoint=subpoint, resolution=resolution)
|
|
13
|
-
self.Tempfile = []
|
|
14
|
-
def hsdBlock(self, srcHSDfiles, tmppath, fillvalue=65535) :
|
|
15
|
-
''' 对H8、H9的HSD文件进行解析、拼接成NOM '''
|
|
16
|
-
# HS_H09_20230115_0400_B01_FLDK_R10_S0110.DAT.bz2
|
|
17
|
-
BandID, BlockIDMin, BlockIDMax, SegmentTotal = self.setHSDInfo(srcHSDfiles)
|
|
18
|
-
outdata = None
|
|
19
|
-
BlockIDs = []
|
|
20
|
-
with tqdm(total=len(srcHSDfiles), iterable='iterable',
|
|
21
|
-
desc = '正在进行第%i波段块合成' %(BandID), mininterval=1) as pbar:
|
|
22
|
-
for hsdname in srcHSDfiles :
|
|
23
|
-
if not os.path.isfile(hsdname):
|
|
24
|
-
print('文件不存在【%s】' %(hsdname))
|
|
25
|
-
pbar.update(1)
|
|
26
|
-
continue
|
|
27
|
-
|
|
28
|
-
# 获取文件名信息
|
|
29
|
-
nameinfo = self.getHSDNameInfo(hsdname)
|
|
30
|
-
if nameinfo is None :
|
|
31
|
-
pbar.update(1)
|
|
32
|
-
continue
|
|
33
|
-
SegmentNum = nameinfo['SegmemtID']
|
|
34
|
-
|
|
35
|
-
# print('正在解压bz2文件【%s】' %(hsdname))
|
|
36
|
-
self._unzipped = self.unzip_file(hsdname, tmppath)
|
|
37
|
-
if self._unzipped:
|
|
38
|
-
self.is_zipped = True
|
|
39
|
-
filename = self._unzipped
|
|
40
|
-
|
|
41
|
-
self.Tempfile.append(filename)
|
|
42
|
-
else:
|
|
43
|
-
filename = hsdname
|
|
44
|
-
|
|
45
|
-
if filename.endswith('.bz2') :
|
|
46
|
-
print('解压bz2文件失败【%s】' %(filename))
|
|
47
|
-
pbar.update(1)
|
|
48
|
-
continue
|
|
49
|
-
|
|
50
|
-
# 根据块号对数据进行拼接
|
|
51
|
-
data = self.readhsd(filename, SegmentNum)
|
|
52
|
-
if data is None :
|
|
53
|
-
pbar.update(1)
|
|
54
|
-
continue
|
|
55
|
-
|
|
56
|
-
if outdata is None :
|
|
57
|
-
line, pixel = data.shape
|
|
58
|
-
outdata = np.full(shape=(line*SegmentTotal, pixel),
|
|
59
|
-
fill_value=fillvalue, dtype=np.uint16)
|
|
60
|
-
|
|
61
|
-
data[np.isnan(data)] = fillvalue/100.0
|
|
62
|
-
outdata[(SegmentNum-BlockIDMin)*line:(SegmentNum-BlockIDMin+1)*line, :] \
|
|
63
|
-
= np.array(data*100.0, dtype=np.uint16)
|
|
64
|
-
BlockIDs.append(SegmentNum)
|
|
65
|
-
pbar.update(1)
|
|
66
|
-
pbar.close()
|
|
67
|
-
self.__del__()
|
|
68
|
-
return outdata
|
|
69
|
-
|
|
70
|
-
def setHSDInfo(self, filelist):
|
|
71
|
-
|
|
72
|
-
BandID = None
|
|
73
|
-
BlockIDs = []
|
|
74
|
-
for filename in filelist :
|
|
75
|
-
nameinfo = self.getHSDNameInfo(filename)
|
|
76
|
-
if nameinfo is None :
|
|
77
|
-
continue
|
|
78
|
-
|
|
79
|
-
if BandID is None :
|
|
80
|
-
BandID = nameinfo['BandID']
|
|
81
|
-
elif BandID != nameinfo['BandID'] :
|
|
82
|
-
raise Exception('输入的文件列表中有多个波段的块数据文件【%s】' %(filename))
|
|
83
|
-
BlockIDs.append(nameinfo['SegmemtID'])
|
|
84
|
-
|
|
85
|
-
BlockIDMin = np.nanmin(BlockIDs)
|
|
86
|
-
BlockIDMax = np.nanmax(BlockIDs)
|
|
87
|
-
|
|
88
|
-
SegmentTotal = int(BlockIDMax-BlockIDMin+1)
|
|
89
|
-
|
|
90
|
-
return BandID, BlockIDMin, BlockIDMax, SegmentTotal
|
|
91
|
-
|
|
92
|
-
def getHSDNameInfo(self, filename):
|
|
93
|
-
|
|
94
|
-
basename = os.path.basename(filename)
|
|
95
|
-
basename = basename.split('.')[0]
|
|
96
|
-
if len(basename) != 39 :
|
|
97
|
-
print('非标准文件名,需要输入文件名【HS_H09_YYYYMMDD_HHMM_BXX_FLDK_R20_S0810】')
|
|
98
|
-
return None
|
|
99
|
-
|
|
100
|
-
nameinfo = {}
|
|
101
|
-
namelist = basename.split('_')
|
|
102
|
-
|
|
103
|
-
nameinfo['SatID'] = namelist[1]
|
|
104
|
-
nameinfo['StartTime'] = datetime.datetime.strptime('%s %s' %(namelist[2], namelist[3]), '%Y%m%d %H%M')
|
|
105
|
-
nameinfo['BandID'] = int(namelist[4][1:]) # 2-digit band number (varies from "01" to "16");
|
|
106
|
-
nameinfo['ObsType'] = namelist[5]
|
|
107
|
-
nameinfo['Resolution'] = float(namelist[6][1:])/10.0/100 # spatial resolution ("05": 0.5km, "10": 1.0km, "20": 2.0km);
|
|
108
|
-
nameinfo['SegmemtID'] = int(namelist[7][1:3])
|
|
109
|
-
nameinfo['SegmemtTotal'] = int(namelist[7][3:5]) # total number of segments (fixed to "10")
|
|
110
|
-
|
|
111
|
-
return nameinfo
|
|
112
|
-
|
|
113
|
-
def __del__(self):
|
|
114
|
-
# pass
|
|
115
|
-
for filename in self.Tempfile :
|
|
116
|
-
if os.path.isfile(filename) :
|
|
117
|
-
try:
|
|
118
|
-
os.remove(filename)
|
|
119
|
-
except BaseException as e :
|
|
120
|
-
time.sleep(1)
|
|
121
|
-
try:
|
|
122
|
-
fp = open(filename, 'r')
|
|
123
|
-
fp.close()
|
|
124
|
-
os.remove(filename)
|
|
125
|
-
except BaseException as e :
|
|
126
|
-
pass
|