pysteps.cascade.decomposition.decomposition_fft

pysteps.cascade.decomposition.decomposition_fft

pysteps.cascade.decomposition.decomposition_fft(field, bp_filter, **kwargs)

Decompose a two-dimensional input field into multiple spatial scales by using the Fast Fourier Transform (FFT) and a set of bandpass filters.

Parameters
field: array_like

Two-dimensional array containing the input field. All values are required to be finite.

bp_filter: dict

A filter returned by a method implemented in pysteps.cascade.bandpass_filters.

Returns
out: ndarray

A dictionary described in the module documentation. The number of cascade levels is determined from the filter (see pysteps.cascade.bandpass_filters).

Other Parameters
fft_method: str or tuple

A string or a (function,kwargs) tuple defining the FFT method to use (see pysteps.utils.interface.get_method()). Defaults to “numpy”. This option is not used if input_domain and output_domain are both set to “spectral”.

normalize: bool

If True, normalize the cascade levels to zero mean and unit variance. Requires that compute_stats is True. Implies that compute_stats is True. Defaults to False.

mask: array_like

Optional mask to use for computing the statistics for the cascade levels. Pixels with mask==False are excluded from the computations. This option is not used if output domain is “spectral”.

input_domain: {“spatial”, “spectral”}

The domain of the input field. If “spectral”, the input is assumed to be in the spectral domain. Defaults to “spatial”.

output_domain: {“spatial”, “spectral”}

If “spatial”, the output cascade levels are transformed back to the spatial domain by using the inverse FFT. If “spectral”, the cascade is kept in the spectral domain. Defaults to “spatial”.

compute_stats: bool

If True, the output dictionary contains the keys “means” and “stds” for the mean and standard deviation of each output cascade level. Defaults to False.

compact_output: bool

Applicable if output_domain is “spectral”. If set to True, only the parts of the Fourier spectrum with non-negligible filter weights are stored. Defaults to False.