Julia package of https://github.com/mfuderer/colorResources
Note
This package goes with the guideline paper published in Magnetic Resonance in Medicine, https://doi.org/10.1002/mrm.30290, Color-map recommendation for MR relaxomtry maps .
It is to be used to display quantitative maps of T1, R1, T2, R2, T2*, R2*, T1rho and R1rho.
The package contains the Lipari color map, the Navia color map as well as the logarithm-processing as referred to by the publication.
In essence, following commands are needed:
# x is the image to be displayed, (loLev,upLev) is the relevant range of the values
# we take T1 as an example
using QMRIColors
cmap,imClip = relaxationColorMap("T1",x,loLev,upLev)
VIEW(imClip,c=cmap,clim=(loLev,upLev))
Important
In the above, "VIEW" refers to your favorite viewing software (CairoMakie, Plots, PyPlot, ...) : [see example](https://magneticresonanceimaging.github.io/QMRIColors.jl/dev/example/)
The relaxationColorMap() function selects the appropriate color-map (in this case, Lipari for T1) and applies logarithm-processing on that color map. The resulting colormap is returned.
The "imclip" image should retain the distinction between the "invalid" value of 0 (which is to be displayed as black) and "valid" values that are below loLev (these should be shown in a dark shade of blue). Depending on the calibration of your screen, you may or may not see that distinction.
A simplified processing can be done by
using QMRIColors
cmap = relaxationColorMap("T1")
VIEW(x,c=cmap) # Do not use range limits here!
Again, "VIEW" refers to your favorite viewing software.