Application: gvSIG desktop » gvSIG raster

<?xml version="1.0" encoding="UTF-8"?>

<projectDescription>

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

	<comment></comment>

	<projects>

	</projects>

	<buildSpec>

		<buildCommand>

			<name>org.eclipse.m2e.core.maven2Builder</name>

			<arguments>

			</arguments>

		</buildCommand>

	</buildSpec>

	<natures>

		<nature>org.eclipse.m2e.core.maven2Nature</nature>

	</natures>

</projectDescription>

multispectral=Raster multiespectral

tasseledcab=Tasseled Cab

use_roi=Usar las regiones de inter?s de la capa de entrada

multispectral=Multispectral raster

tasseledcab=Tasseled Cab

use_roi=Use regions of interest of the input layer

multispectral=Raster multiespectral

tasseledcab=Tasseled Cab

tasseledcap=Tasseled Cab

use_roi=Usar las regiones de inter?s

stadistics=Estad?sticas

calculating=Calculando

multispectral=Multispectral raster

tasseledcab=Tasseled Cab

tasseledcap=Tasseled Cab

use_roi=Use regions of interest

stadistics=Stadistics

calculating=Calculating

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

<help>

	<element name="DESCRIPTION"

		text=""

		description="Descripción" type="0">

		<image description="" file="">

		</image>

	</element>

	<element name="ADDITIONAL_INFO" text=""

		description="Información adicional" type="0">

	</element>

	<element name="EXTENSION_AUTHOR" text=""

		description="Algoritmo creado por" type="0">

	</element>

	<element name="HELP_AUTHOR" text="" description="Ayuda creada por"

		type="0">

	</element>

	<element name="USER_NOTES" text="" description="Notas de usuario"

		type="0">

	</element>

</help>

<?xml version='1.0' encoding='ISO-8859-1' standalone='yes' ?>

<help>

	<element name="DESCRIPTION"

		text=""

		description="Descripción" type="0">

		<image description="" file="">

		</image>

	</element>

	<element name="ADDITIONAL_INFO" text=""

		description="Información adicional" type="0">

	</element>

	<element name="EXTENSION_AUTHOR" text=""

		description="Algoritmo creado por" type="0">

	</element>

	<element name="HELP_AUTHOR" text="" description="Ayuda creada por"

		type="0">

	</element>

	<element name="USER_NOTES" text="" description="Notas de usuario"

		type="0">

	</element>

</help>

org.gvsig.raster.tasseledcap.TasseledCapSextanteLibrary

package org.gvsig.raster.tasseledcap;

import java.util.List;

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

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

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

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

import org.gvsig.geoprocess.lib.sextante.AbstractSextanteGeoProcess;

import org.gvsig.geoprocess.lib.sextante.dataObjects.FLyrRasterIRasterLayer;

import org.gvsig.raster.algorithm.RasterBaseAlgorithmLibrary;

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

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

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

import org.gvsig.raster.fmap.layers.FLyrRaster;

import org.gvsig.raster.tasseledcap.algorithm.TasseledCapAlgorithmLibrary;

import org.gvsig.raster.tasseledcap.algorithm.TasseledCapProcess;

import es.unex.sextante.core.AnalysisExtent;

import es.unex.sextante.core.Sextante;

import es.unex.sextante.dataObjects.IRasterLayer;

import es.unex.sextante.exceptions.GeoAlgorithmExecutionException;

import es.unex.sextante.exceptions.RepeatedParameterNameException;

/**

 * ...

 */

public class TasseledCapSextanteAlgorithm extends AbstractSextanteGeoProcess implements IProcessActions {

    public static final String RESULT            = "RESULT";

    public static final String LAYER             = TasseledCapProcess.RASTER_STORE;

    public static final String USE_ROI           = "USE_ROI";

    public static final String TYPE              = "TYPE";

    public static String[]     LANDSAT_TYPES     = new String[]{"Landsat MS" , "Landsat TM",  "Landsat ETM"};

    private boolean            useROI            = false;

    private DataProcess        taskStats         = null;

    private DataProcess        taskPC            = null;

    private int                type              = 0;

