Application: gvSIG desktop

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional and Generalitat Valenciana.

 *

 * 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., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Iba?ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package org.gvsig.remotesensing.processtest;

import junit.framework.TestCase;

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

import org.gvsig.raster.RasterLibrary;

import org.gvsig.raster.buffer.BufferFactory;

import org.gvsig.raster.buffer.RasterBufferInvalidException;

import org.gvsig.raster.dataset.IBuffer;

import org.gvsig.raster.grid.Grid;

import org.gvsig.remotesensing.principalcomponents.PCImageProcess;

import Jama.Matrix;

import com.iver.cit.gvsig.exceptions.layers.LoadLayerException;

/**

* Este test prueba el proceso de construcci?n de la imagen resultante del proceso

* de analisis de componentes principales de la imagen pc_CreateImageTest.tif de 

* dimensiones 4x4 y tres bandas. 

* 

* ** @author Alejandro Mu?oz Sanchez (alejandro.munoz@uclm.es)

* */

public class TPCImageProcess extends TestCase {

        private String baseDir = "./test-images/";

        private String path1 = baseDir + "pc_CreateImageTest.tif";

        FLyrRasterSE lyr = null;

        Grid dataGrid=null;

        IBuffer   result= null;

        static{

                RasterLibrary.wakeUp();        

        }

        public void start() {

                this.setUp();

                this.testStack();

        }

        public void setUp() {

                try {

                        lyr = FLyrRasterSE.createLayer(

                                        path1,

                                        path1,

                                        null

                                        );

                        BufferFactory ds1 = new BufferFactory(lyr.getDataSource());

                        dataGrid= new Grid(ds1);

                } catch (LoadLayerException e) {

                        System.out.print("Error en la construcci?n de la capa");

                } catch (RasterBufferInvalidException e) {

                        e.printStackTrace();

                }

        }

        public void testStack() {

                PCImageProcess iProcess = new PCImageProcess ();

                iProcess.addParam("layer",lyr);

                iProcess.addParam("grid",dataGrid);

                iProcess.addParam("autovalores",new double[]{425.104,1753.227,18714.369});

                iProcess.addParam("autovectores", new Matrix(new double[][]{{-0.922766631036,0.2990630039783, 0.2430289371596},

                                                                                                                                   {-0.030678162029,-0.685664096085, 0.7272713370632},

                                                                           {0.3841361672896,0.6636460404393, 0.6418826512605}}));

                iProcess.addParam("compselected",new boolean[]{true,true,true});

                iProcess.run();

                result=iProcess.getBufferResult();

                compare();

        }

        private void compare() {

                // Mismo numero de bandas

                assertEquals(result.getBandCount(),3);

                //Comparaci?n de valores

                assertEquals(result.getElemFloat(0, 0,0),1977.6109,1);

                assertEquals(result.getElemFloat(0, 1,0),1864.9603,1);

                assertEquals(result.getElemFloat(0, 2,0),1896.0991,1);

                assertEquals(result.getElemFloat(0, 3,0),1768.2855,1);        

                assertEquals(result.getElemFloat(1, 0,0),2081.7544,1);

                assertEquals(result.getElemFloat(1, 1,0),2062.1047,1);

                assertEquals(result.getElemFloat(1, 2,0),1989.8268,1);

                assertEquals(result.getElemFloat(1, 3,0),1777.2758,1);        

                // Completar con resultado de otras dos bandas

        }

}