Skip to main content

Ctrl+K

GitHub
PyPI
Openradar Discourse

GitHub
PyPI
Openradar Discourse

GitHub
PyPI
Openradar Discourse

Getting Started

Installation
wradlib 2.0 migration
Docker
Python IDE’s
Jupyter Notebooks

User Guide

Community
Core Features
Tutorials and Examples
RADOLAN Guide
Library Reference
Bibliography
Web resources

Developers Guide

Contributing Guide
Release Notes

Repository
Suggest edit
Open issue

.md

Clutter detection using the Gabella approach

Contents

Read the data
Apply filter
Plot results

Clutter detection using the Gabella approach#

import warnings

import matplotlib.pyplot as plt
import numpy as np
import wradlib as wrl
import wradlib_data
import xarray as xr

warnings.filterwarnings("ignore")

/home/docs/checkouts/readthedocs.org/user_builds/wradlib/conda/2.9.0/lib/python3.13/site-packages/tqdm/auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html
  from .autonotebook import tqdm as notebook_tqdm

Read the data#

filename = wradlib_data.DATASETS.fetch("hdf5/2014-08-10--182000.ppi.mvol")
swp = xr.open_dataset(filename, engine="gamic", group="sweep_0")
swp = swp.set_coords("sweep_mode")
swp = swp.wrl.georef.georeference()
display(swp)

<xarray.Dataset> Size: 33MB
Dimensions:            (azimuth: 360, range: 1000)
Coordinates: (12/15)
  * azimuth            (azimuth) float64 3kB 0.5054 1.522 2.527 ... 358.5 359.5
    elevation          (azimuth) float64 3kB 1.505 1.505 1.505 ... 1.505 1.505
    time               (azimuth) datetime64[ns] 3kB ...
  * range              (range) float32 4kB 50.0 150.0 ... 9.985e+04 9.995e+04
    x                  (azimuth, range) float64 3MB 0.4409 1.323 ... -866.6
    y                  (azimuth, range) float64 3MB 49.98 149.9 ... 9.988e+04
    ...                 ...
    bins               (azimuth, range) float32 1MB 50.0 150.0 ... 9.995e+04
    sweep_mode         <U20 80B 'azimuth_surveillance'
    longitude          float64 8B 7.072
    latitude           float64 8B 50.73
    altitude           float64 8B 99.5
    crs_wkt            int64 8B 0
Data variables: (12/17)
    KDP                (azimuth, range) float32 1MB ...
    PHIDP              (azimuth, range) float32 1MB ...
    DBZH               (azimuth, range) float32 1MB ...
    DBZV               (azimuth, range) float32 1MB ...
    RHOHV              (azimuth, range) float32 1MB ...
    DBTH               (azimuth, range) float32 1MB ...
    ...                 ...
    ZDR                (azimuth, range) float32 1MB ...
    sweep_number       int64 8B ...
    prt_mode           <U7 28B ...
    follow_mode        <U7 28B ...
    sweep_fixed_angle  float64 8B ...
    nyquist_velocity   object 8B ...
Attributes:
    source:       gamic
    pulse_width:  2

xarray.Dataset

Dimensions:
- azimuth: 360
- range: 1000

Coordinates: (15)

azimuth
(azimuth)
float64
0.5054 1.522 2.527 ... 358.5 359.5
standard_name :
ray_azimuth_angle
long_name :
azimuth_angle_from_true_north
units :
degrees
axis :
radial_azimuth_coordinate
```
array([  0.505371,   1.521606,   2.526855, ..., 357.50885 , 358.514099,
       359.502869], shape=(360,))
```
elevation
(azimuth)
float64
1.505 1.505 1.505 ... 1.505 1.505
standard_name :
ray_elevation_angle
long_name :
elevation_angle_from_horizontal_plane
units :
degrees
axis :
radial_elevation_coordinate
```
array([1.505127, 1.505127, 1.505127, ..., 1.505127, 1.505127, 1.505127],
      shape=(360,))
```
time
(azimuth)
datetime64[ns]
...
standard_name :
time
```
[360 values with dtype=datetime64[ns]]
```
range
(range)
float32
50.0 150.0 ... 9.985e+04 9.995e+04
units :
meters
standard_name :
projection_range_coordinate
long_name :
range_to_measurement_volume
axis :
radial_range_coordinate
meters_between_gates :
100.0
spacing_is_constant :
true
meters_to_center_of_first_gate :
50.0
```
array([5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04, 9.995e+04],
      shape=(1000,), dtype=float32)
```