    public void defineCharacteristics() {

        setName(getTranslation("tasseled_cap"));

        setGroup(getTranslation("multispectral"));

        try {

            m_Parameters.addInputRasterLayer(LAYER, getTranslation("Input_layer"), true);

            m_Parameters.addBoolean(USE_ROI, getTranslation("use_roi"), false);

            m_Parameters.addSelection(TYPE, getTranslation("type"), LANDSAT_TYPES);

        } catch (RepeatedParameterNameException e) {

            Sextante.addErrorToLog(e);

        }

        addOutputRasterLayer(RESULT, getTranslation("tasseled_cap"));

    }

    public boolean processAlgorithm() throws GeoAlgorithmExecutionException {

    	if(existsOutPutFile(TasseledCapSextanteAlgorithm.RESULT, 0)) {

    		throw new GeoAlgorithmExecutionException(getTranslation("file_exists"));

    	}

    	IRasterLayer input = m_Parameters.getParameterValueAsRasterLayer(LAYER);

		useROI = m_Parameters.getParameterValueAsBoolean(USE_ROI);

		String value = m_Parameters.getParameterValueAsString(TYPE);

		type = 0;

		for (int j = 0; j < LANDSAT_TYPES.length; j++) {

			if(LANDSAT_TYPES[j].equals(value))

				type = j;

		}

    	FLyrRaster lyrRaster = ((FLyrRaster)input.getBaseDataObject());

    	IRasterLayer output = null;

    	output = getNewRORasterLayer(

    			RESULT,

    			Sextante.getText("tasseledcap_description"),

    			input.getDataType(),

    			input.getBandsCount());

    	String fileName = ((FLyrRasterIRasterLayer)output).getName();

    	try {

			setProgressText(getTranslation("calculating"));

			taskPC = createTasseledCapProcess(lyrRaster.getDataStore(), fileName, type);

			taskPC.execute();

			((FLyrRasterIRasterLayer)output).setBaseDataObject(fileName);

		} catch (ProcessInterruptedException e) {

			Sextante.addErrorToLog(e);

		} catch (ProcessException e) {

			Sextante.addErrorToLog(e);

		}

		if(getTaskMonitor().isCanceled())

			return false;

        return true;

    }

    private DataProcess createTasseledCapProcess(RasterDataStore inputStore, String fileName, int type) throws ProcessException {

    	DataProcess taskPC = RasterBaseAlgorithmLibrary.getManager().createRasterTask(TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

    	taskPC.setActions(this);

    	List<String> params = taskPC.getRasterTaskInputParameters(TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

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

    		String paramName = params.get(i);

    		Class<?> paramType = taskPC.getParameterTypeByProcess(TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL, paramName);

    		if(paramType == RasterDataStore.class) {

    			taskPC.addParam(paramName, (RasterDataStore)inputStore);

    		}

			if(paramType == Boolean[].class) {

    			boolean[] bands = new boolean[type == 0 ? 4 : 6];

    			for (int j = 0; j < bands.length; j++) {

					bands[j] = true;

				}

    			taskPC.addParam(paramName, bands);

			}

			if(paramType == Integer.class) {

				taskPC.addParam(paramName, type);

			}

			if(paramName.equals("PATH")) {

				taskPC.addParam(paramName, fileName);

			}

			if(paramName.equals("ROI_EPSG") && useROI) {

				taskPC.addParam(paramName, "EPSG:4326");

    		}

			if(paramName.equals("WINDOW")) {

    			AnalysisExtent ext = getAnalysisExtent();

    			Extent bbox = RasterLocator.getManager().getDataStructFactory().createExtent(

    					ext.getXMin(), ext.getYMax(), ext.getXMax(), ext.getYMin());

    			Extent inputBbox = inputStore.getExtent();

    			if(bbox.getULX() != inputBbox.getULX() ||

    				bbox.getULY() != inputBbox.getULY() ||

    				bbox.getLRX() != inputBbox.getLRX() ||

    				bbox.getLRY() != inputBbox.getLRY()) {

    				taskPC.addParam(paramName, bbox);

    			}

    		}

    	}

    	return taskPC;

    }

	public void interrupted() {

	}

	public void end(Object param) {

	}

	public void updateProgress(int current, int total) {

		boolean cancelled = setProgress(current, total);

		if(!cancelled) {

			if(taskStats != null)

				taskStats.actionCanceled(null);

			if(taskPC != null)

				taskPC.actionCanceled(null);

		}

	}

    /*

     * TODO: Customized panels

    @Override

    public Class<? extends GeoAlgorithmParametersPanel> getCustomParametersPanelClass() {

        return PrincipalComponentsParametersPanel.class;

    }*/

}

package org.gvsig.raster.tasseledcap;

import org.gvsig.geoprocess.algorithm.base.core.AlgorithmAbstractLibrary;

import org.gvsig.i18n.Messages;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of <code>PrincipalComponentsSextanteLibrary</code> library.

