Application: gvSIG desktop » gvSIG raster

org.gvsig.raster.algorithm.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

 * 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.algorithm;

import org.gvsig.i18n.Messages;

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.extensionpoint.ExtensionPoint;

import org.gvsig.tools.extensionpoint.ExtensionPointManager;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of RasterBaseAlgorithmLibrary library.

 * 

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

 */

public class RasterBaseAlgorithmLibrary extends AbstractLibrary {

	public static final String         REGISTER_PROCESS_LABEL             = "RasterProcess";

	public static final String         REGISTER_INPUT_PARAMETERS_LABEL    = "RasterProcessInputParameters";

	public static final String         REGISTER_OUTPUT_PARAMETERS_LABEL   = "RasterProcessOutputParameters";

	/**

	 * Registers a raster process

	 */

	@SuppressWarnings("unchecked")

	public static void register(String processLabel, Class process) {

		ExtensionPointManager extensionPoints = ToolsLocator.getExtensionPointManager();

		ExtensionPoint point = extensionPoints.add(REGISTER_PROCESS_LABEL);

		point.append(processLabel, "", process);

	}

	/**

	 *  Registers input parameters of a raster process

	 */

	@SuppressWarnings("unchecked")

	public static String registerInputParameter(String parameterLabel, Class parameterClass) {

		return registerParameter(parameterLabel, parameterClass, REGISTER_INPUT_PARAMETERS_LABEL);

	}

	/**

	 *  Registers output parameters of a raster process

	 */

	@SuppressWarnings("unchecked")

	public static String registerOutputParameter(String parameterLabel, Class parameterClass) {

		return registerParameter(parameterLabel, parameterClass, REGISTER_OUTPUT_PARAMETERS_LABEL);

	}

	/**

	 * Registers a raster process

	 */

	@SuppressWarnings("unchecked")

	private static String registerParameter(String parameterLabel, Class parameterClass, String type) {

		ExtensionPointManager extensionPoints = ToolsLocator.getExtensionPointManager();

		ExtensionPoint point = extensionPoints.add(type);

		int i = 0;

		//Si no hay un par?metro registrado con ese nombre se registra

		if(point.get(parameterLabel) == null) {

			point.append(parameterLabel, "", parameterClass);

			return parameterLabel;

		} else {

			//Si hay un par?metro registrado con ese nombre pero los tipos son iguales no se hace nada

			if(point.get(parameterLabel).getExtension() == parameterClass)

				return parameterLabel;

			//Buscamos un nombre no registrado para el par?metro a?adiendo sufijos

			String newName = parameterLabel + "_" + i;

			while (point.get(newName) != null && point.get(newName).getExtension() != parameterClass) {

				i++;

				newName = parameterLabel + "_" + i;

			}

			point.append(parameterLabel + "_" + i, "", parameterClass);

		}

		return parameterLabel + "_" + i;

	}

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

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

            RasterBaseAlgorithmLibrary.class.getClassLoader(), 

            RasterBaseAlgorithmLibrary.class.getClass().getName());

        //registerGeoProcess(new RasterReprojectAlgorithmLibrary());

    }

    /**

     * Returns the instance of the manager

     * @return

     */

    public static RasterBaseAlgorithmManager getManager() {

    	return RasterBaseAlgorithmManager.getInstance();

    }

}

/* 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.algorithm;

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

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

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

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.extensionpoint.ExtensionPoint;

import org.gvsig.tools.extensionpoint.ExtensionPointManager;

/**

 * Default {@link RasterManager} implementation.

 * 

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

 * @version $Id$

 */

public class RasterBaseAlgorithmManager {

	private static RasterBaseAlgorithmManager internalInstance  = new RasterBaseAlgorithmManager();

	/**

	 * Gets an instance of this object for internal use.

	 * @return DefaultRasterManager

	 */

	public static RasterBaseAlgorithmManager getInstance() {

		return internalInstance;

	}

	/**

	 * Returns a new instance of a process.

	 * @param processLabel

	 *        Label to identify the process

	 * @return

	 */

