(57)【要約】【課題】少ない手間で画像処理に適用されるパラメーターを適切に設定する。【解決手段】医用動画像が取得される。医用動画像は、体内の構造物を撮影することにより得られ、第１のフレーム画像及び第２のフレーム画像を有する。第１のパラメーターが取得される。第１のパラメーター値を適用して第１のフレーム画像に画像処理が行われる。第２のパラメーター値が第１のパラメーター値から導出される。第２のパラメーター値は、一般的には、第１のパラメーター値と異なる。第２のパラメーター値を適用して第２のフレーム画像に画像処理が行われる。【選択図】図１
2015167719.tif 000002
【特許請求の範囲】
【請求項１】
体内の構造物を撮影することにより得られ少なくとも１個の第１のフレーム画像及び第２のフレーム画像を有する医用動画像を取得する取得部と、
少なくとも１個の第１のパラメーター値を取得するパラメーター取得部と、
前記少なくとも１個の第１のパラメーター値を適用して前記少なくとも１個の第１のフレーム画像に画像処理を行う第１の画像処理部と、
前記少なくとも１個の第１のパラメーター値と異なる第２のパラメーター値を前記少なくとも１個の第１のパラメーター値から導出するパラメーター導出部と、
前記第２のパラメーター値を適用して前記第２のフレーム画像に前記画像処理を行う第２の画像処理部と、
を備える画像処理装置。
【請求項２】
前記少なくとも１個の第１のパラメーター値及び前記第２のパラメーター値が画素値であり、
前記パラメーター導出部は、
画素値を階級とし画素数を度数とするヒストグラムを画像についてのヒストグラムであるとした場合に、前記少なくとも１個の第１のフレーム画像について少なくとも１個の第１のヒストグラムを作成し、前記第２のフレーム画像について第２のヒストグラムを作成し、前記少なくとも１個の第１のヒストグラムの特徴点と前記少なくとも１個の第１のパラメーター値との関係が前記第２のヒストグラムの特徴点及び前記第２のパラメーター値において再現されるように前記第２のパラメーター値を導出する
請求項１の画像処理装置。
【請求項３】
前記少なくとも１個の第１のフレーム画像が２個以上の第１のフレーム画像であり、
前記少なくとも１個の第１のパラメーター値が２個以上の第１のパラメーター値であり、
前記少なくとも１個の第１のヒストグラムが２個以上の第１のヒストグラムであり、
前記パラメーター取得部は、
前記２個以上の第１のフレーム画像の各々に適用される第１のパラメーター値を取得することにより前記２個以上の第１のパラメーター値を取得し、
前記パラメーター導出部は、
前記２個以上の第１のフレーム画像の各々について第１のヒストグラムを作成することにより前記２個以上の第１のヒストグラムを作成し、前記２個以上のヒストグラムを平均化することにより平均化された第１のヒストグラムを得、前記２個以上の第１のパラメーター値を平均化することにより平均化された第１のパラメーターを得、前記平均化された第１のヒストグラムの特徴点と前記平均化された第１のパラメーター値との関係が前記第２のヒストグラムの特徴点と前記第２のパラメーター値との関係において再現されるように前記第２のパラメーター値を導出する
請求項２の画像処理装置。
【請求項４】
前記少なくとも１個の第１のフレーム画像が２個以上の第１のフレーム画像であり、
前記少なくとも１個の第１のパラメーター値が２個以上の第１のパラメーター値であり、
前記パラメーター取得部は、
前記２個以上の第１のフレーム画像の各々に適用される第１のパラメーター値を取得することにより前記２個以上の第１のパラメーター値を取得し、
前記パラメーター導出部は、
前記２個以上の第１のパラメーター値を補間又は外挿することにより前記第２のパラメーター値を導出する
請求項１の画像処理装置。
【請求項５】
前記パラメーター導出部は、
前記２個以上の第１のパラメーター値からパラメーター値の周期的な変化を推定し、パラメーター値が推定した周期的な変化をすると仮定して前記第２のパラメーター値を導出する
請求項４の画像処理装置。
【請求項６】
前記構造物が肺及び心臓の両方又は片方である
請求項１から５までのいずれかの画像処理装置。
【請求項７】
前記画像処理が階調変換である
請求項１から６までのいずれかの画像処理装置。
【請求項８】
前記階調変換がオフセット補正及びゲイン補正の両方又は片方である
請求項７の画像処理装置。
【請求項９】
体内の構造物を撮影し医用動画像を生成する撮影装置と、
請求項１から８までのいずれかの画像処理装置と、
を備える画像撮影システム。
【請求項１０】
体内の構造物を撮影することにより得られ少なくとも１個の第１のフレーム画像及び第２のフレーム画像を有する医用動画像を取得する第１の工程と、
少なくとも１個の第１のパラメーター値を取得する第２の工程と、
前記少なくとも１個の第１のパラメーター値を適用して前記少なくとも１個の第１のフレーム画像に画像処理を行う第３の工程と、
前記少なくとも１個の第１のパラメーター値と異なる第２のパラメーター値を前記少なくとも１個の第１のパラメーター値から導出する第４の工程と、
前記第２のパラメーター値を適用して前記第２のフレーム画像に前記画像処理を行う第５の工程と、
をコンピューターに実行させる画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、体内の構造物を撮影することにより得られる医用動画像の処理に関する。
【背景技術】
【０００２】
胸部Ｘ線動画像の撮影においては、例えば、Ｘ線管によりＸ線が発生させられ、発生させられたＸ線が人体を透過させられ、人体を透過させられたＸ線が平面Ｘ線検出器（ＦＰＤ）により検出される。胸部Ｘ線動画像は、例えば、医師による診断を支援するために撮影される。着目する構造物を視認しやすくする、着目する構造物を抽出しやすくする等の目的で、胸部Ｘ線動画像に画像処理が行われる場合もある。
【０００３】
特許文献１に記載された技術においては、処理対象画像の画素ブロックと参照画像の画素ブロックとの類似度が求められ、類似度からフィルター係数が求められ、フィルター係数にしたがって処理対象画像と参照画像とが重み付け加算される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−１２６５３０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
胸部Ｘ線動画像は、多数のフレーム画像を有する。このため、胸部Ｘ線動画像に画像処理が行われる場合は、画像処理に適用されるパラメーターが多数のフレーム画像の各々に設定される。しかし、多数のフレーム画像の各々に個別のパラメーターを指定することは面倒な作業である。一方で、多数のフレーム画像の各々に同じパラメーターが設定された場合は、画像処理が適切に行われない。心臓、肺等の動きにより、適切なパラメーターがフレーム画像ごとに異なるためである。
【０００６】
この問題は、Ｘ線撮影以外により医用動画像が撮影される場合にも生じる。
【０００７】
本発明は、この問題を解決するためになされる。本発明の目的は、少ない手間で画像処理に適用されるパラメーターを適切に設定することである。
【課題を解決するための手段】
【０００８】
医用動画像が取得される。医用動画像は、体内の構造物を撮影することにより得られ、第１のフレーム画像及び第２のフレーム画像を有する。第１のパラメーターが取得される。第１のパラメーター値を適用して第１のフレーム画像に画像処理が行われる。第２のパラメーター値が第１のパラメーター値から導出される。第２のパラメーター値は、一般的には、第１のパラメーター値と異なる。第２のパラメーター値を適用して第２のフレーム画像に画像処理が行われる。
【発明の効果】
【０００９】
少ない手間で画像処理に適用されるパラメーターが適切に設定される。
【００１０】
これらの及びこれら以外の本発明の目的、特徴、局面及び利点は、添付図面とともに考慮されたときに下記の本発明の詳細な説明によってより明白となる。
【図面の簡単な説明】
【００１１】
【図１】画像撮影システムのブロック図である。
【図２】画像処理の内容を示す模式図である。
【図３】被補正領域の設定例を示す模式図である。
【図４】手動補正前のフレーム画像についてのヒストグラムの例を示す模式図である。
【図５】自動補正前のフレーム画像についてのヒストグラムの例を示す模式図である。
【図６】手動補正後のフレーム画像についてのヒストグラムの例を示す模式図である。
【図７】自動補正後のフレーム画像についてのヒストグラムの例を示す模式図である。
【図８】画像処理の内容を示す模式図である。
【図９】画像処理の内容を示す模式図である。
【図１０】下画素値の導出例を示す模式図である。
【図１１】コンピューターのブロック図である。
【発明を実施するための形態】
【００１２】
（１）画像撮影システム
図１のブロック図は、画像撮影システムの模式図である。
【００１３】
図１に示される画像撮影システム１０００は、撮影装置１０１０及び画像処理装置１０１１を備える。
【００１４】
撮影装置１０１０は、Ｘ線撮影により医用動画像を生成する。撮影装置１０１０は、Ｘ線管によりＸ線を発生し、発生したＸ線に人体を透過させ、人体を透過したＸ線を平面Ｘ線検出器（ＦＰＤ）により検出する。これにより、撮影装置１０１０は、人の体内の構造物を撮影し、人の体内の構造物を描写したフレーム画像を生成する。撮影装置１０１０は、２回以上の撮影を行い、２個以上のフレーム画像を有する医用動画像を生成する。医用動画像は、２次元動画像であり、デジタルデータである。撮影装置１０１０がＸ線撮影以外により医用動画像を生成してもよい。例えば、撮影装置１０１０が超音波撮影、核磁気共鳴画像法（ＭＲＩ）、コンピューター断層撮影（ＣＴ）等により医用動画像を生成してもよい。撮影装置１０１０が人以外の動物の体内の構造物を撮影してもよい。医用動画像が３次元動画像であってもよい。画像処理装置１０１１は、動きに応じてＸ線透過量が変化する構造物が撮影された医用動画像の画像処理に適する。例えば、画像処理装置１０１１は、肺及び心臓の両方又は片方が撮影された医用動画像の画像処理に適する。しかし、肺及び心臓以外の構造物を撮影した医用動画像を画像処理装置１０１１が処理してもよい。
【００１５】
（２）画像処理
図２の模式図は、画像処理の内容を示す。
【００１６】
図２に示される画像処理においては、n個のフレーム画像Img(1),・・・,Img(n)のうちの１個のフレーム画像Img(i)に手動のオフセット／ゲイン補正（手動補正）が行われ、n個以上のフレーム画像のうちの残余のn-1個のフレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々に自動のオフセット／ゲイン補正（自動補正）が行われる。
【００１７】
オフセット／ゲイン補正以外の階調変換が行われてもよい。例えば、ガンマ補正が行われてもよい。階調変換以外の画像処理が行われてもよい。例えば、シャープニング等が行われてもよい。
【００１８】
画像処理においては、手動補正に適用される下画素値L(i)から自動補正に適用される下画素値L(j)(j=1,・・・,i-1,i+1,・・・,n)が導出され、手動補正に適用される上画素値U(i)から自動補正に適用される上画素値U(j)(j=1,・・・,i-1,i+1,・・・,n)が導出される。
【００１９】
オフセット／ゲイン補正に適用されるパラメーターが下画素値及び上画素値以外であってもよい。例えば、オフセット／ゲイン補正に適用されるパラメーターがオフセット値、ゲイン値等であってもよい。一般的に言って、画像処理に適用されるパラメーターは、画像処理の結果に影響を与える処理条件を指定する。画像処理がガンマ補正である場合はガンマ補正に適用されるパラメーターはガンマ値等である。画像処理がシャープニングである場合はシャープニングに適用されるパラメーターは強度、半径、閾値等である。
【００２０】
（３）画像処理装置
画像処理装置１０１１は、撮影装置１０１０により撮影された医用動画像を取得し、取得した医用動画像に画像処理を行う。
【００２１】
画像処理装置１０１１は、取得部１０２０、手動設定部１０２１、自動設定部１０２２、パラメーター取得部１０２３、パラメーター導出部１０２４、手動補正部１０２５、自動補正部１０２６等を備える。
【００２２】
取得部１０２０は、通信を行い、撮影装置１０１０から直接的に又は撮影装置１０１０から画像サーバー等の撮影装置１０１０以外の装置を経由して医用動画像を取得する。医用動画像が記録された光ディスク等の記録媒体を読み取ることにより取得部１０２０が医用動画像を取得してもよい。
【００２３】
手動設定部１０２１は、操作者の操作を検出し、当該操作に従って手動補正が行われるフレーム画像Img(i)に被補正領域を設定する。被補正領域は、着目する構造物を描写した像を含むように設定される。例えば、横隔膜に着目する場合は、図３の模式図に示されるように、横隔膜を描写した像を含む矩形状の被補正領域１０３０がフレーム画像Img(i)に設定される。被補正領域が矩形状でなくてもよい。例えば、被補正領域が円形状、楕円形状等であってもよい。着目する構造物を描写した像が画像処理により抽出され、抽出された像を含むように被補正領域が設定されてもよい。被補正領域が設定されず、フレーム画像Img(i)の全体に手動補正が行われてもよい。
【００２４】
自動設定部１０２２は、フレーム画像Img(i)に設定される被補正領域に含まれる像に対応する像を含むように自動補正が行われるフレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々に被補正領域を設定する。フレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々に設定される被補正領域は、フレーム画像Img(i)に設定される被補正領域をトラッキングすることにより設定されてもよいし、フレーム画像Img(i)に設定される被補正領域に描画されている画像をテンプレート画像とするテンプレートマッチングにより設定されてもよい。フレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の全部について着目する構造物を描写した像がフレーム画像Img(j)に設定される被補正領域に含まれる場合は、フレーム画像Img(j)に設定される被補正領域の位置がフレーム画像Img(i)に設定される被補正領域の位置と一致させられてもよい。フレーム画像Img(i)に被補正領域が設定されない場合は、フレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々にも被補正領域が設定されず、フレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々の全体に自動補正が行われる。
【００２５】
フレーム画像Img(i)に被補正領域が設定され、その後にフレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々に被補正領域が設定される。これが望ましいのは、手動補正に適用される下画素値L(i)及び上画素値U(i)は、フレーム画像Img(i)に設定される被補正領域に着目して決定されることが多いからである。
【００２６】
パラメーター取得部１０２３は、操作者の操作を検出し、当該操作に従って手動補正に適用される下画素値L(i)及び上画素値U(i)を取得する。下画素値L(i)及び上画素値U(i)は、着目する構造物を描写した像を視認しやすくなるように決定される。着目する構造物を描写した像を画像処理で抽出しやすくなるように下画素値L(i)及び上画素値U(i)が決定されてもよい。
【００２７】
パラメーター導出部１０２４は、下画素値L(i)から自動補正に適用される下画素値L(j)(j=1,・・・,i-1,i+1,・・・,n)を導出し、上画素値U(i)から自動補正に適用される上画素値U(j)(j=1,・・・,i-1,i+1,・・・,n)を導出する。下画素値L(j)(j=1,・・・,i-1,i+1,・・・,n)の各々は、偶然の一致はありえるが、一般的には下画素値L(i)と異なる。上画素値U(j)(j=1,・・・,i-1,i+1,・・・,n)の各々は、偶然の一致はありうるが、一般的には上画素値U(i)と異なる。これにより、手動補正に適用される下画素値L(i)及び上画素値U(i)を設定するだけで下画素値L(1),・・・,L(n)及び上画素値U(1),・・・,U(n)が適切に設定される。
【００２８】
下画素値L(j)が導出される場合は、フレーム画像Img(i)についてヒストグラムが作成され、フレーム画像Img(j)についてもヒストグラムが作成される。ヒストグラムにおいては、画素値が階級であり、画素数が度数である。下画素値L(j)は、フレーム画像Img(i)についてのヒストグラムの特徴点と下画素値L(i)との関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び下画素値L(j)において再現されるように導出される。上画素値U(j)も、フレーム画像Img(i)についてのヒストグラムの特徴点と上画素値U(i)との関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び上画素値U(j)において再現されるように導出される。
【００２９】
図４の模式図及び図５の模式図の各々は、フレーム画像についてのヒストグラムの例を示す。図４に示されるヒストグラムは、手動補正前のフレーム画像Img(i)についてのものである。図５に示されるヒストグラムは、自動補正前のフレーム画像Img(j)についてのものである。
【００３０】
図４に示されるようにフレーム画像Img(i)についてのヒストグラムの特徴点である画素数が０でない最も小さい画素値に下画素値L(i)が一致する場合は、図５に示されるようにフレーム画像Img(j)についてのヒストグラムの特徴点である画素数が０でない最も小さい画素値に下画素値L(j)が一致するように下画素値L(j)が導出される。これにより、フレーム画像Img(i)についてのヒストグラムの特徴点に下画素値L(i)が一致するという関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び下画素値L(j)において再現される。また、図４に示されるようにフレーム画像Img(i)についてのヒストグラムの特徴点である画素数が極小となる画素値に上画素値U(i)が一致する場合は、図５に示されるようにフレーム画像Img(j)についてのヒストグラムの特徴点である画素数が極小になる画素値に上画素値U(j)が一致するように上画素値U(j)が導出される。これにより、フレーム画像Img(i)についてのヒストグラムの特徴点に上画素値U(i)が一致するという関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び上画素値U(j)において再現される。
【００３１】
特徴点と画素値との関係は、一致することでなくてもよい。例えば、特徴点と画素値との関係が特徴点と画素値との画素値差が基準範囲内であることであってもよい。
【００３２】
図６の模式図及び図７の模式図の各々は、フレーム画像についてのヒストグラムの例を示す。図６に示されるヒストグラムは、手動補正後のフレーム画像Img'(i)についてのものである。図７に示されるヒストグラムは、自動補正後のフレーム画像Img'(j)についてのものである。
【００３３】
手動補正部１０２５は、フレーム画像Img(i)における下画素値L(i)がフレーム画像Img'(i)における最小画素値MINとなり、フレーム画像Img(i)における上画素値U(i)がフレーム画像Img'(i)における最大画素値MAXとなるように、オフセット補正及びゲイン補正の両方又は片方を行う。
【００３４】
自動補正部１０２６は、フレーム画像Img(j)(j=1,・・・,i-1,i+1,・・・,n)の各々について、フレーム画像Img(j)における下画素値L(j)がフレーム画像Img'(j)における最小画素値MINとなり、フレーム画像Img(j)における上画素値U(j)がフレーム画像Img'(j)における最大画素値MAXとなるように、オフセット補正及びゲイン補正の両方又は片方を行う。
【００３５】
これにより、フレーム画像Img(1),・・・,Img(n)の変化を追うことや一のフレーム画像と他のフレーム画像とを対比することが容易になる。
【００３６】
（４）２個以上のフレーム画像の手動補正（その１）
図８の模式図は、画像処理の内容を示す。図８は、n個のフレーム画像Img(1),・・・,Img(n)のうちの連続する２個のフレーム画像Img(i)及びImg(i+1)の各々に手動補正が行われ、n個のフレーム画像Img(1),・・・,Img(n)のうちの残余のn-2個のフレーム画像Img(j)(j=1,・・・,i-1,i+2,・・・,n)の各々に自動補正が行われる場合を示す。
【００３７】
２個のフレーム画像Img(i)及びImg(i+1)の各々に手動補正が行われる場合は、パラメーター取得部１０２３は、フレーム画像Img(i)に行われる手動補正に適用される下画素値L(i)及び上画素値U(i)を取得し、フレーム画像Img(i+1)に行われる手動補正に適用される下画素値L(i+1)及び上画素値U(i+1)を取得する。これにより、２個の下画素値L(i)及びL(i+1)が取得され、２個の上画素値U(i)及びU(i+1)が取得される。
【００３８】
パラメーター導出部１０２４は、２個の下画素値L(i)及びL(i+1)から自動補正に適用される下画素値L(j)(j=1,・・・,i-1,i+2,・・・,n)を導出し、２個の上画素値U(i)及びU(i+1)から自動補正に適用される上画素値U(j)(j=1,・・・,i-1,i+1,・・・,n)を導出する。下画素値L(j)及び上画素値U(j)が導出される場合は、フレーム画像Img(i)及びImg(i+1)の各々についてヒストグラムが作成され２個のヒストグラムが作成され、フレーム画像Img(j)についてヒストグラムが作成される。２個の下画素値L(i)及びL(i+1)が平均化され、２個の上画素値U(i)及びU(i+1)が平均化され、２個のヒストグラムが平均化される。これにより、平均化された下画素値、平均化された上画素値及び平均化されたヒストグラムが得られる。下画素値L(j)は、平均化されたヒストグラムの特徴点と平均化された下画素値との関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び下画素値L(j)において再現されるように導出される。上画素値U(j)も、平均化されたヒストグラムの特徴点と平均化された上画素値との関係が、フレーム画像Img(j)についてのヒストグラムの特徴点及び上画素値U(j)において再現されるように導出される。
【００３９】
３個以上のフレーム画像に手動補正が行われてもよく、手動補正が行われるフレーム画像が連続していなくてもよい。
【００４０】
２個以上のフレーム画像の各々に被補正領域が設定され手動補正が行われる場合は、手動補正が行われるフレーム画像に設定された被補正領域の位置を補間又は外挿することにより自動補正が行われるフレーム画像に設定される被補正領域の位置が決定されてもよい。位置に代えて又は位置に加えて形状が補間されてもよい。
【００４１】
（５）２個以上のフレーム画像への手動補正（その２）
図９の模式図は、画像処理の内容を示す。図９は、n個のフレーム画像Img(1),・・・,Img(n)のうちの連続しない２個のフレーム画像Img(i)及びImg(i+3)の各々に手動補正が行われ、n個のフレーム画像Img(1),・・・,Img(n)のうちの残余のn-2個のフレーム画像Img(j)(j=1,・・・,i-1,i+1,i+2,i+4,・・・,n)の各々に自動補正が行われる場合を示す。
【００４２】
２個のフレーム画像Img(i)及びImg(i+3)に手動補正が行われる場合は、パラメーター取得部１０２３は、フレーム画像Img(i)に行われる手動補正に適用される下画素値L(i)及び上画素値U(i)を取得し、フレーム画像Img(i+3)に行われる手動補正に適用される下画素値L(i+3)及び上画素値U(i+3)を取得する。これにより、２個の下画素値L(i)及びL(i+3)が取得され、２個の上画素値U(i)及びU(i+3)が取得される。
【００４３】
パラメーター導出部１０２４は、２個の下画素値L(i)及びL(i+3)から自動補正に適用される下画素値L(j)(j=1,・・・,i-1,i+1,i+2,i+4,・・・,n)を導出し、２個の上画素値U(i)及びU(i+3)から自動補正に適用される上画素値U(j)(j=1,・・・,i-1,i+1,i+2,i+4,・・・,n)を導出する。フレーム画像Img(i)とフレーム画像Img(i+3)との間に撮影されたフレーム画像Img(j)(j=1+1,i+2)の各々に行われる自動補正に適用される下画素値L(j)及び上画素値U(j)が導出される場合は、２個の下画素値L(i)及びL(i+3)を補間することにより下画素値L(j)が導出され、２個の上画素値U(i)及びU(i+3)を補間することにより上画素値U(j)が導出される。フレーム画像Img(i)とフレーム画像Img(i+3)との間に撮影されていないフレーム画像Img(j)(j=1,・・・,i-1,1+4・・・,n)に行われる自動補正に適用される下画素値L(j)及び上画素値U(j)が導出される場合は、２個の下画素値L(i)及びL(i+3)を外挿することにより下画素値L(j)が導出され、２個の上画素値U(i)及びU(i+3)を外挿することにより上画素値L(j)が導出される。
【００４４】
手動補正が行われるフレーム画像の間に撮影されるフレーム画像の数が２個から増減されてもよい。手動補正が行われるフレーム画像が２個から増加させられてよい。補間は、線形補間であってもよいし、非線形補間であってもよい。外挿も、線形外挿であってもよいし、非線形外挿であってもよい。
【００４５】
図１０の模式図は、画像処理に適用される下画素値の導出例を示す。図１０は、手動補正に適用される下画素値L(1),L(3)及びL(6)から自動補正に適用される下画素値L(2),L(4),L(5),L(7),L(8)及びL(9)を下画素値の周期的な変化を利用して導出する場合を示す。
【００４６】
図１０に示される導出においては、パラメーター導出部１０２４は、下画素値L(1),L(3)及びL(6)から下画素値の周期的な変化Lを推定し、下画素値が推定した周期的な変化Lをすると仮定して下画素値L(2),L(4),L(5),L(7),L(8)及びL(9)を導出する。手動補正が行われるフレーム画像の数が増減されてもよく、自動補正が行われるフレーム画像の数が増減されてもよい。手動補正に適用される上画素値から自動補正に適用される上画素値を導出することも、同様に行うことができる。
【００４７】
補間が線形補間である場合又は外挿が線形外挿である場合は、下画素値及び上画素値が単調に増減する場合にのみ自動補正に適用される下画素値及び上画素値が適切に導出されるが、下画素値及び上画素値が不規則に変化する場合は自動補正に適用される下画素値及び上画素値が適切に導出されない。例えば、下画素値及び上画素値の変化が極値を有する場合は、自動補正に適用される下画素値及び上画素値が適切に導出されない。これに対して、周期的な変化を推定し推定した周期的な変化をすることを仮定する場合は、下画素値及び上画素値が適切に導出される可能性が高まる。
【００４８】
望ましくは、撮影されたものの周期的な変化のモデルを活用して下画素値の周期的な変化が推定される。例えば、肺が撮影されている場合は、呼吸周期と同じ周期を有する周期的な変化が推定される。
【００４９】
（６）コンピューター
図１１は、コンピューターのブロック図である。
【００５０】
図１１に示されるコンピューター１０４０は、ＣＰＵ１０５０、メモリー１０５１、ハードディスクドライブ１０５２等を備え、画像処理装置１０１１となる。ハードディスクドライブ１０５２が他の種類の補助記憶装置に置き換えられてもよい。ハードディスクドライブ１０５２には、画像処理プログラム１０６０がインストールされる。画像処理装置１０１１の機能は、ＣＰＵ１０５０が画像処理プログラム１０６０をメモリー１０５１にロードして実行することにより実現される。画像処理装置１０１１の機能の全部又は一部がソフトウェアを伴わないハードウェアにより実現されてもよい。画像処理装置１０１１の機能の全部又は一部が２台以上のコンピューターにより実現されてもよい。画像処理プログラム１０６０は、コンピューター１０４０の出荷時にハードディスクドライブ１０５２にプリインストールされてもよいし、コンピューター１０４０の出荷後にハードディスクドライブ１０５２にインストールされてもよい。画像処理プログラム１０６０のインストールは、画像処理プログラム１０６０が記録された光ディスク等の記録媒体１０７０を読み取ることにより行われてもよし、ネットワーク１０７１を経由してダウンロードすることにより行われてもよい。
【００５１】
本発明は詳細に示され記述されたが、上記の記述は全ての局面において例示であって限定的ではない。したがって、本発明の範囲からはずれることなく無数の修正及び変形が案出されうると解される。
【符号の説明】
【００５２】
１０００  画像撮影システム
Img(1),・・・,Img(n)  フレーム画像
L(1),・・・,L(n)  下画素値
U(1),・・・,U(n)  上画素値
【図１】
2015167719.tif 000003
【図２】
2015167719.tif 000004
【図３】
2015167719.tif 000005
【図４】
2015167719.tif 000006
【図５】
2015167719.tif 000007
【図６】
2015167719.tif 000008
【図７】
2015167719.tif 000009
【図８】
2015167719.tif 000010
【図９】
2015167719.tif 000011
【図１０】
2015167719.tif 000012
【図１１】
2015167719.tif 000013