shancx 1.9.33.226__tar.gz → 1.9.33.227__tar.gz

{shancx-1.9.33.226 → shancx-1.9.33.227}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: shancx
-Version: 1.9.33.226
+Version: 1.9.33.227
 Summary: A simple timer decorator
 Home-page: https://gitee.com/shancx
 Author: shancx

{shancx-1.9.33.226 → shancx-1.9.33.227}/setup.py

@@ -3,7 +3,7 @@
 from setuptools import setup, find_packages
 setup(
     name="shancx",
-    version="1.9.33.226",
+    version="1.9.33.227",
     packages=find_packages(),
     description="A simple timer decorator",
     long_description=open("README.md").read(),

{shancx-1.9.33.226 → shancx-1.9.33.227}/shancx/Read.py

@@ -1,8 +1,14 @@
+
 import xarray as xr
 import numpy as np
 import gzip
 import tempfile
 import os
+backend_kwargs = {
+    'errors': 'ignore',  
+    # 'errors': 'raise'   
+    # 'errors': 'warn'    
+}
 def getGrib(filepath):
     tmp_path = filepath 
     if filepath.endswith('.gz'):
@@ -11,7 +17,7 @@ def getGrib(filepath):
                 tmp.write(f.read())
                 tmp_path = tmp.name   
     try:
-        ds = xr.open_dataset(tmp_path, engine='cfgrib')
+        ds = xr.open_dataset(tmp_path, engine='cfgrib', backend_kwargs=backend_kwargs)
         var_name = list(ds.data_vars)[0]
         data_var = ds[var_name] 
         data_values = np.where(data_var.values < -900, np.nan, data_var.values)
@@ -32,16 +38,15 @@ def getGrib(filepath):
                     lons = var.longitude.values
                     break
             else:
-                raise ValueError("未找到经纬度坐标")
-        lat_range = (float(np.nanmin(lats)), float(np.nanmax(lats)))
-        lon_range = (float(np.nanmin(lons)), float(np.nanmax(lons)))        
+                raise ValueError("未找到经纬度坐标")       
         return {
-            'data': data,
+            'data': data_values,
             'lats': lats,
             'lons': lons 
         }        
     finally:
         os.unlink(tmp_path)
+
 """
 input_Path = './Composite_00.50_20250701-010044.grib2.gz'
 result = getGrib(input_Path)
@@ -50,15 +55,17 @@ result = getGrib(input_Path)
 import pygrib
 import numpy as np
 import pandas as pd
-def readGrib(file_path, target_param=None):
+import gzip
+import tempfile
+def readGrib(filepath, target_param=None):
+    filepath= filepath  
+    if filepath.endswith('.gz'):
+        with gzip.open(filepath, 'rb') as f:
+            with tempfile.NamedTemporaryFile(suffix='.grib2', delete=False) as tmp:
+                tmp.write(f.read())
+                filepath = tmp.name  
     try:
-        with pygrib.open(file_path) as grbs:    
-            file_path= file_path  
-            if filepath.endswith('.gz'):
-                with gzip.open(filepath, 'rb') as f:
-                    with tempfile.NamedTemporaryFile(suffix='.grib2', delete=False) as tmp:
-                        tmp.write(f.read())
-                        file_path = tmp.name       
+        with pygrib.open(filepath) as grbs:         
             if target_param:
                 try:
                     grb = grbs.select(shortName=target_param)[0]
@@ -99,6 +106,7 @@ def readGrib(file_path, target_param=None):
     except Exception as e:
         print(f"GRIB读取错误: {str(e)}")
         return None
+"""
 if __name__ == "__main__":
     path = "/mnt/wtx_weather_forecast/CMA_DATA/NAFP/EC/C1D/2024/2024112720/ECMFC1D_PRTY_1_2024112712_GLB_1_2.grib2"
     result = readGrib(path)    
@@ -107,5 +115,6 @@ if __name__ == "__main__":
         print("经度范围:", np.min(result['lons']), "~", np.max(result['lons']))
         print("纬度范围:", np.min(result['lats']), "~", np.max(result['lats']))
         print("参数单位:", result['metadata']['units'])
+"""
   
 

{shancx-1.9.33.226 → shancx-1.9.33.227}/shancx.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: shancx
-Version: 1.9.33.226
+Version: 1.9.33.227
 Summary: A simple timer decorator
 Home-page: https://gitee.com/shancx
 Author: shancx