최근 MEMS (Micro Electro Mechanical System) 기술이 발전함에 따라 크기가 작고 저가의 IMMU (Inertial and Magnetic Measurement Unit)가 차량, 무인 잠수정, 수중 로봇, 스마트폰 등 3차원 자세 정보를 필요로 하는 다양한 분야에서 사용되고 있다.^1-5 여기서, IMMU는 3축 가속도계(Accelerometer), 3축 자이로 센서(Gyroscope)와 3축 마그네토미터(Magnetometer)로 구성된 센서모듈로서, 물체의 3차원 자세추정 정확성을 높이기 위해서는 센서 사용 전에 보정(Calibration) 과정이 필요하다. 그 중에서도 방위각(Azimuth) 정보를 제공하는 마그네토미터는 다른 두 센서에 비해 주변환경의 영향을 크게 받기 때문에 반드시 보정 과정을 거쳐야 한다.^6-8

3축 마그네토미터는 지구자기장의 크기와 방향을 측정하여 방위각 정보를 얻기 위해 사용된다. 하지만 센서 제조 과정에서 생기는 오차와 주변환경에 의해 발생하는 오차로 인하여 신호 왜곡이 발생할 수 있다. 여기서, 센서 자체가 가지는 오차는 민감도(Scale Factor), 비직교성(Non-Orthogonality)과 오프셋(Offset)이 있고, 주변환경에 의해 발생하는 오차는 강성 왜곡(Hard-Iron Distortion)과 약성 왜곡(Soft-Iron Distortion)이 있다. 이러한 오차들로 인해 보정되지 않은 마그네토미터를 이용하여 방위각 추정시 정확성이 크게 저하될 수 있다.⁹

마그네토미터를 통한 방위각 추정 정확성을 높이기 위하여 다양한 보정 기법들이 연구되었다. 가장 기본적인 보정 기법으로는 스윙 방식(Swing Method)이 있다. 스윙 방식은 센서의 각 축을 중심으로 회전시켜 얻은 데이터를 통해 오프셋과 강성 왜곡에 해당하는 오차를 보정하는 기법으로 민감도, 비직교성과 약성 왜곡으로 인한 오차는 보정할 수 없다. 또한, 센서를 회전시키는 과정에서 회전축과 센서축을 완벽하게 일치시켜야 한다는 단점이 있다.¹⁰

스윙 방식과 달리 (1) 오프셋과 강성 왜곡 외에도 다양한 오차 보정이 가능하며, (2) 일정회전축에 의한 데이터가 아닌 임의 회전에 의한 데이터가 사용가능한 타원체 피팅(Ellipsoid-Fitting) 방식이 연구되어 왔다.^11-17 타원체 피팅 방식의 기본 원리는 다음과 같다. 이론적으로 3축 마그네토미터를 임의로 회전시키며 일정한 지구자기장을 측정하여 공간상에 출력하면 구 형태를 형성해야 한다. 하지만 마그네토미터 오차 요인들로 인하여 편향된 타원체 형태를 형성한다. 이에 따라 측정 데이터를 타원체 방정식으로 표현할 수 있으며, 이때 타원체 방정식의 계수는 마그네토미터 오차 파라미터(Parameter)들로 이루어진 함수이다. 즉, 3차원상에서 회전하여 계측된 3축 마그네토미터 측정값을 이용하여 타원체 방정식 계수 및 오차 파라미터를 결정할 수 있다.

