Probably the easiest thing would be to just extend the code in Images and submit a PR (the code is not very complicated).
However, you can do very fancy things with MapInfo objects. This is untested, but it should be close: immutable TwoColormap <: MapInfo colormap1 colormap2 end function map!(dest, mapi::TwoColormap, src::Tuple{AbstractArray,AbstractArray}) img1, img2 = src for I in eachindex(dest) dest[I] = clamp(RGBmapi.colormap1[img1[I]] + mapi.colormap2[img2[I]]) end dest end --Tim On Thursday, November 12, 2015 07:23:07 AM Tobias Knopp wrote: > Hi, > > I am using the OverlayImage type from the Images.jl package to overlay two > different grayscale images (tomographic data). > If I understand it correctly OverlayImage is restricted to colormaps that > go from black to a certain RGB value. Has anybody an idea how this could be > extended to Colormaps provided by Colors.jl? > > So my need is: > Input: two 3D datasets (FloatingPoint) + two Colormaps + WindowWidth > WindowLevel for each > Output: Combined 3D dataset as RGBA values. > > Thanks > > Tobias