Application: gvSIG desktop » gvSIG raster

/*

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 gvSIG Association.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.raster.tools.algorithm.layerdatatype;

import org.gvsig.i18n.Messages;

import org.gvsig.raster.tools.algorithm.base.RasterBaseAlgorithmLibrary;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of MaskthresholdAlgorithmLibrary library.

 * 

 * @author <a href="mailto:nachobrodin@gmail.com">Nacho Brodin</a>

 */

public class LayerDatatypeAlgorithmLibrary extends AbstractLibrary {

	public static final String         PROCESS_LABEL   = "MaskthresholdProcess";

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

    	//Registers the process and its parameters

    	RasterBaseAlgorithmLibrary.register(PROCESS_LABEL, LayerDatatypeProcess.class);

    	LayerDatatypeProcess.registerParameters();

        Messages.addResourceFamily(

            "org.gvsig.raster.tools.algorithm.maskthreshold", 

            LayerDatatypeAlgorithmLibrary.class.getClassLoader(), 

            LayerDatatypeAlgorithmLibrary.class.getClass().getName());

        //registerGeoProcess(new RasterReprojectAlgorithmLibrary());

    }

}

/* gvSIG. Geographic Information System of the Valencian Government

*

* Copyright (C) 2007-2008 Infrastructures and Transports Department

* of the Valencian Government (CIT)

* 

* This program is free software; you can redistribute it and/or

* modify it under the terms of the GNU General Public License

* as published by the Free Software Foundation; either version 2

* of the License, or (at your option) any later version.

* 

* This program is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

* GNU General Public License for more details.

* 

* You should have received a copy of the GNU General Public License

* along with this program; if not, write to the Free Software

* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, 

* MA  02110-1301, USA.

* 

*/

package org.gvsig.raster.tools.algorithm.layerdatatype;

import org.gvsig.raster.tools.algorithm.base.process.ProcessException;

/**

 * This exception is thrown if happen problems processing data with this algorithm.

 * 

 * @author Nacho Brodin (nachobrodin@gmail.com)

 *

 */

public class LayerDatatypeException extends ProcessException {

	private static final long serialVersionUID = -3022090543908771484L;

	public LayerDatatypeException(String msg){

		super(msg);

	}

