(57)【要約】  本発明は、複数のスペクトル画像のセットを変換するための方法に関する。当該方法は、当該セットにおける複数画像を等しく配置された複数区域に分割するステップと、当該複数区域の各々に関して、当該複数画像のセット中の予め定められた特徴を計算するステップと、当該複数画像の各々に関して、当該複数区域の各々における複数の輝度値を当該区域の当該予め定められた特徴の関数で正規化するステップとを含む。本発明はまた、対応するコンピュータプログラムプロダクトおよび対応する画像処理システムに関する。
2018521418.tif 000002
【特許請求の範囲】
【請求項１】
複数の二次元スペクトル画像のセットを変換するための、コンピュータによって実現される方法であって、
−  前記セットにおける複数画像を等しく配置された複数区域に分割するステップ（２１０）と、
−  前記複数区域の各々に関して、前記複数画像のセット中の予め定められた特徴を計算するステップ（２２０）と、
−  前記複数画像の各々に関して、前記複数区域の各々における複数の輝度値を前記区域の前記予め定められた特徴の関数で正規化するステップ（２３０）とを含む、コンピュータによって実現される方法。
【請求項２】
前記予め定められた特徴は平均輝度であり、前記正規化するステップは、前記複数区域の各々における前記複数の輝度値を前記平均輝度値に対して正規化するステップを含む、請求項１に記載の、コンピュータによって実現される方法。
【請求項３】
前記複数区域は個々の画素に対応する、請求項１または２に記載の、コンピュータによって実現される方法。
【請求項４】
前記複数区域は、別個の波長帯域を表わすそれぞれの複数の画素を含む複数の矩形ブロックに対応する、請求項１または２に記載の、コンピュータによって実現される方法。
【請求項５】
前掲の請求項のいずれかに記載の方法のステップを実行するようプロセッサに指示するように構成されたコード手段を含むコンピュータプログラムプロダクト。
【請求項６】
複数のスペクトル画像を変換するためのシステム（３００）であって、
−  複数の二次元スペクトル画像のセットを受取るように適合された入力手段（３１０）と、
−  処理手段（３２０）とを含み、前記処理手段（３２０）は、
−  前記セットにおける複数画像を等しく配置された複数区域に分割し、
−  前記複数区域の各々に関して、前記複数画像のセット中の予め定められた特徴を計算し、
−  前記複数画像の各々に関して、前記複数区域の各々における複数の輝度値を前記区域の前記予め定められた特徴の関数で正規化するように構成されており、前記システム（３００）はさらに、
−  前記処理手段（３２０）によって処理されるとおりに前記複数画像のセットを出力するように適合された出力手段（３３０）を含む、システム（３００）。
【請求項７】
前記予め定められた特徴は平均輝度であり、前記正規化することは、前記複数区域の各々における前記複数の輝度値を前記平均輝度値に対して正規化することを含む、請求項６に記載のシステム。
【請求項８】
前記複数区域は個々の画素に対応する、請求項７に記載のシステム。
【請求項９】
前記複数区域は、別個の波長帯域を表わすそれぞれの複数の画素を含む複数の矩形ブロックに対応する、請求項７に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
発明の分野
本発明は、たとえば、航空映像において得られる画像に適用されるべき画像処理の分野に関する。より特定的には、本発明はスペクトル画像を変換するための方法およびシステムに関する。
【背景技術】
【０００２】
背景
マルチスペクトル撮像およびハイパースペクトル撮像の分野においては、近年、画像取得のための従来の「プッシュブルーム（pushbroom）」法の代替案が提案されてきた。この代替案は、センサのさまざまな部分をさまざまな波長に対して高感度にする薄膜スペクトルフィルタを備えた二次元センサを用いることを含む。送信された波長の変化は、典型的には、センサが背景全体にわたって動かされる方向（すなわち、センサ構成が搭載された航空機または宇宙航空機の場合には飛行方向）に沿って徐々にまたは段階的に起こる。したがって、このようなセンサによって得られる各々の画像は実際にはモザイク画であって、画像のうち別々の部分が別々のそれぞれの波長帯域の放射において見られるように地形の対応する部分を表現している。「スペクトル画像」という語は、以下では、撮像された放射をそれぞれ別個の波長で表現する別個の領域を含む画像を表わすのに用いられる。
【０００３】
このタイプのセンサは、ＩＭＥＣの名義で「スペクトル撮像システムのための集積回路（Integrated Circuit for Spectral Imaging System）」と題された国際特許出願公報ＷＯ２０１１／０６４４０３Ａ１に開示されている。
【０００４】
フィルタベースのマルチスペクトルセンサの使用に関してそれ自体が呈する深刻な問題は、取得したマルチスペクトル画像を幾何学的に参照および位置合わせすることである。
【０００５】
スペクトル画像を参照する際の問題に対処する一方法は、スペクトル画像に対して既知の一定した関係をもつ一連のパンクロマティック画像を同時に取得することである。このような同時取得のために適合されたセンサ構成が、本出願人の名義で「マルチスペクトルデータの幾何学的参照（Geometric Referencing of Multi-Spectral Data）」と題された国際特許出願公報ＷＯ２０１１／０７３４３０Ａ１において詳細に開示されている。この解決策は、検知装置によって生成されるデータ量を増加させてしまう。このデータは、検知装置からどうにかして画像処理インフラストラクチャに送信し返さなければならないものである。このような送信は多大な費用および／または多大な時間が掛かる可能性がある。
【０００６】
したがって、追加のパンクロマティック画像情報に依拠しないマルチスペクトル画像およびハイパースペクトル画像を参照して位置合わせする方法が必要とされる。
【発明の概要】
【課題を解決するための手段】
【０００７】
発明の概要
本発明の一局面に従うと、複数のスペクトル画像のセットを変換するための、コンピュータによって実現される方法が提供される。当該方法は、上記セットにおける複数画像を等しく配置された複数区域に分割するステップと、上記複数区域の各々に関して、上記複数画像のセット中の予め定められた特徴を計算するステップと、上記複数画像の各々に関して、上記複数区域の各々における複数の輝度値を上記区域の上記予め定められた特徴の関数で正規化するステップとを含む。
【０００８】
上記複数区域に関して、複数の輝度値についての１つ以上の代表的な特徴を計算することができる。１つのこのような特徴として、１区域にわたる平均輝度値が挙げられる。別の有用な特徴として、測定されるべきコントラストを示す輝度値の標準偏差が挙げられる。より一般的には、輝度値の分布を計算し、より大きな１セットの特徴で表現することができるだろう。１区域当たりに得られる１セットの特徴は正規化係数として用いることができる。これらの特徴を用いて正規化を適用した後、これらの特徴の値は、結果として得られた画像におけるさまざまな区域にわたって均一になる。
【０００９】
正規化係数を決定するための手順は、画像内容の影響を平均化するために十分に大きな１セットの複数画像を平均化することによって行なわれる。その後、同じ複数画像について、または、同じ機器を用いて同じ方法で取得された他の複数画像について、設定された係数を用いて正規化を行なうことができる。この手順により、すべての新しいセットの画像についての新しい係数を計算する必要がなくなるため、作業方法が簡素化される。
【００１０】
本発明は、特に、フィルタがあるせいでさまざまな波長帯域に反応してしまうセンサ内のさまざまな部分によって物理的特徴が取得されているので、これら物理的特徴が連続取得された異なる複数画像において異なる態様で表示されるという事実により、さまざまなスペクトル画像の相互位置合わせがより困難になるという発明者の見識に基づいている。
【００１１】
本発明は、特に、さまざまな波長帯域で物理的特徴を表現しており取得順が同じであるさまざまなスペクトル画像の間の所与の物理的特徴の輝度差に対して２つの成分が存在する、という発明者の見識にさらに基づいている。（１）物理的特徴はさまざまな波長帯域で異なる反射率を有し得る。（２）センサはさまざまな波長帯域で異なる感度を有し得る。第２の要因は、それぞれの部分を表わす平均値に対して複数画像のうちさまざまな部分を正規化することによって補償することができる。
【００１２】
第１の要因を補償することは不可能であるが、発明者らは、意外にも、第２の要因だけを補償した後に位置合わせアルゴリズムの効率が既に大幅に改善していることを発見した。この効果は、現実世界の物理的オブジェクトが、典型的には、対象となるスペクトルの大部分にわたる波長の関数で徐々に変化する反射率を呈するという事実によるものと考えられている。
【００１３】
本発明に従った方法の一実施形態においては、予め定められた特徴は平均輝度であり、正規化することは、複数区域の各々における輝度値を平均輝度値に対して正規化することを含む。
【００１４】
本発明に従った方法の一実施形態においては、複数区域は個々の画素に対応する。
この実施形態の利点は、個々の画素・フィルタの組合せの感度が変化した理由を、（フィルタの製造公差または不純物を含む）このような変化の原因にかかわらず、説明することができるように、センサが１画素単位で有効に較正される点である。これによりアーティファクトが最大限に抑制されることとなる。画素・フィルタの組合せに光学系を追加することにより、完全な撮像システムが得られる。光学系によってもたらされる感度変化を含めるようにしてこれら感度変化を補正することを選択することもでき、または、感度変化を排除することでシステムをさまざまな光学系にとって汎用的なままにしておくことを選択することもできる。
【００１５】
本発明に従った方法の一実施形態においては、複数区域は、別個の波長帯域を表わす複数のそれぞれの画素を含む複数の矩形ブロックに対応する。
【００１６】
複数画素のブロックごとに正規化を実行できることがこの実施形態の利点である。この場合、ブロックは、典型的には、矩形細片のセンサまたは複数の矩形区域の組合せを表わしている。
【００１７】
本発明の一局面に従うと、上述の方法のステップを実行するようプロセッサに指示するように構成されたコード手段を含むコンピュータプログラムプロダクトが提供される。
【００１８】
本発明の一局面に従うと、複数のスペクトル画像を変換するためのシステムが提供される。当該システムは、複数のスペクトル画像のセットを受取るように適合された入力手段と、処理手段とを含む。処理手段は、上記セットにおける複数画像を等しく配置された複数区域に分割し、上記複数区域の各々に関して、上記複数画像のセット中の予め定められた特徴を計算し、上記複数画像の各々に関して、上記複数区域の各々における複数の輝度値を上記区域の上記予め定められた特徴の関数で正規化するように構成されている。当該システムはさらに、上記処理手段によって処理されるように上記複数画像のセットを出力するように適合された出力手段を含む。
【００１９】
本発明に従ったシステムの一実施形態においては、予め定められた特徴は平均輝度であり、正規化することは、複数区域の各々における複数の輝度値を平均輝度値に対して正規化することを含む。
【００２０】
本発明に従ったシステムの一実施形態においては、複数区域は個々の画素に対応する。
本発明に従ったシステムの一実施形態においては、複数区域は、別個の波長帯域を表わすそれぞれの複数の画素を含む複数の矩形ブロックに対応する。
【００２１】
本発明に従ったコンピュータプログラムプロダクトおよびシステムについての実施形態の技術的な作用および利点は、本発明に従った方法の対応する実施形態についての作用および利点に必要な変更を加えたものに対応している。
【００２２】
図の簡単な説明
本発明の実施形態のこれらおよび他の技術的な作用および利点を、ここで、添付の図面に関連付けてより詳細に説明する。
【図面の簡単な説明】
【００２３】
【図１】マルチスペクトルセンサ、特にハイパースペクトルセンサの連続取得によって撮像された領域を示す斜視図である。
【図２】本発明に従った方法の一実施形態を示すフローチャートである。
【図３】本発明に従ったシステムの一実施形態を概略的に示す図である。
【図４】本発明に従った方法の別の実施形態を示すフローチャートである。
【発明を実施するための形態】
【００２４】
実施形態の説明
図１は、マルチスペクトルセンサ、特にハイパースペクトルセンサの連続取得によって撮像された領域の斜視図を示す。ハイパースペクトルセンサの一例が、本出願人の名義で前述の国際特許出願公報ＷＯ２０１１／０７３４３０Ａ１において開示されており、第２の（視程）センサと連携して動作する「第１のセンサ」として記載されている。ＷＯ２０１１／０７３４３０Ａ１の「第１のセンサ」は、本発明を明確にするために参照されているが、本発明がこれに限定されないものと理解されるべきである。
【００２５】
検知素子のさまざまな部分がさまざまな波長に対して高感度であることは、このようなハイパースペクトルセンサにとっては特有のことである。この効果は、検知素子に、当該検知素子の表面にわたって変化する波長応答を有するフィルタリング層を備えることによって得られる可能性がある。したがって、このようなハイパースペクトルセンサによって得られた各画像は実際にはモザイク画であって、画像のうち別々の部分が、それぞれ別々の波長帯域の放射に見られるように地形の対応する部分を表現している。１つの特定の波長帯域の放射に見られるような任意の所与の区域の画像を得るために、多数のこのようなモザイク画のうち関連する部分同士を継ぎ合わせなければならない。すべての関連するスペクトル帯域において充分な空間カバレッジを確実にするために、これらのハイパースペクトルセンサが（センサ移動の速度に応じて、非常に高いフレームレートを必要とする可能性がある）密集配置された画像を必要とすることは明らかである。
【００２６】
発明者らは、任意の所与の波長帯域で画像同士を「継ぎ合わせる」のに必要とされるこれらのスペクトル画像を幾何学的に参照する効率が予備的な正規化ステップを実行することによって大幅に改善されることを発見した。
【００２７】
図２は、本発明に従った方法の一実施形態のフローチャートを示す。
第１のステップにおいて、複数のスペクトル画像のセットを変換するための方法２００の例示的な実施形態は、上記セットにおける複数画像を等しく配置された複数区域に分割するステップ（２１０）を含む。これら区域は、個々の画素、同じ波長帯域に対して感度の高い画素の細片、または、画素の他のいずれかのグループ化に対応し得る。
【００２８】
複数区域の各々に関して、複数画像のセットにわたって平均輝度値が計算される（２２０）。検出器応答は、露光時にセンサのそれぞれの画素によって検知されるような（たとえば、ＣＣＤセンサまたはＣＭＯＳセンサ上に光子が入射した結果蓄積された電荷の量で表わされるような）画素値であり、かつ、通常、画素の輝度を表わすスカラー（グレイスケール）値である。「平均」は、演算平均もしくは幾何学平均、適切に選択されたパーセンタイル値、またはメジアン値などの輝度値の適切な統計値であってもよい。このステップは１回の動作として示されているが、当業者であれば、この動作が、たとえば反復してまたは並列に実行され得ることを理解するだろう。
【００２９】
複数画像の各々に関して、複数区域の各々における輝度値が当該区域の平均に対して正規化される（２３０）。この正規化は、対象となっている画素が属する区域の平均値で各スカラー画素値を割ることを含んでもよい。このステップは１回の動作として示されているが、当業者であれば、この動作が、たとえば反復してまたは並列に実行され得ることを理解するだろう。
【００３０】
正規化ステップの結果、完全にフラットなスペクトル反応を有する（すなわち、対象となるスペクトルにわたって真っ白である）仮想オブジェクトが正規化された画像の各々において等しい輝度で表現され得る可能性がある。これにより、結果として得られる画像はシミュレートされたパンクロマティック画像として作用し得る。
【００３１】
上記セットの複数画像が正規化されると、これら画像は、以下に例示されるように、任意の所望の処理アルゴリズム（図示せず）で入力として用いられ得る。
【００３２】
図３は、本発明に従ったシステムの一実施形態を概略的に示す。図示されたシステム３００は、複数のスペクトル画像のセットを受取るように適合された入力手段３１０を含む。
【００３３】
システム３００は処理手段３２０をさらに含む。処理手段３２０は、
−  当該セットにおける複数画像を等しく配置された複数区域に分割し、
−  複数区域の各々に関して、複数画像のセットにわたる平均輝度値を計算し、
−  複数画像の各々に関して、複数区域の各々における複数の輝度値を当該区域の平均に対して正規化するように構成されている。
【００３４】
処理手段３２０は、専用ハードウェア（たとえばＡＳＩＣ）、構成可能ハードウェア（たとえばＦＰＧＡ）、プログラマブルな構成要素（たとえば、ＤＳＰもしくは適切なソフトウェアを備えた汎用プロセッサ）、またはこれらのいずれかの組合せで実現され得る。同じ構成要素が他の機能を含んでもよい。システム３００はさらに、処理手段３２０によって処理されるように複数画像のセットを出力するように適合された出力手段３３０を含む。
【００３５】
「入力手段」および「出力手段」という語は、それぞれ、データを提供して受取ることのできる別のエンティティと通信するのに必要なハードウェアおよびソフトウェアを示している。好ましくは、このようなハードウェアおよびソフトウェアは、一般に容認されている業界基準に従って機能する。したがって、インターフェイス手段の物理およびデータリンク層の局面は、ＩＥＥＥ  Ｓｔｄ  ８０２．３（イーサネット（登録商標））、ＩＥＥＥ  Ｓｔｄ  ８０２．１１（無線ＬＡＮ）、ＵＳＢなどの規格に従って機能し得る。インターフェイス手段のネットワークおよびトランスポート層の局面はＴＣＰ／ＩＰプロトコルスタックに従って機能し得る。この明細書中において言及されたさまざまなインターフェイス（３１０および３２０）はハードウェアおよび／またはソフトウェアを共有し得る。
【００３６】
上述の正規化プロセスの後、本発明に従った方法は、正規化された画像に対して実行されるべきさらなるステップを含み得る。本発明に従ったシステムについて図示された実施形態においては、これらのステップは、処理手段３２０によって実行されてもよく、または出力手段３３０に接続された追加のシステムによって実行されてもよい。
【００３７】
図４は、本発明の一実施形態に従った、正規化後に有利に実行され得る例示的な写真測量動作を表わすフローチャートを示す。当該方法は、それぞれ別個の波長で撮像された放射を表わす別個の区域を含む一連の画像において撮像されたオブジェクトの写真測量３Ｄ復元を実行する。これらの画像は、最初に、典型的には上述のとおりハイパースペクトルセンサで取得される（４１０）。センサは航空機に搭載されるものであってもよい。正規化動作は、図２に関連付けて上述したステップに対応する前処理ステップ２００として示されている。当該方法は、次いで、上記一連の画像から複数のサブセットを選択するステップ（４２０）を含む。複数のサブセットの各々は複数の画像を含んでおり、その各々の画像は、同じサブセットにおける少なくとも１つの他の画像の視野と部分的に重複する視野を表わしている。好ましくは、複数のサブセットは互いにばらばらであり、上記サブセットを結合したものは上記一連の画像と一致している。次に、複数サブセットのそれぞれにおける複数画像に対して写真測量３Ｄ復元を実行することによって、中間３Ｄモデルのセットが生成される（４３０）。これらの中間３Ｄモデルは、次いで、当該３Ｄモデルのセットから複合３Ｄモデルへと再結合される（４４０）。マルチスペクトル画像またはハイパースペクトル画像の場合、３Ｄ復元は、たとえば、位置の移動と、さまざまな画像に現われるのと同じ特徴を有する関連する視野角差とを検出することによって、単一の波長帯域画像に対して同じ様に進めることができる。
【００３８】
上述の実施形態は写真測量３Ｄ復元に関するものであるが、本発明の正規化方法は、有利には、他の任意のアルゴリズムのための前処理ステップとして用いられてもよい。他の任意のアルゴリズムは、連続したさまざまなスペクトル画像で撮像された同じ物理的特徴を認識することに依存しているとともに、これらさまざまなスペクトル画像が異なる波長感度を持つ複数帯域に位置しているせいでこれらさまざまなスペクトル画像においてさまざまな異なる態様で現われる。
【００３９】
本発明はまた、コンピューティング装置上で実行されたときに本発明に従った方法のうちのいずれかの機能を提供するコンピュータプログラムプロダクトを含む。このようなコンピュータプログラムプロダクトは、プログラマブルプロセッサによって実行されるべき機械読取り可能コードを担持するキャリア媒体において有形的に具体化することができる。このため、本発明は、コンピューティング手段上で実行されたときに上述の方法のいずれかを実行するための命令を提供するコンピュータプログラムプロダクトを担持するキャリア媒体に関する。「キャリア媒体」という語は、実行すべき命令をプロセッサに提供することに関与する如何なる媒体をも指している。このような媒体は、不揮発性媒体および伝送媒体を含むがこれらに限定されない多くの形態を取り得る。不揮発性媒体は、たとえば、大容量記憶装置の一部である記憶装置などの光ディスクまたは磁気ディスクを含む。コンピュータ読取り可能媒体の一般的な形態は、ＣＤ−ＲＯＭ、ＤＶＤ、フレキシブルディスクもしくはフロッピー（登録商標）ディスク、テープ、メモリチップもしくはカートリッジ、または、コンピュータが読取り可能な他の如何なる媒体をも含む。コンピュータ読取り可能媒体のさまざまな形態は、プロセッサに対する実行すべき１つ以上の命令の１つ以上のシーケンスを担持するのに関与し得る。コンピュータプログラムプロダクトはまた、ＬＡＮ、ＷＡＮまたはインターネットなどのネットワークの搬送波を介して送信することができる。伝送媒体は、電波通信および赤外線データ通信中に生成されるような音波または光波の形態を取ることができる。伝送媒体は、コンピュータ内におけるバスを含むワイヤを含め、同軸ケーブル、銅線および光ファイバーを含む。
【００４０】
特定の実施形態に関連付けて本発明を上述してきたが、これは本発明を明確にするためになされたものであって、本発明を限定するためになされたものではない。当業者であれば、開示された特徴のさまざまな変更例およびさまざまな組合せが本発明の範囲から逸脱することなく実現可能であることを認識するだろう。
【図１】
2018521418.tif 000003
【図２】
2018521418.tif 000004
【図３】
2018521418.tif 000005
【図４】
2018521418.tif 000006
【手続補正書】
【提出日】平成30年3月9日(2018.3.9)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
それぞれ別個の波長で撮像された放射を表わす別個の複数領域を含む複数の二次元画像のセットを変換するための、コンピュータによって実現される方法であって、
−  前記セットにおける複数画像を等しく配置された複数区域に分割するステップ（２１０）と、
−  前記複数区域の各々に関して、前記複数画像のセット中の予め定められた特徴を計算するステップ（２２０）と、
−  前記複数画像の各々に関して、前記複数区域の各々における複数の輝度値を前記区域の前記予め定められた特徴の関数で正規化するステップ（２３０）とを含む、コンピュータによって実現される方法。
【請求項２】
前記予め定められた特徴は平均輝度であり、前記正規化するステップは、前記複数区域の各々における前記複数の輝度値を前記平均輝度値に対して正規化するステップを含む、請求項１に記載の、コンピュータによって実現される方法。
【請求項３】
前記複数区域は個々の画素に対応する、請求項１または２に記載の、コンピュータによって実現される方法。
【請求項４】
前記複数区域は、別個の波長帯域を表わすそれぞれの複数の画素を含む複数の矩形ブロックに対応する、請求項１または２に記載の、コンピュータによって実現される方法。
【請求項５】
前掲の請求項のいずれかに記載の方法のステップを実行するようプロセッサに指示するように構成されたコード手段を含むコンピュータプログラムプロダクト。
【請求項６】
それぞれ別個の波長で撮像された放射を表わす別個の複数領域を含む画像を変換するためのシステム（３００）であって、
−  それぞれ別個の波長で撮像された放射を表わす別個の複数領域を含む複数の二次元画像のセットを受取るように適合された入力手段（３１０）と、
−  処理手段（３２０）とを含み、前記処理手段（３２０）は、
−  前記セットにおける複数画像を等しく配置された複数区域に分割し、
−  前記複数区域の各々に関して、前記複数画像のセット中の予め定められた特徴を計算し、
−  前記複数画像の各々に関して、前記複数区域の各々における複数の輝度値を前記区域の前記予め定められた特徴の関数で正規化するように構成されており、前記システム（３００）はさらに、
−  前記処理手段（３２０）によって処理されるとおりに前記複数画像のセットを出力するように適合された出力手段（３３０）を含む、システム（３００）。
【請求項７】
前記予め定められた特徴は平均輝度であり、前記正規化することは、前記複数区域の各々における前記複数の輝度値を前記平均輝度値に対して正規化することを含む、請求項６に記載のシステム。
【請求項８】
前記複数区域は個々の画素に対応する、請求項７に記載のシステム。
【請求項９】
前記複数区域は、別個の波長帯域を表わすそれぞれの複数の画素を含む複数の矩形ブロックに対応する、請求項７に記載のシステム。
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図２
【補正方法】変更
【補正の内容】
【図２】
2018521418.tif 000007
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図４
【補正方法】変更
【補正の内容】
【図４】
2018521418.tif 000008
【国際調査報告】
2018521418.tif 000009
2018521418.tif 000010