GPU accelerated image processing for everyone
By Robert Haase, Kisha Sivanathan
Takes two label maps, and counts for every label in label map 1 how many labels overlap with it in label map 2.
The resulting map is generated from the label map 1 by replacing the labels with the respective count.
Categories: Labels, Measurements, Visualisation
Availability: Available in Fiji by activating the update sites clij and clij2. This function is part of clij2_-2.5.0.1.jar.
Ext.CLIJ2_labelOverlapCountMap(Image label_map1, Image label_map2, Image overlap_count_map_destination);
// init CLIJ and GPU import net.haesleinhuepf.clij2.CLIJ2; import net.haesleinhuepf.clij.clearcl.ClearCLBuffer; CLIJ2 clij2 = CLIJ2.getInstance(); // get input parameters ClearCLBuffer label_map1 = clij2.push(label_map1ImagePlus); ClearCLBuffer label_map2 = clij2.push(label_map2ImagePlus); overlap_count_map_destination = clij2.create(label_map1);
// Execute operation on GPU clij2.labelOverlapCountMap(label_map1, label_map2, overlap_count_map_destination);
// show result overlap_count_map_destinationImagePlus = clij2.pull(overlap_count_map_destination); overlap_count_map_destinationImagePlus.show(); // cleanup memory on GPU clij2.release(label_map1); clij2.release(label_map2); clij2.release(overlap_count_map_destination);
% init CLIJ and GPU clij2 = init_clatlab(); % get input parameters label_map1 = clij2.pushMat(label_map1_matrix); label_map2 = clij2.pushMat(label_map2_matrix); overlap_count_map_destination = clij2.create(label_map1);
% Execute operation on GPU clij2.labelOverlapCountMap(label_map1, label_map2, overlap_count_map_destination);
% show result overlap_count_map_destination = clij2.pullMat(overlap_count_map_destination) % cleanup memory on GPU clij2.release(label_map1); clij2.release(label_map2); clij2.release(overlap_count_map_destination);
// init CLIJ and GPU importClass(net.haesleinhuepf.clicy.CLICY); importClass(Packages.icy.main.Icy); clij2 = CLICY.getInstance(); // get input parameters label_map1_sequence = getSequence(); label_map1 = clij2.pushSequence(label_map1_sequence); label_map2_sequence = getSequence(); label_map2 = clij2.pushSequence(label_map2_sequence); overlap_count_map_destination = clij2.create(label_map1);
// Execute operation on GPU clij2.labelOverlapCountMap(label_map1, label_map2, overlap_count_map_destination);
// show result overlap_count_map_destination_sequence = clij2.pullSequence(overlap_count_map_destination) Icy.addSequence(overlap_count_map_destination_sequence); // cleanup memory on GPU clij2.release(label_map1); clij2.release(label_map2); clij2.release(overlap_count_map_destination);