Remove LaueOps::generateMisorientationColors due to licensing.

Original author LPGL the original code which was copied into this repository
by accident.

Signed-off-by: Michael Jackson <mike.jackson@bluequartz.net>

Author: imikejackson

Source/EbsdLib/LaueOps/CubicLowOps.cpp

@@ -1135,128 +1135,6 @@ EbsdLib::UInt8ArrayType::Pointer CubicLowOps::generateIPFTriangleLegend(int imag
   return image;
 }
 
-// -----------------------------------------------------------------------------
-//
-// -----------------------------------------------------------------------------
-EbsdLib::Rgb CubicLowOps::generateMisorientationColor(const QuatD& q, const QuatD& refFrame) const
-{
-  throw std::out_of_range("generateMisorientationColor::generateMisorientationColor NOT Implemented");
-
-  // double n1, n2, n3, w;
-  double x, x1, x2, x3, x4;
-  double y, y1, y2, y3, y4;
-  double z, z1, z2, z3, z4;
-  double k, h, s, v, c, r, g, b;
-
-  QuatD q1 = q;
-  QuatD q2 = refFrame;
-
-  OrientationD axisAngle = calculateMisorientation(q1, q2);
-
-  // eq c7.1
-  k = tan(axisAngle[3] / 2.0);
-  x = axisAngle[0];
-  y = axisAngle[1];
-  z = axisAngle[2];
-  OrientationType rod(x, y, z, k);
-  rod = getMDFFZRod(rod);
-  x = rod[0];
-  y = rod[1];
-  z = rod[2];
-  k = rod[3];
-
-  // eq c7.2
-  k = atan2(y, x);
-  if(k < M_PI_4)
-  {
-    x1 = z;
-    y1 = x;
-    z1 = y;
-  }
-  else
-  {
-    x1 = y;
-    y1 = z;
-    z1 = x;
-  }
-
-  // eq c7.3
-  // 3 rotation matricies (in paper) can be multiplied into one (here) for simplicity / speed
-  // g1*g2*g3 = {{sqrt(2/3), 0, 1/sqrt(3)},{-1/sqrt(6), 1/sqrt(2), 1/sqrt(3)},{-1/sqrt(6), 1/sqrt(2), 1/sqrt(3)}}
-  x2 = x1 * sqrt(2.0f / 3.0) - (y1 + z1) / sqrt(6.0);
-  y2 = (y1 - z1) / sqrt(2.0);
-  z2 = (x1 + y1 + z1) / sqrt(3.0);
-
-  // eq c7.4
-  k = (sqrt(3.0) * y2 + x2) / (2.0f * pow(x2 * x2 + y2 * y2, 1.5f));
-  x3 = x2 * (x2 + sqrt(3.0) * y2) * (x2 - sqrt(3.0) * y2) * k;
-  y3 = y2 * (y2 + sqrt(3.0) * x2) * (sqrt(3.0) * x2 - y2) * k;
-  z3 = z2 * sqrt(3.0);
-
-  // eq c7.5 these hsv are from 0 to 1 in cartesian coordinates
-  x4 = -x3;
-  y4 = -y3;
-  z4 = z3;
-
-  // convert to traditional hsv (0-1)
-  h = fmod(atan2(y4, x4) + 2.0f * M_PI, 2.0f * M_PI) / (2.0f * M_PI);
-  s = sqrt(x4 * x4 + y4 * y4);
-  v = z4;
-  if(v > 0)
-  {
-    s = s / v;
-  }
-
-  // hsv to rgb (from wikipedia hsv/hsl page)
-  c = v * s;
-  k = c * (1 - fabs(fmod(h * 6, 2) - 1)); // x in wiki article
-  h = h * 6;
-  r = 0;
-  g = 0;
-  b = 0;
-
-  if(h >= 0)
-  {
-    if(h < 1)
-    {
-      r = c;
-      g = k;
-    }
-    else if(h < 2)
-    {
-      r = k;
-      g = c;
-    }
-    else if(h < 3)
-    {
-      g = c;
-      b = k;
-    }
-    else if(h < 4)
-    {
-      g = k;
-      b = c;
-    }
-    else if(h < 5)
-    {
-      r = k;
-      b = c;
-    }
-    else if(h < 6)
-    {
-      r = c;
-      b = k;
-    }
-  }
-
-  // adjust lumosity and invert
-  r = (r + (v - c));
-  g = (g + (v - c));
-  b = (b + (v - c));
-
-  return EbsdLib::RgbColor::dRgb(static_cast<int32_t>(r * 255), static_cast<int32_t>(g * 255), static_cast<int32_t>(b * 255), 255);
-}
-
 // -----------------------------------------------------------------------------
 CubicLowOps::Pointer CubicLowOps::NullPointer()
 {

Source/EbsdLib/LaueOps/CubicLowOps.h

@@ -205,14 +205,6 @@ class EbsdLib_EXPORT CubicLowOps : public LaueOps
    */
   EbsdLib::Rgb generateRodriguesColor(double r1, double r2, double r3) const override;
 
-  /**
-   * @brief generateMisorientationColor Generates a color based on the method developed by C. Schuh and S. Patala.
-   * @param q Quaternion representing the direction
-   * @param refDir The sample reference direction
-   * @return Returns the ARGB Quadruplet EbsdLib::Rgb
-   */
-  EbsdLib::Rgb generateMisorientationColor(const QuatD& q, const QuatD& refFrame) const override;
-
   /**
    * @brief generatePoleFigure This method will generate a number of pole figures for this crystal symmetry and the Euler
    * angles that are passed in.