Software Gimbals for your compass

If 12 factor calibration is used on both accelerometer and magnetometer data, and the Rodrigues formula is applied do the magnetic dip statistics improve?

Code to apply 12 factor calibrations ( Using a LSM9DS1 IMU )

//Use 12 element calibration on raw count magnetometer data
function cal8G(S){
 var i;
 var mm=[0,0,0];
 var scale=8;
 var mres=[
     [0.970972836,0.171154738,-0.000344029],
     [-0.095223453,0.956653319,-0.000125418],
     [0.017179224,0.005190024,1],
     ];
 var moff=[1057.176483,120.7727841,549.551664];
 for(i=0;i<3;i++) S.m[i] -=moff[i];
 for(i=0;i<3;i++)
  mm[i]=S.m[0]*mres[i][0]+S.m[1]*mres[i][1]+S.m[2]*mres[i][2];
  for(i=0;i<3;i++)S.m[i]=mm[i]*scale/32768;
}//end cal8G(S)
//Use 12 element calibration on raw count accelerometer data
function calA2(S){
 var i;
 var mm=[0,0,0];
 var scale=2;
 var mres=[
   [1.005708319,-0.007635345,-1.55429E-05],
   [0.001530038,1.006717582,-9.04159E-05],
   [0.017299441,0.097632586,1]
   ];
 var moff=[-126.0318135,-350.2409885,730.2274842];

 for(i=0;i<3;i++) S.a[i] -=moff[i];
 for(i=0;i<3;i++)
  mm[i]=S.a[0]*mres[i][0]+S.a[1]*mres[i][1]+S.a[2]*mres[i][2];
  for(i=0;i<3;i++)S.a[i]=mm[i]*scale/32768;
}//end cal8G(S)
////////////////////
function readall(W){
if(W.magAvailable()){
  W.readMag();
  W.readAccel();
  cal8G(W);
 calA2(W);
console.log(W.a,',',W.m);
//console.log(W.a);
}//endif
}//end readall

Data were collected and pasted into the following code that is used to determine the signs of the magnetometer axis that produce the best magnetic dip angle statistics.

function test(){
 FF=require("Rodrigues").connect();
 var i,j,k,l;
 var d="";
  var t=[0,0,0];
  console.log("mmm,mmp,mpm,mpp,pmm,pmp,ppm,ppp");
 for(i=0;i<Adata.length;i++){
  d=" "+i;
 A=FF.Setup(Adata[i],[0,0,1]);
 Theta=FF.theta;
 Kn=FF.Kn;
 A=FF.An;
for(j=-1;j<2;j=j+2){
 for(k=-1;k<2;k=k+2){
  for(l=-1;l<2;l=l+2){
   t[0]=Mdata[i][0]*j;
   t[1]=Mdata[i][1]*k;
   t[2]=Mdata[i][2]*l;
   B=FF.RotVector(t);
   An=FF.An;
 dip=pirate*Math.atan(B[2]/Math.sqrt(B[0]*B[0]+B[1]*B[1]));
     d=d+","+dip;
  }//nextl
 }//nextk
}//nextj
console.log(d);
}//next i
}

This produces output that looks like the following:

This output is copied into a CSV file and the CSV file is opened using Excel. Descriptive Statistics are then applied to the columns of data

The mpp and pmm columns produce the smallest coffoecoent of variation (COV) values. The following code uses the +X, -Y, -Z (pmm) combination.

function test(){
 FF=require("Rodrigues").connect();
 var i;
 for(i=0;i<Adata.length;i++){
 A=FF.Setup(Adata[i],[0,0,1]);
   Theta=FF.theta;
   Kn=FF.Kn;
 A=FF.An;
   t[0]=Mdata[i][0]*1;
   t[1]=Mdata[i][1]*-1;
   t[2]=Mdata[i][2]*-1;
  B=FF.RotVector(t);
   An=FF.An;
  heading=pirate*Math.atan2(-B[1],B[0]);
 if(heading<0)heading+=360;
 dip=pirate*Math.atan(B[2]/Math.sqrt(B[0]*B[0]+B[1]*B[1]));
console.log(i,',',A,',',An,',',B,',',
 Theta,',',heading,',',dip);
}//next i
}

This produces the following: