Microsoft MVP성태의 닷넷 이야기
Graphics: 16. 3D 공간에서 두 점이 이루는 각도 구하기 [링크 복사], [링크+제목 복사]
조회: 3643
글쓴 사람
홈페이지
첨부 파일

3D 공간에서 두 점이 이루는 각도 구하기

예를 들어, 다음과 같이 3D 공간에서 면과 점이 있다고 가정해 보겠습니다.

angle_between_two_0.png
(오른손 좌표계 시스템의 x-y-z 축을 사용하고 있습니다.)

이때 점 C와 점 A가 이루는 각도는 어떻게 될까요? 간단하게 점과 점의 위치를 빼서 벡터를 만들고,

점 C와 카메라의 벡터 = 점 C - 점 A
                    = (-1.4, -2.9, -5)

그 벡터가 이루는 각을 3차원에 따라 각 축 별로 3개의 각도를 구할 수 있습니다.

  • x-축을 기준으로 y-z 평면의 각
  • y-축을 기준으로 x-z 평면의 각
  • z-축을 기준으로 x-y 평면의 각

역시 시각화를 해보는 게 좋겠지요? ^^ GeoGebra 앱에서 "Home" 모양의 아이콘을 누르면,

angle_between_two_1.png

차례대로 "y축을 기준으로 x-z 평면", "z축을 기준으로 x-y 평면", "x축을 기준으로 y-z 평면" 상에서의 물체를 볼 수 있는 기능이 제공됩니다.

자... 멋진 툴이 있으니, 그럼 우선 첫 번째 아이콘인 "y축을 기준으로 x-z 평면"을 다음과 같이 볼 수 있습니다. ^^

angle_between_two_2_y.png

보는 바와 같이 x-z 평면에서 벡터가 이루는 각을 쉽게 확인할 수 있습니다. 점 A에서 x-축으로 수선의 발을 그렸을 때,

angle_between_two_3.png

직각 삼각형이 되고 이때 우리가 아는 벡터 A의 x, z 성분의 값으로 인해 아크 탄젠트 함수를 사용하면 점 A에서 이루는 각도를 구할 수 있습니다.

각 A = arc_tangent(vx, vz);

C# 코드로 구해 보면 다음과 같습니다.

Vector3 v3 = new Vector3(-1.4f, -2.9f, -5f);

{
    double yAngle = Math.Atan2(v3.X, v3.Z);
    Console.WriteLine(RadianToDegree(yAngle));
}

실행해 보면, -164.357753796137도가 나오는데, 벡터가 음의 방향이어서 그런 것일 뿐 180도를 빼면 약 15.6도가 나옵니다. 동일한 각도를 Geogebra 기하를 이용해서도 구해 볼까요? ^^

2차원 좌표계에서 x, z의 성분으로 점을 표시한 후 그와 연결한 벡터를 그리고 그 사잇각을 다음과 같이 구할 수 있습니다.

angle_between_two_4.png

잘 맞아 들어가는군요. ^^




마찬가지 방식으로 "z축을 기준으로 x-y 평면"과,

angle_between_two_2_z.png

"x축을 기준으로 y-z 평면" 상의 모습을 확인하며,

angle_between_two_2_x.png

각도를 구해 보면 이렇습니다.

Vector3 v3 = new Vector3(-1.4f, -2.9f, -5f);
Console.WriteLine("xAngle");
{
    double xAngle = Math.Atan2(v3.Y, v3.Z);
    Console.WriteLine(RadianToDegree(xAngle));
}

Console.WriteLine("yAngle");
{
    double yAngle = Math.Atan2(v3.X, v3.Z);
    Console.WriteLine(RadianToDegree(yAngle));
}

Console.WriteLine("zAngle");
{
    double zAngle = Math.Atan2(v3.Y, v3.X);
    Console.WriteLine(RadianToDegree(zAngle));
}

/*
출력 결과

xAngle
-149.886266031276
yAngle
-164.357753796137
zAngle
-115.769326504636
*/

(첨부 파일은 이 글에서 예제로 사용한 geogebra 파일입니다.)




참고로, 2개의 벡터가 이루는 각은 아크코사인을 이용해 구할 수 있습니다.

vector v1 = ...;
vector v2 = ...;

vector nv1 = normalize(v1);
vector nv2 = normalize(v2);

float angle = acos(mul(nv1, nv2));

일반적으로 acos의 정의역이 -1 ~ 1이기 때문에 정규화된 벡터의 내적 연산으로 NaN 값을 반환하는 경우는 없습니다. 문제는 0 벡터인 경우 normalize 하는 과정에서 NaN이 나온다는 것입니다.

vector v = vector(0, 0, 0, 0);
vector normal = normalize(v); // normal == NaN

왜냐하면 정규화라는 것이 벡터 크기로 나누는 것인데,

${
\hat v = \left ( \frac{ v_x }{ \Vert v \Vert }, \frac{ v_y }{ \Vert v \Vert }, \frac{ v_z }{ \Vert v \Vert } \right )
}$