x

(azimuth, range)

float64

0.4409 1.323 ... -865.7 -866.6

axis :: X
long_name :: Easting
standard_name :: projection_x_coordinate
units :: metre

array([[ 4.40856362e-01,  1.32256868e+00,  2.20428045e+00, ...,
         8.79196725e+02,  8.80077772e+02,  8.80958818e+02],
       [ 1.32722219e+00,  3.98166534e+00,  6.63610685e+00, ...,
         2.64686983e+03,  2.64952227e+03,  2.65217471e+03],
       [ 2.20359591e+00,  6.61078567e+00,  1.10179727e+01, ...,
         4.39461573e+03,  4.39901959e+03,  4.40342345e+03],
       ...,
       [-2.17247790e+00, -6.51743169e+00, -1.08623828e+01, ...,
        -4.33255731e+03, -4.33689899e+03, -4.34124066e+03],
       [-1.29608506e+00, -3.88825396e+00, -6.48042127e+00, ...,
        -2.58477326e+03, -2.58736347e+03, -2.58995368e+03],
       [-4.33668669e-01, -1.30100560e+00, -2.16834200e+00, ...,
        -8.64862358e+02, -8.65729041e+02, -8.66595722e+02]],
      shape=(360, 1000))

y

(azimuth, range)

float64

49.98 149.9 ... 9.978e+04 9.988e+04

axis :: Y
long_name :: Northing
standard_name :: projection_y_coordinate
units :: metre

array([[4.99802109e+01, 1.49940586e+02, 2.49900900e+02, ...,
        9.96751812e+04, 9.97750661e+04, 9.98749509e+04],
       [4.99645306e+01, 1.49893546e+02, 2.49822499e+02, ...,
        9.96439101e+04, 9.97437636e+04, 9.98436171e+04],
       [4.99335559e+01, 1.49800621e+02, 2.49667625e+02, ...,
        9.95821374e+04, 9.96819291e+04, 9.97817207e+04],
       ...,
       [4.99349195e+01, 1.49804712e+02, 2.49674443e+02, ...,
        9.95848567e+04, 9.96846511e+04, 9.97844454e+04],
       [4.99653480e+01, 1.49895998e+02, 2.49826585e+02, ...,
        9.96455402e+04, 9.97453954e+04, 9.98452505e+04],
       [4.99802738e+01, 1.49940775e+02, 2.49901214e+02, ...,
        9.96753066e+04, 9.97751916e+04, 9.98750766e+04]],
      shape=(360, 1000))

z

(azimuth, range)

float64

100.8 103.4 ... 3.309e+03 3.313e+03

standard_name :: height_above_ground
units :: metre

array([[ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       ...,
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169]], shape=(360, 1000))

gr

(azimuth, range)

float64

49.98 149.9 ... 9.978e+04 9.988e+04

standard_name :: distance_from_radar
units :: metre

