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)> Size: 369kB
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)
  * range       (range) float64 1kB 500.0 1.5e+03 ... 1.265e+05 1.275e+05
  * azimuth     (azimuth) float64 3kB 0.5 1.5 2.5 3.5 ... 357.5 358.5 359.5
    elevation   (azimuth) float64 3kB 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
    longitude   float64 8B 0.0
    latitude    float64 8B 0.0
    altitude    float64 8B 0.0
    ...          ...
    y           (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    z           (azimuth, range) float64 369kB 0.0147 0.1323 ... 940.8 955.7
    gr          (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    rays        (azimuth, range) float64 369kB 0.5 0.5 0.5 ... 359.5 359.5 359.5
    bins        (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    crs_wkt     int64 8B 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)

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.])

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])
```

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

standard_name :: east_west_distance_from_radar
units :: meters

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

standard_name :: north_south_distance_from_radar
units :: meters

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

standard_name :: height_above_ground
units :: meters

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

standard_name :: distance_from_radar
units :: meters

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.]])

crs_wkt
()
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["Azimuthal Equidistant",ID["EPSG",1125]],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)

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))

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))

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)> Size: 46kB
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)
  * range       (range) float64 1kB 500.0 1.5e+03 ... 1.265e+05 1.275e+05
  * azimuth     (azimuth) float64 3kB 0.5 1.5 2.5 3.5 ... 357.5 358.5 359.5
    elevation   (azimuth) float64 3kB 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
    longitude   float64 8B 0.0
    latitude    float64 8B 0.0
    altitude    float64 8B 0.0
    ...          ...
    y           (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    z           (azimuth, range) float64 369kB 0.0147 0.1323 ... 940.8 955.7
    gr          (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    rays        (azimuth, range) float64 369kB 0.5 0.5 0.5 ... 359.5 359.5 359.5
    bins        (azimuth, range) float64 369kB 500.0 1.5e+03 ... 1.275e+05
    crs_wkt     int64 8B 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)

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.])

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])
```

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

standard_name :: east_west_distance_from_radar
units :: meters

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

standard_name :: north_south_distance_from_radar
units :: meters

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

standard_name :: height_above_ground
units :: meters

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

standard_name :: distance_from_radar
units :: meters

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.]])

crs_wkt
()
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["Azimuthal Equidistant",ID["EPSG",1125]],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)

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))

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))

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.