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 Ibez, 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;

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

import org.gvsig.remotesensing.principalcomponents.gui.PrincipalComponentPanel;

import com.iver.andami.PluginServices;

import com.iver.andami.plugins.Extension;

import com.iver.cit.gvsig.fmap.MapContext;

import com.iver.cit.gvsig.fmap.layers.FLayers;

import com.iver.cit.gvsig.project.documents.view.IProjectView;

import com.iver.cit.gvsig.project.documents.view.gui.View;

/**

 * Extensi?n para el c?lculo de Componentes Principales.

 *

 * @author Diego Guerrero Sevilla (diego.guerrero@uclm.es)

 *

 */

public class PrincipalComponentsExtension extends Extension {

        private static final double NODATA = -99999999999999.;

        public void initialize() {

                //PrincipalComponentListManager.register();

        }

        public void execute(String actionCommand) {

                if (actionCommand.equals("principal_components")){

                        com.iver.andami.ui.mdiManager.IWindow activeWindow = PluginServices.getMDIManager().getActiveWindow();

                        //si la ventana activa es de tipo Vista

                        if (activeWindow instanceof View) {

                                PrincipalComponentPanel pcPanel = new PrincipalComponentPanel ((View)activeWindow);

                                PluginServices.getMDIManager().addWindow(pcPanel);

                        }

                }

        }

        public boolean isEnabled() {

                com.iver.andami.ui.mdiManager.IWindow f = PluginServices.getMDIManager().getActiveWindow();

                if (f == null) {

                        return false;

                }

                if (f.getClass() == View.class) {

                        View vista = (View) f;

                        IProjectView model = vista.getModel();

                        MapContext mapa = model.getMapContext();

                        FLayers layers = mapa.getLayers();

                        for (int i = 0; i < layers.getLayersCount(); i++)

                                if (layers.getLayer(i) instanceof FLyrRasterSE)

                                        return true;

                }

                return false;

        }

        public boolean isVisible() {

                com.iver.andami.ui.mdiManager.IWindow f = PluginServices.getMDIManager().getActiveWindow();

                if (f == null) {

                        return false;

                }

                if (f instanceof View) {

                        View vista = (View) f;

                        IProjectView model = vista.getModel();

                        MapContext mapa = model.getMapContext();

                        return mapa.getLayers().getLayersCount() > 0;

                } else {

                        return false;

                }

        }

        private static double covariance(double[] xx, double x, double[] yy, double y) {

                double dSum = 0;

                int iValues = 0;

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

                        if (xx[i] != NODATA && yy[i] != NODATA){

                                dSum += (xx[i] - x) * (yy[i] - y);

                                iValues++;

                        }

                }

                if (iValues > 1){

                        return dSum / (double) (iValues - 1);

                }

                else{

                        return NODATA;

                }

        }

        private static double covariance(byte[][] xx, double x, byte[][] yy, double y) {

                double dSum = 0;

                int iValues = 0;

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

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

                                if (xx[i][j] != NODATA && yy[i][j] != NODATA){

                                        dSum += (xx[i][j] - x) * (yy[i][j] - y);

                                        iValues++;

                                }

                        }

                }

                if (iValues > 1){

                        return dSum / (double) (iValues - 1);

                }

                else{

                        return NODATA;

                }

        }

        //Covariance of two 1D arrays of doubles, xx and yy

        public static double covariance(byte[][] xx, byte[][] yy){

        int nx = xx.length;

        int ny = xx[0].length;

        if(nx!=yy.length)throw new IllegalArgumentException("length of x variable array, " + nx + " and length of y array, " + yy.length + " are different");

        double sumx=0.0D, meanx=0.0D;

        double sumy=0.0D, meany=0.0D;

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

                for(int j=0; j<ny; j++){

                sumx+=xx[i][j];

                sumy+=yy[i][j];

                }

        }

        meanx=sumx/((double)nx*(double)ny);

        meany=sumy/((double)nx*(double)ny);

        double sum=0.0D;

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

                for(int j=0; j<ny; j++)

                sum+=(xx[i][j]-meanx)*(yy[i][j]-meany);

        }

        return sum/((double)(nx*ny-1));

        }

}