example_iris_identify.cpp

Go to the documentation of this file.

 
#include "sfe_toolkit/sfe_core.h"
#include "sfe_toolkit/sfe_iris.h"
#include <regex>
#include <vector>
 
#include <iomanip>
#include <sstream>
 
#include "annotate.hpp"
#include "solvers.h"
#include "utils.hpp"
#include <unordered_map>
 
std::string image_match = "assets/iris/iris1_1.png";
std::string image_probe = "assets/iris/iris1_0.png";
std::string image_non_match = "assets/iris/iris0.png";
 
std::string solver_iris_detect = SOLVER_IRIS_DETECT;
std::string solver_iris_template = SOLVER_IRIS_TEMPLATE;
 
const auto SOLVER_PARAMETERS = std::vector<SFESolverParameter>{};
 
float identification_threshold = 0.6f;
 
SFEIrisTemplate extract_template(std::string image_path,
                                 SFESolver &detector_solver,
                                 SFESolver &template_solver) {
 
  SFEImage image{};
  DEFER(sfeImageFree(image));
  SFEError error{};
  { // STAGE 1: Load image
    // Load image data from file
    auto image_data = utils::readFile(image_path);
 
    // Decode image from data
    error = sfeImageDecode(image_data.data(), image_data.size(), &image);
    utils::checkError(error);
  }
 
  SFEDetection detected_iris = {};
  size_t iris_count = 1;
  { // STAGE 2: Detect iris in the image
 
    // Detect iris in the image
    error = sfeDetect(detector_solver, image, 0.1f, &detected_iris, &iris_count);
    utils::checkError(error);
 
    if (iris_count == 0 || detected_iris.confidence < 0.5) {
      throw std::runtime_error("No iris detected in the probe image.");
    }
 
    std::cout << "Found iris in the image. Extracting template..." << std::endl;
  }
 
  SFEIrisTemplate iris_template;
  { // STAGE 3 Extract probe iris template
    // Use template extraction solver to create new iris
    // template
    error = sfeIrisTemplateExtract(template_solver, image, &detected_iris,
                                   &iris_template);
    utils::checkError(error);
  }
  return iris_template;
}
 
void printHelp() {
  std::cout << "Help: Usage of the program." << std::endl;
  std::cout << "Options:" << std::endl;
  std::cout << "-h: Display help." << std::endl;
  std::cout << "-p: Probe image file." << std::endl;
  std::cout << "-g: Matching image file. Our \"iris database\"." << std::endl;
  std::cout << "-d: Path to detector solver." << std::endl;
  std::cout << "-e: Path to extraction solver." << std::endl;
  std::cout << "-t: Detection threshold. <0,1>" << std::endl;
  std::cout << "-i: Identification threshold. <0,1>" << std::endl;
}
 
