001/*******************************************************************************
002 * This software is provided as a supplement to the authors' textbooks on digital
003 * image processing published by Springer-Verlag in various languages and editions.
004 * Permission to use and distribute this software is granted under the BSD 2-Clause
005 * "Simplified" License (see http://opensource.org/licenses/BSD-2-Clause).
006 * Copyright (c) 2006-2023 Wilhelm Burger, Mark J. Burge. All rights reserved.
007 * Visit https://imagingbook.com for additional details.
008 ******************************************************************************/
009package Ch21_Geometric_Operations;
010
011import ij.ImagePlus;
012import ij.plugin.filter.PlugInFilter;
013import ij.process.ImageProcessor;
014import imagingbook.common.geometry.basic.Pnt2d;
015import imagingbook.common.geometry.mappings.Mapping2D;
016import imagingbook.common.ij.DialogUtils;
017import imagingbook.common.ij.IjUtils;
018import imagingbook.common.image.ImageMapper;
019import imagingbook.common.image.interpolation.InterpolationMethod;
020import imagingbook.core.jdoc.JavaDocHelp;
021import imagingbook.sampleimages.GeneralSampleImage;
022
023/**
024 * <p>
025 * ImageJ plugin, applies a non-linear "tapestry" transformation to the current image. See Sec. 21.1.7 (Exercise 21.9
026 * and Fig. 21.18) of [1] for details. Optionally opens a sample image if no image is currently open.
027 * </p>
028 * <p>
029 * [1] W. Burger, M.J. Burge, <em>Digital Image Processing &ndash; An Algorithmic Introduction</em>, 3rd ed, Springer
030 * (2022).
031 * </p>
032 *
033 * @author WB
034 * @version 2022/11/28
035 * @see ImageMapper
036 * @see Mapping2D
037 */
038public class Map_Nonlinear_Tapestry implements PlugInFilter, JavaDocHelp {
039
040        private static final double a = 5.0;
041        private static final double tx = 30;
042        private static final double ty = 30;
043        
044                
045        /**
046         * Constructor, asks to open a predefined sample image if no other image
047         * is currently open.
048         */
049        public Map_Nonlinear_Tapestry() {
050                if (IjUtils.noCurrentImage()) {
051                        DialogUtils.askForSampleImage(GeneralSampleImage.Flower);
052                }
053        }
054        
055        @Override
056        public int setup(String arg, ImagePlus imp) {
057                return DOES_ALL;
058        }
059
060        @Override
061        public void run(ImageProcessor ip) {
062                
063                final double xc = 0.5 * ip.getWidth();
064                final double yc = 0.5 * ip.getHeight();
065                
066                Mapping2D imap = new Mapping2D() {
067                        final double bx = 2 * Math.PI / tx;
068                        final double by = 2 * Math.PI / ty;
069                        
070                        @Override
071                        public Pnt2d applyTo(Pnt2d pnt) {
072                                double x = pnt.getX();
073                                double y = pnt.getY();
074                                double xx = x + a * Math.sin((x - xc) * bx);    // yes, both are sin()!
075                                double yy = y + a * Math.sin((y - yc) * by);
076                                return Pnt2d.from(xx,  yy);
077                        }
078                };
079                
080                new ImageMapper(imap, null, InterpolationMethod.Bicubic).map(ip);
081        }
082
083}