root / trunk / libraries / libRaster / src-test / org / gvsig / fmap / dataaccess / TestGdalClassHistogram.java @ 10740
History | View | Annotate | Download (5.59 KB)
1 |
/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana
|
---|---|
2 |
*
|
3 |
* Copyright (C) 2007 IVER T.I. 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 |
package org.gvsig.fmap.dataaccess; |
20 |
|
21 |
import junit.framework.TestCase; |
22 |
|
23 |
import org.gvsig.raster.dataaccess.DataSource; |
24 |
import org.gvsig.raster.driver.FileNotOpenException; |
25 |
import org.gvsig.raster.driver.IBuffer; |
26 |
import org.gvsig.raster.driver.NotSupportedExtensionException; |
27 |
import org.gvsig.raster.driver.RasterDataset; |
28 |
import org.gvsig.raster.driver.RasterDriverException; |
29 |
import org.gvsig.raster.shared.DataClass; |
30 |
import org.gvsig.raster.shared.DataClassList; |
31 |
import org.gvsig.raster.shared.RasterLibrary; |
32 |
|
33 |
/**
|
34 |
* Test a un histograma usando clases. El test se realiza sobre un raster en
|
35 |
* coma flotante de 28x25 pixels. Se crean clases de tres tipo de intervalos distintos y
|
36 |
* se testean los valores obtenidos en el histograma resultante. El ?ltimo
|
37 |
* test comprueba la generaci?n autom?tica de clases.
|
38 |
*
|
39 |
* @author Nacho Brodin (nachobrodin@gmail.com)
|
40 |
*/
|
41 |
public class TestGdalClassHistogram extends TestCase{ |
42 |
|
43 |
private String baseDir = "./test-images/"; |
44 |
private String path = baseDir + "miniRaster28x25F32.tif"; |
45 |
private RasterDataset f = null; |
46 |
private DataSource ds = null; |
47 |
|
48 |
public void setUp() { |
49 |
System.out.println("TestGdalClassHistogram running..."); |
50 |
} |
51 |
|
52 |
static{
|
53 |
RasterLibrary.wakeUp(); |
54 |
} |
55 |
|
56 |
public void testStack(){ |
57 |
int[] drawableBands = {0, 1, 2}; |
58 |
try {
|
59 |
f = RasterDataset.openFile(null, path);
|
60 |
} catch (NotSupportedExtensionException e) {
|
61 |
return;
|
62 |
} catch (RasterDriverException e) {
|
63 |
return;
|
64 |
} |
65 |
ds = new DataSource(f); |
66 |
ds.addDrawableBands(drawableBands); |
67 |
ds.setAreaOfInterest(0, 0, 28, 25); |
68 |
//printData();
|
69 |
DataClassList classList = new DataClassList();
|
70 |
DataClass classInterval1 = new DataClass(0D, 1000D, 0, ""); |
71 |
DataClass classInterval2 = new DataClass(1000D, 2000D, 1, ""); |
72 |
classList.addClass(classInterval1); |
73 |
classList.addClass(classInterval2); |
74 |
long[][] histogram = null; |
75 |
try {
|
76 |
histogram = ds.getHistogram(classList); |
77 |
} catch (FileNotOpenException e) {
|
78 |
e.printStackTrace(); |
79 |
} catch (RasterDriverException e) {
|
80 |
e.printStackTrace(); |
81 |
} |
82 |
testHist1(histogram); |
83 |
|
84 |
classList.clear(); |
85 |
classInterval1 = new DataClass(0D, 1200D, 0, ""); |
86 |
classInterval2 = new DataClass(1200D, 2000D, 1, ""); |
87 |
classList.addClass(classInterval1); |
88 |
classList.addClass(classInterval2); |
89 |
try {
|
90 |
histogram = ds.getHistogram(classList); |
91 |
} catch (FileNotOpenException e) {
|
92 |
e.printStackTrace(); |
93 |
} catch (RasterDriverException e) {
|
94 |
e.printStackTrace(); |
95 |
} |
96 |
testHist2(histogram); |
97 |
|
98 |
classList.clear(); |
99 |
classInterval1 = new DataClass(0D, 1150D, 0, ""); |
100 |
classInterval2 = new DataClass(115D, 1175D, 1, ""); |
101 |
DataClass classInterval3 = new DataClass(1175D, 1200D, 2, ""); |
102 |
DataClass classInterval4 = new DataClass(1200D, 1210D, 3, ""); |
103 |
DataClass classInterval5 = new DataClass(1210D, 2000D, 4, ""); |
104 |
classList.addClass(classInterval1); |
105 |
classList.addClass(classInterval2); |
106 |
classList.addClass(classInterval3); |
107 |
classList.addClass(classInterval4); |
108 |
classList.addClass(classInterval5); |
109 |
try {
|
110 |
histogram = ds.getHistogram(classList); |
111 |
} catch (FileNotOpenException e) {
|
112 |
e.printStackTrace(); |
113 |
} catch (RasterDriverException e) {
|
114 |
e.printStackTrace(); |
115 |
} |
116 |
testHist3(histogram); |
117 |
|
118 |
try {
|
119 |
histogram = ds.getHistogram(null);
|
120 |
} catch (FileNotOpenException e) {
|
121 |
e.printStackTrace(); |
122 |
} catch (RasterDriverException e) {
|
123 |
e.printStackTrace(); |
124 |
} |
125 |
testHist4(histogram); |
126 |
//print(histogram);
|
127 |
} |
128 |
|
129 |
private void testHist1(long[][] histogram){ |
130 |
assertEquals(histogram[0][0], 0); |
131 |
assertEquals(histogram[0][1], 700); |
132 |
} |
133 |
|
134 |
private void testHist2(long[][] histogram){ |
135 |
assertEquals(histogram[0][0], 180); |
136 |
assertEquals(histogram[0][1], 520); |
137 |
} |
138 |
|
139 |
private void testHist3(long[][] histogram){ |
140 |
assertEquals(histogram[0][0], 4); |
141 |
assertEquals(histogram[0][1], 38); |
142 |
assertEquals(histogram[0][2], 138); |
143 |
assertEquals(histogram[0][3], 100); |
144 |
assertEquals(histogram[0][4], 420); |
145 |
} |
146 |
|
147 |
private void testHist4(long[][] histogram){ |
148 |
assertEquals(histogram[0][0], 10); |
149 |
assertEquals(histogram[0][1], 17); |
150 |
assertEquals(histogram[0][2], 36); |
151 |
assertEquals(histogram[0][3], 72); |
152 |
assertEquals(histogram[0][4], 119); |
153 |
assertEquals(histogram[0][5], 129); |
154 |
assertEquals(histogram[0][6], 98); |
155 |
assertEquals(histogram[0][7], 93); |
156 |
assertEquals(histogram[0][8], 76); |
157 |
assertEquals(histogram[0][9], 49); |
158 |
} |
159 |
|
160 |
private void print(long[][] histogram){ |
161 |
for (int i = 0; i < histogram.length; i++) { |
162 |
for (int j = 0; j < histogram[i].length; j++){ |
163 |
System.out.print(histogram[i][j] + " "); |
164 |
} |
165 |
System.out.println(""); |
166 |
} |
167 |
} |
168 |
|
169 |
/**
|
170 |
* Imprime todos los pixels de la fuente de datos en RGB
|
171 |
*/
|
172 |
private void printData(){ |
173 |
IBuffer raster = ds.getRasterBuf(); |
174 |
for(int line = 0; line < raster.getHeight(); line++){ |
175 |
for(int col = 0; col < raster.getWidth(); col++) |
176 |
System.out.print("(" + raster.getElemFloat(line, col, 0) + ")"); |
177 |
System.out.println();
|
178 |
} |
179 |
} |
180 |
} |