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
Tutorials and Examples
Library Reference
Bibliography
Web resources

Developers Guide

Contributing Guide
Release Notes

Repository
Suggest edit
Open issue

.ipynb

Clutter detection using the Gabella approach

Contents

Read the data
Apply filter
Plot results

Clutter detection using the Gabella approach#

[1]:

import matplotlib.pyplot as plt
import numpy as np
import wradlib as wrl
import warnings

warnings.filterwarnings("ignore")
try:
    get_ipython().run_line_magic("matplotlib inline")
except:
    plt.ion()
import numpy as np

Read the data#

[2]:

filename = wrl.util.get_wradlib_data_file("misc/polar_dBZ_fbg.gz")
data = np.loadtxt(filename)
data = wrl.georef.create_xarray_dataarray(data, rf=0.001).wrl.georef.georeference()
data

Downloading file 'misc/polar_dBZ_fbg.gz' from 'https://github.com/wradlib/wradlib-data/raw/pooch/data/misc/polar_dBZ_fbg.gz' to '/home/runner/work/wradlib/wradlib/wradlib-data'.

[2]:

<xarray.DataArray (azimuth: 360, range: 128)>
array([[  7.01,   8.28,   5.46, ..., -10.  , -10.  , -10.  ],
       [  7.01,   7.01,   7.01, ..., -10.  , -10.  , -10.  ],
       [  7.01,   7.67,   5.46, ..., -10.  , -10.  , -10.  ],
       ...,
       [  7.01,   8.28,   5.46, ...,  14.87,  17.76,  15.51],
       [  7.67,   7.67,   5.46, ...,   0.64,   8.83,   4.53],
       [  8.83,   7.67,   4.53, ...,  -8.99,  -4.17,  -8.99]])