분모가 0이 되므로 normalize의 결과가 NaN이 나옵니다. 이 때문에 mul 연산도 NaN이 되고, acos까지 NaN으로 끝납니다.




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]

[연관 글]





[최초 등록일: ]
[최종 수정일: 8/1/2018 ]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer@outlook.com

비밀번호

댓글 쓴 사람
 




... 16  17  18  19  20  21  22  23  24  25  26  27  28  29  [30]  ...
NoWriterDateCnt.TitleFile(s)
11397정성태12/12/20173959.NET Framework: 703. 양자 컴퓨팅을 위한 마이크로소프트의 Q# 언어
11396정성태12/8/201712197개발 환경 구성: 343. Visual Studio - 리눅스 용 프로젝트의 인텔리센스를 위한 헤더 파일 처리 방법 [3]
11395정성태12/8/20172898오류 유형: 441. 이벤트 로그 - Time Provider NtpClient: No valid response has been received from domain controller
11394정성태12/8/20172814개발 환경 구성: 342. 비주얼 스튜디오에서 실행하던 ASP.NET Core (.NET Framework) 응용 프로그램을 명령행에서 실행하는 방법
11393정성태12/7/20174315Windows: 145. 윈도우 10 빌드 17046부터 WSL에서 백그라운드 작업 지원 [3]
11392정성태12/7/20172216개발 환경 구성: 341. openSUSE에 닷넷 코어 설치
11391정성태12/7/20172745개발 환경 구성: 340. WSL을 이용해 윈도우 PC 1대에서 openSUSE 응용 프로그램을 Visual Studio로 개발하는 방법 [1]
11390정성태12/7/20177781개발 환경 구성: 339. WSL을 이용해 윈도우 PC 1대에서 Linux 응용 프로그램을 Visual Studio로 개발하는 방법 [3]
11389정성태12/7/20172634오류 유형: 440. .NET Core 오류 - 0x80131620 Unable to load DLL 'libuv'
11388정성태10/4/20184040개발 환경 구성: 338. WSL 또는 Ubuntu에 닷넷 코어 설치 [2]
11387정성태12/6/20173953오류 유형: 439. 이벤트 로그 - Data Sharing Service 서비스의 %%3239247874 오류 메시지
11386정성태12/5/20172661오류 유형: 438. Hyper-V - '...' failed to add device 'Virtual CD/DVD Disk'
11385정성태12/14/20176790VC++: 121. DXGI를 이용한 윈도우 화면 캡처 소스 코드(Visual C++) [3]파일 다운로드1
11384정성태12/5/20173338오류 유형: 437. Visual C++ - Cannot open include file: 'SDKDDKVer.h'
11383정성태9/18/20193823디버깅 기술: 110. 비동기 코드 실행 중 예외로 인한 ASP.NET 프로세스 비정상 종료 현상 [1]
11382정성태12/4/20174553오류 유형: 436. System.Data.SqlClient.SqlException (0x80131904): Connection Timeout Expired 예외 발생 시 "[Pre-Login] initialization=48; handshake=1944;" 값의 의미
11381정성태11/30/20172804.NET Framework: 702. 한글이 포함된 바이트 배열을 나눈 경우 한글이 깨지지 않도록 다시 조합하는 방법(두 번째 이야기)파일 다운로드1
11380정성태11/30/20173353디버깅 기술: 109. windbg - (x64에서의 인자 값 추적을 이용한) Thread.Abort 시 대상이 되는 스레드를 식별하는 방법
11379정성태11/30/20172657오류 유형: 435. System.Web.HttpException - Session state has created a session id, but cannot save it because the response was already flushed by the application.
11378정성태11/29/20174086.NET Framework: 701. 한글이 포함된 바이트 배열을 나눈 경우 한글이 깨지지 않도록 다시 조합하는 방법 [1]파일 다운로드1
11377정성태11/29/20173705.NET Framework: 700. CommonOpenFileDialog 사용 시 사용자가 선택한 파일 목록을 구하는 방법 [3]파일 다운로드1
11376정성태11/28/20173936VS.NET IDE: 123. Visual Studio 편집기의 \r\n (crlf) 개행을 \n으로 폴더 단위로 설정하는 방법
11375정성태11/28/20172665오류 유형: 434. Visual Studio로 ASP.NET 디버깅 중 System.Web.HttpException - Could not load type 오류
11374정성태8/9/20184492사물인터넷: 14. 라즈베리 파이 - (윈도우의 NT 서비스처럼) 부팅 시 시작하는 프로그램 설정
11373정성태7/16/20183334오류 유형: 433. Raspberry Pi/Windows 다중 플랫폼 지원 컴파일 관련 오류 기록
11372정성태10/2/20185880사물인터넷: 13. 윈도우즈 사용자를 위한 라즈베리 파이 제로 W 모델을 설정하는 방법 [3]
... 16  17  18  19  20  21  22  23  24  25  26  27  28  29  [30]  ...