 */

public class TasseledCapSextanteLibrary extends AlgorithmAbstractLibrary {

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

        Messages.addResourceFamily(

            "org.gvsig.raster.tasseledcap.sextante.i18n.text", TasseledCapSextanteLibrary.class

                .getClassLoader(), TasseledCapSextanteLibrary.class.getClass().getName());

        registerGeoProcess(new TasseledCapSextanteAlgorithm());

    }

}

<?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/xsd/maven-4.0.0.xsd">

	<modelVersion>4.0.0</modelVersion>

	<artifactId>org.gvsig.raster.tasseledcap.toolbox.algorithm</artifactId>

	<packaging>jar</packaging>

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

	<parent>

		<groupId>org.gvsig</groupId>

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

		<version>2.2.246</version>

	</parent>

	<dependencies>

		<dependency>

		    <groupId>org.gvsig</groupId>

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

            <scope>compile</scope>

   		</dependency>

		<dependency>

		    <groupId>org.gvsig</groupId>

   			<artifactId>org.gvsig.geoprocess.algorithm.base</artifactId>

            <scope>compile</scope>

   		</dependency>

		<dependency>

			<groupId>org.gvsig</groupId>

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

			<scope>compile</scope>

		</dependency>

		<dependency>

			<groupId>org.gvsig</groupId>

			<artifactId>org.gvsig.geoprocess.lib.sextante</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.ui</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.fmap.mapcontext.api</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

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

            <scope>runtime</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.metadata.lib.basic.api</artifactId>

            <scope>compile</scope>

        </dependency>

	</dependencies>

</project>

<?xml version="1.0" encoding="UTF-8"?>

