Statistics
| Revision:

root / trunk / libraries / libRaster / src / org / gvsig / raster / grid / filter / convolution / ConvolutionDoubleFilter.java @ 21803

History | View | Annotate | Download (3.47 KB)

1
/* gvSIG. Sistema de Informaci?n Geogrfica de la Generalitat Valenciana
2
 *
3
 * Copyright (C) 2006 Instituto de Desarrollo Regional and Generalitat Valenciana.
4
 *
5
 * This program is free software; you can redistribute it and/or
6
 * modify it under the terms of the GNU General Public License
7
 * as published by the Free Software Foundation; either version 2
8
 * of the License, or (at your option) any later version.
9
 *
10
 * This program is distributed in the hope that it will be useful,
11
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13
 * GNU General Public License for more details.
14
 *
15
 * You should have received a copy of the GNU General Public License
16
 * along with this program; if not, write to the Free Software
17
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.
18
 *
19
 * For more information, contact:
20
 *
21
 *  Generalitat Valenciana
22
 *   Conselleria d'Infraestructures i Transport
23
 *   Av. Blasco Iba?ez, 50
24
 *   46010 VALENCIA
25
 *   SPAIN
26
 *
27
 *      +34 963862235
28
 *   gvsig@gva.es
29
 *      www.gvsig.gva.es
30
 *
31
 *    or
32
 *
33
 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)
34
 *   Campus Universitario s/n
35
 *   02071 Alabacete
36
 *   Spain
37
 *
38
 *   +34 967 599 200
39
 */
40
package org.gvsig.raster.grid.filter.convolution;
41

    
42
import org.gvsig.raster.buffer.RasterBuffer;
43
import org.gvsig.raster.dataset.IBuffer;
44

    
45
/**
46
 * Filtro de convolucion para Buffer de tipo Double
47
 * @author Alejandro Mu?oz        <alejandro.munoz@uclm.es> 
48
 * @author Diego Guerrero Sevilla  <diego.guerrero@uclm.es> 
49
 * */
50

    
51
public class ConvolutionDoubleFilter extends ConvolutionFilter {
52

    
53
        public ConvolutionDoubleFilter(){
54
                super();
55
        }
56
        /** 
57
         * @param Kernel a aplicar. En caso de que no se trate de un kernel definido en ConvolutionFilter, se puede pasar como 
58
         * parametro el kernel se pretende aplicar.
59
         * **/
60
        public ConvolutionDoubleFilter(Kernel k){
61
                super();
62
                super.kernel=k;
63
        }
64
        
65
        
66
        public void pre(){
67
                super.pre();
68
                rasterResult = RasterBuffer.getBuffer(IBuffer.TYPE_DOUBLE, raster.getWidth(), raster.getHeight(), raster.getBandCount(), true);
69
        }
70

    
71
        
72
        /** Aplicacion del filtro para el pixel de la posicion line, col */
73
        public void process(int col, int line) {
74
                
75
                int lado= kernel.getLado();
76
                int semiLado = (lado - 1) >> 1;
77
                double ventana[][]= new double[lado][lado];
78
                double resultConvolution=0;
79
                for (int band = 0; band < raster.getBandCount(); band++) {        
80
                        if ((col - semiLado >= 0) && (line - semiLado >= 0) && (col + semiLado < width) && (line + semiLado < height)){
81
                                for (int j = -semiLado; j <= semiLado; j++)
82
                                        for (int i = -semiLado; i <= semiLado; i++)
83
                                                ventana[i + semiLado][j + semiLado] = raster.getElemDouble(line + j, col + i, band);
84
                                Kernel Kventana= new Kernel(ventana);
85
                                resultConvolution=(double)kernel.convolution(Kventana);
86
                                if (resultConvolution<0)
87
                                        resultConvolution = 0;
88
                                rasterResult.setElem(line, col, band,resultConvolution);        
89
                        }
90
                        else rasterResult.setElem(line, col,band,raster.getElemDouble(line,col,band));
91
                }        
92
        }
93
        
94
        /**
95
         * @return  tipo de dato del buffer de entrada
96
         * */
97
        public int getInRasterDataType() {
98
                return RasterBuffer.TYPE_DOUBLE;
99
        }
100

    
101
        /**
102
         * @return  tipo de dato del buffer de salida
103
         * */                
104
        public int getOutRasterDataType() {
105
                return RasterBuffer.TYPE_DOUBLE;
106
        }
107

    
108
        /**
109
         * @return  buffer resultante tras aplicar el filtro
110
         * */
111
        public Object getResult(String name) {
112
                if (name.equals("raster"))
113
                        return this.rasterResult;
114
                return null;
115
        }
116
        
117
}