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 java.util.Arrays;

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

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

import org.gvsig.raster.impl.DefaultRasterManager;

import org.gvsig.raster.impl.buffer.RasterBuffer;

/**

 * Filtro de Mediana que se aplica en la imagen. Toma como entrada la imagen y

 * el lado de la ventana de filtrado.

 *

 * @author Diego Guerrero Sevilla <diego.guerrero@uclm.es>

 */

public class MedianIntegerFilter extends MedianFilter {

        private int[] window = null;

        public MedianIntegerFilter() {

                super();

        }

        public void pre() throws FilterAddException {

                super.pre();

                window = new int[sizeWindow];

                rasterResult = DefaultRasterManager.getInstance().createBuffer(Buffer.TYPE_INT, raster.getWidth(), raster.getHeight(), raster.getBandCount(), true);

        }

        public void process(int col, int line) {

                // Obtener el vector con la ventanas de muestras

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

                        int k = 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.getElemInt(line + j, col + i, band) & 0xff;

                                                k++;

                                        }

                                }

                        // Ordenar los valores de las ventanas, se supone que usa quickSort.

                        Arrays.sort(window, 0, k);

                        // Extraer los elementos centrales y asignarselos al pixel (x,y)

                        rasterResult.setElem(line, col, band, (int) window[k >> 1]);

                }

        }

        public int getInRasterDataType() {

                return RasterBuffer.TYPE_INT;

        }

        public int getOutRasterDataType() {

                return RasterBuffer.TYPE_INT;

        }

}