(57)【特許請求の範囲】
【請求項１】
表示パラメータ値の算出に用いる係数パラメータ値を医用画像に関する撮像部位、撮像条件、画像関連情報、ユーザ識別情報のうちの２以上の項目に対応付けて格納し、
表示する医用画像に対応する係数パラメータ値を読み出し、
読み出された係数パラメータ値に基づいて表示パラメータ値を算出し設定する医用画像の表示パラメータ値設定方法であって、
ユーザによる前記表示パラメータ値の修正を受け付け、
前記表示パラメータ値がユーザにより修正される度にその修正回数を反映する収束関数で前記係数パラメータ値を更新することを特徴とする医用画像の表示パラメータ値設定方法。
【請求項２】
請求項１に記載の医用画像の表示パラメータ値設定方法において、
前記医用画像はＭＲ画像であり、
前記表示パラメータ値はＷＷ値またはＷＬ値であることを特徴とする医用画像の表示パラメータ値設定方法。
【請求項３】
請求項１または請求項２に記載の医用画像の表示パラメータ値設定方法において、
係数パラメータ値の修正回数をＮｔとし、更新前の係数パラメータ値をＰとし、ユーザによる修正後の表示パラメータ値をＷとし、医用画像から得られる統計量をＳとするとき、更新後の係数パラメータ値Ｐ’を、
Ｐ’＝（Ｎｔ×Ｐ＋Ｗ／Ｓ）／（Ｎｔ＋１）
なる収束関数で算出することを特徴とする医用画像の表示パラメータ値設定方法。
【請求項４】
ＷＷ値またはＷＬ値の算出に用いる係数パラメータ値をＭＲ画像に関する撮像部位、撮像条件、画像関連情報、ユーザ識別情報のうちの２以上の項目に対応付けて格納する係数パラメータ値格納手段と、
表示するＭＲ画像に対応する係数パラメータ値を前記係数パラメータ値格納手段から読み出してその係数パラメータ値に基づいてＷＷ値またはＷＬ値を算出するＷＷ／ＷＬ値算出手段と、
算出されたＷＷ値またはＷＬ値により画像上に前記ＭＲ画像を表示するＭＲ画像表示手段と、
ユーザによる前記ＷＷ値またはＷＬ値の修正を受け付けるＷＷ／ＷＬ値修正受付手段と、
前記ＷＷ値またはＷＬ値がユーザにより修正される度にその修正回数を反映する収束関数で前記ＷＷ／ＷＬ値を更新するＷＷ／ＷＬ更新手段とを具備したことを特徴とするＭＲＩ装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、医用画像の表示パラメータ値設定方法、医用画像診断装置およびＭＲＩ装置に関し、さらに詳しくは、医用画像に適合した表示パラメータ値を容易に設定できる医用画像の表示パラメータ値設定方法、医用画像診断装置およびＭＲＩ装置に関する。特に、ＭＲ（Magnetic  Resonance）画像のＷＷ（Window Width）値またはＷＬ（Window Level）値を設定するのに有用である。
【０００２】
【従来の技術】
一般に、ＭＲ画像やＣＴ画像（Computed Tomography）などの医用画像の表示に際して、ＷＷ値やＷＬ値などの表示パラメータを設定する必要がある。ＷＷ値やＷＬ値の設定によって画像の見え方が全く異なるため、医用画像を表示する度に、ユーザ（医師や技師）は、適正なＷＷ値やＷＬ値を設定する。
しかし、画像ごとに設定を行うのは作業負担が大きいので、設定を支援して作業負担を軽減した次の従来技術が知られている。
（１）特開平９−２９９３４９号公報には、図７，図８に示すように、トレーニングシリーズと呼ばれる一連の断層画像上で適正に設定されたＷＷ／ＷＬ値をテーブル曲線ＷＷ（Ｘ’），ＷＬ（Ｘ’）として持ち、設定対象シリーズと呼ばれる一連の断層画像の画素値の統計量（標準偏差／平均値）を基にして算出されたＷＷオフセット量／ＷＬオフセット量で前記テーブル曲線ＷＷ（Ｘ’），ＷＬ（Ｘ’）を修正して出力テーブル曲線ＷＷA（Ｘ’），ＷＬA（Ｘ’）を作成し、その出力テーブル曲線ＷＷA（Ｘ’），ＷＬA（Ｘ’）を用いて、前記設定対象シリーズの各断層画像のスライス位置Ｘ’（頭下部，頭中部，頭上部）対応のＷＷ値／ＷＬ値を出力するようにした技術が開示されている。
（２）図９に示すように、係数パラメータテーブル５０に、撮像部位ごとに係数パラメータＰww／ＰwLを予め格納しておき、表示する画像の撮像部位が指定されると、前記係数パラメータテーブル５０から該当する係数パラメータＰww／ＰwLを読み出し、それらを用いて、
ＷＷ＝Ｐww×Ｓ１
ＷＬ＝ＰwL×Ｓ２
（Ｓ１，Ｓ２はＭＲ画像から得られる統計量）
によりＷＷ／ＷＬ値を算出し、ＷＷ値／ＷＬ値を自動設定する技術が知られている。なお、前記Ｓ１は、例えばＭＲ画像の画素値の標準偏差である。前記Ｓ２は、例えばＭＲ画像の画素値の平均値である。
【０００３】
【発明が解決しようとする課題】
上記（１）の従来技術では、設定対象シリーズの断層画像とトレーニングシリーズの断層画像の種類が基本的に一致している必要があり、両者の種類が異なる場合（例えば一方の断層画像の撮像角がアキシャルで，他方の断層撮像角がサジタルのとき）には、適用し難い。
上記（２）の従来技術では、撮像部位に係数パラメータＰww，ＰwLを対応付けるだけなので、ＭＲ画像としての適合性に影響を与える他の要因（例えば撮像条件）を十分に反映することが出来ない。
そこで、本発明の目的は、医用画像に適合した表示パラメータ値を容易に設定できる医用画像の表示パラメータ値設定方法、医用画像診断装置およびＭＲＩ装置を提供することにある。
【０００４】
【課題を解決するための手段】
第１の観点では、本発明は、表示パラメータ値の算出に用いる係数パラメータ値を医用画像に関する撮像部位，撮像条件，画像関連情報，ユーザ識別情報のうちの２以上の項目に対応付けて格納し、表示する医用画像に対応する係数パラメータ値を読み出し、読み出された係数パラメータ値に基づいて表示パラメータ値を算出し設定することを特徴とする医用画像の表示パラメータ値設定方法を提供する。
上記構成において、医用画像の撮像部位としては、例えば、頭部（Head），脊椎（Spine）等の種別がある。撮像条件としては、例えば撮像法の種別や，撮像角の種別がある。ＭＲ画像の場合、撮像法としては、ＳＥ（Spin Echo）法や，ＧＲＥ（GRadient Echo）法や，ＩＲ（Inversion Recovery）法などがある。撮像角としては、アキシャル（Axial），サジタル（Sagital），コロナル（Coronal）などがある。画像関連情報としては、医用画像のＳＮＲ（Signal to Noise Ratio）がある。ユーザ識別情報としては、例えば医師や技師の識別情報がある。上記第１の観点による表示パラメータ値設定方法では、係数パラメータ値を医用画像に関する撮像部位，撮像条件，画像関連情報，ユーザ識別情報のうちの２以上の項目に対応付けて格納するので、撮像条件や画像関連情報により微妙に異なる画像特性やユーザの好みを忠実に反映した係数パラメータ値を容易に得ることが出来る。したがって、適正値に近い表示パラメータ値を容易に設定することが可能となり、ユーザの作業負担を軽減することが出来る。
【０００５】
第２の観点では、本発明は、上記第１の観点の医用画像の表示パラメータ値設定方法において、前記医用画像はＭＲ画像であり、前記表示パラメータ値はＷＷ値またはＷＬ値であることを特徴とする医用画像の表示パラメータ値設定方法を提供する。
上記第２の観点による表示パラメータ値設定方法では、ＭＲ画像のＷＷ値やＷＬ値を設定する。したがって、適正値に近いＷＷ値やＷＬ値を容易に設定することが可能となり、ユーザの作業負担を軽減することが出来る。
【０００６】
第３の観点では、本発明は、上記第１の観点または上記第２の観点の医用画像の表示パラメータ値設定方法において、ユーザによる前記表示パラメータ値の修正を受け付け、前記表示パラメータ値がユーザにより修正される度にその修正回数を反映する収束関数で前記係数パラメータ値を更新することを特徴とする医用画像の表示パラメータ値設定方法を提供する。
上記第３の観点による表示パラメータ値設定方法では、画面上に表示された医用画像を見ながら、ユーザが表示パラメータ値を容易に修正できる。
また、表示パラメータ値がユーザにより修正される度にその修正回数を反映する収束関数で係数パラメータ値を更新する。すなわち、修正回数が少ない時点ほど、大きな重みで係数パラメータ値を更新する。そして、修正回数が多くなると、小さな重みで係数パラメータ値を更新する。したがって、修正を繰り返すのにつれて、係数パラメータ値が収束する。
この結果、臨床上の目的やユーザごとの好みに十分に合致した医用画像を得ることが出来る。
【０００７】
第４の観点では、本発明は、上記第３の観点の表示パラメータ値設定方法において、係数パラメータ値の修正回数をＮｔとし、更新前の係数パラメータ値をＰとし、ユーザによる修正後の表示パラメータ値をＷとし、医用画像から得られる統計量をＳとするとき、更新後の係数パラメータ値Ｐ’を、Ｐ’＝（Ｎｔ×Ｐ＋Ｗ／Ｓ）／（Ｎｔ＋１）なる収束関数で算出することを特徴とする医用画像の表示パラメータ値設定方法を提供する。
本願の発明者の研究によれば、上式により更新後の係数パラメータ値Ｐ’を算出することで、修正回数Ｎｔに対して適切な重みで更新値を算出できることが判った。また、設定対象の医用画像から得られる統計量Ｓ（医用画像の画素値の標準偏差や平均値など）を反映することで、更新値をいっそう適正に算出できることが判った。
そこで、上記第４の観点による表示パラメータ値設定方法では、表示パラメータ値をより容易に適正値に近づけることが可能となり、ユーザ好みの医用画像を容易に得ることが出来る。
【０００８】
第５の観点では、本発明は、表示パラメータ値の算出に用いる係数パラメータ値を医用画像に関する撮像部位，撮像条件，画像関連情報，ユーザ識別情報のうちの２以上の項目に対応付けて格納する係数パラメータ値格納手段と、表示する医用画像に対応する係数パラメータ値を前記係数パラメータ値格納手段から読み出してその係数パラメータ値に基づいて表示パラメータ値を算出し設定する表示パラメータ値算出設定手段とを具備したことを特徴とする医用画像診断装置を提供する。
上記第５の観点による医用画像診断装置は、上記第１の観点による表示パラメータ値設定方法を好適に実施することが出来るため、適正値に近い表示パラメータ値を容易に設定することが可能となり、ユーザの作業負担を軽減することが出来る。
【０００９】
第６の観点では、本発明は、ＷＷ値またはＷＬ値の算出に用いる係数パラメータ値をＭＲ画像に関する撮像部位，撮像条件，画像関連情報，ユーザ識別情報のうちの２以上の項目に対応付けて格納する係数パラメータ値格納手段と、表示するＭＲ画像に対応する係数パラメータ値を前記係数パラメータ値格納手段から読み出してその係数パラメータ値に基づいてＷＷ値またはＷＬ値を算出するＷＷ／ＷＬ値算出手段と、算出されたＷＷ値またはＷＬ値により画面上に前記ＭＲ画像を表示するＭＲ画像表示手段と、ユーザによる前記ＷＷ値またはＷＬ値の修正を受け付けるＷＷ／ＷＬ値修正受付手段と、前記ＷＷ値またはＷＬ値がユーザにより修正される度にその修正回数を反映する収束関数で前記ＷＷ／ＷＬ値を更新するＷＷ／ＷＬ値更新手段とを具備したことを特徴とするＭＲＩ装置を提供する。上記第６の観点によるＭＲＩ装置では、ＷＷ値またはＷＬ値がユーザにより修正される度にその修正回数を反映する収束関数で係数パラメータ値を更新する。すなわち、修正回数が少ない時点ほど、大きな重みで係数パラメータ値を更新する。そして、修正回数が多くなると、小さな重みで係数パラメータ値を更新する。したがって、修正を繰り返すのにつれて、係数パラメータ値が収束する。
この結果、臨床上の目的やユーザの好みにいっそう合致したＭＲ画像を容易に得ることが出来る。
【００１０】
【発明の実施の形態】
以下、図に示す実施形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
図１は、本発明の一実施形態にかかるＭＲＩ装置を示すブロック図である。
このＭＲＩ装置１００は、撮像装置１と、処理装置２と、記憶装置３と、表示装置４と、入力装置５とを具備している。
前記撮像装置１は、被検体Ｋをスキャンし、得られたスキャンデータを処理装置２に渡す。
前記処理装置２は、前記撮像装置１からスキャンデータを取得するデータ取得部２ａと、スキャンデータからＭＲ画像の画像データを生成する画像データ生成部２ｂと、ＷＷ値，ＷＬ値を設定するためのＷＷ／ＷＬ値設定部２ｃとを含んでいる。
前記記憶装置３は、前記スキャンデータや前記画像データを記憶する。
前記表示装置４は、前記ＭＲ画像や各種の情報を表示する。
前記入力装置５は、操作者が入力する情報を受け付ける。
【００１１】
図２は、前記ＷＷ／ＷＬ値設定部２ｃによりＷＷ値およびＷＬ値を設定する処理を示すフロー図である。
ステップＳＴ１では、修正回数Ｎｔ＝０に初期化する。
ステップＳＴ２では、記憶装置３から、一連のＭＲ画像（これをトレーニングシリーズのＭＲ画像という）のうち、一枚目のＭＲ画像を指定する。一連のＭＲ画像は、例えば、被検体の脊椎に沿って一定間隔で撮像した一群のマルチスライスイメージである。
【００１２】
ステップＳＴ３では、指定されたＭＲ画像を読み込む。
ステップＳＴ４では、ＭＲ画像に関するヘッダ情報を読み取る。一般に、ヘッダ情報は、撮像部位（例えば頭部，脊椎）や、撮像法（例えばＳＥ法，ＧＲＥ法，ＩＲ法）や、コントラストの種類（例えば造影剤でコントラスト強調したＴ１強調画像）や、撮像角（例えばアキシャル，サジタル，コロナル）や、スキャンパラメータ（繰り返し時間ＴＲ，エコー時間ＴＥ，反転時間ＴＩ）を示すフラグを含む。
ステップＳＴ５では、上記ヘッダ情報に基づいて、ＭＲ画像の種類を判別する。ＭＲ画像の種類は、例えば、撮像部位，撮像法，コントラストの種類，撮像角により記述される。具体例を示せば、Head／SE／T1WI／Ax（頭部，ＳＥ法，Ｔ１強調画像，アキシャル）や，Spine／SE／T2WI／Sag（脊椎，ＳＥ法，Ｔ２強調画像，サジタル）である。
【００１３】
ステップＳＴ６では、次式により、ＭＲ画像のＳＮＲを算出し、整数化する。
ＳＮＲ＝(ＦＯＶx／Ｍx)×(ＦＯＶy／Ｍy)×Ｔhk×√｛ＮEX／ＢW｝×100
ただし、ＦＯＶx，ＦＯＶyは、ｘ方向，ｙ方向の撮像領域長〔ｍｍ〕である。Ｍx，Ｍyは、撮像マトリクス数である。Ｔhkは、スライス厚〔ｍｍ〕である。ＮEXは、励起回数である。ＢWは、受信バンド幅〔Ｈｚ〕である。例えば、ＦＯＶx×ＦＯＶy＝２２０×２２０〔ｍｍ〕，Ｍx×Ｍy＝２５６×２５６，Ｔhk＝５〔ｍｍ〕，ＮEX＝２，ＢW＝１６０００〔Ｈｚ〕の場合、ＳＮＲ≒４．１となるので、四捨五入によりＳＮＲ＝４に整数化される。
【００１４】
ステップＳＴ７では、ＭＲ画像の種類に対応した係数パラメータテーブルから係数パラメータＰww，ＰwLを読み出す。
図３に、ＭＲ画像の種類がHead／Axに対応する係数パラメータテーブル１０ａ’を例示する。図４に、ＭＲ画像の種類がHead／Ax以外に対応する係数パラメータテーブル１０ｂ’を例示する。
前記係数パラメータテーブル１０ａ’，１０ｂ’内には、ユーザの好みを入れていない標準的な係数パラメータＰww，ＰwLが予め格納されている。また、初期状態では、ユーザによるＷＷ値，ＷＬ値の修正（後述のステップＳＴ１１）を１回も受け付けていないから、修正回数Ｎｔとして、すべて“０”が格納されている。
例えば、Head／SE／T1WI／Ax，ＳＮＲ＝１１，スライス位置（特開平９−２９９３４９号公報に開示のように、頭部像面積および関心領域内平均画素値を基にして、スライス位置を決定する）が脳下部ＬのＭＲ画像の場合、図３の係数パラメータテーブル１０ａ’からＳＮＲ＝１１，脳内のスライス位置＝脳下部Ｌに対応する係数パラメータＰww＝４．７９，ＰwL＝１．１４を読み出す。また、例えば、Spine／SE／T2WI／Sag，ＳＮＲ＝８のＭＲ画像の場合、図４の係数パラメータテーブル１０ｂ’からＳＮＲ＝８に対応する係数パラメータＰww＝６．１９，ＰwL＝１．４６を読み出す。
【００１５】
ステップＳＴ８では、読み出した係数パラメータＰww，ＰwLを用いて、次式により、ＷＷ値，ＷＬ値を算出する。
ＷＷ＝Ｐww×ＳＤ
ＷＬ＝ＰwL×Mean
ただし、ＳＤは、ＭＲ画像の画素値の標準偏差である。Meanは、ＭＲ画像の画素値の平均値である。なお、前記標準偏差ＳＤおよび前記平均値Meanは、ＭＲ画像の関心領域について算出することが好ましい。関心領域の自動算出方法として、特開平９−２９９３４９号公報には、ＭＲ画像の像輪郭を円で近似し、その円の半径を所定値（例えば２０ｍｍ）だけ縮めた円の後ろ半分を関心領域として算出する方法が開示されている。
【００１６】
ステップＳＴ９では、上記ステップＳＴ６で算出されたＷＷ値，ＷＬ値により、ＭＲ画像を表示装置４の画面上に表示する。
ステップＳＴ１０では、ユーザ（技師や医師など）は、表示されたＭＲ画像を見て、ＷＷ値，ＷＬ値が適正か否か判定する。ＷＷ値と，ＷＬ値の両方が適正であればステップＳＴ１４に進み、ＷＷ値と，ＷＬ値の少なくとも一方が適正でなければステップＳＴ１１に進む。
ステップＳＴ１１では、ユーザは、入力装置５により、ＷＷ値，ＷＬ値を適正値に修正する。修正後のＷＷ値をＷＷｍとし、修正後のＷＬ値をＷＬｍとする。
【００１７】
ステップＳＴ１２では、修正回数Ｎｔを１だけインクリメントする。
ステップＳＴ１３では、次式により、係数パラメータＰww，ＰwLを更新し、更新後の係数パラメータＰ’ww，Ｐ’wLおよび修正回数Ｎｔを、係数パラメータテーブル１０ａ（図３）または１０ｂ（図４）の該当箇所に格納する。
Ｐ’ww＝（Ｎｔ×Ｐww＋ＷＷｍ／ＳＤ）／（Ｎｔ＋１）
Ｐ’wL＝（Ｎｔ×ＰwL＋ＷＬｍ／Mean）／（Ｎｔ＋１）
【００１８】
ステップＳＴ１４では、トレーニングシリーズ内のＭＲ画像をすべて表示したか否か判定し、すべて表示したらＷＷ／ＷＬ値設定処理を終了し、そうでなければステップＳＴ１５に進む。
ステップＳＴ１５では、トレーニングシリーズ内で、次のＭＲ画像を指定する。その後、上記ステップＳＴ３に戻る。
【００１９】
図５，図６は、上記ＷＷ／ＷＬ値設定処理により、係数パラメータＰww，ＰwLおよび学習回数Ｎｔを更新後の係数パラメータテーブル１０ａ，１０ｂを示す。図５の係数パラメータテーブル１０ａでは、例えばＳＮＲ＝１１，スライス位置＝脳下部Ｌに対応して、Ｐww＝４．６０，ＰwL＝１．１４，Ｎｔ＝２４が格納されている。図６の係数パラメータテーブル１０ｂでは、例えばＳＮＲ＝８に対応して、係数パラメータＰww＝５．１８，ＰwL＝１．３２，Ｎｔ＝６が格納されている。
【００２０】
以上のＭＲＩ装置１００によれば、ＷＷ値またはＷＬ値をユーザが修正する度に、その修正回数Ｎｔを反映した重みで係数パラメータテーブル１０ａ，１０ｂ内の係数パラメータ値Ｐww，ＰwLを更新する。したがって、前記更新後には、ＭＲ画像に対応した係数パラメータ値Ｐww，ＰwLを係数パラメータテーブル１０ａ，１０ｂから読み出してＷＷ値，ＷＬ値を算出すれば、ユーザ好みのＭＲ画像を表示し得るＷＷ値，ＷＬ値を自動設定することが出来る。
【００２１】
なお、以上のＭＲＩ装置１００の構成を、次のように変更してもよい。
（１）上記ＷＷ／ＷＬ値設定処理（図２参照）により更新された係数パラメータ値Ｐww，ＰwLをそのまま用いて複数のＭＲ画像を切り換えて表示すると、撮像位置が隣接する画像間で視覚的に不自然な印象を与える（ＷＷ値，ＷＬ値に過大な差異が生じる）場合がある。これを軽減するため、表示対象のＭＲ画像と，その近傍のＭＲ画像のＷＷ値，ＷＬ値の平均値を算出し、それを設定値として採用してもよい。
（２）係数パラメータテーブル１０ａ，１０ｂ（図５，図６参照）を、ユーザ識別子に対応して作成し、ユーザごとに係数パラメータＰww，ＰwLを管理してもよい。この場合には、ＭＲ画像の好みがユーザごとに大きく異なる場合でも、各ユーザに適したＷＷ／ＷＬ値を設定することが出来る。
【００２２】
【発明の効果】
本発明の医用画像の表示パラメータ値設定方法、医用画像診断およびＭＲＩ装置によれば、装置係数パラメータ値を医用画像に関する撮像部位，撮像条件，画像関連情報，ユーザ識別情報のうちの２以上の項目に対応付けて格納するので、係数パラメータ値をきめ細かく管理することができ、適正な表示パラメータ値を容易に設定することが出来る。また、表示パラメータ値がユーザにより修正される度にその修正回数を反映する収束関数で係数パラメータ値を更新することで、係数パラメータ値を迅速に適正値に収束させることができ、設定効率を向上することが出来る。
特に、ＭＲ画像のＷＷ値またはＷＬ値を設定すれば、ユーザ好みのＭＲ画像を容易に表示可能となり、臨床上の有用性が高い。
【図面の簡単な説明】
【図１】本発明の一実施形態にかかるＭＲＩ装置を示すブロック図である。
【図２】図１のＭＲＩ装置によるＷＷ／ＷＬ値設定処理を示すフロー図である。
【図３】ＭＲ画像の種類がHead／Axに対応する係数パラメータテーブルの初期状態を示す説明図である。
【図４】ＭＲ画像の種類がHead／Ax以外に対応する係数パラメータテーブルの初期状態を示す説明図である。
【図５】ＭＲ画像の種類がHead／Axialに対応する更新後の係数パラメータテーブルを示す説明図である。
【図６】ＭＲ画像の種類がHead／Axial以外に対応する更新後の係数パラメータテーブルを示す説明図である。
【図７】出力ＷＷテーブル曲線の説明図である。
【図８】出力ＷＬテーブル曲線の説明図である。
【図９】従来のＷＷ／ＷＬ値設定方法に用いる係数パラメータテーブルを示す説明図である。
【符号の説明】
１００              ＭＲＩ装置
２                  処理装置
２ａ                データ取得部
２ｂ                画像データ生成部
２ｃ                ＷＷ／ＷＬ値設定部
４                  表示装置
５                  入力装置
１０ａ，１０ｂ      係数パラメータテーブル（初期状態）
１０ａ’，１０ｂ’  係数パラメータテーブル（更新後）
【図１】
0003746158.tif 000002
【図２】
0003746158.tif 000003
【図３】
0003746158.tif 000004
【図４】
0003746158.tif 000005
【図５】
0003746158.tif 000006
【図６】
0003746158.tif 000007
【図７】
0003746158.tif 000008
【図８】
0003746158.tif 000009
【図９】
0003746158.tif 000010