image =	localEnergy3dImg (inImg3d,inHalfKnlSizeX,inHalfKnlSizeY,inHalfKnlSizeZ)
image =	localEnergy3dImg (inImg3d,inHalfKnlSizeX,inHalfKnlSizeY,inHalfKnlSizeZ,inOptNbClasses)

Detailed Description

local energy on a 3d image

This algorithm computes for each pixel of the output image associated local energy on a box neighbourhood of the input image.

Given an input image $InImg$ , box kernel half sizes $N_x = InHalfKnlSizeX$ , $N_y = InHakfKnlSizeY$ and $N_z = InHalfKnlSizeZ$ and a number of classes $n$ , the output image values are given by :

$OutImg[x, y, z] = \sum_{i=0}^{n}{p^2(i)}$

where:

equals to , with the number of pixels in the kernel whose values belong to the interval of the bin of the histogram computed on the box kernel, with the following parameters:
- number of classes = ,
- $min_{histo}$ = minimum value of the whole input image,
- $max_{histo}$ = maximum value of the whole input image.

The number of classes is an optional parameter of the algorithm. Its default value equals to 16. If the number of classes specified by the user exceeds the maximum allowed number of classes given the input image data type and dynamic range, it is automatically adjusted. For instance:

a number of classes initialized to value 16 will be automatically adjusted to value 2 for a binary input image,
a number of classes initialized to value 256 will be automatically adjusted to value 10 for an integral input image whose dynamic range equals to .

The borders of the input image are handled by padding pixels with a mirror reflection of the border pixels in the input image (see Border policy for more details).

Here is an example of an output image computed from the LocalEnergy3dImg algorithm on a 8-bits grey level, with a kernel of size 3x4x5 pixels and 8 classes:

Example of Python code :

Example imports

import PyIPSDK

import PyIPSDK.IPSDKIPLStats as stats

Code Example

    # opening of input image
    inImg = PyIPSDK.loadTiffImageFile(inputImgPath)
    
    # local entropy 2d image computation
    outImg = stats.localEnergy3dImg(inImg, 3, 4, 5, 8)

Example of C++ code :

Example informations

Header file

#include <IPSDKIPL/IPSDKIPLStats/Processor/LocalEnergy3dImg/LocalEnergy3dImg.h>

Code Example

// compute local energy on input image

ImagePtr pOutImg = localEnergy3dImg(pInImg, inHalfKnlSizeX, inHalfKnlSizeY, inHalfKnlSizeZ, nbClasses);