Application: gvSIG desktop » gvSIG raster

/* 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.impl.grid.filter.correction;

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

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

/**

 * Proceso que aplica el filtro de Moda a un raster de tipo double

 * 

 * 23/07/2008

 * @author Nacho Brodin nachobrodin@gmail.com

 */

public class ModeDoubleFilter extends ModeFilter {

        private double[]                   window = null;

        private double                     tempValue       = 0;

        public ModeDoubleFilter() {

                super();

        }

        public void pre() throws FilterAddException {

                super.pre();

                window = new double[sizeWindow];

        }

        public void process(int col, int line) {

                for (int band = 0; band < numberOfBandsToProcess(); band++) {

                        int k = 0;

                        count = 0;

                        for (int i = -halfSide; i <= halfSide; i++) {

                                for (int j = -halfSide; j <= halfSide; j++) {

                                        if ((col + i >= 0) && (line + j >= 0) && (col + i < width) && (line + j < height)) {

                                                window[k] = raster.getElemDouble(line + j, col + i, band);

                                                if(i == -halfSide && j == -halfSide) 

                                                        tempValue = window[k];

                                                if(tempValue == window[k])

                                                        count ++;

                                                k++;

                                        }

                                }

                        }

                        if(count > sizeWindow)

                                return;

                        k = 1;

                        while(k < window.length) {

                                int auxCount = 0;

                                for (int i = k; i < window.length; i++) {

                                        if(window[i] == window[k])

                                                auxCount ++;

                                }

                                if(auxCount > count) {

                                        count = auxCount;

                                        tempValue = window[k];

                                        if(count > sizeWindow)

                                                return;

                                }

                                k++;

                        }

                        rasterResult.setElem(line, col, band, (double) tempValue);

                }

        }

        public int getInRasterDataType() {

                return Buffer.TYPE_DOUBLE;

        }

        public int getOutRasterDataType() {

                return Buffer.TYPE_DOUBLE;

        }

}