Coordinates: (12/14)
  * azimuth      (azimuth) float64 0.5 1.5 2.5 3.5 ... 356.5 357.5 358.5 359.5
  * range        (range) float64 500.0 1.5e+03 2.5e+03 ... 1.265e+05 1.275e+05
    elevation    (azimuth) float64 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0 0.0
    longitude    float64 0.0
    latitude     float64 0.0
    altitude     float64 0.0
    ...           ...
    y            (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    z            (azimuth, range) float64 0.0147 0.1323 0.3675 ... 940.8 955.7
    gr           (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    rays         (azimuth, range) float64 0.5 0.5 0.5 0.5 ... 359.5 359.5 359.5
    bins         (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    spatial_ref  int64 0

xarray.DataArray

azimuth: 360
range: 128

7.01 8.28 5.46 7.01 5.46 6.28 ... -10.0 -10.0 -10.0 -8.99 -4.17 -8.99

array([[  7.01,   8.28,   5.46, ..., -10.  , -10.  , -10.  ],
       [  7.01,   7.01,   7.01, ..., -10.  , -10.  , -10.  ],
       [  7.01,   7.67,   5.46, ..., -10.  , -10.  , -10.  ],
       ...,
       [  7.01,   8.28,   5.46, ...,  14.87,  17.76,  15.51],
       [  7.67,   7.67,   5.46, ...,   0.64,   8.83,   4.53],
       [  8.83,   7.67,   4.53, ...,  -8.99,  -4.17,  -8.99]])

Coordinates: (14)

azimuth
(azimuth)
float64
0.5 1.5 2.5 ... 357.5 358.5 359.5
```
array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5])
```

range

(range)

float64

500.0 1.5e+03 ... 1.275e+05

array([   500.,   1500.,   2500.,   3500.,   4500.,   5500.,   6500.,   7500.,
         8500.,   9500.,  10500.,  11500.,  12500.,  13500.,  14500.,  15500.,
        16500.,  17500.,  18500.,  19500.,  20500.,  21500.,  22500.,  23500.,
        24500.,  25500.,  26500.,  27500.,  28500.,  29500.,  30500.,  31500.,
        32500.,  33500.,  34500.,  35500.,  36500.,  37500.,  38500.,  39500.,
        40500.,  41500.,  42500.,  43500.,  44500.,  45500.,  46500.,  47500.,
        48500.,  49500.,  50500.,  51500.,  52500.,  53500.,  54500.,  55500.,
        56500.,  57500.,  58500.,  59500.,  60500.,  61500.,  62500.,  63500.,
        64500.,  65500.,  66500.,  67500.,  68500.,  69500.,  70500.,  71500.,
        72500.,  73500.,  74500.,  75500.,  76500.,  77500.,  78500.,  79500.,
        80500.,  81500.,  82500.,  83500.,  84500.,  85500.,  86500.,  87500.,
        88500.,  89500.,  90500.,  91500.,  92500.,  93500.,  94500.,  95500.,
        96500.,  97500.,  98500.,  99500., 100500., 101500., 102500., 103500.,
       104500., 105500., 106500., 107500., 108500., 109500., 110500., 111500.,
       112500., 113500., 114500., 115500., 116500., 117500., 118500., 119500.,
       120500., 121500., 122500., 123500., 124500., 125500., 126500., 127500.])

elevation

(azimuth)

float64

0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0.])

longitude
()
float64
0.0
```
array(0.)
```
latitude
()
float64
0.0
```
array(0.)
```
altitude
()
float64
0.0
```
array(0.)
```

sweep_mode

()

<U20

'azimuth_surveillance'

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

x

(azimuth, range)

float64

4.363 13.09 ... -1.113e+03

array([[ 4.36326774e+00,  1.30898031e+01,  2.18163381e+01, ...,
         1.09510071e+03,  1.10382533e+03,  1.11254992e+03],
       [ 1.30884741e+01,  3.92654221e+01,  6.54423689e+01, ...,
         3.28496856e+03,  3.31113976e+03,  3.33731087e+03],
       [ 2.18096937e+01,  6.54290804e+01,  1.09048465e+02, ...,
         5.47383577e+03,  5.51744558e+03,  5.56105525e+03],
       ...,
       [-2.18096937e+01, -6.54290804e+01, -1.09048465e+02, ...,
        -5.47383577e+03, -5.51744558e+03, -5.56105525e+03],
       [-1.30884741e+01, -3.92654221e+01, -6.54423689e+01, ...,
        -3.28496856e+03, -3.31113976e+03, -3.33731087e+03],
       [-4.36326774e+00, -1.30898031e+01, -2.18163381e+01, ...,
        -1.09510071e+03, -1.10382533e+03, -1.11254992e+03]])

y

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   499.98096096,   1499.94286904,   2499.90473565, ...,
        125486.1123282 , 126485.85478572, 127485.59374664],
       [   499.82866191,   1499.48597191,   2499.14324045, ...,
        125447.88804277, 126447.32596906, 127446.7603998 ],
       [   499.52411022,   1498.57231683,   2497.62048201, ...,
        125371.45111541, 126370.28007198, 127369.10553514],
       ...,
       [   499.52411022,   1498.57231683,   2497.62048201, ...,
        125371.45111541, 126370.28007198, 127369.10553514],
       [   499.82866191,   1499.48597191,   2499.14324045, ...,
        125447.88804277, 126447.32596906, 127446.7603998 ],
       [   499.98096096,   1499.94286904,   2499.90473565, ...,
        125486.1123282 , 126485.85478572, 127485.59374664]])

z

(azimuth, range)

float64

0.0147 0.1323 ... 940.8 955.7

array([[1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       ...,
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02]])

gr

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       ...,
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225]])

rays

(azimuth, range)

float64

0.5 0.5 0.5 ... 359.5 359.5 359.5

array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],
       [  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],
       [  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],
       ...,
       [357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],
       [358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],
       [359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]])

bins

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       ...,
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.]])

spatial_ref
()
int64
0
crs_wkt :
PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6378137,0,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Modified Azimuthal Equidistant",ID["EPSG",9832]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",0,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 :
6378137.0
inverse_flattening :
0.0
reference_ellipsoid_name :
unknown
longitude_of_prime_meridian :
0.0
prime_meridian_name :
Greenwich
geographic_crs_name :
unknown
horizontal_datum_name :
unknown
projected_crs_name :
unknown
grid_mapping_name :
azimuthal_equidistant
latitude_of_projection_origin :
0.0
longitude_of_projection_origin :
0.0
false_easting :
0.0
false_northing :
0.0
```
array(0)
```

Indexes: (2)

azimuth

PandasIndex

PandasIndex(Index([  0.5,   1.5,   2.5,   3.5,   4.5,   5.5,   6.5,   7.5,   8.5,   9.5,
       ...
       350.5, 351.5, 352.5, 353.5, 354.5, 355.5, 356.5, 357.5, 358.5, 359.5],
      dtype='float64', name='azimuth', length=360))

range

PandasIndex

PandasIndex(Index([   500.0,   1500.0,   2500.0,   3500.0,   4500.0,   5500.0,   6500.0,
         7500.0,   8500.0,   9500.0,
       ...
       118500.0, 119500.0, 120500.0, 121500.0, 122500.0, 123500.0, 124500.0,
       125500.0, 126500.0, 127500.0],
      dtype='float64', name='range', length=128))

