001/**
002 * Copyright (c) 2011, The University of Southampton and the individual contributors.
003 * All rights reserved.
004 *
005 * Redistribution and use in source and binary forms, with or without modification,
006 * are permitted provided that the following conditions are met:
007 *
008 *   *  Redistributions of source code must retain the above copyright notice,
009 *      this list of conditions and the following disclaimer.
010 *
011 *   *  Redistributions in binary form must reproduce the above copyright notice,
012 *      this list of conditions and the following disclaimer in the documentation
013 *      and/or other materials provided with the distribution.
014 *
015 *   *  Neither the name of the University of Southampton nor the names of its
016 *      contributors may be used to endorse or promote products derived from this
017 *      software without specific prior written permission.
018 *
019 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
020 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
021 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
022 * DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR
023 * ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
024 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
025 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
026 * ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
027 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
028 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
029 */
030package org.openimaj.workinprogress;
031
032import java.io.File;
033import java.io.IOException;
034
035import org.openimaj.image.DisplayUtilities;
036import org.openimaj.image.FImage;
037import org.openimaj.image.ImageUtilities;
038import org.openimaj.image.MBFImage;
039import org.openimaj.image.processing.algorithm.FilterSupport;
040import org.openimaj.image.processing.algorithm.MedianFilter;
041import org.openimaj.image.processing.resize.ResizeProcessor;
042import org.openimaj.video.xuggle.XuggleVideo;
043
044public class BackgroundSubtractor {
045        public static void main(String[] args) throws IOException {
046                final XuggleVideo xv = new XuggleVideo(new File("/Users/jon/Desktop/merlin/tunnel.mp4"));
047                final FImage bg = ResizeProcessor.halfSize(
048                                ImageUtilities.readF(new File("/Users/jon/Desktop/merlin/tunnel-background.png"))
049                                );
050
051                // final XuggleVideoWriter xvw = new
052                // XuggleVideoWriter("/Users/jon/Desktop/merlin/tunnel-proc.mp4",
053                // bg.width,
054                // bg.height, xv.getFPS());
055                for (final MBFImage frc : xv) {
056                        final FImage fr = ResizeProcessor.halfSize(frc.flatten());
057                        final MBFImage diff = diff(bg, fr);
058
059                        // xvw.addFrame(diff);
060                        DisplayUtilities.displayName(diff, "");
061                }
062                // xvw.close();
063        }
064
065        static MBFImage diff(FImage bg, FImage fg) {
066                final FImage df = new FImage(bg.getWidth(), bg.getHeight());
067                final float[][] dff = df.pixels;
068
069                final float[][] bgfr = bg.pixels;
070                final float[][] fgfr = fg.pixels;
071
072                for (int y = 0; y < df.getHeight(); y++) {
073                        for (int x = 0; x < df.getWidth(); x++) {
074                                final float dr = bgfr[y][x] - fgfr[y][x];
075                                final float ssd = dr * dr;
076
077                                if (ssd < 0.03) {
078                                        dff[y][x] = 0;
079                                } else {
080                                        dff[y][x] = 1;
081                                }
082                        }
083                }
084
085                // Dilate.dilate(df, 1);
086                // Erode.erode(df, 2);
087                df.processInplace(new MedianFilter(FilterSupport.createBlockSupport(3, 3)));
088                df.processInplace(new MedianFilter(FilterSupport.createBlockSupport(3, 3)));
089
090                return df.toRGB();
091        }
092}