技術特征：

1.一種極光測量儀器控制方法，其特征在于，包括：

接收衛(wèi)星廣播的定位數(shù)據(jù)；

根據(jù)所述定位數(shù)據(jù)中的時間參數(shù)計算當前時刻地球磁場的磁矩三分量，基于所述磁矩三分量計算得到由地理坐標系到地磁坐標系的坐標轉換參數(shù)；

根據(jù)所述定位數(shù)據(jù)中的地理坐標、所述坐標轉換參數(shù)計算地磁坐標，在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作。

2.根據(jù)權利要求1所述的極光測量儀器控制方法，其特征在于，所述根據(jù)所述定位數(shù)據(jù)中的時間參數(shù)計算當前時刻地球磁場的磁矩三分量包括：

根據(jù)所述定位數(shù)據(jù)中的時間參數(shù)計算當前時刻相對于磁矩更新起始時間的年數(shù)；

根據(jù)磁矩三分量的初始值、磁矩三分量的年變化值以及所述年數(shù)計算當前時刻地球磁場的磁矩三分量。

3.根據(jù)權利要求2所述的極光測量儀器控制方法，其特征在于，計算當前時刻相對于磁矩更新起始時間的年數(shù)的計算式為：

ΔT＝(D_GPS+T_GPS/1000.0/86400.0-EPOCH)/365.25；

其中，ΔT表示相對于磁矩更新起始時間的年數(shù)，定位數(shù)據(jù)中的時間參數(shù)D_GPS為J2000坐標系的天數(shù)、T_GPS為J2000坐標系的毫秒數(shù)，EPOCH表示磁矩更新起始時間，365.25表示在一個更新周期內每年的平均天數(shù)；

根據(jù)以下計算式計算當前時刻地球磁場的磁矩三分量：

G10＝G100+ΔG10×ΔT；

G11＝G110+ΔG11×ΔT；

H11＝H110+ΔH11×ΔT；

其中，G10、G11、H11表示當前時刻的磁矩三分量，G100、G110、H110表示在地理坐標系中磁矩三分量的初始值，ΛG10、ΛG11、ΛH11表示在地理坐標系中磁矩三分量的年變化值；

根據(jù)以下計算式計算所述坐標轉換參數(shù)：

SQ＝G11²+H11²；

$SQQ=\sqrt{SQ};$

$SQR=\sqrt{G{10}^2+SQ};$

SL0＝-H11/SQQ；

CL0＝-G11/SQQ；

ST0＝SQQ/SQR；

CT0＝-G10/SQR；

STCL＝ST0*CL0；

STSL＝ST0*SL0；

CTSL＝CT0*SL0；

CTCL＝CT0*CL0；

所述根據(jù)所述定位數(shù)據(jù)中的地理坐標、所述坐標轉換參數(shù)計算地磁坐標包括：根據(jù)以下計算式計算地磁坐標：

X_M＝X_G×CTCL+Y_G×CTSL-Z_G×ST0；

Y_M＝Y_G×CL0-X_G×SL0；

Z_M＝X_G×STCL+Y_G×STSL+Z_G×CT0；

其中，X_M、Y_M、Z_M表示地磁坐標，X_G、Y_G、Z_G表示地理坐標。

4.根據(jù)權利要求3所述的極光測量儀器控制方法，其特征在于，所述在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作包括：

根據(jù)如下計算式計算磁緯角度：

$\mathrm{\θ}=\arcsin \sqrt{\frac{{Z_M}^2}{{X_M}^2+{Y_M}^2+{Z_M}^2}};$

當所述磁緯角度滿足預設條件時觸發(fā)所述極光測量儀器開始工作。

5.根據(jù)權利要求3所述的極光測量儀器控制方法，其特征在于，所述在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作包括：

計算1/sin²θ＝(X_M/Z_M)²+(Y_M/Z_M)²+1，當(X_M/Z_M)²+(Y_M/Z_M)²滿足預設條件時觸發(fā)所述極光測量儀器開始工作。

6.一種極光測量儀器控制系統(tǒng)，其特征在于，包括：

接收模塊，用于接收衛(wèi)星廣播的定位數(shù)據(jù)；

