Transform image intensities based on histogram mapping. — histogramWarpImageIntensities • ANTsR

Apply B-spline 1-D maps to an input image for intensity warping.

histogramWarpImageIntensities(
  image,
  breakPoints = c(0.25, 0.5, 0.75),
  displacements = NULL,
  clampEndPoints = c(FALSE, FALSE),
  sdDisplacements = 0.05,
  transformDomainSize = 20
)

Arguments

image: input image.
breakPoints: parametric points at which the intensity transform displacements are specified between [0, 1]. Alternatively, a single number can be given and the sequence is linearly spaced in [0, 1].
displacements: displacements to define intensity warping. Length must be equal to the breakPoints. Alternatively, if NULL random displacements are chosen (random normal: mean = 0, sd = sdDisplacements).
clampEndPoints: specify non-zero intensity change at the ends of the histogram.
sdDisplacements: characterize the randomness of the intensity displacement.
transformDomainSize: Defines the sampling resolution of the B-spline warping.

Value

warped intensity image

Author

Tustison NJ

Examples


library( ANTsR )
image <- antsImageRead( getANTsRData( "r16" ) )
transformedImage <- histogramWarpImageIntensities( image, transformDomainSize = 10 )
rm(image); gc()
#>            used  (Mb) gc trigger  (Mb) limit (Mb) max used  (Mb)
#> Ncells  3479516 185.9    5814731 310.6         NA  5814731 310.6
#> Vcells 11378691  86.9   22785357 173.9     102400 22784809 173.9
rm(transformedImage); gc()
#>            used  (Mb) gc trigger  (Mb) limit (Mb) max used  (Mb)
#> Ncells  3479421 185.9    5814731 310.6         NA  5814731 310.6
#> Vcells 11378569  86.9   22785357 173.9     102400 22784809 173.9