	public LayerDatatypeException(String msg, Throwable e){

		super(msg, e);

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

*

* Copyright (C) 2007-2008 Infrastructures and Transports Department

* of the Valencian Government (CIT)

* 

* This program is free software; you can redistribute it and/or

* modify it under the terms of the GNU General Public License

* as published by the Free Software Foundation; either version 2

* of the License, or (at your option) any later version.

* 

* This program is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

* GNU General Public License for more details.

* 

* You should have received a copy of the GNU General Public License

* along with this program; if not, write to the Free Software

* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, 

* MA  02110-1301, USA.

* 

*/

package org.gvsig.raster.tools.algorithm.layerdatatype;

import java.util.HashMap;

import javax.swing.SwingUtilities;

import org.gvsig.fmap.dal.coverage.RasterLocator;

import org.gvsig.fmap.dal.coverage.dataset.Buffer;

import org.gvsig.fmap.dal.coverage.datastruct.NoData;

import org.gvsig.fmap.dal.coverage.exception.InvalidSetViewException;

import org.gvsig.fmap.dal.coverage.exception.ProcessInterruptedException;

import org.gvsig.fmap.dal.coverage.exception.RasterDriverException;

import org.gvsig.fmap.dal.coverage.store.RasterDataStore;

import org.gvsig.fmap.dal.coverage.store.RasterQuery;

import org.gvsig.i18n.Messages;

import org.gvsig.raster.tools.algorithm.base.RasterBaseAlgorithmLibrary;

import org.gvsig.raster.tools.algorithm.base.process.RasterProcess;

/**

 * Process to change a <code>RasterDataStore</code> of data type

 *

 * 10/12/2007

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public class LayerDatatypeProcess extends RasterProcess {

	public static String      RASTER_STORE1     = "RasterStore1";

	public static String      PATH              = "Path";

	public static String      FILENAME          = "FileName";

	public static String      DATATYPE          = "Datatype";

	public static String      ADJUST_DEC2INT    = "AdjustDec2Int";

	public static String      ADJUST_BIG2SMALL  = "AdjustBig2Small";

	public static String      TIME              = "Time";

	public static String[]    DEC2INT_OPTIONS   = new String[]{"Trunk", "Round", 

        "Ceil", "Floor"};

	public static String[]    BIG2SMALL_OPTIONS = new String[]{"Trunk", "Maxvalue", 

        "NoData"};

	public static String      RESULT_DSTORE   	= "RESULT_DSTORE";

	private RasterDataStore   store             = null;

	private String            filename          = null;

	private int               datatype          = 0;

	private int               dec2int           = 0;

	private int               big2smal          = 0;

	private long              milis             = 0;

	private int               percent           = 0;

	private NoData            byteNoData        = null;

	private NoData            shortNoData       = null;

	private NoData            intNoData         = null;

	private NoData            floatNoData       = null;

	public static void registerParameters() {

		RASTER_STORE1 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_STORE1, RasterDataStore.class);

		PATH = RasterBaseAlgorithmLibrary.registerInputParameter(PATH, String.class);

		DATATYPE = RasterBaseAlgorithmLibrary.registerInputParameter(DATATYPE, Integer.class);

		ADJUST_DEC2INT = RasterBaseAlgorithmLibrary.registerInputParameter(ADJUST_DEC2INT, Integer.class);

		ADJUST_BIG2SMALL = RasterBaseAlgorithmLibrary.registerInputParameter(ADJUST_BIG2SMALL, Integer.class);

		RESULT_DSTORE = RasterBaseAlgorithmLibrary.registerOutputParameter(RESULT_DSTORE, RasterDataStore.class);

		FILENAME = RasterBaseAlgorithmLibrary.registerOutputParameter(FILENAME, String.class);

		TIME = RasterBaseAlgorithmLibrary.registerOutputParameter(TIME, Long.class);

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.rastertools.RasterProcess#init()

	 */

	public void init() {

		store = getParam(RASTER_STORE1) != null ? (RasterDataStore)getParam(RASTER_STORE1) : null;

		filename = getStringParam(PATH);

		datatype = (Integer)getParam(DATATYPE);

		dec2int = (Integer)getParam(ADJUST_DEC2INT);

		big2smal = (Integer)getParam(ADJUST_BIG2SMALL);

		byteNoData = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(

				store.getBandCount(), Buffer.TYPE_BYTE);

		shortNoData = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(

				store.getBandCount(), Buffer.TYPE_SHORT);

		intNoData = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(

				store.getBandCount(), Buffer.TYPE_INT);

		floatNoData = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(

				store.getBandCount(), Buffer.TYPE_FLOAT);

	}

	public void process() throws ProcessInterruptedException {

		long t1 = new java.util.Date().getTime();

		insertLineLog(Messages.getText("layer_datatype"));

		try {

			if (store == null)

				throw new LayerDatatypeException("intput_not_valid");

			int inDataType = store.getDataType()[0];

			Buffer buf = RasterLocator.getManager().createBuffer(

					datatype, 

					(int)store.getWidth(), 

					(int)store.getHeight(), 

					store.getBandCount(), 

					true);

			RasterQuery query = RasterLocator.getManager().createQuery();

			query.setAllDrawableBands();

			query.setAreaOfInterest();

			Buffer sourceBuffer;

			try {

				sourceBuffer = store.query(query);

			} catch (RasterDriverException e) {

				throw new LayerDatatypeException("");

			} catch (InvalidSetViewException e) {

				throw new LayerDatatypeException("");

			}

			double v = 0;

			switch (inDataType) {

			case Buffer.TYPE_BYTE:

				for (int row = 0; row < buf.getHeight(); row++) {

					for (int col = 0; col < buf.getWidth(); col++) {

						for (int iBand = 0; iBand < buf.getBandCount(); iBand++) {

							writePixelFromByte(row, col, iBand, sourceBuffer.getElemByte(row, col, iBand), buf);

						}

					}

					percent = (int)((row * 100) / buf.getHeight());

				}

				break;

			case Buffer.TYPE_SHORT:

				for (int row = 0; row < buf.getHeight(); row++) {

					for (int col = 0; col < buf.getWidth(); col++) {

						for (int iBand = 0; iBand < buf.getBandCount(); iBand++) {

							writePixelFromShort(row, col, iBand, sourceBuffer.getElemShort(row, col, iBand), buf);

						}

					}

					percent = (int)((row * 100) / buf.getHeight());

				}

				break;

			case Buffer.TYPE_INT:

				for (int row = 0; row < buf.getHeight(); row++) {

					for (int col = 0; col < buf.getWidth(); col++) {

						for (int iBand = 0; iBand < buf.getBandCount(); iBand++) {

							writePixelFromInt(row, col, iBand, sourceBuffer.getElemInt(row, col, iBand), buf);

						}

					}

					percent = (int)((row * 100) / buf.getHeight());

				}

				break;

			case Buffer.TYPE_FLOAT:

				for (int row = 0; row < buf.getHeight(); row++) {

					for (int col = 0; col < buf.getWidth(); col++) {

						for (int iBand = 0; iBand < buf.getBandCount(); iBand++) {

							v = getDecimalValue(sourceBuffer.getElemFloat(row, col, iBand), inDataType);

							writePixelFromDouble(row, col, iBand, v, buf);

						}

					}

					percent = (int)((row * 100) / buf.getHeight());

				}

				break;

			case Buffer.TYPE_DOUBLE:

				for (int row = 0; row < buf.getHeight(); row++) {

					for (int col = 0; col < buf.getWidth(); col++) {

						for (int iBand = 0; iBand < buf.getBandCount(); iBand++) {

							v = getDecimalValue(sourceBuffer.getElemDouble(row, col, iBand), inDataType);

							writePixelFromDouble(row, col, iBand, v, buf);

						}

					}

					percent = (int)((row * 100) / buf.getHeight());

				}

				break;

			}

			super.exportRaster(filename, 

					buf, 

					sourceBuffer.getStore().getCellSize(), 

					sourceBuffer.getStore().getExtent().getULX(), 

					sourceBuffer.getStore().getExtent().getULY());

			long t2 = new java.util.Date().getTime();

			milis = t2 - t1;

			SwingUtilities.invokeLater(new Runnable() {

				public void run() {

					if (externalActions != null) {

						HashMap<String, Object> map = new HashMap<String, Object>();

						map.put(FILENAME, filename);

						map.put(TIME, new Long(milis));

						externalActions.end(map);

					}

				}

			});

		} catch (LayerDatatypeException e) {

			if (incrementableTask != null)

				incrementableTask.processFinalize();

			messageBoxError("error_changing_the_datatype", this, e);

		}

	}

	public double getDecimalValue(double value, int inDataType) {

		if(inDataType == Buffer.TYPE_FLOAT || inDataType == Buffer.TYPE_DOUBLE) {

			switch (dec2int) {

			case 0:return (double)((long)value);

			case 1:return Math.round(value);

			case 2:return Math.ceil(value);

			case 3:return Math.floor(value);

			}

		} 

		return value;

	}

	private void writePixelFromByte(int row, int col, int band, byte value, Buffer dst) {

		switch (dst.getDataType()) {

		case Buffer.TYPE_BYTE: dst.setElem(row, col, band, value);break;

		case Buffer.TYPE_SHORT: dst.setElem(row, col, band, (short)value);break;

		case Buffer.TYPE_INT: dst.setElem(row, col, band, (int)value);break;

		case Buffer.TYPE_FLOAT: dst.setElem(row, col, band, (float)value);break;

		case Buffer.TYPE_DOUBLE: dst.setElem(row, col, band, (double)value);break;

		}

	}

	private void writePixelFromShort(int row, int col, int band, short value, Buffer dst) {

		switch (dst.getDataType()) {

		case Buffer.TYPE_BYTE: 

			dst.setElem(row, col, band, (byte)value);

			if(value > Byte.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Byte.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Byte)byteNoData.getValue());break;

				}

			} 

