http://www.labbookpages.co.uk/software/imgProc/otsuThreshold.html

TM 좌표를 WGS84로 변환하는 엑셀 수식입니다.

국토지리정보원에서 제공하는 표준 정보를 이용해서 만든 자료입니다.

링크:  https://autoc.tistory.com/315

https://modoocode.com/224

https://doitnow-man.tistory.com/215

https://neokido.tistory.com/entry/TRACE-C-Debug%EB%A5%BC-%EC%9C%84%ED%95%9C-TRACE-%EC%B0%8D%EA%B8%B0

https://seujeum.tistory.com/13

std::vector의 push_back() 함수를 사용하면 충돌현상이 발생함

1. 따라서 push_back()이 발생할 때 아래처럼 처리하면 push_back()의 효과를 낼 수 있음

output.resize(input.rows);
int k = 0;

#pragma omp parallel for shared(k, input)
for( int i = 0; i < input.rows; i++ )
{
    if(IsGoodMatch(input[I])
    { 
        Newvalues newValues;
        ...
        // ! prevent other threads to modify k !
        output[k] = newValues;
        k++;
        // ! allow other threads to modify k again !
    }
} 

output.resize(k);

2. 또는 멀티 쓰레드 처리를 위해 개발된 concurrency vector의 push_back()을 사용하면 됨

#include <concurrent_vector.h>
Concurrency::concurrent_vector<int> // in c++11.

https://stackoverflow.com/questions/19268217/stdvector-push-back-fails-when-used-in-a-parallel-for-loop

void func1() 
{
    #define update_leftmost(i) \
    { \
        const vec2_t &p  = polygon[i];                                   \
        const vec2_t &lm = polygon[iLeftMost];                           \
        if (p.x < lm.x || (p.x == lm.x && p.y < lm.y)) iLeftMost = i;    \
    }
    
    #define init_vertex(i, iPrev, iNext) \
        vertices[i].index = i;                                           \
        vertices[i].prev = &vertices[iPrev];                             \
        vertices[i].next = &vertices[iNext];                             \
        vertices[i].prev->index = iPrev;                                 \
        vertices[i].next->index = iNext;                                 \
        vertices[i].winding_value = winding_value(polygon, vertices[i]); \
        vertices[i].reflex = false;

    // initialize vertices

    init_vertex(0    , N - 1, 1);
    init_vertex(N - 1, N - 2, 0);

    update_leftmost(N - 1);

    for (int i = 1; i < N - 1; i++)
    {
        init_vertex(i, i - 1, i + 1);
        update_leftmost(i);
    }
}

Distance from a point to a line/segment

// 일반 함수
/* 점과 직선 사이의 거리 계산 */
static std::array<int, 3> cross(const std::array<int, 3> &a,
	const std::array<int, 3> &b)
{
	std::array<int, 3> result;
	result[0] = a[1] * b[2] - a[2] * b[1];
	result[1] = a[2] * b[0] - a[0] * b[2];
	result[2] = a[0] * b[1] - a[1] * b[0];
	return result;
}

static double point_to_line_distance(const cv::Point &p, const cv::Vec4i &line)
{
	std::array<int, 3> pa{ { line[0], line[1], 1 } };
	std::array<int, 3> pb{ { line[2], line[3], 1 } };
	std::array<int, 3> l = cross(pa, pb);
	return std::abs((p.x * l[0] + p.y * l[1] + l[2])) * 1.0 /
		std::sqrt(double(l[0] * l[0] + l[1] * l[1]));
}

// 또는

// 템플릿 함수화
template<typename point3_t>
point3_t cross(const point3_t vec1_pos, const point3_t vec2_pos)
{
	return point3_t(
		vec1_pos.y * vec2_pos.z - vec1_pos.z * vec2_pos.y,
		vec1_pos.z * vec2_pos.x - vec1_pos.x * vec2_pos.z,
		vec1_pos.x * vec2_pos.y - vec1_pos.y * vec2_pos.x);
}
    
/* 점과 직선 사이의 거리 계산 */
template <typename point2_t, typename vec4_t>
double point_to_line_distance(const point2_t &p, const vec4_t &line)
{
	struct point3_d { double x, y, z; point3_d() {}; point3_d(double a, double b, double c) { x = a; y = b; z = c; }; };
	point3_d line_s(line[0], line[1], 1);
	point3_d line_e(line[2], line[3], 1);
	point3_d line_parameters = cross(line_s, line_e);	// line_parameters[3] mean {a, b, c} of line equation(ax + by + c = 0).
	return std::abs((line_parameters.x * p.x + line_parameters.y * p.y + line_parameters.z)) * 1.0 / std::sqrt(double(line_parameters.x * line_parameters.x + line_parameters.y * line_parameters.y));
}

https://stackoverflow.com/questions/12132352/distance-from-a-point-to-a-line-segment

C++ 07.13 - 동적으로 배열 할당하기 (Dynamically allocating arrays)
https://boycoding.tistory.com/205

헤더파일에 작성한 문법에 이상이 없는 함수 코드가 빌드되지 않고 error LNK2005 가 발생할 때는,

static 이나 inline 함수로 선언한다.

Some possible solutions are:

1. Declare the functions as static, which implies internal linkage
2. Inline the functions
3. Move the implementation out of the header and into a c/c++ file

Read more here: What is external linkage and internal linkage?

https://stackoverflow.com/questions/6964819/function-already-defined-error-in-c

3차원 영상(공간?) 도메인에서 3차원 배열의 인덱스를 1차원으로 변환하는 식은 아래와 같음

public int to1D( int x, int y, int z ) {
    return (z * xMax * yMax) + (y * xMax) + x;
}

public int[] to3D( int idx ) {
    final int z = idx / (xMax * yMax);
    idx -= (z * xMax * yMax);
    final int y = idx / xMax;
    final int x = idx % xMax;
    return new int[]{ x, y, z };
}

https://stackoverflow.com/questions/7367770/how-to-flatten-or-index-3d-array-in-1d-array