(57)【要約】  対象物のボリュームデータ画像を再構成する方法は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、受信された円状投影データから前処理された円状投影データを生成し、対象物に関する直線状経路に沿って検出器によって収集された直線状投影データを受信し、受信された直線状投影データから前処理された直線状投影データを生成し、前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて対象物に関する再構成された円状経路ボリューム画像を生成し、前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて対象物に関する再構成された直線状経路ボリューム画像を生成し、再構成された円状経路ボリューム画像と再構成された直線状経路ボリューム画像とを結像して、対象物のボリューム画像を生成する。装置とコンピュータプログラム製品もまた記述される。【選択図】図１
2010510833.tif 000002
【特許請求の範囲】
【請求項１】
対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成するデータ収集部と、
前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成する円状投影データ再構成部と、
前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成する直線状投影データ再構成部と、
前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する画像ボリューム処理部と、
を具備するボリューム画像再構成装置。
【請求項２】
前記直線状経路は、前記円状経路の平面に対して垂直でない、請求項１記載の装置。
【請求項３】
前記円状投影データ再構成部は、バックプロジェクション部と方向性フィルタリング部とを具備し、
前記バックプロジェクション部は、フィルタリングされた検出器空間データに基づいて前記再構成された円状経路ボリューム画像を生成し、
前記方向性フィルタリング部は、
前記前処理された円状投影データを、予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、リビニングされたデータを生成するリサンプリング部と、
前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、
前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記フィルタリングされた検出器空間データを生成する逆リビニング部と、を備える、
請求項１記載の装置。
【請求項４】
前記円状投影データ再構成部は、
前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、
前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、
前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
請求項１記載の装置。
【請求項５】
前記円状投影データ再構成部は、
前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、
前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するハイブリッドフィルタリング部と、
前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
請求項１記載の装置。
【請求項６】
前記データ収集部は、さらに、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された円状投影データを受信し、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、請求項１記載の装置。
【請求項７】
前記直線状投影データ再構成部は、
前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、
前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、
前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
請求項１記載の装置。
【請求項８】
前記直線状投影データ再構成部は、
前記前処理された直線状投影データを、２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、ダウンサンプリングされたデータを生成するダウンサンプリング部と、
前記ダウンサンプリングされたデータの微分係数を計算し、微分係数データを生成する微分器と、
前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、
前記フィルタリングされた検出器空間データをバックプロジェクションし、ダウンサンプリングされた直線状経路ボリューム画像を生成するバックプロジェクション部と、
前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングし、前記再構成された直線状経路ボリューム画像を生成するアップサンプリング部と、を備える、
請求項１記載の装置。
【請求項９】
前記対象物に電流に対応する放射線を照射するものであって、前記電流は２０ｍＡよりも低い又は等しいＸ線管と、
前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する検出器と、をさらに備える、
請求項１記載の装置。
【請求項１０】
前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整する寝台不正確さ補償部をさらに備え、
前記直線状投影データ再構成部は、さらに前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、
請求項１記載の装置。
【請求項１１】
前記直線状投影データ再構成部は、
前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、
前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
請求項１記載の装置。
【請求項１２】
前記対象物を放射線で制御可能に照射する線源と、
前記対象物が前記放射線で照射された場合、前記直線状投影データを生成する検出器と、
前記線源を制御して前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射する動的コリメーション部と、をさらに備える、
請求項１記載の装置。
【請求項１３】
前記データ収集部は、さらに、異なる捕獲時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる捕獲時間にそれぞれ対応する複数の前処理された円状投影データを生成し、
前記画像ボリューム処理部は、さらに、前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成するものであり、
前記対象物の前記複数のボリューム画像を前記画像ボリューム処理部から統合し、前記複数のボリューム画像の統合されたデータを生成する統合部と、
前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する合成部と、をさらに備える、
請求項１記載の装置。
【請求項１４】
前記対象物は繰り返し位相を示すものであり、
前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、
前記合成部は、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、
請求項１３記載の装置。
【請求項１５】
前記対象物は繰り返し位相を示すものであり、
前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、
前記データ収集部は、さらに、異なる獲得時間に円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、
請求項１記載の装置。
【請求項１６】
前記データ収集部は、さらに、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、
前記画像ボリューム処理部は、さらに、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記対象物の前記部分の前記複数のボリューム画像を合成し、前記対象物の前記ボリューム画像を生成する、
請求項１記載の装置。
【請求項１７】
対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、
前記受信された円状投影データから前処理された円状投影データを生成し、
前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、
前記受信された直線状投影データから前処理された直線状投影データを生成し、
前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、
前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、
前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、
を特徴とする対象物のボリューム画像を再構成する方法。
【請求項１８】
前記直線状経路は、前記円状経路の平面に対して垂直でないこと、を特徴とする請求項１７記載の方法。
【請求項１９】
前記再構成された円状経路ボリューム画像を生成することは、さらに
リビニングされたデータを生成するために、前記前処理された円状投影データを予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、
フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、
前記フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションすること、
を特徴とする請求項１７記載の方法。
【請求項２０】
前記再構成された円状経路ボリューム画像を生成することは、さらに
プレ重みづけされたデータを得るために、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、
フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、
前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションすること、
を特徴とする請求項１７記載の方法。
【請求項２１】
前記再構成された円状経路ボリューム画像を生成することは、さらに、
プレ重みづけされたデータを得るために、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、
フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってハイブリッドフィルタリングし、
前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションすること、
を特徴とする請求項１７記載の方法。
【請求項２２】
前記円状投影データを受信することは、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された前記円状投影データを受信し、
前記前処理された円状データを生成することは、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成すること、
を特徴とする請求項１７記載の方法。
【請求項２３】
前記再構成された直線状経路ボリューム画像を生成することは、さらに、
微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、
リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
フィルタリングされたデータを生成するために、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、
フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションすること、
を特徴とする請求項１７記載の方法。
【請求項２４】
さらに、
ダウンサンプリングされたデータを生成するために、前記前処理された直線状投影データを２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、
微分係数データを生成するために、前記ダウンサンプリングされたデータの微分係数を計算し、
リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、
フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
ダウンサンプリングされた直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションし、
前記再構成された直線状経路ボリューム画像を生成するために、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングすること、
を特徴とする請求項１７記載の方法。
【請求項２５】
さらに、
前記対象物に、Ｘ線管内の電流であって２０ｍＡよりも低い又は等しい前記電流に対応する放射線を照射し、
前記対象物が前記放射線に照射された場合、前記直線状投影データを生成すること、
を特徴とする請求項１７記載の方法。
【請求項２６】
さらに、
前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整することによって寝台不正確さ補償し、
前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成すること、
を特徴とする請求項１７記載の方法。
【請求項２７】
前記直線状経路ボリューム画像を生成することは、さらに、
微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、
リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、
前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションすること、
を特徴とする請求項１７記載の方法。
【請求項２８】
さらに、
前記対象物を放射線で制御可能に照射し、
前記対象物が前記放射線で照射された場合、前記直線状投影データを生成し、
前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射するように前記照射することを制御すること、
を特徴とする請求項１７記載の方法。
【請求項２９】
前記円状投影データを前記受信することは、異なる獲得時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記結合することは、前記前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成することであり、
さらに、
前記対象物の前記複数のボリューム画像を統合し、
前記複数のボリューム画像の統合されたデータを生成し、
前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成すること、
を特徴とする請求項１７記載の方法。
【請求項３０】
前記対象物は繰り返し位相を示すものであり、
さらに、
前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、
前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成すること、
を特徴とする請求項２９記載の方法。
【請求項３１】
前記対象物は繰り返し位相を示すものであり、
さらに、
前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、
前記円状経路の複数部分に沿って、異なる獲得時間に前記検出器によって収集された前記円状投影データを受信し、
異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、
前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成すること、
を特徴とする請求項１７記載の方法。
【請求項３２】
前記円状投影データを前記受信することは、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、
前記前処理された円状投影データを前記生成することは、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、
前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、
前記合成することは、前記対象物の前記ボリューム画像を生成するために、前記対象物の前記部分の前記複数のボリューム画像を合成する、
ことを特徴とする請求項１７記載の方法。
【請求項３３】
コンピュータにより実行された場合に、
対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、
前記受信された円状投影データから前処理された円状投影データを生成し、
前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、
前記受信された直線状投影データから前処理された直線状投影データを生成し、
前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、
前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、
前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、
を具備するステップをコンピュータに実行させるコンピュータプログラムインストラクションを記憶するコンピュータ読み取り可能媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概して医用画像の再構成に関する。より詳細には、本発明は、医用画像の質と円状及び直線状軌道スキャンからの投影データに基づく３次元コーンビームコンピュータ断層撮影（ＣＴ）再構成の効率との向上のための新しい方法及びシステム装置に関する。
【０００２】
本発明は、下記の参考文献のリスト内で参照され記載されている様々な技術の使用を含む。これら参考文献は、括弧内の対応する参照番号により明細書中で引用され、参照によりここに本質的に組込まれる。
【表１】
2010510833.tif 000003
【背景技術】
【０００３】
コンピュータ断層撮影（ＣＴ）装置により生成された再構成画像の質と効率とは、ＣＴ装置の全般的な有効性にとって重要である。画像を再構成する際に用いられるアルゴリズムは、質と効率とに影響を与える。
【０００４】
それに加え、円状収集（circular acquisition）における１回転で全ての人体器官を覆うＣＴスキャナに対する臨床的な要求がある。円状軌道は、心臓／冠状動脈撮像やパーフュージョン等のような動的検査においてヘリカルよりも優位である。例えば、人の心臓を覆うことは、アキシャル方向に関して１０〜１２センチメートルの検出器サイズを要求することを明記しておく。
【０００５】
たった１回のスキャンでほとんどの人体器官をスキャンするために、大ボリューム範囲を備えるスキャナが要求される。このような大ボリューム範囲は、医用ＣＴにおいて広く用いられているフェルドカンプ（Feldkamp）、デイビス（Davis）、クレス（Kress）のアルゴリズム（ＦＤＫ）に基づく慣習的な近似又は非厳密な再構成法に対する挑戦である。
【０００６】
ＦＤＫは、高い達成可能空間解像度（high achievable spatial resolution）を有する円状軌道Ｘ線ＣＴデータを再構成するための効率的な方法である。一方、他のアルゴリズムは、リビニング（ribinning）のような空間解像度を低下させる幾つかのステップを含む。しかしながら、ＦＤＫのコーンビームジオメトリは、円状軌道の理論的な不完全性、特に体軸方向に関する急な減弱変化がおきる場所、に起因するコーンビームアーチファクトを受ける。
【０００７】
一方、厳密なフィルタードバックプロジェクション（ＦＢＰ）再構成アルゴリズムが提案されている［１１］。Ｘ線軌跡が単一の円を構成している場合、コーンビームの厳密な再構成は不可能であることが知られている。なぜならば、軌跡に交差しない多くの平面が存在し、これら平面の間のデータが計測できないからである。換言すれば、例えば図１に示すように、円状データ１２は、空間１０を覆う不完全なラドン空間を与える。この制限を解消するため、図２に示すように、円状スキャン軌跡２２を直線状スキャン軌跡２０で補うことが提案され、様々な円及び直線状アルゴリズムが提案されている［６〜１０］。しかしながら、［６］〜［１０］で提案されたアルゴリズムは、商業的なＣＴスキャナには全く適していなかった。すなわち、２つの要因：（ａ）それらアルゴリズムの非能率な構成に起因する画像再構成のために長時間を要求する（有効なshift-invariantなフィルタードバックプロジェクション型ではなかった）、（ｂ）製造に実用的でない、とても大きい検出器を要求する、すなわち、同等には、いわゆる“長体問題（long object problem）”（これは、人体スキャンにおいてよくあることであり、検出器の外側にまで亘る対象物（体軸方向に関して打ち切られる）は、再構成可能できないことを意味する）を受ける。
【０００８】
カツェビッチ（Kathevich）の厳密な直線状＋円状再構成法は、円状及び直線状データの両方に対するshift-invariant（ピクセル依存なしに独立して行なわれた、再構成された対象物のコンボルーション）と、フィルタードバックプロジェクション（ＦＢＰ）とを含む。
【０００９】
しかしながら、カツェビッチによるアルゴリズムの厳密な方法は、現状のＣＴスキャナのハードウェアに適用するのが容易でない。なぜならば、ヒルベルトコンボルーションは周波数領域で適用され、ヒルベルトコンボルーションは検出されたデータに対する微分を要求するからである。これは、空間解像度の低下を起こしうる。それに加え、ヒルベルト変換は、リビニングを要求する。これは、計算上の効率を相当に低下させ、解像度を多少低下させる。このように、円状及び直線状データの両方に対するそのようなアルゴリズムは、多くの時間を消費しうり、付加的な高価な処理資源を要求し、そして商業的なＣＴスキャナに適さない。
【発明の開示】
【００１０】
本発明に係るある目的は、再構成画像の質を向上することとコーンビームアーチファクトを除去することとである。次節において我々は、コーンビームアーチファクトを除去するための、円状スキャンと直線状スキャンとを組合せる方法を記述する。
【００１１】
従って、上述の技術における再構成アルゴリズムの問題を解消するため、本発明は、方法、システム装置、及びコンピュータプログラム製品を提供する。
【００１２】
本発明のある目的は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成するデータ収集部と、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成する円状投影データ再構成部と、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成する直線状投影データ再構成部と、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する画像ボリューム処理部と、を具備する新規のボリューム画像再構成装置を提供することにある。
【００１３】
本発明の他の目的は、前記直線状経路が前記円状経路の平面に対して垂直でない、新規の装置を提供することにある。
【００１４】
本発明の他の目的は、前記円状投影データ再構成部が、バックプロジェクション部と方向性フィルタリング部とを具備し、前記バックプロジェクション部が、フィルタリングされた検出器空間データに基づいて前記再構成された円状経路ボリューム画像を生成し、前記方向性フィルタリング部が、前記前処理された円状投影データを、予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、リビニングされたデータを生成するリサンプリング部と、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記フィルタリングされた検出器空間データを生成する逆リビニング部と、を備える、新規の装置を提供することにある。
【００１５】
本発明の他の目的は、前記円状投影データ再構成部が、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
【００１６】
本発明の他の目的は、前記円状投影データ再構成部が、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するハイブリッドフィルタリング部と、前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
【００１７】
本発明の他の目的は、前記データ収集部が、さらに、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された円状投影データを受信し、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、新規の装置を提供することにある。
【００１８】
本発明の他の目的は、前記直線状投影データ再構成部が、前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
【００１９】
本発明の他の目的は、前記直線状投影データ再構成部が、前記前処理された直線状投影データを、２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、ダウンサンプリングされたデータを生成するダウンサンプリング部と、前記ダウンサンプリングされたデータの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、前記フィルタリングされた検出器空間データをバックプロジェクションし、ダウンサンプリングされた直線状経路ボリューム画像を生成するバックプロジェクション部と、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングし、前記再構成された直線状経路ボリューム画像を生成するアップサンプリング部と、を備える新規の装置を提供することにある。
【００２０】
本発明の他の目的は、前記対象物に電流に対応する放射線を照射するものであって、前記電流は２０ｍＡよりも低い又は等しいＸ線管と、前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する検出器と、をさらに備える新規の装置を提供することにある。
【００２１】
本発明の他の目的は、前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整する寝台不正確さ補償部をさらに備え、前記直線状投影データ再構成部が、さらに前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、新規の装置を提供することにある。
【００２２】
本発明の他の目的は、前記直線状投影データ再構成部が、前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
【００２３】
本発明の他の目的は、前記対象物を放射線で制御可能に照射する線源と、前記対象物が前記放射線で照射された場合、前記直線状投影データを生成する検出器と、前記線源を制御して前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射する動的コリメーション部と、をさらに備える新規の装置を提供することにある。
【００２４】
本発明の他の目的は、前記データ収集部が、さらに、異なる捕獲時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる捕獲時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記画像ボリューム処理部が、さらに、前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成するものであり、前記対象物の前記複数のボリューム画像を前記画像ボリューム処理部から統合し、前記複数のボリューム画像の統合されたデータを生成する統合部と、前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する合成部と、をさらに備える新規の装置を提供することにある。
【００２５】
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、前記合成部が、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、新規の装置を提供することにある。
【００２６】
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、前記データ収集部が、さらに、異なる獲得時間に円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、新規の装置を提供することにある。
【００２７】
本発明の他の目的は、前記データ収集部が、さらに、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、前記画像ボリューム処理部は、さらに、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記対象物の前記部分の前記複数のボリューム画像を合成し、前記対象物の前記ボリューム画像を生成する、新規の装置を提供することにある。
【００２８】
本発明の他の目的は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成し、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する、対象物のボリューム画像を再構成する新規の方法を提供することにある。
【００２９】
本発明の他の目的は、前記直線状経路が、前記円状経路の平面に対して垂直でない、新規の方法を提供することにある。
【００３０】
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらに、リビニングされたデータを生成するために、前記前処理された円状投影データを予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、前記フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションする、新規の方法を提供することにある。
【００３１】
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらにプレ重みづけされたデータを得るために、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションする、新規の方法を提供することにある。
【００３２】
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらに、プレ重みづけされたデータを得るために、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってハイブリッドフィルタリングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションする、新規の方法を提供することにある。
【００３３】
本発明の他の目的は、前記円状投影データを受信することは、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状データを生成することは、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、新規の方法を提供することにある。
【００３４】
本発明の他の目的は、前記再構成された直線状経路ボリューム画像を生成することは、さらに、微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションする、新規の方法を提供することにある。
【００３５】
本発明の他の目的は、さらに、ダウンサンプリングされたデータを生成するために、前記前処理された直線状投影データを２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、微分係数データを生成するために、前記ダウンサンプリングされたデータの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、ダウンサンプリングされた直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するために、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングする、新規の方法を提供することにある。
【００３６】
本発明の他の目的は、さらに、前記対象物に、Ｘ線管内の電流であって２０ｍＡよりも低い又は等しい前記電流に対応する放射線を照射し、前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する、新規の方法を提供することにある。
【００３７】
本発明の他の目的は、さらに、前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整することによって寝台不正確さ補償し、前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、新規の方法を提供することにある。
【００３８】
本発明の他の目的は、前記直線状経路ボリューム画像を生成することは、さらに、微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションする、新規の方法を提供することにある。
【００３９】
本発明の他の目的は、さらに、前記対象物を放射線で制御可能に照射し、前記対象物が前記放射線で照射された場合、前記直線状投影データを生成し、前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射するように前記照射することを制御する、新規の方法を提供することにある。
【００４０】
本発明の他の目的は、前記円状投影データを前記受信することは、異なる獲得時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記結合することは、前記前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成することであり、さらに、前記対象物の前記複数のボリューム画像を統合し、前記複数のボリューム画像の統合されたデータを生成し、前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する、新規の方法を提供することにある。
【００４１】
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、さらに、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、新規の方法を提供することにある。
【００４２】
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、さらに、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、前記円状経路の複数部分に沿って、異なる獲得時間に前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、新規の方法を提供することにある。
【００４３】
本発明の他の目的は、前記円状投影データを前記受信することは、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記合成することは、前記対象物の前記ボリューム画像を生成するために、前記対象物の前記部分の前記複数のボリューム画像を合成する、新規の方法を提供することにある。
【００４４】
本発明の他の目的は、コンピュータにより実行された場合に、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成し、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、を具備するステップをコンピュータに実行させるコンピュータプログラムインストラクションを記憶するコンピュータ読み取り可能媒体を提供することにある。
【００４５】
ハードウェア開発の影響とコストとを最小にするため、本発明は、現行の装置に対する最小の変化で実行できうる円状データのランプ再構成を含む。それに加え、ランプフィルタに基づく再構成アルゴリズム（例えば、ＦＤＫ、またはザミャチン（Zamyatin）、タグチ（Taguchi）、シルバー（Silver）（ＺＴＳ）ランプフィルタに基づく）は、慣習的なＣＴ医用撮像装置でよく用いられうるので、それぞれが撮像対象の体部分に基づく特定のアプリケーションのために用意されているランプ関数の確立されたライブラリーが存在する。このように、本発明は、このランプ関数の現存するライブラリーからのランプ関数を利用可能であり、異なるアプローチに関するフィルタ関数に特有なアプリケーションの開発（例えば、ヒルベルト関数に対する体部分に特有のフィルタリングアプローチに基づくアプローチ）に関連づけられる不確定とコストとを避ける。
【００４６】
直線状データは、不完全な円状軌道に起因するコーンビームアーチファクトを相殺する追加項を提供する。
【００４７】
本アプローチは、擬人ファントムのシミュレーションデータと実測のスキャンデータとに基づいて実現可能であると証明されている。本アプローチは、スキャン中における造影剤注入、体動、及びミスアライメントに関して安定している証明されている。直線状スキャンに起因して患者が照射される余分な線量もまた、円状スキャンのそれに比して有利に少ない。本発明は、コーン角がとても大きい場合でさえも、コーンビームアーチファクトの影響のない再構成をも可能とし、さらに検出器列の数が２５６を超える場合であっても使用可能である。
【００４８】
本発明の実施形態によれば、円状データのランプフィルタに基づく再構成（例えば、ＦＤＫ又はＺＴＣのような）が行なわれる。さらに、低線量の直線状スキャン（ヒルベルト変換を含む方法を用いて再構成される）が追加される。直線状データは、再構成を速めるため、及びさらなる処理改善（例えば、メモリ削減や処理容量）を得るため、２から４の係数によってダウンサンプリングされうる。このように、本発明は、商業的なＣＴスキャナに適している。
【００４９】
幾つかの臨床応用は、円状軌道に対して直線状部分を傾かせるガントリチルトを要求する。換言すれば、その直線状部分は、円状軌道の平面に対して垂直でない。このように、本発明の実施形態は、ガントリが傾いている場合、例えば、脳撮像期間中に適している。この場合、患者の目と患者の耳の中心とを結ぶ眼窩耳孔（ＯＭ）線に対して再構成面が平行するように、スキャン前にガントリが傾かれる。このＯＭ線を含む面は、基準解剖面として供される。また、ガントリを傾けることは、患者の目のレンズへの照射線量を最小にするのに役立つ。
【００５０】
本発明の他の方法、システム、及びコンピュータプログラム製品は、以下の図面及び好ましい実施形態の詳細な記述の検討のもと、一般人あるいは当業者にとって明らかになろう。
【００５１】
本発明のより完全な理解及びそれの多くの付随する利点は、以下の詳細な記述を、以下の添付図面に関連づけて考慮して参照することによってより容易に理解できるようになるであろう。
【図面の簡単な説明】
【００５２】
【図１】図１は、ラドン空間における円状データ範囲の一例の図である。
【図２】図２は、直線状スキャン軌跡と円状スキャン軌跡との間の関係の一例の等角投影図である。
【図３】図３は、本発明の実施のブロック図である。
【図４】図４は、本発明の実施に係る物理座標系のブロック図である。
【図５】本発明の実施に係る画像再構成の実行方法のフローチャートである。
【図６】本発明の実施に係る直線状ボリューム画像の再構成の実行方法のフローチャートである。
【図７Ａ】図７Ａは、円状スキャン中において収集されたデータに対する、検出器上のフィルタリングラインの一例である。
【図７Ｂ】図７Ｂは、直線状スキャン中において収集されたデータに対する、検出器上のフィルタリングラインの一例である。
【図８】図８は、直線状スキャンの範囲の一例のブロック図である。
【図９】図９は、平面及び湾曲した検出器のジオメトリを考慮するためのさらなる用語を有するブロック図。
【図１０】図１０は、本実施形態の一例に係るｗ対ρのグラフ上におけるＢＰＪエリアのプロットである。
【図１１Ａ】図１１Ａは、傾いたガントリを含む本発明の実施形態の一例の等角投影図である。
【図１１Ｂ】図１１Ｂは、円状軌道に関連付けられた、傾いた座標の等角投影図である。
【図１２Ａ】図１２Ａは、本発明の実施形態に係る座標の等角投影図である。
【図１２Ｂ】図１２Ｂは、本発明の実施形態に係る座標の他の図である。
【図１３Ａ】図１３Ａは、慣習的な円状ＦＤＫ法を用いた、シミュレーションデータの再構成画像の一例である。
【図１３Ｂ】図１３Ｂは、直線状データ再構成に組合わされたランプに基づく円状データ再構成を有する本発明の実施形態を用いた、シミュレーションデータの再構成画像の一例である。
【図１４Ａ】図１４Ａは、ゲート制御された円状スキャンデータをのみを用いて画像が再構成された場合の心像動の影響を示す再構成画像である。
【図１４Ｂ】図１４Ｂは、本発明の実施形態に係るゲート制御された円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された場合の心臓動の影響を示す再構成画像である。
【図１４Ｃ】図１４Ｃは、図１４Ａと図１４Ｂとの間の差分の画像である。
【図１４Ｄ】図１４Ｄは、ゲート制御されていない円状スキャンデータのみを用いて画像が再構成された場合のシンチレータ心臓動の影響を示す再構成画像である。
【図１４Ｅ】図１４Ｅは、ゲート制御されていない円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された場合の心臓動の影響を示す再構成画像である。
【図１４Ｆ】図１４Ｆは、図１４Ｄと図１４Ｅとの間の差分の画像である。
【図１５】図１５は、擬人ファントムの画像である。
【図１６Ａ】図１６Ａは、ＦＤＫ円状スキャンデータのみを用いる再構成画像の一例である。
【図１６Ｂ】図１６Ｂは、直線状再構成データの一例である。
【図１６Ｃ】図１６Ｃは、本発明の実施形態を用いる再構成画像の一例である。
【図１７Ａ】図１７Ａは、ＦＤＫ円状データのみを用いて再構成された画像の一例である。
【図１７Ｂ】図１７Ｂは、本発明の実施形態に係る、管電流５０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。
【図１７Ｃ】図１７Ｃは、図１７Ｂの例において用いられた直線状データのみの一例である。
【図１７Ｄ】図１７Ｄは、本発明の実施形態に係る、管電流３０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。
【図１７Ｅ】図１７Ｅは、図１７Ｄの例において用いられた直線状データのみの一例である。
【図１７Ｆ】図１７Ｆは、本発明の実施形態に係る、管電流５０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。
【図１７Ｇ】図１７Ｇは、図１７Ｆの例において用いられた直線状データのみの一例である。
【図１８Ａ】図１８Ａは、ＦＤＫ円状データ再構成のみを用いて再構成された画像である。
【図１８Ｂ】図１８Ｂは、本発明の実施形態に係る、正確な寝台位置情報で再構成された画像である。
【図１８Ｃ】図１８Ｃは、不正確な寝台位置情報で再構成された画像である。
【図１９Ａ】図１９Ａは、ＦＤＫ円状スキャンのみで得られたＦＯＶの図である。
【図１９Ｂ】図１９Ｂは、仮想的な拡張検出器を用いて得られた拡張ＦＯＶの図である。
【図２０】図２０は、本発明の実施形態に係る直線状スキャンコリメーションを用いる直線状スキャンの照射線量を低減する方法におけるＦＯＶの図である。
【図２１】図２１は、動的コリメーションに対するコリメーション窓の一例である。
【図２２】図２２は、本発明に係るジオメトリの更なる図である。
【図２３Ａ】図２３Ａは、直線状スキャンデータ中のフィルタラインのグラフである。
【図２３Ｂ】図２３Ｂは、直線状スキャンデータ中のフィルタラインの数のグラフである。
【図２４】図２４は、本発明の他の実施形態に係る直線状及び円状アルゴリズムの方法のブロック図である。
【図２５Ａ】図２５Ａは、ショートスキャンコンボルーションにおいて用いられるカーブのファミリーの一例である。
【図２５Ｂ】図２５Ｂは、ショートスキャンコンボルーションにおいて用いられるカーブの他のファミリーの一例である。
【図２５Ｃ】図２５Ｃは、図２５Ａにおけるカーブのファミリーの他の表現である。
【図２５Ｄ】図２５Ｄは、図２５Ｂにおけるカーブのファミリーの他の表現である。
【図２６Ａ】図２６Ａは、複数のファミリーからのラインにおけるフルスキャンコンボルーションの一例である。
【図２６Ｂ】図２６Ｂは、図２６Ａの２つのファミリーのラインの重なりの１次元図である。
【図２７Ａ】図２７Ａは、本発明の実施形態に係る逆（inverse）リビニングカーブの一例である。
【図２７Ｂ】図２７Ｂは、逆（reverse）リビニング処理の一例である。
【図２８Ａ】図２８Ａは、本発明の実施形態に係る拡張再構成関心ボリューム（ＶＯＩ）の一例である。
【図２８Ｂ】図２８Ｂは、ＦＤＫ―ＶＯＩを拡張させたＦＤＫ―ＶＯＩの形状の線図である。
【図２９Ａ】図２９Ａは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【図２９Ｂ】図２９Ｂは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【図２９Ｃ】図２９Ｃは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【図２９Ｄ】図２９Ｄは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【図２９Ｅ】図２９Ｅは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【図２９Ｆ】図２９Ｆは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
【発明を実施するための最良の形態】
【００５３】
では図面、より詳細には図３を参照せよ。種々の図面の全体に亘って、同一の参照符号は、同一又は対応する部分を示している。図３は、本発明に係るＸ線コンピュータ断層撮影装置の実施形態のブロック図である。撮像装置は、略錐形状のＸ線束のコーンビームを発生するＸ線源３を有するガントリ１と、２次元状に配列された複数の検出素子（すなわち、一列に配列された複数の素子が複数列に積まれている）を含む２次元アレイ型のＸ線検出器５とを含む。Ｘ線源３と２次元アレイ型のＸ線検出器５とは、回転リング２に設置され、スライド制の寝台６に載置されている被検体に対向して向き合っている。２次元アレイ型のＸ線検出器５は、回転リング２上に搭載されている。各検出素子は、検出されたデータの１チャンネルに対応する。Ｘ線源３からのＸ線は、Ｘ線フィルタ４を介して被検体に向けられる。被検体を透過したＸ線は、２次元アレイ型のＸ線検出器５によって電気信号として検出される。
【００５４】
Ｘ線制御装置８は、高電圧発生器７にトリガ信号を供給する。高電圧発生器７は、トリガ信号を受信したタイミングで高電圧をＸ線源３に印加し、Ｘ線をＸ線源３から放射させる。ガントリ／寝台制御器９は、ガントリ１の回転リング２の回転と寝台６のスライドシートのスライドとを同期制御する。システム制御器１０は、全システムのコントロールセンターを構成し、被検体から見てＸ線源３が円状スキャンと直線状スキャンとを実行するようにＸ線制御器８とガントリ／寝台制御器９とを制御する。
【００５５】
本実施形態に係る円状スキャン期間中、回転リング２（線源３と検出器５とを含む）は、スライド制の寝台６が所定の位置に留まったまま、所定の角速度で連続的に回転し、Ｘ線は、所定角度おきに連続的または間欠的にＸ線源３から放射される。このように、寝台６が静止している時、円状検出器５は、スキャンされている対象物に関する真の円状経路に沿って移動し、螺旋経路上を移動しない。直線状スキャン中、寝台６は、Ｘ線が連続的または間欠的にＸ線源３（静的な位置に保持されている、例えば、直接的に寝台６の上に）から放射されたまま、直線的にガントリ内を所定速度で移動される。
【００５６】
上記の直線状スキャンの実施形態において、撮像される患者又は対象物は、線源と検出器とが所定位置に保持されている間、ガントリ内を移動される。しかしながら、本発明は、対象物と検出器との両方が移動する実施形態と同様に、対象物が所定位置に保持され、検出器が対象物に沿って移動する実施形態をも含む。円状スキャンと直線状スキャンとは、あらゆる順番で交互に行なわれる。
【００５７】
上記の円状スキャンの実施形態において、検出器５の経路は、円であると示された。しかしながら、本発明は、円状経路に限定されず、対象物の周りを囲む他の湾曲された経路をも含む。
【００５８】
２次元配列型Ｘ線検出器５の出力信号は、投影データを生成するために、特に直線状スキャン中においては直線状投影データを生成し、円状スキャン中においては円状投影データを生成するために、チャンネル毎にデータ収集部１１により増幅され、デジタル信号に変換される。投影データは、当業者に知られている方法を用いて前処理され、前処理されたデータは、データ収集部１１から再構成処理部１２に出力される。再構成処理部１２は、前処理されたデータを用いて、スキャンされた対象物に関するボリューム画像を生成する。
【００５９】
再構成されたボリューム画像は、ディスプレイ装置１４に送信され、そこで３次元画像または断層像として視覚的に表示される。
【００６０】
この明細書の目的のため、我々は、ボリューム画像をスキャンされている身体物の描写であると定義するものとする。その画像は、幾つかの撮像技術によって発生される。撮像技術の例は、テレビジョン、ＣＣＤカメラ、又はＸ線、音波、あるいは超音波撮像装置を含みうる。画像が記録される初期媒体は、電子的固体装置、画像フィルム、又は輝尽性蛍光体のような幾つかの他の装置が可能である。記録された画像は、その後、電子的（ＣＣＤ信号の場合のように）又は機械的／光学的方法（写真フィルムをデジタル化する場合や輝尽性蛍光体からのデータをデジタル化する場合のように）の組合せによって、デジタル形式に変換され得る。
【００６１】
本実施形態に用いられる検出器５は、１２０ミリメートルのボリューム範囲を有する、２５６列の医用ＣＴスキャナ検出器を含みうる［１−４］。当業者は、本発明が他の慣習的なＣＴスキャナ検出器にも適用するということを認識するであろう。本実施形態に係る補足的なスキャナパラメータは、表１に示される。  
【表２】
2010510833.tif 000004
【００６２】
図４は、本撮像装置の物理座標システムの定義の一例のブロック図である。円状軌跡上の線源位置は、ｙＣ（β）＝（Ｒｃｏｓβ，Ｒｓｉｎβ，０）である。ここで、Ｒは円状軌跡の半径、βはビュー角である。直線状軌跡上の線源位置は、ｙＬ（ｈ）＝（Ｒｃｏｓθ，Ｒｓｉｎθ，ｈ）、ここでθは所定の直線角度、ｈはｚ軸に沿う直線状軌跡上の線源の垂直座標である。慣習的な完全アルゴリズム［１１］に対する方程式は、方程式（１）に示すような、円及び直線状データの積分にそれぞれ対応する２つの積分項の加算を含む。  
【数１】
2010510833.tif 000005
【００６３】
方程式（１）において、ｘは再構成点、γはファン角、νは垂直検出器座標である。ＨｓｉｎγＤｉ，ｉ＝１，２，は、カーネル１／ｓｉｎγを有する１次元方向ヒルベルトフィルタである。線及び円状スキャンに対するフィルタリング方向は、異なる。
【００６４】
一方、本発明によれば、方程式（１）の厳密なアルゴリズムは簡約されている。第１に、湾曲した検出器上のフィルタリング方向は、円状データに対して略水平である（それらは平面検出器上において厳密に水平）。従って、本発明によれば、円状データに対してリビニングが適用されない。第２に、γ微分係数に適用されるヒルベルトコンボルーションは、ランプコンボルーションに置き換え可能であることが提案されている［１２］。しかしながら、本発明は、さらに進んでおり、ヒルベルト再構成に比して速度と効率とを改善するためだけでなく、解像度を改善するために方程式（１）の第１項の全てをＦＤＫ再構成に置き換える（比較のため［１２］を参照見よ）。
【００６５】
本発明によれば、円状スキャンは、ＦＤＫ再構成アルゴリズム又はＺＴＳ再構成アルゴリズムのような、ランプフィルタを有する再構成アルゴリズムを用いて再構成される。さらに、本発明によれば、直線状スキャンデータは、再構成された円状スキャンデータを補足するための低周波数補正項の計算のために再構成される。
【００６６】
図５は、本発明の実施形態に係る画像再構成を実行する方法のフローチャートである。図５において、その方法は、例えば、前処理器又は前処理器を含む検出器からの前処理された円状投影データが受信されるステップＳ５００を含む。ステップＳ５０２において、受信されたデータは、円状投影ボリューム画像を生成するために再構成される。ステップＳ５０１において、前処理された直線状投影データは受信され、そしてステップＳ５０３において、前処理された直線状スキャン投影データは、直線状投影ボリューム画像を生成するため、フィルタードバックプロジェクション（ＦＢＰ）ヒルベルト再構成を用いて再構成される。ステップＳ５０４において、円状投影ボリューム画像と直線状投影ボリューム画像とは、検出された対象物の最終及び全ボリューム画像を生成するために組合される。
【００６７】
図６は、本発明に係る、例えば、図５におけるステップＳ５０３を実行する、直線状ボリューム画像の再構成を実行する方法のフローチャートである。図６において、ステップＳ６０１は、受信された前処理後の直線状投影データの微分を実行する。ステップＳ６０２において、微分されたデータは、フィルタリングラインに沿ってリビニングされ、ステップＳ６０３において、リビニングされたデータは、フィルタリングラインに沿ってヒルベルトフィルタでフィルタリングされる。ステップＳ６０４において、ヒルベルトフィルタリングされたデータは、検出器空間又は格子点に逆リビニングされ、ステップ６０５において、リビニングされたデータは、直線状投影ボリューム画像を生成するため、３次元バックプロジェクションされる。
【００６８】
図７Ａは、円状スキャン中に収集されたデータに対する検出器上のフィルタリングラインの一例を示す。図７Ｂは、直線状スキャン中に収集されたデータに対する検出器上のフィルタリングラインの一例を示す。円状データにおいて、フィルタリング方向は略水平であり、直線状データに対しては、フィルタリング方向は傾いている。
【００６９】
直線状データのコンボルーション。初めに、ＰＩラインの概念を紹介する。再構成ピクセルをｘとすれば、ＰＩラインは、円状軌跡上の数点と直線状軌道上の数点とを接続し、ピクセルｘに交差する直線である。ＰＩラインは、各再構成ピクセルｘに対して存在することが［１１］に示されている。平面検出器上のＰＩ窓（軌跡の直線状部分と円状部分との間をＰＩラインによって追跡された）は、放物線：νｆｌ＝−｛（ｈｕ２）／（２Ｒ２）｝−ｈ／２によって与えられる。ここで、ｕ及びνｆｌは、それぞれ水平及び垂直平面検出器座標である。放物線の頂点は、（０，−ｈ／２）に位置する。フィルタリングラインのパラメータとしてρ＝ｔａｎγを用い、そして平面検出器上のラインνｆｌ＝ｍｕ＋ｂが湾曲検出器上のカーブν＝ｂｃｏｓγ＋ｍＲｓｉｎγに投影されることを考慮すると、我々は、図７Ｂに示すように（図７Ｂの垂直軸はνに対応し、水平軸はγに対応する）、コンボルーションラインの下記の方程式を得る。  
【数２】
2010510833.tif 000006
【００７０】
フィルタリングファミリーは、検出器エリアの二重範囲（double coverage）を提供する。すなわち、ＰＩ窓内のあらゆるピクセルに対して、そのピクセルを通過する２つのフィルタリングラインを見つけられる。フルスキャンモードにおいては、両方のラインからの寄与が加えられる。ショートスキャンモードにおいては、１つのサブファミリーのみが用いられる。
【００７１】
微分は、線源パラメータ（ｄ／ｄｈ）の方向に関して実行され、２つの連続した投影間の単純投影減算として実行され、それらの直線距離で割り算される。あるいは、３点以上の有限差分式で実行できる。
【００７２】
さらにフィルタリングラインへのリビニングは、直線状データのヒルベルト再構成の一部として実行されうる。リビニングは、カーテシアン検出器座標から非カーテシアンフィルタリング座標へのリサンプリングとして実行されることを明記しておく。通常２次元内挿が要求されるが、我々はν方向に関する１次元内挿のみを用いてリサンプリングを適用し、そして我々はγ方向に関してサンプリングを同じに保つ。
【００７３】
直線状データのバックプロジェクション。直線状データのバックプロジェクションは、直接的にフィルタリングラインから実行されうる。このように、ＰＩ窓内のデータのみが用いられる（すなわち、ν＞νπ（γ）、直線状データの冗長性が考慮されるように）。他の方法は、検出器座標に逆リビニングを適用することである。この場合、直線状検出器データは、ＰＩ窓でマスクされなければならず、方程式（３）によって与えられる。  
【数３】
2010510833.tif 000007
【００７４】
図８は、直線状スキャンの範囲の一例のブロック図である。アイソセンタにおける検出器の幅をＷとする。そのとき、直線状スキャンに必要な長さは２Ｗ、すなわち、直線状スキャンの長さは、アイソセンタにおける検出器サイズの２倍、又は概ね物理検出器（physical detector）のサイズでなければならない。この場合、直線状データは、ＦＤＫによって正確に再構成されうる六角形領域（図８において影付けされた領域）を完全に覆う。すなわち、直線状データは、ＦＤＫフィルタ又はＺＴＳフィルタ［１５］又は［１６］のようなランプに基づくフィルタを含むランプフィルタに基づくアルゴリズムに応じたフィールドオブビューの厳密な再構成に十分足りている。
【００７５】
図９は、平面及び湾曲検出器のジオメトリを考慮するためのさらなる用語を有する図である。Ｗは、検出器幅の半分（例えば、Ｗ＝６４ｍｍ）、Ｈは１方向に関する直線状スキャンの長さ、そしてｈは直線の線源位置、ｈ＜Ｈである。ここで我々は、ｈ＞０の場合のみ記述し、そしてｈ＜０の場合は同様に扱いうるとする。（ν，γ）は、円柱状検出器座標である。νは垂直座標であり、γはファン角である。（ｗ，ρ）は、平面検出器座標である。ｗは、垂直座標であり、ｗ＝ν／ｃｏｓγである。ρは水平座標であり、ρ＝ｔａｎγである。図７Ａ及び図７Ｂのそれぞれの垂直軸はρに対応し、水平軸はγに対応する。
【００７６】
平面検出器。初めに、平面検出器に対する方程式が示される。そして、それら方程式は、円柱状検出器に適用されるために、平面検出器上のラインｗ＝ｍρ＋ｂが円柱状検出器上のカーブν＝ｂｃｏｓγ＋ｍｓｉｎγに対応することを示す関係を用いて変換される。
【００７７】
ＰＩ窓の投影境界は、方程式（４）により与えられる。  
【数４】
2010510833.tif 000008
【００７８】
上記方程式によれば、検出器上のコンボルーションラインは、放物線ｗＰＩ（ρ，ｈ）の接線である。放物線上の点（ｗ，ρ）の接線の傾きは、  
【数５】
2010510833.tif 000009
【００７９】
によって与えられる。
【００８０】
このように、点（ｘ０，ｙ０）を通り傾きｍを有するラインの方程式は、  
【数６】
2010510833.tif 000010
【００８１】
によって与えられる。
【００８２】
放物線ｗＰＩ（ρ，ｈ）上の与えられた点（ｗ０，ρ０）に対して、この点を通る接線の傾きはｍ＝−ｈρ０であり、ｗ０＝−ｈ（１＋ρ０２）／２である。この時、この点を通る接線の方程式は、ｗ＝−ｈ（１＋ρ０２）／２−ｈρ０（ρ−ρ０）である。この方程式を解くことによって、我々は平面検出器上の接線の方程式を得る。  
【数７】
2010510833.tif 000011
【００８３】
次に、ρ０の極限、すなわちρｍｉｎ及びρｍａｘは、２つの場合、すなわちフルスキャン及びハーフスキャンのそれぞれにおいて、フィルタリングラインが検出器上のＢＰＪエリアの全体を覆うように選択される。
【００８４】
フルスキャンの場合において、フィルタリングファミリーは、対称的である、すなわち、ρｍｉｎ＝−ρｍａｘである。従って、ρｍａｘを探せば十分である。パラメータρｍａｘは、ＦＤＫボリューム投影のピーク点を覆うフィルタリングラインによって決定される（臨界状態）。このように、点Ａの座標は、（ρ＝０,ｗ＝Ｗ−ｈ）である。これら座標を方程式（７）に代入すると、我々は、ρ０がρｍａｘに対応するということが分かる。  
【数８】
2010510833.tif 000012
【００８５】
これから、  
【数９】
2010510833.tif 000013
【００８６】
であり、ここで、  
【数１０】
2010510833.tif 000014
【００８７】
である。
【００８８】
式（９）〜（１０）は、直線状スキャンの長さＨに制限を課す。座標ρＸは、負であってはならない。これは、条件：ｈ≦２Ｗを強要する。従って、ＦＤＫボリューム再構成のために、直線状スキャンの長さは、検出器範囲の厳密に２倍、Ｈ＝２Ｗでなければならない。
【００８９】
図１０は、Ｗ＝６４ｍｍ，ｈ＝１００ｍｍ，及びρＸ＝０．２８における本実施形態の一例に係る、ｗ対ρのグラフ上のＢＰＪエリアのプロットである。
【００９０】
ハーフスキャンエリアにおいて、フィルタリングファミリーは、ＢＰＪは１つの側のみを覆うべきであるので、対称でない。我々が右からＢＰＪを覆う場合、ρｍａｘは、方程式（１０）で与えられる。
【００９１】
パラメータρｍｉｎは、ＢＰＪエリアによって決定される。ρ１＝ｔａｎγｍａｘとする。このとき、  
【数１１】
2010510833.tif 000015
【００９２】
ＰＩ窓の投影境界は、  
【数１２】
2010510833.tif 000016
【００９３】
によって与えられる。
【００９４】
コンボルーションライン。方程式（７）から我々は、  
【数１３】
2010510833.tif 000017
【００９５】
であることが分かる。（今我々は、簡単のため“ρ０”の代わりに“ρ”と書いた）、そして我々は、平面検出器上のフィルタリングカーブの方程式を得る。  
【数１４】
2010510833.tif 000018
【００９６】
円柱状検出器上のフィルタリングファミリーは、平面検出器上のファミリーと同じ指標付きパラメータを用いることを明記しておく。従って我々は、ρに対する極限を見つけるために方程式（１０）〜（１１）を用いることができる。
【００９７】
図１１Ａは、傾いたガントリを含む本発明の実施形態の等角投影図である。図１１Ｂは、円状軌道に関連する傾いた座標の等角投影図である。傾いたガントリにおいて、アルゴリズムの構造は上記の傾いていない場合と同様のままであり、画像再構成は傾いた座標系で実行され、このことは円状データの再構成がガントリの傾きがないのと同様の状態であることを意味する。従って、この節において、我々は、傾いたガントリにおける直線状データからの再構成のみを記述する。
【００９８】
本実施形態によれば、図１１Ａに示すように、ξをガントリのチルト角とし、図１２Ａにおいて、μバー（μバーは、方程式上ではμの上に―を付したものとする）は、直線状軌跡μの単位方向ベクトルとする。ここで、図１２Ａ上の影付けされた領域は、再構成ボリュームを表す（すなわち、データ外挿が必要ない最大領域）。円と直線との交点ｙ０バー（ｙ０バー、方程式上ではｙ０の上に―を付したものとする）が座標（Ｒ，０，０）を有するように、軸ｘ１，ｘ２を回転すると便利である。
【００９９】
本実施形態において、ベクトルμバーは、座標  
【数１５】
2010510833.tif 000019
【０１００】
を有し、直線状線源の座標は、  
【数１６】
2010510833.tif 000020
【０１０１】
により与えられる。
【０１０２】
方程式（１８）は、以下に示すような、平面検出器上のフィルタリングラインのための、方程式（１９）は、湾曲検出器上のフィルタリングラインのためのものである。
【数１７】
2010510833.tif 000021
【０１０３】
方程式（１８）及び（１９）においてρＬは、方程式ρＬ＝−ｃｏｔ（λ／２）に基づくフィルタリングパラメータである。ここでλは、フィルタリング面の接点への角度である。
【０１０４】
直線状スキャンの範囲は、  
【数１８】
2010510833.tif 000022
【０１０５】
によって定義される。ガントリの傾きがゼロの場合において、方程式（２０）が上記のようにｈｍａｘ＝２Ｗにまとめられることが容易にわかる。
【０１０６】
本発明の全実施形態は、便宜的に、慣習的で汎用なコンピュータ又は、当事者にとって明らかなように、本発明の教示に応じてプログラムされたマイクロプロセッサーを用いて実行される。適切なソフトウェアは、ソフトウェア業界において通常のスキルを有する者にとって明らかなように、本開示の教示に基づく通常のスキルを備えたプログラマによって容易に用意されうる。特に、コンピュータの筐体は、ＣＰＵ、メモリ（例えば、ＤＲＡＭや、ＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、ＳＲＡＭ、ＳＤＲＡＭ、フラッシュＲＡＭ）、及び他のオプションの専用論理装置（例えばＡＳＩＣＳ）又は適合性のある論理デバイス（例えばＧＡＬ、及び再プログラム可能なＦＰＧＡ）を備えるマザーボードを収容する。コンピュータは、複数の入力装置（例えばキーボードやマウス）、及びモニターを制御するためのディスプレイカードもまた含む。さらにコンピュータは、フロッピー（登録商標）ディスクドライブ；他の着脱可能な媒体装置（例えばコンパクトディスク、テープ、及び着脱可能な光磁気媒体）；及びハードディスク又は他の固定された高密度メディアドライブを含んでいても良く、適切なデバイスバス（例えばＳＣＳＩバス、拡張ＩＤＥバス、又はウルトラＤＭＡバス）を用いて接続されていてもよい。コンピュータは、さらにコンパクトディスク読取装置、コンパクトディスク読み取り／書き込みユニット、又はコンパクトディスク・ジュークボックスを含んでもよい。それらは同じデバイスバス又は他のデバイスバスに接続されてもよい。
【０１０７】
本発明に関連したコンピュータ読取り可能な媒体の例は、コンパクトディスク、ハードディスク、フロッピー（登録商標）ディスク、テープ、光磁気ディスク、ＰＲＯＭ（例えばＥＰＲＯＭ、ＥＥＰＲＯＭ、フラッシュＥＰＲＯＭ）、ＤＲＡＭ、ＳＲＡＭ、ＳＤＲＡＭ等を含む。任意の一つ又はこれらコンピュータ読取り可能な媒体の組合せに格納されて、本発明は、コンピュータの両方のハードウェアを制御するため、及びコンピュータが人間のユーザーと対話できるようにするためのソフトウェアを含む。そのようなソフトウェアは、開発ツールのようなデバイス・ドライバや、オペレーティングシステム、ユーザーアプリケーションを含んでいてもよいが、これらに制限されない。本発明のコンピュータプログラム製品は、コンピュータによって実行された時にコンピュータに本発明の方法を行なわせるコンピュータプログラムインストラクション（例えばコンピュータコード装置）を格納するあらゆるコンピュータ読み取り可能媒体を含んでいる。本発明のコンピュータコード装置は、スクリプトや、インタープリタ、ダイナミックリンクライブラリ（ＤＬＬ）、Ｊａｖａ（登録商標）クラス、完全実行可能プログラムを含むあらゆる解釈可能又は実行可能なコードメカニズムでもよいが、これらに制限されない。さらに、本発明の処理部分は、よりよい性能、信頼性、及び／又はコストのために分配（例えば（１）多重ＣＰＵ間、（２）少なくとも１つのＣＰＵと少なくとも１つの配置可能な論理装置との間で）されてもよい。例えば、輪郭（outline）又は画像は、第1のコンピュータ上で選択され、遠隔診断のための第２のコンピュータへ送信されてもよい。
【０１０８】
その発明は、当業者に容易に明らかなように、特定用途向け集積回路の準備によって又は慣習的な構成回路の適切なネットワークを相互に接続することによって実行されてもよい。
【０１０９】
本発明への画像データの源は、Ｘ線装置またはＣＴ装置のような任意の適切な画像収集装置が可能である。さらに、収集されたデータは、既にデジタル形式でなければデジタル化されてもよい。あるいは、収集され処理される画像データの源は、画像収集装置によって生成されたデータを格納するメモリでもよい。そのメモリは、局部メモリであってもよいし、又は遠隔メモリであってもよく、その場合、ＰＡＣＳ（Picture Archiving Computer System）のようなデータ通信網は、本発明に係わる処理のために画像データにアクセスするために用いられてもよい。
【０１１０】
本発明の実施形態の実行は、コーンビームアーチファクトの完全な削除を含み、シミュレートされた、ｚ方向に関して急な変化を有する背骨を含むノイズのない胴部を用いて評価される。
【０１１１】
図１３Ａは、慣習的な円状ＦＤＫ方法を用いてシミュレートされたデータの再構成画像の一例である。
【０１１２】
図１３Ｂは、上記のように、直線状データ再構成と組合された、ランプに基づく円状データ再構成を有する本発明の実施形態を用いてシミュレートされたデータの再構成画像の一例である。図１３Ａ及び図１３Ｂにおけるディスプレイウィンドウは、Ｌ／Ｗ＝０／２００である。
【０１１３】
本発明の実施形態の実行は、人体の心臓動を現実的に模した動的な動態ファントムを用いて評価された。
【０１１４】
図１４Ａ〜１４Ｆは、シミュレートされた心臓動の影響の結果の一例である。
【０１１５】
図１４Ａは、ゲート制御された円状スキャンデータのみを用いて画像が再構成された時の心臓動の影響を示す再構成画像である。
【０１１６】
図１４Ｂは、本発明の実施形態における、ゲート制御された円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
【０１１７】
図１４Ｃは、図１４Ａと図１４Ｂとの差分を示す。
【０１１８】
図１４Ｄは、ゲート制御されていない円状スキャンデータのみを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
【０１１９】
図１４Ｅは、ゲート制御されていない円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
【０１２０】
図１４Ｆは、図１４Ｄと図１４Ｅとの差分を示す。
【０１２１】
画像例は、直線状スキャンデータのゲート制御を要求することなく円状スキャンデータをゲート制御することによって、再構成が改善され、コーンビームアーチファクトが除去されたことを示す。さらに、コーンビームアーチファクトの低減は、円状スキャンデータにおける器官動によって影響されない。
【０１２２】
ゲート制御された円状データは、心臓が比較的動いていないとき、その円状データが同一の特定心位相中に獲得された（captured）画像のみを含むように制限する効果を有する。ゲート制御は、所望の心位相中のスキャンデータのみを生成する（例えば、その位相に対応する時間のみに線源で患者を放射する）、又は複数の心位相中にデータを獲得して所望の心周期に対応するデータをのみを抽出することによって実行されうる。
【０１２３】
擬人ファントムの投影データは、円状及び直線状軌跡に沿って２５６スライススキャナで収集された。図１５は、擬人ファントムの画像である。幾つかのデータセットは、直線放射線量管電流（ｍＡ）の設定と体部分すなわち頭、胸、及び腹も依存する性能を調査するために収集された。この例に対するスキャンパラメータは、表１に与えられている。この例において、最大コーン角は、１２度を越える。直線状スキャンのための管電流は、円状スキャンのそれに比してかなり低いことを明記しておく。
【０１２４】
図１６Ａは、ＦＤＫ円状スキャンデータのみを用いた再構成画像の一例を示す。図１６Ｂは、直線状再構成データの一例である。図１６Ｃは、上記のような、本発明の実施形態を用いた再構成画像の一例である。図１６Ａ〜１６Ｃの例において、表示ウィンドウは、Ｌ／Ｗ＝４０／２００である。この例は、本発明のアプローチは、目に見えるコーンビームアーチファクトを消去することをさらに例証している。
【０１２５】
さらに図１６Ｂは、直線状スキャンデータがコーンビームアーチファクト補正のみに寄与し、他のいかなる解剖学的構造を含んでいないことを例証している。直線状スキャンデータは、ｚ方向のみに関して振動数が高いことを表し、ｘｙ平面に関してとても振動数が低いことを示している。
【０１２６】
さらに、とても低い管電流（ｍＡ）に応じて、とても低い放射線量を用いて直線状スキャンを実行することが可能である。上記の本発明の実施形態の性能評価は、コーンビームアーチファクトが５０ｍＡより低い、好ましくは２０ｍＡより低い、又は約１０ｍＡの管電流を用いて低減されうるとういうことを示す。
【０１２７】
図１７Ａは、ＦＤＫ円状データのみを用いて再構成された画像の一例を示す。図１７Ｂは、５０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｃは、図１７Ｂの例において用いられた直線状データのみの例である。図１７Ｄは、３０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｅは、図１７Ｄの例において用いられた直線状データのみの例である。図１７Ｆは、５０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｇは、図１７Ｆの例において用いられた直線状データのみの例である。従って、この例は、患者への直線状スキャンの余分な線量が最小であり、それでもコーンビームアーチファクトが低減されていることを示している。
【０１２８】
さらに上記の実施形態は、寝台位置に関するデータの不正確さに敏感でありうる。例えば、図１８Ａは、ＦＤＫ円状データ再構成のみを用いて再構成された画像である。図１８Ｂは、正確な寝台位置情報を用いた、上記の本発明の実施形態により再構成された画像である。図１８Ｃは、不正確な寝台位置情報、特に、＋／−０．３ｍｍの不正確さを有する寝台位置情報を用いた、上記の本発明の実施形態により再構成された画像である。図１８Ｃの例は、不正確な寝台位置情報に起因するストリークアーチファクトを示す。
【０１２９】
本発明のさらなる実施形態は、直線状再構成又は円状再構成の１つ又は両方を生成するために用いられる再構成パラメータを調整することによる、寝台に関する不正確を補償することを含む。その再構成パラメータは、円及び直線状軌道との共通点の付近にある円状投影データと円状投影データとの間の相関に基づいて決定される（この点において両方の軌道からの投影は、一致していなければならない）。
【０１３０】
直線状スキャンの再構成ボリューム範囲は、仮想的な拡張検出器（virtual expanded detector）を用いることによってさらに低減されうる。ＦＤＫにおける円状経路のみの再構成において、ＦＯＶ、図１９Ａの例に示すように、６角形形状により制限される（例えば、１９Ａにおける６角形形状に形作られたＲＯＩ領域）。ＦＤＫ―ＦＯＶを増加させるために、仮想的な拡張検出器が用いられ得る。拡張検出器列における非計測データは、境界の検出器例における計測データを重複させることによって、又は外挿によって得られる。図１９Ｂは、仮想的な拡張検出器を用いて得られた、拡張されたＦＯＶを示す。ＦＤＫでのＦＯＶの外側にある拡張領域における再構成は、ＦＤＫでのＦＯＶの内側よりも正確さに欠けることを明記しておく。
【０１３１】
図２０の（Ａ）〜（Ｆ）は、直線状スキャンコリメーションを用いる直線状スキャンの照射線量を低減させる方法における複数の線図である。各線図において、ｈは直線状スキャン経路上のｚ座標を示し、影付けされた領域は直線状スキャンにおけるバックプロジェクションステップにおいて用いられたデータを示す。さらに、（Ｂ）〜（Ｆ）は、ｈ＞０（直線状スキャンの上部）の場合のみを示しているが、ｈ＜０の場合は対称的であり、対応する方法で扱われうる。このように、本実施形態によれば、直線状スキャンに患者によって受けられる放射線量は、影付けされた領域の外側の放射を避けるため、動的なコリメーションを用いて低減される。
【０１３２】
図２１は、動的コリメーションに対するコリメーション窓の一例を示す。この例によれば、方程式（９）、（１０）、（１３）、（１５）を用いて我々は、垂直トップ（vertical top）とボトム境界(bottom boundaries)とのために、νＴＯＰ（ｈ）とνＢＯＴ（ｈ）とをそれぞれ得る。  
【数１９】
2010510833.tif 000023
【０１３３】
さらに、方程式（２１）は、ｈ＞５．８１ｍｍの時におけるνＴＯＰを計算するためにのみ用いられる。ｈ＜５．８１ｍｍの時、νＴＯＰ＝Ｗである。
【０１３４】
図２２は、本発明の係るジオメトリのさらなる図である。本発明によれば、以下のパラメータは、システムオペレータによってスキャン条件に応じて決定される。  
【数２０】
2010510833.tif 000024
【０１３５】
Ｚｓｔａｒｔ  直線状スキャンの開始ｚ位置  
Ｚｅｎｄ  直線状スキャンの終了ｚ位置  
Ｚｃｉｒｃ  円のｚ位置  
ＴｏｔａｌＶｉｅｗｓ  直線ビューの全数  
以下のパラメータは、本発明の実施形態に応じて計算される。  
【数２１】
2010510833.tif 000025
【０１３６】
Ｎｖｉｅｗ  再構成のためのビュー数  
ＶｉｅｗＪｕｍｐ  再構成のための最初のビュー  
Δｈ  直線ビューピッチ  
ｈ  直線ビュー位置、円に相対的
さらに  
【数２２】
2010510833.tif 000026
【０１３７】
又は、  
【数２３】
2010510833.tif 000027
【０１３８】
によれば、方程式（２７）が満足される。  
【数２４】
2010510833.tif 000028
【０１３９】
そして、そのような定義は、Ｚｅｎｄ＞Ｚｓｔａｒｔ、又はＺｅｎｄ＜Ｚｓｔａｒｔの場合、独立である。さらに、  
【数２５】
2010510833.tif 000029
【０１４０】
図２３Ａは、直線状スキャンデータにおけるフィルタラインのグラフである。
【０１４１】
図２３Ｂは、直線状スキャンデータにおけるフィルタラインの数のグラフである。詳細には、この例によれば、フィルタリングラインの数は、直線状部分ＡＣ＝Ｗ−ｈ／２の長さに比例する。好ましくは、フィルタリングラインＮｆｌの数は、以下のように計算される。  
【数２６】
2010510833.tif 000030
【０１４２】
図２３Ｂは、フィルタラインＮｆｌの数のグラフである。
【０１４３】
図２４は、本発明の他の実施形態に係る直線及び円状アルゴリズムの方法のブロック図である。この場合において、収集された直線状データにおけるサンプル（チャンネル）の数はファクタＫＤＳによって低減されうる。ここで、例としてダウンサンプリング係数ＫＤＳが２、３、４、又は５に一致しうる。
【０１４４】
図２５Ａと図２５Ｂとは、ショートスキャンコンボルーションにおいて用いられるカーブの２つのファミリーの一例を示す。図２５Ｃと図２５Ｄとは、カーブの２つのファミリーのさらなる表現を示す図である。フィルタリング平面（又は湾曲）は、フィルタリング平面と検出器面（平坦又は円筒形。検出器のタイプによる）との交差によって得られる。特に、図２５Ａに示されるカーブは、図２５Ｃにおける弧Ａに接するフィルタリング平面に対応し、図２５Ｂにおけるカーブは、図２５Ｄにおける弧Ｂに接するフィルタリング平面に対応する。
【０１４５】
本発明の実施形態によれば、ショートスキャン（これは円弧の部分のみに沿ってスキャンする）は、円状スキャンの一部が用いられることに依存する、ラインの１つのファミリーのみを用いる（すなわち、図２５Ａ又は図２５Ｂの１つ）。すなわち、各検出器ピクセルに対して、１つのフィルタリングラインのみが存在する。
【０１４６】
図２６Ａは、複数のファミリー（すなわち、図２６Ａにおける全ライン）からのラインが用いられる場合における、フルスキャンコンボルーションの一例を示す図である。このファミリーは、検出器を２回覆う。すなわち、各検出器ピクセルに対して、２つのフィルタリングラインが存在する。各フィルタリングに対応する値を得るために、２つの値の差分を得る。  
【数２７】
2010510833.tif 000031
【０１４７】
図２６Ｂは、図２６Ａの２つのファミリーのラインの重なりの１次元図である。
【０１４８】
図２７Ａと図２７Ｂとは、本発明の実施形態に係る逆リビニングの一例を示す。図２７Ａは、チャンネル１０に対する逆リビニングカーブ２７０、チャンネル４４７（例えば、センターチャンネル）に対する逆リビニングカーブ２７１、及びチャンネル８８７に対する逆リビニングカーブ２７２示す。
【０１４９】
図２７Ｂは、逆リビニング処理の一例を示す。この例によれば、初めに逆リビニングカーブ２７４とセグメント座標Ｓｅｇとの交点のρＡとρＢとを得る。次に、これらρＡとρＢとにおけるＣＯＮＶ値を内挿し、差分をとる。  
【数２８】
2010510833.tif 000032
【０１５０】
図２８Ａは、本発明の実施形態に係る、拡張再構成関心ボリューム（ＶＯＩ）の一例を示す。物理検出器２８２は、領域２８４に対応するデータを受信する。しかしながら、ｚ軸に関する長さ２Ｗを有する仮想的に拡張された検出器２８６は、領域２８８に対応するデータを受信する。このように、表示ボリューム２９０内において、ＦＤＫ―ＶＯＩ２９２は物理検出器に対応し、ＦＤＫ―ＶＯＩ２９４は仮想的に拡張された検出器に対応する。図２８Ｂは、分かりやすさのため、ＦＤＫ―ＶＯＩ２９２とＦＤＫ―ＶＯＩ２９４との形状の別々の線図である。
【０１５１】
複数の画像、大きいスキャン対象物の差分部分のそれぞれは、発明的なシステムによって連続して獲得されうり、全体的な対象物の１つの画像を生成するために組合わされうる。特に、複数の円状スキャンからのデータは、１つ又はそれ以上の直線状スキャンからのデータに組合わされうる。
【０１５２】
図２９Ａ〜Ｆは、１つ又はそれ以上の直線状スキャンを有する複数の円状スキャンからのデータの結合例を示す。図２９Ａ〜Ｆのそれぞれにおいて、各分離された水平線３００は、単一の直線状スキャンに対応し、各円３０１は、単一の円状スキャン経路（これは、この経路を通過する１つ以上の円状スキャンを含みうる）。さらに、様々な直線及び円状スキャンがあらゆる適当な順番において実行されうる。
【０１５３】
それに加え、前処理された円状データは、受信された円状投影データ中の一重に計測された及び二重に計測されたレイに対して別個の重みが与えられうる。例として、完全回転以下のデータを用いた場合、（パーカー（Parker）、シルバー（Silver））［１７］、［１８］により記述されたような冗長性重みづけが用いられうる。
【０１５４】
ｇ（β，γ，α）を、２つの角（γ，α）によって与えられる方向での線源位置ａ（β）からのレイに沿う線積分を表すものとする。ここで、γはファン角を、αはコーン角を表す。角γ及びαは、物理検出器によって制限される。すなわち、−γｍ≦γ≦γｍ、及び−αｍ≦α≦αｍである。例えば、幾つかのスキャナは、−３０°＜γ＜３０°、−６°＜α＜６°を有する。ファンビームジオメトリの場合には、ｇ（β，γ）＝ｇ（β＋π＋２γ，−γ）であることを明記しておく。従って、我々が全回転分のファンビームデータを得る場合（いわゆるフルスキャン）、各レイは、２回計数される。そして我々は、ｗ（β，γ）＝１／２でデータを重みづけしうる。画像を再構成するために、全回転に亘ってスキャンする必要はない。ｇ（β，γ）＝ｇ（β＋π＋２γ，−γ）という関係から、我々は、π＋２γｍの再構成範囲のみが厳密なファンビーム再構成に十分であることが分かる。ここでγｍは、検出器が取り得る最大ファン角である。
【０１５５】
パーカー（Parker）［１７］において、不連続性ができるだけ均一に分布されるように、
最小な完全データセット（π＋２γｍ）においてデータを重みづけすることが提案されている。そして以下の重みづけ関数が提案された。  
【数２９】
2010510833.tif 000033
【０１５６】
ここから、定義されていない場合、重みはゼロであるとする。例えば、β＜０、あるいはβ＞π＋２γｍの場合、ｗＰ（β，γ）＝０である。仮想的にγｍを増加させることによって、我々は、より大きい再構成範囲を得られ、従って良好なノイズ特性が得られることも明記しておく（シルバー（Silver））［１８］。ｗＰ（β，γ）に対する方程式において、物理的な最大ファン角γｍを仮想的なΓ（Γ≧γｍ）で置き換えることにより、我々は、他の重み関数を得る。  
【数３０】
2010510833.tif 000034
【０１５７】
パーカー（Parker）の重み関数は、ＭＨＳ重みづけの特殊な事例であり、従ってｗＰ（β，γ）の代わりにｗＭＨＳ（β，γ）を用いることが好まれる。
【０１５８】
他の重み関数がノー（Noo）によって提案された［１９］。これは、任意の再構成範囲Λ＝（β０，β１）を許可する。ここでβ０及びβ１はそれぞれ、再構成角範囲の開始及び終了点である。これは、  
【数３１】
2010510833.tif 000035
【０１５９】
によって与えられる。ここで、  
【数３２】
2010510833.tif 000036
【０１６０】
である。そして関数ｃ（β）は、  
【数３３】
2010510833.tif 000037
【０１６１】
によって与えられる。ここで、Δβは、所定数（２０°）、又は全体の再構成範囲（β０−β１）の割合（但し５０％を超えない）として選ばれ得る平滑化間隔である。
【０１６２】
本発明の多数の改良及び変形が上記の技術に基づいて可能である。従って、添付された請求項の範囲内において、ここに特別に記述された以外のやり方で本発明が実施可能であるということが認識されるべきである。
【図１】
2010510833.tif 000038
【図２】
2010510833.tif 000039
【図３】
2010510833.tif 000040
【図４】
2010510833.tif 000041
【図５】
2010510833.tif 000042
【図６】
2010510833.tif 000043
【図７Ａ】
2010510833.tif 000044
【図７Ｂ】
2010510833.tif 000045
【図８】
2010510833.tif 000046
【図９】
2010510833.tif 000047
【図１０】
2010510833.tif 000048
【図１１Ａ】
2010510833.tif 000049
【図１１Ｂ】
2010510833.tif 000050
【図１２Ａ】
2010510833.tif 000051
【図１２Ｂ】
2010510833.tif 000052
【図１３Ａ】
2010510833.tif 000053
【図１３Ｂ】
2010510833.tif 000054
【図１４Ａ】
2010510833.tif 000055
【図１４Ｂ】
2010510833.tif 000056
【図１４Ｃ】
2010510833.tif 000057
【図１４Ｄ】
2010510833.tif 000058
【図１４Ｅ】
2010510833.tif 000059
【図１４Ｆ】
2010510833.tif 000060
【図１５】
2010510833.tif 000061
【図１６Ａ】
2010510833.tif 000062
【図１６Ｂ】
2010510833.tif 000063
【図１６Ｃ】
2010510833.tif 000064
【図１７Ａ】
2010510833.tif 000065
【図１７Ｂ】
2010510833.tif 000066
【図１７Ｃ】
2010510833.tif 000067
【図１７Ｄ】
2010510833.tif 000068
【図１７Ｅ】
2010510833.tif 000069
【図１７Ｆ】
2010510833.tif 000070
【図１７Ｇ】
2010510833.tif 000071
【図１８Ａ】
2010510833.tif 000072
【図１８Ｂ】
2010510833.tif 000073
【図１８Ｃ】
2010510833.tif 000074
【図１９Ａ】
2010510833.tif 000075
【図１９Ｂ】
2010510833.tif 000076
【図２０】
2010510833.tif 000077
【図２１】
2010510833.tif 000078
【図２２】
2010510833.tif 000079
【図２３Ａ】
2010510833.tif 000080
【図２３Ｂ】
2010510833.tif 000081
【図２４】
2010510833.tif 000082
【図２５Ａ】
2010510833.tif 000083
【図２５Ｂ】
2010510833.tif 000084
【図２５Ｃ】
2010510833.tif 000085
【図２５Ｄ】
2010510833.tif 000086
【図２６Ａ】
2010510833.tif 000087
【図２６Ｂ】
2010510833.tif 000088
【図２７Ａ】
2010510833.tif 000089
【図２７Ｂ】
2010510833.tif 000090
【図２８Ａ】
2010510833.tif 000091
【図２８Ｂ】
2010510833.tif 000092
【図２９Ａ】
2010510833.tif 000093
【図２９Ｂ】
2010510833.tif 000094
【図２９Ｃ】
2010510833.tif 000095
【図２９Ｄ】
2010510833.tif 000096
【図２９Ｅ】
2010510833.tif 000097
【図２９Ｆ】
2010510833.tif 000098
【国際調査報告】
2010510833.tif 000099