Revision 3578
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/RasterFilterStack.java | ||
|---|---|---|
| 24 | 24 |
package org.cresques.io.raster; |
| 25 | 25 |
|
| 26 | 26 |
import java.awt.Image; |
| 27 |
|
|
| 28 | 27 |
import java.util.Vector; |
| 29 | 28 |
|
| 29 |
import org.cresques.geo.ViewPortData; |
|
| 30 |
import org.cresques.px.Extent; |
|
| 30 | 31 |
|
| 32 |
|
|
| 31 | 33 |
/** |
| 32 | 34 |
* Esta clase representa la pila de filtros que debe ser manejada desde el |
| 33 | 35 |
* RasterFilterStackManager. |
| ... | ... | |
| 35 | 37 |
* |
| 36 | 38 |
*/ |
| 37 | 39 |
public class RasterFilterStack {
|
| 38 |
private Image image = null; |
|
| 39 |
private RasterBuf rasterBuf = null; |
|
| 40 |
private RasterFilter rasterFilter = null; |
|
| 41 |
private int[] order = null; |
|
| 42 |
private RasterStats stats = null; |
|
| 43 |
Vector stack = new Vector(); |
|
| 44 |
private int typeFilter = -2; |
|
| 40 |
private Image image = null; |
|
| 41 |
private RasterBuf rasterBuf = null; |
|
| 42 |
private RasterFilter rasterFilter = null; |
|
| 43 |
private int[] order = null; |
|
| 44 |
private RasterStats stats = null; |
|
| 45 |
private Vector stack = new Vector(); |
|
| 46 |
private int typeFilter = -2; |
|
| 47 |
private ViewPortData viewPortData = null; |
|
| 48 |
private Extent extent = null; |
|
| 45 | 49 |
|
| 46 | 50 |
/** |
| 47 | 51 |
* Constructor |
| ... | ... | |
| 238 | 242 |
public RasterFilter get(int i) {
|
| 239 | 243 |
return ((Filter) stack.get(i)).filter; |
| 240 | 244 |
} |
| 245 |
|
|
| 246 |
/** |
|
| 247 |
* Obtiene el filtro apilado que corresponde con el tipo |
|
| 248 |
* @param type Tipo de filtro |
|
| 249 |
* @return Filtro |
|
| 250 |
*/ |
|
| 251 |
public RasterFilter getByType(int type) {
|
|
| 252 |
for(int i=0;i<lenght();i++){
|
|
| 253 |
if(((Filter) stack.get(i)).type == type) |
|
| 254 |
return ((Filter) stack.get(i)).filter; |
|
| 255 |
} |
|
| 256 |
return null; |
|
| 257 |
} |
|
| 241 | 258 |
|
| 242 | 259 |
/** |
| 243 | 260 |
* Obtiene el tipo del filtro de la pila de la posici?n i |
| ... | ... | |
| 343 | 360 |
} else {
|
| 344 | 361 |
filter.addParam("raster", rasterBuf);
|
| 345 | 362 |
} |
| 346 |
|
|
| 363 |
|
|
| 364 |
filter.setExtent(extent); |
|
| 365 |
filter.setViewPortData(viewPortData); |
|
| 347 | 366 |
filter.execute(); |
| 348 | 367 |
|
| 349 | 368 |
if (filter.getResult("stats") != null) {
|
| ... | ... | |
| 381 | 400 |
this.rasterBuf = rasterBuf; |
| 382 | 401 |
executeFilterByDataType(RasterBuf.TYPE_SHORT); |
| 383 | 402 |
} |
| 384 |
|
|
| 403 |
|
|
| 385 | 404 |
/** |
| 386 | 405 |
*Muestra el contenido de la pila de filtros para depuraci?n |
| 387 | 406 |
*/ |
| ... | ... | |
| 412 | 431 |
this.type = type; |
| 413 | 432 |
} |
| 414 | 433 |
} |
| 434 |
|
|
| 435 |
/** |
|
| 436 |
* @param viewPortData The viewPortData to set. |
|
| 437 |
*/ |
|
| 438 |
public void setViewPortData(ViewPortData viewPortData) {
|
|
| 439 |
this.viewPortData = viewPortData; |
|
| 440 |
} |
|
| 441 |
/** |
|
| 442 |
* @param extent The extent to set. |
|
| 443 |
*/ |
|
| 444 |
public void setExtent(Extent extent) {
|
|
| 445 |
this.extent = extent; |
|
| 446 |
} |
|
| 415 | 447 |
} |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/ColorSpaceConversion.java | ||
|---|---|---|
| 1 |
/* |
|
| 2 |
* Cresques Mapping Suite. Graphic Library for constructing mapping applications. |
|
| 3 |
* |
|
| 4 |
* Copyright (C) 2004-5. |
|
| 5 |
* |
|
| 6 |
* This program is free software; you can redistribute it and/or |
|
| 7 |
* modify it under the terms of the GNU General Public License |
|
| 8 |
* as published by the Free Software Foundation; either version 2 |
|
| 9 |
* of the License, or (at your option) any later version. |
|
| 10 |
* |
|
| 11 |
* This program is distributed in the hope that it will be useful, |
|
| 12 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of |
|
| 13 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|
| 14 |
* GNU General Public License for more details. |
|
| 15 |
* |
|
| 16 |
* You should have received a copy of the GNU General Public License |
|
| 17 |
* along with this program; if not, write to the Free Software |
|
| 18 |
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,USA. |
|
| 19 |
* |
|
| 20 |
* For more information, contact: |
|
| 21 |
* |
|
| 22 |
* cresques@gmail.com |
|
| 23 |
*/ |
|
| 24 |
package org.cresques.io.raster; |
|
| 25 |
|
|
| 26 |
|
|
| 27 |
|
|
| 28 |
/** |
|
| 29 |
* |
|
| 30 |
* @author Nacho Brodin (brodin_ign@gva.es) |
|
| 31 |
*/ |
|
| 32 |
|
|
| 33 |
public class ColorSpaceConversion {
|
|
| 34 |
public static int H = 0, R = 0; |
|
| 35 |
public static int S = 1, G = 1; |
|
| 36 |
public static int L = 2, B = 2; |
|
| 37 |
|
|
| 38 |
|
|
| 39 |
/** |
|
| 40 |
* Convierte HSL a RGB |
|
| 41 |
* @param h hue |
|
| 42 |
* @param s saturation |
|
| 43 |
* @param l ligthness |
|
| 44 |
* @return |
|
| 45 |
*/ |
|
| 46 |
public static int[] HSLtoRGB(int h, int s, int i){
|
|
| 47 |
double red; /* the red band output */ |
|
| 48 |
double red255; /* the red band output */ |
|
| 49 |
double green; /* the green band output */ |
|
| 50 |
double green255; /* the green band output */ |
|
| 51 |
double blue; /* the blue band output */ |
|
| 52 |
double blue255; /* the blue band output */ |
|
| 53 |
double m1; /* value used for determining RGB */ |
|
| 54 |
double m2; /* value used for determining RGB */ |
|
| 55 |
double scalei; /* intensity value */ |
|
| 56 |
double scales; /* saturation value */ |
|
| 57 |
double hue; /* hue */ |
|
| 58 |
double savehue; /* save the hue for future processing */ |
|
| 59 |
int[] res = new int[3]; |
|
| 60 |
|
|
| 61 |
if(h == s && h == 0){
|
|
| 62 |
res[0] = res[1] = res[2] = i; |
|
| 63 |
return res; |
|
| 64 |
} |
|
| 65 |
|
|
| 66 |
red = green = blue = 0.0; |
|
| 67 |
scalei = (double) (i / 255.0); |
|
| 68 |
scales = (double) (s / 255.0); |
|
| 69 |
m2 = 0.0; |
|
| 70 |
|
|
| 71 |
if (scalei <= 0.50) |
|
| 72 |
m2 = scalei*(1.0 + scales); |
|
| 73 |
else if (scalei > 0.50) |
|
| 74 |
m2 = scalei+scales-(scalei*scales); |
|
| 75 |
m1 = 2.0 * scalei - m2; |
|
| 76 |
|
|
| 77 |
hue = (double) 360.0 * h / 255.0; |
|
| 78 |
|
|
| 79 |
if (scales == 0.0){
|
|
| 80 |
if (hue == -1.0){
|
|
| 81 |
red = scalei; |
|
| 82 |
green = scalei; |
|
| 83 |
blue = scalei; |
|
| 84 |
} |
|
| 85 |
}else{
|
|
| 86 |
savehue = hue + 120.0; |
|
| 87 |
if (savehue > 360.0) |
|
| 88 |
savehue -= 360.0; |
|
| 89 |
if (savehue < 0.0) |
|
| 90 |
savehue += 360.0; |
|
| 91 |
if (savehue < 60.0) |
|
| 92 |
red = m1 + (m2-m1) * savehue/60.0; |
|
| 93 |
else if (savehue < 180.0) |
|
| 94 |
red = m2; |
|
| 95 |
else if (savehue < 240.0) |
|
| 96 |
red = m1 + (m2-m1) * (240.0-savehue)/60.0; |
|
| 97 |
else |
|
| 98 |
red = m1; |
|
| 99 |
|
|
| 100 |
savehue = hue; |
|
| 101 |
if (savehue > 360.0) |
|
| 102 |
savehue -= 360.0; |
|
| 103 |
if (savehue < 0.0) |
|
| 104 |
savehue += 360.0; |
|
| 105 |
if (savehue < 60.0) |
|
| 106 |
green = m1 + (m2-m1) * savehue/60.0; |
|
| 107 |
else if (savehue < 180.0) |
|
| 108 |
green = m2; |
|
| 109 |
else if (savehue < 240.0) |
|
| 110 |
green = m1 + (m2-m1) * (240.0-savehue)/60.0; |
|
| 111 |
else |
|
| 112 |
green = m1; |
|
| 113 |
|
|
| 114 |
|
|
| 115 |
savehue = hue - 120.0; |
|
| 116 |
if (savehue > 360.0) |
|
| 117 |
savehue -= 360.0; |
|
| 118 |
if (savehue < 0.0) |
|
| 119 |
savehue += 360.0; |
|
| 120 |
if (savehue < 60.0) |
|
| 121 |
blue = m1 + (m2-m1) * savehue/60.0; |
|
| 122 |
else if (savehue < 180.0) |
|
| 123 |
blue = m2; |
|
| 124 |
else if (savehue < 240.0) |
|
| 125 |
blue = m1 + (m2-m1) * (240.0-savehue)/60.0; |
|
| 126 |
else |
|
| 127 |
blue = m1; |
|
| 128 |
} |
|
| 129 |
|
|
| 130 |
red255 = red*255.0; |
|
| 131 |
green255 = green*255.0; |
|
| 132 |
blue255 = blue*255.0; |
|
| 133 |
if (red255 > 255.0) |
|
| 134 |
red = 255.0; |
|
| 135 |
else |
|
| 136 |
red = red255; |
|
| 137 |
if (green255 > 255.0) |
|
| 138 |
green = 255.0; |
|
| 139 |
else |
|
| 140 |
green = green255; |
|
| 141 |
if (blue255 > 255.0) |
|
| 142 |
blue = 255.0; |
|
| 143 |
else |
|
| 144 |
blue = blue255; |
|
| 145 |
|
|
| 146 |
if (red > 254.5) red = 254.5; |
|
| 147 |
if (red < 0.0) red = 0.0; |
|
| 148 |
if (green > 254.5) green = 254.5; |
|
| 149 |
if (green < 0.0) green = 0.0; |
|
| 150 |
if (blue > 254.5) blue = 254.5; |
|
| 151 |
if (blue < 0.0) blue = 0.0; |
|
| 152 |
|
|
| 153 |
res[0] = (int) (red + 0.5); |
|
| 154 |
res[1] = (int) (green +0.5); |
|
| 155 |
res[2] = (int) (blue + 0.5); |
|
| 156 |
return res; |
|
| 157 |
} |
|
| 158 |
|
|
| 159 |
/** |
|
| 160 |
* Convierte RGB a HSL |
|
| 161 |
* @param r red |
|
| 162 |
* @param g green |
|
| 163 |
* @param b blue |
|
| 164 |
* @return |
|
| 165 |
*/ |
|
| 166 |
public static double[] RGBtoHSL(int rojo, int verde, int azul){
|
|
| 167 |
double scaler, scaleg, scaleb; |
|
| 168 |
double sat; |
|
| 169 |
double red, green, blue; |
|
| 170 |
double low, high, intens; |
|
| 171 |
double hue=0.0; |
|
| 172 |
double[] res1 = new double[3]; |
|
| 173 |
|
|
| 174 |
scaler = (double) (rojo / 255.0); |
|
| 175 |
scaleg = (double) (verde / 255.0); |
|
| 176 |
scaleb = (double) (azul / 255.0); |
|
| 177 |
|
|
| 178 |
high = scaler; |
|
| 179 |
|
|
| 180 |
if (scaleg > high) |
|
| 181 |
high = scaleg; |
|
| 182 |
if (scaleb > high) |
|
| 183 |
high = scaleb; |
|
| 184 |
low = scaler; |
|
| 185 |
if (scaleg < low) |
|
| 186 |
low = scaleg; |
|
| 187 |
if (scaleb < low) |
|
| 188 |
low = scaleb; |
|
| 189 |
|
|
| 190 |
intens = ((high + low) / 2.0); |
|
| 191 |
if (high == low){
|
|
| 192 |
sat = 0.0; |
|
| 193 |
hue = 0.0; |
|
| 194 |
res1[0] = (double) hue; |
|
| 195 |
res1[2] = (double) intens; |
|
| 196 |
res1[1] = (double) sat; |
|
| 197 |
}else if (high != low){
|
|
| 198 |
if (intens <= 0.5) |
|
| 199 |
sat = (high-low)/(high+low); |
|
| 200 |
else |
|
| 201 |
sat = (high-low)/(2 - high - low); |
|
| 202 |
|
|
| 203 |
red = (high-scaler)/(high-low); |
|
| 204 |
green = (high-scaleg)/(high-low); |
|
| 205 |
blue =(high-scaleb)/(high-low); |
|
| 206 |
|
|
| 207 |
if (scaler == high) |
|
| 208 |
hue = blue - green; |
|
| 209 |
else |
|
| 210 |
if (scaleg == high) |
|
| 211 |
hue = 2 + red - blue; |
|
| 212 |
else |
|
| 213 |
if (scaleb == high) |
|
| 214 |
hue = 4 + green - red; |
|
| 215 |
|
|
| 216 |
hue *= 60.0; |
|
| 217 |
|
|
| 218 |
if (hue < 0.0) |
|
| 219 |
hue += 360.0; |
|
| 220 |
|
|
| 221 |
//res[0] = (byte)(255.0 * hue / 360.0 + 0.5); |
|
| 222 |
//res[1] = (byte) (intens*255. + 0.5); |
|
| 223 |
//res[2] = (byte) (sat*255. + 0.5); |
|
| 224 |
res1[0] = (double)(hue); |
|
| 225 |
res1[2] = (double) (intens); |
|
| 226 |
res1[1] = (double) (sat); |
|
| 227 |
} |
|
| 228 |
return res1; |
|
| 229 |
} |
|
| 230 |
} |
|
| 0 | 231 | |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/SharpeningFilter.java | ||
|---|---|---|
| 1 |
/* |
|
| 2 |
* Cresques Mapping Suite. Graphic Library for constructing mapping applications. |
|
| 3 |
* |
|
| 4 |
* Copyright (C) 2004-5. |
|
| 5 |
* |
|
| 6 |
* This program is free software; you can redistribute it and/or |
|
| 7 |
* modify it under the terms of the GNU General Public License |
|
| 8 |
* as published by the Free Software Foundation; either version 2 |
|
| 9 |
* of the License, or (at your option) any later version. |
|
| 10 |
* |
|
| 11 |
* This program is distributed in the hope that it will be useful, |
|
| 12 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of |
|
| 13 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|
| 14 |
* GNU General Public License for more details. |
|
| 15 |
* |
|
| 16 |
* You should have received a copy of the GNU General Public License |
|
| 17 |
* along with this program; if not, write to the Free Software |
|
| 18 |
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,USA. |
|
| 19 |
* |
|
| 20 |
* For more information, contact: |
|
| 21 |
* |
|
| 22 |
* cresques@gmail.com |
|
| 23 |
*/ |
|
| 24 |
package org.cresques.io.raster; |
|
| 25 |
|
|
| 26 |
import java.awt.Image; |
|
| 27 |
|
|
| 28 |
import org.cresques.io.GeoRasterFile; |
|
| 29 |
|
|
| 30 |
|
|
| 31 |
/** |
|
| 32 |
* Clase base para los filtros de sharpening en sus diferentes tipos |
|
| 33 |
* de datos. |
|
| 34 |
* @author Nacho Brodin (brodin_ign@gva.es) |
|
| 35 |
* |
|
| 36 |
*/ |
|
| 37 |
public abstract class SharpeningFilter extends RasterFilter {
|
|
| 38 |
|
|
| 39 |
/** |
|
| 40 |
* Posici?n en la lista de GeoRasterFiles de la imagen de refinado. |
|
| 41 |
*/ |
|
| 42 |
protected int posPancromatica = 0; |
|
| 43 |
|
|
| 44 |
/** |
|
| 45 |
* Lista de GeoRasterFiles que componen los datos de entrada |
|
| 46 |
*/ |
|
| 47 |
protected GeoRasterFile[] files = null; |
|
| 48 |
|
|
| 49 |
/** |
|
| 50 |
* Alto de la imagen pancrom?tica |
|
| 51 |
*/ |
|
| 52 |
protected int heightPancr = 0; |
|
| 53 |
|
|
| 54 |
/** |
|
| 55 |
* Ancho de la imagen pancrom?tica |
|
| 56 |
*/ |
|
| 57 |
protected int widthPancr = 0; |
|
| 58 |
|
|
| 59 |
/** |
|
| 60 |
* Alto de las bandas multiespectrales |
|
| 61 |
*/ |
|
| 62 |
protected int heightMultiespec = 0; |
|
| 63 |
|
|
| 64 |
/** |
|
| 65 |
* Ancho de las bandas multiespectrales |
|
| 66 |
*/ |
|
| 67 |
protected int widthMultiespec = 0; |
|
| 68 |
|
|
| 69 |
/** |
|
| 70 |
* Image que corresponde a la imagen de refinado |
|
| 71 |
*/ |
|
| 72 |
protected Image imagePancr = null; |
|
| 73 |
|
|
| 74 |
/** |
|
| 75 |
* Image que corresponde a la imagen de entrada |
|
| 76 |
*/ |
|
| 77 |
protected Image imageMultiespec = null; |
|
| 78 |
|
|
| 79 |
/** |
|
| 80 |
* N?mero de bandas de la operaci?n incluido la de refinado |
|
| 81 |
*/ |
|
| 82 |
protected int nBands = 0; |
|
| 83 |
|
|
| 84 |
/** |
|
| 85 |
* Relaci?n entre la pancromatica y la multiespectral en X |
|
| 86 |
*/ |
|
| 87 |
protected int relX = 0; |
|
| 88 |
|
|
| 89 |
/** |
|
| 90 |
* Relaci?n entre la pancromatica y la multiespectral en Y |
|
| 91 |
*/ |
|
| 92 |
protected int relY = 0; |
|
| 93 |
|
|
| 94 |
/** |
|
| 95 |
* Vector de 3 elementos que corresponden a R, G y B respectivamente. El valor contenido |
|
| 96 |
* en cada uno de ellos corresponde al n?mero de banda de fichero que se visualiza en R, G y B. |
|
| 97 |
*/ |
|
| 98 |
protected int[] bandOrder = null; |
|
| 99 |
|
|
| 100 |
/** |
|
| 101 |
* Valor de alpha aplicado a toda la imagen |
|
| 102 |
*/ |
|
| 103 |
protected int alpha = 0; |
|
| 104 |
|
|
| 105 |
/** |
|
| 106 |
* M?todo de calculo del refinado |
|
| 107 |
*/ |
|
| 108 |
protected String method = "ihs"; |
|
| 109 |
|
|
| 110 |
/** |
|
| 111 |
* Coeficiente |
|
| 112 |
*/ |
|
| 113 |
protected double coef = 0.15; |
|
| 114 |
|
|
| 115 |
/** |
|
| 116 |
* Constructor |
|
| 117 |
* |
|
| 118 |
*/ |
|
| 119 |
public SharpeningFilter() {
|
|
| 120 |
super(); |
|
| 121 |
} |
|
| 122 |
|
|
| 123 |
/** |
|
| 124 |
* |
|
| 125 |
*/ |
|
| 126 |
public void pre() {
|
|
| 127 |
|
|
| 128 |
} |
|
| 129 |
|
|
| 130 |
} |
|
| 0 | 131 | |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/RasterFilter.java | ||
|---|---|---|
| 24 | 24 |
package org.cresques.io.raster; |
| 25 | 25 |
|
| 26 | 26 |
import java.awt.Image; |
| 27 |
|
|
| 28 | 27 |
import java.util.Hashtable; |
| 29 | 28 |
|
| 29 |
import org.cresques.geo.ViewPortData; |
|
| 30 |
import org.cresques.px.Extent; |
|
| 30 | 31 |
|
| 32 |
|
|
| 31 | 33 |
/** |
| 32 | 34 |
* Filtro para raster. Ancestro de todos los filtros. |
| 33 | 35 |
* @author Luis W. Sevilla (sevilla_lui@gva.es) |
| 34 | 36 |
*/ |
| 35 | 37 |
public abstract class RasterFilter implements IRasterFilter {
|
| 36 |
protected RasterBuf raster = null; |
|
| 37 |
protected RasterStats stats = null; |
|
| 38 |
protected Image image = null; |
|
| 39 |
protected int[] px = new int[4]; |
|
| 40 |
protected int height = 0; |
|
| 41 |
protected int width = 0; |
|
| 42 |
protected Hashtable params = new Hashtable(); |
|
| 43 |
protected int incX = 1; |
|
| 44 |
protected int incY = 1; |
|
| 45 |
protected boolean exec = true; //Si est? a false no se ejecuta el filtro |
|
| 38 |
protected RasterBuf raster = null; |
|
| 39 |
protected RasterStats stats = null; |
|
| 40 |
protected Image image = null; |
|
| 41 |
protected int[] px = new int[4]; |
|
| 42 |
protected int height = 0; |
|
| 43 |
protected int width = 0; |
|
| 44 |
protected Hashtable params = new Hashtable(); |
|
| 45 |
protected int incX = 1; |
|
| 46 |
protected int incY = 1; |
|
| 47 |
protected ViewPortData viewPortData = null; |
|
| 48 |
protected Extent extent = null; |
|
| 49 |
|
|
| 50 |
/** |
|
| 51 |
* Variable que control la aplicaci?n o no del filtro. Si est? a false aunque est? en |
|
| 52 |
* la pila el filtro no se ejecutar?. |
|
| 53 |
*/ |
|
| 54 |
protected boolean exec = true; |
|
| 46 | 55 |
|
| 47 | 56 |
/** |
| 57 |
* @author Nacho Brodin (brodin_ign@gva.es) |
|
| 58 |
*/ |
|
| 59 |
public static class Kernel{
|
|
| 60 |
public float[][] kernel = null; |
|
| 61 |
public Kernel(float[][] k){
|
|
| 62 |
this.kernel = k; |
|
| 63 |
} |
|
| 64 |
|
|
| 65 |
public float kernelOperation(Kernel k){
|
|
| 66 |
float res = 0; |
|
| 67 |
for(int i=0;i<3;i++) |
|
| 68 |
for(int j=0;j<3;j++) |
|
| 69 |
res += kernel[i][j] * k.kernel[i][j]; |
|
| 70 |
return res; |
|
| 71 |
} |
|
| 72 |
|
|
| 73 |
} |
|
| 74 |
|
|
| 75 |
/** |
|
| 48 | 76 |
* Constructor |
| 49 | 77 |
*/ |
| 50 | 78 |
public RasterFilter() {
|
| ... | ... | |
| 78 | 106 |
|
| 79 | 107 |
/** |
| 80 | 108 |
* A?ade un par?metro al filtro |
| 109 |
* @param name Clave del par?metro |
|
| 110 |
* @param param Objeto pasado como par?metro |
|
| 81 | 111 |
*/ |
| 82 | 112 |
public void addParam(String name, Object param) {
|
| 83 | 113 |
params.put(name, param); |
| 84 | 114 |
} |
| 85 | 115 |
|
| 86 | 116 |
/** |
| 117 |
* Elimina un par?metro del filtro |
|
| 118 |
* @param name Clave del par?metro a eliminar |
|
| 119 |
*/ |
|
| 120 |
public void removeParam(String name){
|
|
| 121 |
params.remove(name); |
|
| 122 |
} |
|
| 123 |
|
|
| 124 |
/** |
|
| 125 |
* Obtiene un par?metro a partir de la clave |
|
| 126 |
* @param name Par?metro |
|
| 127 |
* @return Par?metro |
|
| 128 |
*/ |
|
| 129 |
public Object getParam(String name){
|
|
| 130 |
return params.get(name); |
|
| 131 |
} |
|
| 132 |
|
|
| 133 |
/** |
|
| 87 | 134 |
* Funci?n que contiene el c?digo a ejecutar antes de aplicar el filtro |
| 88 | 135 |
*/ |
| 89 | 136 |
abstract public void pre(); |
| ... | ... | |
| 118 | 165 |
* @param name clave para obtener un objeto resultado del filtro. |
| 119 | 166 |
*/ |
| 120 | 167 |
abstract public Object getResult(String name); |
| 168 |
|
|
| 169 |
/** |
|
| 170 |
* @param viewPortData The viewPortData to set. |
|
| 171 |
*/ |
|
| 172 |
public void setViewPortData(ViewPortData viewPortData) {
|
|
| 173 |
this.viewPortData = viewPortData; |
|
| 174 |
} |
|
| 175 |
/** |
|
| 176 |
* @param extent The extent to set. |
|
| 177 |
*/ |
|
| 178 |
public void setExtent(Extent extent) {
|
|
| 179 |
this.extent = extent; |
|
| 180 |
} |
|
| 121 | 181 |
} |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/SharpeningImageFilter.java | ||
|---|---|---|
| 1 |
/* |
|
| 2 |
* Cresques Mapping Suite. Graphic Library for constructing mapping applications. |
|
| 3 |
* |
|
| 4 |
* Copyright (C) 2004-5. |
|
| 5 |
* |
|
| 6 |
* This program is free software; you can redistribute it and/or |
|
| 7 |
* modify it under the terms of the GNU General Public License |
|
| 8 |
* as published by the Free Software Foundation; either version 2 |
|
| 9 |
* of the License, or (at your option) any later version. |
|
| 10 |
* |
|
| 11 |
* This program is distributed in the hope that it will be useful, |
|
| 12 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of |
|
| 13 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|
| 14 |
* GNU General Public License for more details. |
|
| 15 |
* |
|
| 16 |
* You should have received a copy of the GNU General Public License |
|
| 17 |
* along with this program; if not, write to the Free Software |
|
| 18 |
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,USA. |
|
| 19 |
* |
|
| 20 |
* For more information, contact: |
|
| 21 |
* |
|
| 22 |
* cresques@gmail.com |
|
| 23 |
*/ |
|
| 24 |
package org.cresques.io.raster; |
|
| 25 |
|
|
| 26 |
import java.awt.Dimension; |
|
| 27 |
import java.awt.Graphics2D; |
|
| 28 |
import java.awt.Image; |
|
| 29 |
import java.awt.image.BufferedImage; |
|
| 30 |
import java.awt.image.PixelGrabber; |
|
| 31 |
|
|
| 32 |
import org.cresques.geo.ViewPortData; |
|
| 33 |
import org.cresques.io.GeoRasterFile; |
|
| 34 |
import org.cresques.px.PxRaster; |
|
| 35 |
|
|
| 36 |
|
|
| 37 |
/** |
|
| 38 |
* Filtro de sharpening que se aplica sobre la vista. Tiene como entradas el objeto Image |
|
| 39 |
* con las bandas del raster visualizadas en la Vista y la banda pancrom?tica. A partir de |
|
| 40 |
* estas entradas genera un Image de salida como la que hay pero de resoluci?n igual a la |
|
| 41 |
* de la pancrom?tica. |
|
| 42 |
* @author Nacho Brodin (brodin_ign@gva.es) |
|
| 43 |
* |
|
| 44 |
*/ |
|
| 45 |
public class SharpeningImageFilter extends SharpeningFilter {
|
|
| 46 |
|
|
| 47 |
//Variable que dice para cada pixel de la imagen de entrada cuantos de salida hay. |
|
| 48 |
private int nOutputPixelsPerInputPixel; |
|
| 49 |
//Vector con los desplazamientos de los pixeles de salida dentro del vector |
|
| 50 |
private int[] outKernel; |
|
| 51 |
//Kernel aplicado a la pancromatica |
|
| 52 |
private Kernel kernel = null; |
|
| 53 |
private PixelGrabber pgAct = null, pgPrev = null, pgNext = null; |
|
| 54 |
int[] bufferPrev = null, bufferAct = null, bufferNext = null; |
|
| 55 |
double coef = 0.20; |
|
| 56 |
double coefBrovey = 0; |
|
| 57 |
|
|
| 58 |
/** |
|
| 59 |
* Constructor |
|
| 60 |
* |
|
| 61 |
*/ |
|
| 62 |
public SharpeningImageFilter() {
|
|
| 63 |
super(); |
|
| 64 |
} |
|
| 65 |
|
|
| 66 |
/** |
|
| 67 |
* Validamos que haya alguna otra banda adem?s de la pancrom?tica y que la pancrom?tica |
|
| 68 |
* sea de mayor resoluci?n que las dem?s. |
|
| 69 |
*/ |
|
| 70 |
private void checkInput(){
|
|
| 71 |
exec = true; |
|
| 72 |
if(heightMultiespec >= heightPancr || widthMultiespec >= widthPancr || heightMultiespec == 0 || widthMultiespec == 0) |
|
| 73 |
exec = false; |
|
| 74 |
for(int i=0;i<files.length;i++){
|
|
| 75 |
if(i != posPancromatica){
|
|
| 76 |
if(files[i].getHeight() != heightMultiespec || files[i].getWidth() != widthMultiespec){
|
|
| 77 |
exec = false; |
|
| 78 |
break; |
|
| 79 |
} |
|
| 80 |
} |
|
| 81 |
} |
|
| 82 |
} |
|
| 83 |
|
|
| 84 |
/** |
|
| 85 |
* Carga los par?metros pasados al filtro en las variables de |
|
| 86 |
* instancia que corresponde, |
|
| 87 |
*/ |
|
| 88 |
private void loadParam(){
|
|
| 89 |
//Carga de parametros |
|
| 90 |
this.image = (Image) params.get("raster");
|
|
| 91 |
posPancromatica = ((Integer) params.get("pancromatica")).intValue();
|
|
| 92 |
files = (GeoRasterFile[]) params.get("files");
|
|
| 93 |
bandOrder = (int[]) params.get("order");
|
|
| 94 |
alpha = ((Integer) params.get("alpha")).intValue();
|
|
| 95 |
method = (String) params.get("method");
|
|
| 96 |
coef = ((Double) params.get("coef")).doubleValue();
|
|
| 97 |
coefBrovey = ((Integer) params.get("coefBrovey")).intValue();
|
|
| 98 |
|
|
| 99 |
//Asignamos el nombre de la banda pancrom?tica como par?metro para que cuando |
|
| 100 |
//volvamos a abrir el dialogo pueda recuperarse y seleccionarse en el cuadro |
|
| 101 |
this.addParam("pancrName",files[posPancromatica].getName().substring(
|
|
| 102 |
files[posPancromatica].getName().lastIndexOf("/")+1,
|
|
| 103 |
files[posPancromatica].getName().length())); |
|
| 104 |
|
|
| 105 |
height = image.getHeight(null); |
|
| 106 |
width = image.getWidth(null); |
|
| 107 |
heightPancr = files[posPancromatica].getHeight(); |
|
| 108 |
widthPancr = files[posPancromatica].getWidth(); |
|
| 109 |
for(int i=0;i<files.length;i++){
|
|
| 110 |
if(i != posPancromatica){
|
|
| 111 |
heightMultiespec = files[i].getHeight(); |
|
| 112 |
widthMultiespec = files[i].getWidth(); |
|
| 113 |
} |
|
| 114 |
} |
|
| 115 |
|
|
| 116 |
relX = (int)widthPancr/widthMultiespec; |
|
| 117 |
relY = (int)heightPancr/heightMultiespec; |
|
| 118 |
} |
|
| 119 |
|
|
| 120 |
/* (non-Javadoc) |
|
| 121 |
* @see org.cresques.io.raster.IRasterFilter#pre() |
|
| 122 |
*/ |
|
| 123 |
public void pre() {
|
|
| 124 |
|
|
| 125 |
loadParam(); |
|
| 126 |
checkInput(); |
|
| 127 |
|
|
| 128 |
//Creamos el buffer donde se pinta la pancromatica |
|
| 129 |
|
|
| 130 |
imagePancr = new BufferedImage(width, height, BufferedImage.TYPE_INT_ARGB); |
|
| 131 |
Graphics2D graphicsPancr = (Graphics2D)imagePancr.getGraphics(); |
|
| 132 |
PxRaster rasterPancr = new PxRaster(files[posPancromatica], null, files[posPancromatica].getView()); |
|
| 133 |
|
|
| 134 |
rasterPancr.setBand(GeoRasterFile.RED_BAND, 0); |
|
| 135 |
rasterPancr.setBand(GeoRasterFile.GREEN_BAND, 1); |
|
| 136 |
rasterPancr.setBand(GeoRasterFile.BLUE_BAND, 2); |
|
| 137 |
|
|
| 138 |
ViewPortData vp = (ViewPortData)viewPortData.clone(); |
|
| 139 |
vp.zoom(extent); |
|
| 140 |
rasterPancr.draw(graphicsPancr, vp); |
|
| 141 |
rasterPancr = null; |
|
| 142 |
|
|
| 143 |
//Creamos el buffer donde se pinta la imagen de entrada |
|
| 144 |
|
|
| 145 |
imageMultiespec = new BufferedImage(width, height, BufferedImage.TYPE_INT_RGB); |
|
| 146 |
Graphics2D graphicsMultiespec = (Graphics2D)imageMultiespec.getGraphics(); |
|
| 147 |
PxRaster rasterMultiespec = null; |
|
| 148 |
boolean first = true; |
|
| 149 |
for(int i = 0;i<files.length;i++){
|
|
| 150 |
if(first){
|
|
| 151 |
rasterMultiespec = new PxRaster(files[i], null, files[posPancromatica].getView()); |
|
| 152 |
first = false; |
|
| 153 |
}else |
|
| 154 |
rasterMultiespec.addFile(files[i].getName()); |
|
| 155 |
} |
|
| 156 |
|
|
| 157 |
rasterMultiespec.setBand(GeoRasterFile.RED_BAND, bandOrder[0]); |
|
| 158 |
rasterMultiespec.setBand(GeoRasterFile.GREEN_BAND, bandOrder[1]); |
|
| 159 |
rasterMultiespec.setBand(GeoRasterFile.BLUE_BAND, bandOrder[2]); |
|
| 160 |
|
|
| 161 |
rasterMultiespec.draw(graphicsMultiespec, vp); |
|
| 162 |
rasterMultiespec = null; |
|
| 163 |
|
|
| 164 |
super.pre(); |
|
| 165 |
} |
|
| 166 |
|
|
| 167 |
/** |
|
| 168 |
* Aplica la operaci?n de refinamiento sobre el buffer Multiespectral que contiene |
|
| 169 |
* el RGB que es pasado por par?metro utilizando la pancrom?tica. |
|
| 170 |
* @param bufferInput Buffer rgb |
|
| 171 |
* @param length longitud del buffer de la pancromatica utilizado |
|
| 172 |
* @param iLine l?nea leida de la imagen multiespectral |
|
| 173 |
* @return buffer con el resultado de la operaci?n |
|
| 174 |
*/ |
|
| 175 |
private int[] processBrovey(int[] bufferInput, int iLine){
|
|
| 176 |
double[] hsl; |
|
| 177 |
int[] rgb; |
|
| 178 |
PixelGrabber pg = null; |
|
| 179 |
|
|
| 180 |
//Longitud del buffer de salida |
|
| 181 |
int widthDst = width * relX; |
|
| 182 |
//Buffer de salida |
|
| 183 |
int[] bufferPancr = new int[width]; |
|
| 184 |
|
|
| 185 |
pg = new PixelGrabber(imagePancr, 0, iLine, width, 1, bufferPancr , 0, width); |
|
| 186 |
|
|
| 187 |
|
|
| 188 |
try {
|
|
| 189 |
pg.grabPixels(); |
|
| 190 |
} catch (InterruptedException e) {
|
|
| 191 |
e.printStackTrace(); |
|
| 192 |
} |
|
| 193 |
|
|
| 194 |
for(int iElem=0; iElem<width; iElem++){
|
|
| 195 |
int r = ((bufferInput[iElem] >> 16) & 0xff); |
|
| 196 |
int g = ((bufferInput[iElem] >> 8) & 0xff); |
|
| 197 |
int b = (bufferInput[iElem] & 0xff); |
|
| 198 |
byte i = (byte)((bufferPancr[iElem] >> 16) & 0xff); |
|
| 199 |
double scale = (3.0*(i+coefBrovey))/(r+g+b+1.0); |
|
| 200 |
r *= scale;g *= scale;b *= scale; |
|
| 201 |
bufferPancr[iElem] = ((alpha << 24) & 0xff000000) | |
|
| 202 |
((r << 16) & 0x00ff0000)| |
|
| 203 |
((g << 8) & 0x0000ff00) | |
|
| 204 |
(b & 0x000000ff); |
|
| 205 |
} |
|
| 206 |
return bufferPancr; |
|
| 207 |
|
|
| 208 |
} |
|
| 209 |
|
|
| 210 |
/** |
|
| 211 |
* Aplica la operaci?n de refinamiento sobre el buffer Multiespectral que contiene |
|
| 212 |
* el RGB que es pasado por par?metro utilizando la pancrom?tica. |
|
| 213 |
* @param bufferInput Buffer rgb |
|
| 214 |
* @param length longitud del buffer de la pancromatica utilizado |
|
| 215 |
* @param iLine l?nea leida de la imagen multiespectral |
|
| 216 |
* @return buffer con el resultado de la operaci?n |
|
| 217 |
*/ |
|
| 218 |
private int[] processKernel(int[] bufferInput, int iLine){
|
|
| 219 |
int[] bufferDst = new int[width]; |
|
| 220 |
|
|
| 221 |
if(iLine != 0){
|
|
| 222 |
bufferPrev = bufferAct; |
|
| 223 |
bufferAct = bufferNext; |
|
| 224 |
}else{
|
|
| 225 |
bufferPrev = null; |
|
| 226 |
bufferAct = new int[width]; |
|
| 227 |
pgAct = new PixelGrabber(imagePancr, 0, iLine, width, 1, bufferAct , 0, width); |
|
| 228 |
try {
|
|
| 229 |
pgAct.grabPixels(); |
|
| 230 |
} catch (InterruptedException e) {e.printStackTrace();}
|
|
| 231 |
} |
|
| 232 |
|
|
| 233 |
if(iLine != (height - 1)){
|
|
| 234 |
bufferNext = new int[width]; |
|
| 235 |
pgNext = new PixelGrabber(imagePancr, 0, iLine + 1, width, 1, bufferNext , 0, width); |
|
| 236 |
try {
|
|
| 237 |
pgNext.grabPixels(); |
|
| 238 |
} catch (InterruptedException e) {e.printStackTrace();}
|
|
| 239 |
}else |
|
| 240 |
bufferNext = null; |
|
| 241 |
|
|
| 242 |
|
|
| 243 |
for(int iElem=0; iElem<width; iElem++){
|
|
| 244 |
int a = ((bufferInput[iElem] >> 24) & 0xff); |
|
| 245 |
int r = ((bufferInput[iElem] >> 16) & 0xff); |
|
| 246 |
int g = ((bufferInput[iElem] >> 8) & 0xff); |
|
| 247 |
int b = (bufferInput[iElem] & 0xff); |
|
| 248 |
|
|
| 249 |
float[][] op = new float[3][3]; |
|
| 250 |
if(bufferPrev != null){
|
|
| 251 |
if(iElem != 0)op[0][0] = (bufferPrev[iElem - 1] & 0x000000ff); |
|
| 252 |
op[0][1] = (bufferPrev[iElem] & 0x000000ff); |
|
| 253 |
if(iElem != (width -1))op[0][2] = (bufferPrev[iElem + 1] & 0x000000ff); |
|
| 254 |
} |
|
| 255 |
if(iElem != 0)op[1][0] = (bufferAct[iElem - 1] & 0x000000ff); |
|
| 256 |
op[1][1] = (bufferAct[iElem] & 0x000000ff); |
|
| 257 |
if(iElem != (width -1))op[1][2] = (bufferAct[iElem + 1] & 0x000000ff); |
|
| 258 |
if(bufferNext != null){
|
|
| 259 |
if(iElem != 0)op[2][0] = (bufferNext[iElem - 1] & 0x000000ff); |
|
| 260 |
op[2][1] = (bufferNext[iElem] & 0x000000ff); |
|
| 261 |
if(iElem != (width -1))op[2][2] = (bufferNext[iElem + 1] & 0x000000ff); |
|
| 262 |
} |
|
| 263 |
Kernel operando = new Kernel(op); |
|
| 264 |
double i = kernel.kernelOperation(operando) * 0.15; |
|
| 265 |
r += i; |
|
| 266 |
g += i; |
|
| 267 |
b += i; |
|
| 268 |
bufferDst[iElem] = ((a << 24) & 0xff000000) | |
|
| 269 |
((r << 16) & 0x00ff0000)| |
|
| 270 |
((g << 8) & 0x0000ff00) | |
|
| 271 |
(b & 0x000000ff); |
|
| 272 |
} |
|
| 273 |
return bufferDst; |
|
| 274 |
|
|
| 275 |
} |
|
| 276 |
|
|
| 277 |
/** |
|
| 278 |
* Aplica la operaci?n de refinamiento sobre el buffer Multiespectral que contiene |
|
| 279 |
* el RGB que es pasado por par?metro utilizando la pancrom?tica. |
|
| 280 |
* @param bufferInput Buffer rgb |
|
| 281 |
* @param length longitud del buffer de la pancromatica utilizado |
|
| 282 |
* @param iLine l?nea leida de la imagen multiespectral |
|
| 283 |
* @return buffer con el resultado de la operaci?n |
|
| 284 |
*/ |
|
| 285 |
private int[] processIHS(int[] bufferInput, int iLine){
|
|
| 286 |
double[] hsl; |
|
| 287 |
int[] rgb; |
|
| 288 |
PixelGrabber pg = null; |
|
| 289 |
|
|
| 290 |
//Buffer de salida |
|
| 291 |
int[] bufferPancr = new int[width]; |
|
| 292 |
|
|
| 293 |
pg = new PixelGrabber(imagePancr, 0, iLine, width, 1, bufferPancr , 0, width); |
|
| 294 |
|
|
| 295 |
try {
|
|
| 296 |
pg.grabPixels(); |
|
| 297 |
} catch (InterruptedException e) {
|
|
| 298 |
e.printStackTrace(); |
|
| 299 |
} |
|
| 300 |
|
|
| 301 |
int[] uvw , tmp = new int[3]; |
|
| 302 |
double[] xyz; |
|
| 303 |
for(int iElem=0; iElem<width; iElem++){
|
|
| 304 |
xyz = ColorSpaceConversion.RGBtoHSL( (bufferInput[iElem] >> 16) & 0x000000ff, |
|
| 305 |
(bufferInput[iElem] >> 8) & 0x000000ff, |
|
| 306 |
bufferInput[iElem] & 0x000000ff); |
|
| 307 |
|
|
| 308 |
xyz[2] = ((bufferPancr[iElem] & 0x000000ff)/255.0) + coef; |
|
| 309 |
|
|
| 310 |
tmp[0] = (int)(255.0 * xyz[0] / 360.0 + 0.5); |
|
| 311 |
tmp[2] = (int) (xyz[2]*255. + 0.5); |
|
| 312 |
tmp[1] = (int) (xyz[1]*255. + 0.5); |
|
| 313 |
|
|
| 314 |
uvw = ColorSpaceConversion.HSLtoRGB(tmp[ColorSpaceConversion.H], |
|
| 315 |
tmp[ColorSpaceConversion.S], |
|
| 316 |
tmp[ColorSpaceConversion.L]); |
|
| 317 |
bufferPancr[iElem] = ((alpha << 24) & 0xff000000) | |
|
| 318 |
((uvw[0] << 16) & 0x00ff0000)| |
|
| 319 |
((uvw[1] << 8) & 0x0000ff00) | |
|
| 320 |
(uvw[2] & 0x000000ff); |
|
| 321 |
} |
|
| 322 |
return bufferPancr; |
|
| 323 |
|
|
| 324 |
} |
|
| 325 |
|
|
| 326 |
/** |
|
| 327 |
* Aplica el filtro sobre el raster pasado pixel a pixel |
|
| 328 |
*/ |
|
| 329 |
public void execute() {
|
|
| 330 |
pre(); |
|
| 331 |
|
|
| 332 |
if (exec) {
|
|
| 333 |
int[] pRGBArrayMultiesp = new int[width]; |
|
| 334 |
int[] pRGBArrayPancr = null; |
|
| 335 |
PixelGrabber pg = null; |
|
| 336 |
int widthDst = width * relX; |
|
| 337 |
|
|
| 338 |
//Para cada linea leemos los valores RGB del image. Aplicamos el algoritmo |
|
| 339 |
//y escribimos el resultado. |
|
| 340 |
|
|
| 341 |
if(method.equals("hsl")){
|
|
| 342 |
float[][] k = {{-1F, -1F, -1F},{-1F, 8F, -1F},{-1F, -1F, -1F}};
|
|
| 343 |
kernel = new Kernel(k); |
|
| 344 |
for(int iLine=0;iLine<height ;iLine++){
|
|
| 345 |
pg = new PixelGrabber(imageMultiespec, 0, iLine, width, 1, pRGBArrayMultiesp, 0, width); |
|
| 346 |
try {
|
|
| 347 |
pg.grabPixels(); |
|
| 348 |
} catch (InterruptedException e) {e.printStackTrace();}
|
|
| 349 |
pRGBArrayPancr = processIHS(pRGBArrayMultiesp, iLine); |
|
| 350 |
((BufferedImage)imagePancr).setRGB(0, iLine, width, 1, pRGBArrayPancr, 0, width); |
|
| 351 |
} |
|
| 352 |
}else{
|
|
| 353 |
for(int iLine=0;iLine<height ;iLine++){
|
|
| 354 |
pg = new PixelGrabber(imageMultiespec, 0, iLine, width, 1, pRGBArrayMultiesp, 0, width); |
|
| 355 |
try {
|
|
| 356 |
pg.grabPixels(); |
|
| 357 |
} catch (InterruptedException e) {e.printStackTrace();}
|
|
| 358 |
pRGBArrayPancr = processBrovey(pRGBArrayMultiesp, iLine); |
|
| 359 |
((BufferedImage)imagePancr).setRGB(0, iLine, width, 1, pRGBArrayPancr, 0, width); |
|
| 360 |
} |
|
| 361 |
} |
|
| 362 |
} |
|
| 363 |
|
|
| 364 |
post(); |
|
| 365 |
} |
|
| 366 |
|
|
| 367 |
/* (non-Javadoc) |
|
| 368 |
* @see org.cresques.io.raster.IRasterFilter#process(int, int) |
|
| 369 |
*/ |
|
| 370 |
public void process(int x, int y) {
|
|
| 371 |
|
|
| 372 |
} |
|
| 373 |
|
|
| 374 |
/* (non-Javadoc) |
|
| 375 |
* @see org.cresques.io.raster.IRasterFilter#getInRasterDataType() |
|
| 376 |
*/ |
|
| 377 |
public int getInRasterDataType() {
|
|
| 378 |
return RasterBuf.TYPE_IMAGE; |
|
| 379 |
} |
|
| 380 |
|
|
| 381 |
/* (non-Javadoc) |
|
| 382 |
* @see org.cresques.io.raster.IRasterFilter#getOutRasterDataType() |
|
| 383 |
*/ |
|
| 384 |
public int getOutRasterDataType() {
|
|
| 385 |
return RasterBuf.TYPE_IMAGE; |
|
| 386 |
} |
|
| 387 |
|
|
| 388 |
/* (non-Javadoc) |
|
| 389 |
* @see org.cresques.io.raster.IRasterFilter#getResult(java.lang.String) |
|
| 390 |
*/ |
|
| 391 |
public Object getResult(String name) {
|
|
| 392 |
if (name.equals("raster")) {
|
|
| 393 |
return (Object) this.imagePancr; |
|
| 394 |
} else {
|
|
| 395 |
return null; |
|
| 396 |
} |
|
| 397 |
} |
|
| 398 |
|
|
| 399 |
/* (non-Javadoc) |
|
| 400 |
* @see org.cresques.io.raster.RasterFilter#processLine(int) |
|
| 401 |
*/ |
|
| 402 |
public void processLine(int y) {
|
|
| 403 |
} |
|
| 404 |
|
|
| 405 |
/* (non-Javadoc) |
|
| 406 |
* @see org.cresques.io.raster.IRasterFilter#post() |
|
| 407 |
*/ |
|
| 408 |
public void post() {
|
|
| 409 |
} |
|
| 410 |
} |
|
| 411 |
|
|
| 0 | 412 | |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/raster/RasterFilterStackManager.java | ||
|---|---|---|
| 24 | 24 |
package org.cresques.io.raster; |
| 25 | 25 |
|
| 26 | 26 |
import java.lang.reflect.Constructor; |
| 27 |
|
|
| 28 | 27 |
import java.util.ArrayList; |
| 29 | 28 |
import java.util.Hashtable; |
| 30 | 29 |
import java.util.StringTokenizer; |
| 31 | 30 |
import java.util.regex.Matcher; |
| 32 | 31 |
import java.util.regex.Pattern; |
| 33 | 32 |
|
| 33 |
import org.cresques.io.GeoRasterFile; |
|
| 34 | 34 |
|
| 35 |
|
|
| 35 | 36 |
/** |
| 36 | 37 |
* Esta clase es de la parte cliente y es la encargada de la gesti?n |
| 37 | 38 |
* de la pila de filtros. Es la que conoce el orden en que se deben apilar |
| ... | ... | |
| 59 | 60 |
typeFilters.put("computeminmax", new Integer(2));
|
| 60 | 61 |
typeFilters.put("tail", new Integer(3));
|
| 61 | 62 |
typeFilters.put("removebands", new Integer(4));
|
| 63 |
typeFilters.put("sharpening", new Integer(5));
|
|
| 62 | 64 |
init(); |
| 63 | 65 |
this.filterStack.setOrder(order); |
| 64 | 66 |
} |
| ... | ... | |
| 77 | 79 |
*/ |
| 78 | 80 |
protected void init() {
|
| 79 | 81 |
order = new int[typeFilters.size()]; |
| 80 |
order[0] = ((Integer) typeFilters.get("computeminmax")).intValue();
|
|
| 81 |
order[1] = ((Integer) typeFilters.get("tail")).intValue();
|
|
| 82 |
order[2] = ((Integer) typeFilters.get("enhanced")).intValue();
|
|
| 83 |
order[3] = ((Integer) typeFilters.get("transparency")).intValue();
|
|
| 84 |
order[4] = ((Integer) typeFilters.get("removebands")).intValue();
|
|
| 82 |
order[0] = ((Integer) typeFilters.get("sharpening")).intValue();
|
|
| 83 |
order[1] = ((Integer) typeFilters.get("computeminmax")).intValue();
|
|
| 84 |
order[2] = ((Integer) typeFilters.get("tail")).intValue();
|
|
| 85 |
order[3] = ((Integer) typeFilters.get("enhanced")).intValue();
|
|
| 86 |
order[4] = ((Integer) typeFilters.get("transparency")).intValue();
|
|
| 87 |
order[5] = ((Integer) typeFilters.get("removebands")).intValue();
|
|
| 85 | 88 |
} |
| 86 | 89 |
|
| 87 | 90 |
/** |
| ... | ... | |
| 134 | 137 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("transparency")).intValue())) {
|
| 135 | 138 |
case RasterBuf.TYPE_IMAGE: |
| 136 | 139 |
filtro = new TransparencyImageFilter(); |
| 137 |
|
|
| 138 | 140 |
break; |
| 139 |
|
|
| 140 | 141 |
case RasterBuf.TYPE_SHORT: |
| 141 | 142 |
case RasterBuf.TYPE_USHORT: |
| 142 | 143 |
case RasterBuf.TYPE_INT: |
| 143 | 144 |
filtro = new TransparencyShortFilter(); |
| 144 |
|
|
| 145 | 145 |
break; |
| 146 | 146 |
} |
| 147 | 147 |
|
| ... | ... | |
| 250 | 250 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("removebands")).intValue())) {
|
| 251 | 251 |
case RasterBuf.TYPE_IMAGE: |
| 252 | 252 |
filtro = new RemoveBandsImageFilter(); |
| 253 |
|
|
| 254 | 253 |
break; |
| 255 |
|
|
| 256 | 254 |
case RasterBuf.TYPE_SHORT: |
| 257 | 255 |
case RasterBuf.TYPE_USHORT: |
| 258 | 256 |
case RasterBuf.TYPE_INT: |
| 259 | 257 |
filtro = new RemoveBandsShortFilter(); |
| 260 |
|
|
| 261 | 258 |
break; |
| 262 | 259 |
} |
| 263 | 260 |
|
| ... | ... | |
| 286 | 283 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("tail")).intValue())) {
|
| 287 | 284 |
case RasterBuf.TYPE_IMAGE: |
| 288 | 285 |
filtro = new PercentTailTrimImageFilter(); |
| 289 |
|
|
| 290 | 286 |
break; |
| 291 |
|
|
| 292 | 287 |
case RasterBuf.TYPE_SHORT: |
| 293 | 288 |
case RasterBuf.TYPE_USHORT: |
| 294 | 289 |
case RasterBuf.TYPE_INT: |
| 295 | 290 |
filtro = new PercentTailTrimShortFilter(); |
| 296 |
|
|
| 297 | 291 |
break; |
| 298 | 292 |
} |
| 299 | 293 |
|
| ... | ... | |
| 323 | 317 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("enhanced")).intValue())) {
|
| 324 | 318 |
case RasterBuf.TYPE_IMAGE: |
| 325 | 319 |
filtro = new LinearEnhancementImageFilter(); |
| 326 |
|
|
| 327 | 320 |
break; |
| 328 |
|
|
| 329 | 321 |
case RasterBuf.TYPE_SHORT: |
| 330 | 322 |
case RasterBuf.TYPE_USHORT: |
| 331 | 323 |
case RasterBuf.TYPE_INT: |
| 332 | 324 |
filtro = new LinearEnhancementShortFilter(); |
| 333 |
|
|
| 334 | 325 |
break; |
| 335 | 326 |
} |
| 336 | 327 |
|
| ... | ... | |
| 371 | 362 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("computeminmax")).intValue())) {
|
| 372 | 363 |
case RasterBuf.TYPE_IMAGE: |
| 373 | 364 |
filtro = new ComputeMinMaxImageFilter(); |
| 374 |
|
|
| 375 | 365 |
break; |
| 376 |
|
|
| 377 | 366 |
case RasterBuf.TYPE_SHORT: |
| 378 | 367 |
case RasterBuf.TYPE_USHORT: |
| 379 | 368 |
case RasterBuf.TYPE_INT: |
| 380 | 369 |
filtro = new ComputeMinMaxShortFilter(); |
| 381 |
|
|
| 382 | 370 |
break; |
| 383 | 371 |
} |
| 384 | 372 |
|
| ... | ... | |
| 391 | 379 |
} |
| 392 | 380 |
|
| 393 | 381 |
/** |
| 382 |
* A?ade un filtro ComputeMinMax |
|
| 383 |
* @param files Ficheros que componen la imagen |
|
| 384 |
* @param i Posici?n de la imagen de refinado dentro de la lista de GeoRasterFile |
|
| 385 |
* @param order Orden de visualizado de las bandas RGB |
|
| 386 |
* @param alpha Alpha aplicado a toda la imagen |
|
| 387 |
* @param method M?todo de calculo del refinado |
|
| 388 |
* @param coef Coeficiente |
|
| 389 |
*/ |
|
| 390 |
public void addSharpeningFilter(GeoRasterFile[] files, int i, int[] order, int alpha, String method, double coef, int coefBrovey) {
|
|
| 391 |
//Si ya hay un filtro de sharpening nos lo cargamos |
|
| 392 |
if (filterStack.isActive(((Integer) typeFilters.get("sharpening")).intValue())) {
|
|
| 393 |
filterStack.removeFilter(((Integer) typeFilters.get("sharpening")).intValue());
|
|
| 394 |
} |
|
| 395 |
|
|
| 396 |
RasterFilter filtro = null; |
|
| 397 |
|
|
| 398 |
switch (filterStack.getDataTypeInFilter(((Integer) typeFilters.get("sharpening")).intValue())) {
|
|
| 399 |
case RasterBuf.TYPE_IMAGE: |
|
| 400 |
filtro = new SharpeningImageFilter(); |
|
| 401 |
break; |
|
| 402 |
case RasterBuf.TYPE_SHORT: |
|
| 403 |
case RasterBuf.TYPE_USHORT: |
|
| 404 |
case RasterBuf.TYPE_INT: |
|
| 405 |
//filtro = new SharpeningShortFilter(); |
|
| 406 |
break; |
|
| 407 |
} |
|
| 408 |
|
|
| 409 |
filtro.addParam("stats", filterStack.getStats());
|
|
| 410 |
filtro.addParam("pancromatica", new Integer(i));
|
|
| 411 |
filtro.addParam("files", files);
|
|
| 412 |
filtro.addParam("order", order);
|
|
| 413 |
filtro.addParam("alpha", new Integer(alpha));
|
|
| 414 |
filtro.addParam("method", method);
|
|
| 415 |
filtro.addParam("coef", new Double(coef));
|
|
| 416 |
filtro.addParam("coefBrovey", new Integer(coefBrovey));
|
|
| 417 |
filterStack.addFilter(((Integer) typeFilters.get("sharpening")).intValue(),filtro);
|
|
| 418 |
|
|
| 419 |
this.controlTypes(); |
|
| 420 |
} |
|
| 421 |
|
|
| 422 |
/** |
|
| 423 |
* Obtiene de la pila el primer filtro del tipo solicitado si existe sino devuelve null |
|
| 424 |
* @param type Tipo de filtro |
|
| 425 |
* @return Filtro |
|
| 426 |
*/ |
|
| 427 |
public RasterFilter getFilter(String type){
|
|
| 428 |
int intType = getTypeFilter(type); |
|
| 429 |
return filterStack.getByType(intType); |
|
| 430 |
|
|
| 431 |
} |
|
| 432 |
|
|
| 433 |
/** |
|
| 394 | 434 |
* Obtiene el tipo de filtro a partir del objeto RasterFilter |
| 395 | 435 |
* @param rasterFilter Objeto RasterFilter del cual se quiere saber que tipo de filtro contiene |
| 396 | 436 |
* @return Tipo de filtro seg?n las constantes contenidas en RasterFilterStackManager |
| ... | ... | |
| 431 | 471 |
if ((i - 1) >= 0) {
|
| 432 | 472 |
//Para el primer filtro comprobamos con el tipo de dato de entrada a la pila |
| 433 | 473 |
if (i == filterStack.lenght()) {
|
| 434 |
if (filterStack.getInitDataType() != filterStack.get(i - 1) |
|
| 435 |
.getInRasterDataType()) {
|
|
| 436 |
String oldClass = filterStack.get(i - 1).getClass() |
|
| 437 |
.toString().substring(filterStack.get(i - |
|
| 438 |
1) |
|
| 439 |
.getClass() |
|
| 440 |
.toString() |
|
| 441 |
.lastIndexOf(".") +
|
|
| 442 |
1, |
|
| 443 |
filterStack.get(i - |
|
| 444 |
1) |
|
| 445 |
.getClass() |
|
| 446 |
.toString() |
|
| 447 |
.length()); |
|
| 448 |
Pattern p = Pattern.compile(RasterBuf.typesToString(filterStack.get(i - |
|
| 449 |
1) |
|
| 450 |
.getInRasterDataType())); |
|
| 474 |
if (filterStack.getInitDataType() != |
|
| 475 |
filterStack.get(i - 1).getInRasterDataType()) {
|
|
| 476 |
String oldClass = filterStack.get(i - 1). |
|
| 477 |
getClass(). |
|
| 478 |
toString(). |
|
| 479 |
substring(filterStack.get(i - 1). |
|
| 480 |
getClass(). |
|
| 481 |
toString(). |
|
| 482 |
lastIndexOf(".") +1,
|
|
| 483 |
filterStack.get(i - 1). |
|
| 484 |
getClass(). |
|
| 485 |
toString(). |
|
| 486 |
length()); |
|
| 487 |
Pattern p = Pattern.compile(RasterBuf.typesToString(filterStack.get(i - 1).getInRasterDataType())); |
|
| 451 | 488 |
Matcher m = p.matcher(oldClass); |
| 452 | 489 |
String newClass = m.replaceAll(RasterBuf.typesToString(filterStack.getInitDataType())); |
| 453 | 490 |
|
| 454 |
//System.out.println("==>"+oldClass+" "+newClass);
|
|
| 455 | 491 |
try {
|
| 456 | 492 |
Class filterClass = Class.forName("org.cresques.io.raster." +
|
| 457 | 493 |
newClass); |
| ... | ... | |
| 466 | 502 |
} |
| 467 | 503 |
|
| 468 | 504 |
//Desde el filtro 2 en adelante se compara la salida de uno con la entrada del siguiente |
| 469 |
} else if (filterStack.get(i).getOutRasterDataType() != filterStack.get(i - |
|
| 470 |
1) |
|
| 471 |
.getInRasterDataType()) {
|
|
| 472 |
String oldClass = filterStack.get(i - 1).getClass() |
|
| 473 |
.toString().substring(filterStack.get(i - |
|
| 474 |
1) |
|
| 475 |
.getClass() |
|
| 476 |
.toString() |
|
| 477 |
.lastIndexOf(".") +
|
|
| 478 |
1, |
|
| 479 |
filterStack.get(i - |
|
| 480 |
1) |
|
| 481 |
.getClass() |
|
| 482 |
.toString() |
|
| 483 |
.length()); |
|
| 484 |
Pattern p = Pattern.compile(RasterBuf.typesToString(filterStack.get(i - |
|
| 485 |
1) |
|
| 486 |
.getInRasterDataType())); |
|
| 505 |
} else if (filterStack.get(i).getOutRasterDataType() != |
|
| 506 |
filterStack.get(i - 1).getInRasterDataType()) {
|
|
| 507 |
String oldClass = filterStack.get(i - 1).getClass(). |
|
| 508 |
toString().substring(filterStack.get(i - 1). |
|
| 509 |
getClass(). |
|
| 510 |
toString(). |
|
| 511 |
lastIndexOf(".") + 1,
|
|
| 512 |
filterStack.get(i - 1). |
|
| 513 |
getClass(). |
|
| 514 |
toString(). |
|
| 515 |
length()); |
|
| 516 |
Pattern p = Pattern.compile(RasterBuf.typesToString(filterStack.get(i - 1).getInRasterDataType())); |
|
| 487 | 517 |
Matcher m = p.matcher(oldClass); |
| 488 |
String newClass = m.replaceAll(RasterBuf.typesToString(filterStack.get(i) |
|
| 489 |
.getOutRasterDataType())); |
|
| 518 |
String newClass = m.replaceAll(RasterBuf.typesToString(filterStack.get(i).getOutRasterDataType())); |
|
| 490 | 519 |
|
| 491 |
//System.out.println("==>"+oldClass+" "+newClass);
|
|
| 492 | 520 |
try {
|
| 493 | 521 |
Class filterClass = Class.forName("org.cresques.io.raster." +
|
| 494 | 522 |
newClass); |
| ... | ... | |
| 660 | 688 |
getValue(fil).equals("true")) {
|
| 661 | 689 |
filters.remove(filteri); |
| 662 | 690 |
|
| 663 |
boolean remove = false; |
|
| 664 |
|
|
| 665 | 691 |
for (int propFilter = 0; propFilter < filters.size(); |
| 666 | 692 |
propFilter++) {
|
| 667 | 693 |
String elem = (String) filters.get(propFilter); |
| ... | ... | |
| 711 | 737 |
getValue(fil).equals("true")) {
|
| 712 | 738 |
filters.remove(filteri); |
| 713 | 739 |
|
| 714 |
boolean remove = false; |
|
| 715 |
|
|
| 716 | 740 |
for (int propFilter = 0; propFilter < filters.size(); |
| 717 | 741 |
propFilter++) {
|
| 718 | 742 |
String elem = (String) filters.get(propFilter); |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/io/GdalWriter.java | ||
|---|---|---|
| 23 | 23 |
*/ |
| 24 | 24 |
package org.cresques.io; |
| 25 | 25 |
|
| 26 |
import java.io.BufferedOutputStream; |
|
| 27 |
import java.io.DataOutputStream; |
|
| 28 |
import java.io.File; |
|
| 29 |
import java.io.FileOutputStream; |
|
| 30 |
import java.io.IOException; |
|
| 31 |
|
|
| 32 |
import org.cresques.cts.IProjection; |
|
| 33 |
import org.cresques.px.PxRaster; |
|
| 34 |
|
|
| 26 | 35 |
import es.gva.cit.jecwcompress.EcwException; |
| 27 | 36 |
import es.gva.cit.jgdal.Gdal; |
| 28 | 37 |
import es.gva.cit.jgdal.GdalBuffer; |
| ... | ... | |
| 32 | 41 |
import es.gva.cit.jgdal.GeoTransform; |
| 33 | 42 |
import es.gva.cit.jogr.OGRSpatialReference; |
| 34 | 43 |
|
| 35 |
import org.cresques.px.PxRaster; |
|
| 36 | 44 |
|
| 37 |
import java.io.IOException; |
|
| 38 |
|
|
| 39 |
|
|
| 40 | 45 |
/** |
| 41 | 46 |
* Driver para la escritura a trav?s de Gdal. |
| 42 | 47 |
* Puede exportar un fichero de un formato a otro desde un GeoRasterFile |
| ... | ... | |
| 55 | 60 |
public class GdalWriter extends GeoRasterWriter {
|
| 56 | 61 |
static {
|
| 57 | 62 |
GeoRasterWriter.registerWriterExtension("tif", GdalWriter.class);
|
| 63 |
GeoRasterWriter.registerWriterExtension("jpg", GdalWriter.class);
|
|
| 64 |
GeoRasterWriter.registerWriterExtension("jpeg", GdalWriter.class);
|
|
| 58 | 65 |
} |
| 59 | 66 |
|
| 60 | 67 |
public final int windowSizeX = 386; |
| ... | ... | |
| 66 | 73 |
private Gdal dset_destino = null; |
| 67 | 74 |
private GdalRasterBand rband = null; |
| 68 | 75 |
private String drvType = null; //Tipo de driver |
| 69 |
private String[] supportedDrv = { "GTiff" }; //Lista de drivers de escritura que soporta
|
|
| 76 |
private String[] supportedDrv = { "GTiff", "JPEG" }; //Lista de drivers de escritura que soporta
|
|
| 70 | 77 |
private GeoTransform geot = null; //Datos de georeferenciaci?n |
| 71 | 78 |
private OGRSpatialReference oSRS; //Datos de proyecci?n |
| 72 | 79 |
private GdalBuffer buf = null; //Buffer de origen de gdal |
| ... | ... | |
| 86 | 93 |
public GdalWriter(String drvType) {
|
| 87 | 94 |
this.support = new GdalSupportOptions(drvType); |
| 88 | 95 |
this.ident = this.drvType = drvType; |
| 89 |
this.driver = "tif"; |
|
| 96 |
if(drvType.equals("GTiff"))
|
|
| 97 |
this.driver = "tif"; |
|
| 98 |
else if(drvType.equals("JPEG"))
|
|
| 99 |
this.driver = "jpg"; |
|
| 90 | 100 |
this.support.setBlockSize(64); |
| 91 | 101 |
this.support.setPhotometric("RGB");
|
| 92 | 102 |
this.support.setInterleave("BAND");
|
| ... | ... | |
| 108 | 118 |
String drvType) throws GdalException, IOException {
|
| 109 | 119 |
this.support = new GdalSupportOptions(drvType); |
| 110 | 120 |
this.ident = this.drvType = drvType; |
| 111 |
this.driver = "tif"; |
|
| 112 |
|
|
| 121 |
if(drvType.equals("GTiff"))
|
|
| 122 |
this.driver = "tif"; |
|
| 123 |
else if(drvType.equals("JPEG"))
|
|
| 124 |
this.driver = "jpg"; |
|
| 125 |
|
|
| 113 | 126 |
this.outfilename = outfilename; |
| 114 | 127 |
this.infilename = infilename; |
| 115 | 128 |
this.currentRaster = raster; |
| ... | ... | |
| 134 | 147 |
|
| 135 | 148 |
nBands = currentRaster.getBandCount(); |
| 136 | 149 |
|
| 137 |
this.support.setBlockSize(currentRaster.getBlockSize());
|
|
| 150 |
this.support.setBlockSize(64/*currentRaster.getBlockSize()*/);
|
|
| 138 | 151 |
|
| 139 | 152 |
if ((sizeWindowX < 0) || (sizeWindowY < 0)) {
|
| 140 | 153 |
throw new IOException("Tama?o del fichero de salida erroneo.");
|
| ... | ... | |
| 187 | 200 |
throws GdalException, IOException {
|
| 188 | 201 |
this.support = new GdalSupportOptions(drvType); |
| 189 | 202 |
this.ident = this.drvType = drvType; |
| 190 |
this.driver = "tif"; |
|
| 203 |
if(drvType.equals("GTiff"))
|
|
| 204 |
this.driver = "tif"; |
|
| 205 |
else if(drvType.equals("JPEG"))
|
|
| 206 |
this.driver = "jpg"; |
|
| 191 | 207 |
|
| 192 | 208 |
this.dataWriter = dataWriter; |
| 193 | 209 |
this.outfilename = outFilename; |
| ... | ... | |
| 286 | 302 |
} |
| 287 | 303 |
|
| 288 | 304 |
//Obtenemos el driver y creamos el dataset del destino |
| 305 |
drv = Gdal.getDriverByName(drvType); |
|
| 289 | 306 |
|
| 290 |
drv = Gdal.getDriverByName(drvType); |
|
| 291 |
|
|
| 307 |
/*try{
|
|
| 308 |
drv = Gdal.getDriverByName("JPEG");
|
|
| 309 |
Gdal src = new Gdal(); |
|
| 310 |
src.open("/home/nacho/imagenes/q101866.tif", 0);
|
|
| 311 |
drv.createCopy( "/home/nacho/outgvGIG.jpg", src, false ); |
|
| 312 |
}catch(IOException e){}
|
|
| 313 |
return;*/ |
|
| 314 |
|
|
| 292 | 315 |
//System.out.println("++++++>"+params.length);
|
| 293 | 316 |
//for(int i=0;i<params.length;i++)System.out.println("=====>"+params[i]);
|
| 317 |
|
|
| 294 | 318 |
if (dset_destino != null) {
|
| 295 | 319 |
dset_destino.close(); |
| 296 | 320 |
dset_destino = null; |
| ... | ... | |
| 513 | 537 |
} |
| 514 | 538 |
|
| 515 | 539 |
/** |
| 540 |
* Realiza una copia en el formato especificado. |
|
| 541 |
* @throws IOException |
|
| 542 |
*/ |
|
| 543 |
public static void createCopy(GdalDriver driverDst, String dst, String src, boolean bstrict, String[] params, IProjection proj) throws IOException, GdalException {
|
|
| 544 |
if (dst == null || src == null) {
|
|
| 545 |
throw new IOException("No se ha asignado un fichero de entrada.");
|
|
| 546 |
} |
|
| 547 |
|
|
| 548 |
GdalFile gdalFile = new GdalFile(proj, src); |
|
| 549 |
driverDst.createCopy(dst, gdalFile.file, bstrict, params); |
|
| 550 |
if(dst.endsWith(".jpg") || dst.endsWith(".jpeg"))
|
|
| 551 |
GdalWriter.createWorldFile(dst, gdalFile); |
|
| 552 |
gdalFile.close(); |
|
| 553 |
} |
|
| 554 |
|
|
| 555 |
/** |
|
| 556 |
* Crea un fichero de georeferenciaci?n |
|
| 557 |
* @param img |
|
| 558 |
* @param name |
|
| 559 |
* @return |
|
| 560 |
* @throws IOException |
|
| 561 |
*/ |
|
| 562 |
private static void createWorldFile(String name, GdalFile gdalFile) throws IOException{
|
|
| 563 |
File tfw = null; |
|
| 564 |
|
|
| 565 |
String extWorldFile = ".wld"; |
|
| 566 |
if(name.endsWith("tif"))
|
|
| 567 |
extWorldFile = ".tfw"; |
|
| 568 |
if(name.endsWith("jpg") || name.endsWith("jpeg"))
|
|
| 569 |
extWorldFile = ".jpgw"; |
|
| 570 |
|
|
| 571 |
tfw = new File(name.substring(0, name.lastIndexOf(".")) + extWorldFile);
|
|
| 572 |
|
|
| 573 |
//Generamos un world file para gdal |
|
| 574 |
DataOutputStream dos = new DataOutputStream( new BufferedOutputStream(new FileOutputStream(tfw)) ); |
|
| 575 |
dos.writeBytes((gdalFile.getExtent().getMax().getX() - gdalFile.getExtent().getMin().getX())/gdalFile.getWidth()+"\n"); |
|
| 576 |
dos.writeBytes("0.0\n");
|
|
| 577 |
dos.writeBytes("0.0\n");
|
|
| 578 |
dos.writeBytes((gdalFile.getExtent().getMax().getY() - gdalFile.getExtent().getMin().getY())/gdalFile.getHeight()+"\n"); |
|
| 579 |
dos.writeBytes(""+gdalFile.getExtent().getMin().getX()+"\n");
|
|
| 580 |
dos.writeBytes(""+gdalFile.getExtent().getMin().getY()+"\n");
|
|
| 581 |
dos.close(); |
|
| 582 |
} |
|
| 583 |
|
|
| 584 |
|
|
| 585 |
/** |
|
| 516 | 586 |
* Realiza la escritura de datos con los datos que le pasa el cliente. |
| 517 | 587 |
* @throws IOException |
| 518 | 588 |
*/ |
| trunk/libraries/libCq CMS for java.old/src/org/cresques/px/PxRaster.java | ||
|---|---|---|
| 213 | 213 |
geoFile = listFiles; |
| 214 | 214 |
} else {
|
| 215 | 215 |
System.err.println("PxRaster.addFile(): Imagen no cargada.");
|
| 216 |
|
|
| 217 | 216 |
return; |
| 218 | 217 |
} |
| 219 | 218 |
} |
| ... | ... | |
| 778 | 777 |
stackManager.addTailFilter(this.percentFilterInit, 0D, false); |
| 779 | 778 |
} |
Also available in: Unified diff