File: combinatorialthresholder.cpp

package info (click to toggle)
aoflagger 3.1.0-2
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 4,868 kB
  • sloc: cpp: 52,164; python: 152; sh: 60; makefile: 17
file content (53 lines) | stat: -rw-r--r-- 1,854 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
#include "../../structures/image2d.h"

#include "combinatorialthresholder.h"
#include "thresholdtools.h"

void CombinatorialThresholder::HorizontalVarThreshold(const Image2D* input,
                                                      Mask2D* mask,
                                                      size_t length,
                                                      num_t threshold) {
  size_t width = input->Width() - length + 1;
  for (size_t y = 0; y < input->Height(); ++y) {
    for (size_t x = 0; x < width; ++x) {
      bool flag = true;
      for (size_t i = 0; i < length; ++i) {
        if (input->Value(x + i, y) < threshold &&
            input->Value(x + i, y) > -threshold) {
          flag = false;
          break;
        }
      }
      if (flag) {
        for (size_t i = 0; i < length; ++i) mask->SetValue(x + i, y, true);
      }
    }
  }
}

void CombinatorialThresholder::VerticalVarThreshold(const Image2D* input,
                                                    Mask2D* mask, size_t length,
                                                    num_t threshold) {
  size_t height = input->Height() - length + 1;
  for (size_t y = 0; y < height; ++y) {
    for (size_t x = 0; x < input->Width(); ++x) {
      bool flag = true;
      for (size_t i = 0; i < length; ++i) {
        if (input->Value(x, y + i) <= threshold &&
            input->Value(x, y + i) >= -threshold) {
          flag = false;
          break;
        }
      }
      if (flag) {
        for (size_t i = 0; i < length; ++i) mask->SetValue(x, y + i, true);
      }
    }
  }
}

void CombinatorialThresholder::VarThreshold(const Image2D* input, Mask2D* mask,
                                            size_t length, num_t threshold) {
  HorizontalVarThreshold(input, mask, length, threshold);
  VerticalVarThreshold(input, mask, length, threshold);
}