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 Ch11_Circle_Ellipse_Fitting;
010
011
012import ij.ImagePlus;
013import ij.gui.GenericDialog;
014import ij.gui.NewImage;
015import ij.gui.PointRoi;
016import ij.plugin.PlugIn;
017import ij.process.ImageProcessor;
018import imagingbook.common.color.sets.BasicAwtColor;
019import imagingbook.common.geometry.basic.Pnt2d;
020import imagingbook.common.geometry.circle.GeometricCircle;
021import imagingbook.common.geometry.fitting.circle.utils.CircleSampler;
022import imagingbook.common.ij.DialogUtils;
023import imagingbook.common.ij.DialogUtils.DialogLabel;
024import imagingbook.common.ij.DialogUtils.DialogStringColumns;
025import imagingbook.common.ij.RoiUtils;
026import imagingbook.common.ij.overlay.ColoredStroke;
027import imagingbook.common.ij.overlay.ShapeOverlayAdapter;
028import imagingbook.common.util.ParameterBundle;
029import imagingbook.core.jdoc.JavaDocHelp;
030
031import static imagingbook.common.ij.DialogUtils.addToDialog;
032
033/**
034 * Samples points on a given (ideal) circle and creates a new image with the
035 * sample points marked and also contained in a {@link PointRoi}. Image size,
036 * circle parameters and noise can be specified. The result can be used as a
037 * test image for circle fitting. Note that the resulting image has an inverted
038 * LUT, i.e., the background value is 0 and marked points have value 255.
039 * 
040 * @author WB
041 * @version 2022/10/03
042 */
043public class Circle_Make_Random implements PlugIn, JavaDocHelp {
044        
045        public static class Parameters implements ParameterBundle<Circle_Make_Random> {
046                
047                @DialogLabel("image title")@DialogStringColumns(12)
048                public String Title = "RandomCircle"; // Circle_Make_Random.class.getSimpleName();
049                
050                @DialogLabel("image width")
051                public int W = 400;
052                @DialogLabel("image height")
053                public int H = 400;
054                
055                @DialogLabel("number of points")
056                public int n = 20;
057                
058                @DialogLabel("circle center (xc)")
059                public double xc = 200;
060                
061                @DialogLabel("circle center (yc)")
062                public double yc = 190;
063                
064                @DialogLabel("circle radius (r)")
065                public double r = 150;
066                
067                @DialogLabel("start angle (deg)")
068                public double angle0 = 0;
069                
070                @DialogLabel("stop angle (deg)")
071                public double angle1 = 180;
072                
073                @DialogLabel("x/y noise (sigma)")
074                public double sigma = 5.0; //2.0;
075                @DialogLabel("Random seed (0 = none)")
076                public int seed = 0;
077                
078                @DialogLabel("show real circle")
079                public boolean ShowRealCurve = true;
080                @DialogLabel("circle color")
081                public BasicAwtColor StrokeColor = BasicAwtColor.Green;
082                @DialogLabel("stroke width")
083                public double StrokeWidth = 1.0;
084        };
085        
086        private static Parameters params = new Parameters();
087        
088        @Override
089        public void run(String arg) {
090                if (!runDialog()) {
091                        return;
092                }
093                
094                GeometricCircle realCircle = new GeometricCircle(params.xc, params.yc, params.r);
095                CircleSampler sampler = new CircleSampler(realCircle, params.seed);
096                Pnt2d[] points =
097                                sampler.getPoints(params.n, Math.toRadians(params.angle0), Math.toRadians(params.angle1), params.sigma);
098
099                ImagePlus im = NewImage.createByteImage(params.Title, params.W, params.H, 1, NewImage.FILL_BLACK);
100                im.setRoi(RoiUtils.toPointRoi(points));
101                
102                ImageProcessor ip = im.getProcessor();
103                for (Pnt2d p : points) {
104                        int u = (int) Math.rint(p.getX());
105                        int v = (int) Math.rint(p.getY());
106                        ip.putPixel(u, v, 255);
107                }
108                ip.invertLut();
109                
110                if (params.ShowRealCurve) {
111                        ShapeOverlayAdapter ola = new ShapeOverlayAdapter();
112                        ColoredStroke circleStroke = new ColoredStroke(params.StrokeWidth, params.StrokeColor.getColor());
113                        ola.addShapes(realCircle.getShapes(3), circleStroke);
114                        im.setOverlay(ola.getOverlay());
115                }
116                
117                im.show();
118        }
119        
120        // ------------------------------------------
121        
122        private boolean runDialog() {
123                GenericDialog gd = new GenericDialog(this.getClass().getSimpleName());
124                gd.addHelp(getJavaDocUrl());
125                gd.addMessage(DialogUtils.formatText(50,
126                                "This plugin samples points on a given (ideal) circle and",
127                                "creates a new image with the sample points marked and also",
128                                "contained in a ROI (float coordinates)."
129                                ));
130                
131                addToDialog(params, gd);
132
133                gd.showDialog();
134                if (gd.wasCanceled())
135                        return false;
136
137                return params.validate();
138        }
139
140}