	public RasterProcess createRasterTask(String processLabel) throws ProcessException {

		ExtensionPointManager extensionPoints = ToolsLocator.getExtensionPointManager();

		ExtensionPoint point = extensionPoints.add(RasterBaseAlgorithmLibrary.REGISTER_PROCESS_LABEL);

		try {

			return (RasterProcess)point.create(processLabel);

		} catch (InstantiationException e) {

			throw new ProcessException("", e);

		} catch (IllegalAccessException e) {

			throw new ProcessException("", 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.algorithm.util;

/**

 * Interface for operations

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public interface Operation {

	public boolean compare(double a, double b);

	public void invoke();

}

/* 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.algorithm.util;

/**

 * Interface for operations

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public abstract class AbstractOperation implements Operation {

	public boolean compare(double a, double b) {

		return false;

	}

	public void invoke() {

		//do nothing

	}

}

/* 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.algorithm.util;

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

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

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

/**

 * Calculates a pixel value using a interpolation method

 * @author Nacho Brodin nachobrodin@gmail.com

 * @author Victor Olaya

 */

public class Interpolation {

	private Buffer        buffer    = null;

	private double        nodata    = 0;

	public Interpolation(Buffer buf) {

		this.buffer = buf;

		NoData nodata = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(

				1, Buffer.TYPE_DOUBLE);

		this.nodata = nodata.getValue().doubleValue();

	}

	private double[] getKernel(int x, int y, int band) {

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

			return getKernelByte(x, y, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_DOUBLE) {

			return getKernelByte(x, y, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_FLOAT) {

			return getKernelByte(x, y, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_SHORT) {

			return getKernelByte(x, y, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_INT) {

			return getKernelByte(x, y, band);

		}

		return null;

	}

	public double getNearestNeighbour(double x, double y, int band) {

		int dy = (int)Math.round(y);

		int dx = (int)Math.round(x);

		dy = dy < buffer.getHeight() ? dy : buffer.getHeight() - 1;

		dx = dx < buffer.getWidth() ? dx : buffer.getWidth() - 1;

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

			return (double)buffer.getElemByte(dy, dx, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_DOUBLE) {

			return (double)buffer.getElemDouble(dy, dx, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_FLOAT) {

			return (double)buffer.getElemFloat(dy, dx, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_SHORT) {

			return (double)buffer.getElemShort(dy, dx, band);

		}

		if(buffer.getDataType() == Buffer.TYPE_INT) {

			return (double)buffer.getElemInt(dy, dx, band);

		}

		return nodata;

	}

	/**

	 * Calcula los valores N y Z para el m?todo bilinear y obtiene el valor del pixel como

	 * Z / N

	 * @param dx distancia en X desde el centro del pixel hasta el punto. Es un valor entre 0 y 1

	 * @param dy distancia en Y desde el centro del pixel hasta el punto. Es un valor entre 0 y 1

	 * @param kernel valor del pixel y alrededor 

	 * @return valor del pixel

	 */

	public double getBilinearValue(double x, double y, int band) {

		double[] kernel = getKernel((int)x, (int)y, band);

		double dx = x - ((int) x);

		double dy = y - ((int) y);

		double z = 0.0, n = 0.0, d;

		d = (1.0 - dx) * (1.0 - dy);

		z += d * kernel[0];

		n += d;

		d = dx * (1.0 - dy);

		z += d * kernel[1]; 

		n += d;

		d = (1.0 - dx) * dy;

		z += d * kernel[2]; 

		n += d;

		d = dx * dy;

		z += d * kernel[3]; 

		n += d;

		double b = 0;

		if(n > 0.0)

			b = (z / n);

		return b;

	}

	/**

	 * Calcula los valores N y Z para el m?todo de distancia inversa y calcula el valor del

	 * pixel como Z / N.

	 * @param dx distancia en X desde el centro del pixel hasta el punto. Es un valor entre 0 y 1

	 * @param dy distancia en Y desde el centro del pixel hasta el punto. Es un valor entre 0 y 1

	 * @param kernel valor del pixel y alrededor 

	 * @return valor del pixel

	 */

	public double getInverseDistance(double x, double y, int band) {

		double[] kernel = getKernel((int)x, (int)y, band);

		double dx = x - ((int) x);

		double dy = y - ((int) y);

		double z = 0.0, n = 0.0, d;

		d = 1.0 / Math.sqrt(dx * dx + dy * dy);

		z += d * kernel[0];

		n += d;

		d = 1.0 / Math.sqrt((1.0 - dx) * ( 1.0 - dx) + dy * dy);

		z += d * kernel[1]; 

		n += d;

		d = 1.0 / Math.sqrt(dx*dx + (1.0-dy)*(1.0-dy));

		z += d * kernel[2]; 

		n += d;

		d = 1.0 / Math.sqrt((1.0 - dx) *( 1.0 - dx) + (1.0 - dy) * (1.0 - dy));

		z += d * kernel[3]; 

		n += d;

		double b = 0;

		if(n > 0.0)

			b = (z / n);

		return b;

	}

	/**

	 * Obtiene un kernel de cuatro elemento que corresponden a los pixeles (x, y), (x + 1, y),

	 * (x, y + 1), (x + 1, y + 1). Si los pixeles x + 1 o y + 1 se salen del raster de origen

	 * se tomar? x e y. 

	 * @param x Coordenada X del pixel inicial

	 * @param y Coordenada Y del pixel inicial

	 * @param band N?mero de banda.

	 * @return Kernel solicitado en forma de array.

	 */

	private double[] getKernelByte(int x, int y, int band) {

		double[] d = new double[4];

		d[0] = (buffer.getElemByte(y, x, band) & 0xff);

		int nextX = ((x + 1) >= buffer.getWidth()) ? x : (x + 1);

		int nextY = ((y + 1) >= buffer.getHeight()) ? y : (y + 1);

		d[1] = (buffer.getElemByte(y, nextX, band) & 0xff);

		d[2] = (buffer.getElemByte(nextY, x, band) & 0xff);

		d[3] = (buffer.getElemByte(nextY, nextX, band) & 0xff);

		return d;

	}

	/**

	 * Obtiene un kernel de cuatro elemento que corresponden a los pixeles (x, y), (x + 1, y),

	 * (x, y + 1), (x + 1, y + 1). Si los pixeles x + 1 o y + 1 se salen del raster de origen

	 * se tomar? x e y. 

	 * @param x Coordenada X del pixel inicial

	 * @param y Coordenada Y del pixel inicial

	 * @param band N?mero de banda.

	 * @return Kernel solicitado en forma de array.

	 */

	@SuppressWarnings("unused")

	private double[] getKernelShort(int x, int y, int band) {

		double[] d = new double[4];

		d[0] = (buffer.getElemShort(y, x, band) & 0xffff);

		int nextX = ((x + 1) >= buffer.getWidth()) ? x : (x + 1);

		int nextY = ((y + 1) >= buffer.getHeight()) ? y : (y + 1);

		d[1] = (buffer.getElemShort(y, nextX, band) & 0xffff);

		d[2] = (buffer.getElemShort(nextY, x, band) & 0xffff);

		d[3] = (buffer.getElemShort(nextY, nextX, band) & 0xffff);

		return d;

	}

	/**

	 * Obtiene un kernel de cuatro elemento que corresponden a los pixeles (x, y), (x + 1, y),

	 * (x, y + 1), (x + 1, y + 1). Si los pixeles x + 1 o y + 1 se salen del raster de origen

	 * se tomar? x e y. 

	 * @param x Coordenada X del pixel inicial

	 * @param y Coordenada Y del pixel inicial

	 * @param band N?mero de banda.

	 * @return Kernel solicitado en forma de array.

	 */

	@SuppressWarnings("unused")

	private double[] getKernelInt(int x, int y, int band) {

		double[] d = new double[4];

		d[0] = (buffer.getElemInt(y, x, band) & 0xffffffff);

		int nextX = ((x + 1) >= buffer.getWidth()) ? x : (x + 1);

		int nextY = ((y + 1) >= buffer.getHeight()) ? y : (y + 1);

		d[1] = (buffer.getElemInt(y, nextX, band) & 0xffffffff);

		d[2] = (buffer.getElemInt(nextY, x, band) & 0xffffffff);

		d[3] = (buffer.getElemInt(nextY, nextX, band) & 0xffffffff);

		return d;

	}

	/**

	 * Obtiene un kernel de cuatro elemento que corresponden a los pixeles (x, y), (x + 1, y),

	 * (x, y + 1), (x + 1, y + 1). Si los pixeles x + 1 o y + 1 se salen del raster de origen

	 * se tomar? x e y. 

	 * @param x Coordenada X del pixel inicial

	 * @param y Coordenada Y del pixel inicial

	 * @param band N?mero de banda.

	 * @return Kernel solicitado en forma de array.

	 */

	@SuppressWarnings("unused")

	private double[] getKernelFloat(int x, int y, int band) {

		double[] d = new double[4];

		d[0] = buffer.getElemFloat(y, x, band);

		int nextX = ((x + 1) >= buffer.getWidth()) ? x : (x + 1);

		int nextY = ((y + 1) >= buffer.getHeight()) ? y : (y + 1);

		d[1] = buffer.getElemFloat(y, nextX, band);

		d[2] = buffer.getElemFloat(nextY, x, band);

		d[3] = buffer.getElemFloat(nextY, nextX, band);

		return d;

	}

	/**

	 * Obtiene un kernel de cuatro elemento que corresponden a los pixeles (x, y), (x + 1, y),

	 * (x, y + 1), (x + 1, y + 1). Si los pixeles x + 1 o y + 1 se salen del raster de origen

	 * se tomar? x e y. 

	 * @param x Coordenada X del pixel inicial

	 * @param y Coordenada Y del pixel inicial

	 * @param band N?mero de banda.

	 * @return Kernel solicitado en forma de array.

	 */

	@SuppressWarnings("unused")

	private double[] getKernelDouble(int x, int y, int band) {

		double[] d = new double[4];

		d[0] = buffer.getElemDouble(y, x, band);

		int nextX = ((x + 1) >= buffer.getWidth()) ? x : (x + 1);

		int nextY = ((y + 1) >= buffer.getHeight()) ? y : (y + 1);

		d[1] = buffer.getElemDouble(y, nextX, band);

		d[2] = buffer.getElemDouble(nextY, x, band);

		d[3] = buffer.getElemDouble(nextY, nextX, band);

		return d;

	}

}

/* 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.algorithm.process;

/**

 * Acciones que puede ejecutar un proceso sobre un objeto externo. Un proceso

 * debe ser independiente de objetos que sean dependientes de su funcionalidad. Es

 * decir, cualquier extensi?n debe poder usar un proceso y este ejecutar metodos

 * de la funcionalidad. Esto se hace ha trav?s de este interfaz.

 * 

 * @version 26/09/2007

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

 *

 */

public interface IProcessActions {

	/**

	 * El proceso comunica que ha sido interrumpido

	 *

	 */

	public void interrupted();

	/**

	 * Acciones de finalizaci?n del proceso

	 */

	public void end(Object param);

}

/* 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.algorithm.process;

/**

 * This exception is thrown when a problem inside a raster process

 * is detected. 

 * 

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

 */

public class ProcessException extends Exception {

	private static final long serialVersionUID = 666908550965442025L;

	public ProcessException(String msg, Throwable e) {

		super(msg, e);

	}

	public ProcessException(String msg){

		super(msg);

	}

}

/* 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.algorithm.process;

import java.awt.Component;

import java.awt.geom.AffineTransform;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import java.util.ArrayList;

import java.util.Hashtable;

import javax.swing.JOptionPane;

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

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

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

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

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

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

import org.gvsig.fmap.dal.coverage.process.CancelEvent;

import org.gvsig.fmap.dal.coverage.process.TaskEventManager;

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

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

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

import org.gvsig.fmap.dal.coverage.store.props.ColorInterpretation;

import org.gvsig.gui.beans.buttonspanel.ButtonsPanel;

import org.gvsig.gui.beans.incrementabletask.IIncrementable;

import org.gvsig.gui.beans.incrementabletask.IncrementableEvent;

import org.gvsig.gui.beans.incrementabletask.IncrementableListener;

import org.gvsig.gui.beans.incrementabletask.IncrementableTask;

import org.gvsig.i18n.Messages;

import org.gvsig.tools.dispose.Disposable;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

/**

 * Clase base de todos los procesos raster. En ella se genstionan todas las

 * funciones comunes como incremento de la tarea, gesti?n de eventos a la tarea, 

 * par?metros de la tarea, etc ...

 * 

 * 18/12/2007

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public abstract class RasterProcess implements IIncrementable, IncrementableListener, Runnable, Disposable {

	protected IncrementableTask incrementableTask = null;

	protected volatile Thread   blinker           = null;

	protected TaskEventManager  taskEventManager  = null;

	protected IProcessActions   externalActions   = null;

	protected Hashtable<String, Object>

	                            taskParams        = new Hashtable<String, Object>();

	private String              log               = "";

	private String              lastLine          = "";

	private long                time              = 0;

	private boolean             progressActive    = true;

	private IProcessActions     queueActions      = null;

	private Logger              logger            = LoggerFactory.getLogger(RasterProcess.class.toString());

	protected int               percent           = 0;

	/**

	 * Crea la ventana de IncrementableTask

	 */

	public IncrementableTask getIncrementableTask() {

		if (incrementableTask == null) {

			incrementableTask = new IncrementableTask(this);

			incrementableTask.addIncrementableListener(this);

		}

		return incrementableTask;

	}

	/**

	 * Define si se puede cancelar el proceso. Por defecto es true.

	 * @param enabled

	 */

	public void setCancelable(boolean enabled) {

		getIncrementableTask().getButtonsPanel().setEnabled(ButtonsPanel.BUTTON_CANCEL, enabled);

	}

	/**

	 * Muestra la ventana de IncrementableTask

	 */

	public void showIncrementableWindow() {

		if (progressActive) {

			getIncrementableTask().showWindow();

			getIncrementableTask().start();

		}

	}

	/**

	 * Arranca el proceso de recorte de un layer

	 */

	public void start() {

		showIncrementableWindow();

		if(blinker == null)

			blinker = new Thread(this);

		blinker.start();

	}

	/**

	 * Proceso de carga de par?metros. Puede no hacerse una carga de par?metros 

	 * explicita y obtenerlos directamente de la tabla Hash cuando vayan a usarse. 

	 * Esto queda a elecci?n del programador. En caso de hacerse una carga de par?metros

	 * explicita esta llamada deber?a hacerse justo despues del constructor de la clase

	 * que contendr? el proceso.

	 */

	public abstract void init();

	/**

	 * Proceso

	 * @throws InterruptedException

	 */

	public abstract void process() throws ProcessInterruptedException, ProcessException;

	/**

	 * Obtenci?n de un objeto de resultado. Este objeto deber? ser el mismo que el que se

	 * pasa por par?metro en el "end" de IProcessActions. Este m?todo es util cuando la tarea no se

	 * lanza en un Thread (ejecutamos process directamente) y queremos recoger el resultado

	 * sin registrarnos al evento de finalizaci?n de tarea.

	 * @return objeto resultado de la tarea

	 */

	public Object getResult() {

		return null;

	}

	/**

	 * Proceso

	 * @throws ProcessException 

	 */

	public void execute() throws ProcessInterruptedException, ProcessException {

		init();

		process();

	}

	/**

	 * M?todo donde se ejecutar? el Thread. Este har? las acciones globales para 

	 * cualquier tarea y llamar? al m?todo execute especifico de una tarea.

	 */

	public void run() {

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

		try {

			taskEventManager = RasterLocator.getManager().createRasterTask(this);

			execute();

		} catch (ProcessInterruptedException e) {

			if (externalActions != null)

				externalActions.interrupted();

			Thread.currentThread().interrupt();

		} catch (ProcessException e) {

			if (progressActive) {

				if (incrementableTask != null) {

					getIncrementableTask().processFinalize();

					incrementableTask = null;

				}

			}

			logger.warn(RasterLocator.getManager().getRasterUtils().getTrace(e));

			if(e.getMessage() != null && e.getMessage().compareTo("") != 0)

				messageBoxError(e.getMessage(), this);

			else

				messageBoxError("error_processing", this);

			queueActions = null;

		} catch (Exception e) {

			if (progressActive) {

				if (incrementableTask != null) {

					getIncrementableTask().processFinalize();

					incrementableTask = null;

				}

			}

			logger.warn(RasterLocator.getManager().getRasterUtils().getTrace(e));

			if(e.getMessage() != null && e.getMessage().compareTo("") != 0)

				messageBoxError(e.getMessage(), null);

			else

				messageBoxError("error_processing", null);

			queueActions = null;

		} finally {

			taskEventManager.removeTask();

			if (progressActive) {

				if (incrementableTask != null)

					getIncrementableTask().processFinalize();

			}

			if (queueActions != null)

				queueActions.end(this);

			time = new java.util.Date().getTime() - t1;

			blinker = null;

		}

	}

	/**

	 * Activa o desactiva el interfaz de progreso en el lanzamiento de la tarea como un thread.

	 * @param active true para activarlo o false para desactivarlo

	 */

	public void setProgressActive(boolean active) {

		this.progressActive = active;

	}

	/**

	 * Obtiene el tiempo que tard? en ejecutarse la tarea 

	 * la ?ltima vez que se proces?

	 */

	public long getTime() {

		return time;

	}

	/**

	 * Obtiene el objeto para ejecutar acciones externar.

	 * @param IProcessActions

	 */

	public IProcessActions getActions() {

		return externalActions;

	}

	/**

	 * Asigna el objeto para ejecutar acciones externar.

	 * @param IProcessActions

	 */

	public void setActions(IProcessActions actions) {

		this.externalActions = actions;

	}

	/**

	 * Obtiene el objeto para ejecutar acciones de la cola de procesos de ejecuci?n exclusiva.

	 * @param IProcessActions

	 */

	public IProcessActions getUniqueProcessActions() {

		return queueActions;

	}

	/**

	 * Asigna el objeto para ejecutar acciones externar.

	 * @param IProcessActions

	 */

	public void setUniqueProcessActions(IProcessActions actions) {

		this.queueActions = actions;

	}

	/**

	 * Inserta una nueva l?nea en el log del cuadro de incremento de tarea

	 * @param line

	 */

	protected void insertLineLog(String line) {

		lastLine = line;

		log = log + line + "\n";

	}

	/**

	 * Obtiene la ?ltima l?nea introducida en el log del cuadro de incremento.

	 */

	public String getLabel() {

		return lastLine;

	}

	/**

	 * Obtiene el texto de log del cuadro de incremento completo.

	 */

	public String getLog() {

		return log;

	}

	/**

	 * Un evento de cancelado es enviado a la tarea cuando actionCanceled es activado. Para

	 * ello se crear? un objeto CancelEvent y se asignar? a la tarea en ejecuci?n. Esta lo

	 * procesar? cuando pueda e interrumpir? el proceso.

	 */

	public void actionCanceled(IncrementableEvent e) {

		taskEventManager.setEvent(new CancelEvent(this));

	}

	/**

	 * A?ade un par?metro a la tarea

	 * @param name Clave del par?metro

	 * @param param Objeto pasado como par?metro

	 */

	public void addParam(String name, Object param) {

		if (param != null)

			taskParams.put(name, param);

		else

			taskParams.remove(name);

	}

	/**

	 * Elimina un par?metro de la tarea

	 * @param name Clave del par?metro a eliminar

	 */

	public void removeParam(String name) {

		taskParams.remove(name);

	}

	/**

	 * Obtiene un par?metro a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public Object getParam(String name) {

		return taskParams.get(name);

	}

	/**

	 * Obtiene un par?metro String a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public String getStringParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof String) ? (String)value : null;

	}

	/**

	 * Obtiene un par?metro byte a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public byte getByteParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Byte) ? ((Byte)value).byteValue() : 0;

	}

	/**

	 * Obtiene un par?metro float a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public float getFloatParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Float) ? ((Float)value).floatValue() : 0F;

	}

	/**

	 * Obtiene un par?metro double a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public double getDoubleParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Double) ? ((Double)value).doubleValue() : 0D;

	}

	/**

	 * Obtiene un par?metro entero a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public int getIntParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Integer) ? ((Integer)value).intValue() : 0;

	}

	/**

	 * Obtiene un par?metro boolean a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public boolean getBooleanParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Boolean) ? ((Boolean)value).booleanValue() : false;

	}

	/**

	 * Obtiene un par?metro int[] a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public int[] getIntArrayParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof int[]) ? ((int[])value) : null;

	}

	/**

	 * Obtiene un par?metro double[] a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public double[] getDoubleArrayParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof double[]) ? ((double[])value) : null;

	}

	/**

	 * Obtiene un par?metro extent a partir de la clave

	 * @param name Par?metro

	 * @return Par?metro

	 */

	public Extent getExtentParam(String name) {

		Object value = taskParams.get(name);

		return (value != null && value instanceof Extent) ? ((Extent)value) : null;

	}

	/**

	 * Updates the percentaje of progress and manages the cancelation

	 * @param parcial

	 * @param total

	 * @throws ProcessInterruptedException

	 */

	public void updatePercent(int parcial, int total) throws ProcessInterruptedException {

		if (taskEventManager != null && taskEventManager.getEvent() != null)

			taskEventManager.manageEvent(taskEventManager.getEvent());

		percent = (int)((parcial * 100) / total);

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IncrementableListener#actionResumed(org.gvsig.gui.beans.incrementabletask.IncrementableEvent)

	 */

	public void actionResumed(IncrementableEvent e) {	

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.gui.beans.incrementabletask.IncrementableListener#actionSuspended(org.gvsig.gui.beans.incrementabletask.IncrementableEvent)

	 */

	public void actionSuspended(IncrementableEvent e) {

	}

	/*

	 * (non-Javadoc)

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

	 */

	public boolean isCancelable() {

		return true;

	}

	/*

	 * (non-Javadoc)

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

	 */

	public boolean isPausable() {

		return false;

	}

	/**

	 * Registra un mensaje de error en el log de gvSIG

	 * @param msg Mensaje a guardar en el log

	 * @param parent Objeto que hizo disparar el mensaje

	 * @param exception Excepcion que ha sido recogida

	 */

	public static void debug(String msg, Object parent, Exception exception) {

		if(parent != null)

		    LoggerFactory

            .getLogger(parent.getClass()).debug(Messages.getText(msg), exception);

	}

	public void dispose() {

		try {

			finalize();

		} catch (Throwable e) {

		}

	}

	/*

	 * (non-Javadoc)

	 * @see java.lang.Object#finalize()

	 */

	@Override

	protected void finalize() throws Throwable {

		incrementableTask           = null;

		blinker                     = null;

		externalActions             = null;

		if(taskParams != null) {

			taskParams.clear();

			taskParams = null;

		}

		log                         = null;

		lastLine                    = null;

		queueActions                = null;

		logger                      = null;

		super.finalize();

	}

	/**

	 * Shows a error dialog with a text and a accept button 

	 * @param msg Message to show in the dialog

	 * @param parentWindow Parent window

	 */

	public static void messageBoxError(String msg, Object parentWindow){

		String string = Messages.getText("accept");

		Object[] options = {string};

		JOptionPane.showOptionDialog((Component)parentWindow,

					"<html>" + Messages.getText(msg).replaceAll("\n", "<br>") + "</html>",

					Messages.getText("confirmacion"),

					JOptionPane.OK_OPTION,

					JOptionPane.ERROR_MESSAGE,

					null,

					options,

					string);

	}

	/**

	 * Muestra un dialogo de error con un texto y un bot?n de aceptar. El error es

	 * registrado en el log de gvSIG con la excepcion que se le pase por parametro

	 * @param msg Mensaje a mostrar en el dialogo.

	 * @param parentWindow Ventana desde la que se lanza el dialogo

	 * @param exception Excepcion que ha sido recogida

	 */

	public static void messageBoxError(String msg, Object parentWindow, Exception exception) {

		debug(msg, parentWindow, exception);

		messageBoxError(msg, parentWindow);

	}

	/**

	 * Muestra un dialogo de error con un texto y un bot?n de aceptar. Se le pasa como ?ltimo par?metros

	 * una lista de excepciones que ser?n guardadas en el log

	 * @param msg Mensaje a mostrar en el dialogo.

	 * @param parentWindow Ventana desde la que se lanza el dialogo

	 * @param exception Excepcion que ha sido recogida

	 */

	public static void messageBoxError(String msg, Object parentWindow, ArrayList<Exception> exception) {

		for (int i = 0; i < exception.size(); i++) 

			debug(msg, parentWindow, exception.get(i));

		messageBoxError(msg, parentWindow);

	}

	/**

	 * Exports a raster buffer to disk

	 * @param sFilename

	 * @param buf

	 * @param cellsize

	 * @param minX

	 * @param minY

	 * @return

	 */

	public boolean exportRaster(final String sFilename, 

			Buffer buf, 

			Buffer alphaBuffer, 

			double cellsize, 

			double minX, 

			double minY) {

		try {

			RasterManager manager = RasterLocator.getManager();

			final DataServerWriter writerBufferServer = manager.createDataServerWriter();

			writerBufferServer.setBuffer(buf, -1);

			int nBands = buf.getBandCount();

			ColorInterpretation colorInterpretation = null;

			if(alphaBuffer != null) {

				writerBufferServer.setAlphaBuffer(alphaBuffer);

				nBands ++;

				colorInterpretation = RasterLocator.getManager().getDataStructFactory().createColorInterpretation(

						new String[] { 

								ColorInterpretation.RED_BAND, 

								ColorInterpretation.GREEN_BAND, 

								ColorInterpretation.BLUE_BAND,

								ColorInterpretation.ALPHA_BAND});

			} else {

				if(nBands == 1)

					colorInterpretation = RasterLocator.getManager().getDataStructFactory().createColorInterpretation(

							new String[] { ColorInterpretation.GRAY_BAND });

			}

			final Params params = manager.createWriterParams(sFilename);

			final AffineTransform affineTransform =

				new AffineTransform(cellsize, 0, 0,

						-cellsize, minX, minY);

			final RasterWriter writer = manager.createWriter(

						writerBufferServer, 

						sFilename,

						nBands, 

						affineTransform, 

						buf.getWidth(),

						buf.getHeight(), 

						buf.getDataType(), 

						params, 

						null);

			if(colorInterpretation != null)

				writer.setColorBandsInterpretation(colorInterpretation.getValues());

			writer.dataWrite();

			writer.writeClose();

		} catch (final Exception e) {

			e.printStackTrace();

			return false;

		}

		return true;

	}

	/**

	 * Exports a raster buffer to disk

	 * @param sFilename

	 * @param buf

	 * @param cellsize

	 * @param minX

	 * @param minY

	 * @return

	 */

	public boolean exportRaster(final String sFilename, 

			Buffer buf, 

			double cellsize, 

			double minX, 

			double minY) {

        return exportRaster(sFilename, buf, null, cellsize, minX, minY);

    }

	/**

	 * Gets a list of rectangles which represents the pixel coordinates of each DataStore.

	 * This rectangle is the area that intersects with the other DataStores in the 	list.

	 * @param dataStoreList

	 * @return

	 */

	protected Rectangle2D[] getIntersectionInPxCoords(RasterDataStore[] dataStoreList) {

		Extent extentIntersect = null;

		int nRectangles = 0;

		for (int i = 0; i < dataStoreList.length; i++) {

			if(dataStoreList[i] != null) {

				if(extentIntersect == null)

					extentIntersect = dataStoreList[i].getExtent();

				else

					extentIntersect = extentIntersect.intersection(dataStoreList[i].getExtent());

				nRectangles ++;

			}

		}

		Rectangle2D[] result = new Rectangle2D[nRectangles];

		Point2D p1 = new Point2D.Double(extentIntersect.getULX(), extentIntersect.getULY());

		Point2D p2 = new Point2D.Double(extentIntersect.getLRX(), extentIntersect.getLRY());

		int cont = 0;

		for (int i = 0; i < dataStoreList.length; i++) {

			if(dataStoreList[i] != null) {

				Point2D init = dataStoreList[i].worldToRaster(p1);

				Point2D end = dataStoreList[i].worldToRaster(p2);

				result[cont] = new Rectangle2D.Double(

						init.getX(), 

						init.getY(), 

						Math.abs(end.getX() - init.getX()), 

						Math.abs(end.getY() - init.getY()));

				cont++;

			}

		}

		return result;

	}

}

<?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.algorithm</artifactId>

	<packaging>jar</packaging>

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

	<parent>

		<groupId>org.gvsig</groupId>

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

		<version>2.1.0-SNAPSHOT</version>

	</parent>

	<dependencies>

		<dependency>

			<groupId>org.gvsig</groupId>

			<artifactId>org.gvsig.raster.lib.api</artifactId>

            <scope>compile</scope>

		</dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.raster.lib.impl</artifactId>

            <scope>compile</scope>

        </dependency>

     	 <dependency>

          <groupId>org.gvsig</groupId>

          <artifactId>org.gvsig.metadata.lib.basic.impl</artifactId>

          <scope>runtime</scope>

      	</dependency>

      	<dependency>

          <groupId>org.gvsig</groupId>

          <artifactId>org.gvsig.metadata.swing.basic.impl</artifactId>

          <scope>runtime</scope>

      	</dependency>

        <dependency>

			<groupId>org.gvsig</groupId>

			<artifactId>org.gvsig.raster.cache.lib.api</artifactId>

            <scope>compile</scope>

		</dependency>

		<dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.tools.lib</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.i18n</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.ui</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.projection.api</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.projection.cresques.impl</artifactId>

            <scope>runtime</scope>

        </dependency>

         <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.compat.api</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.compat.se</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.dal.api</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.dal.impl</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.dal.file.lib</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.dal.spi</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.geometry.api</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.geometry.impl</artifactId>

            <scope>runtime</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.fmap.geometry.operation</artifactId>

            <scope>runtime</scope>

        </dependency>

	</dependencies>

    <properties>

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

    </properties>

</project>

org.gvsig.fmap.dal.coverage.RasterLibrary

/* 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