최근 들어 다음과 같은 타원체 피팅 방식에 기반한 다양한 방법들이 발표되었다: 타원체 방정식의 대각화 분해(Eigen-Decomposition)를 통한 선형 최소자승법(Least Square)을 사용하는 방법(이하 Method A),¹¹ 비선형(Nonlinear) 최소자승법을 적용한 방법(이하 Method B),¹² 타원체 방정식의 계수를 구하는 단계와 계수로부터 파라미터를 구하는 단계를 명시적으로 구분하여 적용한 2단계 알고리즘(이하 Method C),¹³ 그리고 Method B와 보정 모델과 방정식구성에 차이를 보이는 또다른 비선형 최소자승법 적용방법 (이하 Method D).¹⁴ 이외에도 최대가능도법(Maximum Likelihood Estimation),¹⁵ 격차 진화 알고리즘(Differential Evolution Algorithm)¹⁶ 그리고 확장 칼만 필터(Extended Kalman Filter)¹⁷ 등을 사용한 기법들이 연구되었다.

위에 기술한 바와 같이 다양한 타원체 피팅 방식 보정기법이 발표되었으나, 현재까지 이들 기법에 대한 성능비교와 장단점 분석은 이루어 지지 않았다. 본 논문에서는 타원체 피팅을 기반으로 하는 보정 기법 Method A-D의 방위각 추정성능을 비교하고자 한다. 이때, Method A와 B에서는 마그네토미터 오차들이 모두 고려된 반면, Method C와 D에서는 약성 왜곡이 고려되지 않았다. 따라서 Method C와 D에서 약성 왜곡 미고려에 따른 영향을 제거하기 위하여 Method C와 D에 Method B의 오차모델을 적용한 Method C′과 D′을 추가로 비교한다. 본 논문은 이러한 비교를 통해 각 방법의 장단점과 특징을 분석하고, 보정성능의 우열을 가늠해 보고자 한다.

Table 1은 각 마그네토미터 보정 기법의 특성을 나타낸 표이다. Method C와 D에서는 약성 왜곡뿐 아니라 강성 왜곡도 고려되지 않았다. 하지만, 강성 왜곡은 오프셋과 동일한 형태의 오차로 센서 각 축의 출력값을 증감시키므로, Method C와 D에서도 강성 왜곡이 보정된다고 할 수 있겠다.

검증 실험에는 GY-87 센서모듈과 광학식 모션 캡쳐 시스템(Optitrack Flex 13)을 사용하였다. 여기서 GY-87 센서모듈은 6축 관성센서(MPU6050)와 3축 마그네토미터(HMC5883L)로 구성된 센서모듈이며, Arduino Uno R3에 연결하여 센서 측정값을 PC로 전송하였다. 광학식 모션 캡쳐 시스템은 방위각 추정 성능을 비교하기 위한 참조값(Reference)을 얻기 위해 사용되었다. 센서 측정값과 참조값을 동일 좌표계에서 얻기 위해 육면체 플라스틱 상자에 GY-87 센서모듈, Arduino Uno R3와 마커를 부착하여 실험하였다 (Fig. 1 참조). 또한, 광학식 모션 캡쳐 시스템과 GY-87 센서모듈의 계측 주기는 100 Hz로 동일하게 설정하였다.

Fig. 1 
Experimental setup

보정 성능 비교를 위해 무작위로 회전시켜 얻은 보정 데이터와 2가지 자세 실험을 실행하였다. 여기서 보정 데이터는 파라미터를 구하기 위해 사용되며, 2300개의 데이터를 사용하였다. 자세 실험은 구한 파라미터를 사용하여 방위각 추정 성능을 비교하기 위한 실험이다. 또한, 다양한 주변환경에 따른 성능 비교를 위해 보정 및 실험 데이터에 추가적으로 임의의 신호 왜곡을 가하였다.

방위각 추정 성능 비교를 위한 자세 실험은 다음과 같다.

(1) Test 1: 1분동안 센서의 Z축을 중심으로 약 45^o씩 회전시킨 후 5초간 정적 상태를 유지하였다.

(2) Test 2: 1분동안 센서를 무작위로 자세를 변경시킨 후 5초간 정적 상태를 유지하는 것을 반복하였다.

추가적인 신호 왜곡을 가하기 위하여 다음과 같이 마그네토미터 오차 파라미터를 임의로 설정하였다.

