oafuncs 0.0.98.11__py3-none-any.whl → 0.0.98.13__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.
- oafuncs/_script/netcdf_merge.py +7 -35
- oafuncs/oa_tool.py +6 -1
- {oafuncs-0.0.98.11.dist-info → oafuncs-0.0.98.13.dist-info}/METADATA +1 -1
- {oafuncs-0.0.98.11.dist-info → oafuncs-0.0.98.13.dist-info}/RECORD +7 -7
- {oafuncs-0.0.98.11.dist-info → oafuncs-0.0.98.13.dist-info}/WHEEL +0 -0
- {oafuncs-0.0.98.11.dist-info → oafuncs-0.0.98.13.dist-info}/licenses/LICENSE.txt +0 -0
- {oafuncs-0.0.98.11.dist-info → oafuncs-0.0.98.13.dist-info}/top_level.txt +0 -0
oafuncs/_script/netcdf_merge.py
CHANGED
@@ -3,7 +3,7 @@ from typing import List, Optional, Union
|
|
3
3
|
|
4
4
|
import numpy as np
|
5
5
|
import xarray as xr
|
6
|
-
|
6
|
+
|
7
7
|
from oafuncs import pbar
|
8
8
|
|
9
9
|
|
@@ -25,6 +25,7 @@ def merge_nc(file_list: Union[str, List[str]], var_name: Optional[Union[str, Lis
|
|
25
25
|
merge(file_list, var_name=['u', 'v'], dim_name='time', target_filename='merged.nc')
|
26
26
|
merge(file_list, var_name=None, dim_name='time', target_filename='merged.nc')
|
27
27
|
"""
|
28
|
+
from oafuncs._script.netcdf_write import save_to_nc
|
28
29
|
|
29
30
|
if target_filename is None:
|
30
31
|
target_filename = "merged.nc"
|
@@ -51,7 +52,7 @@ def merge_nc(file_list: Union[str, List[str]], var_name: Optional[Union[str, Lis
|
|
51
52
|
# 初始化合并数据字典
|
52
53
|
merged_data = {}
|
53
54
|
|
54
|
-
for i, file in pbar(enumerate(file_list), description="Reading files",
|
55
|
+
for i, file in pbar(enumerate(file_list), description="Reading files", total=len(file_list)):
|
55
56
|
with xr.open_dataset(file) as ds:
|
56
57
|
for var in var_names:
|
57
58
|
data_var = ds[var]
|
@@ -64,53 +65,24 @@ def merge_nc(file_list: Union[str, List[str]], var_name: Optional[Union[str, Lis
|
|
64
65
|
else:
|
65
66
|
merged_data[var] = data_var.fillna(0)
|
66
67
|
|
67
|
-
for var in pbar(merged_data, description="Merging variables"
|
68
|
+
for var in pbar(merged_data, description="Merging variables"):
|
68
69
|
if isinstance(merged_data[var], list):
|
69
70
|
# 判断类型,时间类型用NaT填充
|
70
71
|
if np.issubdtype(merged_data[var][0].dtype, np.datetime64):
|
71
72
|
merged_data[var] = xr.concat(merged_data[var], dim=dim_name).fillna(np.datetime64("NaT"))
|
72
73
|
else:
|
73
74
|
merged_data[var] = xr.concat(merged_data[var], dim=dim_name).fillna(0)
|
74
|
-
# print(f"Variable '{var}' merged: min={merged_data[var].min().values:.3f}, max={merged_data[var].max().values:.3f}, mean={merged_data[var].mean().values:.3f}")
|
75
75
|
|
76
|
-
# 修改写入数据部分,支持压缩并设置基数和比例因子
|
77
|
-
# print("\nWriting data to file ...")
|
78
76
|
if os.path.exists(target_filename):
|
79
77
|
print("Warning: The target file already exists. Removing it ...")
|
80
78
|
os.remove(target_filename)
|
81
|
-
|
82
|
-
|
83
|
-
encoding = {}
|
84
|
-
for var in merged_dataset.data_vars:
|
85
|
-
data = merged_dataset[var].values
|
86
|
-
# print(f"Variable '{var}' ready for writing: min={data.min():.3f}, max={data.max():.3f}, mean={data.mean():.3f}")
|
87
|
-
if data.dtype.kind in {"i", "u", "f"}: # 仅对数值型数据进行压缩
|
88
|
-
data_range = data.max() - data.min()
|
89
|
-
if data_range > 0: # 避免范围过小导致的精度问题
|
90
|
-
scale_factor = data_range / (2**16 - 1)
|
91
|
-
add_offset = data.min()
|
92
|
-
encoding[var] = {
|
93
|
-
"zlib": True,
|
94
|
-
"complevel": 4,
|
95
|
-
"dtype": "int16",
|
96
|
-
"scale_factor": scale_factor,
|
97
|
-
"add_offset": add_offset,
|
98
|
-
"_FillValue": -32767,
|
99
|
-
}
|
100
|
-
else:
|
101
|
-
encoding[var] = {"zlib": True, "complevel": 4} # 范围过小时禁用缩放
|
102
|
-
else:
|
103
|
-
encoding[var] = {"zlib": True, "complevel": 4} # 非数值型数据不使用缩放
|
104
|
-
|
105
|
-
# 确保写入时不会因编码问题导致数据丢失
|
106
|
-
# merged_dataset.to_netcdf(target_filename, encoding=encoding)
|
107
|
-
delayed_write = merged_dataset.to_netcdf(target_filename, encoding=encoding, compute=False)
|
108
|
-
with ProgressBar():
|
109
|
-
delayed_write.compute()
|
79
|
+
|
80
|
+
save_to_nc(target_filename, merged_data)
|
110
81
|
|
111
82
|
print(f'\nFile "{target_filename}" has been successfully created.')
|
112
83
|
|
113
84
|
|
85
|
+
|
114
86
|
# Example usage
|
115
87
|
if __name__ == "__main__":
|
116
88
|
files_to_merge = ["file1.nc", "file2.nc", "file3.nc"]
|
oafuncs/oa_tool.py
CHANGED
@@ -2,6 +2,7 @@ import datetime
|
|
2
2
|
import logging
|
3
3
|
from typing import Any, Callable, Dict, Iterable, List, Optional, Tuple, Union
|
4
4
|
|
5
|
+
|
5
6
|
from rich import print
|
6
7
|
|
7
8
|
from ._script.parallel import ParallelExecutor
|
@@ -174,6 +175,8 @@ def pbar(
|
|
174
175
|
"""
|
175
176
|
from ._script.cprogressbar import ColorProgressBar
|
176
177
|
import random
|
178
|
+
|
179
|
+
# number = random.randint(1, 999)
|
177
180
|
|
178
181
|
def _generate_random_color_hex():
|
179
182
|
"""Generate a random color in hexadecimal format."""
|
@@ -186,9 +189,11 @@ def pbar(
|
|
186
189
|
color = _generate_random_color_hex()
|
187
190
|
|
188
191
|
style = f"bold {color if color != 'None' else 'green'}"
|
192
|
+
# print(f"[{style}]~*^* {description} *^*~ -> {number:03d}[/{style}]")
|
189
193
|
print(f"[{style}]~*^* {description} *^*~[/{style}]")
|
190
194
|
|
191
|
-
description=''
|
195
|
+
# description=f'{number:03d}'
|
196
|
+
description = ""
|
192
197
|
|
193
198
|
return ColorProgressBar(
|
194
199
|
iterable=iterable,
|
@@ -7,12 +7,12 @@ oafuncs/oa_file.py,sha256=j9gXJgPOJsliu4IOUc4bc-luW4yBvQyNCEmMyDVjUwQ,16404
|
|
7
7
|
oafuncs/oa_help.py,sha256=_4AZgRDq5Or0vauNvq5IDDHIBoBfdOQtzak-mG1wwAw,4537
|
8
8
|
oafuncs/oa_nc.py,sha256=lLEPjj4qgdEw1al0r1nKGZUnMP_ejT8A2NKD4lrs2kc,11936
|
9
9
|
oafuncs/oa_python.py,sha256=NkopwkYFGSEuVljnTBvXCl6o2CeyRNBqRXSsUl3euEE,5192
|
10
|
-
oafuncs/oa_tool.py,sha256=
|
10
|
+
oafuncs/oa_tool.py,sha256=rpPkLqWhqMmqlCc5wjL8qMTg3gThCkSrYJckbX_0iJc,8631
|
11
11
|
oafuncs/_data/hycom.png,sha256=MadKs6Gyj5n9-TOu7L4atQfTXtF9dvN9w-tdU9IfygI,10945710
|
12
12
|
oafuncs/_data/oafuncs.png,sha256=o3VD7wm-kwDea5E98JqxXl04_78cBX7VcdUt7uQXGiU,3679898
|
13
13
|
oafuncs/_script/cprogressbar.py,sha256=UIgGcLFs-6IgWlITuBLaQqrpt4OAK3Mst5RlCiNfZdQ,15772
|
14
14
|
oafuncs/_script/email.py,sha256=lL4HGKrr524-g0xLlgs-4u7x4-u7DtgNoD9AL8XJKj4,3058
|
15
|
-
oafuncs/_script/netcdf_merge.py,sha256=
|
15
|
+
oafuncs/_script/netcdf_merge.py,sha256=vWvrFUmzwZgTCn2E_9iBKmPLf93ekvD3Ud2NOQgJxGQ,3829
|
16
16
|
oafuncs/_script/netcdf_modify.py,sha256=sGRUYNhfGgf9JV70rnBzw3bzuTRSXzBTL_RMDnDPeLQ,4552
|
17
17
|
oafuncs/_script/netcdf_write.py,sha256=iO1Qv9bp6RLiw1D8Nrv7tX_8X-diUZaX3Nxhk6pJ5Nw,8556
|
18
18
|
oafuncs/_script/parallel.py,sha256=T9Aie-e4LcbKlFTLZ0l4lhEN3SBVa84jRcrAsIm8s0I,8767
|
@@ -37,8 +37,8 @@ oafuncs/oa_sign/__init__.py,sha256=QKqTFrJDFK40C5uvk48GlRRbGFzO40rgkYwu6dYxatM,5
|
|
37
37
|
oafuncs/oa_sign/meteorological.py,sha256=8091SHo2L8kl4dCFmmSH5NGVHDku5i5lSiLEG5DLnOQ,6489
|
38
38
|
oafuncs/oa_sign/ocean.py,sha256=xrW-rWD7xBWsB5PuCyEwQ1Q_RDKq2KCLz-LOONHgldU,5932
|
39
39
|
oafuncs/oa_sign/scientific.py,sha256=a4JxOBgm9vzNZKpJ_GQIQf7cokkraV5nh23HGbmTYKw,5064
|
40
|
-
oafuncs-0.0.98.
|
41
|
-
oafuncs-0.0.98.
|
42
|
-
oafuncs-0.0.98.
|
43
|
-
oafuncs-0.0.98.
|
44
|
-
oafuncs-0.0.98.
|
40
|
+
oafuncs-0.0.98.13.dist-info/licenses/LICENSE.txt,sha256=rMtLpVg8sKiSlwClfR9w_Dd_5WubTQgoOzE2PDFxzs4,1074
|
41
|
+
oafuncs-0.0.98.13.dist-info/METADATA,sha256=S7BtUR9AaROiwbgyU3-wFsbVo8ySvRNUR3M9x383uA8,4273
|
42
|
+
oafuncs-0.0.98.13.dist-info/WHEEL,sha256=SmOxYU7pzNKBqASvQJ7DjX3XGUF92lrGhMb3R6_iiqI,91
|
43
|
+
oafuncs-0.0.98.13.dist-info/top_level.txt,sha256=bgC35QkXbN4EmPHEveg_xGIZ5i9NNPYWqtJqaKqTPsQ,8
|
44
|
+
oafuncs-0.0.98.13.dist-info/RECORD,,
|
File without changes
|
File without changes
|
File without changes
|