LSST Applications g042eb84c57+730a74494b,g04e9c324dd+8c5ae1fdc5,g134cb467dc+1f1e3e7524,g199a45376c+0ba108daf9,g1fd858c14a+fa7d31856b,g210f2d0738+f66ac109ec,g262e1987ae+83a3acc0e5,g29ae962dfc+d856a2cb1f,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+a1e0c9f713,g47891489e3+0d594cb711,g4d44eb3520+c57ec8f3ed,g4d7b6aa1c5+f66ac109ec,g53246c7159+8c5ae1fdc5,g56a1a4eaf3+fd7ad03fde,g64539dfbff+f66ac109ec,g67b6fd64d1+0d594cb711,g67fd3c3899+f66ac109ec,g6985122a63+0d594cb711,g74acd417e5+3098891321,g786e29fd12+668abc6043,g81db2e9a8d+98e2ab9f28,g87389fa792+8856018cbb,g89139ef638+0d594cb711,g8d7436a09f+80fda9ce03,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+a8a29bda4b,g99822b682c+e3ec3c61f9,g9d5c6a246b+0d5dac0c3d,ga41d0fce20+9243b26dd2,gbf99507273+8c5ae1fdc5,gd7ef33dd92+0d594cb711,gdab6d2f7ff+3098891321,ge410e46f29+0d594cb711,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.38
LSST Data Management Base Package
|
Classes | |
class | MultiResolutionSupport |
Functions | |
np.ndarray | bspline_convolve (np.ndarray image, int scale) |
int | get_starlet_scales (Sequence[int] image_shape, int|None scales=None) |
np.ndarray | starlet_transform (np.ndarray image, int|None scales=None, int generation=2, Callable|None convolve2d=None) |
np.ndarray | multiband_starlet_transform (np.ndarray image, int|None scales=None, int generation=2, Callable|None convolve2d=None) |
np.ndarray | starlet_reconstruction (np.ndarray starlets, int generation=2, Callable|None convolve2d=None) |
np.ndarray | multiband_starlet_reconstruction (np.ndarray starlets, int generation=2, Callable|None convolve2d=None) |
MultiResolutionSupport | get_multiresolution_support (np.ndarray image, np.ndarray starlets, np.floating sigma, float sigma_scaling=3, float epsilon=1e-1, int max_iter=20, str image_type="ground") |
np.ndarray | apply_wavelet_denoising (np.ndarray image, np.floating|None sigma=None, float sigma_scaling=3, float epsilon=1e-1, int max_iter=20, str image_type="ground", bool positive=True) |
np.ndarray lsst.scarlet.lite.wavelet.apply_wavelet_denoising | ( | np.ndarray | image, |
np.floating | None | sigma = None, | ||
float | sigma_scaling = 3, | ||
float | epsilon = 1e-1, | ||
int | max_iter = 20, | ||
str | image_type = "ground", | ||
bool | positive = True ) |
Apply wavelet denoising Uses the algorithm and notation from Starck et al. 2011, section 4.1 Parameters ---------- image: The image to denoise sigma: The standard deviation of the image sigma_scaling: The threshold in units of sigma to declare a coefficient significant epsilon: Convergence criteria for determining the support max_iter: The maximum number of iterations. This applies to both finding the support and the denoising loop. image_type: The type of image that is being used. This should be "ground" for ground based images with wide PSFs or "space" for images from space-based telescopes with a narrow PSF. positive: Whether or not the expected result should be positive Returns ------- result: The resulting denoised image after `max_iter` iterations.
Definition at line 336 of file wavelet.py.
np.ndarray lsst.scarlet.lite.wavelet.bspline_convolve | ( | np.ndarray | image, |
int | scale ) |
Convolve an image with a bspline at a given scale. This uses the spline `h1d = np.array([1.0 / 16, 1.0 / 4, 3.0 / 8, 1.0 / 4, 1.0 / 16])` from Starck et al. 2011. Parameters ---------- image: The 2D image or wavelet coefficients to convolve. scale: The wavelet scale for the convolution. This sets the spacing between adjacent pixels with the spline. Returns ------- result: The result of convolving the `image` with the spline.
Definition at line 36 of file wavelet.py.
MultiResolutionSupport lsst.scarlet.lite.wavelet.get_multiresolution_support | ( | np.ndarray | image, |
np.ndarray | starlets, | ||
np.floating | sigma, | ||
float | sigma_scaling = 3, | ||
float | epsilon = 1e-1, | ||
int | max_iter = 20, | ||
str | image_type = "ground" ) |
Calculate the multi-resolution support for a dictionary of starlet coefficients. This is different for ground and space based telescopes. For space-based telescopes the procedure in Starck and Murtagh 1998 iteratively calculates the multi-resolution support. For ground based images, where the PSF is much wider and there are no pixels with no signal at all scales, we use a modified method that estimates support at each scale independently. Parameters ---------- image: The image to transform into starlet coefficients. starlets: The starlet dictionary used to reconstruct `image` with dimension (scales+1, Ny, Nx). sigma: The standard deviation of the `image`. sigma_scaling: The multiple of `sigma` to use to calculate significance. Coefficients `w` where `|w| > K*sigma_j`, where `sigma_j` is standard deviation at the jth scale, are considered significant. epsilon: The convergence criteria of the algorithm. Once `|new_sigma_j - sigma_j|/new_sigma_j < epsilon` the algorithm has completed. max_iter: Maximum number of iterations to fit `sigma_j` at each scale. image_type: The type of image that is being used. This should be "ground" for ground based images with wide PSFs or "space" for images from space-based telescopes with a narrow PSF. Returns ------- M: Mask with significant coefficients in `starlets` set to `True`.
Definition at line 243 of file wavelet.py.
int lsst.scarlet.lite.wavelet.get_starlet_scales | ( | Sequence[int] | image_shape, |
int | None | scales = None ) |
Get the number of scales to use in the starlet transform. Parameters ---------- image_shape: The 2D shape of the image that is being transformed scales: The number of scales to transform with starlets. The total dimension of the starlet will have `scales+1` dimensions, since it will also hold the image at all scales higher than `scales`. Returns ------- result: Number of scales, adjusted for the size of the image.
Definition at line 80 of file wavelet.py.
np.ndarray lsst.scarlet.lite.wavelet.multiband_starlet_reconstruction | ( | np.ndarray | starlets, |
int | generation = 2, | ||
Callable | None | convolve2d = None ) |
Reconstruct a multiband image. See `starlet_reconstruction` for a description of the remainder of the parameters.
Definition at line 220 of file wavelet.py.
np.ndarray lsst.scarlet.lite.wavelet.multiband_starlet_transform | ( | np.ndarray | image, |
int | None | scales = None, | ||
int | generation = 2, | ||
Callable | None | convolve2d = None ) |
Perform a starlet transform of a multiband image. See `starlet_transform` for a description of the parameters.
Definition at line 161 of file wavelet.py.
np.ndarray lsst.scarlet.lite.wavelet.starlet_reconstruction | ( | np.ndarray | starlets, |
int | generation = 2, | ||
Callable | None | convolve2d = None ) |
Reconstruct an image from a dictionary of starlets Parameters ---------- starlets: The starlet dictionary used to reconstruct the image with dimension (scales+1, Ny, Nx). generation: The generation of the starlet transform (either ``1`` or ``2``). convolve2d: The filter function to use to convolve the image with starlets in 2D. Returns ------- image: The 2D image reconstructed from the input `starlet`.
Definition at line 183 of file wavelet.py.
np.ndarray lsst.scarlet.lite.wavelet.starlet_transform | ( | np.ndarray | image, |
int | None | scales = None, | ||
int | generation = 2, | ||
Callable | None | convolve2d = None ) |
Perform a starlet transform, or 2nd gen starlet transform. Parameters ---------- image: The image to transform into starlet coefficients. scales: The number of scale to transform with starlets. The total dimension of the starlet will have `scales+1` dimensions, since it will also hold the image at all scales higher than `scales`. generation: The generation of the transform. This must be `1` or `2`. convolve2d: The filter function to use to convolve the image with starlets in 2D. Returns ------- starlet: The starlet dictionary for the input `image`.
Definition at line 105 of file wavelet.py.