1 /**
2 * Copyright (c) 2011, The University of Southampton and the individual contributors.
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without modification,
6 * are permitted provided that the following conditions are met:
7 *
8 * * Redistributions of source code must retain the above copyright notice,
9 * this list of conditions and the following disclaimer.
10 *
11 * * Redistributions in binary form must reproduce the above copyright notice,
12 * this list of conditions and the following disclaimer in the documentation
13 * and/or other materials provided with the distribution.
14 *
15 * * Neither the name of the University of Southampton nor the names of its
16 * contributors may be used to endorse or promote products derived from this
17 * software without specific prior written permission.
18 *
19 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
20 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
21 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
22 * DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR
23 * ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
24 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
25 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
26 * ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
27 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
28 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
29 */
30 package org.openimaj.image.processing.transform;
31
32 import org.openimaj.image.FImage;
33 import org.openimaj.image.MBFImage;
34 import org.openimaj.math.geometry.point.Point2dImpl;
35 import org.openimaj.math.geometry.shape.Rectangle;
36 import org.openimaj.math.geometry.shape.Shape;
37
38 /**
39 * @author Jonathon Hare (jsh2@ecs.soton.ac.uk)
40 * @author Sina Samangooei (ss@ecs.soton.ac.uk)
41 *
42 * Perform a set of matrix transforms on a set of images and construct a single image containing all the pixels (or a window of the pixels)
43 * in the projected space.
44 */
45 public class MBFProjectionProcessor extends ProjectionProcessor<Float[], MBFImage> {
46 @Override
47 public MBFImage performProjection(int windowMinC , int windowMinR , MBFImage output) {
48 final FImage[] bands = new FImage[output.numBands()];
49 for(int i = 0; i < bands.length; i++){
50 bands[i] = output.getBand(i);
51 }
52 final FImage[][] input = new FImage[this.projectedShapes.size()][];
53 for(int i = 0; i < input.length; i++){
54 MBFImage inputMBF = this.images.get(i);
55 input[i] = new FImage[inputMBF.numBands()];
56 for(int j = 0; j < input[i].length; j++){
57 input[i][j] = inputMBF.getBand(j);
58 }
59 }
60 for(int y = 0; y < output.getHeight(); y++)
61 {
62 for(int x = 0; x < output.getWidth(); x++){
63 Point2dImpl realPoint = new Point2dImpl(windowMinC + x,windowMinR + y);
64 int i = 0;
65 for(int k = 0; k < this.projectedRectangles.size(); k++){
66 Rectangle r = this.projectedRectangles.get(k);
67 Shape s = this.projectedShapes.get(k);
68 if( r.isInside(realPoint) && s.isInside(realPoint)){
69 double[][] transform = this.transformsInverted.get(i).getArray();
70
71 float xt = (float)transform[0][0] * realPoint.x + (float)transform[0][1] * realPoint.y + (float)transform[0][2];
72 float yt = (float)transform[1][0] * realPoint.x + (float)transform[1][1] * realPoint.y + (float)transform[1][2];
73 float zt = (float)transform[2][0] * realPoint.x + (float)transform[2][1] * realPoint.y + (float)transform[2][2];
74
75 xt /= zt;
76 yt /= zt;
77
78 for(int j = 0; j < bands.length; j++){
79 FImage in = input[i][j];
80 FImage out = bands[j];
81 out.pixels[y][x] = in.getPixelInterpNative(xt, yt,out.pixels[y][x]);
82 }
83
84 }
85 i++;
86 }
87 }
88 }
89 return output;
90 }
91 }