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.demos.sandbox; 031 032import java.io.IOException; 033 034import org.openimaj.image.DisplayUtilities; 035import org.openimaj.image.FImage; 036import org.openimaj.image.ImageUtilities; 037import org.openimaj.image.MBFImage; 038import org.openimaj.image.colour.RGBColour; 039import org.openimaj.image.model.asm.ActiveShapeModel.IterationResult; 040import org.openimaj.image.model.asm.MultiResolutionActiveShapeModel; 041import org.openimaj.image.model.asm.datasets.ShapeModelDataset; 042import org.openimaj.image.model.asm.datasets.ShapeModelDatasets; 043import org.openimaj.image.model.landmark.FNormalLandmarkModel; 044import org.openimaj.image.pixel.sampling.FLineSampler; 045import org.openimaj.math.geometry.line.Line2d; 046import org.openimaj.math.geometry.point.Point2d; 047import org.openimaj.math.geometry.point.PointList; 048import org.openimaj.math.geometry.point.PointListConnections; 049import org.openimaj.math.geometry.shape.PointDistributionModel; 050import org.openimaj.math.geometry.transforms.TransformUtilities; 051import org.openimaj.math.matrix.algorithm.pca.PrincipalComponentAnalysis.NumberComponentSelector; 052 053import Jama.Matrix; 054 055public class ASMPlayground { 056 /** 057 * @param args 058 * @throws IOException 059 */ 060 public static void main(String[] args) throws IOException { 061 final ShapeModelDataset<FImage> dataset = ShapeModelDatasets.loadASFDataset("/Users/jsh2/Downloads/imm_face_db", 062 ImageUtilities.FIMAGE_READER); 063 final PointListConnections conns = dataset.getConnections(); 064 065 final float scale = 0.02f; 066 final FNormalLandmarkModel.Factory factory = new FNormalLandmarkModel.Factory(conns, 067 FLineSampler.INTERPOLATED_DERIVATIVE, 5, 9, scale); 068 final MultiResolutionActiveShapeModel<FImage> asm = MultiResolutionActiveShapeModel.trainModel(3, 069 new NumberComponentSelector(19), dataset, new PointDistributionModel.BoxConstraint(3), factory); 070 071 final Matrix pose = TransformUtilities.translateMatrix(300, 300).times(TransformUtilities.scaleMatrix(70, 70)); 072 PointList shape = asm.getPDM().getMean().transform(pose); 073 // PointList shape = ASFDataset.readASF(new File(dir, 074 // "01-1m.asf")).firstObject(); 075 final FImage img = dataset.get(0).secondObject(); 076 // PointList shape = ASFDataset.readASF(new File(dir, 077 // "16-6m.asf")).firstObject(); 078 // FImage img = ASFDataset.readASF(new File(dir, 079 // "16-6m.asf")).secondObject(); 080 081 final MBFImage image = img.toRGB(); 082 image.drawLines(conns.getLines(shape), 1, RGBColour.RED); 083 084 final long t1 = System.currentTimeMillis(); 085 final IterationResult newData = asm.fit(img, shape); 086 final long t2 = System.currentTimeMillis(); 087 088 shape = newData.shape; 089 090 System.out.println(newData.fit); 091 System.out.println(t2 - t1); 092 093 image.drawLines(conns.getLines(shape), 1, RGBColour.GREEN); 094 095 final float shapeScale = shape.computeIntrinsicScale(); 096 for (final Point2d pt : shape) { 097 final Line2d normal = conns.calculateNormalLine(pt, shape, scale * shapeScale); 098 if (normal != null) 099 image.drawLine(normal, 1, RGBColour.BLUE); 100 } 101 102 DisplayUtilities.display(image); 103 } 104}