<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/xsd/maven-4.0.0.xsd">

	<modelVersion>4.0.0</modelVersion>

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

	<packaging>pom</packaging>

	<version>2.2.246</version>

	<name>${project.artifactId}</name>

	<description />

	<inceptionYear>2011</inceptionYear>

	<parent>

      <groupId>org.gvsig</groupId>

      <artifactId>org.gvsig.desktop</artifactId>

      <version>2.0.445</version>

  </parent>

        <properties>

            <!-- El plugin versions:use-latest-versions falla con scope import -->

            <!-- asi que toca usar el versions:update-properties que si que funciona -->

            <org.gvsig.raster.version>2.2.265</org.gvsig.raster.version>

            <org.gvsig.geoprocess.version>2.2.272</org.gvsig.geoprocess.version>

        </properties>

    <repositories>

      <repository>

        <id>gvsig-public-http-repository</id>

        <name>gvSIG maven public HTTP repository</name>

        <url>http://devel.gvsig.org/m2repo/j2se</url>

        <releases>

          <enabled>true</enabled>

          <updatePolicy>daily</updatePolicy>

          <checksumPolicy>warn</checksumPolicy>

        </releases>

        <snapshots>

          <enabled>true</enabled>

          <updatePolicy>daily</updatePolicy>

          <checksumPolicy>warn</checksumPolicy>

        </snapshots>

      </repository>

    </repositories>

    <scm>

        <connection>scm:svn:https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.tasseledcap/tags/org.gvsig.raster.tasseledcap-2.2.246</connection>

        <developerConnection>scm:svn:https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.tasseledcap/tags/org.gvsig.raster.tasseledcap-2.2.246</developerConnection>

        <url>https://devel.gvsig.org/redmine/projects/gvsig-raster/repository/show/org.gvsig.raster.tasseledcap/tags/org.gvsig.raster.tasseledcap-2.2.246</url>

    </scm>

    <build>

        <plugins>

            <plugin>

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-release-plugin</artifactId>

                <configuration>

                    <tagBase>https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.tasseledcap/tags/</tagBase>

                    <goals>deploy</goals>

                </configuration>

            </plugin>

        </plugins>

    </build>

	<dependencyManagement>

		<dependencies>

                        <dependency>

                                <groupId>org.gvsig</groupId>

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

                                <version>${org.gvsig.raster.version}</version>

                                <type>pom</type>

                                <scope>import</scope>

                        </dependency>

                        <dependency>

                                <groupId>org.gvsig</groupId>

                                <artifactId>org.gvsig.geoprocess</artifactId>

                                <version>${org.gvsig.geoprocess.version}</version>

                                <type>pom</type>

                                <scope>import</scope>

                        </dependency>

			<dependency>

				<groupId>org.gvsig.legacy</groupId>

				<artifactId>jama</artifactId>

				<version>1.0.3</version>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

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

				<version>2.2.1-SNAPSHOT</version>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

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

				<version>2.2.1-SNAPSHOT</version>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster.tasseledcap.swing.api</artifactId>

				<version>2.2.246</version>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster.tasseledcap.swing.impl</artifactId>

				<version>2.2.246</version>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

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

				<version>2.2.246</version>

			</dependency>

   			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster.tasseledcap.toolbox.algorithm</artifactId>

				<version>2.2.246</version>

			</dependency>

		</dependencies>

	</dependencyManagement>

	<modules>

		<module>org.gvsig.raster.tasseledcap.algorithm</module>

		<module>org.gvsig.raster.tasseledcap.toolbox.algorithm</module>

		<module>org.gvsig.raster.tasseledcap.swing</module>

		<module>org.gvsig.raster.tasseledcap.app.tasseledcapclient</module>

	</modules>

</project>

org.gvsig.raster.tasseledcap.algorithm.TasseledCapAlgorithmLibrary

layers_dont_intersect=Las capas no intersectan

roi_not_inside=La zona de estudio no esta contenida en la extensi?n de las capas

preparing_buffers=Preparando buffers

writting_in_buffer=Escribiendo en buffer de salida

processing_tc=Procesando Tasseled Cap

tasseled_cap=Tasseled Cap

type=Tipo

layers_dont_intersect=Layers do not intersect

roi_not_inside=The region of interest are not inside the extension of the layers

preparing_buffers=Preparing buffers

writting_in_buffer=Writing in the output buffer

processing_tc=Procesando Tasseled Cap

tasseled_cap=Tasseled Cap

type=Type

package org.gvsig.raster.tasseledcap.algorithm;

import org.gvsig.i18n.Messages;

import org.gvsig.raster.algorithm.RasterBaseAlgorithmLibrary;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.LibraryException;

/**

 * Initialization of PrincipalComponentsAlgorithmLibrary library.

 */

public class TasseledCapAlgorithmLibrary extends AbstractLibrary {

	public static final String         TASSELEDCAP_PROCESS_LABEL         = "TasseledCapProcess";

    @Override

    protected void doInitialize() throws LibraryException {

        // Nothing to do

    }

    @Override

    protected void doPostInitialize() throws LibraryException {

    	//Registers the process and its parameters

    	RasterBaseAlgorithmLibrary.register(TASSELEDCAP_PROCESS_LABEL, TasseledCapProcess.class);

    	TasseledCapProcess.registerParameters();

        Messages.addResourceFamily(

            "org.gvsig.raster.tasseledcap.algorithm", 

            TasseledCapAlgorithmLibrary.class.getClassLoader(), 

            TasseledCapAlgorithmLibrary.class.getClass().getName());

        //registerGeoProcess(new RasterReprojectAlgorithmLibrary());

        Messages.addResourceFamily("org.gvsig.raster.tasseledcap.algorithm.i18n.text",

        		TasseledCapAlgorithmLibrary.class.getClassLoader(),

        		TasseledCapAlgorithmLibrary.class.getClass().getName()); 

    }

}