Attributes: (0)

Apply filter#

[3]:

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

[3]:

<xarray.DataArray 'filter_gabella' (azimuth: 360, range: 128)>
array([[False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       [False, False, False, ..., False, False, False],
       ...,
       [False, False, False, ...,  True,  True,  True],
       [False, False, False, ...,  True,  True,  True],
       [False, False, False, ..., False, False, False]])
Coordinates: (12/14)
  * azimuth      (azimuth) float64 0.5 1.5 2.5 3.5 ... 356.5 357.5 358.5 359.5
  * range        (range) float64 500.0 1.5e+03 2.5e+03 ... 1.265e+05 1.275e+05
    elevation    (azimuth) float64 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0 0.0
    longitude    float64 0.0
    latitude     float64 0.0
    altitude     float64 0.0
    ...           ...
    y            (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    z            (azimuth, range) float64 0.0147 0.1323 0.3675 ... 940.8 955.7
    gr           (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    rays         (azimuth, range) float64 0.5 0.5 0.5 0.5 ... 359.5 359.5 359.5
    bins         (azimuth, range) float64 500.0 1.5e+03 ... 1.265e+05 1.275e+05
    spatial_ref  int64 0

xarray.DataArray

'filter_gabella'

azimuth: 360
range: 128

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, ...,  True,  True,  True],
       [False, False, False, ...,  True,  True,  True],
       [False, False, False, ..., False, False, False]])

Coordinates: (14)

azimuth
(azimuth)
float64
0.5 1.5 2.5 ... 357.5 358.5 359.5
```
array([  0.5,   1.5,   2.5, ..., 357.5, 358.5, 359.5])
```

range

(range)

float64

500.0 1.5e+03 ... 1.275e+05

array([   500.,   1500.,   2500.,   3500.,   4500.,   5500.,   6500.,   7500.,
         8500.,   9500.,  10500.,  11500.,  12500.,  13500.,  14500.,  15500.,
        16500.,  17500.,  18500.,  19500.,  20500.,  21500.,  22500.,  23500.,
        24500.,  25500.,  26500.,  27500.,  28500.,  29500.,  30500.,  31500.,
        32500.,  33500.,  34500.,  35500.,  36500.,  37500.,  38500.,  39500.,
        40500.,  41500.,  42500.,  43500.,  44500.,  45500.,  46500.,  47500.,
        48500.,  49500.,  50500.,  51500.,  52500.,  53500.,  54500.,  55500.,
        56500.,  57500.,  58500.,  59500.,  60500.,  61500.,  62500.,  63500.,
        64500.,  65500.,  66500.,  67500.,  68500.,  69500.,  70500.,  71500.,
        72500.,  73500.,  74500.,  75500.,  76500.,  77500.,  78500.,  79500.,
        80500.,  81500.,  82500.,  83500.,  84500.,  85500.,  86500.,  87500.,
        88500.,  89500.,  90500.,  91500.,  92500.,  93500.,  94500.,  95500.,
        96500.,  97500.,  98500.,  99500., 100500., 101500., 102500., 103500.,
       104500., 105500., 106500., 107500., 108500., 109500., 110500., 111500.,
       112500., 113500., 114500., 115500., 116500., 117500., 118500., 119500.,
       120500., 121500., 122500., 123500., 124500., 125500., 126500., 127500.])

elevation

(azimuth)

float64

0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0.])

longitude
()
float64
0.0
```
array(0.)
```
latitude
()
float64
0.0
```
array(0.)
```
altitude
()
float64
0.0
```
array(0.)
```

sweep_mode

()

<U20

'azimuth_surveillance'

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

x

(azimuth, range)

float64

4.363 13.09 ... -1.113e+03

array([[ 4.36326774e+00,  1.30898031e+01,  2.18163381e+01, ...,
         1.09510071e+03,  1.10382533e+03,  1.11254992e+03],
       [ 1.30884741e+01,  3.92654221e+01,  6.54423689e+01, ...,
         3.28496856e+03,  3.31113976e+03,  3.33731087e+03],
       [ 2.18096937e+01,  6.54290804e+01,  1.09048465e+02, ...,
         5.47383577e+03,  5.51744558e+03,  5.56105525e+03],
       ...,
       [-2.18096937e+01, -6.54290804e+01, -1.09048465e+02, ...,
        -5.47383577e+03, -5.51744558e+03, -5.56105525e+03],
       [-1.30884741e+01, -3.92654221e+01, -6.54423689e+01, ...,
        -3.28496856e+03, -3.31113976e+03, -3.33731087e+03],
       [-4.36326774e+00, -1.30898031e+01, -2.18163381e+01, ...,
        -1.09510071e+03, -1.10382533e+03, -1.11254992e+03]])

