001/*******************************************************************************
002 * Permission to use and distribute this software is granted under the BSD 2-Clause
003 * "Simplified" License (see http://opensource.org/licenses/BSD-2-Clause).
004 * Copyright (c) 2016-2023 Wilhelm Burger. All rights reserved.
005 * Visit https://imagingbook.com for additional details.
006 ******************************************************************************/
007package imagingbook.calibration.zhang.util;
008
009import ij.ImagePlus;
010import ij.gui.Plot;
011import ij.gui.PlotWindow;
012import ij.process.FloatProcessor;
013import imagingbook.calibration.zhang.Camera;
014import imagingbook.common.math.Matrix;
015import org.apache.commons.math3.linear.RealMatrix;
016
017import java.awt.Color;
018
019
020/**
021 * This class defines ImageJ-related utility methods used for camera calibration.
022 *
023 * @author WB
024 */
025public class IjUtil {
026
027        /**
028         * Creates a float-type image of the supplied matrix. Matrix rows/columns correspond to image rows/columns. Typical
029         * usage: makeImage("my matrix", M).show()
030         *
031         * @param title the image title
032         * @param M the matrix to be displayed
033         * @return the resulting image
034         */
035        public static ImagePlus makeImage(String title, RealMatrix M) {
036                float[][] fdata = Matrix.toFloat(M.transpose().getData());
037                ImagePlus im = new ImagePlus(title, new FloatProcessor(fdata));
038                im.show();
039                return im;
040        }
041        
042        // -------------------------------------------------------------------
043        
044        
045        public static PlotWindow plotLensDistortionFunction(Camera cam, double rmax) {
046                int n = 100;
047//              double[] K = cam.getK();
048                double[] xVals = new double[n];
049                double[] yVals = new double[n];
050                for (int i = 0; i < n; i++) {
051                        double r = i * (rmax / n);
052                        //double d = cam.distFun2(r * r); //getRadialDistortion(r, K);
053                        double d = cam.D(r);
054                        xVals[i] = r;
055                        yVals[i] = d;
056                }
057                Plot plot = new Plot("d(r)", "x", "f(x)");
058                plot.setColor(Color.BLUE);
059                plot.setLimits(0, 1.25, -0.1, 0.1);
060                plot.addPoints(xVals, yVals, Plot.LINE);
061                plot.draw();
062                return plot.show();
063        }
064        
065        // ---------------------------------------------------------------
066        
067//      public void showJacobian(MultivariateMatrixFunction jacobianFun, RealVector point) {
068//              double[][] J = jacobianFun.value(point.toArray());
069//              FloatProcessor fp = new FloatProcessor(J[0].length, J.length);
070//              for (int i = 0; i < J.length; i++) {
071//                      for (int j = 0; j < J[i].length; j++) {
072//                              fp.setf(j, i, (float) J[i][j]); 
073//                      }
074//              }
075//              (new ImagePlus("Jacobian", fp)).show();
076//      }
077}