ENH: Add DASK_GAUSSIAN scaling and fix ITK_GAUSSIAN #40
Conversation
|
Example output for comparison
|
Changes: - Adds `DASK_GAUSSIAN` scaling option which combines dask_image Gaussian blurring with ITK's ResampleImageFilter for downsampling - Resolves issue in ITK_GAUSSIAN where kernel radius was improperly set from local block metadata, resulting severely different images from DASK_GAUSSIAN blurring results - Adds tests and update baseline images
There was a problem hiding this comment.
@tbirdso awesome!!
Can we use something from dask-image for the resampling so we do not have the itk dependency for this method? E.g. affine_transform, map_coordinates
The `DASK_GAUSSIAN` method previously relied on ITK for calculating output metadata at each scale and resampling after dask_image blurring. This commit updates the `DASK_GAUSSIAN` procedure to manually compute output metadata without ITK and to use the dask_image `affine_transform` method for downsampling.
|
@thewtex Please see updated changes substituting Does |
|
That's really cool to see. Thanks for sharing and working on this Matt and Tom! 😄 |
Changes:
DASK_GAUSSIANscaling option which combines dask_image Gaussianblurring with ITK's ResampleImageFilter for downsampling
resulting in severely different images from DASK_GAUSSIAN blurring results
Note that Gaussian kernels with same sigma value are still generated slightly differently between ITK and dask_image, so baseline image magnitudes differ by a small amount.