mergechannels 0.1.1

No project description provided

mergechannels

This project was originally conceived because I often find myself wanting to apply and blend colormaps to images while working from Python, and for historical reasons many of my favorite colormaps are distributed as FIJI lookup tables. I also care about things likes speed and memory usage (e.g., not casting large arrays to floating point dtypes, not creating multiple whole arrays just to add them together), so I was interested in seeing if I could at least match matplotlib's colormapping performance with my own hand-rolled colorizer in Rust.

Installation

Install pre-compiled binaries from PyPI:

pip install mergechannels

Build from source on your own machine:

pip install git+https://github.com/zacswider/mergechannels.git

Usage

NOTE: scikit-image is not a dependency of this project, but is used in the examples below to fetch images.

colorize a single image

from skimage import data
import mergechannels as mc

img = data.camera()
colorized = mc.apply_color_map(img, 'Red/Green')
fig, axes = plt.subplots(1, 2, figsize=(6, 3), dpi=150)
for ax in axes: ax.axis('off')
(a, b) = axes
a.imshow(img, cmap='gray')
b.imshow(colorized)
plt.show()
print(colorized.shape, colorized.dtype)
>> (512, 512, 3) uint8

apply a different colormap to each channel

from skimage import data
import matplotlib.pyplot as plt
import mergechannels as mc

cells, nuclei = data.cells3d().max(axis=0)
cells = to_uint8(cells)     # normalize your own arrays, mergechannels doesn't currently handle this
nuclei = to_uint8(nuclei)   # normalize your own arrays, mergechannels doesn't currently handle this

fig, axes = plt.subplots(2, 2, figsize=(6, 6), dpi=150)
for ax in axes.ravel(): ax.axis('off')
(a, b, c, d) = axes.ravel()
a.imshow(cells, cmap='gray')
b.imshow(nuclei, cmap='gray')
c.imshow(
    mc.merge(
        [cells, nuclei],
        ['Orange Hot', 'Cyan Hot'],  # maximum blending is the default
    ),
)
d.imshow(
    mc.merge(
        [cells, nuclei],
        ['I Blue', 'I Forest'],
        blending='min',     # use minimum blending with inverted colormaps
    ),
)
fig.tight_layout()
plt.show()

Roadmap

mergechannels is currently incredibly simple. It can apply one or more colormaps to one or more 2D 8-bit images and that's it.

• Add support for any numerical dtype
• Add support for 3D images
• Add option to return any colormap as a matplotlib colormap
• Add support for directly passing matplotlib colormaps instead of colormap names
• Parallelize colormap application on large images (if it's helpful)
• Add option to overlay binary or instance masks onto colorized images

Acknowledgements

There are other great colormapping libraries available (e.g., microfilm, cmap) that are more feature-rich than this one, but which don't address my goals. My hope is that this project can fill an un-met niche and otherwise maintain full compatibility with these and similar libraries.

This project incorporates a number of colormaps that were hand-crafted by Christophe Leterrier and were originally distributed here under the MIT license: https://github.com/cleterrier/ChrisLUTs