여기서, C는 민감도, 비직교성과 약성 왜곡에 의한 오차에 해당하는 행렬이며, b는 오프셋과 강성 왜곡에 의한 오차로 단위는 [mGauss]이다. 오프셋과 강성 왜곡만이 존재할 경우와 모든 오차 성분이 존재할 경우의 보정 기법 성능 비교를 위해 보정을 위한 데이터와 자세 실험 데이터에 다음과 같이 추가적인 왜곡을 가하였다.

여기서, B_orig는 추가적인 왜곡을 가하기 전의 원래(Original) 보정 및 실험 데이터이다. B_dist.1과 B_dist.2는 보정 및 실험 데이터에 임의의 왜곡을 가한 후(Distorted)로 각각 (1) 오프셋과 강성 왜곡만 왜곡을 가한 경우와 (2) 모든 오차 성분에 왜곡을 가한 경우의 데이터이다. 이에 맞추어, 자세실험 데이터 역시, B_dist.1를 적용한 (즉, b만 가한) Test 1-1과 Test 2-1과 B_dist.2를 적용한 (즉, b와 C를 가한) Test 1-2와 Test 2-2로 구분된다. 보정 기법간 방위각 추정 성능 비교를 위해 평균 제곱근 오차(Root Mean Square Error, RMSE)를 사용하였다.

Table 2는 각 실험에 따른 방위각 추정 RMSE 결과이다. 여기서, Method C′과 Method D′은 Method B에서 사용된 오차 모델링을 적용한 후의 결과이다. 또한, Figs. 2와 3은 Method A-D를 사용하였을 때 Test 1-2와 Test 2-2에서의 방위각 추정 결과 중 일부 데이터를 발췌하여 나타낸 그래프이다.

Fig. 2 
Results of azimuth estimation in Test 1-2

Fig. 3 
Results of azimuth estimation in Test 2-2

오프셋과 강성 왜곡만이 존재하는 Test 1-1과 Test 2-1에서 결과를 확인해보면, 보정 이전의 데이터를 사용하여 방위각을 추정한 경우 RMSE는 각각 88.61°, 104.09°로 정확성이 매우 저하되었다. Test 1-1에서 Method A를 사용한 경우 방위각 추정 RMSE가 1.11^o로 다른 보정 기법들에 비해 성능이 우수하였다. 반면, Method D에서 RMSE가 1.77°로 가장 큰 오차를 보였다. Test 2-1에서는 Method A-D′ 모두 방위각 추정 RMSE가 0.88^o 또는 0.89^o로 성능 차이를 보이지 않았다.

모든 오차가 존재하는 Test 1-2와 Test 2-2에서도 보정 이전의 데이터를 사용하여 방위각을 추정한 경우 RMSE가 96.34^o 이상으로 방위각 추정 성능이 매우 저하되었다. 모든 마그네토미터 오차를 고려한 Method A와 B는 Test 1-2에서는 1.97^o, Test 2-2에서는 1.53^o 이하의 오차를 보였다. 이에 비해 Method C는 방위각 추정 RMSE가 Test 1-2에서는 7.64^o, Test 2-2에서는 4.32^o로 이전과 비교하여 오차가 크게 증가하였다. Method C의 경우 오차 모델링 변경 이후에 Test 1-2와 Test 2-2에서 RMSE가 각각 3.62^o, 2.27^o로 이전과 비교하여 방위각 추정 성능이 향상되었다. 하지만 Method A와 B보다는 여전히 큰 오차를 보였다. Method D는 가장 큰 오차를 보였으며, 오차 모델링 변경 후에도 방위각 추정 성능이 향상되지 않았다. 이는 Method D에서 파라미터를 구하는데 있어 식(25)를 사용하여 발생하는 문제로 판단된다. 식(25)는 데이터마다 생기는 오차의 합으로 이루어져 있어, 이를 통해 파라미터를 구하게 될 경우 각 데이터의 오차는 최소화되지 않을 수 있다. 즉, 각 데이터의 오차는 남아있지만 오차를 더해주는 과정에서 상쇄되어 오차의 합이 최소가 되는 파라미터가 구해질 수 있다.

