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

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

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

/**

 * Filtro de ecualizaci?n de histograma para tipo de datos short.

 *

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

 */

public class EqualizationShortFilter extends EqualizationFilter {

        public void process(int col, int line) {

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

                        short p = raster.getElemShort(line, col, iBand);

                        if(!equalizationActive(iBand)) {

                                rasterResult.setElem(line, col, iBand, p);

                                continue;

                        }

                        /*if (p > maxBandValue[renderBands[iBand]])

                                p = (short) maxBandValue[renderBands[iBand]];

                        else if (p < minBandValue[renderBands[iBand]])

                                p = (short) minBandValue[renderBands[iBand]];*/

                        int pos = (int)(((p + dto[iBand]) * nClasses) / distance[iBand]);

                        int ecualizationPositive = (int)(lahe[iBand][pos]);

                        int ecualizationNegative = (int)(lahe[iBand][nElements - pos]);

                        double value = ((nElements - ecualizationNegative) + ecualizationPositive) / 2;

                        rasterResult.setElem(line, col, iBand, (short)value);

                }

        }

        double[] distance = null;

        double[] dto = null;

        public void pre() throws FilterAddException {

                super.pre();

                distance = new double[raster.getBandCount()];

                dto = new double[raster.getBandCount()];

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

                        distance[i] = maxBandValue[renderBands[i]] - minBandValue[renderBands[i]];

                        dto[i] = -minBandValue[renderBands[i]];

                }

        }

        public int getInRasterDataType() {

                return Buffer.TYPE_SHORT;

        }

        public int getOutRasterDataType() {

                return Buffer.TYPE_SHORT;

        }

}