등고선 그리기

이번 작업을 하면서 등고선 그리는 일이 필요하게 되었습니다.
며칠동안 인터넷을 뒤져 가면서 등고선 그리는 루틴( 소스이면 좋겠다고 생각 함)을 찾아봤지만 
이거다 싶은 것을 발견하지 못했습니다.
그나마 겨우 찾아낸것이 Kriging 알고리즘으로 등고선을 그린다는 것이었습니다.
Kriging 알고리즘은 그리 복잡하지 않은 몇개의 식으로 표현되어 있었는데 이것을 어떻게 프로그래밍해야 할지 알수가 없었습니다.
다시 인터넷을 뒤지던 중 나의 상황에 맞는 루틴을 코드구루에서 찾을 수 있었습니다.
이것은 아주 단순했습니다.
저의 경우는 1mm 간격으로 놓여 있는 그리드의 측정점(x , y , z)들을 미리 알고 있는 상태였습니다.
아래 이미지는 측정 데이터들의 일부분을 캡쳐한 것입니다.

TOP에서 봤을때

ISO로 봤을때

이 측정점들로 삼각망들을 형성할 수 있습니다.
min z ~ max z의 간격을 7등분하였고 각 등분된 간격에 무지개 색상을 각각 할당하였습니다.
등분된 간격에 들어가는 삼각망을 할당된 색상으로 칠하기만 하면 되는 것입니다.
좀더 자세히 생각해 보면, 최초의 삼각망은 LEVEL( min z ~ max z를 등분하는 값)에 의해 분할이 일어나게 됩니다.
이렇게 분할이 일어날때 삼각망은 2개의 다각형으로 분할되는데
이 다각형을 다시 삼각형으로 나누어줍니다.
$N$개의 정점을 가지는 다각형은 $N-2$개의 삼각형으로 나눌 수 있습니다.
$N-1,i,i+1$의 3개의 정점을 가지고 삼각형을 생성한다고 생각하면 됩니다.
$(i = 0,...,N-3)$

이러한 작업은 모든 삼각망에 대해 모든 LEVEL에 걸쳐 수행하게 됩니다.
결과로 생기는 모든 삼각망들은 모든 LEVEL에 의해 나누어지게 됩니다.
삼각망들을 2D로 내려 할당된 색상으로 칠해주면 됩니다.

여담으로 소스포지에 보면 FMesh라는 프로그램이 있는데
이 프로그램을 분석해보면 3D 좌표를 가지고 먼저 삼각망을 구성하고 나서
앞서 이야기한 내용대로 삼각망을 분할하여 CONTOUR 곡면을 생성하고 있었습니다.

저의 경우는 이미 삼각망이 구성되었는 상황이니 쉽게 작업을 진행할 수 있었습니다.
아래는 최종 출력 파일입니다.