			break;

		case Buffer.TYPE_SHORT: dst.setElem(row, col, band, value);break;

		case Buffer.TYPE_INT: dst.setElem(row, col, band, (int)value);break;

		case Buffer.TYPE_FLOAT: dst.setElem(row, col, band, (float)value);break;

		case Buffer.TYPE_DOUBLE: dst.setElem(row, col, band, (double)value);break;

		}

	}

	private void writePixelFromInt(int row, int col, int band, int value, Buffer dst) {

		switch (dst.getDataType()) {

		case Buffer.TYPE_BYTE: 

			dst.setElem(row, col, band, (byte)value);

			if(value > Byte.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Byte.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Byte)byteNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_SHORT:

			dst.setElem(row, col, band, (short)value);

			if(value > Short.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Short.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Short)shortNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_INT: dst.setElem(row, col, band, (int)value);break;

		case Buffer.TYPE_FLOAT: dst.setElem(row, col, band, (float)value);break;

		case Buffer.TYPE_DOUBLE: dst.setElem(row, col, band, (double)value);break;

		}

	}

	private void writePixelFromDouble(int row, int col, int band, double value, Buffer dst) {

		switch (dst.getDataType()) {

		case Buffer.TYPE_BYTE: 

			dst.setElem(row, col, band, (byte)value);

			if(value > Byte.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Byte.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Byte)byteNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_SHORT:

			dst.setElem(row, col, band, (short)value);

			if(value > Short.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Short.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Short)shortNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_INT:

			dst.setElem(row, col, band, (int)value);

			if(value > Integer.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Integer.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Integer)intNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_FLOAT:

			dst.setElem(row, col, band, (float)value);

			if(value > Float.MAX_VALUE) {

				switch (big2smal) {

				case 1:dst.setElem(row, col, band, Float.MAX_VALUE);break;

				case 2:dst.setElem(row, col, band, (Float)floatNoData.getValue());break;

				}

			}

			break;

		case Buffer.TYPE_DOUBLE: dst.setElem(row, col, band, (double)value);break;

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.raster.tools.app.basic.raster.process.RasterProcess#getResult()

	 */

	public Object getResult() {

		HashMap<String, Object> map = new HashMap<String, Object>();

		map.put(FILENAME, filename);

		map.put(TIME, new Long(milis));

		return map;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getPercent()

	 */

	public int getPercent() {

		return percent;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getTitle()

	 */

	public String getTitle() {

		return Messages.getText("layer_datatype");

	}

}

org.gvsig.raster.tools.algorithm.layerdatatype.LayerDatatypeAlgorithmLibrary

<?xml version="1.0" encoding="ISO-8859-1"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/maven-v4_0_0.xsd">

	<modelVersion>4.0.0</modelVersion>

	<artifactId>org.gvsig.raster.tools.algorithm.layerdatatype</artifactId>

	<packaging>jar</packaging>

	<name>org.gvsig.raster.tools.algorithm.layerdatatype</name>

	<parent>

		<groupId>org.gvsig</groupId>

		<artifactId>org.gvsig.raster.tools.algorithm</artifactId>

		<version>2.0.0</version>

	</parent>

    <properties>

        <build-dir>${basedir}/../../../build</build-dir>

    </properties>

    <dependencies>

		<dependency>

			<groupId>org.gvsig</groupId>

			<artifactId>org.gvsig.raster.tools.algorithm.base</artifactId>

            <scope>compile</scope>

		</dependency>

	</dependencies>

</project>

/*

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 gvSIG Association.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.raster.tools.algorithm.maskthreshold;

import org.gvsig.i18n.Messages;

import org.gvsig.raster.tools.algorithm.base.RasterBaseAlgorithmLibrary;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of MaskthresholdAlgorithmLibrary library.

 * 

 * @author <a href="mailto:nachobrodin@gmail.com">Nacho Brodin</a>

 */

public class MaskthresholdAlgorithmLibrary extends AbstractLibrary {

	public static final String         PROCESS_LABEL   = "MaskthresholdProcess";

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    public void doPostInitialize() throws LibraryException {

    	//Registers the process and its parameters

    	RasterBaseAlgorithmLibrary.register(PROCESS_LABEL, MaskthresholdProcess.class);

    	MaskthresholdProcess.registerParameters();

        Messages.addResourceFamily(

            "org.gvsig.raster.tools.algorithm.maskthreshold", 

            MaskthresholdAlgorithmLibrary.class.getClassLoader(), 

            MaskthresholdAlgorithmLibrary.class.getClass().getName());

        //registerGeoProcess(new RasterReprojectAlgorithmLibrary());

    }

}

/* gvSIG. Geographic Information System of the Valencian Government

*

* Copyright (C) 2007-2008 Infrastructures and Transports Department

* of the Valencian Government (CIT)

* 

* This program is free software; you can redistribute it and/or

* modify it under the terms of the GNU General Public License

* as published by the Free Software Foundation; either version 2

* of the License, or (at your option) any later version.

* 

* This program is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

* GNU General Public License for more details.

* 

* You should have received a copy of the GNU General Public License

* along with this program; if not, write to the Free Software

* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, 

* MA  02110-1301, USA.

* 

*/

package org.gvsig.raster.tools.algorithm.maskthreshold;

/**

 * This exception is thrown if happen problems processing data with this algorithm.

 * 

 * @author Nacho Brodin (nachobrodin@gmail.com)

 *

 */

public class MaskthresholdException extends Exception {

	private static final long serialVersionUID = -3022090543908771484L;

	public MaskthresholdException(String msg){

		super(msg);

	}

	public MaskthresholdException(String msg, Throwable e){

		super(msg, e);

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

*

* Copyright (C) 2007-2008 Infrastructures and Transports Department

* of the Valencian Government (CIT)

* 

* This program is free software; you can redistribute it and/or

* modify it under the terms of the GNU General Public License

* as published by the Free Software Foundation; either version 2

* of the License, or (at your option) any later version.

* 

* This program is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

* GNU General Public License for more details.

* 

* You should have received a copy of the GNU General Public License

* along with this program; if not, write to the Free Software

* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, 

* MA  02110-1301, USA.

* 

*/

package org.gvsig.raster.tools.algorithm.maskthreshold;

import java.awt.geom.Rectangle2D;

import java.util.HashMap;

import javax.swing.SwingUtilities;

import org.gvsig.fmap.dal.coverage.RasterLocator;

import org.gvsig.fmap.dal.coverage.dataset.Buffer;

import org.gvsig.fmap.dal.coverage.datastruct.NoData;

import org.gvsig.fmap.dal.coverage.exception.InvalidSetViewException;

import org.gvsig.fmap.dal.coverage.exception.ProcessInterruptedException;

import org.gvsig.fmap.dal.coverage.exception.RasterDriverException;

import org.gvsig.fmap.dal.coverage.store.RasterDataStore;

import org.gvsig.fmap.dal.coverage.store.RasterQuery;

import org.gvsig.i18n.Messages;

import org.gvsig.raster.tools.algorithm.base.RasterBaseAlgorithmLibrary;

import org.gvsig.raster.tools.algorithm.base.process.ProcessException;

import org.gvsig.raster.tools.algorithm.base.process.RasterProcess;

import org.gvsig.raster.tools.algorithm.base.util.Operation;

/**

 * Process to create a mask using a threshold. This threshold could be a value

 * of a pixel or a fixed value.

 * 

 * Possible use cases:

 * <UL>

 * <LI><b>Fixed-oplayer:</b>  A input layer, an operator layer and a fixed value if the condition is satisfied.

 * If the condition is not satisfied the value is set to NoData</LI>

 * 

 * <LI><b>Fixed-threshold:</b>  A input layer, a threshold and a fixed value if the condition is satisfied.

 * If the condition is not satisfied the value is set to NoData</LI>

 *

 * <LI><b>Replicate-oplayer:</b>  A input layer, an operator layer and a output layer. If the condition is satisfied.

 * the result layer replicates the output bands. If the condition is not satisfied the values is set to NoData</LI>

 * 

 * <LI><b>Replicate-threshold:</b>  A input layer, a threshold layer and a output layer. If the condition is satisfied.

 * the result layer replicates the output bands. If the condition is not satisfied the values is set to NoData</LI>

 * 

 * <LI><b>Transparent-oplayer:</b> A input layer, an operator layer and a RGB output layer. If the condition is satisfied

 * the layer result replicates the output RGB with transparency values. If the condition is not satisfied 

 * the layer result only replicates the output RGB.</LI>

 * 

 * <LI><b>Transparent-threshold:</b> A input layer, a threshold and a RGB output layer. If the condition is satisfied

 * the layer result replicates the output RGB with transparency values. If the condition is not satisfied 

 * the layer result only replicates the output RGB.</LI>

 * 

 * <LI><b>Opacity-oplayer:</b> A input layer, an operator layer, a RGB output layer. If the condition is satisfied

 * the layer result replicates the output layer values and if not this pixels are set to transparent </LI>

 * 

 * <LI><b>Opacity-threshold:</b> A input layer, an operator layer, a threshold. If the condition is satisfied

 * the layer result replicates the output layer values and if not this pixels are set to transparent </LI>

 * </UL>

 * 

 * <table>

 * <tr>

 * <th>band1</th><th>band2</th><th>threshold</th><th>output</th>

 * </tr>

 * <tr>

 * <td>-1</td><td>-1</td><td>-</td><td>el n?mero de bandas de la capa de operaci?n debe 

 *                                   ser mayor o igual a la de entrada. Compara la banda 

 *                                   i de la capa de entrada con la banda i de la de operaci?n</td>

 * </tr>

 * <tr>

 * <td>-1</td><td>0</td><td>-</td><td>Compara todas las bandas de la capa de entrada con la 

 *                                   banda cero de la capa de operaci?n</td>

 * </tr>

 * <tr>

 * <td>-1</td><td>-</td><td>value</td><td>Compara todas las bandas de la capa de entrada con el valor 

 *                                   de threshold</td>

 * </tr>

 * <tr>

 * <td>0</td><td>-1</td><td>-</td><td>Caso no permitido</td>

 * </tr>

 * <tr>

 * <td>0</td><td>-</td><td>value</td><td>Compara la banda cero de la capa de entrada con el valor

 *                                   de threshold</td>

 * </tr>                                  

 * </table>  

 * 

 * Input parameters

 * <UL>

 * <LI><b>RASTER_STORE1:</b>(Store method) Input store.</LI>

 * <LI><b>RASTER_STORE2:</b>(Store method) Operation store. If this parameter is null then THRESHOLD is mandatory</LI>

 * <LI><b>RASTER_STORE3:</b>(Store method) Output store. If this parameter is null then FIXED_VALUE is mandatory</LI>

 * <LI><b>RASTER_BUFFER1:</b>(Buffer method) Input buffer</LI>

 * <LI><b>RASTER_BUFFER2:</b>(Buffer method) Operation buffer. If this parameter is null then THRESHOLD is mandatory</LI>

 * <LI><b>RASTER_BUFFER3:</b>(Buffer method) Output buffer. If this parameter is null then FIXED_VALUE is mandatory</LI>

 * <LI><b>BAND1:</b> RASTER_STORE1 band. Band to compare with the threshold or operation band. If BAND1 is lesser than 0 then all

 * bands will be compared. If BAND1 is lesser than zero and BAND2 is lesser than zero then the store </LI>

 * <LI><b>BAND2:</b> RASTER_STORE2 band. Band to compare with the input store band</LI>

 * <LI><b>BAND3:</b> RASTER:STORE3 band. Band to put in the output when the method is "layer bands selecteds"</LI>

 * <LI><b>THRESHOLD:</b> Threshold to compare with the input data</LI>

 * <LI><b>FIXED_VALUE:</b> Value to put in the output when the method is "Fixed Value" or "NoData"</LI>

 * <LI><b>OPERATION:</b> Operation. Index in list MASK_OPERATIONS</LI>

 * <LI><b>SATISFY_METHOD:</b> Output method. Index in list OUTPUT_METHODS</LI>

 * <LI><b>PATH:</b></LI>

 * </UL>

 * 

 * Possible output parameters

 * <UL>

 * <LI><b>BUFFER:</b></LI>

 * <LI><b>FILENAME:</b></LI>

 * <LI><b>TIME:</b></LI>

 * </UL>

 * 10/12/2007

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public class MaskthresholdProcess extends RasterProcess {

	//Possible input parameters

	public static String      RASTER_STORE1    = "RasterStore1"; //Input store

	public static String      RASTER_STORE2    = "RasterStore2"; //Operation store if exists (threshold)

	public static String      RASTER_STORE3    = "RasterStore3"; //Output store if exists (fixed value)

	public static String      RASTER_BUFFER1   = "Buffer1";      //Intput buffer

	public static String      RASTER_BUFFER2   = "Buffer2";      //Operation buffer if exists (threshold)

	public static String      RASTER_BUFFER3   = "Buffer3";      //Output buffer if exists (fixed value)

	public static String      BAND1            = "Band1";       

	public static String      BAND2            = "Band2";

	public static String      BAND3            = "Band3";

	public static String      THRESHOLD        = "Threshold";

	public static String      FIXED_VALUE      = "FixedValue";

	public static String      OPERATION        = "Operation";

	public static String      SATISFY_METHOD   = "Satisfy_Method";

	public static String      PATH             = "Path";

	//Possible output parameters

	public static String      BUFFER           = "Buffer";

	public static String      FILENAME         = "FileName";

	public static String      TIME             = "Time";

	public static String[]    MASK_OPERATIONS  = new String[]{"Greater", "Lower", 

		                                                      "GreaterOrEquals", "LoweOrEquals"};

	public static String[]    OUTPUT_METHODS   = new String[]{"NoData", "FixedValue", 

		                                                      "Transparency", "Opacity", 

		                                                      "LayerAndBandsSelected", "LayerAndBandsSelectedInverse"};

	private RasterDataStore   inputStore       = null;

	private RasterDataStore   operationStore   = null;

	private RasterDataStore   outputStore      = null;

	private Buffer            inputBuffer      = null;

	private Buffer            operationBuffer  = null;

	private Buffer            outputBuffer     = null;

	private int               band1            = 0;

	private int               band2            = 0;

	private int               band3            = 0; //-1 All

	private int               resultNumBands   = 0;

	private String            filename         = null; //If is null only returns a buffer

	private double            threshold        = 0;

	private double            fixedValue       = 0;

	private NoData            nodata           = null;

	private AbstractOperation operation        = null;

	private AbstractOperation opSatisfy        = null;

	private AbstractOperation opNotSatisfy     = null;

	private long              milis            = 0;

	private int               percent          = 0;

	private Buffer            resultBuffer     = null;

	public static void registerParameters() {

		RASTER_STORE1 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_STORE1, RasterDataStore.class);

		RASTER_STORE2 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_STORE2, RasterDataStore.class);

		RASTER_STORE3 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_STORE3, RasterDataStore.class);

		RASTER_BUFFER1 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_BUFFER1, Buffer.class);

		RASTER_BUFFER2 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_BUFFER2, Buffer.class);

		RASTER_BUFFER3 = RasterBaseAlgorithmLibrary.registerInputParameter(RASTER_BUFFER3, Buffer.class);

		BAND1 = RasterBaseAlgorithmLibrary.registerInputParameter(BAND1, Integer.class);

		BAND2 = RasterBaseAlgorithmLibrary.registerInputParameter(BAND2, Integer.class);

		BAND3 = RasterBaseAlgorithmLibrary.registerInputParameter(BAND3, Integer.class);

		SATISFY_METHOD = RasterBaseAlgorithmLibrary.registerInputParameter(SATISFY_METHOD, Integer.class);

		THRESHOLD = RasterBaseAlgorithmLibrary.registerInputParameter(THRESHOLD, Double.class);

		FIXED_VALUE = RasterBaseAlgorithmLibrary.registerInputParameter(FIXED_VALUE, Double.class);

		OPERATION = RasterBaseAlgorithmLibrary.registerInputParameter(OPERATION, Integer.class);

		PATH = RasterBaseAlgorithmLibrary.registerInputParameter(PATH, String.class);

		FILENAME = RasterBaseAlgorithmLibrary.registerOutputParameter(FILENAME, String.class);

		BUFFER = RasterBaseAlgorithmLibrary.registerOutputParameter(BUFFER, Buffer.class);

		TIME = RasterBaseAlgorithmLibrary.registerOutputParameter(TIME, Long.class);

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.rastertools.RasterProcess#init()

	 */

	public void init() {

		inputStore = getParam(RASTER_STORE1) != null ? (RasterDataStore)getParam(RASTER_STORE1) : null;

		operationStore = getParam(RASTER_STORE2) != null ? (RasterDataStore)getParam(RASTER_STORE2) : null; //Default Null

		outputStore = getParam(RASTER_STORE3) != null ? (RasterDataStore)getParam(RASTER_STORE3) : null; //Default Null

		inputBuffer = getParam(RASTER_BUFFER1) != null ? (Buffer)getParam(RASTER_BUFFER1) : null;

		operationBuffer = getParam(RASTER_BUFFER2) != null ? (Buffer)getParam(RASTER_BUFFER2) : null;

		outputBuffer = getParam(RASTER_BUFFER3) != null ? (Buffer)getParam(RASTER_BUFFER3) : null;

		band1 = (Integer)getIntParam(BAND1); //Default 0

		band2 = (Integer)getIntParam(BAND2); //Default 0

		band3 = (Integer)getIntParam(BAND3); //Default 0

		//Caso no permitido. No se puede comparar una banda de la de entrada 

		//con todas de la capa de operaci?n. Evitamos el pete

		if(band2 < 0 && band1 >= 0) 

			band2 = 0;

		int satisfy = (Integer)getIntParam(SATISFY_METHOD); //Default 0

		if(satisfy == 0) {

			opSatisfy = new NoDataOp();

			opNotSatisfy = new FixedValueOp();

			resultNumBands = 1;

		}

		if(satisfy == 1) {

			opSatisfy = new FixedValueOp();

			opNotSatisfy = new NoDataOp();

			resultNumBands = 1;

		}

		if(satisfy == 2) {

			opSatisfy = new TransparencyOp();

			opNotSatisfy = new OpacityOp();

			resultNumBands = 4;

		}

		if(satisfy == 3) {

			opSatisfy = new OpacityOp();

			opNotSatisfy = new TransparencyOp();

			resultNumBands = 4;

		}

		if(satisfy == 4) {

			opSatisfy = new LayerOp();

			opNotSatisfy = new NoDataOp();

		}

		if(satisfy == 5) {

			opSatisfy = new NoDataOp();

			opNotSatisfy = new LayerOp();

		}

		if(satisfy == 4 || satisfy == 5) {

			if(outputStore != null)

				resultNumBands = band3 >= 0 ? band3 + 1 : outputStore.getBandCount();

			if(outputBuffer != null)

				resultNumBands = band3 >= 0 ? band3 + 1 : outputBuffer.getBandCount();

		}

		int op = (Integer)getIntParam(OPERATION); //Default 0

		if(op == 0)

			operation = new Greater();

		if(op == 1)

			operation = new Lower();

		if(op == 2)

			operation = new GreaterOrEquals();

		if(op == 3)

			operation = new LowerOrEquals();

		threshold = (Double)getDoubleParam(THRESHOLD); //Default 0

		fixedValue = (Double)getDoubleParam(FIXED_VALUE); //Default 0

		filename = getStringParam(PATH);

	}

	/**

	 * M?todo donde se ejecutar? el Thread, aqu? se reproyecta el raster.

	 */

	public void process() throws ProcessInterruptedException, ProcessException {

		long t1 = new java.util.Date().getTime();

		insertLineLog(Messages.getText("maskthreshold"));

		try {

			if(inputStore != null)

				checkStore();

			else

				checkBuffer();

			Rectangle2D[] pxBBox = null;

			if(inputStore != null) {

				inputBuffer = readROBuffer(inputStore, band1);

				RasterDataStore[] storeList = new RasterDataStore[]{inputStore, operationStore, outputStore};

				pxBBox = getIntersectionInPxCoords(storeList);

			} else {

				Rectangle2D r = new Rectangle2D.Double(0, 0, inputBuffer.getWidth(), inputBuffer.getHeight());

				pxBBox = new Rectangle2D[]{r, r, r};

			}

			if(operationStore != null) 

				operationBuffer = readROBuffer(operationStore, band2);

			int dataType = Buffer.TYPE_DOUBLE;

			if(outputStore != null) { 

				outputBuffer = readROBuffer(outputStore, band3);

			}

			if(outputBuffer != null) {

				dataType = outputBuffer.getDataType();

			}

			resultBuffer = RasterLocator.getManager().createBuffer(

					dataType, 

					(int)pxBBox[0].getWidth(), 

					(int)pxBBox[0].getHeight(), 

					resultNumBands, true);

			write(resultBuffer, pxBBox);

			if(filename != null)

				exportRaster(filename, resultBuffer, inputStore.getCellSize(), 

						inputStore.getExtent().getULX(), inputStore.getExtent().getULY());

			long t2 = new java.util.Date().getTime();

			milis = t2 - t1;

			SwingUtilities.invokeLater(new Runnable() {

				public void run() {

					if (externalActions != null) {

						HashMap<String, Object> map = new HashMap<String, Object>();

						map.put(FILENAME, filename);

						map.put(BUFFER, resultBuffer);

						map.put(TIME, new Long(milis));

						externalActions.end(map);

					}

				}

			});

		} catch (MaskthresholdException e) {

			if (incrementableTask != null)

				incrementableTask.processFinalize();

			messageBoxError("error_processing_maskthreshold", null, e);

		} finally {

			if (incrementableTask != null) {

				incrementableTask.processFinalize();

				setProgressActive(false);

			}

		}

	}

	/**

	 * Checks input parameters for DataStores

	 * @throws MaskthresholdException

	 */

	private void checkStore() throws MaskthresholdException {

		if(opSatisfy instanceof TransparencyOp || opSatisfy instanceof OpacityOp && 

				(outputStore == null || outputStore.getBandCount() < 3 || outputStore.getDataType()[0] != Buffer.TYPE_BYTE)) {

			throw new MaskthresholdException("Not valid parameters. Transparency method have 3 output bands");

		}

		if((opSatisfy instanceof LayerOp || opNotSatisfy instanceof LayerOp) && 

				outputStore == null)

			throw new MaskthresholdException("Not valid parameters. Layer method have output data store");

		if (inputStore == null)

			throw new MaskthresholdException("Layer not valid");

		if(operationStore != null) {

			if(!inputStore.getExtent().intersects(operationStore.getExtent()))

				throw new MaskthresholdException("Extents don't intersect");

			if(inputStore.getCellSize() != operationStore.getCellSize())

				throw new MaskthresholdException("Cell sizes are not equal");

			if(band1 < 0 && band2 < 0) { // Si se comparan todas las bandas, la capa dos debe tener m?s bandas que la uno o las mismas

				if(operationStore.getBandCount() < inputStore.getBandCount())

					throw new MaskthresholdException("If exists operation layer and all bands have to be compared then the operation layer must have more or the same number of bands that the input layer.");

			}

		}

	}

	/**

	 * Checks input parameters for Buffers

	 * @throws MaskthresholdException

	 */

	private void checkBuffer() throws MaskthresholdException {

		if(opSatisfy instanceof TransparencyOp || opSatisfy instanceof OpacityOp && 

				(outputBuffer == null || outputBuffer.getBandCount() < 3 || outputBuffer.getDataType() != Buffer.TYPE_BYTE)) {

			throw new MaskthresholdException("Not valid parameters. Transparency method have 3 output bands");

		}

		if((opSatisfy instanceof LayerOp || opNotSatisfy instanceof LayerOp) && 

				outputBuffer == null)

			throw new MaskthresholdException("Not valid parameters. Layer method have output data store");

		if (inputBuffer == null)

			throw new MaskthresholdException("Layer not valid");

		if(outputBuffer != null) {

			if ((inputBuffer.getWidth() != outputBuffer.getWidth()) || 

					(inputBuffer.getHeight() != outputBuffer.getHeight()))

				throw new MaskthresholdException("Buffer sizes not valid");

		}

		if(operationBuffer != null) {

			if(band1 < 0 && band2 < 0) { // Si se comparan todas las bandas, la capa dos debe tener m?s bandas que la uno o las mismas

				if(operationBuffer.getBandCount() < inputBuffer.getBandCount())

					throw new MaskthresholdException("If exists operation layer and all bands have to be compared then the operation layer must have more or the same number of bands that the input layer.");

			}

		}

	}

	/**

	 * Reads a band of a buffer

	 * @param store

	 * @param band

	 * @return

	 * @throws MaskthresholdException

	 * @throws ProcessInterruptedException

	 */

	private Buffer readROBuffer(RasterDataStore store, int band) throws MaskthresholdException, ProcessInterruptedException {

		RasterQuery query = RasterLocator.getManager().createQuery();

		if(band >= 0 && band < store.getBandCount())

			query.setDrawableBands(new int[]{band});

		else

			query.setAllDrawableBands();

		query.setReadOnly(true);

		query.setAreaOfInterest();

		try {

			return store.query(query);

		} catch (RasterDriverException e) {

			throw new MaskthresholdException("Error reading input raster.");

		} catch (InvalidSetViewException e) {

			throw new MaskthresholdException("Error reading input raster.");

		}

	}

	private void write(Buffer result, Rectangle2D[] pxBbox) {

		if(nodata == null)

			nodata = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(1, result.getDataType());

		int colSource = (int)pxBbox[0].getX();

		int rowSource = (int)pxBbox[0].getY();

		int colOperation = colSource;

		int rowOperation = rowSource;

		if(pxBbox.length > 1) {

			colOperation = (int)pxBbox[1].getX();

			rowOperation = (int)pxBbox[1].getY();

		}

		opSatisfy.fixedValue = opNotSatisfy.fixedValue = fixedValue;

		opSatisfy.result = opNotSatisfy.result = result;

		opSatisfy.outBuffer = opNotSatisfy.outBuffer = outputBuffer;

		opSatisfy.band = opNotSatisfy.band = band3;

		int inputInitBand = band1 < 0 ? 0 : band1;

		int inputEndBand = band1 < 0 ? inputBuffer.getBandCount() : band1 + 1;

		for (int row = 0; row < result.getHeight(); row++) {

			if(pxBbox.length > 1) {

				colOperation = (int)pxBbox[1].getX();

				colSource = (int)pxBbox[0].getX();

			}

			colSource = 0;

			for (int col = 0; col < result.getWidth(); col++) {

				boolean compare = true;

				for (int iBands = inputInitBand; iBands < inputEndBand; iBands++) {

					double op1 = getValue(inputBuffer, rowSource, colSource, iBands);

					double op2 = getValue(operationBuffer, rowOperation, colOperation, getBandToCompare(operationBuffer, iBands));

					if(!operation.compare(op1, op2)) {

						compare = false;

						break;

					}

				}

				if(compare) {

					opSatisfy.row = row;

					opSatisfy.col = col;

					opSatisfy.invoke();

				} else {

					opNotSatisfy.row = row;

					opNotSatisfy.col = col;

					opNotSatisfy.invoke();

				}

				colSource ++;

				colOperation++;

			}

			rowSource ++;

			rowOperation ++;

			percent = (int)((row * 100) / result.getHeight());

		}

	}

	/**

	 * Gets the band to operate. If the buffer is null then it does not 

	 * matter and it returns zero.

	 * @param buf

	 * @param inputBand

	 * @return

	 */

	private int getBandToCompare(Buffer buf, int inputBand) {

		if(buf == null)

			return 0;

		if(band2 < 0)

			return inputBand;

		else

			return band2;

	}

	/**

	 * Gets the operation value in using a double value. If there is a operation data store then

	 * returns the value of the pixel in this position else it returns a threshold value

	 * @param operationBuffer

	 * @param rowOperation

	 * @param colOperation

	 * @return

	 */

	private double getValue(Buffer buf, int rowOperation, int colOperation, int band) {

		if(buf == null)

			return threshold;

		if(buf.getDataType() == Buffer.TYPE_BYTE)

			return buf.getElemByte(rowOperation, colOperation, band);

		if(buf.getDataType() == Buffer.TYPE_FLOAT)

			return buf.getElemFloat(rowOperation, colOperation, band);

		if(buf.getDataType() == Buffer.TYPE_DOUBLE)

			return buf.getElemDouble(rowOperation, colOperation, band);

		if(buf.getDataType() == Buffer.TYPE_SHORT)

			return buf.getElemShort(rowOperation, colOperation, band);

		if(buf.getDataType() == Buffer.TYPE_INT)

			return buf.getElemInt(rowOperation, colOperation, band);

		return 0;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.raster.tools.app.basic.raster.process.RasterProcess#getResult()

	 */

	public Object getResult() {

		HashMap<String, Object> map = new HashMap<String, Object>();

		map.put(FILENAME, filename);

		map.put(BUFFER, resultBuffer);

		map.put(TIME, new Long(milis));

		return map;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getPercent()

	 */

	public int getPercent() {

		return percent;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getTitle()

	 */

	public String getTitle() {

		return Messages.getText("maskthreshold");

	}

	abstract class AbstractOperation implements Operation {

		public int       row          = 0;

		public int       col          = 0;

		public Buffer    result       = null;

		public double    fixedValue   = 0;

		public Buffer    outBuffer   = null;

		public int       band         = 0;

		public void invoke() {

		}

		public boolean compare(double a, double b) {

			return false;

		}

	}

	class LayerOp extends AbstractOperation {

		public void invoke() {

			if(result.getDataType() == Buffer.TYPE_BYTE) {

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						byte vSource = outBuffer.getElemByte(row, col, i);

						result.setElem(row, col, i, vSource);

					}

				} else {

					byte vSource = outBuffer.getElemByte(row, col, band);

					result.setElem(row, col, band, vSource);

				}

			} else if(result.getDataType() == Buffer.TYPE_FLOAT) {

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						float vSource = outBuffer.getElemFloat(row, col, i);

						result.setElem(row, col, i, vSource);

					}

				} else {

					float vSource = outBuffer.getElemFloat(row, col, band);

					result.setElem(row, col, band, vSource);

				}

			} else if(result.getDataType() == Buffer.TYPE_DOUBLE) {

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						double vSource = outBuffer.getElemDouble(row, col, i);

						result.setElem(row, col, i, vSource);

					}

				} else {

					double vSource = outBuffer.getElemDouble(row, col, band);

					result.setElem(row, col, band, vSource);

				}

			} else if(result.getDataType() == Buffer.TYPE_SHORT) {

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						float vSource = outBuffer.getElemShort(row, col, i);

						result.setElem(row, col, i, vSource);

					}

				} else {

					float vSource = outBuffer.getElemShort(row, col, band);

					result.setElem(row, col, band, vSource);

				}

			} else if(result.getDataType() == Buffer.TYPE_INT) {

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						int vSource = outBuffer.getElemInt(row, col, i);

						result.setElem(row, col, i, vSource);

					}

				} else {

					int vSource = outBuffer.getElemInt(row, col, band);

					result.setElem(row, col, band, vSource);

				}

			}

		}

	}

	class TransparencyOp extends AbstractOperation {

		public void invoke() {

			if(result.getDataType() == Buffer.TYPE_BYTE) {

				for (int i = 0; i < 3; i++) {

					byte vSource = outBuffer.getElemByte(row, col, i);

					result.setElem(row, col, i, vSource);

				}

				result.setElem(row, col, 4, (byte)0);

			}

		}

	}

	class OpacityOp extends AbstractOperation {

		public void invoke() {

			if(result.getDataType() == Buffer.TYPE_BYTE) {

				for (int i = 0; i < 3; i++) {

					byte vSource = outBuffer.getElemByte(row, col, i);

					result.setElem(row, col, i, vSource);

				}

				result.setElem(row, col, 4, (byte)255);

			}

		}

	}

	class FixedValueOp extends AbstractOperation {

		public void invoke() {

			if(result.getDataType() == Buffer.TYPE_BYTE)

				result.setElem(row, col, 0, (byte)fixedValue);

			else if(result.getDataType() == Buffer.TYPE_FLOAT)

				result.setElem(row, col, 0, (float)fixedValue);

			else if(result.getDataType() == Buffer.TYPE_DOUBLE)

				result.setElem(row, col, 0, fixedValue);

			else if(result.getDataType() == Buffer.TYPE_SHORT)

				result.setElem(row, col, 0, (short)fixedValue);

			else if(result.getDataType() == Buffer.TYPE_INT)

				result.setElem(row, col, 0, (int)fixedValue);

		}

	}

	class NoDataOp extends AbstractOperation {

		public void invoke() {

			if(result.getDataType() == Buffer.TYPE_BYTE)

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						result.setElem(row, col, i, nodata.getValue().byteValue());

					}

				} else {

					result.setElem(row, col, band, nodata.getValue().byteValue());

				}

			else if(result.getDataType() == Buffer.TYPE_FLOAT)

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						result.setElem(row, col, i, nodata.getValue().floatValue());

					}

				} else {

					result.setElem(row, col, band, nodata.getValue().floatValue());

				}

			else if(result.getDataType() == Buffer.TYPE_DOUBLE)

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						result.setElem(row, col, i, nodata.getValue().doubleValue());

					}

				} else {

					result.setElem(row, col, band, nodata.getValue().doubleValue());

				}

			else if(result.getDataType() == Buffer.TYPE_SHORT)

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						result.setElem(row, col, i, nodata.getValue().shortValue());

					}

				} else {

					result.setElem(row, col, band, nodata.getValue().shortValue());

				}

			else if(result.getDataType() == Buffer.TYPE_INT)

				if(band < 0) {

					for (int i = 0; i < result.getBandCount(); i++) {

						result.setElem(row, col, i, nodata.getValue().intValue());

					}

				} else {

					result.setElem(row, col, band, nodata.getValue().intValue());

				}

		}

	}

	class Greater extends AbstractOperation {

		public boolean compare(double a, double b) {

			return a > b;

		}

	}

	class Lower extends AbstractOperation {

		public boolean compare(double a, double b) {

			return a < b;

		}

	}

	class GreaterOrEquals extends AbstractOperation {

		public boolean compare(double a, double b) {

			return a >= b;

		}

	}

	class LowerOrEquals extends AbstractOperation {

		public boolean compare(double a, double b) {

			return a <= b;

		}

	}

}

org.gvsig.raster.tools.algorithm.maskthreshold.MaskthresholdAlgorithmLibrary

<?xml version="1.0" encoding="ISO-8859-1"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/maven-v4_0_0.xsd">

	<modelVersion>4.0.0</modelVersion>

	<artifactId>org.gvsig.raster.tools.algorithm.maskthreshold</artifactId>

	<packaging>jar</packaging>

	<name>org.gvsig.raster.tools.algorithm.maskthreshold</name>

	<parent>

		<groupId>org.gvsig</groupId>

		<artifactId>org.gvsig.raster.tools.algorithm</artifactId>

		<version>2.0.0</version>

	</parent>

    <properties>

        <build-dir>${basedir}/../../../build</build-dir>

    </properties>

    <dependencies>

		<dependency>

			<groupId>org.gvsig</groupId>

			<artifactId>org.gvsig.raster.tools.algorithm.base</artifactId>

            <scope>compile</scope>

		</dependency>

	</dependencies>

</project>

org.gvsig.raster.tools.algorithm.base.RasterBaseAlgorithmLibrary

/*

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2012 gvSIG Association.

 *

 * This program is free software; you can redistribute it and/or

 * modify it under the terms of the GNU General Public License

 * as published by the Free Software Foundation; either version 2

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License