| gvSIG desktop 1

package org.gvsig.fmap.raster.grid.roi;

import junit.framework.TestCase;

import org.gvsig.raster.buffer.RasterBufferInvalidException;

import org.gvsig.raster.dataset.IBuffer;

import org.gvsig.raster.grid.Grid;

import org.gvsig.raster.grid.GridExtent;

import com.iver.cit.gvsig.fmap.core.GeneralPathX;

import com.iver.cit.gvsig.fmap.core.IGeometry;

import com.iver.cit.gvsig.fmap.core.ShapeFactory;

/**Test calculo de estadisticas avanzadas (Matriz de varianza covarianza) para una ROI*/

public class TestRoiStatisticAdvanced extends TestCase {

        Grid inputGrid = null;

        GridExtent extent = null;

        int cellSize = 1;

        double xMin = 0;

        double xMax = 3;

        double yMin = 0;

        double yMax = 3;

        int bands[] ={0,1};

        VectorialROI roi = null;

        GeneralPathX path = null;

        IGeometry geometry = null;

        static{

        //        RasterLibrary.wakeUp();

        }

        public void start() {

                this.setUp();

                this.testStack();

        }

        public void setUp() {

                int bands[] = {0,1};

                extent = new GridExtent();

                extent.setCellSize(cellSize);

                extent.setXRange(xMin, xMax);

                extent.setYRange(yMin, yMax);

                try {

                        inputGrid = new Grid(extent, extent, IBuffer.TYPE_BYTE, bands);

                        //inputGrid.getRasterBuf().setLineInBandByte(new byte[]{94,3,50}, 0, 0);

                        inputGrid.getRasterBuf().setLineInBandByte(new byte[]{68,66,64}, 1, 0);

                        inputGrid.getRasterBuf().setLineInBandByte(new byte[]{74,66,62}, 2, 0);

                        //inputGrid.getRasterBuf().setLineInBandByte(new byte[]{48,40,96}, 0, 1);

                        inputGrid.getRasterBuf().setLineInBandByte(new byte[]{84,80,78}, 1, 1);

                        inputGrid.getRasterBuf().setLineInBandByte(new byte[]{90,84,78}, 2, 1);

                } catch (RasterBufferInvalidException e) {

                        // TODO Auto-generated catch block

                        e.printStackTrace();

                }

        }

        public void testStack(){

                VectorialROI r= new VectorialROI(inputGrid);

                path = new GeneralPathX();

                path.moveTo(3,3);

                path.lineTo(0,3);

                path.lineTo(0,0);

                path.lineTo(3,0);

                path.closePath();

                geometry=ShapeFactory.createPolygon2D(path);

                r.addGeometry(geometry);

                double[][] resultado ={{14.2222,15.11111},{15.1111,17.8888}};

                double matrix[][]=r.getVarCovMatrix();

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

                        for(int j=0; j<matrix[0].length;j++){

                                assertEquals(matrix[i][j],resultado[i][j],0.01);

                        }

        }

}