(57)【要約】【課題】  放射線画像に対する画像処理条件を、ユーザによる調整結果を利用することにより、自動的に最適化する。【解決手段】  画像処理条件に応じたパラメータを記憶するパラメータ記憶手段と、前記パラメータ記憶手段で記憶されている画像処理条件に応じたパラメータで前記放射線画像に画像処理を行う画像処理手段と、前記画像処理手段によって画像処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整手段と、前記画像処理条件調整手段で調整された画像処理条件毎のパラメータ調整量を記憶するパラメータ調整量記憶手段と、前記パラメータ調整量に基づいて前記パラメータ記憶手段で記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更手段と、前記パラメータの変更を所定の条件で停止させるための判断を行うパラメータ修正許可判定手段と、を有する。【選択図】      図１
2005111249.tif 000002
【特許請求の範囲】
【請求項１】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理方法であって、
１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理ステップと、
前記画像処理ステップによって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整ステップと、
前記画像処理条件調整ステップで調整された前記画像処理条件に含まれるパラメータの調整量を記憶するパラメータ調整量記憶ステップと、
前記記憶されたパラメータの調整量に基づいて、前記パラメータ記憶ステップに記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更ステップと、
を有することを特徴とする画像処理方法。
【請求項２】
複数の画像処理条件を有するとともに適切な画像処理条件を前記画像処理ステップに供給する画像処理条件設定ステップを更に有し、
前記パラメータ調整量記憶ステップは、複数の各画像処理条件毎にパラメータの調整量を記憶するとともに、
前記パラメータ変更ステップは、前記パラメータ調整量に基づいて、各画像処理条件毎にパラメータの変更を行う、
ことを特徴とする請求項１記載の画像処理方法。
【請求項３】
前記パラメータ調整量記憶ステップは、前記画像処理条件調整ステップで調整された前記画像処理条件毎のパラメータ調整量に関して、該画像処理条件毎に複数回の画像処理についてのパラメータ調整量を記憶し、
前記パラメータ変更ステップは、前記記憶されたパラメータ調整量のうち、調整量の絶対値が大きい所定回分の調整量を除外してパラメータの変更を行う、
ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法。
【請求項４】
前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、
ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法。
【請求項５】
所定の条件を満たす場合に、前記パラメータ変更ステップにおけるパラメータの変更を停止するパラメータ修正許可判定ステップをさらに有する、
ことを特徴とする請求項１記載の画像処理方法。
【請求項６】
前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、
ことを特徴とする請求項５記載の画像処理方法。
【請求項７】
前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、
ことを特徴とする請求項５または請求項６のいずれかに記載の画像処理方法。
【請求項８】
前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、
Ｔｈｄ2＞Ｔｈｄ1、
となる関係が成立する、
ことを特徴とする請求項７記載の画像処理方法。
【請求項９】
前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について
１−Ｔｈｄ3＜Rp＜１＋Ｔｈｄ3、
の条件を満たす場合はパラメータの変更を停止する、
ことを特徴とする請求項５記載の画像処理方法。
【請求項１０】
前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について
１−Ｔｈｄ4＞Rp、または１＋Ｔｈｄ4＜Rp、
の条件を満たす場合はパラメータの変更を再開する、
ことを特徴とする請求項５または請求項９のいずれかに記載の画像処理方法。
【請求項１１】
前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、
パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、
Ｔｈｄ4＞Ｔｈｄ3、
となる関係が成立する、
ことを特徴とする請求項１０記載の画像処理方法。
【請求項１２】
前記パラメータの調整量について、前記パラメータ調整量記憶ステップに記憶するか否かを判定するパラメータ調整量記憶許可判定ステップをさらに有する、
ことを特徴とする請求項１乃至請求項１１のいずれかに記載の画像処理方法。
【請求項１３】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理方法であって、
１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理ステップと、
前記画像処理ステップによって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整ステップと、
前記画像処理条件調整ステップで調整された前記画像処理条件に含まれるパラメータの調整後のパラメータの値を記憶する修正パラメータ記憶ステップと、
前記調整後のパラメータに基づいて、前記パラメータ記憶ステップに記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更ステップとを有する、
ことを特徴とする画像処理方法。
【請求項１４】
前記パラメータ変更ステップは、離散的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回の調整後のパラメータ値の最頻値で前記変更するパラメータを決定する、
ことを特徴とする請求項１３記載の画像処理方法。
【請求項１５】
前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整後のパラメータの平均値と当該パラメータの値との差が減少するようにパラメータを変更する、
ことを特徴とする請求項１３記載の画像処理方法。
【請求項１６】
前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、
ことを特徴とする請求項１３または請求項１４のいずれかに記載の画像処理方法。
【請求項１７】
所定の条件を満たす場合に、パラメータの変更を停止するパラメータ修正許可判定ステップをさらに有する、
ことを特徴とする請求項１３記載の画像処理方法。
【請求項１８】
前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、
ことを特徴とする請求項１７記載の画像処理方法。
【請求項１９】
前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、
ことを特徴とする請求項１７または請求項１８のいずれかに記載の画像処理方法。
【請求項２０】
前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、
Ｔｈｄ2＞Ｔｈｄ1、
となる関係が成立する、
ことを特徴とする請求項１９記載の画像処理方法。
【請求項２１】
前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について、
１−Ｔｈｄ3＜Rp＜１＋Ｔｈｄ3、
の条件を満たす場合はパラメータの変更を停止する、
ことを特徴とする請求項１７記載の画像処理方法。
【請求項２２】
前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について
１−Ｔｈｄ4＞Rp、または１＋Ｔｈｄ4＜Rp、
の条件を満たす場合はパラメータの変更を再開する、
ことを特徴とする請求項１７または請求項２１のいずれかに記載の画像処理方法。
【請求項２３】
前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、
Ｔｈｄ4＞Ｔｈｄ3、
となる関係が成立する、
ことを特徴とする請求項２２に記載の画像処理方法。
【請求項２４】
前記調整後のパラメータについて、前記修正パラメータ記憶ステップに記憶するか否かを判定する修正パラメータ記憶許可判定ステップをさらに有する、
ことを特徴とする請求項１３乃至請求項２２のいずれかに記載の画像処理方法。
【請求項２５】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理装置であって、
１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理手段と、
前記画像処理手段によって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整手段と、
前記画像処理条件調整手段で調整された前記画像処理条件に含まれるパラメータの調整量を記憶するパラメータ調整量記憶手段と、
前記記憶されたパラメータの調整量に基づいて、前記パラメータ記憶手段に記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更手段と、
を有することを特徴とする画像処理装置。
【請求項２６】
複数の画像処理条件を有するとともに適切な画像処理条件を前記画像処理手段に供給する画像処理条件設定手段を更に有し、
前記パラメータ調整量記憶手段は、複数の各画像処理条件毎にパラメータの調整量を記憶するとともに、
前記パラメータ変更手段は、前記パラメータ調整量に基づいて、各画像処理条件毎にパラメータの変更を行う、
ことを特徴とする請求項２５記載の画像処理装置。
【請求項２７】
前記パラメータ調整量記憶手段は、前記画像処理条件調整手段で調整された前記画像処理条件毎のパラメータ調整量に関して、該画像処理条件毎に複数回の画像処理についてのパラメータ調整量を記憶し、
前記パラメータ変更手段は、前記記憶されたパラメータ調整量のうち、調整量の絶対値が大きい所定回分の調整量を除外してパラメータの変更を行う、
ことを特徴とする請求項２５または請求項２６のいずれかに記載の画像処理装置。
【請求項２８】
前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、
ことを特徴とする請求項２５または請求項２６のいずれかに記載の画像処理装置。
【請求項２９】
所定の条件を満たす場合に、前記パラメータ変更手段におけるパラメータの変更を停止するパラメータ修正許可判定手段をさらに有する、
ことを特徴とする請求項２５記載の画像処理装置。
【請求項３０】
前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、
ことを特徴とする請求項２９記載の画像処理装置。
【請求項３１】
前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、
ことを特徴とする請求項２９または請求項３０のいずれかに記載の画像処理装置。
【請求項３２】
前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、
Ｔｈｄ2＞Ｔｈｄ1、
となる関係が成立する、
ことを特徴とする請求項３１記載の画像処理装置。
【請求項３３】
前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について、
１−Ｔｈｄ3＜Rp＜１＋Ｔｈｄ3、
の条件を満たす場合はパラメータの変更を停止する、
ことを特徴とする請求項２９記載の画像処理装置。
【請求項３４】
前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について、
１−Ｔｈｄ4＞Rp、または１＋Ｔｈｄ4＜Rp、
の条件を満たす場合はパラメータの変更を再開する、
ことを特徴とする請求項２９または請求項３３のいずれかに記載の画像処理装置。
【請求項３５】
前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、
パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、
Ｔｈｄ4＞Ｔｈｄ3、
となる関係が成立する、
ことを特徴とする請求項３４記載の画像処理装置。
【請求項３６】
前記パラメータの調整量について、前記パラメータ調整量記憶手段に記憶するか否かを判定するパラメータ調整量記憶許可判定手段をさらに有する、
ことを特徴とする請求項２５乃至請求項３５のいずれかに記載の画像処理装置。
【請求項３７】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理装置であって、
１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理手段と、
前記画像処理手段によって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整手段と、
前記画像処理条件調整手段で調整された前記画像処理条件に含まれるパラメータの調整後のパラメータの値を記憶する修正パラメータ記憶手段と、
前記調整後のパラメータに基づいて、前記パラメータ記憶手段に記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更手段とを有する、
ことを特徴とする画像処理装置。
【請求項３８】
前記パラメータ変更手段は、離散的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回の調整後のパラメータ値の最頻値で前記変更するパラメータを決定する、
ことを特徴とする請求項３７記載の画像処理装置。
【請求項３９】
前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整後のパラメータの平均値と当該パラメータの値との差が減少するようにパラメータを変更する、
ことを特徴とする請求項３７記載の画像処理装置。
【請求項４０】
前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、
ことを特徴とする請求項３７または請求項３８のいずれかに記載の画像処理装置。
【請求項４１】
所定の条件を満たす場合に、パラメータの変更を停止するパラメータ修正許可判定手段をさらに有する、
ことを特徴とする請求項３７記載の画像処理装置。
【請求項４２】
前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、
とを特徴とする請求項４１記載の画像処理装置。
【請求項４３】
前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、
ことを特徴とする請求項４１または請求項４２のいずれかに記載の画像処理装置。
【請求項４４】
前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、
Ｔｈｄ2＞Ｔｈｄ1、となる関係が成立する、
ことを特徴とする請求項４３記載の画像処理装置。
【請求項４５】
前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について
１−Ｔｈｄ3＜Rp＜１＋Ｔｈｄ3、
の条件を満たす場合はパラメータの変更を停止する、
ことを特徴とする請求項４１記載の画像処理装置。
【請求項４６】
前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Rpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について
１−Ｔｈｄ4＞Rp、または１＋Ｔｈｄ4＜Rp、
の条件を満たす場合はパラメータの変更を再開する、
ことを特徴とする請求項４１または請求項４５のいずれかに記載の画像処理装置。
【請求項４７】
前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、
Ｔｈｄ4＞Ｔｈｄ3、
となる関係が成立する、
ことを特徴とする請求項４６に記載の画像処理装置。
【請求項４８】
前記調整後のパラメータについて、前記修正パラメータ記憶手段に記憶するか否かを判定する修正パラメータ記憶許可判定手段をさらに有する、
ことを特徴とする請求項３７乃至請求項４６のいずれかに記載の画像処理装置。
【請求項４９】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対して画像処理を実行する画像処理プログラムであって、
請求項１乃至請求項２４の画像処理を実行するために、各ステップを実行するルーチンを有することを特徴とする画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は放射線画像を処理する際の画像処理方法および画像処理装置ならびに画像処理プログラムに関し、さらに詳しくは、ユーザが加えた画像濃度やコントラスト等の画像処理の調整量に基づいて自動的にユーザの好みに合わせた画像処理条件を安定した状態に学習設定できる画像処理方法および画像処理装置ならびに画像処理プログラムに関する。
【背景技術】
【０００２】
近年、放射線画像を直接デジタル画像として撮影できる装置が開発されている。たとえば、被写体に照射された放射線量を検出し、その検出量に対応して形成される放射線画像を電気信号として得る装置としては、輝尽性蛍光体を用いたディテクタを用いる方法が特開昭55-12429号公報、特開昭63-189853号公報など、多数開示されている。
【０００３】
このような装置では、シート状の基板に輝尽性蛍光体を塗布、あるいは蒸着等によって固着したディテクタに、いったん被写体を透過した放射線を照射して輝尽性蛍光体に放射線を吸収させる。
【０００４】
その後、この輝尽性蛍光体を光または熱エネルギーで励起することにより、この輝尽性蛍光体が上記吸収によって蓄積している放射線エネルギーを蛍光として放射させ、この蛍光を光電変換して画像信号を得るようにしている。
【０００５】
一方、照射された放射線の強度に応じた電荷を光導電層に生成し、生成された電荷を二次元的に配列された複数のコンデンサに蓄積し、それら蓄積された電荷を取り出すことにより得られる放射線画像検出装置が提案されている。
【０００６】
このような放射線画像検出装置では、フラットパネルディテクタ（ＦＰＤ）と呼ばれるものを使用している。この種のＦＰＤは、特開平9-90048号公報に記載されているように、照射された放射線強度に応じた蛍光を発する蛍光体と、蛍光体から発する蛍光を直接または縮小光学系を介して受光して光電変換を行うフォトダイオードやＣＣＤのような光電変換素子の組み合わせによって実現されるものが知られている。また特開平6-342098号公報に記載されているように、照射された放射を直接電荷に変換するものも知られている。
【０００７】
これらの放射線画像検出装置においては、取得した画像に対して、診断に適した画像となるよう、階調変換処理やエッジ強調処理などの画像処理を行うことが一般的に行われている。
【０００８】
一方、これらの画像処理について、画像取得後に一々手動で行うことは煩雑である。特に撮影数の多い病院などでは、このように画像取得後に手動で画像処理の設定を行うとすると、放射線技師の負担が大きくなりすぎる。
【０００９】
そこでこれらの装置では、頭部や、胸部、腹部など撮影部位や撮影方向に応じて、最適と思われる画像処理条件を予め設定しておき、画像取得時に条件を指定して、自動的に画像処理を行う構成となっている。
【００１０】
このような自動的に画像処理を行うことが、以下の特許文献１乃至特許文献１０に記載されている。
また、複数の画像処理装置がネットワークを介して接続されたシステムにおいて、任意の画像処理装置で画像処理を行っても、他の画像処理装置で同じ画像に対し、同じ画像処理を行った場合と同一の結果が得られるよう画像処理条件を共通管理する手段を設け、その共通管理された画像処理条件を修正することで、各画像処理装置毎に修正する必要が無く、操作性を向上させたシステムが、以下の特許文献１１に開示されている。
【００１１】
一方、実際にこれら装置で取得された画像を用いて診断を行う医師の好みなどにより、最適な画像処理条件は施設毎に異なる。そのため、これらの装置を各施設に導入する際に、医師等の希望に沿って画像処理条件を調節する必要がある。
【００１２】
このような画像処理条件の調整は、撮影部位及び撮影方向毎に手動で行わなければならず、煩雑で時間を要する作業を必要とする。そこで、以上のような画像処理条件の調整を行う手段としては、以下の特許文献１２や特許文献１３に記載された記述が存在している。
【特許文献１】特開2000-83938号公報（第１頁、図１）
【特許文献２】特開2000-23952号公報（第１頁、図１）
【特許文献３】特開昭63-262141号公報（第１頁、図１）
【特許文献４】特開平4-341246号公報（第１頁、図１）
【特許文献５】特開平8-62751号公報（第１頁、図１）
【特許文献６】特開200１−29335号公報（第１頁、図１）
【特許文献７】特開平7-178076号公報（第１頁、図１）
【特許文献８】特開平8-62751号公報（第１頁、図１）
【特許文献９】特開2002-183726号公報（第１頁、図１）
【特許文献１０】特開2002-183727号公報（第１頁、図１）
【特許文献１１】特開2002-279395号公報（第１頁、図１）
【特許文献１２】特開平1１−88678号公報(第１頁、図１)
【特許文献１３】特開2000-175035号公報(第１頁、図１)
【発明の開示】
【発明が解決しようとする課題】
【００１３】
〈第１の課題〉
ところが、実際にこれら装置で取得された画像を用いて診断を行う医師の好みなどにより、最適な画像処理条件は施設毎に異なる。そのため、これらの装置を各施設に導入する際に、医師等の希望に沿って画像処理条件を調節する必要がある。
【００１４】
従来、このような画像処理条件の調整は、撮影部位及び撮影方向毎に手動で行わなければならず、煩雑で時間を要する作業を必要とした。そのため、このようなカスタマイズ作業を自動的に行える装置が望まれている。
【００１５】
〈第２の課題〉
また、そのようなカスタマイズ作業を自動的に行わせるとしても、常にカスタマイズ作業を続けているとすると、自動的なパラメータ値の変更が続くため、安定した動作が望めなくなる可能性も生じる。
【００１６】
〈第３の課題〉
また、複数の画像処理装置がネットワークを介して接続されたシステムにおいて、画像処理条件を自動的に調整する際、もし、各装置で個別に調整を行われると、同じ撮影部位、同じ撮影体位の画像に対し、画像処理装置毎に処理結果が変わってしまう。これは、安定した診断画像を提供できなくなるため、望ましくない。そこで、全ての画像処理装置で同じ画像処理条件となるよう、画像処理条件の調整が行われることが望ましい。
【００１７】
また、これら画像処理装置を複数の放射線技師などが使用する場合、経験を積んだ放射線技師は、診断に適した画像を安定的に供給できるよう、画像調整を行うことが可能である。しかし、未だ経験の浅い放射線技師の場合、画像調整の結果が安定的ではないことも考えられる。このような場合には、特定の使用者の調整結果のみを用いて、画像処理条件を調整することが望ましい。
【００１８】
本発明は以上のような課題に鑑みてなされたものであって、放射線画像に対する画像処理条件を、ユーザによる調整結果を利用することにより、自動的に最適化することが可能な画像処理方法および画像処理装置ならびに画像処理プログラムを実現することを第１の目的とする。
【００１９】
また、本発明は以上のような課題に鑑みてなされたものであって、放射線画像に対する画像処理条件を、ユーザによる調整結果を利用することにより、自動的にかつ安定した状態で最適化することが可能な画像処理方法および画像処理装置ならびに画像処理プログラムを実現することを第２の目的とする。
【００２０】
さらに、本発明は以上のような課題に鑑みてなされたものであって、放射線画像に対する画像処理条件を、ユーザによる調整結果を利用することにより、自動的に最適化する際、画像調整結果を反映する対象を特定することにより、複数の画像処理装置や複数のユーザが使用する場合でも簡便かつ安定的に画像処理条件の調整が行われる画像処理方法および画像処理装置ならびに画像処理プログラムを実現することを第３の目的とする。
【課題を解決するための手段】
【００２１】
すなわち、前記した課題を解決は、以下に列記する発明により解決される。
（１）請求項１記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理方法であって、１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理ステップと、前記画像処理ステップによって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整ステップと、前記画像処理条件調整ステップで調整された前記画像処理条件に含まれるパラメータの調整量を記憶するパラメータ調整量記憶ステップと、前記記憶されたパラメータの調整量に基づいて、前記パラメータ記憶ステップに記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更ステップと、を有することを特徴とする画像処理方法である。
【００２２】
また、請求項２５記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理装置であって、１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理手段と、前記画像処理手段によって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整手段と、前記画像処理条件調整手段で調整された前記画像処理条件に含まれるパラメータの調整量を記憶するパラメータ調整量記憶手段と、前記記憶されたパラメータの調整量に基づいて、前記パラメータ記憶手段に記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更手段と、を有することを特徴とする画像処理装置である。
【００２３】
（２）請求項２記載の発明は、複数の画像処理条件を有するとともに適切な画像処理条件を前記画像処理ステップに供給する画像処理条件設定ステップを更に有し、前記パラメータ調整量記憶ステップは、複数の各画像処理条件毎にパラメータの調整量を記憶するとともに、前記パラメータ変更ステップは、前記パラメータ調整量に基づいて、各画像処理条件毎にパラメータの変更を行う、ことを特徴とする請求項１記載の画像処理方法である。
【００２４】
また、請求項２６記載の発明は、複数の画像処理条件を有するとともに適切な画像処理条件を前記画像処理手段に供給する画像処理条件設定手段を更に有し、前記パラメータ調整量記憶手段は、複数の各画像処理条件毎にパラメータの調整量を記憶するとともに、前記パラメータ変更手段は、前記パラメータ調整量に基づいて、各画像処理条件毎にパラメータの変更を行う、ことを特徴とする請求項２５記載の画像処理装置である。
【００２５】
（３）請求項３記載の発明は、前記パラメータ調整量記憶ステップは、前記画像処理条件調整ステップで調整された前記画像処理条件毎のパラメータ調整量に関して、該画像処理条件毎に複数回の画像処理についてのパラメータ調整量を記憶し、前記パラメータ変更ステップは、前記記憶されたパラメータ調整量のうち、調整量の絶対値が大きい所定回分の調整量を除外してパラメータの変更を行う、ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法である。
【００２６】
また、請求項２７記載の発明は、前記パラメータ調整量記憶手段は、前記画像処理条件調整手段で調整された前記画像処理条件毎のパラメータ調整量に関して、該画像処理条件毎に複数回の画像処理についてのパラメータ調整量を記憶し、前記パラメータ変更手段は、前記記憶されたパラメータ調整量のうち、調整量の絶対値が大きい所定回分の調整量を除外してパラメータの変更を行う、ことを特徴とする請求項２５または請求項２６のいずれかに記載の画像処理装置である。
【００２７】
（４）請求項４記載の発明は、前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法である。
【００２８】
また、請求項２８記載の発明は、前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、ことを特徴とする請求項２５または請求項２６のいずれかに記載の画像処理装置である。
【００２９】
（５）請求項５記載の発明は、所定の条件を満たす場合に、前記パラメータ変更ステップにおけるパラメータの変更を停止するパラメータ修正許可判定ステップをさらに有する、ことを特徴とする請求項１記載の画像処理方法である。
【００３０】
また、請求項２９記載の発明は、所定の条件を満たす場合に、前記パラメータ変更手段におけるパラメータの変更を停止するパラメータ修正許可判定手段をさらに有する、ことを特徴とする請求項２５記載の画像処理装置である。
【００３１】
（６）請求項６記載の発明は、前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、ことを特徴とする請求項５記載の画像処理方法である。
【００３２】
また、請求項３０記載の発明は、前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、ことを特徴とする請求項２９記載の画像処理装置である。
【００３３】
（７）請求項７記載の発明は、前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、ことを特徴とする請求項５または請求項６のいずれかに記載の画像処理方法である。
【００３４】
また、請求項３１記載の発明は、前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、ことを特徴とする請求項２９または請求項３０のいずれかに記載の画像処理装置である。
【００３５】
（８）請求項８記載の発明は、前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、Ｔｈｄ2＞thd1、となる関係が成立する、ことを特徴とする請求項７記載の画像処理方法である。
【００３６】
また、請求項３２記載の発明は、前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、Ｔｈｄ2＞Ｔｈｄ1、となる関係が成立する、ことを特徴とする請求項３１記載の画像処理装置である。
【００３７】
（９）請求項９記載の発明は、前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について１−Ｔｈｄ3＜Ｒp＜１＋Ｔｈｄ3、の条件を満たす場合はパラメータの変更を停止する、ことを特徴とする請求項５記載の画像処理方法である。
【００３８】
また、請求項３３記載の発明は、前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について１−Ｔｈｄ3＜Ｒp＜１＋Ｔｈｄ3、の条件を満たす場合はパラメータの変更を停止する、ことを特徴とする請求項２９記載の画像処理装置である。
【００３９】
（１０）請求項１０記載の発明は、前記パラメータ修正許可判定ステップは、前記パラメータ調整量記憶ステップに記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について１−Ｔｈｄ4＞Ｒp、または１＋Ｔｈｄ4＜Ｒp、の条件を満たす場合はパラメータの変更を再開する、ことを特徴とする請求項５または請求項９のいずれかに記載の画像処理方法である。
【００４０】
また、請求項３４記載の発明は、前記パラメータ修正許可判定手段は、前記パラメータ調整量記憶手段に記憶されたパラメータ調整量から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について１−Ｔｈｄ4＞Ｒp、または１＋Ｔｈｄ4＜Ｒp、の条件を満たす場合はパラメータの変更を再開する、ことを特徴とする請求項２９または請求項３３のいずれかに記載の画像処理装置である。
【００４１】
（１１）請求項１１記載の発明は、前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、Ｔｈｄ4＞Ｔｈｄ3、となる関係が成立する、ことを特徴とする請求項１０記載の画像処理方法である。
【００４２】
また、請求項３５記載の発明は、前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、Ｔｈｄ4＞Ｔｈｄ3、となる関係が成立する、ことを特徴とする請求項３４記載の画像処理装置である。
【００４３】
（１２）請求項１２記載の発明は、前記パラメータの調整量について、前記パラメータ調整量記憶ステップに記憶するか否かを判定するパラメータ調整量記憶許可判定ステップをさらに有する、ことを特徴とする請求項１乃至請求項１１のいずれかに記載の画像処理方法である。
【００４４】
また、請求項３６記載の発明は、前記パラメータの調整量について、前記パラメータ調整量記憶手段に記憶するか否かを判定するパラメータ調整量記憶許可判定手段をさらに有する、ことを特徴とする請求項２５乃至請求項３５のいずれかに記載の画像処理装置である。
【００４５】
（１３）請求項１３記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理方法であって、１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理ステップと、前記画像処理ステップによって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整ステップと、前記画像処理条件調整ステップで調整された前記画像処理条件に含まれるパラメータの調整後のパラメータの値を記憶する修正パラメータ記憶ステップと、前記調整後のパラメータに基づいて、前記パラメータ記憶ステップに記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更ステップとを有する、ことを特徴とする画像処理方法である。
【００４６】
また、請求項３７記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理装置であって、１以上のパラメータを含む予め設定された画像処理条件で、前記放射線画像の画像処理を行う画像処理手段と、前記画像処理手段によって処理が施された処理済み画像に対し、さらに画像処理条件の調整を行う画像処理条件調整手段と、前記画像処理条件調整手段で調整された前記画像処理条件に含まれるパラメータの調整後のパラメータの値を記憶する修正パラメータ記憶手段と、前記調整後のパラメータに基づいて、前記パラメータ記憶手段に記憶されている前記画像処理条件に応じたパラメータを変更するパラメータ変更手段とを有する、ことを特徴とする画像処理装置である。
【００４７】
（１４）請求項１４記載の発明は、前記パラメータ変更ステップは、離散的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回の調整後のパラメータ値の最頻値で前記変更するパラメータを決定する、ことを特徴とする請求項１３記載の画像処理方法である。
【００４８】
また、請求項３８記載の発明は、前記パラメータ変更手段は、離散的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回の調整後のパラメータ値の最頻値で前記変更するパラメータを決定する、ことを特徴とする請求項３７記載の画像処理装置である。
【００４９】
（１５）請求項１５記載の発明は、前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整後のパラメータの平均値と当該パラメータの値との差が減少するようにパラメータを変更する、ことを特徴とする請求項１３記載の画像処理方法である。
【００５０】
また、請求項３９記載の発明は、前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整後のパラメータの平均値と当該パラメータの値との差が減少するようにパラメータを変更する、ことを特徴とする請求項３７記載の画像処理装置である。
【００５１】
（１６）請求項１６記載の発明は、前記パラメータ変更ステップは、連続的な値を有するパラメータについては、前記パラメータ記憶ステップで記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、ことを特徴とする請求項１３または請求項１４のいずれかに記載の画像処理方法である。
【００５２】
また、請求項４０記載の発明は、前記パラメータ変更手段は、連続的な値を有するパラメータについては、前記パラメータ記憶手段で記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更する、ことを特徴とする請求項３７または請求項３８のいずれかに記載の画像処理装置である。
【００５３】
（１７）請求項１７記載の発明は、所定の条件を満たす場合に、パラメータの変更を停止するパラメータ修正許可判定ステップをさらに有する、ことを特徴とする請求項１３記載の画像処理方法である。
【００５４】
また、請求項４１記載の発明は、所定の条件を満たす場合に、パラメータの変更を停止するパラメータ修正許可判定手段をさらに有する、ことを特徴とする請求項３７記載の画像処理装置である。
【００５５】
（１８）請求項１８記載の発明は、前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、とを特徴とする請求項１７記載の画像処理方法である。
【００５６】
また、請求項４２記載の発明は、前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止する、ことを特徴とする請求項４１記載の画像処理装置である。
【００５７】
（１９）請求項１９記載の発明は、前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、ことを特徴とする請求項１７または請求項１８のいずれかに記載の画像処理方法である。
【００５８】
また、請求項４３記載の発明は、前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開する、ことを特徴とする請求項４１または請求項４２のいずれかに記載の画像処理装置である。
【００５９】
（２０）請求項２０記載の発明は、前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、Ｔｈｄ2＞thd1、となる関係が成立する、ことを特徴とする請求項１９記載の画像処理方法である。
【００６０】
また、請求項４４記載の発明は、前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ1と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ2との間に、Ｔｈｄ2＞Ｔｈｄ1、となる関係が成立する、ことを特徴とする請求項４３記載の画像処理装置である。
【００６１】
（２１）請求項２１記載の発明は、前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について１−Ｔｈｄ3＜Ｒp＜１＋Ｔｈｄ3、の条件を満たす場合はパラメータの変更を停止する、ことを特徴とする請求項１７記載の画像処理方法である。
【００６２】
また、請求項４５記載の発明は、前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について１−Ｔｈｄ3＜Ｒp＜１＋Ｔｈｄ3、の条件を満たす場合はパラメータの変更を停止する、ことを特徴とする請求項４１記載の画像処理装置である。
【００６３】
（２２）請求項２２記載の発明は、前記パラメータ修正許可判定ステップは、前記修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について１−Ｔｈｄ4＞Ｒp、または１＋Ｔｈｄ4＜Ｒp、の条件を満たす場合はパラメータの変更を再開する、ことを特徴とする請求項１７または請求項２１のいずれかに記載の画像処理方法である。
【００６４】
また、請求項４６記載の発明は、前記パラメータ修正許可判定手段は、前記修正パラメータ記憶手段に記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について１−Ｔｈｄ4＞Ｒp、または１＋Ｔｈｄ4＜Ｒp、の条件を満たす場合はパラメータの変更を再開する、ことを特徴とする請求項４１または請求項４５のいずれかに記載の画像処理装置である。
【００６５】
（２３）請求項２３記載の発明は、前記パラメータ修正許可判定ステップにおいて、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、Ｔｈｄ4＞Ｔｈｄ3、となる関係が成立する、ことを特徴とする請求項２２に記載の画像処理方法である。
【００６６】
また、請求項４７記載の発明は、前記パラメータ修正許可判定手段において、パラメータの変更を停止する前記所定の閾値Ｔｈｄ3と、パラメータの変更を再開する前記所定の閾値Ｔｈｄ4との間に、Ｔｈｄ4＞Ｔｈｄ3、となる関係が成立する、ことを特徴とする請求項４６に記載の画像処理装置である。
【００６７】
（２４）請求項２４記載の発明は、前記調整後のパラメータについて、前記修正パラメータ記憶ステップに記憶するか否かを判定する修正パラメータ記憶許可判定ステップをさらに有する、ことを特徴とする請求項１３乃至請求項２２のいずれかに記載の画像処理方法である。
【００６８】
また、請求項４８記載の発明は、前記調整後のパラメータについて、前記修正パラメータ記憶手段に記憶するか否かを判定する修正パラメータ記憶許可判定手段をさらに有する、ことを特徴とする請求項３７乃至請求項４６のいずれかに記載の画像処理装置である。
【００６９】
（２５）なお、請求項４９記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対して画像処理を実行する画像処理プログラムであって、上記（１）〜（２４）の画像処理の各ステップを実行するルーチンを有することを特徴とする画像処理プログラムである。
【発明の効果】
【００７０】
以上、説明したように、本発明によれば、以下のような効果が得られる。
（１）請求項１，２５記載の発明では、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し診断に適した画像を得るための画像処理を行う際に、設定された画像処理条件毎に該画像処理条件に応じたパラメータを前もって記憶しておき、記憶されている画像処理条件に応じたパラメータで放射線画像に画像処理を行うようにしており、画像処理が施された処理済み画像に対しさらに画像処理条件の調整が行われると、調整されたパラメータ調整量を記憶し、記憶されているパラメータ調整量によってパラメータを変更する。
【００７１】
この結果、ユーザが調整した調整量によってパラメータが徐々に変更されていくため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定できる。
【００７２】
（２）請求項２，２６記載の発明では、複数の各画像処理条件毎にパラメータ調整量を記憶しており、画像処理条件毎に記憶された調整量に基づいて各画像処理条件毎にパラメータの変更を行うようにしている。この結果、撮影部位または撮影方向に基づいて設定される複数の画像処理条件を各々最適化することができる。
【００７３】
（３）請求項３，２７記載の発明では、調整された画像処理条件毎のパラメータ調整量に関して該画像処理毎に複数回の画像処理についてのパラメータ調整量を記憶し、パラメータ調整量のうち絶対値が大きい所定回分のパラメータ調整量を除外してパラメータの変更を行うようにしている。この結果、誤った処理を行った場合の画像処理条件の調整量を用いて、最適でない画像処理条件の設定を行うといった危険性を減少させることができる。
【００７４】
（４）請求項４，２８記載の発明では、連続的な値を有するパラメータについては、記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更している。この結果、画像処理条件を定めるパラメータが連続的な場合でも、最適な画像処理条件を決定できる。
【００７５】
（５）請求項５，２９記載の発明では、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し診断に適した画像を得るための画像処理を行う際に、設定された画像処理条件毎に該画像処理条件に応じたパラメータを前もって記憶しておき、記憶されている画像処理条件に応じたパラメータで放射線画像に画像処理を行うようにしており、画像処理が施された処理済み画像に対しさらに画像処理条件の調整が行われると、調整されたパラメータ調整量を記憶し、記憶されているパラメータ調整量によってパラメータを変更する。そして、このパラメータの変更について、所定の条件を満たした場合には停止させるための判断（パラメータ修正許可判定）を行うようにしている。
【００７６】
したがって、ユーザが調整した調整量に応じてパラメータが徐々に変更されていくため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定できる（変更許可状態＝学習機能ＯＮ）。この際に、ユーザが調整した調整量が小さくなり所定の条件を満たす場合には、最適な画像処理条件に近づいたと考えられるため、以上のパラメータの変更を停止させる（変更不許可状態＝学習機能ＯＦＦ）。
【００７７】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の差が所定量未満になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。
【００７８】
（６）請求項６、３０記載の発明では、以上の（５）において、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合を所定の条件を満たしたとして、パラメータの変更を停止させる。
【００７９】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更していき、かつ、調整が所定の条件で収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。
【００８０】
（７）請求項７，３１記載の発明では、以上の（６）において所定の条件を満たしたとして、パラメータの変更を停止させた後、さらに、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合を前記所定の条件の解除としてパラメータの変更を再開させる。
【００８１】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【００８２】
（８）請求項８，３２記載の発明は、以上の（７）においてパラメータ変更について自動停止と自動再開とを行う場合に、再開のための閾値Ｔｈｄ2を停止のための閾値Ｔｈｄ1より大きく設定している。
【００８３】
この結果、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【００８４】
（９）請求項９，３３記載の発明は、以上の（５）において、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との比Ｒｐを求め、該比Ｒｐが、１−Ｔｈｄ3＜Ｒｐ＜１＋Ｔｈｄ3の場合を所定の条件を満たしたとして、パラメータの変更を停止させる。
【００８５】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の比が所定の範囲になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【００８６】
（１０）請求項１０，３４記載の発明は、以上の（９）において所定の条件を満たしたとして、パラメータの変更を停止させた後、さらに、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との比Ｒｐを求め、該比Ｒｐが、１−Ｔｈｄ4＞Ｒｐ，または、１＋Ｔｈｄ4＜Ｒｐ，のいずれかの場合に、前記所定の条件の解除としてパラメータの変更を再開させる。
【００８７】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【００８８】
（１１）請求項１１，３５記載の発明は、以上の（１０）においてパラメータ変更について自動停止と自動再開とを行う場合に、再開のための閾値Ｔｈｄ4を停止のための閾値Ｔｈｄ3より大きく設定している。
【００８９】
この結果、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【００９０】
（１２）請求項１２，３６記載の発明では、（１）〜（１１）でのパラメータの調整量について、パラメータ調整量記憶をするか否かを判定するようにしている。すなわち、パラメータの学習に使用するデータを制限する手段を設けることにより、本来学習に用いるべきでないデータを学習から除外できるので、効率よくパラメータの学習を行うことができる。この結果、自動的に最適化する際、画像調整結果を反映する対象を特定することにより、複数の画像処理装置や複数のユーザが使用する場合でも簡便かつ安定的に画像処理条件の調整が行われるようになる。
【００９１】
（１３）請求項１３，３７記載の発明では、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し診断に適した画像を得るための画像処理を行う際に、設定された画像処理条件毎に該画像処理条件に応じたパラメータを前もって記憶しておき、記憶されている画像処理条件に応じたパラメータで放射線画像に画像処理を行うようにしており、画像処理が施された処理済み画像に対しさらに画像処理条件の調整が行われると、調整後のパラメータを記憶し、この調整後のパラメータによってパラメータを変更する。
【００９２】
この結果、ユーザが調整した調整後パラメータによってそれまでのパラメータが徐々に変更されていくため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定できる。
【００９３】
（１４）請求項１４，３８記載の発明では、離散的な値を有するパラメータについては、記憶された複数回の調整後のパラメータ値の最頻値で決定するようにしている。この結果、画像処理条件を定めるパラメータが離散的な場合でも、最適な画像処理条件を決定できる。
【００９４】
（１５）請求項１５，３９記載の発明では、連続的な値を有するパラメータについては、記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更している。この結果、画像処理条件を定めるパラメータが連続的な場合でも、最適な画像処理条件を決定できる。
【００９５】
（１６）請求項１６，４０記載の発明では、連続的な値を有するパラメータについては、記憶された複数回のパラメータ調整量の平均値が減少するようにパラメータを変更している。この結果、画像処理条件を定めるパラメータが連続的な場合でも、最適な画像処理条件を決定できる。
【００９６】
（１７）請求項１７，４１記載の発明では、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し診断に適した画像を得るための画像処理を行う際に、設定された画像処理条件毎に該画像処理条件に応じたパラメータを前もって記憶しておき、記憶されている画像処理条件に応じたパラメータで放射線画像に画像処理を行うようにしており、画像処理が施された処理済み画像に対しさらに画像処理条件の調整が行われると、調整後のパラメータを記憶し、この調整後のパラメータによってパラメータを変更する。そして、このパラメータの変更について、所定の条件を満たした場合には停止させるための判断（パラメータ修正許可判定）を行うようにしている。
【００９７】
したがって、ユーザが調整した調整量に応じてパラメータが徐々に変更されていくため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定できる（変更許可状態＝学習機能ＯＮ）。この際に、ユーザが調整した調整量が小さくなり所定の条件を満たす場合には、最適な画像処理条件に近づいたと考えられるため、以上のパラメータの変更を停止させる（変更不許可状態＝学習機能ＯＦＦ）。
【００９８】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の差が所定量未満になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。
【００９９】
（１８）請求項１８，４２記載の発明では、以上の（１７）において、修正パラメータ記憶ステップに記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ1未満の場合はパラメータの変更を停止させる。
【０１００】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更していき、かつ、調整が所定の条件で収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。
【０１０１】
（１９）請求項１９，４３記載の発明では、以上の（１７）（１８）において、修正パラメータ記憶で記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の差を求め、該差の絶対値が所定の閾値Ｔｈｄ2以上の場合はパラメータの変更を再開させる。
【０１０２】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【０１０３】
（２０）請求項２０，４４記載の発明は、以上の（１９）において、パラメータ変更について自動停止と自動再開とを行う場合に、再開のための閾値Ｔｈｄ2を停止のための閾値Ｔｈｄ1より大きく設定している。
【０１０４】
この結果、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【０１０５】
（２１）請求項２１，４５記載の発明は、以上の（１７）において、修正パラメータ記憶で記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ3について１−Ｔｈｄ3＜Ｒp＜１＋Ｔｈｄ3、の条件を満たす場合はパラメータの変更を停止させる。
【０１０６】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の比が所定の範囲になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１０７】
（２２）請求項２２，４６記載の発明は、以上の（１７）、（２１）において、修正パラメータ記憶で記憶された調整後のパラメータの値から計算した変更後のパラメータ値と、当該時点でのパラメータ値の比Ｒpを求め、該比の絶対値が所定の閾値Ｔｈｄ4について１−Ｔｈｄ4＞Ｒp、または１＋Ｔｈｄ4＜Ｒp、の条件を満たす場合はパラメータの変更を再開させる。
【０１０８】
この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１０９】
（２３）請求項２３，４７記載の発明は、以上の（２２）において、パラメータ変更について自動停止と自動再開とを行う場合に、再開のための閾値Ｔｈｄ4を停止のための閾値Ｔｈｄ3より大きく設定している。
【０１１０】
この結果、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１１１】
（２４）請求項２４，４８記載の発明では、（１３）〜（２３）でのパラメータの調整量について、パラメータ調整量記憶をするか否かを判定するようにしている。すなわち、パラメータの学習に使用するデータを制限する手段を設けることにより、本来学習に用いるべきでないデータを学習から除外できるので、効率よくパラメータの学習を行うことができる。この結果、自動的に最適化する際、画像調整結果を反映する対象を特定することにより、複数の画像処理装置や複数のユーザが使用する場合でも簡便かつ安定的に画像処理条件の調整が行われるようになる。
【発明を実施するための最良の形態】
【０１１２】
以下、図面を参照して本発明の実施をするための最良の形態を詳細に説明する。
本発明の実施をするための最良の形態の画像処理方法および画像処理装置ならびに画像処理プログラムの好適な実施の形態について説明する。なお、これにより本発明が限定されるものではない。
【０１１３】
なお、本実施例の形態例の各手段は、ハードウェアやファームウェア、またはソフトウェアで構成することが可能である。このため、各ステップ，各手段，各ルーチンの処理手順に沿った機能ブロック図としての図１を示す。
【０１１４】
また、この機能ブロック図である図１は、画像処理方法および画像処理装置ならびに画像処理プログラムの各実施形態を理解するためのフローチャートとしても用いることができる。
【０１１５】
以下、本実施形態の構成および動作（処理）について、（１）〜（９）に分類しつつ詳細に説明する。なお、ここでは、図５に示すフローチャートをも参照して詳細な説明を行う。
【０１１６】
〈全体構成および処理の流れ〉
全体構成：
放射線画像入力手段１０、画像処理条件設定手段２０、パラメータ記憶手段３０、画像処理条件調整手段４０、パラメータ調整量記憶手段５０、パラメータ変更手段６０、パラメータ修正許可判定手段７０、制御手段８０、および、パラメータ調整量記憶許可判定手段９０が、図１に示すように構成される。
【０１１７】
なお、以上の図１における各部（各手段）は、画像処理装置を構成する構成要素であるが、画像処理方法の各ステップ、画像処理プログラムの各ルーチンを構成するものでもある。また、以上の画像処理装置は、ＣＰＵやメモリなどと処理プログラムを組み合わせて構成することも可能であるが、プログラマブルなゲートアレイなどを用いて構成することも可能である。
【０１１８】
処理の流れ：
・放射線画像入力手段１０では、放射線撮影によって得られた画像データを取得する。取得された画像データは、制御手段８０へ送信される。
・制御手段８０では、取得された画像データから画像処理用の縮小画像、表示用の縮小画像を作成する。そして画像処理用の縮小画像を画像処理条件設定手段２０へ、表示用の縮小画像を画像処理条件調整手段４０へ送る。また、制御手段８０は、撮影部位または撮影方向等の情報をユーザインターフェースなどから取得し、これらの情報も画像処理条件設定手段２０へ送る。また、制御手段８０は画像処理手段を構成しており、後述する画像処理条件調整手段４０で決定された画像処理条件を用いて、画像データに対する画像処理を施す。画像処理が施された処理済み画像は、ネットワーク等を介して出力され、フィルム出力されたり、診断用モニタ表示される。一方、制御手段８０は、画像処理条件調整手段４０で決定された画像処理条件、および画像処理条件設定手段２０で設定された画像処理条件を、パラメータ調整量記憶許可判定手段９０へ送信する。
・画像処理条件設定手段２０では、制御手段８０より得た撮影部位または撮影方向等の情報より、必要な画像処理条件を特定して、パラメータ記憶手段３０から呼び出す。そして当該画像処理条件及び上記画像処理用縮小画像に基づいて、画像データ毎に最適な階調変換曲線等の固有の画像処理条件を決定する。決定した画像処理条件は、画像処理条件調整手段４０へ送信される。
・画像処理条件調整手段４０では、前記表示用の縮小画像に対し、前記固有の画像処理条件を用いて画像処理を施し、当該処理済み画像を表示する。その処理済み画像を確認しながら、ユーザは画像の濃度やコントラスト、エッジ強調度など、画像処理に関するパラメータを調整できる。ユーザがパラメータ調整を行うと、デフォルトの画像処理条件と、変更後の画像に対応する画像処理条件との比較を行い、パラメータの調整量及び変更後のパラメータ値を決定する。決定したパラメータは、制御手段８０へ送信される。
・パラメータ調整量記憶手段５０では、前記で求めたパラメータの調整量及び変更後のパラメータ値を、各画像処理条件毎に保存する。
・パラメータ変更手段６０では、各画像処理条件毎に、所定回数の画像処理が行われる度に画像処理条件を定めるパラメータを再計算する。計算したパラメータは、画像処理パラメータ記憶手段３０に送られ、次回以降当該パラメータを用いて画像処理が行われる。
・パラメータ修正許可判定手段７０では、各画像処理条件毎に、所定回数の画像処理が行われる度に画像処理条件を定めるパラメータの変更を継続する（許可）か、停止する（不許可）か、または既に停止している場合に再開するか、について修正許可判定を行い決定する。
・パラメータ調整量記憶許可判定手段９０では、制御手段８０からの画像処理条件（パラメータ）について、パラメータの再計算に用いるパラメータか否かを内部情報によって判別し、パラメータの再計算に用いるパラメータと判断した場合は、当該パラメータをパラメータ調整量記憶手段５０に送信する。
【０１１９】
〈各手段、各ステップ、各ルーチンの詳細内容〉
（１）放射線画像入力：
放射線画像入力ステップに用いられる放射線画像入力手段１０では、被写体を透過した放射線を検知して、画像信号として取得する。具体的な構成例としては、輝尽性蛍光体プレートを用いたものとして、特開平１１-１４２９９８号公報や特開２００２-１５６７１６号公報に記載されたものがある。また、フラットパネルディテクタ（FPD）を入力装置として用いるものには、特開平６-３４２０９８号公報に記載された、検出したX線を直接電荷に変換し、画像信号として取得するものや、特開平９-９００４８号公報に記載された、検出したX線を一旦光に変換した後、その光を受光して電荷に変換する、間接方式のものがある。
【０１２０】
なお、取得される画像データは、１画素当たりの大きさが実寸で０.０４mm〜０.２mm程度である。また本実施例では階調数を１２bit、即ち０〜４０９５の値をとるものとして扱う。さらに、これらの装置から取得した画像データを、ネットワークを介して取得するものもある。取得された画像データは、制御手段８０へ送信される。
【０１２１】
（２）制御：
各種制御を行う制御手段８０では、まず取得された画像データから画像処理用の縮小画像、表示用の縮小画像を作成する。これは、後の処理を高速に行い、かつ演算に必要なハードウェアのリソースを減らすためである。
【０１２２】
なお縮小画像では、１画素当たり被写体の実寸で０.５mm〜５mm程度に相当することが望ましい。０.５mmより小さいと、縮小画像に含まれる画素数が多くなりすぎ、高速化等の目的を達成することができない。また表示画素数の多いモニタを使用しなければ、画面上に画像全体を一度に表示することが困難となる。また５mmより大きくなると、被写体の構造がよく分からなくなり、最適な画像処理条件を求めることが困難になる。本実施例では、画像処理用、表示用とも１画素当たり実寸で１mmの大きさとする。
【０１２３】
次に、撮影部位または撮影方向等の情報をユーザインターフェースなどから取得する。これらの情報は、ユーザが撮影部位等を特定することで行う。例えば、図２に示すように、表示部とタッチパネルとを兼ね備えた当該画像処理装置のユーザインタフェース４０１で、図２（ａ）の『腹部』と表示されたボタン（ハッチング部分）を押すことや、図２（ｂ）で一度『胸部』というボタンを押した後、撮影方向を選択するポップアップ画面４０２から、『正面』というボタンを押すこと等により、各ボタンに対応付けられた画像処理条件と１対１に対応する被写体番号bnが取得される。
【０１２４】
次に、画像処理条件設定手段２０により、画像処理条件が決定されると、制御手段８０は、当該画像処理条件設定手段２０から、当該画像処理条件に相当するパラメータを取得し、制御手段８０内に準備されたキャッシュメモリのような一時保持手段に保持する。
【０１２５】
また、制御手段８０は、表示用縮小画像に対して、当該画像処理条件を用いて画像処理を施す。そして、画像処理が施された表示用縮小画像は、制御手段８０から、画像処理条件調整手段４０へ送信される。
【０１２６】
その後、画像処理条件調整手段４０においてパラメータに修正が加えられる度、制御手段８０は、画像処理条件設定手段２０より修正後の画像処理条件に相当するパラメータを受け取って、再度表示用縮小画像に対して当該修正後の画像処理条件に相当するパラメータを適用して画像処理を実行する（図５Ｓ１）。なお、ここで、制御手段８０で保持する前記画像処理条件については、後述する。
【０１２７】
画像処理条件の調整が終了し、画像処理条件調整手段４０より終了通知を受け取ると、制御手段８０は、その時点で保持している画像処理条件を用いて、原画像データに対する画像処理を施して出力する。また、制御手段８０は、保持しているパラメータ一式を、パラメータ調整量記憶許可判定手段９０を介して、パラメータ調整量記憶手段５０へ送信する。
【０１２８】
（３）画像処理条件設定：
画像処理条件の設定を行う画像処理条件設定手段２０では、制御手段８０から得た被写体番号bnに基づき、特定の画像処理用パラメータセットを画像処理パラメータ記憶手段３０から読み出す。
【０１２９】
当該パラメータセットの呼び出しは、次のように行われる。画像処理条件設定手段２０では、被写体番号bnと、対応する画像処理条件に関連する一連のパラメータが記憶されている画像処理パラメータ記憶手段３０上のメモリアドレスとを対応付けるテーブルを有する。そして、被写体番号bnが取得されると前記テーブルを参照して画像処理パラメータ記憶手段３０から一連のパラメータセットが読み出される。
【０１３０】
このパラメータには、以下に示すようなものが含まれる。
3-a）画像中から照射野領域を検出するためのアルゴリズムで使用するパラメータ、
3-b）照射野内に撮影されている被写体から、特定の構造を見つけるための複数種類のアルゴリズムを選択するためのパラメータ、
3-c）当該特定構造を見つけるアルゴリズムで使用するパラメータ、
3-d）当該特定構造部について、画像処理後に所定の濃度となるよう設定する基準信号値パラメータ、
である。
【０１３１】
具体的には、以下の通りである。
3-d-1)関心領域内の信号値の累積ヒストグラムにおける累積度数の比率h1,h2、
3-d-2)h1に対応する信号値が出力画像でとる信号値D1、
3-d-3)h2に対応する信号値が出力画像でとる信号値D2、
3-d-4)コントラスト調整係数k、
3-e）基準となる階調変換曲線stdLUT、
具体的には、入力信号値をs[s＝０,１,..,４０９５]としたときの出力信号値を要素値とする１次元配列stdLUT[s]である。
3-f）当該画像に対するエッジの強調度を決めるパラメータ、
具体的には、
・強調タイプtyp、
・各周波数帯毎の強調係数β0,β1,..,βn、
が考えられる。
【０１３２】
以上のような各種パラメータが画像処理パラメータ記憶手段３０から読み出されると、これらに基づいて画像処理条件設定手段２０において、画像処理条件の設定が以下のように行われる。
【０１３３】
なお、具体的手順は以下の通りである。
3-f-i）縮小画像に対して、分割パターン認識処理を実行し、特定の照射野領域を含む領域が抽出される。次に当該領域に対して照射野認識処理等が実行され、照射野領域だけを示す情報が得られる。なお、分割パターン認識処理としては、特開２０００-８３９３８号公報等に記載された方法があり、照射野認識処理としては、特開２０００-２３９５２号公報に記載された方法等がある。
【０１３４】
3-f-ii）前記照射野認識処理で得られた、照射野領域情報に基づき、画像診断上注目すべき領域に関心領域（ＲＯＩ）を設定する。ＲＯＩの設定は、被写体の解剖学的構造に基づき、縮小画像を解析することで設定する。ＲＯＩの設定方法としては、特開昭６３-２６２１４１号公報、特開平４-３４１２４６号公報、特開平８-６２７５１号公報に記載された方法や、特開２００１-２９３３５号公報に記載された方法などがある。
【０１３５】
3-f-iii）ＲＯＩが設定されると、ＲＯＩ内部の信号値に関する統計的情報に基づいて、基準信号値が決定される。統計的情報とは、ＲＯＩに含まれる画素の平均値や最頻値、信号値の累積ヒストグラムの所定割合に相当する信号値などがある。基準信号値が決定すると、予め準備された階調変換曲線を当該信号値に基づいてカスタマイズし、当該入力画像に適した階調変換曲線を得る。
【０１３６】
このような階調変換方法としては、特開平７−１７８０７６号公報、特開平８−６２７５１号公報に記載された方法等を利用することができる。
例えば、具体的には以下の手順で階調変換曲線LUT[]を計算する。
【０１３７】
まず、関心領域内の信号値を元に、累積ヒストグラムを算出する。
そして、信号値の低い方から累積度数を求め、関心領域内の全画素数Nroiに対する所定割合h1となる信号値std1を求める。
【０１３８】
次に、信号値の低い方から累積度数を求め、関心領域内の全画素数Nroiに対する所定割合h2となる信号値std2を求める。
一方、基準となる階調変換曲線stdLUT[]から、std1がとるべき出力信号値D1、及びstd2がとるべき出力信号値D2、に対応する入力信号値s'1、s'2を求める。
【０１３９】
そして以下の（１）式と（２）とにより、変換係数aおよび変換係数bを求める。
a＝（std2−std1）／（s'2−s'1）  …（１），
b＝（std1−s'1）                …（２），
そして、画像毎に階調変換曲線を最適化するために、これらa,b及びコントラスト調整係数kを用いて、s＝０〜４０９５に対して階調変換曲線の入力信号値を規格化する。その規格化後の入力信号値s'を以下の（３）式により計算する。
【０１４０】
s'＝k・a・s＋b    …（３），
そして上記stdLUT[]の要素番号値sをs'に置き換えた階調変換曲線LUT[s']を作成する。ここで、s'＜０ 及びs'＞４０９５（＝２12−１）となる場合は無視する（入力信号値s'は０〜４０９５の値しか取り得ないため）とともに、s'が１ずつインクリメントされるよう、補間処理を行う。
【０１４１】
さらに、h2を指定しない場合（例えばh2＜０ のときは参照しないというルールとする）は、上記と同様にstd1を求めた後、
b＝（std1−s'1）    …（４），
を計算する。そして、
s'＝k・s＋b    …（５），
として、stdLUT[s]の要素番号値sをs'に置き換えた階調変換曲線LUT[s']を作ることができる。
【０１４２】
同様に、h1を指定しない場合、std2を求めた後、
b＝（std2−s'2）    …（６），
を計算する。そして、
s'＝k・s＋b    …（７），
として、stdLUT[s]の要素番号値sをs'に置き換えた階調変換曲線LUT[s']を作ることができる。
【０１４３】
3-f-iv）また、上記階調変換処理とは別に、エッジを強調して細部構造を見易くするような強調処理が行われる。強調処理としては、特開２００２−１８３７２６号公報、特開２００２−１８３７２７号公報に記載されたピラミッドアルゴリズムによる方法がある。
【０１４４】
ピラミッドアルゴリズムでは、例えばラプラシアンフィルタを掛けて作成した非鮮鋭画像及び、原画像との差分画像の作成を繰り返すことで、周波数帯域の異なる複数の画像を得ることができ、各周波数帯域の画像毎に変換処理を行った後に加算処理を行うことで、きめ細かな強調処理を行うことができる。
【０１４５】
この変換処理としては、各周波数帯域を表すi番目の差分画像毎に、si'＝βj・si のように、各差分画像の画素信号値siを変換して画素信号値si'を計算する方法がある。そして、このβjの組合せを変えることにより、周波数帯域毎に強調度を変えることができる。
【０１４６】
得られた階調変換曲線、前記信号値D１およびD２、コントラスト調整係数k、エッジ強調タイプTyp、ならびに、エッジ強調度βjなどの画像処理条件は、制御手段８０へ送信される。
【０１４７】
一方、画像処理条件調整手段４０にて画像処理条件の調整がされると、画像処理条件設定手段２０は、その調整に対応する基準信号値の調整値Δs、コントラストの調整値Δr等を画像処理条件調整手段４０より受信する。
【０１４８】
そして、画像処理条件設定手段２０では、上記パラメータD1,D2等に対し、以下の式の処理がなされる。
D1'＝D1＋Δs    …（８），
D2'＝D2＋Δs    …（９），
k'＝k＋Δr・k  …（１０），
として、再度（１）〜（７）式が計算され、階調変換曲線LUTが再作成される。階調変換曲線が再作成されると、当該変更後の階調変換曲線LUT及び、その作成に使用したパラメータD1'、D2'、及びk'が制御手段８０へ送信される。
【０１４９】
同様に、画像処理条件調整手段４０にてエッジ強調度の調整量Δbを受信すると、このΔbが画像処理条件設定手段２０へ送信される。そして、画像処理条件設定手段２０は、下記式によって強調係数の修正を行い、制御手段８０へ修正後の強調係数βj'を送る。
【０１５０】
βj'＝βj＋Δb    …（１１），
さらに、周波数強調タイプ番号typが画像処理条件調整手段４０にて変更されると、画像処理条件設定手段２０は、当該画像処理条件調整手段４０から受信した前記強調タイプ番号typに基づき、パラメータ記憶手段３０から対応するβjの組を呼び出して、当該強調タイプ番号typとともに、制御手段８０へ送信する。
【０１５１】
なお、制御手段８０では、被写体番号bn、初期設定されているパラメータD１,D２,k,βj,typ及び最後の画像処理条件調整後の当該パラメータD１',D２',k',βj',typ'を一時保持手段により保持する。
【０１５２】
（４）パラメータ記憶：
パラメータ記憶を行うパラメータ記憶手段３０には、前記各被写体番号bn毎に、一組の画像処理条件を記憶する。ここで、各組の画像処理条件には、少なくとも以下のパラメータが含まれる。
【０１５３】
4-a）画像中から照射野領域を検出するためのアルゴリズムで使用するパラメータ、
4-b）照射野内に撮影されている被写体から、特定の構造を見つけるための複数種類のアルゴリズムを選択するためのパラメータ、
4-c）当該特定構造を見つけるアルゴリズムで使用するパラメータ、
4-d）当該特定構造部について、画像処理後に所定の濃度となるよう設定する基準信号値パラメータ、
である。
【０１５４】
具体的には、関心領域内の信号値の累積ヒストグラムにおける累積度数の比率h1,h2,及びh1に対応する信号値が出力画像でとる信号値D1,h2に対応する信号値が出力画像でとる信号値D2、である。
【０１５５】
4-e）基準となる階調変換曲線stdLUT、
具体的には、入力信号値をs[s＝０,１,..,４０９５]としたときの出力信号値を要素値とする１次元配列stdLUT[s]、
4-f）当該画像に対するエッジの強調度を決めるパラメータ
具体的には、以下の通りである。
【０１５６】
・強調タイプtyp、
・各周波数帯毎の強調係数β0,β1,..,βn、
である。
【０１５７】
なお、各パラメータは、パラメータ変更手段６０からの修正情報に基づいて、適宜修正される。また、これらと別に、強調タイプ毎の初期設定強調係数の組β0,β1,..,βnが記憶されている。
【０１５８】
（５）画像処理条件調整：
画像処理条件の調整を行う画像処理条件調整手段４０は、画像処理が施された縮小画像を表示する、画像表示手段４１と、画像の濃度、コントラスト等を修正入力する修正手段４２とによって構成される（図３参照）。
【０１５９】
まず、制御手段８０より取得した縮小画像を画像表示手段４１に表示する。そして同時に、基準信号値、または基準信号値から一意に求まる表示用の値を修正手段４２上に表示する。この様子を図３に示す。
【０１６０】
ユーザは、画像表示手段４１上に表示された画像または修正手段４２上に表示された値を見ながら、自分の好みの画像が得られるよう、濃度やコントラスト、エッジ強調度などを変更する（図５Ｓ２）。
【０１６１】
ここで修正手段４２は図３に示すような表示部とタッチパネル等とから構成されたユーザインタフェースとして構成されている。
ここで、ユーザは濃度変更ボタン４２１やコントラスト変更ボタン４２２を、押すことにより、画像の濃度やコントラストを変更することができる。
【０１６２】
なお、濃度変更ボタン４２１等を『押す』という動作は、マウスやタッチペンのようなポインティングデバイスを利用して、画面上の当該部分でクリック動作をすることで実行できる。また画像処理条件調整手段４０がタッチパネルを有する場合には、画面上の当該部分を指で触れることなどにより、当該動作を実現できる。
【０１６３】
パラメータの変更は、例えば上記濃度変更ボタン４２１を一度押すことにより、基準となる入力信号値に対する出力画像の信号値が一定値|Δs|だけ増減させることで行う。上記濃度変更ボタン４２１のうち、『Up』が押されると、Δsは正の値となり、逆に『Dw』が押されるとΔsは負の値となる。
【０１６４】
また同様に、上記コントラスト変更ボタン４２２を一度押すと、階調変換曲線の傾きを一定割合|Δr|だけ傾けることで行う。この際、上記コントラスト変更ボタン４２２のうち、『Up』が押されると、Δrは負の値となり、逆に『Dw』が押されるとΔrは正の値となる。
【０１６５】
同様に、修正手段４２上の強調度変更ボタン４２３を押すと、各周波数帯域毎の画像に対する強調係数βj（i＝０,１,..,n）それぞれに一定値|Δb|が増減される。上記強調度変更ボタン４２３のうち、『Up』が押されると、Δbは正の値となり、逆に『Down』が押されるとΔbは負の値となる。
【０１６６】
このようにして、パラメータの調整が行われる度に、パラメータの増減値Δs、Δr、Δbは、画像処理条件調整手段４０から、パラメータ修正量として、画像処理条件設定手段２０へ送信される。
【０１６７】
さらに、修正手段４２上の強調タイプ変更ボタン４２４を押すと、各βjのバランスを変更することができる。このβjのバランスは、たとえば、図４のように、type１、type２、type３のように、強調する周波数帯域が異なる複数のタイプが予め準備され、その中からタイプを選択することで行う。
【０１６８】
ここでは、強調タイプ変更ボタン４２４が一度押される度に、type１→type２→type３→type１のように、強調タイプが順次切り替わる構成になっていることが、操作性の点で望ましい。
【０１６９】
具体的には、強調タイプ変更ボタン４２４が押される度に、各タイプと１対１に対応する強調タイプ番号typ（本実施例では、type１＝１, type２＝２, type３＝３とする）が画像処理条件設定手段２０へ送信される。
【０１７０】
そして、前述したように、画像処理条件設定手段２０では、階調変換曲線LUTを再作成したり、強調係数の再計算を行い、当該再計算後の階調変換曲線等を制御手段８０へ送信する（図５Ｓ３）。
【０１７１】
その後、再作成された階調変換曲線に基づいて画像処理が行われた表示用縮小画像を再び画像処理条件調整手段４０で受信すると、当該再処理された表示用縮小画像が画像表示手段４１上に表示され、ユーザはパラメータの変更結果を確認することができるようになる。なお、最後に、ユーザによるパラメータの調整が終了した時点では、ユーザが終了ボタン４３を押すことにより、終了通知が制御手段８０へ通知される。
【０１７２】
（６）パラメータ調整量記憶：
パラメータ調整量の記憶を行うパラメータ調整量記憶手段５０では、上記画像処理条件調整手段４０でユーザが加えたパラメータ調整量の総量を、放射線画像入力手段１０で得た画像データ毎に保存する（図５Ｓ５）。
【０１７３】
また、以上のパラメータ調整量の総量の保存と合わせて、前回パラメータ変更手段６０にてパラメータ変更が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num1も保存する。さらに、以上のパラメータ調整量の総量とnum1との保存に合わせて、前回パラメータ修正許可判定手段７０にて修正許可判定が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num2も保存する。
【０１７４】
また、前回パラメータ変更手段６０が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num1及び前回パラメータ修正許可判定手段７０が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num2も一緒に記憶する。
【０１７５】
さらに、パラメータの調整を許可するか、不許可とするかを判別するパラメータ修正許可判定フラグとして、doLearningが記憶される。
doLearningの初期値は‘０’であり、この場合、パラメータの調整が可能となる。一方、doLearningの値が‘１'のときには、パラメータの調整は行われず、一定のパラメータが保持される。
【０１７６】
画像処理条件調整手段４０での調整が終了すると、制御手段８０より、被写体番号bn、初期設定されているパラメータD１,D２,k,βj,typ及び画像処理条件調整後の当該パラメータD１',D２',k',βj',typ'が、パラメータ調整量記憶許可判定手段９０を介してパラメータ調整量記憶手段５０へ送られる。
【０１７７】
そしてパラメータ調整量記憶手段５０では、各パラメータの調整総量ΔD１,ΔD２,Δk,Δβ,Tが以下のように求められる。
ΔD1＝D1'−D1    …（１２），
ΔD2＝D2'−D2    …（１３），
Δk＝k'／k    …（１４），
Δβ＝β0'／β0    …（１５），
T＝typ'[typ'＝|１（type１）,２（type２）,３（type３）]    …（１６），
ここで調整量は、被写体番号bnで特定されるパラメータセット毎に分類される。例えば記憶される調整量はC[bn]＝[Cbn１,Cbn２,...,Cbn２t]として表される。なお、各Cbni（i＝１,２,..,２t）は、（i−１）回前のパラメータ調整量の組を表し、以下のように表される。
【０１７８】
Cbni＝（ΔD1i, ΔD2i, Δki, Δβ, T）、
また、これらパラメータの調整量を求めるとともに、前回パラメータ変更手段６０が呼び出されてから、当該被写体番号bnに対応する画像処理が行われた回数num1について、
num1＝num1＋１、
のように１増加させる（図５Ｓ７）。
【０１７９】
同様に、前回パラメータ修正許可判定手段７０が呼び出されてから、当該被写体番号bnに対応する画像処理が行われた回数num2についても、
num2＝num2＋１、
のように１増加させる。
【０１８０】
ここで、あるパラメータセットについて、所定回数の閾値をＴｈｄn1とした場合（本実施形態においては、当該所定回数の閾値Ｔｈｄn1を１０とする）、num1＞Ｔｈｄn1、かつ、dolearning＝0のとき、パラメータ変更手段６０へC[bn]およびパラメータD1,D2,k,βj,typが送信され、パラメータの調整が開始される（図５Ｓ６、Ｓ８〜）。
【０１８１】
なお、各回のパラメータ調整量Cbniが所定数MaxNum以上記憶された場合、パラメータ調整量記憶手段５０では所定数MaxNumだけ各回のパラメータ調整量Cbniを保存し、これを超えた場合は、古い方から削除される。
【０１８２】
なお、dolearning＝1の場合、パラメータ変更手段６０は呼び出されず、ただ、num1＝0と入力される（図５Ｓ９）。
また、あるパラメータセットについて、所定回数の閾値をＴｈｄn2とした場合（本実施形態においては、当該所定回数の閾値Ｔｈｄn2を２０とする）、num2＞Ｔｈｄn2のとき、パラメータ修正許可判定手段７０へC[bn]およびパラメータD1,D2,k,βj,typ，dolearningが送信され、パラメータ変更を継続するか停止するか、既に停止している場合には再開するか停止を継続するか、の判定（パラメータ修正許可判定）が開始される（図５Ｓ１３、Ｓ１５〜）。
【０１８３】
また、以上の具体例では、パラメータ調整を開始させる（学習機能を働かせる）画像処理回数のＴｈｄn1の閾値を１０、修正許可判定を行わせる画像処理回数の閾値Ｔｈｄn2を２０とするものを示したが、
Ｔｈｄn2＞Ｔｈｄn1、
の条件を満たせば、修正許可判定の周期内に学習機能が働くため望ましい。また、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【０１８４】
また、以上の場合において、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との比Ｒｐを求め、該比Ｒｐが、１−Ｔｈｄ3＜Ｒｐ＜１＋Ｔｈｄ3の場合を所定の条件を満たしたとして、パラメータの変更を停止させるようにしてもよい。この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の比が所定の範囲になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１８５】
また、以上の場合において所定の条件を満たしたとして、パラメータの変更を停止させた後、さらに、パラメータ調整量から計算された調整後のパラメータ値と、当該時点でのパラメータ値との比Ｒｐを求め、該比Ｒｐが、１−Ｔｈｄ4＞Ｒｐ，または、１＋Ｔｈｄ4＜Ｒｐ，のいずれかの場合に、前記所定の条件の解除としてパラメータの変更を再開させるようにしてもよい。この結果、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１８６】
なお、パラメータ変更について自動停止と自動再開とを行う場合に、再開のための閾値Ｔｈｄ4を停止のための閾値Ｔｈｄ3より大きく設定することが望ましい。この結果、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。なお、この場合、差ではなく比によって判定しているため、パラメータの絶対値の変動が大きい場合であっても適切な判定が可能になる。
【０１８７】
また、パラメータ変更手段６０の呼び出しは、前回のパラメータ変更後、Cbniが１回更新されるごとに行ってもよいし、n回（ただしn≦t）更新される毎に行ってもよい。ここで、更新の回数が多ければ、迅速にパラメータの最適化が進むようになる。
【０１８８】
なお、上述したパラメータの調整開始に関する閾値Ｔｈｄn１を、変更する手段を別途設けることが望ましい。例えば、本実施形態にかかる画像処理装置をメンテナンスするための特別のモードで上述した画像の調整を行う場合は、Ｔｈｄn1＝１と設定すれば、１画像に対して画像処理条件の調整が行われる度に、画像処理パラメータの修正が行われるので、経験を積んだ放射線技師などが操作する場合には、早期に最適な画像処理パラメータへ修正することが可能となる。
【０１８９】
ここで画像処理装置をメンテナンスするための特別のモードを設定するためには、通常のモードか、特別のモードかを判別するモード判別フラグを制御手段中に準備しておき、当該画像処理装置の起動時に、ログイン操作を行わせ、その操作時に所定のパスワードを入力することにより、上述したモード判別フラグを特別のモードであることを示す値とすることにより実現できる。
【０１９０】
（７）パラメータ変更：
パラメータ変更手段６０は、パラメータ調整量C[bn]を受け取ると、当該パラメータセットについて、パラメータの最適化を行うため、調整を行う。
【０１９１】
具体的には、Cbni（i＝１,２,..,t）の各要素毎に平均値ΔD1av,ΔD2av,Δkav,Δβavを計算する。
ここで各要素の値で、最も絶対値の大きいものは、上記平均値の計算から除外する。これは、階調変換曲線の作成の際に、例えば照射野認識処理や関心領域設定処理で誤認識が発生することなどにより、本来予定されているのと異なる階調変換曲線を作成してしまうおそれがあり、このような誤った階調変換曲線を元に調整されたパラメータの悪影響を排除するためである。一方、Cbniの要素Tより、最頻値Tmodを求める。
【０１９２】
次に、各平均値ΔD1av, ΔD2av ,Δkavから、パラメータの修正量を算出する。基準信号値std1,std2の一方しか指定されない場合（h1,h2の一方のみ指定されている場合）、指定されている方の出力信号値D1またはD2について、
D1”＝ηΔD1av＋D1    …（１７），
もしくは
D2”＝ηΔD2av＋D2    …（１８），
が計算され、求められたD1”、D2”が新たなパラメータとなる。ここで、ηは、学習係数であり、学習の早さを決めるものである。ηの値は、０〜１のいずれかであり、数値が大きいほど早くパラメータが変更される。また、η＝０のときは、学習が行われない。
【０１９３】
同様に、コントラストの修正量Δkavから、コントラスト調整パラメータkが以下のように変更される。
k”＝ηΔkav ・ k    …（１９），
基準信号値std1,std2が両方指定される場合は、上記（１７）（１８）式を両方計算し、D1,D2双方とも変更する。
【０１９４】
次に、強調タイプのうち、最も使用されたものを特定するため、強調タイプの最頻値Tmodの値を調べる。Tmodが従前に設定された強調タイプtypと等しい場合、新たな各強調係数βj”は、以下のように求められる。
【０１９５】
βj”＝ηΔβav＋βj    …（２０），
一方Tmod≠typの場合、まずTmodの値に対応する強調係数βjのセットがパラメータ記憶手段３０から呼び出される。その後、呼び出された各βjに対して上記（２０）式が計算される。
【０１９６】
さらに、パラメータ調整量C[bn]に記憶されている各調整量について、変更後のパラメータを基準とした調整量に修正する。具体的には以下の通りである。
ΔD1i'＝ΔD1i−ηΔD1av    …（２２），
ΔD2i'＝ΔD2i−ηΔD2av    …（２３），
Δki'＝Δki／ηΔkav    …（２４），
Δβji'＝Δβji−ηΔβav 、（但し、Tmod＝typのとき），
Δβji'＝０、（但し、Tmod ≠ typのとき）    …（２５），
Ti'＝Ti    …（２６），
ここで、変更後のパラメータセットは、パラメータ変更手段６０からパラメータ記憶手段３０へ送られ、以後変更後のパラメータセットにより、画像処理条件が決定されることとなる（図５Ｓ１０，Ｓ１１）。
【０１９７】
そして、最後に、パラメータ変更手段６０は、パラメータ調整量記憶手段５０にて記憶されている、前回のパラメータ調整からの当該画像処理条件の呼び出し回数num1について、
num1＝0、
と修正する（図５Ｓ１２）。
【０１９８】
（８）パラメータ修正許可判定：
パラメータ修正許可判定手段７０では、パラメータ変更手段６０と同様にして、（１７）〜（２０）式に従って、画像処理条件変更後のパラメータ値D1'、D2'、k'等を算出する。
【０１９９】
(8-1)パラメータ変更不許可：
そして、これら画像処理条件変更後のパラメータ値と、画像処理条件調整前のパラメータ値とを比較し、全ての調整対象パラメータについて、各閾値ＴｈｄCrrctD、ＴｈｄCrrctK、ＴｈｄCrrctB、を定めて、下記条件を満たす場合、前記パラメータ修正許可判定フラグdolearningを1に変更する（図５Ｓ１５，１６）。すなわち、変更不許可とする。
｜D1'−D1｜＜ＴｈｄCrrectD、
｜D2'−D2｜＜ＴｈｄCrrectD、
｜k'−k｜＜ＴｈｄCrrectK、
｜βj'−βj｜＜ＴｈｄCrrectB、
このようにすることで、ユーザによるパラメータ調整の結果から算出したパラメータが、元のパラメータと一定の差異しかないような場合に、自動的にパラメータ変更動作を収束させることができる。一旦収束すれば、その後の撮影においては、ユーザによって最適なパラメータを用いて画像処理を施すことが可能である。
【０２００】
(8-2)パラメータ変更許可：
一方、これら画像処理条件変更後のパラメータ値と、画像処理条件調整前のパラメータ値とを比較し、各閾値ＴｈｄUnCrrctD、ＴｈｄUnCrrctK、ＴｈｄUnCrrctBを定めて、いずれか一つでも下記条件を満たす場合、前記パラメータ修正許可判定フラグdolearningを0に変更する（図５Ｓ１５，Ｓ１７）。すなわち、変更許可＝変更再開とする。
｜D1'−D1｜＞ＴｈｄUnCrrectD、
｜D2'−D2｜＞ＴｈｄUnCrrectD、
｜k'−k｜＞ＴｈｄUnCrrectK、
｜βj'−βj｜＞ＴｈｄUnCrrectB、
このようにすることで、パラメータ変更動作が一旦収束しているとしても、ユーザによるパラメータ調整量が大きくなった場合に、再度パラメータ変更を行うように自動的に判別することによって、ユーザにとって最適なパラメータに変更していくことが可能である。
【０２０１】
(8-3)パラメータ修正許可判定時の各種条件：
パラメータ修正許可判定では、以上の各閾値ＴｈｄCrrctD、ＴｈｄCrrctK、ＴｈｄCrrctB、ＴｈｄUnCrrctD、ＴｈｄUnCrrctK、ＴｈｄUnCrrctB、について、以下の条件を満たすことが望ましい。
ＴｈｄCrrctD＜ＴｈｄUnCrrctD、
ＴｈｄCrrctK＜ＴｈｄUnCrrctK、
ＴｈｄCrrctB＜ＴｈｄUnCrrctB、
以上の不許可のための閾値と許可（変更再開）のための閾値とを上記条件とすることで、パラメータ変更の停止・再開を頻繁に繰り返すことによる不安定さを解消することができ、安定性を高めることが可能になる。
【０２０２】
なお、本件出願の発明者が本実施形態について具体的な実験を行ったところ以下の値の範囲で望ましい結果が得られた。
ＴｈｄCrrctD＝０．０３〜０．１、
ＴｈｄCrrctK＝０．１〜０．３、
ＴｈｄCrrctB＝０．０５〜０．１、
ＴｈｄUnCrrctD＝０．１〜０．２、
ＴｈｄUnCrrctK＝０．３〜０．６、
ＴｈｄUnCrrctB＝０．１〜０．２、
なお、この値は発明者が実験を行った条件で望ましい結果が得られた具体例であり、この値に限定されるものではない。
【０２０３】
(９)パラメータ調整量記憶許可判定：
パラメータ調整量記憶許可判定手段９０では、制御手段８０より受信したパラメータ一式をパラメータ調整量記憶手段５０へ送信するか否かを、記憶許可判定フラグpermitfの値により判定する（図５Ｓ４）。この記憶許可判定を行う実施形態の具体例としては、以下のものがある。
【０２０４】
(9-1)パラメータ調整量記憶許可判定の具体例：
特定のユーザによる画像調整の結果だけをパラメータ修正に反映させたい場合、パラメータ調整量記憶許可判定手段９０では、ユーザのID番号に対応した記憶許可判定フラグpermitfを保持しておく。なお、ユーザID番号に対応した記憶許可判定フラグpermitfの値は、以下のような周知の手段により、取得できる。
【０２０５】
まず、本実施形態にかかる画像処理装置の使用をできるユーザを予め登録したユーザに限定するとともに、各ユーザ毎にユニークなID番号を設定する。そのユーザ登録時において、パラメータの修正を行うか否かの情報を、ユーザインタフェースを通じて、ユーザに選択させる。
【０２０６】
例えば、『パラメータの修正を認める』というチェックボックスをユーザインタフェース上に準備し、チェックがなされた場合は、パラメータ修正を許可するという意味で‘１'の値を、そうでない場合はパラメータ修正を許可しないという意味で‘０’の値をユーザID番号に関連付けて記憶しておく。
【０２０７】
そして、当該画像処理装置の使用を開始する際、『ログイン』操作のように、ユーザID番号をユーザインタフェースから当該画像処理装置に入力することにより、使用開始の許可を与える操作をユーザに行わせる。その操作によって入力されたID番号と、ユーザ情報記録手段に記憶されている全てのユーザID番号とを照合し、一致するユーザID番号が存在する場合は、上記ユーザ情報記録手段から、当該ユーザID番号に対応する記憶許可判定フラグの値を読み込む。
【０２０８】
そして、当該操作に伴い、制御手段８０を介して、パラメータ調整量記憶許可判定手段９０に、記憶許可判定フラグの値が送信される。そして、パラメータ調整量記憶許可判定手段９０では、当該記憶許可判定フラグの値をフラグ保持メモリ９１（図示せず）に保管する。
【０２０９】
制御手段８０より、パラメータ一式を受信すると、パラメータ調整量記憶許可判定手段９０では、一時記憶メモリ９２（図示せず）などに当該パラメータ一式を保持する。その後、フラグ保持メモリ９１から記憶許可判定フラグpermitfを読み出し、値を参照する。
【０２１０】
ここで、permitfの値が‘１'の場合は、パラメータ調整量記憶許可判定手段９０は、制御手段８０から受け取った上記パラメータ一式を、パラメータ調整量記憶手段５０へ送信する。
【０２１１】
一方、permitfの値が‘０’の場合は、パラメータ調整量記憶許可判定手段９０は、制御手段８０から受け取った上記一時記憶メモリ９２に保持されたパラメータ一式を消去する。
【０２１２】
このような構成とすることで、熟練した放射線技師のように、より有益な情報を与えるユーザのみを対象としてパラメータの修正を行わせることが可能となり、診断に適した画像を得るためのパラメータへ、少ない学習回数で到達することが可能となる。
【０２１３】
(9-2)パラメータ調整量記憶許可判定の具体例：
パラメータ調整量記憶許可判定手段９０において、ユーザの使用開始操作後、記憶許可判定フラグの値をフラグ保持メモリ９１に保存するまでは、以上の具体例と同一の構成となる。
【０２１４】
一方、制御手段８０より、パラメータ一式を受信すると、パラメータ調整量記憶許可判定手段９０では、フラグ保持メモリ９１から読み出した記憶許可判定フラグpermitfについて、当該パラメータ一式に追加してパラメータ調整量記憶手段５０へ送信する。
【０２１５】
この場合、パラメータ調整量記憶手段５０では、当該パラメータ一式に加えて、その都度記憶許可判定フラグの値も記憶する。
また、パラメータ調整量記憶手段５０では、前回パラメータ変更手段６０が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num１及び前回パラメータ修正許可判定手段７０が実行されてから、当該被写体番号bnに対応する画像処理が行われた回数num２については、それぞれパラメータ調整量記憶許可判定手段９０にて追加された記憶許可判定フラグpermitfの値が‘１'の場合だけ、それぞれ更新し、そうでない場合は直前の値を維持する。
【０２１６】
パラメータ変更手段６０、パラメータ修正許可判定手段７０では、上記パラメータ一式とともに記憶されている記憶許可判定フラグpermitfの値が‘１'の場合だけ、当該パラメータを計算に使用し、permitfの値が‘０’の場合は、当該パラメータの値を無視することにする。
【０２１７】
このような構成とすることで、パラメータ調整量記憶手段５０に記憶されたパラメータ一式を、例えばファイルの形として別途取得することにより、パラメータの修正が行われた否かに関わらず追跡調査することができ、装置のメンテナンス上有用な情報を得ることができる。
【０２１８】
(9-3)パラメータ調整量記憶許可判定の具体例：
上述のような画像処理装置が複数ネットワークに接続されているシステムにおいて、画像処理のパラメータを共有して使用する場合、特定の画像処理装置でのみパラメータ修正を行う。パラメータの修正が行われると、共有されている画像処理パラメータ保持手段（以上の実施形態のパラメータ記憶手段３０に相当）へ修正済みのパラメータを送信し、パラメータの更新がなされる。
【０２１９】
なお、パラメータの共有を行うシステムは、特開２００２−２７９３９５号公報に記載されたもの等がある。パラメータ修正の手順は、以下の通りである。
(9-3i)パラメータの学習を行う画像処理装置では、上述したように、パラメータ調整量を記憶しておくとともに、所定の手順に従って修正パラメータを計算する。
【０２２０】
(9-3ii)修正パラメータが求まると、パラメータの学習を行う画像処理装置から、画像処理パラメータ保持手段へパラメータ修正要求を出す。
(9-3iii)パラメータ修正要求を受けた画像処理パラメータ保持手段では、一時記憶手段に修正前の画像処理パラメータを複写する。
【０２２１】
(9-3iv)一時記憶手段への複写が完了し、ネットワークに接続された他の画像処理装置からのアクセスが途絶えると、上記パラメータ修正要求を出した画像処理装置へ、パラメータ修正許可を与える。
【０２２２】
(9-3v)パラメータ修正許可を受けた画像処理装置は、修正後のパラメータを画像処理パラメータ保持手段へ送信する。
(9-vi)画像処理パラメータ保持手段では、保持している画像処理パラメータを上記修正後のパラメータで置換する。
【０２２３】
(9-vii)置換が正常に終了すると、上述した一時記憶手段に保存されている画像処理パラメータを消去し、パラメータの更新が完了する。このとき、置換が正常に終了したか否かは、画像処理パラメータ保持手段が受信した、修正後の画像処理パラメータのデータサイズが所定値か否かで判断できる。また、パラメータの置換が不正と判断された場合は、一時記憶手段に保存されている修正前の画像処理パラメータを元のメモリに複写することで画像処理パラメータを元に戻す。
【０２２４】
(9-viii)画像処理パラメータの更新が正常に終了すると、ネットワーク上に接続された各画像処理装置に対し、パラメータ更新通知を送信する。
(9-ix)各画像処理装置は、画像処理パラメータ保持手段から、更新された画像処理パラメータを取得する。
【０２２５】
なお、画像処理パラメータ保持手段への、上述した修正要求は、所定のタイミングに限定することも可能である。例えば、画像処理パラメータの学習を行った画像処理装置の起動時、終了時の手順中で行ったり、ユーザが指定した時間に行うようにすることも可能である。
【０２２６】
この場合、画像処理パラメータの学習が行われると、パラメータ修正要求を一時的にパラメータ修正要求保持手段に保存しておき、所定のタイミングにてパラメータ修正要求保持手段を参照し、パラメータ修正要求が存在すれば画像処理パラメータ保持手段へパラメータ修正要求を送信するということで実現できる。
【０２２７】
このような形態とすることで、特定の画像処理装置を、熟練した放射線技師が使用することにより、当該技師による使用の結果、修正されたパラメータを他の装置に反映させることができ、システム全体の画像処理パラメータについて、診断に最適な画像を得るためのパラメータへ、少ない学習回数で到達することが可能となる。
【０２２８】
〈本発明の実施形態の処理例とそれによる効果〉
ユーザが調整した調整量あるいは調整値によってパラメータが徐々に変更されていくため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定でき、自動的に最適化することが可能になる。
【０２２９】
また、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整前後の差が所定量未満になるように収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。また、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更していき、かつ、調整が所定の条件で収束した場合には変更を停止させることで、自動的にかつ安定した状態で最適化することが可能になる。そして、放射線画像に対する画像処理条件をユーザによる調整結果を利用することにより自動的に変更する際に、調整が収束した場合には変更を一時停止させているが、調整前後の差が所定量を超えるようになった場合には設定されている状態とユーザが望む状態とに乖離が生じているので変更を自動的に再開させることで、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。さらに、停止の閾値と再開の閾値を設定することで、停止と再開とを不用意に繰り返すことがなくなり、自動的にかつ安定した状態を維持しつつ最適化することが可能になる。
【０２３０】
また、以上のパラメータの調整量について、パラメータ調整量記憶をするか否かを判定するようにしている。すなわち、パラメータの学習に使用するデータを制限する手段を設けることにより、本来学習に用いるべきでないデータを学習から除外できるので、効率よくパラメータの学習を行うことができる。この結果、自動的に最適化する際、画像調整結果を反映する対象を特定することにより、複数の画像処理装置や複数のユーザが使用する場合でも簡便かつ安定的に画像処理条件の調整が行われるようになる。
【０２３１】
〈その他の実施形態（１）〉
以上の実施形態の説明では、パラメータ調整量を用いて各種の処理を行うようにしていたが、パラメータ調整量の代わりに、調整後のパラメータ（修正パラメータ）を記憶し、上記処理と同等の処理を行うようにしても、同様の効果を得ることができる。
【０２３２】
〈その他の実施形態（２）〉
本発明においては、画像処理のパラメータを自動的に最適化する装置について説明してきたが、場合によっては、画像処理のパラメータを初期値に戻したい場合も考えられる。そのため、画像処理パラメータを初期値へ戻す手段を設けておくことが望ましい。
【０２３３】
かかる手段としては、周知の手段により実現できる。例えば、予め画像処理パラメータの初期値をパラメータ記憶手段３０中に保存しておき、ユーザインタフェースにパラメータ初期化指示手段を設け、当該手段によって初期化指示がなされた場合に、パラメータ記憶手段３０に保存されている画像処理パラメータの初期値で、画像処理パラメータを置換することにより、実現できる。
【０２３４】
〈その他の実施形態（３）〉
以上の実施形態において、パラメータ調整量記憶許可判定手段９０とパラメータ修正許可判定手段７０とを設けない構成であっても、ユーザが調整した調整量によってパラメータが徐々に変更されていく。このため、記憶されているパラメータがユーザが望んでいる画像処理の状態に近づいていくことになり、ユーザが手動で画像処理条件の事前調整を行うことなく、最適な画像処理条件を決定できる効果が得られ、上述した第１の課題を解決することができる。
【０２３５】
なお、パラメータ修正許可判定手段７０を設けることで、自動的なパラメータ値の変更が停止・再開可能になるため、安定した動作が望めるようになり、第２の課題を解決することができるようになる。
【０２３６】
さらに、パラメータ修正許可判定手段７０とパラメータ調整量記憶許可判定手段９０とを設けることで、自動的なパラメータ値の変更が停止・再開可能になり、さらに許可された場合に学習機能が働くため、特定の画像処理装置を、熟練した放射線技師が使用することにより、当該技師による使用の結果、修正されたパラメータを他の装置に反映させることができ、システム全体の画像処理パラメータについて、診断に最適な画像を得るためのパラメータへ、少ない学習回数で到達することが可能となり、第３の課題をも解決することができるようになる。
【図面の簡単な説明】
【０２３７】
【図１】本発明の一実施形態の全体構成あるいは全体処理の流れを示す説明図である。
【図２】本発明の実施形態における被写体情報選択の様子を示す概要図である。
【図３】本発明の実施形態における画像処理条件調整の様子を示す概要図である。
【図４】本発明の実施形態における各強調タイプ毎の強調係数の様子を示す概要図である。
【図５】本発明の実施形態における動作例を示すフローチャートである。
【符号の説明】
【０２３８】
１０  放射線画像入力手段
２０  画像処理条件設定手段
３０  パラメータ記憶手段
４０  画像処理条件調整手段
５０  パラメータ調整量記憶手段
６０  パラメータ変更手段
７０  パラメータ修正許可判定手段
８０  制御手段（画像処理手段を含む）
９０  パラメータ調整量記憶許可判定手段
【図１】
2005111249.tif 000003
【図２】
2005111249.tif 000004
【図３】
2005111249.tif 000005
【図４】
2005111249.tif 000006
【図５】
2005111249.tif 000007