磁気測定システム
【課題】 動揺する環境下においても南北方向、東西方向、鉛直方向の磁界を精度良く測定する磁気測定システムを提供する。
【解決手段】 3軸磁気センサと、この３軸磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、3軸磁気センサのオフセット、直行度、感度及び姿勢センサと3軸磁気センサとの軸平行度のずれを、南北方向の磁界の変動と全磁力の変動が最小となるように補正し、これらの補正値を用いて3軸磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変更する。
【解決手段】 3軸磁気センサと、この３軸磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、3軸磁気センサのオフセット、直行度、感度及び姿勢センサと3軸磁気センサとの軸平行度のずれを、南北方向の磁界の変動と全磁力の変動が最小となるように補正し、これらの補正値を用いて3軸磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変更する。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、海上における船舶の磁気測定や環境磁気測定において、磁気センサが動揺していても、南北方向、東西方向、鉛直方向の磁界を精度良く測定することが可能な磁気測定システムに関する。
【背景技術】
【０００２】
従来、海上等の動揺する環境下において、南北方向、東西方向、鉛直方向の磁界を測定する場合、３軸磁気センサの動揺を姿勢センサを用いて測定し、この動揺情報を用いて、演算処理により南北方向、東西方向、鉛直方向に変換することができる変換式は次式のとおりである（例えば、非特許文献１参照）。
【数１】
【０００３】
ここで、Ｈxin：磁気センサＸ軸出力、ＨYin：磁気センサＹ軸出力、ＨZin：磁気センサＺ軸出力、ＨＸ：南北方向磁界、ＨＹ：東西方向磁界、ＨＺ：鉛直方向磁界、Ｐ：姿勢センサPitch角出力、Ｒ：姿勢センサRoll角出力、Ｙ：姿勢センサYaw角出力である。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】茂原正道著、「宇宙工学入門」、培風館、１９９４年１０月７日、ｐ．６２−６５
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来技術では、３軸磁気センサの各軸出力を南北方向、東西方向、鉛直方向の磁界を測定する場合、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれにより、数１をそのまま適用しても、東西方向、南北方向、鉛直方向の磁界を精度良く求めることは困難であった。
【０００６】
本発明は、そのような問題点に着目してなされたものであって、動揺する環境下においても南北方向、東西方向、鉛直方向の磁界を精度良く測定する磁気測定システムを提供することを目的としている。
【課題を解決するための手段】
【０００７】
前記課題を解決するために、この発明の磁気測定システムは、3軸磁気センサと、この磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正し、これらの補正値を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換することを特徴とする。
【０００８】
この磁気測定システムによれば、従来、動揺下における３軸磁気センサ出力と姿勢センサの出力を用いて、南北方向、東西方向、鉛直方向の磁界を求めても、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサの軸平行度のずれにより、南北方向、東西方向、鉛直方向の磁界を精度良く求めることが困難であったが、これらを求め、補正することにより、南北方向、東西方向、鉛直方向の磁界を精度良く演算処理により求めることができる。
【０００９】
具体的には、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを求め、これらを補正した後に上記数１を適用することにより、精度良く東西方向、南北方向、鉛直方向の磁界を求める。磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサの軸平行度のずれは、例えば南北方向の磁界の変動と全磁力の変動が最小となるように最急降下法やＧＡ（遺伝的アルゴリズム）等の最適パラメータ探査法にて求める。ここで、評価方法として、南北方向の磁界の変動を一例として挙げたが、東西方向や鉛直方向を評価しても良いし、それぞれを組み合わせて評価しても良い。
【発明の効果】
【００１０】
この発明によれば、従来困難であった動揺する環境下における南北方向、東西方向、鉛直方向の磁気を精度良く測定することが可能である。
【図面の簡単な説明】
【００１１】
【図１】実施形態に係る磁気測定システムの構成例の概略構成図である。
【図２】同磁気測定システムのセンサ筐体の内部の構成を示すブロック図である。
【図３】同磁気測定システムの処理概要を示すフロー図である。
【図４】同磁気測定システムに用いる３軸磁気センサの出力を示す図である。
【図５】同磁気測定システムに用いる３軸磁気センサの出力を数１を用いて南北方向、東西方向、鉛直方向の磁界に変換した結果を示す図である。
【図６】同磁気測定システムに用いる磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正するための補正係数をＧＡを用いて求めるフロー図である。
【発明を実施するための形態】
【００１２】
以下、実施の形態により、この発明を更に詳細に説明する。
図１は、この発明の実施形態に係る磁気測定システムの構成例を示す概略図である。センサ筐体１は、このセンサ筐体１にロープ４によって連結されたフロート２とおもり３によって、例えば海中に吊り下げられる。なお、符号５は海面を示している。
【００１３】
センサ筐体１の内部には、図２に示すように、Ｘ軸、Ｙ軸、Ｚ軸の各磁界を測定する３軸磁気センサ１０と、磁気センサ１０の動揺を測定する姿勢センサ１１と、磁気センサ１０及び姿勢センサ１１のアナログ信号をデジタル信号に変換するＡ／Ｄ１２と、例えば図６のフロー図に従って補正係数を算出する機能、求めた補正係数を用いて磁気センサ１０及び姿勢センサ１１の出力に数１を適用して南北方向、東西方向、鉛直方向の磁界を演算する機能等を有するＣＰＵ１３と、磁気センサ１０及び姿勢センサ１１の出力や補正係数等を格納するメモリ１４とが内蔵されている。
【００１４】
図３に処理概要を示す。メモリ１４に格納された磁気センサ１０の出力及び姿勢センサ１１の出力を補正係数計算に使用し（ステップＳＴ２１）、求めた補正係数を元に、ＣＰＵ１３で南北方向、東西方向、鉛直方向の磁界に変換する演算処理を行う（ステップＳＴ２２）。
【００１５】
図４に、３軸磁気センサの出力を示す。この図で明らかなように、３軸磁気センサは、波浪や潮流により動揺し、地磁気中で姿勢や方位が変化するため、Ｘ軸、Ｙ軸、Ｚ軸の各出力が大きく変動する。
【００１６】
図５に、３軸磁気センサの出力を上記数１を用いて南北方向、東西方向、鉛直方向の磁界に変換した結果を示す。数１の変換処理により、磁気センサの出力は、動揺下においても、南北方向、東西方向、鉛直方向の磁界に変換されるので、静止状態と等価の測定結果が得られるが、このままでは、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれにより、南北方向、東西方向、鉛直方向の磁界を精度良く求めることが困難となる。すなわち、図５の変換後のグラフに楕円印で示すように、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれが補正できていない場合、変動が大きくなる。
【００１７】
この変換後の変動を解決するための一例として、図６に、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正するための補正係数をＧＡを用いて求めるフロー図の例を示す。先ず、メモリ１４に格納された磁気センサ１０及び姿勢センサ１１の測定データを読み込む（ステップＳＴ１）。次に、ＧＡにより、次式数2、数3に示す磁気センサ１０のオフセット、直行度、感度及び磁気センサ１０と姿勢センサ１１との軸平行度の補正係数を設定する（ステップＳＴ２）。
【数２】
ここで、
Ｈｘ’〜Ｈｚ’ ＝磁気センサの測定データ
Ｈｘ〜Ｈｚ ＝オフセット、直交度、感度補正後 磁気データ
ａ11、ａ22、ａ33 ＝感度の補正係数
ａ12、ａ13、ａ21、ａ23、ａ31、ａ32 ＝直交度の補正係数
ａ14、ａ24、ａ34 ＝オフセットの補正係数
である。
【００１８】
（数３）
Ｐ＝Ｐ’＋Ｐｈ
Ｐ＝Ｒ’＋Ｒｈ
Ｙ＝Ｙ’＋Ｙｈ
ここで、
Ｐ’、Ｒ’、Ｙ’＝姿勢制御センサのデータ
Ｐｈ、Ｒｈ、Ｙｈ＝軸平行度の補正係数
Ｐ、Ｒ、Ｙ＝補正後の姿勢データ
である。
【００１９】
設定後、補正係数が有用であるかどうか評価する（ステップＳＴ３〜ＳＴ８）。具体的には、ステップＳＴ３にて、Ｘ軸、Ｙ軸、Ｚ軸の各出力から全磁力［（ＨＸ２＋ＨＹ２＋ＨＺ２）１／２］を計算し、ステップＳＴ４にて、全磁力のPeak to Peak値（差値）を求める。次いで、スッテプＳＴ５にて、例えば南北方向の磁界を計算し、ステップＳＴ６にて、南北方向の磁界変換後のPeak to Peak値を求める。続いて、スステップＳＴ７にて、適合度を計算する。すなわち、上記の全磁力のP-P値と南北方向のP-P値とが最小となるように適合度を計算し、ステップＳＴ８にて、収束条件を満たしているかどうか（最小となったかどうか）評価を行う。満たしていない場合は、ステップＳＴ９に進み、各Peak to Peak値が最小となるようにパラメータを変えるために補正係数をＧＡオペレータにて設定し直す。収束条件を満たした場合は、ステップＳＴ１０にて、求めた補正係数をメモリ１４に保存する。
【００２０】
以上の処理により、図３に示すように、補正係数を求め（ステップＳＴ２１）、補正後の補正係数に基づいて姿勢センサ及び磁気センサ出力に数１を適用することにより（ステップＳＴ２２）、南北方向、東西方向、鉛直方向の磁界を精度良く求めることができる。
【００２１】
なお、この例では、補正係数を求めるためにＧＡを用いたが、最急降下法やＳＡ（焼き鈍し法）等、他の最適パラメータ探査法を用いても良い。また、図３や図６のフロー図に示す処理は、センサ筐体１内のＣＰＵ１３で実行してもよいし、例えば海上の船舶内に設けたコンピュータにセンサ筐体１からの信号を取り込み、船舶内のコンピュータで計算を行ってもよい。
【符号の説明】
【００２２】
１ センサ筐体
２ フロート
３ おもり
４ ロープ
５ 海面
１０ 磁気センサ
１１ 姿勢センサ
１２ Ａ／Ｄ
１３ ＣＰＵ
１４ メモリ
【技術分野】
【０００１】
この発明は、海上における船舶の磁気測定や環境磁気測定において、磁気センサが動揺していても、南北方向、東西方向、鉛直方向の磁界を精度良く測定することが可能な磁気測定システムに関する。
【背景技術】
【０００２】
従来、海上等の動揺する環境下において、南北方向、東西方向、鉛直方向の磁界を測定する場合、３軸磁気センサの動揺を姿勢センサを用いて測定し、この動揺情報を用いて、演算処理により南北方向、東西方向、鉛直方向に変換することができる変換式は次式のとおりである（例えば、非特許文献１参照）。
【数１】
【０００３】
ここで、Ｈxin：磁気センサＸ軸出力、ＨYin：磁気センサＹ軸出力、ＨZin：磁気センサＺ軸出力、ＨＸ：南北方向磁界、ＨＹ：東西方向磁界、ＨＺ：鉛直方向磁界、Ｐ：姿勢センサPitch角出力、Ｒ：姿勢センサRoll角出力、Ｙ：姿勢センサYaw角出力である。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】茂原正道著、「宇宙工学入門」、培風館、１９９４年１０月７日、ｐ．６２−６５
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来技術では、３軸磁気センサの各軸出力を南北方向、東西方向、鉛直方向の磁界を測定する場合、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれにより、数１をそのまま適用しても、東西方向、南北方向、鉛直方向の磁界を精度良く求めることは困難であった。
【０００６】
本発明は、そのような問題点に着目してなされたものであって、動揺する環境下においても南北方向、東西方向、鉛直方向の磁界を精度良く測定する磁気測定システムを提供することを目的としている。
【課題を解決するための手段】
【０００７】
前記課題を解決するために、この発明の磁気測定システムは、3軸磁気センサと、この磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正し、これらの補正値を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換することを特徴とする。
【０００８】
この磁気測定システムによれば、従来、動揺下における３軸磁気センサ出力と姿勢センサの出力を用いて、南北方向、東西方向、鉛直方向の磁界を求めても、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサの軸平行度のずれにより、南北方向、東西方向、鉛直方向の磁界を精度良く求めることが困難であったが、これらを求め、補正することにより、南北方向、東西方向、鉛直方向の磁界を精度良く演算処理により求めることができる。
【０００９】
具体的には、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを求め、これらを補正した後に上記数１を適用することにより、精度良く東西方向、南北方向、鉛直方向の磁界を求める。磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサの軸平行度のずれは、例えば南北方向の磁界の変動と全磁力の変動が最小となるように最急降下法やＧＡ（遺伝的アルゴリズム）等の最適パラメータ探査法にて求める。ここで、評価方法として、南北方向の磁界の変動を一例として挙げたが、東西方向や鉛直方向を評価しても良いし、それぞれを組み合わせて評価しても良い。
【発明の効果】
【００１０】
この発明によれば、従来困難であった動揺する環境下における南北方向、東西方向、鉛直方向の磁気を精度良く測定することが可能である。
【図面の簡単な説明】
【００１１】
【図１】実施形態に係る磁気測定システムの構成例の概略構成図である。
【図２】同磁気測定システムのセンサ筐体の内部の構成を示すブロック図である。
【図３】同磁気測定システムの処理概要を示すフロー図である。
【図４】同磁気測定システムに用いる３軸磁気センサの出力を示す図である。
【図５】同磁気測定システムに用いる３軸磁気センサの出力を数１を用いて南北方向、東西方向、鉛直方向の磁界に変換した結果を示す図である。
【図６】同磁気測定システムに用いる磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正するための補正係数をＧＡを用いて求めるフロー図である。
【発明を実施するための形態】
【００１２】
以下、実施の形態により、この発明を更に詳細に説明する。
図１は、この発明の実施形態に係る磁気測定システムの構成例を示す概略図である。センサ筐体１は、このセンサ筐体１にロープ４によって連結されたフロート２とおもり３によって、例えば海中に吊り下げられる。なお、符号５は海面を示している。
【００１３】
センサ筐体１の内部には、図２に示すように、Ｘ軸、Ｙ軸、Ｚ軸の各磁界を測定する３軸磁気センサ１０と、磁気センサ１０の動揺を測定する姿勢センサ１１と、磁気センサ１０及び姿勢センサ１１のアナログ信号をデジタル信号に変換するＡ／Ｄ１２と、例えば図６のフロー図に従って補正係数を算出する機能、求めた補正係数を用いて磁気センサ１０及び姿勢センサ１１の出力に数１を適用して南北方向、東西方向、鉛直方向の磁界を演算する機能等を有するＣＰＵ１３と、磁気センサ１０及び姿勢センサ１１の出力や補正係数等を格納するメモリ１４とが内蔵されている。
【００１４】
図３に処理概要を示す。メモリ１４に格納された磁気センサ１０の出力及び姿勢センサ１１の出力を補正係数計算に使用し（ステップＳＴ２１）、求めた補正係数を元に、ＣＰＵ１３で南北方向、東西方向、鉛直方向の磁界に変換する演算処理を行う（ステップＳＴ２２）。
【００１５】
図４に、３軸磁気センサの出力を示す。この図で明らかなように、３軸磁気センサは、波浪や潮流により動揺し、地磁気中で姿勢や方位が変化するため、Ｘ軸、Ｙ軸、Ｚ軸の各出力が大きく変動する。
【００１６】
図５に、３軸磁気センサの出力を上記数１を用いて南北方向、東西方向、鉛直方向の磁界に変換した結果を示す。数１の変換処理により、磁気センサの出力は、動揺下においても、南北方向、東西方向、鉛直方向の磁界に変換されるので、静止状態と等価の測定結果が得られるが、このままでは、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれにより、南北方向、東西方向、鉛直方向の磁界を精度良く求めることが困難となる。すなわち、図５の変換後のグラフに楕円印で示すように、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれが補正できていない場合、変動が大きくなる。
【００１７】
この変換後の変動を解決するための一例として、図６に、磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正するための補正係数をＧＡを用いて求めるフロー図の例を示す。先ず、メモリ１４に格納された磁気センサ１０及び姿勢センサ１１の測定データを読み込む（ステップＳＴ１）。次に、ＧＡにより、次式数2、数3に示す磁気センサ１０のオフセット、直行度、感度及び磁気センサ１０と姿勢センサ１１との軸平行度の補正係数を設定する（ステップＳＴ２）。
【数２】
ここで、
Ｈｘ’〜Ｈｚ’ ＝磁気センサの測定データ
Ｈｘ〜Ｈｚ ＝オフセット、直交度、感度補正後 磁気データ
ａ11、ａ22、ａ33 ＝感度の補正係数
ａ12、ａ13、ａ21、ａ23、ａ31、ａ32 ＝直交度の補正係数
ａ14、ａ24、ａ34 ＝オフセットの補正係数
である。
【００１８】
（数３）
Ｐ＝Ｐ’＋Ｐｈ
Ｐ＝Ｒ’＋Ｒｈ
Ｙ＝Ｙ’＋Ｙｈ
ここで、
Ｐ’、Ｒ’、Ｙ’＝姿勢制御センサのデータ
Ｐｈ、Ｒｈ、Ｙｈ＝軸平行度の補正係数
Ｐ、Ｒ、Ｙ＝補正後の姿勢データ
である。
【００１９】
設定後、補正係数が有用であるかどうか評価する（ステップＳＴ３〜ＳＴ８）。具体的には、ステップＳＴ３にて、Ｘ軸、Ｙ軸、Ｚ軸の各出力から全磁力［（ＨＸ２＋ＨＹ２＋ＨＺ２）１／２］を計算し、ステップＳＴ４にて、全磁力のPeak to Peak値（差値）を求める。次いで、スッテプＳＴ５にて、例えば南北方向の磁界を計算し、ステップＳＴ６にて、南北方向の磁界変換後のPeak to Peak値を求める。続いて、スステップＳＴ７にて、適合度を計算する。すなわち、上記の全磁力のP-P値と南北方向のP-P値とが最小となるように適合度を計算し、ステップＳＴ８にて、収束条件を満たしているかどうか（最小となったかどうか）評価を行う。満たしていない場合は、ステップＳＴ９に進み、各Peak to Peak値が最小となるようにパラメータを変えるために補正係数をＧＡオペレータにて設定し直す。収束条件を満たした場合は、ステップＳＴ１０にて、求めた補正係数をメモリ１４に保存する。
【００２０】
以上の処理により、図３に示すように、補正係数を求め（ステップＳＴ２１）、補正後の補正係数に基づいて姿勢センサ及び磁気センサ出力に数１を適用することにより（ステップＳＴ２２）、南北方向、東西方向、鉛直方向の磁界を精度良く求めることができる。
【００２１】
なお、この例では、補正係数を求めるためにＧＡを用いたが、最急降下法やＳＡ（焼き鈍し法）等、他の最適パラメータ探査法を用いても良い。また、図３や図６のフロー図に示す処理は、センサ筐体１内のＣＰＵ１３で実行してもよいし、例えば海上の船舶内に設けたコンピュータにセンサ筐体１からの信号を取り込み、船舶内のコンピュータで計算を行ってもよい。
【符号の説明】
【００２２】
１ センサ筐体
２ フロート
３ おもり
４ ロープ
５ 海面
１０ 磁気センサ
１１ 姿勢センサ
１２ Ａ／Ｄ
１３ ＣＰＵ
１４ メモリ
【特許請求の範囲】
【請求項１】
3軸磁気センサと、この磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、
磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正し、これらの補正値を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換することを特徴とする磁気測定システム。
【請求項２】
前記磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれは、南北方向、東西方向、鉛直方向のうちの少なくとも１方向の磁界の変動と全磁力の変動とが最小となるように補正することを特徴とする請求項１記載の磁気測定システム。
【請求項１】
3軸磁気センサと、この磁気センサの動揺を測定する姿勢センサとを備え、姿勢センサにより得られた動揺情報を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換する磁気測定システムにおいて、
磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれを補正し、これらの補正値を用いて磁気センサの測定データを南北方向、東西方向、鉛直方向の磁界に変換することを特徴とする磁気測定システム。
【請求項２】
前記磁気センサのオフセット、直行度、感度及び姿勢センサと磁気センサとの軸平行度のずれは、南北方向、東西方向、鉛直方向のうちの少なくとも１方向の磁界の変動と全磁力の変動とが最小となるように補正することを特徴とする請求項１記載の磁気測定システム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図２】
【図３】
【図４】
【図５】
【図６】
【公開番号】特開２０１１−４７９１０（Ｐ２０１１−４７９１０Ａ）
【公開日】平成２３年３月１０日（２０１１．３．１０）
【国際特許分類】
【出願番号】特願２００９−１９９０２５（Ｐ２００９−１９９０２５）
【出願日】平成２１年８月２８日（２００９．８．２８）
【出願人】（３９００１４３０６）防衛省技術研究本部長 (169)
【出願人】（０００００１９９３）株式会社島津製作所 (3,708)
【Ｆターム（参考）】
【公開日】平成２３年３月１０日（２０１１．３．１０）
【国際特許分類】
【出願日】平成２１年８月２８日（２００９．８．２８）
【出願人】（３９００１４３０６）防衛省技術研究本部長 (169)
【出願人】（０００００１９９３）株式会社島津製作所 (3,708)
【Ｆターム（参考）】
[ Back to top ]