chore: update UAM & swan & Transparent img (#134)

* 优化对风雷产品支持

* 添加UAM产品参数解析

* fix transparent img

* fix:透明图投影

* fix: error to export pyart

* fix: close file after read SWAN

Author: pysoer

cinrad/io/export.py

@@ -88,7 +88,7 @@ def standard_data_to_pyart(f: StandardData, radius: int = 460) -> pyart.core.Rad
     fixed_angle["data"] = np.array(f.el)
 
     fields = {}
-    nscans = f.get_nscans()
+    # nscans = f.get_nscans()
 
     all_var = list()
     for lvl in f.data:
@@ -100,7 +100,7 @@ def standard_data_to_pyart(f: StandardData, radius: int = 460) -> pyart.core.Rad
             name = mapping[mom]
             dic = filemetadata(name)
             dic["_FillValue"] = pyart.config.get_fillvalue()
-            raw_arr = [f.get_raw(nel, radius, mom) for nel in range(nscans)]
+            raw_arr = [f.get_raw(nel, radius, mom) for nel in f.available_tilt(mom)]
             sel_arr = [i if not isinstance(i, tuple) else i[0] for i in raw_arr]
             moment_data = np.ma.vstack(sel_arr)
             dic["data"] = moment_data

cinrad/io/level3.py

@@ -217,7 +217,7 @@ def __init__(self, file: Any, product: Optional[str] = None):
             if not product
             else product
         )
-        for pname in ["CR", "3DREF", "反射率"]:
+        for pname in ["CR", "REF", "反射率"]:
             if pname in self.product_name:
                 self.product_name = "CR"
         start_lon = header["start_lon"][0]
@@ -235,6 +235,7 @@ def __init__(self, file: Any, product: Optional[str] = None):
         else:
             # Leave data unchanged because the scale and offset are unclear
             self.data = np.ma.masked_equal(out, 0)
+        f.close()
 
     def get_data(self) -> Dataset:
         r"""
@@ -421,6 +422,18 @@ def _cappi(buf):
         params["filled"] = filled
         return params
 
+    @staticmethod
+    def _uam(buf):
+        params = {}
+        iRange = np.frombuffer(buf.read(4), "i4")[0]
+        baseProd1 = np.frombuffer(buf.read(4), "i4")[0]
+        dataType1 = np.frombuffer(buf.read(4), "i4")[0]
+        params["Range"] = iRange
+        params["baseProd1"] = baseProd1
+        params["dataType1"] = dataType1
+        params["elevation"] = 0
+        return params
+    
     @staticmethod
     def _empty(buf):
         pass
@@ -434,6 +447,7 @@ def parse(cls, product_type, param_bytes):
             3: cls._cappi,
             51: cls._ppi,
             52: cls._ppi,
+            44: cls._uam,
             18: cls._empty,
         }
         params = {"elevation": 0}
@@ -1019,6 +1033,7 @@ def _parse_uam_fmt(self):
                 "site_longitude": self.stationlon,
                 "site_latitude": self.stationlat,
                 "task": self.task_name,
+                **self.params,
             },
         )
         ds["longitude"] = DataArray(lon[:, 0])

cinrad/visualize/utils.py

@@ -279,8 +279,8 @@ def create_geoaxes(
     from cartopy import __version__
 
     if style == "transparent":
-        ax = fig.add_axes([0, 0, 1, 1], projection=proj)
-        ax.set_aspect("equal")
+        ax = fig.add_axes([0, 0, 1, 1], projection=ccrs.GOOGLE_MERCATOR)
+        ax.set_aspect("auto")
     else:
         ax = fig.add_axes(GEOAXES_POS, projection=proj)
     if __version__ < "0.18":