int main(int argc, char *argv[]) {
 
  { // STAGE 0: Parse command line arguments
    // Map to store argument values
    std::unordered_map<std::string, std::string> args;
 
    // Parse command line arguments
    for (int i = 1; i < argc; ++i) {
      std::string arg = argv[i];
      if (arg[0] == '-') {
        if (arg == "-h") {
          printHelp();
          return 0;
        }
        // Check if there's a next argument and it isn't another option
        if (i + 1 < argc && argv[i + 1][0] != '-') {
          args[arg] = argv[++i];
        } else {
          std::cerr << "Option " << arg << " requires a value." << std::endl;
          return 1;
        }
      } else {
        std::cerr << "Unknown option: " << arg << std::endl;
        return 1;
      }
    }
 
    // Assign values to variables based on parsed arguments
    if (args.count("-p"))
      image_probe = args["-p"];
    if (args.count("-g"))
      image_match = args["-g"];
    if (args.count("-d"))
      solver_iris_detect = args["-d"];
    if (args.count("-e"))
      solver_iris_template = args["-e"];
    if (args.count("-i"))
      identification_threshold = std::stof(args["-i"]);
 
    utils::printFormatted("EXAMPLE PARAMETERS");
    // Print resource names
    std::cout << "Probe image: " << image_probe << std::endl;
    std::cout << "Matching image: " << image_match << std::endl;
 
    // Print solver names
    std::cout << "Iris detection solver: " << solver_iris_detect << std::endl;
    std::cout << "Iris template solver: " << solver_iris_template << std::endl;
 
    // Print other options
    std::cout << "Identification threshold: " << identification_threshold
              << std::endl;
  }
 
  utils::printToolkitInfo();
 
  SFEError error{};
 
  SFESolver detector_solver{};
  DEFER(sfeSolverFree(detector_solver));
  SFESolver landmarks_solver{};
  DEFER(sfeSolverFree(landmarks_solver));
  SFESolver template_solver{};
  DEFER(sfeSolverFree(template_solver));
  { // STAGE 1.1: loading solvers
    utils::printFormatted("LOADING SOLVERS");
    // Initialize Iris detection solver
    error = sfeSolverCreate(
        solver_iris_detect.c_str(), SOLVER_PARAMETERS.data(),
        SOLVER_PARAMETERS.size(), &detector_solver);
    utils::checkError(error);
 
    // Initialize template extraction solver
    error = sfeSolverCreate(
        solver_iris_template.c_str(), SOLVER_PARAMETERS.data(),
        SOLVER_PARAMETERS.size(), &template_solver);
    utils::checkError(error);
  }
 
  SFEIrisTemplate probe_iris_template;
  { // STAGE 1: Extract probe iris template
    utils::printFormatted("PROBE TEMPLATE EXTRACTION");
 
    probe_iris_template =
        extract_template(image_probe, detector_solver, template_solver);
  }
 
  { // STAGE 1.1: (optional) Print template attributes
    utils::printFormatted("PROBE TEMPLATE ATTRIBUTES");
 
    SFEIrisTemplateVersion version = {};
    error = sfeIrisTemplateVersion(&probe_iris_template, &version);
    utils::checkError(error);
 
    std::cout << "Version of the probe template: " << version.version_major
              << '.' << version.version_minor << std::endl;
 
    float quality = {};
    error = sfeIrisTemplateQuality(&probe_iris_template, &quality);
    utils::checkError(error);
 
    std::cout << "Quality of the probe template: " << quality << std::endl;
  }
 
  size_t gallery_size = 10;
  std::vector<SFEIrisTemplate> iris_template_gallery(gallery_size);
  { // STAGE 2: Create iris templates gallery
    utils::printFormatted("IRIS GALLERY EXTRACTION");
 
    auto matching_iris_template =
        extract_template(image_match, detector_solver, template_solver);
 
    auto non_matching_iris_template =
        extract_template(image_non_match, detector_solver, template_solver);
 
    // Fill the gallery with non-matching templates for demonstration purposes
    for (size_t i = 0; i < gallery_size; i++) {
      iris_template_gallery[i] = non_matching_iris_template;
    }
 
    // Add matching template to the gallery at index 5
    iris_template_gallery[5] = matching_iris_template;
  }
 
  size_t candidate_count = 1;
  std::vector<SFETemplateIdentificationCandidate> identification_results(
      candidate_count);
  { // STAGE 3: Identify the probe template
    utils::printFormatted("1:N IDENTIFICATION");
 
    error = sfeIrisTemplateIdentify(
        &probe_iris_template, iris_template_gallery.data(),
        iris_template_gallery.size(), identification_threshold,
        identification_results.data(), &candidate_count, 4);
    utils::checkError(error);
 
    if (candidate_count == 0) {
      std::cout << "No candidates found above identification threshold "
                << identification_threshold << std::endl;
      return 0;
    }
 
    std::cout << "Found " << identification_results.size()
              << " candidates above identification threshold "
              << identification_threshold << std::endl;
    for (auto &result : identification_results)
      std::cout << "Template index: #" << result.index
                << ", score: " << result.score << std::endl;
  }
 
  { // STAGE 4: 1:1 matching with top candidate to showcase 1:1 matching
    utils::printFormatted("1:1 MATCHING");
 
    if (identification_results.size() == 0) {
      std::cout << "No candidates found above identification threshold "
                << identification_threshold << std::endl;
      return 0;
    }
 
    // Since it's sorted vector by score, the best candidate is the first one
    auto match_iris_template =
        iris_template_gallery[identification_results[0].index];
 
    float match_confidence;
    error = sfeIrisTemplateMatch(&probe_iris_template, &match_iris_template,
                                 &match_confidence);
    utils::checkError(error);
 
    std::cout
        << "Matching score of probe template with the best template, score: "
        << match_confidence << std::endl;
  }
 
  utils::printFormatted("FINISHED");
}