Spectral/ColorConvert

// ---------------------------------------------------------------------------
// ---------------------------------------------------------------------------

#include <iostream>

#include <cvb/spectral/spectral.hpp>
#include <cvb/spectral/cube.hpp>
#include <cvb/spectral/arithmetic.hpp>
#include <cvb/spectral/interpolator.hpp>
#include <cvb/spectral/convert_color_space.hpp>

int main()
{
  try
  {
    // Load ENVI-Cubes
    Cvb::String objectHeader(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/leaf.hdr")));
    Cvb::String objectBinary(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/leaf.bin")));
    Cvb::String whiteRefHeader(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/WRef.hdr")));
    Cvb::String whiteRefBinary(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/WRef.bin")));
    Cvb::String blackRefHeader(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/BRef.hdr")));
    Cvb::String blackRefBinary(Cvb::ExpandPath(Cvb::InstallPath() + CVB_LIT("tutorial/Spectral/Images/Leaf/BRef.bin")));

    std::cout << "Load cube from: \n\t" <<
      objectHeader << "\n\t" << objectBinary << std::endl;
    std::cout << "White Reference: \n\t" <<
      whiteRefHeader << "\n\t" << whiteRefBinary << std::endl;
    std::cout << "Black Reference: \n\t" <<
      blackRefHeader << "\n\t" << blackRefBinary << std::endl;

    auto whiteRef = Cvb::Cube::Load(whiteRefHeader, whiteRefBinary);
    auto blackRef = Cvb::Cube::Load(blackRefHeader, blackRefBinary);
    auto imgCube = Cvb::Cube::Load(objectHeader, objectBinary);

    // Normalize
    auto normalizedCube = Cvb::Spectral::Normalization(*imgCube, *whiteRef, *blackRef, Cvb::Spectral::NormalizationMethod::AverageReferences1);
    std::cout << "Normalization done!" << std::endl;

    // Convert to Lab and RGB
    auto interpolator = Cvb::Interpolator::Create(*normalizedCube,
      Cvb::Spectral::StdObserver::CIE2006_2deg,
      Cvb::Spectral::StdIlluminant::D65,
      Cvb::Spectral::InterpolationMethod::Linear);
    auto labImg = Cvb::Spectral::CubeToLab(*normalizedCube, *interpolator);
    std::cout << "Lab calculation done!" << std::endl;
    auto rgbImg = Cvb::Spectral::LabToRGB8(*labImg, *interpolator);
    std::cout << "sRGB calculation done!" << std::endl;

    Cvb::String labFilename(Cvb::ExpandPath(CVB_LIT("%TMP%LabImage.tiff")));
    Cvb::String rgbFilename(Cvb::ExpandPath(CVB_LIT("%TMP%sRGBImage.tiff")));
    labImg->Save(labFilename);
    rgbImg->Save(rgbFilename);
    std::cout << "Saved Lab image to:\n\t" <<
      labFilename << std::endl;
    std::cout << "Saved sRGB image to:\n\t" <<
      rgbFilename << std::endl;
  }
  catch (const std::exception& error)
  {
    std::cout << error.what() << std::endl;
  }
}