?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
Use internal SNAPSHOT couplings again
[trakem2.git] / TrakEM2_ / src / main / java / ini / trakem2 / utils / MCCube.java
blobe94c6b664906da6b7be6567b0541ef3693807a08
1 package ini.trakem2.utils;
2
3 import ij.IJ;
4 import ij3d.Volume;
5
6 import java.awt.Polygon;
7 import java.awt.Rectangle;
8 import java.awt.Shape;
9 import java.awt.geom.Area;
10 import java.awt.geom.PathIterator;
11 import java.util.ArrayList;
12 import java.util.HashMap;
13 import java.util.List;
14
15 import javax.vecmath.Point3f;
16
17 import mpicbg.imglib.container.shapelist.ShapeList;
18
19 public final class MCCube {
20 // vertexes
21 private Point3f[] v;
22
23 // interpolated values
24 private Point3f[] e;
25
26 private MCCube() {
27 this.v = new Point3f[8];
28 for(int i = 0; i < 8; i++)
29 v[i] = new Point3f();
30 this.e = new Point3f[12];
31 for(int i = 0; i < 12; i++)
32 e[i] = new Point3f();
33 }
34
35 /**
36 * initializes a MCCube object
37 * _________ 0______x
38 * /v0 v1/| /|
39 * /________/ | / |
40 * |v3 v2| /v5 y/ |z
41 * |________|/
42 * v7 v6
43 */
44 public void init(int x, int y, int z){
45 v[0].set(x, y, z);
46 v[1].set(x + 1, y, z);
47 v[2].set(x + 1, y + 1, z);
48 v[3].set(x, y + 1, z);
49 v[4].set(x, y, z + 1);
50 v[5].set(x + 1, y, z + 1);
51 v[6].set(x + 1, y + 1, z + 1);
52 v[7].set(x, y + 1, z + 1);
53 }
54
55 /**
56 * computes the interpolated point along a specified whose
57 * intensity equals the reference value
58 * @param v1 first extremity of the edge
59 * @param v2 second extremity of the edge
60 * @param result stores the resulting edge
61 * return the point on the edge where intensity equals the isovalue;
62 * @return false if the interpolated point is beyond edge boundaries
63 */
64 private boolean computeEdge(Point3f v1, int i1,
65 Point3f v2, int i2,
66 Point3f result, final Carrier car) {
67
68 // 30 --- 50 --- 70 : t=0.5
69 // 70 --- 50 --- 30 : t=0.5
70 ///int i1 = car.intensity(v1);
71 ///int i2 = car.intensity(v2);
72 if(i2 < i1)
73 return computeEdge(v2, i2, v1, i1, result, car);
74
75 float t = (car.threshold - i1) / (float) (i2 - i1);
76 if (t >= 0 && t <= 1) {
77 // v1 + t*(v2-v1)
78 result.set(v2);
79 result.sub(v1);
80 result.scale(t);
81 result.add(v1);
82 return true;
83 }
84 result.set(-1, -1, -1);
85 return false;
86 }
87
88 /**
89 * computes interpolated values along each edge of the cube
90 * (null if interpolated value doesn't belong to the edge)
91 */
92 private void computeEdges(final Carrier car) {
93 int i0 = car.intensity(v[0]);
94 int i1 = car.intensity(v[1]);
95 int i2 = car.intensity(v[2]);
96 int i3 = car.intensity(v[3]);
97 int i4 = car.intensity(v[4]);
98 int i5 = car.intensity(v[5]);
99 int i6 = car.intensity(v[6]);
100 int i7 = car.intensity(v[7]);
101
102 this.computeEdge(v[0], i0, v[1], i1, e[0], car);
103 this.computeEdge(v[1], i1, v[2], i2, e[1], car);
104 this.computeEdge(v[2], i2, v[3], i3, e[2], car);
105 this.computeEdge(v[3], i3, v[0], i0, e[3], car);
106
107 this.computeEdge(v[4], i4, v[5], i5, e[4], car);
108 this.computeEdge(v[5], i5, v[6], i6, e[5], car);
109 this.computeEdge(v[6], i6, v[7], i7, e[6], car);
110 this.computeEdge(v[7], i7, v[4], i4, e[7], car);
111
112 this.computeEdge(v[0], i0, v[4], i4, e[8], car);
113 this.computeEdge(v[1], i1, v[5], i5, e[9], car);
114 this.computeEdge(v[3], i3, v[7], i7, e[10], car);
115 this.computeEdge(v[2], i2, v[6], i6, e[11], car);
116 }
117
118 /**
119 * indicates if a number corresponds to an ambigous case
120 * @param n number of the case to test
121 * @return true if the case if ambigous
122 */
123 private static boolean isAmbigous(int n) {
124 boolean result = false;
125 for (int index = 0; index < MCCube.ambigous.length; index++) {
126 result |= MCCube.ambigous[index] == n;
127 }
128 return result;
129 }
130
131 private void getTriangles(List<Point3f> list, final Carrier car){
132 int cn = caseNumber(car);
133 boolean directTable = !(isAmbigous(cn));
134 directTable = true;
135
136 // address in the table
137 int offset = directTable ? cn*15 : (255-cn)*15;
138 for(int index = 0; index < 5; index++){
139 // if there's a triangle
140 if (faces[offset] != -1) {
141 // pick up vertexes of the current triangle
142 list.add(new Point3f(this.e[faces[offset+0]]));
143 list.add(new Point3f(this.e[faces[offset+1]]));
144 list.add(new Point3f(this.e[faces[offset+2]]));
145 }
146 offset += 3;
147 }
148 }
149
150 /**
151 * computes the case number of the cube
152 * @return the number of the case corresponding to the cube
153 */
154 private int caseNumber(final Carrier car) {
155 int caseNumber = 0;
156 for (int index = -1;
157 ++index < v.length;
158 caseNumber +=
159 (car.intensity(v[index]) - car.threshold > 0)
160 ? 1 << index
161 : 0);
162 return caseNumber;
163 }
164
165 /**
166 * An encapsulating class to avoid thread collisions on static fields.
167 */
168 private static final class Carrier {
169 int w, h, d;
170 Volume volume;
171 float threshold;
172
173 final int intensity(final Point3f p) {
174 if(p.x < 0 || p.y < 0 || p.z < 0
175 || p.x >= w || p.y >= h || p.z >= d)
176 return 0;
177 return volume.load((int)p.x, (int)p.y, (int)p.z);
178 }
179 }
180
181 /**
182 * Create a list of triangles from the specified image data and the
183 * given isovalue.
184 * @param volume
185 * @param thresh
186 * @return
187 */
188 public static final List<Point3f> getTriangles(Volume volume, int thresh){
189 List<Point3f> tri = new ArrayList<Point3f>();
190 final Carrier car = new Carrier();
191 car.w = volume.xDim;
192 car.h = volume.yDim;
193 car.d = volume.zDim;
194 car.threshold = thresh + 0.5f;
195 car.volume = volume;
196
197 if (volume instanceof ImgLibVolume && ((ImgLibVolume)volume).getImage().getContainer() instanceof ShapeList) {
198 getShapeListImageTriangles((ImgLibVolume)volume, car, tri);
199 } else {
200 MCCube cube = new MCCube();
201 for(int z = -1; z < car.d+1; z+=1){
202 for(int x = -1; x < car.w+1; x+=1){
203 for(int y = -1; y < car.h+1; y+=1){
204 cube.init(x, y, z);
205 cube.computeEdges(car);
206 cube.getTriangles(tri, car);
207 }
208 }
209 IJ.showProgress(z, car.d-2);
210 }
211 }
212
213 // convert pixel coordinates
214 for(int i = 0; i < tri.size(); i++) {
215 Point3f p = (Point3f)tri.get(i);
216 p.x = (float) (p.x * volume.pw + volume.minCoord.x);
217 p.y = (float) (p.y * volume.ph + volume.minCoord.y);
218 p.z = (float) (p.z * volume.pd + volume.minCoord.z);
219 }
220 return tri;
221 }
222
223 /** Identical to getTriangles, but iterates only the minimal necessary bounding box, by asking the shapes objects. */
224 private static final void getShapeListImageTriangles(final ImgLibVolume volume, final Carrier car, final List<Point3f> tri) {
225 final ShapeList sli = (ShapeList) volume.getImage().getContainer();
226 final ArrayList<ArrayList<Shape>> shapeLists = sli.getShapeLists();
227 final Area[] sectionAreas = new Area[shapeLists.size()];
228 // Create one Area for each section, composed of the addition of all Shape instances
229 {
230 int next = -1;
231 for (final ArrayList<Shape> shapeList : shapeLists) {
232 next++;
233 if (shapeList.isEmpty()) {
234 continue;
235 }
236 final Area a = new Area(shapeList.get(0));
237 for (int i=1; i<shapeList.size(); i++) {
238 a.add(new Area(shapeList.get(i)));
239 }
240 sectionAreas[next] = a;
241 }
242 }
243 // Fuse Area instances for previous and next sections
244 final Area[] scanAreas = new Area[sectionAreas.length];
245 for (int i=0; i<sectionAreas.length; i++) {
246 if (null == sectionAreas[i]) continue;
247 final Area a = new Area(sectionAreas[i]);
248 if (i-1 < 0 || null == sectionAreas[i-1]) {}
249 else a.add(sectionAreas[i-1]);
250 if (i+1 > sectionAreas.length -1 || null == sectionAreas[i+1]) {}
251 else a.add(sectionAreas[i+1]);
252 scanAreas[i] = a;
253 }
254 // Collect the bounds of all subareas in each scanArea:
255 final HashMap<Integer,ArrayList<Rectangle>> sectionBounds = new HashMap<Integer,ArrayList<Rectangle>>();
256 for (int i=0; i<scanAreas.length; i++) {
257 if (null == scanAreas[i]) continue;
258 final ArrayList<Rectangle> bs = new ArrayList<Rectangle>();
259 Polygon pol = new Polygon();
260 final float[] coords = new float[6];
261 for (final PathIterator pit = scanAreas[i].getPathIterator(null); !pit.isDone(); pit.next()) {
262 switch (pit.currentSegment(coords)) {
263 case PathIterator.SEG_MOVETO:
264 case PathIterator.SEG_LINETO:
265 pol.addPoint((int)coords[0], (int)coords[1]);
266 break;
267 case PathIterator.SEG_CLOSE:
268 bs.add(pol.getBounds());
269 pol = new Polygon();
270 break;
271 default:
272 System.out.println("WARNING: unhandled seg type.");
273 break;
274 }
275 }
276 sectionBounds.put(i, bs);
277 }
278
279 // Add Z paddings on top and bottom
280 sectionBounds.put(-1, sectionBounds.get(0));
281 sectionBounds.put(car.d, sectionBounds.get(car.d-1));
282
283 // Scan only relevant areas:
284 final MCCube cube = new MCCube();
285 for (int z = -1; z < car.d + 1; z += 1) {
286 final ArrayList<Rectangle> bs = sectionBounds.get(z);
287 if (null == bs || bs.isEmpty()) continue;
288 for (final Rectangle bounds : bs) {
289 for (int x = bounds.x -1; x < bounds.x + bounds.width +2; x+=1) {
290 for (int y = bounds.y -1; y < bounds.y + bounds.height +2; y+=1) {
291 cube.init(x, y, z);
292 cube.computeEdges(car);
293 cube.getTriangles(tri, car);
294 }
295 }
296 }
297
298 IJ.showProgress(z, car.d-2);
299 }
300 }
301
302 protected static final int ambigous[] = {
303 250,
304 245,
305 237,
306 231,
307 222,
308 219,
309 189,
310 183,
311 175,
312 126,
313 123,
314 95,
315 234,
316 233,
317 227,
318 214,
319 213,
320 211,
321 203,
322 199,
323 188,
324 186,
325 182,
326 174,
327 171,
328 158,
329 151,
330 124,
331 121,
332 117,
333 109,
334 107,
335 93,
336 87,
337 62,
338 61,
339 229,
340 218,
341 181,
342 173,
343 167,
344 122,
345 94,
346 91,
347 150,
348 170,
349 195,
350 135,
351 149,
352 154,
353 163,
354 166,
355 169,
356 172,
357 180,
358 197,
359 202,
360 210,
361 225,
362 165
363 };
364
365 // triangles to be drawn in each case
366 private static final int faces[] =
367 {
368 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
369 0, 8, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
370 0, 1, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
371 1, 8, 3, 9, 8, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
372 1, 2, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
373 0, 8, 3, 1, 2, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1,
374 9, 2, 11, 0, 2, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1,
375 2, 8, 3, 2, 11, 8, 11, 9, 8, -1, -1, -1, -1, -1, -1,
376 3, 10, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
377 0, 10, 2, 8, 10, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1,
378 1, 9, 0, 2, 3, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
379 1, 10, 2, 1, 9, 10, 9, 8, 10, -1, -1, -1, -1, -1, -1,
380 3, 11, 1, 10, 11, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1,
381 0, 11, 1, 0, 8, 11, 8, 10, 11, -1, -1, -1, -1, -1, -1,
382 3, 9, 0, 3, 10, 9, 10, 11, 9, -1, -1, -1, -1, -1, -1,
383 9, 8, 11, 11, 8, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
384 4, 7, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
385 4, 3, 0, 7, 3, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1,
386 0, 1, 9, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
387 4, 1, 9, 4, 7, 1, 7, 3, 1, -1, -1, -1, -1, -1, -1,
388 1, 2, 11, 8, 4, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
389 3, 4, 7, 3, 0, 4, 1, 2, 11, -1, -1, -1, -1, -1, -1,
390 9, 2, 11, 9, 0, 2, 8, 4, 7, -1, -1, -1, -1, -1, -1,
391 2, 11, 9, 2, 9, 7, 2, 7, 3, 7, 9, 4, -1, -1, -1,
392 8, 4, 7, 3, 10, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1,
393 10, 4, 7, 10, 2, 4, 2, 0, 4, -1, -1, -1, -1, -1, -1,
394 9, 0, 1, 8, 4, 7, 2, 3, 10, -1, -1, -1, -1, -1, -1,
395 4, 7, 10, 9, 4, 10, 9, 10, 2, 9, 2, 1, -1, -1, -1,
396 3, 11, 1, 3, 10, 11, 7, 8, 4, -1, -1, -1, -1, -1, -1,
397 1, 10, 11, 1, 4, 10, 1, 0, 4, 7, 10, 4, -1, -1, -1,
398 4, 7, 8, 9, 0, 10, 9, 10, 11, 10, 0, 3, -1, -1, -1,
399 4, 7, 10, 4, 10, 9, 9, 10, 11, -1, -1, -1, -1, -1, -1,
400 9, 5, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
401 9, 5, 4, 0, 8, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1,
402 0, 5, 4, 1, 5, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1,
403 8, 5, 4, 8, 3, 5, 3, 1, 5, -1, -1, -1, -1, -1, -1,
404 1, 2, 11, 9, 5, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1,
405 3, 0, 8, 1, 2, 11, 4, 9, 5, -1, -1, -1, -1, -1, -1,
406 5, 2, 11, 5, 4, 2, 4, 0, 2, -1, -1, -1, -1, -1, -1,
407 2, 11, 5, 3, 2, 5, 3, 5, 4, 3, 4, 8, -1, -1, -1,
408 9, 5, 4, 2, 3, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
409 0, 10, 2, 0, 8, 10, 4, 9, 5, -1, -1, -1, -1, -1, -1,
410 0, 5, 4, 0, 1, 5, 2, 3, 10, -1, -1, -1, -1, -1, -1,
411 2, 1, 5, 2, 5, 8, 2, 8, 10, 4, 8, 5, -1, -1, -1,
412 11, 3, 10, 11, 1, 3, 9, 5, 4, -1, -1, -1, -1, -1, -1,
413 4, 9, 5, 0, 8, 1, 8, 11, 1, 8, 10, 11, -1, -1, -1,
414 5, 4, 0, 5, 0, 10, 5, 10, 11, 10, 0, 3, -1, -1, -1,
415 5, 4, 8, 5, 8, 11, 11, 8, 10, -1, -1, -1, -1, -1, -1,
416 9, 7, 8, 5, 7, 9, -1, -1, -1, -1, -1, -1, -1, -1, -1,
417 9, 3, 0, 9, 5, 3, 5, 7, 3, -1, -1, -1, -1, -1, -1,
418 0, 7, 8, 0, 1, 7, 1, 5, 7, -1, -1, -1, -1, -1, -1,
419 1, 5, 3, 3, 5, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
420 9, 7, 8, 9, 5, 7, 11, 1, 2, -1, -1, -1, -1, -1, -1,
421 11, 1, 2, 9, 5, 0, 5, 3, 0, 5, 7, 3, -1, -1, -1,
422 8, 0, 2, 8, 2, 5, 8, 5, 7, 11, 5, 2, -1, -1, -1,
423 2, 11, 5, 2, 5, 3, 3, 5, 7, -1, -1, -1, -1, -1, -1,
424 7, 9, 5, 7, 8, 9, 3, 10, 2, -1, -1, -1, -1, -1, -1,
425 9, 5, 7, 9, 7, 2, 9, 2, 0, 2, 7, 10, -1, -1, -1,
426 2, 3, 10, 0, 1, 8, 1, 7, 8, 1, 5, 7, -1, -1, -1,
427 10, 2, 1, 10, 1, 7, 7, 1, 5, -1, -1, -1, -1, -1, -1,
428 9, 5, 8, 8, 5, 7, 11, 1, 3, 11, 3, 10, -1, -1, -1,
429 5, 7, 0, 5, 0, 9, 7, 10, 0, 1, 0, 11, 10, 11, 0,
430 10, 11, 0, 10, 0, 3, 11, 5, 0, 8, 0, 7, 5, 7, 0,
431 10, 11, 5, 7, 10, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1,
432 11, 6, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
433 0, 8, 3, 5, 11, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
434 9, 0, 1, 5, 11, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
435 1, 8, 3, 1, 9, 8, 5, 11, 6, -1, -1, -1, -1, -1, -1,
436 1, 6, 5, 2, 6, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
437 1, 6, 5, 1, 2, 6, 3, 0, 8, -1, -1, -1, -1, -1, -1,
438 9, 6, 5, 9, 0, 6, 0, 2, 6, -1, -1, -1, -1, -1, -1,
439 5, 9, 8, 5, 8, 2, 5, 2, 6, 3, 2, 8, -1, -1, -1,
440 2, 3, 10, 11, 6, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1,
441 10, 0, 8, 10, 2, 0, 11, 6, 5, -1, -1, -1, -1, -1, -1,
442 0, 1, 9, 2, 3, 10, 5, 11, 6, -1, -1, -1, -1, -1, -1,
443 5, 11, 6, 1, 9, 2, 9, 10, 2, 9, 8, 10, -1, -1, -1,
444 6, 3, 10, 6, 5, 3, 5, 1, 3, -1, -1, -1, -1, -1, -1,
445 0, 8, 10, 0, 10, 5, 0, 5, 1, 5, 10, 6, -1, -1, -1,
446 3, 10, 6, 0, 3, 6, 0, 6, 5, 0, 5, 9, -1, -1, -1,
447 6, 5, 9, 6, 9, 10, 10, 9, 8, -1, -1, -1, -1, -1, -1,
448 5, 11, 6, 4, 7, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1,
449 4, 3, 0, 4, 7, 3, 6, 5, 11, -1, -1, -1, -1, -1, -1,
450 1, 9, 0, 5, 11, 6, 8, 4, 7, -1, -1, -1, -1, -1, -1,
451 11, 6, 5, 1, 9, 7, 1, 7, 3, 7, 9, 4, -1, -1, -1,
452 6, 1, 2, 6, 5, 1, 4, 7, 8, -1, -1, -1, -1, -1, -1,
453 1, 2, 5, 5, 2, 6, 3, 0, 4, 3, 4, 7, -1, -1, -1,
454 8, 4, 7, 9, 0, 5, 0, 6, 5, 0, 2, 6, -1, -1, -1,
455 7, 3, 9, 7, 9, 4, 3, 2, 9, 5, 9, 6, 2, 6, 9,
456 3, 10, 2, 7, 8, 4, 11, 6, 5, -1, -1, -1, -1, -1, -1,
457 5, 11, 6, 4, 7, 2, 4, 2, 0, 2, 7, 10, -1, -1, -1,
458 0, 1, 9, 4, 7, 8, 2, 3, 10, 5, 11, 6, -1, -1, -1,
459 9, 2, 1, 9, 10, 2, 9, 4, 10, 7, 10, 4, 5, 11, 6,
460 8, 4, 7, 3, 10, 5, 3, 5, 1, 5, 10, 6, -1, -1, -1,
461 5, 1, 10, 5, 10, 6, 1, 0, 10, 7, 10, 4, 0, 4, 10,
462 0, 5, 9, 0, 6, 5, 0, 3, 6, 10, 6, 3, 8, 4, 7,
463 6, 5, 9, 6, 9, 10, 4, 7, 9, 7, 10, 9, -1, -1, -1,
464 11, 4, 9, 6, 4, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1,
465 4, 11, 6, 4, 9, 11, 0, 8, 3, -1, -1, -1, -1, -1, -1,
466 11, 0, 1, 11, 6, 0, 6, 4, 0, -1, -1, -1, -1, -1, -1,
467 8, 3, 1, 8, 1, 6, 8, 6, 4, 6, 1, 11, -1, -1, -1,
468 1, 4, 9, 1, 2, 4, 2, 6, 4, -1, -1, -1, -1, -1, -1,
469 3, 0, 8, 1, 2, 9, 2, 4, 9, 2, 6, 4, -1, -1, -1,
470 0, 2, 4, 4, 2, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
471 8, 3, 2, 8, 2, 4, 4, 2, 6, -1, -1, -1, -1, -1, -1,
472 11, 4, 9, 11, 6, 4, 10, 2, 3, -1, -1, -1, -1, -1, -1,
473 0, 8, 2, 2, 8, 10, 4, 9, 11, 4, 11, 6, -1, -1, -1,
474 3, 10, 2, 0, 1, 6, 0, 6, 4, 6, 1, 11, -1, -1, -1,
475 6, 4, 1, 6, 1, 11, 4, 8, 1, 2, 1, 10, 8, 10, 1,
476 9, 6, 4, 9, 3, 6, 9, 1, 3, 10, 6, 3, -1, -1, -1,
477 8, 10, 1, 8, 1, 0, 10, 6, 1, 9, 1, 4, 6, 4, 1,
478 3, 10, 6, 3, 6, 0, 0, 6, 4, -1, -1, -1, -1, -1, -1,
479 6, 4, 8, 10, 6, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1,
480 7, 11, 6, 7, 8, 11, 8, 9, 11, -1, -1, -1, -1, -1, -1,
481 0, 7, 3, 0, 11, 7, 0, 9, 11, 6, 7, 11, -1, -1, -1,
482 11, 6, 7, 1, 11, 7, 1, 7, 8, 1, 8, 0, -1, -1, -1,
483 11, 6, 7, 11, 7, 1, 1, 7, 3, -1, -1, -1, -1, -1, -1,
484 1, 2, 6, 1, 6, 8, 1, 8, 9, 8, 6, 7, -1, -1, -1,
485 2, 6, 9, 2, 9, 1, 6, 7, 9, 0, 9, 3, 7, 3, 9,
486 7, 8, 0, 7, 0, 6, 6, 0, 2, -1, -1, -1, -1, -1, -1,
487 7, 3, 2, 6, 7, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1,
488 2, 3, 10, 11, 6, 8, 11, 8, 9, 8, 6, 7, -1, -1, -1,
489 2, 0, 7, 2, 7, 10, 0, 9, 7, 6, 7, 11, 9, 11, 7,
490 1, 8, 0, 1, 7, 8, 1, 11, 7, 6, 7, 11, 2, 3, 10,
491 10, 2, 1, 10, 1, 7, 11, 6, 1, 6, 7, 1, -1, -1, -1,
492 8, 9, 6, 8, 6, 7, 9, 1, 6, 10, 6, 3, 1, 3, 6,
493 0, 9, 1, 10, 6, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
494 7, 8, 0, 7, 0, 6, 3, 10, 0, 10, 6, 0, -1, -1, -1,
495 7, 10, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
496 7, 6, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
497 3, 0, 8, 10, 7, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
498 0, 1, 9, 10, 7, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
499 8, 1, 9, 8, 3, 1, 10, 7, 6, -1, -1, -1, -1, -1, -1,
500 11, 1, 2, 6, 10, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
501 1, 2, 11, 3, 0, 8, 6, 10, 7, -1, -1, -1, -1, -1, -1,
502 2, 9, 0, 2, 11, 9, 6, 10, 7, -1, -1, -1, -1, -1, -1,
503 6, 10, 7, 2, 11, 3, 11, 8, 3, 11, 9, 8, -1, -1, -1,
504 7, 2, 3, 6, 2, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
505 7, 0, 8, 7, 6, 0, 6, 2, 0, -1, -1, -1, -1, -1, -1,
506 2, 7, 6, 2, 3, 7, 0, 1, 9, -1, -1, -1, -1, -1, -1,
507 1, 6, 2, 1, 8, 6, 1, 9, 8, 8, 7, 6, -1, -1, -1,
508 11, 7, 6, 11, 1, 7, 1, 3, 7, -1, -1, -1, -1, -1, -1,
509 11, 7, 6, 1, 7, 11, 1, 8, 7, 1, 0, 8, -1, -1, -1,
510 0, 3, 7, 0, 7, 11, 0, 11, 9, 6, 11, 7, -1, -1, -1,
511 7, 6, 11, 7, 11, 8, 8, 11, 9, -1, -1, -1, -1, -1, -1,
512 6, 8, 4, 10, 8, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
513 3, 6, 10, 3, 0, 6, 0, 4, 6, -1, -1, -1, -1, -1, -1,
514 8, 6, 10, 8, 4, 6, 9, 0, 1, -1, -1, -1, -1, -1, -1,
515 9, 4, 6, 9, 6, 3, 9, 3, 1, 10, 3, 6, -1, -1, -1,
516 6, 8, 4, 6, 10, 8, 2, 11, 1, -1, -1, -1, -1, -1, -1,
517 1, 2, 11, 3, 0, 10, 0, 6, 10, 0, 4, 6, -1, -1, -1,
518 4, 10, 8, 4, 6, 10, 0, 2, 9, 2, 11, 9, -1, -1, -1,
519 11, 9, 3, 11, 3, 2, 9, 4, 3, 10, 3, 6, 4, 6, 3,
520 8, 2, 3, 8, 4, 2, 4, 6, 2, -1, -1, -1, -1, -1, -1,
521 0, 4, 2, 4, 6, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1,
522 1, 9, 0, 2, 3, 4, 2, 4, 6, 4, 3, 8, -1, -1, -1,
523 1, 9, 4, 1, 4, 2, 2, 4, 6, -1, -1, -1, -1, -1, -1,
524 8, 1, 3, 8, 6, 1, 8, 4, 6, 6, 11, 1, -1, -1, -1,
525 11, 1, 0, 11, 0, 6, 6, 0, 4, -1, -1, -1, -1, -1, -1,
526 4, 6, 3, 4, 3, 8, 6, 11, 3, 0, 3, 9, 11, 9, 3,
527 11, 9, 4, 6, 11, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1,
528 4, 9, 5, 7, 6, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
529 0, 8, 3, 4, 9, 5, 10, 7, 6, -1, -1, -1, -1, -1, -1,
530 5, 0, 1, 5, 4, 0, 7, 6, 10, -1, -1, -1, -1, -1, -1,
531 10, 7, 6, 8, 3, 4, 3, 5, 4, 3, 1, 5, -1, -1, -1,
532 9, 5, 4, 11, 1, 2, 7, 6, 10, -1, -1, -1, -1, -1, -1,
533 6, 10, 7, 1, 2, 11, 0, 8, 3, 4, 9, 5, -1, -1, -1,
534 7, 6, 10, 5, 4, 11, 4, 2, 11, 4, 0, 2, -1, -1, -1,
535 3, 4, 8, 3, 5, 4, 3, 2, 5, 11, 5, 2, 10, 7, 6,
536 7, 2, 3, 7, 6, 2, 5, 4, 9, -1, -1, -1, -1, -1, -1,
537 9, 5, 4, 0, 8, 6, 0, 6, 2, 6, 8, 7, -1, -1, -1,
538 3, 6, 2, 3, 7, 6, 1, 5, 0, 5, 4, 0, -1, -1, -1,
539 6, 2, 8, 6, 8, 7, 2, 1, 8, 4, 8, 5, 1, 5, 8,
540 9, 5, 4, 11, 1, 6, 1, 7, 6, 1, 3, 7, -1, -1, -1,
541 1, 6, 11, 1, 7, 6, 1, 0, 7, 8, 7, 0, 9, 5, 4,
542 4, 0, 11, 4, 11, 5, 0, 3, 11, 6, 11, 7, 3, 7, 11,
543 7, 6, 11, 7, 11, 8, 5, 4, 11, 4, 8, 11, -1, -1, -1,
544 6, 9, 5, 6, 10, 9, 10, 8, 9, -1, -1, -1, -1, -1, -1,
545 3, 6, 10, 0, 6, 3, 0, 5, 6, 0, 9, 5, -1, -1, -1,
546 0, 10, 8, 0, 5, 10, 0, 1, 5, 5, 6, 10, -1, -1, -1,
547 6, 10, 3, 6, 3, 5, 5, 3, 1, -1, -1, -1, -1, -1, -1,
548 1, 2, 11, 9, 5, 10, 9, 10, 8, 10, 5, 6, -1, -1, -1,
549 0, 10, 3, 0, 6, 10, 0, 9, 6, 5, 6, 9, 1, 2, 11,
550 10, 8, 5, 10, 5, 6, 8, 0, 5, 11, 5, 2, 0, 2, 5,
551 6, 10, 3, 6, 3, 5, 2, 11, 3, 11, 5, 3, -1, -1, -1,
552 5, 8, 9, 5, 2, 8, 5, 6, 2, 3, 8, 2, -1, -1, -1,
553 9, 5, 6, 9, 6, 0, 0, 6, 2, -1, -1, -1, -1, -1, -1,
554 1, 5, 8, 1, 8, 0, 5, 6, 8, 3, 8, 2, 6, 2, 8,
555 1, 5, 6, 2, 1, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1,
556 1, 3, 6, 1, 6, 11, 3, 8, 6, 5, 6, 9, 8, 9, 6,
557 11, 1, 0, 11, 0, 6, 9, 5, 0, 5, 6, 0, -1, -1, -1,
558 0, 3, 8, 5, 6, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1,
559 11, 5, 6, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
560 10, 5, 11, 7, 5, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
561 10, 5, 11, 10, 7, 5, 8, 3, 0, -1, -1, -1, -1, -1, -1,
562 5, 10, 7, 5, 11, 10, 1, 9, 0, -1, -1, -1, -1, -1, -1,
563 11, 7, 5, 11, 10, 7, 9, 8, 1, 8, 3, 1, -1, -1, -1,
564 10, 1, 2, 10, 7, 1, 7, 5, 1, -1, -1, -1, -1, -1, -1,
565 0, 8, 3, 1, 2, 7, 1, 7, 5, 7, 2, 10, -1, -1, -1,
566 9, 7, 5, 9, 2, 7, 9, 0, 2, 2, 10, 7, -1, -1, -1,
567 7, 5, 2, 7, 2, 10, 5, 9, 2, 3, 2, 8, 9, 8, 2,
568 2, 5, 11, 2, 3, 5, 3, 7, 5, -1, -1, -1, -1, -1, -1,
569 8, 2, 0, 8, 5, 2, 8, 7, 5, 11, 2, 5, -1, -1, -1,
570 9, 0, 1, 5, 11, 3, 5, 3, 7, 3, 11, 2, -1, -1, -1,
571 9, 8, 2, 9, 2, 1, 8, 7, 2, 11, 2, 5, 7, 5, 2,
572 1, 3, 5, 3, 7, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1,
573 0, 8, 7, 0, 7, 1, 1, 7, 5, -1, -1, -1, -1, -1, -1,
574 9, 0, 3, 9, 3, 5, 5, 3, 7, -1, -1, -1, -1, -1, -1,
575 9, 8, 7, 5, 9, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1,
576 5, 8, 4, 5, 11, 8, 11, 10, 8, -1, -1, -1, -1, -1, -1,
577 5, 0, 4, 5, 10, 0, 5, 11, 10, 10, 3, 0, -1, -1, -1,
578 0, 1, 9, 8, 4, 11, 8, 11, 10, 11, 4, 5, -1, -1, -1,
579 11, 10, 4, 11, 4, 5, 10, 3, 4, 9, 4, 1, 3, 1, 4,
580 2, 5, 1, 2, 8, 5, 2, 10, 8, 4, 5, 8, -1, -1, -1,
581 0, 4, 10, 0, 10, 3, 4, 5, 10, 2, 10, 1, 5, 1, 10,
582 0, 2, 5, 0, 5, 9, 2, 10, 5, 4, 5, 8, 10, 8, 5,
583 9, 4, 5, 2, 10, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1,
584 2, 5, 11, 3, 5, 2, 3, 4, 5, 3, 8, 4, -1, -1, -1,
585 5, 11, 2, 5, 2, 4, 4, 2, 0, -1, -1, -1, -1, -1, -1,
586 3, 11, 2, 3, 5, 11, 3, 8, 5, 4, 5, 8, 0, 1, 9,
587 5, 11, 2, 5, 2, 4, 1, 9, 2, 9, 4, 2, -1, -1, -1,
588 8, 4, 5, 8, 5, 3, 3, 5, 1, -1, -1, -1, -1, -1, -1,
589 0, 4, 5, 1, 0, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1,
590 8, 4, 5, 8, 5, 3, 9, 0, 5, 0, 3, 5, -1, -1, -1,
591 9, 4, 5, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
592 4, 10, 7, 4, 9, 10, 9, 11, 10, -1, -1, -1, -1, -1, -1,
593 0, 8, 3, 4, 9, 7, 9, 10, 7, 9, 11, 10, -1, -1, -1,
594 1, 11, 10, 1, 10, 4, 1, 4, 0, 7, 4, 10, -1, -1, -1,
595 3, 1, 4, 3, 4, 8, 1, 11, 4, 7, 4, 10, 11, 10, 4,
596 4, 10, 7, 9, 10, 4, 9, 2, 10, 9, 1, 2, -1, -1, -1,
597 9, 7, 4, 9, 10, 7, 9, 1, 10, 2, 10, 1, 0, 8, 3,
598 10, 7, 4, 10, 4, 2, 2, 4, 0, -1, -1, -1, -1, -1, -1,
599 10, 7, 4, 10, 4, 2, 8, 3, 4, 3, 2, 4, -1, -1, -1,
600 2, 9, 11, 2, 7, 9, 2, 3, 7, 7, 4, 9, -1, -1, -1,
601 9, 11, 7, 9, 7, 4, 11, 2, 7, 8, 7, 0, 2, 0, 7,
602 3, 7, 11, 3, 11, 2, 7, 4, 11, 1, 11, 0, 4, 0, 11,
603 1, 11, 2, 8, 7, 4, -1, -1, -1, -1, -1, -1, -1, -1, -1,
604 4, 9, 1, 4, 1, 7, 7, 1, 3, -1, -1, -1, -1, -1, -1,
605 4, 9, 1, 4, 1, 7, 0, 8, 1, 8, 7, 1, -1, -1, -1,
606 4, 0, 3, 7, 4, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1,
607 4, 8, 7, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
608 9, 11, 8, 11, 10, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1,
609 3, 0, 9, 3, 9, 10, 10, 9, 11, -1, -1, -1, -1, -1, -1,
610 0, 1, 11, 0, 11, 8, 8, 11, 10, -1, -1, -1, -1, -1, -1,
611 3, 1, 11, 10, 3, 11, -1, -1, -1, -1, -1, -1, -1, -1, -1,
612 1, 2, 10, 1, 10, 9, 9, 10, 8, -1, -1, -1, -1, -1, -1,
613 3, 0, 9, 3, 9, 10, 1, 2, 9, 2, 10, 9, -1, -1, -1,
614 0, 2, 10, 8, 0, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1,
615 3, 2, 10, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
616 2, 3, 8, 2, 8, 11, 11, 8, 9, -1, -1, -1, -1, -1, -1,
617 9, 11, 2, 0, 9, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1,
618 2, 3, 8, 2, 8, 11, 0, 1, 8, 1, 11, 8, -1, -1, -1,
619 1, 11, 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
620 1, 3, 8, 9, 1, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1,
621 0, 9, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
622 0, 3, 8, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
623 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
624 };
625 }
register/edit/analyze/3D model TB of images