metradar 0.1.5__py3-none-any.whl → 0.1.8__py3-none-any.whl

metradar/__init__.py

@@ -1,7 +1,9 @@
 """
-NMC_MET_RADAR is a python package to read or retrieve radar data
+metradar is a python package to read or retrieve radar data
 
 """
 
 __author__ = "The R & D Center for Weather Forecasting Technology in NMC, CMA"
-__version__ = '0.1.1'
+
+
+__all__ = [s for s in dir() if not s.startswith('_')]

metradar/config.py

@@ -0,0 +1,53 @@
+# _*_ coding: utf-8 _*_
+
+# Copyright (c) 2026 NMC Developers.
+# Distributed under the terms of the GPL V3 License.
+
+"""
+Read configure file.
+"""
+
+
+
+import os
+import configparser
+from pathlib import Path
+
+
+def _get_config_dir():
+    """
+    Get default configuration directory.
+    """
+    config_dir = Path.home() / ".metradar"
+    config_dir.mkdir(parents=True, exist_ok=True)
+    return config_dir
+
+# Global Variables
+CONFIG_DIR = _get_config_dir()
+
+
+def _ConfigFetchError(BaseException):
+    pass
+
+
+def _get_config_from_rcfile(rc='config.ini'):
+    """
+    Get configure information from config_dk_met_io.ini file.
+    """
+
+    if not os.path.exists(rc):
+        print(rc + ' not exists!')
+        return None
+  
+    try:
+        config = configparser.ConfigParser()
+        config.read(rc,encoding='utf-8')
+    except IOError as e:
+        raise _ConfigFetchError(str(e))
+    except BaseException as e:
+        raise _ConfigFetchError(str(e))
+
+    return config
+
+# Global Variables
+CONFIG = _get_config_from_rcfile(CONFIG_DIR / 'config.ini')

metradar/core/__init__.py

@@ -0,0 +1,9 @@
+"""
+METRADAR is a python package to read or retrieve radar data
+
+"""
+
+__author__ = "The R & D Center for Weather Forecasting Technology in NMC, CMA"
+
+
+__all__ = [s for s in dir() if not s.startswith('_')]

metradar/{get_cross_section_from_pyart.py → core/get_cross_section.py}

@@ -1,13 +1,15 @@
+# _*_ coding: utf-8 _*_
+
 '''
 从pyart的radar object中获取任意垂直剖面的数据
 
 朱文剑
-2021.10.14
+
 
 '''
 
 import numpy as np
-from pyart.core import Radar
+from pyart.core import Radar,antenna_to_cartesian
 import math
 # import time
 
@@ -164,6 +166,7 @@ def get_cross_radar(radar:Radar,params):
                 #计算相对雷达的仰角
                 # ele = math.degrees(math.atan2(curz,curx))
                 # 这里必须用math.pow，如果直接用称号会有较大误差
+                # print('R=',R,' curz=',curz,' curx=',curx)
                 ele = math.asin((math.pow((curz+R),2) - math.pow(curx,2) - math.pow(R,2))/2/curx/R)/np.pi*180
 
                 if ele < baseel_R or ele > topel_R:
@@ -197,158 +200,3 @@ def get_cross_radar(radar:Radar,params):
                         nflag+=1
     
     return outdic