array([[4.99827785e+01, 1.49947043e+02, 2.49911245e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827803e+01, 1.49947044e+02, 2.49911246e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827818e+01, 1.49947046e+02, 2.49911248e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       ...,
       [4.99827722e+01, 1.49947036e+02, 2.49911238e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827745e+01, 1.49947039e+02, 2.49911241e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827766e+01, 1.49947041e+02, 2.49911243e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04]],
      shape=(360, 1000))

rays

(azimuth, range)

float64

0.5054 0.5054 ... 359.5 359.5

standard_name :: ray_azimuth_angle
long_name :: azimuth_angle_from_true_north
units :: degrees
axis :: radial_azimuth_coordinate

array([[  0.50537109,   0.50537109,   0.50537109, ...,   0.50537109,
          0.50537109,   0.50537109],
       [  1.52160645,   1.52160645,   1.52160645, ...,   1.52160645,
          1.52160645,   1.52160645],
       [  2.52685547,   2.52685547,   2.52685547, ...,   2.52685547,
          2.52685547,   2.52685547],
       ...,
       [357.5088501 , 357.5088501 , 357.5088501 , ..., 357.5088501 ,
        357.5088501 , 357.5088501 ],
       [358.51409912, 358.51409912, 358.51409912, ..., 358.51409912,
        358.51409912, 358.51409912],
       [359.50286865, 359.50286865, 359.50286865, ..., 359.50286865,
        359.50286865, 359.50286865]], shape=(360, 1000))

bins

(azimuth, range)

float32

50.0 150.0 ... 9.985e+04 9.995e+04

units :: meters
standard_name :: projection_range_coordinate
long_name :: range_to_measurement_volume
axis :: radial_range_coordinate
meters_between_gates :: 100.0
spacing_is_constant :: true
meters_to_center_of_first_gate :: 50.0

array([[5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       ...,
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04]], shape=(360, 1000), dtype=float32)

sweep_mode

()

<U20

'azimuth_surveillance'

array('azimuth_surveillance', dtype='<U20')

longitude
()
float64
7.072
long_name :
longitude
units :
degrees_east
standard_name :
longitude
```
array(7.071663)
```
latitude
()
float64
50.73
long_name :
latitude
units :
degrees_north
positive :
up
standard_name :
latitude
```
array(50.73052)
```
altitude
()
float64
99.5
long_name :
altitude
units :
meters
standard_name :
altitude
```
array(99.5)
```
crs_wkt
()
int64
0
crs_wkt :
PROJCRS["AEQD centered at 50.73052N 7.071663E",BASEGEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Azimuthal Equidistant",ID["EPSG",1125]],PARAMETER["Latitude of natural origin",50.73052,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",7.071663,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]
semi_major_axis :
6378137.0
semi_minor_axis :
6356752.314245179
inverse_flattening :
298.257223563
reference_ellipsoid_name :
WGS 84
longitude_of_prime_meridian :
0.0
prime_meridian_name :
Greenwich
geographic_crs_name :
unknown
horizontal_datum_name :
World Geodetic System 1984
projected_crs_name :
AEQD centered at 50.73052N 7.071663E
grid_mapping_name :
azimuthal_equidistant
latitude_of_projection_origin :
50.73052
longitude_of_projection_origin :
7.071663
false_easting :
0.0
false_northing :
0.0
```
array(0)
```

Data variables: (17)
- KDP
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Specific differential phase HV
  standard_name :
  radar_specific_differential_phase_hv
  units :
  degrees per kilometer
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- PHIDP
  (azimuth, range)
  float32
  ...
  format :
  UV16
  long_name :
  Differential phase HV
  standard_name :
  radar_differential_phase_hv
  units :
  degrees
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- DBZH
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Equivalent reflectivity factor H
  standard_name :
  radar_equivalent_reflectivity_factor_h
  units :
  dBZ
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- DBZV
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Equivalent reflectivity factor V
  standard_name :
  radar_equivalent_reflectivity_factor_v
  units :
  dBZ
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- RHOHV
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Correlation coefficient HV
  standard_name :
  radar_correlation_coefficient_hv
  units :
  unitless
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- DBTH
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Total power H (uncorrected reflectivity)
  standard_name :
  radar_equivalent_reflectivity_factor_h
  units :
  dBZ
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- DBTV
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Total power V (uncorrected reflectivity)
  standard_name :
  radar_equivalent_reflectivity_factor_v
  units :
  dBZ
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- VRADH
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Radial velocity of scatterers away from instrument H
  standard_name :
  radial_velocity_of_scatterers_away_from_instrument_h
  units :
  meters per seconds
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- VRADV
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Radial velocity of scatterers away from instrument V
  standard_name :
  radial_velocity_of_scatterers_away_from_instrument_v
  units :
  meters per second
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- WRADH
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Doppler spectrum width H
  standard_name :
  radar_doppler_spectrum_width_h
  units :
  meters per seconds
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- WRADV
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Doppler spectrum width V
  standard_name :
  radar_doppler_spectrum_width_v
  units :
  meters per second
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- ZDR
  (azimuth, range)
  float32
  ...
  format :
  UV8
  long_name :
  Log differential reflectivity H/V
  standard_name :
  radar_differential_reflectivity_hv
  units :
  dB
  _Undetect :
  0.0
```
[360000 values with dtype=float32]
```
- sweep_number
  ()
  int64
  ...
```
[1 values with dtype=int64]
```
- prt_mode
  ()
  <U7
  ...
```
[1 values with dtype=<U7]
```
- follow_mode
  ()
  <U7
  ...
```
[1 values with dtype=<U7]
```
- sweep_fixed_angle
  ()
  float64
  ...
```
[1 values with dtype=float64]
```
- nyquist_velocity
  ()
  object
  ...
  standard_name :
  nyquist_velocity
  units :
  m s-1
```
[1 values with dtype=object]
```
Attributes: (2)
source :
gamic
pulse_width :
2

Apply filter#

clmap = swp.DBTH.wrl.classify.filter_gabella(
    wsize=5, thrsnorain=0.0, tr1=6.0, n_p=8, tr2=1.3
)
clmap

<xarray.DataArray 'filter_gabella' (azimuth: 360, range: 1000)> Size: 360kB
array([[False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       ...,
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False]], shape=(360, 1000))
Coordinates: (12/15)
  * azimuth     (azimuth) float64 3kB 0.5054 1.522 2.527 ... 357.5 358.5 359.5
    elevation   (azimuth) float64 3kB 1.505 1.505 1.505 ... 1.505 1.505 1.505
    time        (azimuth) datetime64[ns] 3kB ...
  * range       (range) float32 4kB 50.0 150.0 250.0 ... 9.985e+04 9.995e+04
    x           (azimuth, range) float64 3MB 0.4409 1.323 ... -865.7 -866.6
    y           (azimuth, range) float64 3MB 49.98 149.9 ... 9.978e+04 9.988e+04
    ...          ...
    bins        (azimuth, range) float32 1MB 50.0 150.0 ... 9.985e+04 9.995e+04
    sweep_mode  <U20 80B 'azimuth_surveillance'
    longitude   float64 8B 7.072
    latitude    float64 8B 50.73
    altitude    float64 8B 99.5
    crs_wkt     int64 8B 0
Attributes:
    standard_name:  clutter_map
    long_name:      Clutter Map
    short_name:     CMAP
    units:          unitless

xarray.DataArray

'filter_gabella'

azimuth: 360
range: 1000

False False False False False False ... False False False False False

array([[False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       ...,
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False]], shape=(360, 1000))

Coordinates: (15)

azimuth
(azimuth)
float64
0.5054 1.522 2.527 ... 358.5 359.5
standard_name :
ray_azimuth_angle
long_name :
azimuth_angle_from_true_north
units :
degrees
axis :
radial_azimuth_coordinate
```
array([  0.505371,   1.521606,   2.526855, ..., 357.50885 , 358.514099,
       359.502869], shape=(360,))
```
elevation
(azimuth)
float64
1.505 1.505 1.505 ... 1.505 1.505
standard_name :
ray_elevation_angle
long_name :
elevation_angle_from_horizontal_plane
units :
degrees
axis :
radial_elevation_coordinate
```
array([1.505127, 1.505127, 1.505127, ..., 1.505127, 1.505127, 1.505127],
      shape=(360,))
```
time
(azimuth)
datetime64[ns]
...
standard_name :
time
```
[360 values with dtype=datetime64[ns]]
```
range
(range)
float32
50.0 150.0 ... 9.985e+04 9.995e+04
units :
meters
standard_name :
projection_range_coordinate
long_name :
range_to_measurement_volume
axis :
radial_range_coordinate
meters_between_gates :
100.0
spacing_is_constant :
true
meters_to_center_of_first_gate :
50.0
```
array([5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04, 9.995e+04],
      shape=(1000,), dtype=float32)
```

x

(azimuth, range)

float64

0.4409 1.323 ... -865.7 -866.6

axis :: X
long_name :: Easting
standard_name :: projection_x_coordinate
units :: metre

array([[ 4.40856362e-01,  1.32256868e+00,  2.20428045e+00, ...,
         8.79196725e+02,  8.80077772e+02,  8.80958818e+02],
       [ 1.32722219e+00,  3.98166534e+00,  6.63610685e+00, ...,
         2.64686983e+03,  2.64952227e+03,  2.65217471e+03],
       [ 2.20359591e+00,  6.61078567e+00,  1.10179727e+01, ...,
         4.39461573e+03,  4.39901959e+03,  4.40342345e+03],
       ...,
       [-2.17247790e+00, -6.51743169e+00, -1.08623828e+01, ...,
        -4.33255731e+03, -4.33689899e+03, -4.34124066e+03],
       [-1.29608506e+00, -3.88825396e+00, -6.48042127e+00, ...,
        -2.58477326e+03, -2.58736347e+03, -2.58995368e+03],
       [-4.33668669e-01, -1.30100560e+00, -2.16834200e+00, ...,
        -8.64862358e+02, -8.65729041e+02, -8.66595722e+02]],
      shape=(360, 1000))

y

(azimuth, range)

float64

49.98 149.9 ... 9.978e+04 9.988e+04

axis :: Y
long_name :: Northing
standard_name :: projection_y_coordinate
units :: metre

array([[4.99802109e+01, 1.49940586e+02, 2.49900900e+02, ...,
        9.96751812e+04, 9.97750661e+04, 9.98749509e+04],
       [4.99645306e+01, 1.49893546e+02, 2.49822499e+02, ...,
        9.96439101e+04, 9.97437636e+04, 9.98436171e+04],
       [4.99335559e+01, 1.49800621e+02, 2.49667625e+02, ...,
        9.95821374e+04, 9.96819291e+04, 9.97817207e+04],
       ...,
       [4.99349195e+01, 1.49804712e+02, 2.49674443e+02, ...,
        9.95848567e+04, 9.96846511e+04, 9.97844454e+04],
       [4.99653480e+01, 1.49895998e+02, 2.49826585e+02, ...,
        9.96455402e+04, 9.97453954e+04, 9.98452505e+04],
       [4.99802738e+01, 1.49940775e+02, 2.49901214e+02, ...,
        9.96753066e+04, 9.97751916e+04, 9.98750766e+04]],
      shape=(360, 1000))

z

(azimuth, range)

float64

100.8 103.4 ... 3.309e+03 3.313e+03

standard_name :: height_above_ground
units :: metre

array([[ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       ...,
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169],
       [ 100.81346717,  103.44128457,  106.07027939, ..., 3305.14502024,
        3308.94609379, 3312.74837169]], shape=(360, 1000))

gr

(azimuth, range)

float64

49.98 149.9 ... 9.978e+04 9.988e+04

standard_name :: distance_from_radar
units :: metre

array([[4.99827785e+01, 1.49947043e+02, 2.49911245e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827803e+01, 1.49947044e+02, 2.49911246e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827818e+01, 1.49947046e+02, 2.49911248e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       ...,
       [4.99827722e+01, 1.49947036e+02, 2.49911238e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827745e+01, 1.49947039e+02, 2.49911241e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04],
       [4.99827766e+01, 1.49947041e+02, 2.49911243e+02, ...,
        9.96790593e+04, 9.97789481e+04, 9.98788368e+04]],
      shape=(360, 1000))

rays

(azimuth, range)

float64

0.5054 0.5054 ... 359.5 359.5

standard_name :: ray_azimuth_angle
long_name :: azimuth_angle_from_true_north
units :: degrees
axis :: radial_azimuth_coordinate

array([[  0.50537109,   0.50537109,   0.50537109, ...,   0.50537109,
          0.50537109,   0.50537109],
       [  1.52160645,   1.52160645,   1.52160645, ...,   1.52160645,
          1.52160645,   1.52160645],
       [  2.52685547,   2.52685547,   2.52685547, ...,   2.52685547,
          2.52685547,   2.52685547],
       ...,
       [357.5088501 , 357.5088501 , 357.5088501 , ..., 357.5088501 ,
        357.5088501 , 357.5088501 ],
       [358.51409912, 358.51409912, 358.51409912, ..., 358.51409912,
        358.51409912, 358.51409912],
       [359.50286865, 359.50286865, 359.50286865, ..., 359.50286865,
        359.50286865, 359.50286865]], shape=(360, 1000))

bins

(azimuth, range)

float32

50.0 150.0 ... 9.985e+04 9.995e+04

units :: meters
standard_name :: projection_range_coordinate
long_name :: range_to_measurement_volume
axis :: radial_range_coordinate
meters_between_gates :: 100.0
spacing_is_constant :: true
meters_to_center_of_first_gate :: 50.0

array([[5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       ...,
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04],
       [5.000e+01, 1.500e+02, 2.500e+02, ..., 9.975e+04, 9.985e+04,
        9.995e+04]], shape=(360, 1000), dtype=float32)

sweep_mode

()

<U20

'azimuth_surveillance'

array('azimuth_surveillance', dtype='<U20')

longitude
()
float64
7.072
long_name :
longitude
units :
degrees_east
standard_name :
longitude
```
array(7.071663)
```
latitude
()
float64
50.73
long_name :
latitude
units :
degrees_north
positive :
up
standard_name :
latitude
```
array(50.73052)
```
altitude
()
float64
99.5
long_name :
altitude
units :
meters
standard_name :
altitude
```
array(99.5)
```
crs_wkt
()
int64
0
crs_wkt :
PROJCRS["AEQD centered at 50.73052N 7.071663E",BASEGEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Azimuthal Equidistant",ID["EPSG",1125]],PARAMETER["Latitude of natural origin",50.73052,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",7.071663,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]
semi_major_axis :
6378137.0
semi_minor_axis :
6356752.314245179
inverse_flattening :
298.257223563
reference_ellipsoid_name :
WGS 84
longitude_of_prime_meridian :
0.0
prime_meridian_name :
Greenwich
geographic_crs_name :
unknown
horizontal_datum_name :
World Geodetic System 1984
projected_crs_name :
AEQD centered at 50.73052N 7.071663E
grid_mapping_name :
azimuthal_equidistant
latitude_of_projection_origin :
50.73052
longitude_of_projection_origin :
7.071663
false_easting :
0.0
false_northing :
0.0
```
array(0)
```

Attributes: (4)
standard_name :
clutter_map
long_name :
Clutter Map
short_name :
CMAP
units :
unitless

Plot results#

fig = plt.figure(figsize=(12, 8))
ax1 = fig.add_subplot(121)
pm = swp.DBTH.wrl.vis.plot(ax=ax1)
ax1.set_title("Reflectivity")
ax2 = fig.add_subplot(122)
pm = clmap.wrl.vis.plot(ax=ax2)
ax2.set_title("Cluttermap")
fig.tight_layout()

../../../_images/0cc8e14596868d9d1463f33ddbcd3dbbb6375622e313226fc98157738316e23a.png

previous

Clutter and Echo Classification

next

Clutter detection - Satellite Cloud Images

Contents

Read the data
Apply filter
Plot results

By Wradlib Community

© Copyright 2011-2026, wradlib developers.