Ejemplo n.º 1
0
 /**
  * The results of the {@link SR_EELS_CharacterisationPlugin} plugin are parsed.
  *
  * <p>This function extracts the values that describe the width of a spectrum depending on its
  * position on the camera.
  *
  * @return a polynomial that fits the given data points
  */
 private SR_EELS_Polynomial_2D getFunctionWidth() {
   final double[][] vals = new double[importer.size()][3];
   int i = 0;
   for (final float[] point : importer) {
     // The width of the spectrum -> z value
     vals[i][0] = point[8];
     // Coordinate on the energy dispersive axis -> x value
     // The same value for top, centre and bottom.
     vals[i][1] = point[0] - CameraSetup.getFullWidth() / 2;
     // coordinate on the lateral axis -> y value
     // The indices for centre, top and bottom are 2, 4 and 6.
     vals[i][2] = point[2] - CameraSetup.getFullHeight() / 2;
     i++;
   }
   /*
    * Define the orders of the 2D polynomial.
    */
   final int m = 2;
   final int n = 2;
   final SR_EELS_Polynomial_2D func = new SR_EELS_Polynomial_2D(m, n);
   final double[] b_fit = new double[(m + 1) * (n + 1)];
   Arrays.fill(b_fit, 1.);
   final LMA lma = new LMA(func, b_fit, vals);
   lma.fit();
   /*
    * TODO: Output information about the fit using IJ.log
    */
   return new SR_EELS_Polynomial_2D(m, n, b_fit);
 }
Ejemplo n.º 2
0
 /**
  * The results of the {@link SR_EELS_CharacterisationPlugin} plugin are parsed.
  *
  * <p>This function extracts the values that describe the pathway of the borders of a spectrum.
  *
  * @return a polynomial that fits the given data points
  */
 private SR_EELS_Polynomial_2D getFunctionBorders() {
   final double[][] vals = new double[3 * importer.size()][3];
   int i = 0;
   for (final float[] point : importer) {
     for (int j = 0; j < 3; j++) {
       // y coordinate of the fit function at the centre of the
       // image/camera ->
       // z value
       vals[i + j][0] = point[2 + 2 * j] - CameraSetup.getFullWidth() / 2;
       // Coordinate on the energy dispersive axis -> x value
       // The same value for top, centre and bottom.
       vals[i + j][1] = point[0] - CameraSetup.getFullHeight() / 2;
       // coordinate on the lateral axis -> y value
       // The indices for centre, top and bottom are 2, 4 and 6.
       vals[i + j][2] =
           importer.getYInterceptPoint(i / 3)[11 + j] - CameraSetup.getFullHeight() / 2;
     }
     i += 3;
   }
   /*
    * Define the orders of the 2D polynomial.
    */
   final int m = 3;
   final int n = 2;
   final SR_EELS_Polynomial_2D func = new SR_EELS_Polynomial_2D(m, n);
   final double[] a_fit = new double[(m + 1) * (n + 1)];
   Arrays.fill(a_fit, 1.);
   final LMA lma = new LMA(func, a_fit, vals);
   lma.fit();
   /*
    * TODO: Output information about the fit using IJ.log
    */
   return new SR_EELS_Polynomial_2D(m, n, a_fit);
 }