-
-# main function
-def get_cross_radar_old(radar:Radar,params):
-    # startangle=None,startrange=None,endangle=None,endrange=None
-    # startangle = 8.5
-    # startrange = 69.23 #km
-
-    # endangle = 18.6;
-    # endrange =  59.03 #km
-    startangle = params['crs_start_azi']
-    startrange = params['crs_start_range']
-    endangle = params['crs_end_azi']
-    endrange = params['crs_end_range']
-    toph = params['top_height']
-
-    startx = startrange*np.sin(np.radians(startangle))
-    starty = startrange*np.cos(np.radians(startangle))
-    endx   = endrange*np.sin(np.radians(endangle))
-    endy   = endrange*np.cos(np.radians(endangle))
-
-    baseh = 0; # km 垂直剖面的起始高度
-    # toph = 20; # km 垂直剖面的上限高度
-
-    validzlv=[]
-    validvlv=[]
-    validReles=[]
-    validDeles=[]
-    outdic=dict()
-
-    for il in range(radar.nsweeps):
-        # if sweep(il).RGates > 0
-        # if il !=1 and il !=3:
-        #     validzlv.append(il)
-        #     validReles.append(radar.get_elevation(il).mean())
-        # if il !=0 and il !=2:
-        #     validvlv.append(il)
-        #     validDeles.append(radar.get_elevation(il).mean())
-        validzlv.append(il)
-        validReles.append(round(radar.get_elevation(il).mean(),1))
-        validvlv.append(il)
-        validDeles.append(round(radar.get_elevation(il).mean(),1))
-
-    tol_length = np.sqrt(np.power((endx-startx),2) + np.power((endy - starty),2))
-
-    #求方程坐标，这里暂时没考虑X=0的情况，后面要完善这一点
-    [a,b] = linefunc(startx,starty,endx,endy)
-
-    # 获取反射率垂直剖面
-        
-    #构建变量，确定坐标维度
-    xreso = 0.25#km
-    if endx < startx:
-        xreso = -1*xreso
-
-    xnum = int(tol_length  / abs(xreso))
-    yreso = 0.25# km
-    ynum = int((toph - baseh)  / yreso)
-
-    
-
-    beamwidth = radar.get_azimuth(validzlv[-1])
-    beamwidth = np.diff(beamwidth)
-    beamwidth = np.mean(beamwidth[beamwidth>0])
-    topel_R = max(validReles) + beamwidth / 2
-    baseel_R = min(validReles) - beamwidth / 2
-
-    topel_D = max(validDeles) + beamwidth / 2
-    baseel_D = min(validDeles) - beamwidth / 2
-
-    tmp=[]
-    for nn in range(len(validReles)-1):
-        tmp.append(np.mean([validReles[nn],validReles[nn+1]]))
-
-    tmp_newRel = sorted(myconcat((tmp,validReles)))
-
-    tmp=[]
-    for nn in range(len(validDeles)-1):
-        tmp.append(np.mean([validDeles[nn],validDeles[nn+1]]))
-
-    tmp_newDel = sorted(tmp+validDeles)
-
-    newRel = [baseel_R+tmp_newRel+topel_R]
-    newDel = [baseel_D+tmp_newDel+topel_D]
-
-    # for varkey in radar.fields.keys():
-     
-    nflag=0
-    outdic=dict()
-    outdic['vertical_km'] = toph - baseh
-    outdic['horizontal_km'] = tol_length
-    outdic['xreso'] = xreso
-    outdic['yreso'] = yreso
-    
-    for varkey in radar.fields.keys():
-        outdic[varkey]=[]
-        outdic[varkey] = np.zeros([ynum,xnum],dtype='float')*np.nan
-
-    
-    # 纯 for 循环版本 =====================
-    valideles=[]
-    validlvl=[]
-    for varkey in radar.fields.keys():
-        print(varkey)
-        if varkey == 'reflectivity' :
-            valideles = validReles
-            validlvl = validzlv
-        elif varkey == 'velocity':
-            valideles = validDeles
-            validlvl = validvlv 
-        else:
-            # print('暂时仅支持回波强度和径向速度')
-            continue
-        
-        curfd = []
-        for il in range(radar.nsweeps):
-            curfd.append(radar.get_field(il,varkey))#sweep(targetsp).dbz(az,gates);
-
-        for iy in range(ynum):
-            for ix in range(xnum):
-                tmpcurx = startx + (ix * xreso) * math.cos(math.atan(a))
-                
-                tmpcury = a * tmpcurx + b
-                tmpaz = (90 - math.degrees(math.atan2(tmpcury,tmpcurx)) + 360)%360
-                
-                curx = math.sqrt(math.pow(tmpcurx,2) + math.pow(tmpcury,2))#startrange + xreso * ix 
-                cury = baseh + iy * yreso
-                
-                #计算相对雷达的仰角
-                ele = math.degrees(math.atan2(cury,curx))
-                # if ele < baseel_R: ele = baseel_R
-                # if ele > topel_R: continue
-                if ele < baseel_R or ele > topel_R:
-                    continue
-
-                tmp = list(abs(np.array(valideles) - ele))
-                idx = tmp.index(min(tmp))
-                idx_el = idx
-                
-                targetsp = validlvl[idx_el]
-                
-                #获取方位角角标
-                tmp = list(abs(radar.get_azimuth(targetsp) - tmpaz))
-                idx = tmp.index(min(tmp))
-                az = idx
-                
-                gates = math.ceil(curx/math.cos(math.radians(valideles[idx_el]))/(radar.range['meters_between_gates'] / 1000))
-                
-                if gates < radar.ngates:
-                    #disp(sweep(targetsp).dbz(az,gates));
-                    if curfd[targetsp].data[az,gates] > -30:
-                        outdic[varkey][iy,ix] = curfd[targetsp].data[az,gates]
-                        nflag+=1
-    
-    return outdic
-

metradar/{mosaic_merge.py → core/mosaic_merge.py}

@@ -1,5 +1,7 @@
+# _*_ coding: utf-8 _*_
+
 '''
-将不同层次的mosacdata stack在一起
+将不同层次的mosacdata stack在一起，主要用于多进程拼图时，提高效率
 ZhuWJ
 
 '''