Median 3d

bool  =	median3dImg_allow8bitsCompression (inImg,inHalfKnlSizeX,inHalfKnlSizeY,inHalfKnlSizeZ,outImg)
image =	median3dImg (inImg,inHalfKnlSizeX,inHalfKnlSizeY,inHalfKnlSizeZ)

Detailed Description

Median filter on a 3d image.

The median filter is a non-linear filter used to reduce impulsive noise in an image while preserving edges.

The median filter computes, for each pixel of input image, the median of its neighbouring pixels.

On output image values are given by:

where :

• function applied to a set of cardinal ( is an odd positive integer) returns the element of the collection containing the elements of the set sorted by increasing value
• is defined by InHalfKnlSizeX attribute
• is defined by InHalfKnlSizeY attribute
• is defined by InHalfKnlSizeZ attribute

Input and output images must have same size.

The algorithm is optimized if at least one of the following conditions is satisfied:

• and input image pixels are encoded on integers
• the kernel contains up to 125 pixels

If none of the previous conditions are satisfied, the user has all the same the possibility to execute a fast 2d median filter, by setting the attribute "InOutOptAllow8bitsCompression" to value true (this attribute is set to false by default). In this case, he will get a resulting image that will be an approximation of a real median (the resulting image will contain a maximum of 256 different grey level values).

Here are 2 examples of a median filter operation applied to a noised 8-bits grey levels input image:

• with :

• with :

See also

http://en.wikipedia.org/wiki/Median_filter

Example of Python code :

Example imports

import PyIPSDK
import PyIPSDK.IPSDKIPLFiltering as filter

Code Example

    # opening of input images
    inImg = PyIPSDK.loadTiffImageFile(inputImgPath)
    
    # median filter 3d computation
    outImg = filter.median3dImg(inImg, 3, 3, 3)

Example of C++ code :

Example informations

Header file

#include <IPSDKIPL/IPSDKIPLFiltering/Processor/Median3dImg/Median3dImg.h>

Code Example

        // opening input image
        ImagePtr pInImg = loadTiffImageFile(inputImgPath);

        if(pInImg->getBufferType() != inImageBufferType)
            pInImg = util::convertImg(pInImg, inImageBufferType);

        ImageGeometryPtr pOutImageGeometry =
            geometry3d(outImageBufferType, pInImg->getSizeX(), pInImg->getSizeY(), pInImg->getSizeZ());

        boost::shared_ptr<MemoryImage> pOutImg =
            boost::make_shared<MemoryImage>();
        pOutImg->init(*pOutImageGeometry);

        // compute median3d of input image
        median3dImg(pInImg, inHalfKnlSizeX, inHalfKnlSizeY, inHalfKnlSizeZ, pOutImg);