package org.gvsig.raster.tasseledcap.algorithm;

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

/**

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

 */

public class TasseledCapException extends ProcessException {

	private static final long serialVersionUID = -3022090543908771484L;

	public TasseledCapException(String msg){

		super(msg);

	}

	public TasseledCapException(String msg, Throwable e){

		super(msg, e);

	}

}

package org.gvsig.raster.tasseledcap.algorithm;

import java.awt.geom.Rectangle2D;

import java.util.List;

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

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

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

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

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

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

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

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

import org.gvsig.i18n.Messages;

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

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

import org.gvsig.raster.roi.ROI;

/**

 * Tasseled Cab Process 

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

 */

public class TasseledCapProcess extends DataProcess {

	public static int             LANDSAT_MS         = 0;

	public static int             LANDSAT_TM         = 1;

	public static int             LANDSAT_ETM        = 2;

	public static String          RASTER_STORE       = "RASTER_STORE";

	public static String          PATH               = "PATH";

	public static String          FILENAME           = "FILENAME";

	public static String          LANDSAT_TYPE       = "LANDSAT_TYPE";

	public static final String    BUFFERS            = "BUFFERS";

	public static final String    GRAPHIC_DATA       = "GRAPHIC_DATA";

	public static final String    BANDS              = "BANDS";

	private RasterDataStore       store              = null;

	private String                filename           = null;

	private List<ROI>             rois               = null;

	private Extent                extentResult       = null;

	private boolean[]             bands              = null;

	private NoData                nodata             = null;

	private double                matrixParams[][]   = null;

	private int                   type               = LANDSAT_MS;

	// MATRICES DE COEFICIENTES PARA CADA TRANSFORMACI?N

	double LandSatMS[][] = {

			/*   Brillo de suelo       */ 	{0.433, 0.632, 0.586, 0.264   },

			/*	 Indice de Verdor      */	{-0.290, -0.562, 0.600, 0.491 },

			/*	 Indice de Senescencia */   {-0.829,  0.522, -0.039, 0.194},

			/*	 Cuarta                */	{0.223,  0.012,  -0.543, 0.810}		

	};

	double LandSatTM[][] = { 

			/*   Brillo    */ 	{0.33183,  0.33121,  0.55177,  0.42514,  0.48087, 0.25252  		},

			/*	 Verdor    */	{-0.24717,  -0.16263,  -0.40639,  0.85468,  0.05493,  -0.11749	},

			/*	 Tercera   */	{0.13929,  0.22490,  0.40359,  0.25178,  -0.70133, -0.45732  	},	

	};

	double LandSatETM[][] = {

			/*   Brightness */ 	{0.3561,  0.3972,  0.3904,  0.6966,  0.2286,  0.1596     },

			/*	 Greenness 	*/	{-0.3344,  -0.3544,  -0.4556,  0.6966,  -0.0242,  -0.2630},

			/*	 Vetness 	*/	{0.2626,  0.2141,  0.0926,  0.0656,  -0.7629, -0.5388    },	

			/*	 Fourth 	*/	{0.0805,  -0.0498,  0.1950, -0.1327, 0.5752,  -0.7775    },

			/*	 Fifth 		*/	{-0.7252,  -0.0202,  0.6683,  0.0631,  -0.1494,  -0.0274 },

			/*	 Sixth 		*/	{0.4000,  -0.8172,  0.3832,  0.0602,  -0.1095,  -0.0985	 }

	};