y

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   499.98096096,   1499.94286904,   2499.90473565, ...,
        125486.1123282 , 126485.85478572, 127485.59374664],
       [   499.82866191,   1499.48597191,   2499.14324045, ...,
        125447.88804277, 126447.32596906, 127446.7603998 ],
       [   499.52411022,   1498.57231683,   2497.62048201, ...,
        125371.45111541, 126370.28007198, 127369.10553514],
       ...,
       [   499.52411022,   1498.57231683,   2497.62048201, ...,
        125371.45111541, 126370.28007198, 127369.10553514],
       [   499.82866191,   1499.48597191,   2499.14324045, ...,
        125447.88804277, 126447.32596906, 127446.7603998 ],
       [   499.98096096,   1499.94286904,   2499.90473565, ...,
        125486.1123282 , 126485.85478572, 127485.59374664]])

z

(azimuth, range)

float64

0.0147 0.1323 ... 940.8 955.7

array([[1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       ...,
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02],
       [1.46986488e-02, 1.32287843e-01, 3.67466228e-01, ...,
        9.25979202e+02, 9.40793815e+02, 9.55725978e+02]])

gr

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       ...,
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225],
       [   499.99999942,   1499.99998444,   2499.99992798, ...,
        125490.89063679, 126490.67116289, 127490.44819225]])

rays

(azimuth, range)

float64

0.5 0.5 0.5 ... 359.5 359.5 359.5

array([[  0.5,   0.5,   0.5, ...,   0.5,   0.5,   0.5],
       [  1.5,   1.5,   1.5, ...,   1.5,   1.5,   1.5],
       [  2.5,   2.5,   2.5, ...,   2.5,   2.5,   2.5],
       ...,
       [357.5, 357.5, 357.5, ..., 357.5, 357.5, 357.5],
       [358.5, 358.5, 358.5, ..., 358.5, 358.5, 358.5],
       [359.5, 359.5, 359.5, ..., 359.5, 359.5, 359.5]])

bins

(azimuth, range)

float64

500.0 1.5e+03 ... 1.275e+05

array([[   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       ...,
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.],
       [   500.,   1500.,   2500., ..., 125500., 126500., 127500.]])

spatial_ref
()
int64
0
crs_wkt :
PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6378137,0,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["Modified Azimuthal Equidistant",ID["EPSG",9832]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",0,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 :
6378137.0
inverse_flattening :
0.0
reference_ellipsoid_name :
unknown
longitude_of_prime_meridian :
0.0
prime_meridian_name :
Greenwich
geographic_crs_name :
unknown
horizontal_datum_name :
unknown
projected_crs_name :
unknown
grid_mapping_name :
azimuthal_equidistant
latitude_of_projection_origin :
0.0
longitude_of_projection_origin :
0.0
false_easting :
0.0
false_northing :
0.0
```
array(0)
```

Indexes: (2)

azimuth

PandasIndex

PandasIndex(Index([  0.5,   1.5,   2.5,   3.5,   4.5,   5.5,   6.5,   7.5,   8.5,   9.5,
       ...
       350.5, 351.5, 352.5, 353.5, 354.5, 355.5, 356.5, 357.5, 358.5, 359.5],
      dtype='float64', name='azimuth', length=360))

range

PandasIndex

PandasIndex(Index([   500.0,   1500.0,   2500.0,   3500.0,   4500.0,   5500.0,   6500.0,
         7500.0,   8500.0,   9500.0,
       ...
       118500.0, 119500.0, 120500.0, 121500.0, 122500.0, 123500.0, 124500.0,
       125500.0, 126500.0, 127500.0],
      dtype='float64', name='range', length=128))

Attributes: (0)

Plot results#

[4]:

fig = plt.figure(figsize=(12, 8))
ax1 = fig.add_subplot(121)
pm = data.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")

[4]:

Text(0.5, 1.0, 'Cluttermap')

../../_images/notebooks_classify_clutter_gabella_8_1.png

previous

Clutter and Echo Classification

next

Clutter detection by using space-born cloud images

Contents

Read the data
Apply filter
Plot results

By Wradlib Community

© Copyright 2011-2023, wradlib developers.