Maximum.java

package fr.unistra.pelican.algorithms.arithmetic;

import fr.unistra.pelican.Algorithm;
import fr.unistra.pelican.AlgorithmException;
import fr.unistra.pelican.Image;
import fr.unistra.pelican.util.mask.MaskStack;

/**
 * Compute the supremum of two images, for each pixel max(inputImage1,
 * inputImage2). 
 * 
 * Works on double precision. 
 * The outputImage format is the same as inputImage1.
 *
 *	Mask management (by witz) :
 *		- computation occurs as if masked pixels were at 0.
 *		- the the output image mask is inputImage1 mask ORed with inputImage2 mask
 * 
 * @author  ?, Benjamin Perret
 */
public class Maximum extends Algorithm {

	/**
	 * First input image
	 */
	public Image inputImage1;

	/**
	 * Second input image
	 */
	public Image inputImage2;

	/**
	 * Supremum of inputs one and two
	 */
	public Image outputImage;

	/**
	 * Constructor
	 * 
	 */
	public Maximum() {

		super();
		super.inputs = "inputImage1,inputImage2";
		super.outputs = "outputImage";
		

	}

	/*
	 * (non-Javadoc)
	 * 
	 * @see fr.unistra.pelican.Algorithm#launch()
	 */
	public void launch() throws AlgorithmException {
		outputImage = inputImage1.copyImage(false);
		MaskStack mask = new MaskStack( MaskStack.OR );
		mask.push( inputImage1.getMask() );
		mask.push( inputImage2.getMask() );
		outputImage.setMask( mask );
		int size = inputImage1.size();

		boolean isHere1, isHere2;
		for ( int i = 0 ; i < size ; ++i ) {  

			isHere1 = inputImage1.isPresent(i);
			isHere2 = inputImage2.isPresent(i);
			double pixel1 = inputImage1.getPixelDouble(i);
			double pixel2 = inputImage2.getPixelDouble(i);
			if ( isHere1 && isHere2 ) 
				outputImage.setPixelDouble( i, pixel1 > pixel2 ? pixel1 : pixel2 );
			else if ( isHere1 && !isHere2 ) outputImage.setPixelDouble( i, pixel1 );
			else if ( !isHere1 && isHere2 ) outputImage.setPixelDouble( i, pixel2 );
			// if ( !isHere1 && !isHere2 ) do nothing.. 
		}
		
	}
	
	/**
	 * Compute the supremum of two images, for each pixel max(inputImage1,
	 * inputImage2). 
	 * 
	 * Works on double precision. 
	 * The outputImage format is the same as inputImage1.
	 * 
	 * @param inputImage1 First input image
	 * @param inputImage2 Second input image
	 * @return Supremum of inputs one and two
	 */
	public static Image exec(Image inputImage1, Image inputImage2) {
		return (Image)new Maximum().process(inputImage1, inputImage2);
	}
}