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.statistics;

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

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

/**

 * Proceso del filtro de recorte de colas aplicado a im?genes double

 *

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

 */

public class TailTrimDoubleFilter extends TailTrimFilter {

        public TailTrimDoubleFilter() {}

        public void pre() throws FilterAddException {

                super.pre();

                sampleDec = new double[raster.getBandCount()][nSamples];

                result = new double[raster.getBandCount()][2];

        }

        public void process(int col, int line) {

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

                        sampleDec[iBand][count] = raster.getElemDouble(line, col, iBand);

                count++;

        }

        public int getInRasterDataType() {

                return Buffer.TYPE_DOUBLE;

        }

        public int getOutRasterDataType() {

                return Buffer.TYPE_DOUBLE;

        }

        public void post() {

                super.post();

                //Cogemos el minimo y m?ximo para cada banda

                for (int i = 0; i < raster.getBandCount(); i++) {

                        result[i][0] = sampleDec[i][posInit + tailSize];

                        result[i][1] = sampleDec[i][(posInit + nSamples) - tailSize];

                }

                stats.setTailTrimValue(tailPercent, result);

                //Cogemos el minimo y m?ximo para cada banda

                for (int iValue = 0; iValue < tailSizeList.length; iValue++) {

                        double [][] res = new double[raster.getBandCount()][2];

                        for (int i = 0; i < raster.getBandCount(); i++) {

                                res[i][0] = sampleDec[i][posInit + tailSizeList[iValue]];

                                res[i][1] = sampleDec[i][(posInit + nSamples) - tailSizeList[iValue]];

                        }

                        stats.setTailTrimValue(tailPercentList[iValue], res);

                }

        }

}