def quick_cr(r_list: Volume_T, resolution: Optional[tuple] = None) -> Grid: r""" Calculate composite reflectivity Paramters --------- r_list: list of cinrad.datastruct.Radial Returns ------- l2_obj: cinrad.datastruct.Grid composite reflectivity """ r_data = list() for i in r_list: r, x, y = grid_2d(i.data, i.lon, i.lat, resolution=resolution) r_data.append(r) cr = np.max(r_data, axis=0) x, y = np.meshgrid(x, y) l2_obj = Grid( np.ma.array(cr, mask=(cr <= 0)), i.drange, i.reso, i.code, i.name, i.scantime, "CR", i.stp["lon"], i.stp["lat"], x, y, **i.scan_info ) return l2_obj
def get_data(self) -> Union[Grid, Radial]: if self.radial_flag: lon, lat = self.projection() return Radial( self.data, self.max_range, self.el, self.reso, self.code, self.name, self.scantime, self.dtype, self.stationlon, self.stationlat, lon, lat, ) else: return Grid( self.data, self.max_range, self.reso, self.code, self.name, self.scantime, self.dtype, self.stationlon, self.stationlat, self.lon, self.lat, )
def get_data(self, level=0) -> Grid: x, y = np.meshgrid(self.lon, self.lat) dtype = self.product_name if self.data.ndim == 2: ret = self.data else: ret = self.data[level] if self.product_name == '3DREF': dtype = 'CR' grid = Grid(ret, np.nan, np.nan, 'SWAN', 'SWAN', self.data_time, dtype, 0, 0, x, y) return grid
def get_data(self, level=0) -> Grid: x, y = np.meshgrid(self.lon, self.lat) dtype = self.product_name if self.data.ndim == 2: ret = self.data else: ret = self.data[level] if self.product_name == "3DREF": dtype = "CR" grid = Grid(ret, np.nan, np.nan, "SWAN", "SWAN", self.data_time, dtype, 0, 0, x, y) return grid
def __call__(self, step: Number_T) -> Grid: x, y = self._process_grid(step, step) grid = self._map_points(x, y) grid = np.ma.masked_outside(grid, 0.1, 100) return Grid( grid, np.nan, step, "RADMAP", "RADMAP", self.scan_time, self.dtype, 0, 0, x, y, )
def quick_cr(r_list: RList) -> Grid: r''' Calculate composite reflectivity Paramters --------- r_list: list of cinrad.datastruct.Radial Returns ------- l2_obj: cinrad.datastruct.Grid composite reflectivity ''' r_data = list() for i in r_list: r, x, y = grid_2d(i.data, i.lon, i.lat) r_data.append(r) cr = np.max(r_data, axis=0) x, y = np.meshgrid(x, y) l2_obj = Grid(np.ma.array(cr, mask=(cr <= 0)), i.drange, 1, i.code, i.name, i.scantime, 'CR', x, y) return l2_obj
def get_data(self) -> Grid: x, y = np.meshgrid(self.lon, self.lat) grid = Grid(self.data, np.nan, np.nan, 'SWAN', 'SWAN', self.data_time, self.product_name, 0, 0, x, y) return grid