計算模塊，用于根據(jù)所述定位數(shù)據(jù)中的時間參數(shù)計算當前時刻地球磁場的磁矩三分量，基于所述磁矩三分量計算得到由地理坐標系到地磁坐標系的坐標轉換參數(shù)；

控制模塊，用于根據(jù)所述定位數(shù)據(jù)中的地理坐標、所述坐標轉換參數(shù)計算地磁坐標，在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作。

7.根據(jù)權利要求6所述的極光測量儀器控制系統(tǒng)，其特征在于，所述計算模塊包括：

年數(shù)計算單元，用于根據(jù)所述定位數(shù)據(jù)中的時間參數(shù)計算當前時刻相對于磁矩更新起始時間的年數(shù)；

磁矩分量計算單元，用于根據(jù)磁矩三分量的初始值、磁矩三分量的年變化值以及所述年數(shù)計算當前時刻地球磁場的磁矩三分量；

參數(shù)計算單元，用于基于所述磁矩三分量計算得到由地理坐標系到地磁坐標系的坐標轉換參數(shù)。

8.根據(jù)權利要求7所述的極光測量儀器控制系統(tǒng)，其特征在于，

所述年數(shù)計算單元具體用于根據(jù)以下計算式計算當前時刻相對于磁矩更新起始時間的年數(shù)：

ΔT＝(D_GPS+T_GPS/1000.0/86400.0-EPOCH)/365.25；

其中，ΔT表示相對于磁矩更新起始時間的年數(shù)，定位數(shù)據(jù)中的時間參數(shù)D_GPS為J2000坐標系的天數(shù)、T_GPS為J2000坐標系的毫秒數(shù)，EPOCH表示磁矩更新起始時間，365.25表示在一個更新周期內每年的平均天數(shù)；

所述磁矩分量計算單元具體用于根據(jù)以下計算式計算當前時刻地球磁場的磁矩三分量：

G10＝G100+ΔG10×ΔT；

G11＝G110+ΔG11×ΔT；

H11＝H110+ΔH11×ΔT；

其中，G10、G11、H11表示當前時刻的磁矩三分量，G100、G110、H110表示在地理坐標系中磁矩三分量的初始值，ΛG10、ΛG11、ΛH11表示在地理坐標系中磁矩三分量的年變化值；

所述參數(shù)計算單元具體用于根據(jù)以下計算式計算所述坐標轉換參數(shù)：

SQ＝G11²+H11²；

$SQQ=\sqrt{SQ};$

$SQR=\sqrt{G{10}^2+SQ};$

SL0＝-H11/SQQ；

CL0＝-G11/SQQ；

ST0＝SQQ/SQR；

CT0＝-G10/SQR；

STCL＝ST0*CL0；

STSL＝ST0*SL0；

CTSL＝CT0*SL0；

CTCL＝CT0*CL0；

所述控制模塊具體用于根據(jù)以下計算式計算地磁坐標：

X_M＝X_G×CTCL+Y_G×CTSL-Z_G×ST0；

Y_M＝Y_G×CL0-X_G×SL0；

Z_M＝X_G×STCL+Y_G×STSL+Z_G×CT0；

其中，X_M、Y_M、Z_M表示地磁坐標，X_G、Y_G、Z_G表示地理坐標。

9.根據(jù)權利要求8所述的極光測量儀器控制系統(tǒng)，其特征在于，所述控制模塊用于在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作包括：

所述控制模塊具體用于根據(jù)如下計算式計算磁緯角度：

$\mathrm{\θ}=arcsin\sqrt{\frac{{Z_M}^2}{{X_M}^2+{Y_M}^2+{Z_M}^2}};$

當所述磁緯角度滿足預設條件時觸發(fā)所述極光測量儀器開始工作。

10.根據(jù)權利要求8所述的極光測量儀器控制系統(tǒng)，其特征在于，所述控制模塊用于在地磁坐標滿足預設條件時觸發(fā)極光測量儀器開始工作包括：

所述控制模塊具體用于計算1/sin²θ＝(X_M/Z_M)²+(Y_M/Z_M)²+1，當(X_M/Z_M)²+(Y_M/Z_M)²滿足預設條件時觸發(fā)所述極光測量儀器開始工作。