void ThreePhaseDecoder::makePhase() {
int n = width * height;
float i1, i2, i3;
for (int i = 0; i < n; i++) {
i1 = (float) graySequence[0][i];
i2 = (float) graySequence[1][i];
i3 = (float) graySequence[2][i];
#ifdef USE_GAMMA
i1 = Gamma(i1 / 255.0, gamma) * 255.0;
i2 = Gamma(i2 / 255.0, gamma) * 255.0;
i3 = Gamma(i3 / 255.0, gamma) * 255.0;
#endif
range[i] = findRange(i1, i2, i3);
mask[i] = range[i] <= rangeThreshold;
ready[i] = !mask[i];
reflectivity[i] = (byte) ((i1 + i2 + i3) / 3);
if(ready[i])
phase[i] = atan2f(sqrtf(3) * (i1 - i3), 2.f * i2 - i1 - i3) / (float) TWO_PI;
}
#ifdef LINEARIZE_PHASE
if(linearize) {
buildLut();
applyLut();
}
#endif
}
/**
* Matlab-os bwmorph emulálás bridge kapcsolóval
* A bináris képen a nem közvetlen, 1 távolságra lévő bitek összekapcsolása
*
*/
cv::Mat bwMorphBridge(const cv::Mat *image)
{
/* Octave-ból LookUpTable */
cv::Mat lookUpTable =
(cv::Mat_<float>(1, 1024) << 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1);
cv::Mat res;
cv::threshold(*image, res, 0, 1, cv::THRESH_BINARY);
int lastPC = std::numeric_limits<int>::max();
for (int pointCount = cv::countNonZero(res); pointCount < lastPC; pointCount = cv::countNonZero(res))
{
lastPC = pointCount;
res = applyLut(image, lookUpTable);
}
res *= 255;
return res;
}
/*
* 2 LookUpTable alkalmazása a képre (Matlab-os átírás miatt)
*/
cv::Mat applyLut(const cv::Mat *image, const cv::Mat lutOne, const cv::Mat lutTwo)
{
cv::Mat dst = applyLut(image, lutOne);
dst = applyLut(image, lutTwo);
return dst;
}