오프셋과 강성 왜곡만이 존재하는 경우에는 Method A-D′ 모두 방위각 추정 RMSE가 1^o 내외로 우수한 성능을 보이면서 방법 간 성능차이는 크지 않았다. 하지만, 모든 오차가 존재하는 경우 Method A와 B가 다른 방법들에 비해 우수한 추정 성능을 보였다. 특히, 전체적으로 Method A가 비교된 6가지 방법 중에서 가장 우수한 성능을 보였다. 또한, 식(26)의 추가적인 왜곡 성분을 약하게 설정하였을 경우에는 전체적으로 방위각 추정 오차가 감소하였다. 그러나 여전히 Method A와 B가 다른 방법들에 비해 우수한 성능을 보였다.

Method A는 Method B-D′와 달리 평면 또는 평면에 가까운 데이터를 사용하여 마그네토미터 오차 파라미터를 추정하는 경우 추정 정확성이 크게 저하된다. 예를 들어, 한 평면내의 측정값들로 제한된 데이터를 사용하여 파라미터를 추정한 경우, Method B는 Test 1-1에서 3.99^o, Test 2-1에서 4.74^o인 반면, Method A는 Test 1-1에서 12.79^o, Test 2-1에서 13.82^o의 RMSE를 보였다 –본 실험은 데이터에 따른 추정성능 민감도 확인을 위해 별도로 진행됨. 하지만, 마그네토미터를 3차원 공간상에서 적절히 회전시켜 주면서 계측을 하기만 하면 해결되는 문제이므로, 실제적인 운용에 있어 문제는 없다고 하겠다.

Method A는 우수한 추정성능에 더하여 선형해석과 파라미터수에서 다음과 같은 장점이 있다 (Table 1 참조). Method A를 제외한 나머지는 비선형 방정식을 풀어내야 하므로, Newton법 등과 같이 사용자가 지정하는 초기값에서 반복수정을 거치면 해를 구한다(Iterative Approach). 따라서, 지정된 초기값에 따라 수렴속도와 해에 많은 차이가 발생할 소지가 있다. 반면, Method A는 초기값 지정이 필요없는 선형해석과정을 통해 해가 도출되므로 해석속도는 물론이고, 해의 일관성과 정확성에서 보다 우수하다고 할 수 있다. 또한 보정 파라미터의 개수를 보면 Method B는 12개이고 나머지 방법들은 9개로 동일하지만, 고려되는 오차의 종류를 감안하면 Method A만 모든 오차를 고려하면서도 파라미터 개수가 9개인 것을 확인할 수 있다.

본 논문에서는 타원체 피팅 방식을 기반으로 하는 네 가지 마그네토미터 보정 기법 성능을 비교 분석하였다. 보정 성능 비교를 위해 두 가지 자세 실험에 따른 방위각 추정 성능을 비교하였으며, 다음과 같은 결론을 얻었다.

(1) 마그네토미터를 통해 방위각 추정을 함에 있어, 마그네토미터 자체가 가지는 오차뿐만 아니라 강성 왜곡과 약성 왜곡같이 주변환경에 의해 발생하는 오차도 고려해야 한다.

(2) 모든 실험에서 Method A가 가장 우수한 방위각 추정 성능을 보였다.

(3) 추정성능에 더하여 Method A는 파라미터의 개수가 상대적으로 적으며, 선형적인 특성을 지니기 때문에 다른 보정 기법과는 달리 초기값 설정이 필요하지 않다는 장점을 가지고 있다. 특히, 초기값이 불필요하다는 장점은 마그네토미터를 사용하며 수시로 보정이 필요한 경우, 실제 운용의 편의성 측면에서 매우 중요하다고 판단된다.