ImageDev

Global Analysis

This example shows how to extract global measurements from an image.
First, this example applies the IntensityStatistics algorithm on a color input image. It computes several statistics from the input image intensities. Each color channel is processed as a separate image and the example shows how to access to the result of each channel. Results are printed in the standard output.

Computing color intensity statistics...
- Average intensity [R, G, B] = [197.06, 186.25, 181.10]
- Standard Deviation = [R, G, B] = [10.05, 16.00, 7.79]
- Red intensity range = [151, 212]
- Green intensity range = [0, 207]
- Blue intensity range = [124, 197]

Note: Two indexes must be provided to access to each statistics result. The first index refers to the channel in case of a color or multi-channel input, like in this example. The second is unused here and refers to the time step in case of a time series input.

Then, the color image is binarized with the ColorThresholding algorithm. The total area occupied by the phase between thresholds is given by the Area2d algorithm. The number of detected objects is computed with the ObjectCount algorithm.

Color thresholding and binary analysis...
- Total area of detected objects = 46256
- Phase fraction of detected objects = 17.71 %
- Number of detected objects = 65 


#include <ImageDev/ImageDev.h>
#include <ioformat/IOFormat.h>

using namespace imagedev;
using namespace ioformat;
using namespace iolink;

int
main()
{
    // ImageDev library initialization
    if ( isInitialized() == false )
        imagedev::init();

    // Open a tif file to analyze
    auto imageInput = ioformat::readImage( std::string( IMAGEDEVDATA_IMAGES_FOLDER ) + "micro_nucleus.tif" );

    // Extract intensity statistics from an RGB image
    std::cout << "Computing color intensity statistics... " << std::endl;
    auto statistics = intensityStatistics( imageInput, IntensityStatistics::RangeMode::MIN_MAX, { 0, 255 } );
    std::cout << "  - Average intensity [R, G, B] = [" << statistics->mean( 0, 0 ) << ", " << statistics->mean( 1, 0 )
              << ", " << statistics->mean( 2, 0 ) << "]" << std::endl;
    std::cout << "  - Standard Deviation = [R, G, B] = [" << statistics->standardDeviation( 0, 0 ) << ", "
              << statistics->standardDeviation( 1, 0 ) << ", " << statistics->standardDeviation( 2, 0 ) << "]"
              << std::endl;
    std::cout << "  - Red intensity range = [" << statistics->minimum( 0, 0 ) << ", " << statistics->maximum( 0, 0 )
              << "]" << std::endl;
    std::cout << "  - Green intensity range = [" << statistics->minimum( 1, 0 ) << ", " << statistics->maximum( 1, 0 )
              << "]" << std::endl;
    std::cout << "  - Blue intensity range = [" << statistics->minimum( 2, 0 ) << ", " << statistics->maximum( 2, 0 )
              << "]" << std::endl;

    // Threshold image pixels having an intensity lower than 190 in each RGB channel
    std::cout << "Color thresholding and binary analysis... " << std::endl;
    auto imageOutput = colorThresholding( imageInput, { 0, 190 }, { 0, 190 }, { 0, 190 } );

    // Measure the object area
    auto area = area2d( imageOutput );
    std::cout << "  - Total area of detected objects = " << area->area( 0 ) << std::endl;
    std::cout << "  - Phase fraction of detected objects = " << ( 100.0 * area->areaFraction( 0 ) ) << " %"
              << std::endl;

    // Count the number of detected objects
    auto count = objectCount( imageOutput );
    std::cout << "  - Number of detected objects = " << count.outputMeasurement->count( 0 ) << std::endl;
    // Rescale binary the image intensities between 0 and 255 for visualization
    // and save it
    imageOutput = rescaleIntensity( imageOutput,
                                    RescaleIntensity::OutputType::UNSIGNED_INTEGER_8_BIT,
                                    RescaleIntensity::RangeMode::MIN_MAX,
                                    { 2, 98 },
                                    { 0, 255 },
                                    { 0, 255 } );
    ioformat::writeView( imageInput, R"(T04_01_output.png)" );

    // ImageDev library finalization
    imagedev::finish();

    return 0;
}


See also