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.workinprogress;
31
32 import java.io.File;
33 import java.io.IOException;
34
35 import org.openimaj.image.DisplayUtilities;
36 import org.openimaj.image.FImage;
37 import org.openimaj.image.ImageUtilities;
38 import org.openimaj.image.MBFImage;
39 import org.openimaj.image.processing.algorithm.FilterSupport;
40 import org.openimaj.image.processing.algorithm.MedianFilter;
41 import org.openimaj.image.processing.resize.ResizeProcessor;
42 import org.openimaj.video.xuggle.XuggleVideo;
43
44 public class BackgroundSubtractor {
45 public static void main(String[] args) throws IOException {
46 final XuggleVideo xv = new XuggleVideo(new File("/Users/jon/Desktop/merlin/tunnel.mp4"));
47 final FImage bg = ResizeProcessor.halfSize(
48 ImageUtilities.readF(new File("/Users/jon/Desktop/merlin/tunnel-background.png"))
49 );
50
51 // final XuggleVideoWriter xvw = new
52 // XuggleVideoWriter("/Users/jon/Desktop/merlin/tunnel-proc.mp4",
53 // bg.width,
54 // bg.height, xv.getFPS());
55 for (final MBFImage frc : xv) {
56 final FImage fr = ResizeProcessor.halfSize(frc.flatten());
57 final MBFImage diff = diff(bg, fr);
58
59 // xvw.addFrame(diff);
60 DisplayUtilities.displayName(diff, "");
61 }
62 // xvw.close();
63 }
64
65 static MBFImage diff(FImage bg, FImage fg) {
66 final FImage df = new FImage(bg.getWidth(), bg.getHeight());
67 final float[][] dff = df.pixels;
68
69 final float[][] bgfr = bg.pixels;
70 final float[][] fgfr = fg.pixels;
71
72 for (int y = 0; y < df.getHeight(); y++) {
73 for (int x = 0; x < df.getWidth(); x++) {
74 final float dr = bgfr[y][x] - fgfr[y][x];
75 final float ssd = dr * dr;
76
77 if (ssd < 0.03) {
78 dff[y][x] = 0;
79 } else {
80 dff[y][x] = 1;
81 }
82 }
83 }
84
85 // Dilate.dilate(df, 1);
86 // Erode.erode(df, 2);
87 df.processInplace(new MedianFilter(FilterSupport.createBlockSupport(3, 3)));
88 df.processInplace(new MedianFilter(FilterSupport.createBlockSupport(3, 3)));
89
90 return df.toRGB();
91 }
92 }