	public static void registerParameters() {

		registerInputParameter(RASTER_STORE, RasterDataStore.class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

		registerInputParameter(PATH, String.class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

		registerInputParameter(BANDS, Boolean[].class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

		registerInputParameter(LANDSAT_TYPE, Integer.class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

		registerOutputParameter(FILENAME, String.class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

		registerOutputParameter(GRAPHIC_DATA, Double[].class, TasseledCapAlgorithmLibrary.TASSELEDCAP_PROCESS_LABEL);

	}

	public void init() {

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

		filename = getStringParam(PATH);

		bands = (boolean[]) getParam(BANDS);

		type = getIntParam(LANDSAT_TYPE);

	}

	private boolean loadCoefMatrix(int type) {

		switch (type) {

		case 0:

			matrixParams = LandSatMS;

			if(getNumberOfOutputBands(bands) != 4) {

				return false;

			}

			break;	

		case 1:

			matrixParams = LandSatTM;	

			if(getNumberOfOutputBands(bands) != 6) {

				return false;

			}

			break;

		case 2:

			matrixParams = LandSatETM;	

			if(getNumberOfOutputBands(bands) != 6) {

				return false;

			}

			break;

		}	

		return true;

	}

	public void process() throws ProcessInterruptedException, ProcessException {

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

		if (store == null) {

			throw new TasseledCapException(Messages.getText("need_a_input"));

		}

		if(!loadCoefMatrix(type)) {

			throw new TasseledCapException(Messages.getText("wrong_number_of_bands"));

		}

		if(getROIEPSG() != null) {

			try {

				rois = store.getRois(getROIEPSG());

			} catch (ROIException e2) {

				logger.error(Messages.getText("error_getting_rois"), e2);

			}

		}

		try {

			store = ((RasterDataStore)store).cloneDataStore();

		} catch (CloneException e) {

			new TasseledCapException("Error cloning the input DataStore", e);

		}

		extentResult = getExtentResult(getOutputWindow(), rois, store);

		Rectangle2D sourcePxBBox = getSourcePxBox(extentResult, store);

		double cellSize = store.getCellSize();

		nodata = RasterLocator.getManager().getDataStructFactory().createDefaultNoData(1, Buffer.TYPE_DOUBLE);

		double nodataValue = nodata.getValue().doubleValue();

		Buffer outputBuffer = createOutputBuffer((int)sourcePxBBox.getWidth(), (int)sourcePxBBox.getHeight(), matrixParams.length, Buffer.TYPE_FLOAT);

		Buffer sourceBuffer = createSourceBuffer(store, sourcePxBBox, bands);

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

		float resultValue = 0;

		for (int i = 0; i < sourceBuffer.getHeight(); i++) {

			for (int j = 0; j < sourceBuffer.getWidth(); j++) {

				if (isInsideOfROI(j, i, rois, extentResult)) {

					for (int coef = 0; coef < matrixParams.length; coef++) {

						for(int band = 0; band < matrixParams[0].length; band++) {	

							resultValue += getValue(sourceBuffer, i, j, band, coef);

						}

						outputBuffer.setElem(i, j, coef, resultValue);

						resultValue = 0;

					}

				} else {

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

						outputBuffer.setElem(i, j, iBand, (float)nodataValue);

					}

				}

			}

			updatePercent((int)(i * 100 / sourcePxBBox.getHeight()), 100);

		}

        super.exportRaster(filename, outputBuffer, cellSize, extentResult.getULX(),

            extentResult.getULY(), nodata, store.getProjection());

		addOutputValue(FILENAME, filename);

		addOutputValue(GRAPHIC_DATA, filename);

	}

	/**

	 * Gets the number of the selected bands

	 * @param bandsPCs

	 * @return

	 */

	private int getNumberOfOutputBands(boolean[] b) {

		int bandCount = 0;

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

			if(b[i])

				bandCount++;

		}

        return bandCount;

	}

	private float getValue(Buffer sourceBuffer, int i, int j, int band, int coef) {

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

			return (float)((sourceBuffer.getElemByte(i, j, band) & 0xff) * matrixParams[coef][band]);

		}

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

			return (float)(sourceBuffer.getElemShort(i, j, band) * matrixParams[coef][band]);

		}

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

			return (float)(sourceBuffer.getElemInt(i, j, band) * matrixParams[coef][band]);

		}

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

			return (float)(sourceBuffer.getElemFloat(i, j, band) * matrixParams[coef][band]);

		}

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

			return (float)(sourceBuffer.getElemDouble(i, j, band) * matrixParams[coef][band]);

		}