(57)【要約】  本発明は、イメージングプレートの品質を確定するための方法であって、イメージングプレートの露光を実施するステップと、画像を確定するためにイメージングプレートのスキャンを実施するステップと、画像の信号対雑音比を確定するステップまたは／および画像におけるエッジ認識を実施するステップと、画像の信号対雑音比に基づいてまたは／および認識されたエッジ構造に基づいてイメージングプレートの品質値を計算するステップとを有する方法に関する。さらに本発明は、そのような方法を実施するためのイメージングプレートスキャナに関する。
2020537157.tif 000002
【特許請求の範囲】
【請求項１】
イメージングプレートの品質を確定するための方法であって、
ａ）前記イメージングプレートの露光を実施するステップ（Ｓ１）と；
ｂ）画像を確定するために前記イメージングプレートのスキャンを実施するステップ（Ｓ２）と；
ｃ）前記画像の信号対雑音比を確定するステップ（Ｓ４）または／および前記画像におけるエッジ認識を実施するステップ（Ｓ５）と；
ｄ）前記画像の前記信号対雑音比に基づいてまたは／および前記認識されたエッジ構造に基づいて前記イメージングプレートの品質値を計算するステップ（Ｓ６、Ｓ７）と
を有する方法。
【請求項２】
前記露光を実施するステップが、前記プレートと記録デバイスの間の特定の間隔の調整を含む、請求項１に記載の方法。
【請求項３】
前記信号対雑音比を確定するステップが、マスクフィルタの適用を含む（Ｓ３１）、請求項１または２に記載の方法。
【請求項４】
前記信号対雑音比を確定するステップが、前記画像の１つまたは複数の部分に対する局所的な信号対雑音比の確定を含む（Ｓ４１、Ｓ４２）、請求項１〜３のいずれか１項に記載の方法。
【請求項５】
局所的に低い信号対雑音比に基づいて、損傷の数または／およびサイズを決定するステップを有する、請求項１〜４のいずれか１項に記載の方法。
【請求項６】
損傷が、認識されたエッジ構造に基づいて決定される（Ｓ５）、請求項１〜５のいずれか１項に記載の方法。
【請求項７】
前記エッジ認識を実施するステップが、キャニーアルゴリズムの使用を含む、請求項１〜６のいずれか１項に記載の方法。
【請求項８】
前記エッジ認識を実施するステップの後に、モルフォロジー膨張（形態的拡張）または／およびクロージング演算が実施される（Ｓ５４）、請求項１〜７のいずれか１項に記載の方法。
【請求項９】
イメージングプレートの前記品質を決定するため、損傷の数が、低い信号対雑音比に基づく損傷を明らかにする画像と、エッジ認識に基づく損傷を明らかにする画像との間の論理ＯＲ結合を明らかにする画像に基づいて確定される、請求項１〜８のいずれか１項に記載の方法。
【請求項１０】
前記方法の１つを実施するために適応されたイメージングプレートスキャナ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、イメージングプレートの品質を確定するための方法およびそのためのイメージングプレートに関する。
【背景技術】
【０００２】
Ｘ線技術において、とりわけ歯科医療のＸ線技術において、今日では、Ｘ線画像を撮影するためのイメージングプレートが使用される。このイメージングプレートは、透明なマトリクス中に包埋された蛍光体材料を含んでいる。これによりいわゆる蓄積中心が生じ、この蓄積中心は、差し込んでくるＸ線光により、励起された準安定状態になり得る。このようなイメージングプレートをＸ線装置内で、例えば患者の咬翼撮影のために露光すると、イメージングプレートはＸ線潜像を、励起されたおよび励起されていない蓄積中心の形態で内包する。
【０００３】
イメージングプレートを読み取るには、イメージングプレートをスキャン装置内で読み取り光によって逐次走査し、これにより、励起された蓄積中心の準安定状態が、蛍光を発しながら高速緩和する状態になる。この蛍光は検出ユニットによって捕捉でき、これにより相応の評価電子機器で、Ｘ線画像を見えるようにすることができる。読み取り工程には、例えば円筒面に沿ったイメージングプレートを読み取り間隙に通すドラムスキャナが用いられる。
【０００４】
イメージングプレート技術の大きな利点は、イメージングプレートの再利用可能性にある。つまりイメージングプレートは、プレート内で蓄積された画像情報がどのみち消去される読み取り後に、一般的には多数回のさらなる撮影工程および蓄積工程に用いることができる。再利用可能性を制限する様々な老朽化プロセスがある。このような老朽化プロセスは、実質的には、時間が経つにつれて形成されるもう励起させられない、したがって蓄積像において暗いままの局所的な蓄積中心を形成する。実際の運用では機械的負荷も生じる可能性があり、したがって、例えば不適切な取扱いにより、イメージングプレートの表面に引っかき傷または点状の損傷が生じ得る。
【０００５】
実際の運用では、イメージングプレートのその時々の蓄積品質を知ることが望ましい。これに関し従来的には、イメージングプレートを目視して、表面にある引っかき傷の、撮影工程および読み取り工程への影響について判断することができる。さらにイメージングプレートを使った撮影の回数を数えることができ、または製造日に基づいてイメージングプレートの寿命を決定することができる。しかしこれらすべての措置は、イメージングプレートの実際の蓄積能力についての信頼できるデータを示してはいない。
【発明の概要】
【０００６】
本発明の課題は、上述の欠点を回避し、とりわけイメージングプレートの蓄積能力および蓄積品質についての正確なデータを示すことを可能にする、イメージングプレートの品質を確定するための方法を提示することである。
【０００７】
この課題は、独立請求項に基づく方法によって解決される。本発明による方法は以下のステップを有する。
【０００８】
イメージングプレートの露光の実施、画像を確定するためのイメージングプレートのスキャンの実施、画像の信号対雑音比の確定または／および画像におけるエッジ認識の実施、ならびに画像の信号対雑音比に基づくまたは／および認識されたエッジ構造に基づくイメージングプレートの品質値の計算。
【０００９】
本発明による方法は、代替的にまたは組み合わせて用い得る２つの品質特徴を提供する。一方の特徴は、言及した信号対雑音比である。画像の信号対雑音比は、イメージングプレートの露光の実施後に、イメージングプレートが、機能を果たす能力のない蓄積中心の意味における欠陥を有するかどうかの示唆をもたらす。読み取り工程中のイメージングプレートの一様な露光にもかかわらず、蓄積中心が蛍光を発しなければ、信号対雑音比が悪化している。
【００１０】
もう一方の特徴は、エッジ認識の実施である。画像におけるエッジ認識も、欠陥の可能性への示唆をもたらす。イメージングプレートを一様に露光した場合、生じる画像は一様の照度、つまり一様のグレー値になることが望ましい。これに対し、エッジ認識により構造が見つかる場合、これはイメージングプレートの欠陥の可能性を示唆している。
【００１１】
イメージングプレートの露光は、特定の線量で実施することが好ましい。これは、信号対雑音比の決定またはエッジ認識の実施を容易にする。ただし未知の線量でも、露光およびその後の計算を実施することはできる。ただし使用する線量は、イメージングプレートが露光不足にならず、かつ露光過剰にならない線量範囲内であることが望ましい。線量は、例えば曝露時間、Ｘ線電圧、および陽極電圧を通して調整できる。
【００１２】
本発明の有利な一形態では、露光を実施するステップが、プレートと記録デバイスの間の特定の間隔の調整を含んでいる。これは、画像の予測される絶対的グレー値の決定を、したがってイメージングプレートの品質値の計算の再現性を容易にし、しかしこの方法のスムーズな実施に必要な前提条件ではない。イメージングプレートが露光過剰にも露光不足にもならないことが保証されればよい。原則的には、線量および間隔からなる値のペアを相応に選択できる。例えばより大きな間隔はより大きな線量によって相殺できる。
【００１３】
本発明の一実施形態では、イメージングプレートの露光が、通常の診察の枠内で行われる。これは、欠陥認識の方法が、例えば患者のところで露光された、したがって画像を含むイメージングプレートにおいても実施できることを意味する。評価方法、つまり信号対雑音比の確定または／およびエッジ認識の実施は、認識された欠陥を記録して、この欠陥をそのイメージングプレートに割り当てる。このイメージングプレートの使用を繰り返し、改めて欠陥認識を実施した後に再び発生している欠陥を、そのようなものとして登録して特徴づけることができる。発生し得るアーチファクト、例えばＸ線を通さず、撮影で反映される物体（例えば金属の詰め物）は、相応に認識し、評価し、かつそのようなものとして特徴づけることができる。
【００１４】
本発明の一発展形態は、信号対雑音比の確定が、マスクフィルタの適用を含んでいる。このマスクフィルタは、例えばイメージングプレートの周縁、例えば１または２ｍｍを、信号対雑音比の計算から排除できる。これに関しイメージングプレートの周縁は、固定的に設定でき、または操作をする者によって入力可能であり得る。その代わりにまたはそれに加えて、マスクフィルタは周縁だけでなくイメージングプレートのそのほかのマーキングも排除できる。これに関しマーキングは、モルフォロジー（形態的）演算によっても認識でき、かつマスクフィルタに付加できる。
【００１５】
本発明の特に好ましい一実施形態では、信号対雑音比の確定が、画像の１つまたは複数の部分に対する局所的な信号対雑音比の確定を含んでいる。例えば、２１×２１ピクセルのサイズのエリアを、信号対雑音比の計算のための基礎として用いることができる。このような１つのエリアサイズに対して１つの局所的な信号対雑音比を、画像全体で計算できる。これにより、画像の各ピクセルに対して１つの帰属の信号対雑音比の値が生じる（２１×２１ピクセルに相応する周縁帯を除いて）。局所的な信号対雑音比が低い領域は、欠陥を有すると思われる領域に相当する。これに関し、例えば閾値を定めることができる。閾値は、例えば局所的に確定された信号対雑音比の最大値および最小値によって確定できる。このために例えば、最小値分を減らした信号対雑音の値と、最大値と最小値の差との比を確定できる。たった今説明した比が０．２以下である、つまり局所的な信号対雑音比が、有り得る値域の２０％以下である場合に、ピクセルを例えば欠陥としてマーキングできる。
【００１６】
その代わりにまたはそれに加えて一実施形態では、イメージングプレートの損傷を、認識されたエッジ構造に基づいて決定できる。エッジ認識により、欠損、つまり例えばもう機能を果たす能力のない蓄積中心を、そうでなければ一様に露光されるプレート内の構造として認識できる。
【００１７】
好ましい一実施形態では、エッジ認識の実施がキャニー（Ｃａｎｎｙ）アルゴリズムの使用を含んでいる。エッジ認識の実施前に、ぼかしまたはボケ（不鮮明さ）付加を実施できることが好ましい。これは、雑音割合の高いピクセルに基づくアーチファクトの生成を軽減する。さらにエッジ認識の実施後に、好ましくは二値画像を得るために、閾値演算を実施できる。
【００１８】
本発明の一変形形態では、エッジ認識の実施後に、とりわけ閾値演算の実施後に、モルフォロジー膨張（形態的拡張）または／およびクロージング演算を実施できる。モルフォロジー膨張およびクロージング（閉じる）演算により原則的に高いグレー値をもつ値は大きくなり、かつ比較的低いグレー値をもつ値は小さくなるので、不規則さ／損傷を有する単一領域の可視性を改善できる。同時に、そのような領域間の小さな隙間を塞ぐことができる。
【００１９】
本発明の好ましい一実施形態は、信号対雑音の演算から得られた画像と、エッジ認識から生じる画像とを、論理ＯＲ演算によって合わせる。このようにして判明するマーキングされたピクセルの総数に基づき、イメージングプレートの品質を決定できる。これに関しては例えば、このようにしてマーキングされたピクセルの数Ｎマーキングされたを、下式に基づいてイメージングプレートの品質を表す数値に換算できる。
【００２０】
【数１】
2020537157.tif 000003
【００２１】
式中、ｆマーキングされたはイメージングプレートの品質評価であり、Ｎマーキングされたはマーキングされた画像ピクセルの数を表し、Ｎ総数は考察される画像内のピクセルの総数を表し、Ｎマスクは除外マスクによってマスキングされるプレート周縁ピクセルおよびマーキングピクセルの数を表す。したがって、マーキングされたピクセルの数が非常に少ない画像を有するイメージングプレートは高い評価を獲得し、その一方で、その画像が多数のマーキングされた、つまり欠陥のあるピクセルを有するイメージングプレートは低い評価を獲得する。
【００２２】
イメージングプレートのこうして得られた評価を、０〜１の間の範囲内にスケーリングするため、および個々の確定された／マーキングされた特徴の重み付けの適合を可能にするため、評価ｆマーキングされたを以下のように変換できる。
【数２】
2020537157.tif 000004
【００２３】
これは、こうして得られた評価と、そのほかの評価基準、例えば視覚的鑑定およびその後のイメージングプレートの評価によるような評価基準との比較を可能にする。これに関しては例えば、パラメータｗ0およびｗ1は、トレーニングデータセットでの線形回帰により、生じる最終値がプレート品質の専門家の評価と可能な限り一致するように選択できる。
【００２４】
以下に、本発明の例示的実施形態を図面に基づいてより詳しく解説する。
【図面の簡単な説明】
【００２５】
【図１】本発明による方法の第１の実施形態のフロー図である。
【図２】図１の方法の第１の変形形態のフロー図である。
【図３】図１の方法の第２の実施形態のフロー図である。
【図４】図３の方法の第２の代替的または追加的な実施形態のフロー図である。
【図５】本発明によるイメージングプレートスキャナの概略図である。
【発明を実施するための形態】
【００２６】
図１は、本発明による方法の１つの例示的実施形態を、簡略化した概略的なフロー図で示している。この方法は、イメージングプレートを露光するステップ（Ｓ１）を含んでいる。イメージングプレートを露光する際、イメージングプレートを特定の線量で露光する場合が有利である。イメージングプレートが、記録デバイスから予め規定された間隔内に位置決めされる場合が有利と分かった。実際には、ＣＤの直径として調整され得る１２ｃｍが適していると分かった。Ｘ線機器のための調整としては、例えばモジュール式の歯（
2020537157.tif 000005
）のための調整を使用できる。
【００２７】
次のステップでは、イメージングプレートを読み取り、かつデジタル画像を生成する（Ｓ２）。読み取りの後または最中に、イメージングプレートが露光された線量を確定できる。線量が高く選択されすぎたかまたは低く選択されすぎたか、操作者への出力を行うことができる。信号対雑音比は線量によって変化し得るので、場合によっては、操作者はイメージングプレートを改めて露光しなければならない。
【００２８】
次のステップでは、イメージングプレートの品質を決定するために画像を前処理する（Ｓ３）。前処理のステップ（Ｓ３）は以下で詳しく解説する。
【００２９】
前処理（Ｓ３）の後、２つのプロセスブロックを同時に並行して、相次いで、または／および互いに対して代替的に実施することができる。
【００３０】
前処理された画像を、信号対雑音比が特定の閾値より低いピクセルがマーキングされているマスクに変換できる（Ｓ４）。その代わりにまたはそれに加えて、前処理された画像においてエッジ認識を実施でき（Ｓ５）、このエッジ認識に基づいて、このようにして認識されたピクセルがマーキングされているマスクを生成できる。
【００３１】
この両方のマスクにより、マーキングされたピクセルの総数を決定できる（Ｓ６）。このマーキングされたピクセルの数と、（除外マスクのピクセルの数だけ減らした）すべてのピクセルの総数とで比を出せば、この比に基づいて、イメージングプレートの品質を決定できる（Ｓ７）。
【００３２】
以下では、読み取り工程から得られた画像の前処理（Ｓ３）をより詳しく論じる。図２は、前処理（Ｓ３）の工程の可能な一実施形態を示している。読み取られた画像に対し、周縁マスクが提供される（Ｓ３１）。この周縁マスクは、例えば画像の周縁を取り囲む間隔であることができ、この間隔の外側では、画像データが放棄されるかまたは処理のために考慮されない。イメージングプレートのための、例えば口腔のための一般的なサイズ（２×３ｃｍ〜約６×８ｃｍ）の場合、１〜２ｍｍで十分であり、この１〜２ｍｍは、取り囲んでいる周縁として制約のないままであることが望ましい。この間隔は、固定的に調整されていてよく、または操作者によって自由に選択可能であってよい。
【００３３】
その後のステップでは、有り得るマーキングが検出され、このマーキングは、イメージングプレートに付けられており、本来の画像内容には属さない（Ｓ３２）。この検出は、例えばモルフォロジー演算によって行うことができ、例えば画像内容中の公知の構造を検索できる。認識されたマーキングまたは類似の構造は、これをもって除外マスクとなる周縁マスクに付加される（Ｓ３３）。
【００３４】
こうして作成された周縁マスクまたは除外マスクを、チェックまたは／および情報のために操作者に表示してもよい（Ｓ３４）。
【００３５】
こうして作成された除外マスクまたは周縁マスクが、読み取られた画像に適用され（Ｓ３５）、これにより前処理された画像が生じる。この画像が次のステップ（Ｓ４、Ｓ５）のベースとなる。
【００３６】
図３は、信号対雑音比の計算によるイメージングプレートのもう使用できないピクセルの数を確定するための方法（Ｓ４）の一実施形態の幾つかの局面をフロー図で示している。方法（Ｓ４）は、まず第１に、局所的な信号対雑音比（「Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  Ｒａｔｉｏ」、ＳＮＲ）の実施（Ｓ４１）を含んでいる。その際、平均化エリアが定められ、この平均化エリア内では、そこに内包されるすべてのピクセルに対してＳＮＲが１つの値として確定される。こうして確定された値が、そこに内包されるピクセルの１つに相応に割り当てられる。この計算ステップを、画像内の可能なすべての位置で実施する。平均化エリアを事実上ピクセルごとにさらにずらし、新たな場所でそれぞれ１つの局所的ＳＮＲを計算する。これにより、（平均化エリアの選択されたサイズに相応する周縁帯を除いて）画像内のほぼすべてのピクセルに対して１つのＳＮＲが生じる。これらの値は、前処理された画像のＳＮＲ画像として理解でき、かつ示すことができる（Ｓ４２）。実際には、上記のイメージングプレートサイズおよび相応の画像サイズの場合、２１×２１ピクセルの平均化エリアが好適と分かった。そのほかの平均化エリアサイズ、例えば半分だけの大きさまたは二倍の大きさも考えられる。
【００３７】
その後のステップ（Ｓ４３）では、追加的な大きさとしてＳＮＲ範囲を確定する（Ｓ４３）。これに関し、ＳＮＲ範囲の計算には局所的ＳＮＲ値を考慮することができる。例えば最小値および最大値を確定できる。その代わりに、これらの局所的ＳＮＲ値から、または前処理された画像のみに基づいて、大域的ＳＮＲ値を算出できる。
【００３８】
その後のステップ（Ｓ４４）では、ＳＮＲ閾値を決定する。ＳＮＲ閾値は、ＳＮＲ範囲を指針とすることができる。その代わりにまたはそれに加えて１つの局所的ＳＮＲ値をＳＮＲ閾値の決定に加えることができる。ＳＮＲ閾値は、例えば以下のように計算できる。
【００３９】
【数３】
2020537157.tif 000006
式中、ＳＮＲ局所（ｘ，ｙ）は局所的ＳＮＲ値を表し、ＳＮＲ最小は最小ＳＮＲ値を表し、ＳＮＲ最大は最大ＳＮＲ値を表す。この値が２０％以下であれば、対応するピクセルを欠陥があるとマーキングする。ただしそのほかのパーセント値、例えば１０％または３０％も考えられる。
【００４０】
これで、こうして決定されたＳＮＲ閾値により、ＳＮＲ画像内の、ＳＮＲがＳＮＲ閾値より低いピクセルを確定できる（Ｓ４５）。
【００４１】
こうしてマーキングされたピクセルが、その後のステップ（Ｓ４６）で、ＳＮＲに基づいて確定された、イメージングプレートのもう使用できない領域をマーキングするマスクへとまとめられる。
【００４２】
図４は、エッジ認識法に基づくイメージングプレートのもう機能を果たす能力のない領域を決定するための方法（Ｓ５）の代替的な一実施形態をフロー図で示している。この実施形態は、信号対雑音比を確定する方法に対して同時にまたは代替的に用いることができる。
【００４３】
第１のステップでは、前処理された画像にボケを付加する（Ｓ５１）。ボケ付加は、個々の強く雑音を発するピクセルに、エッジ認識アルゴリズムが強く反応しすぎることを回避するのに役立つ。
【００４４】
その後、（Ｓ５２）こうして準備された画像においてエッジ認識を実施する。この場合、例えばキャニーアルゴリズムを用いることができる。ただしそのほかのエッジ認識アルゴリズムも考えられる。エッジ認識アルゴリズムの結果もまた、閾値演算（Ｓ５３）により二値画像（画像値として０または１しか存在しない）に変換される画像である。閾値の選択は、例えば全体的な画像品質、雑音レベル、読み取り工程の質などに依存する。
【００４５】
このようにして変換されたエッジ認識画像においては、検出された構造のより良好な可視性のため、エッジ内の小さな隙間または穴を埋めるためのモルフォロジー膨張およびクロージング演算を実施できる（Ｓ５４）。こうして第２のマスクが生じ、この第２のマスクも、画像内、したがって間接的にイメージングプレート内のもう規定通りに機能しない領域をマーキングする。
【００４６】
このように二種類のやり方で生成された画像マスクは、例えば論理ＯＲ演算により相互に結合させることができる。もううまく機能しないピクセルの形態での損傷の数の取得は、すでにずっと上で詳細に解説したように、下式に基づいて、イメージングプレートの評価に換算できる。
【００４７】
【数４】
2020537157.tif 000007
および
【数５】
2020537157.tif 000008
【００４８】
重み付け係数ｗ0およびｗ1は、品質検査法をそれぞれのイメージングプレートの品質格付けの所望のスケーリングに適合するために使用でき、例えばトレーニングデータセットでの線形回帰により、生じる最終値がプレート品質の専門家の評価と可能な限り一致するように確定できる。
【００４９】
上で示した方法と専門家（歯医者、技術者）の格付け結果とを用いて品質格付けを比較すると、両方の格付けの十分に高い相関が判明した。
【００５０】
図５は、イメージングプレート１２を読み取るためのスキャン装置１０を示しており、イメージングプレート１２は、Ｘ線放射によって励起された準安定の蓄積中心の形態で、Ｘ線潜像を有している。スキャン装置１０は、イメージングプレート１２用の支持装置１４を有している。例えば、イメージングプレート１２を支持装置１４上で、負圧により、一般的に柔軟なイメージングプレート１２が支持面１４に平らに添うように固定することができる。スキャン装置１０は、読み取り光源としてレーザ１６をさらに含んでおり、レーザ１６は、赤色の波長をもつ読み取り光線１８を生成し、この読み取り光線１８により、イメージングプレート１２の準安定の蓄積中心を蛍光に励起できる。この蛍光２０は、典型的には青色である。
【００５１】
スキャン装置１０のこの実施形態ではレーザ１６は、読み取り光線１８を、制御可能な偏向ユニットに向けるように配置されている。制御可能な偏向ユニットは、ここではミラー２２として形成されている。そのほかの偏向ユニット、例えば光学系またはその類似物も考えられる。ミラー２２は、マイクロミラーとして、とりわけＭＥＭＳコンポーネントとして実施でき、こうして、ミラー２２と支持装置１４の間の相対運動なしでまたは少しの相対運動だけで、イメージングプレート１２の面を走査できる。その代わりに、ミラー２２をドラムスキャナ用の回転するミラーとして設けてもよい。この場合には、支持装置１４とミラー２２の間の相対運動は、搬送装置（不図示）によって実現される。イメージングプレート１２の読み取りのさらなる技術は、回転するペンタプリズムの使用を含んでいる。
【００５２】
スキャン装置１０は、図面では破線で示唆されたリフレクタ２４をさらに含むことができ、このリフレクタ２４は、イメージングプレート１２の周りの測定空間全体を、光を通さないように包囲しており、これにより、イメージングプレート１２から出ていく蛍光２０が光検出器２６へと反射される。散乱した読み取り光１８が光検出器２６に達するのを阻止するため、適切な措置、例えばダイクロイックフィルタ材料を用いることができる。読み取り工程を制御するため、スキャン装置１０は制御ユニット２８を含んでおり、制御ユニット２８は、制御機能だけでなく、評価機能または修正機能も果たすことができる。ここで示した実施形態では、制御ユニット２８は、イメージングプレートの品質を確定するための前述の方法の１つを実施するために適応されている。
【図１】
2020537157.tif 000009
【図２】
2020537157.tif 000010
【図３】
2020537157.tif 000011
【図４】
2020537157.tif 000012
【図５】
2020537157.tif 000013
【国際調査報告】
2020537157.tif 000014
2020537157.tif 000015
2020537157.tif 000016
2020537157.tif 000017
2020537157.tif 000018
2020537157.tif 000019