JP2004178344A - 高次モーメントを用いた画像投影方法及び画像処理装置 - Google Patents
高次モーメントを用いた画像投影方法及び画像処理装置 Download PDFInfo
- Publication number
- JP2004178344A JP2004178344A JP2002344974A JP2002344974A JP2004178344A JP 2004178344 A JP2004178344 A JP 2004178344A JP 2002344974 A JP2002344974 A JP 2002344974A JP 2002344974 A JP2002344974 A JP 2002344974A JP 2004178344 A JP2004178344 A JP 2004178344A
- Authority
- JP
- Japan
- Prior art keywords
- projection
- image
- dimensional data
- axis
- order moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000013500 data storage Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 6
- 210000000988 bone and bone Anatomy 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
- Image Generation (AREA)
Abstract
【課題】３次元データから作成した投影画像に、投影軸上の全てのデータの値を反映させる画像投影方法を提供する。
【解決手段】投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とする。
【効果】投影軸上の最大値だけでなく、それ以外のデータの値も投影画像に反映させうる。また、投影軸上に最大値を持つデータがただ１点だけなのか複数点あるのかの情報も投影画像に反映されうる。
【選択図】 図３
【解決手段】投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とする。
【効果】投影軸上の最大値だけでなく、それ以外のデータの値も投影画像に反映させうる。また、投影軸上に最大値を持つデータがただ１点だけなのか複数点あるのかの情報も投影画像に反映されうる。
【選択図】 図３
Description
【０００１】
【発明の属する技術分野】
本発明は、高次モーメントを用いた画像投影方法および画像処理装置に関し、さらに詳しくは、３次元データから作成した投影画像に投影軸上の全てのデータの値を反映することが出来る高次モーメントを用いた画像投影方法および画像処理装置に関する。
【０００２】
【従来の技術】
３次元データから投影画像を作成する画像投影方法として、最大値投影方法（Ｍａｘｉｍｕｍ Ｉｎｔｅｎｓｉｔｙ Ｐｒｏｊｅｃｔｉｏｎ）が知られている。
最大値投影方法は、投影面に対して垂直な軸上での３次元データ値の最大値を、その軸と投影面の交点の画素値とする画像投影方法であり、ＭＲＩ（ＭａｇｎｅｔｉｃＲｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置，Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）装置，超音波診断装置において血管の表示方法として使用されている（「医用画像辞典、株式会社日経メディカル開発発行、日経ＢＰ出版センター発売」による）。
【０００３】
【発明が解決しようとする課題】
上記最大値投影方法では、投影画像には最大値だけが反映され、それ以外のデータの値は全く反映されない問題点がある。また、最大値を持つデータがただ１点だけなのか複数点あるのかの情報も反映されない問題点がある。
そこで、本発明の目的は、投影画像に投影軸上の全てのデータの値を反映することが出来る高次モーメントを用いた画像投影方法および画像処理装置を提供することにある。
【０００４】
【課題を解決するための手段】
第１の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第１の観点による高次モーメントを用いた画像投影方法では、投影面の画素値を上記Ｐに基づいて決めるが、このＰは、（ΣＶｉ）ｒからＶｉｒを除いた値であり、全てのＶｉを含んでいる。よって、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来る。
【０００５】
なお、Ｖｉｒを除く理由は、一つの大きな値のデータが支配的になるのを避けるためである。例えばデータがＶ１，Ｖ２の２個とし、ｒ＝２とすると、（ΣＶｉ）ｒ＝Ｖ１２＋２・Ｖ１・Ｖ２＋Ｖ２２となるが、Ｖ１≫Ｖ２なら、（ΣＶｉ）ｒ≒Ｖ１２になってしまい、Ｖ２が反映されない。しかし、（ΣＶｉ）ｒ−ΣＶｉｒ＝２・Ｖ１・Ｖ２となり、Ｖ１≫Ｖ２でも、Ｖ２が反映される。
【０００６】
ちなみに、Ｐ＝ΣＶｉ／ｎとした場合、全てのＶｉを含んでいるため、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来るように見えるが、実際は一つの大きな値のデータが支配的になってしまう。例えばデータがＶ１，Ｖ２の２個とすると、ΣＶｉ＝Ｖ１＋Ｖ２となるが、Ｖ１≫Ｖ２なら、ΣＶｉ≒Ｖ１になってしまい、Ｖ２が反映されない。
【０００７】
第２の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第２の観点による高次モーメントを用いた画像投影方法は、前記第１の観点による高次モーメントを用いた画像投影方法におけるＰをそのまま画素値Ｇとするものである。
【０００８】
第３の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第３の観点による高次モーメントを用いた画像投影方法は、前記第１の観点による高次モーメントを用いた画像投影方法におけるＰの指数関数の値を画素値Ｇとするものである。
【０００９】
第４の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第４の観点による高次モーメントを用いた画像投影方法では、投影面の画素値を上記Ｐに基づいて決めるが、このＰは全てのＶｉを含んでいる。よって、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来る。
【００１０】
第５の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第５の観点による高次モーメントを用いた画像投影方法は、前記第４の観点による高次モーメントを用いた画像投影方法におけるＰをそのまま画素値Ｇとするものである。
【００１１】
第６の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第６の観点による高次モーメントを用いた画像投影方法は、前記第４の観点による高次モーメントを用いた画像投影方法におけるＰの指数関数の値を画素値Ｇとするものである。
【００１２】
第７の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、２≦ｒ≦１２８であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第７の観点による高次モーメントを用いた画像投影方法では、次数ｒ＝２，３，４，…とすることで、投影画像のコントラストが変わってくるため、目的に応じたコントラストとなる次数ｒを選べばよい。ただし、ｒ＝１２８，１２９，１３０，…としても、投影画像のコントラストはほとんど変わらないため、実用上ｒ＝１２８までで足る。
【００１３】
第８の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、操作者がｒを変更可能であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第８の観点による高次モーメントを用いた画像投影方法では、操作者が次数ｒを変更可能であるため、操作者の所望のコントラストとなる次数ｒを選択できる。
【００１４】
第８の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、２≦ｒ≦１２８であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第８の観点による高次モーメントを用いた画像投影方法では、前記第２の観点による高次モーメントを用いた画像投影方法において説明した理由により、実用上十分な次数ｒの範囲で投影画像のコントラストを選択できる。
【００１５】
第９の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第９の観点による画像処理装置では、前記第１の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１６】
第１０の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１０の観点による画像処理装置では、前記第２の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１７】
第１１の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１１の観点による画像処理装置では、前記第３の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１８】
第１２の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１２の観点による画像処理装置では、前記第４の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１９】
第１３の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１３の観点による画像処理装置では、前記第５の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２０】
第１４の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１４の観点による画像処理装置では、前記第６の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２１】
第１５の観点では、本発明は、上記構成の画像処理装置において、２≦ｒ≦１２８であることを特徴とする画像処理装置を提供する。
上記第１５の観点による画像処理装置では、前記第７の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２２】
第１６の観点では、本発明は、上記構成の画像処理装置において、操作者がｒを指定するための次数指定手段を具備したことを特徴とする画像処理装置を提供する。
上記第１６の観点による画像処理装置では、前記第８の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２３】
【発明の実施の形態】
以下、本発明の実施の形態を説明する。なお、これにより本発明が限定されるものではない。
【００２４】
−第１の実施形態−
図１は、第１の実施形態にかかる医用画像診断装置の構成図である。
この医用画像診断装置１００は、撮像装置１と、画像処理装置２とを具備している。
【００２５】
撮像装置１は、Ｘ線ＣＴ装置やＭＲＩ装置や超音波診断装置などであり、被検体Ｋを撮像して得たデータを画像処理装置２に渡す。
【００２６】
画像処理装置２は、撮像装置１から渡されたデータを基に３次元データを構築し記憶する３次元データ構築・記憶部２ａと、操作者が投影方向を指定するための投影方向指定部２ｂと、操作者が次数ｒを指定するための次数指定部２ｃと、高次モーメントを用いた画像投影演算を行なう投影演算部２ｄと，投影画像をディスプレイ画面に表示する投影画像表示部２ｅとを具備している。
【００２７】
図２は、次数ｒを変更するために操作者が操作するスライド・レバーの外観図である。
スライド・レバーを動かすことにより、次数ｒを「２」から「１２８」まで変更できる。
【００２８】
図３は、画像処理装置２における高次モーメントを用いた画像投影処理を示すフローチャートである。
ステップＳＴ１では、３次元データ構築・記憶部２ａは、撮像装置１から渡されたデータを基に３次元データを構築し、記憶する。
【００２９】
ステップＳＴ２では、投影方向指定部２ｂは、投影方向を指定するために操作者が操作するデバイス（例えば、トラックボール）から投影方向を読み込む。
ステップＳＴ３では、投影演算部２ｄは、投影方向に垂直な投影面を設定する。
【００３０】
ステップＳＴ４では、投影演算部２ｄは、投影面上の一つの画素に着目する。
ステップＳＴ５では、着目した画素に対応する投影軸上のｎ個のデータＶｉを３次元データより取り出す。
ステップＳＴ６では、次数指定部２ｃは、次数ｒを指定するために操作者が操作するデバイス（例えば、図２のスライド・レバー）から次数ｒを読み込む。
【００３１】
ステップＳＴ７では、投影演算部２ｄは、画素値Ｇを次式により算出する。
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
【００３２】
ステップＳＴ８では、投影演算部２ｄは、全画素の画素値Ｇを得るまで、上記ステップＳＴ４〜ＳＴ７を繰り返す。
ステップＳＴ９では、投影画像表示部２ｅは、得られた投影画像をディスプレイ画面に表示する。
ステップＳＴ１０では、操作者が「終了」を指示すれば処理を終了し、「終了」を指示しなければ上記ステップＳＴ２に戻る。
【００３３】
図４，図５は、高次モーメントを用いた画像投影演算の数値例を２次元的に示した説明図である。次数ｒ＝２としている。
図４に、高次モーメントを用いた画像投影により、３次元データＴＤ１を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。投影軸ａ，ｂ，ｃ，ｄは、それぞれ画素Ａ，Ｂ，Ｃ，Ｄに対応する投影軸である。
図５に、高次モーメントを用いた画像投影により、３次元データＴＤ２を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
【００３４】
図６，図７は、最大値投影方法を用いた画像投影演算の数値例を２次元的に示した説明図である。
図６に、最大値投影方法を用いた画像投影により、３次元データＴＤ１を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
図７に、最大値投影方法を用いた画像投影により、３次元データＴＤ２を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
【００３５】
高次モーメントを用いた図４と図５とを比較すれば判るように、投影軸上の最大値以外のデータの値も投影画像に反映されている（例えば、投影軸上の最小値が「１０」か「０」かによって画素値が異なっている）。一方、最大値投影方法を用いた図６と図７とを比較すれば判るように、投影軸上の最大値以外のデータの値は投影画像に全く反映されていない（例えば、投影軸上の最小値が「１０」か「０」かによらず画素値は最大値「７０」になっている）。
このことは、投影方向に例えば骨だけが存在しているのか骨と血管が重なって存在しているのかを、最大値投影方法では投影画像上で判別できないが、本発明では投影画像上で判別できることを意味している。
【００３６】
また、図４と図６とを比較すれば判るように、高次モーメントを用いた図４では、最大値を持つデータがただ１点だけなのか複数点あるのかの情報が投影画像に反映されている（投影軸上に最大値「７０」が存在する個数によって画素値が異なっている）。一方、最大値投影方法を用いた図６では、最大値を持つデータがただ１点だけなのか複数点あるのかの情報が投影画像に全く反映されていない（投影軸上に最大値「７０」が存在する個数によらず画素値は「７０」になっている）。
図５と図７とを比較しても、同じことが判る。
このことは、投影方向に例えば骨が１本あるのか複数本重なって存在しているのかを、最大値投影方法では投影画像上で判別できないが、本発明では投影画像上で判別できることを意味している。
【００３７】
−第２の実施形態−
図３のステップＳＴ７において、投影演算部２ｄは、画素値Ｇを次式により算出してもよい。
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映できる。
【００３８】
−第３の実施形態−
図８は、画像処理装置２における高次モーメントを用いた画像投影処理を示すフロー図である。
このフロー図は、図３のステップＳＴ７がステップＳＴ７’になっている以外は図３のフロー図と同じである。そこで、ステップＳＴ７’のみを説明する。
ステップＳＴ７’において、投影演算部２ｄは、画素値Ｇを次式により算出する。
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映でき、第１の実施形態よりも結果が良い場合があった。
【００３９】
−第４の実施形態−
図８のステップＳＴ７’において、投影演算部２ｄは、画素値Ｇを次式により算出してもよい。
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映でき、第２の実施形態よりも結果が良い場合があった。
【００４０】
なお、第１の実施形態〜第４の実施形態のいずれの画素値Ｇとするか、その他の関数Ｇ（Ｐ）とするかは、画像を作成する目的や好みに応じて適宜選べばよい。
【００４１】
【発明の効果】
本発明の高次モーメントを用いた画像投影方法および画像処理装置によれば、３次元データから作成した投影画像に、投影軸上の全てのデータの値を反映することが出来る。
【図面の簡単な説明】
【図１】本発明にかかる医用画像診断装置の構成図である。
【図２】次数ｒを操作者が変更するために操作するスライド・レバーの外観図である。
【図３】第１の実施形態にかかる高次モーメントを用いた画像投影処理を示すフローチャートである。
【図４】高次モーメントを用いた画像投影演算の数値例を２次元的に示した説明図である。
【図５】高次モーメントを用いた画像投影演算の数値例を２次元的に示した別の説明図である。
【図６】最大値投影方法を用いた画像投影演算の数値例を２次元的に示した説明図である。
【図７】最大値投影方法を用いた画像投影演算の数値例を２次元的に示した別の説明図である。
【図８】第３の実施形態にかかる高次モーメントを用いた画像投影処理を示すフローチャートである。
【符号の説明】
１ 撮像装置
２ 画像処理装置
２ａ ３次元データ構築・記憶部
２ｂ 投影方向指定部
２ｃ 次数指定部
２ｄ 投影演算部
２ｅ 投影画像表示部
１００ 医用画像診断装置
【発明の属する技術分野】
本発明は、高次モーメントを用いた画像投影方法および画像処理装置に関し、さらに詳しくは、３次元データから作成した投影画像に投影軸上の全てのデータの値を反映することが出来る高次モーメントを用いた画像投影方法および画像処理装置に関する。
【０００２】
【従来の技術】
３次元データから投影画像を作成する画像投影方法として、最大値投影方法（Ｍａｘｉｍｕｍ Ｉｎｔｅｎｓｉｔｙ Ｐｒｏｊｅｃｔｉｏｎ）が知られている。
最大値投影方法は、投影面に対して垂直な軸上での３次元データ値の最大値を、その軸と投影面の交点の画素値とする画像投影方法であり、ＭＲＩ（ＭａｇｎｅｔｉｃＲｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置，Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）装置，超音波診断装置において血管の表示方法として使用されている（「医用画像辞典、株式会社日経メディカル開発発行、日経ＢＰ出版センター発売」による）。
【０００３】
【発明が解決しようとする課題】
上記最大値投影方法では、投影画像には最大値だけが反映され、それ以外のデータの値は全く反映されない問題点がある。また、最大値を持つデータがただ１点だけなのか複数点あるのかの情報も反映されない問題点がある。
そこで、本発明の目的は、投影画像に投影軸上の全てのデータの値を反映することが出来る高次モーメントを用いた画像投影方法および画像処理装置を提供することにある。
【０００４】
【課題を解決するための手段】
第１の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第１の観点による高次モーメントを用いた画像投影方法では、投影面の画素値を上記Ｐに基づいて決めるが、このＰは、（ΣＶｉ）ｒからＶｉｒを除いた値であり、全てのＶｉを含んでいる。よって、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来る。
【０００５】
なお、Ｖｉｒを除く理由は、一つの大きな値のデータが支配的になるのを避けるためである。例えばデータがＶ１，Ｖ２の２個とし、ｒ＝２とすると、（ΣＶｉ）ｒ＝Ｖ１２＋２・Ｖ１・Ｖ２＋Ｖ２２となるが、Ｖ１≫Ｖ２なら、（ΣＶｉ）ｒ≒Ｖ１２になってしまい、Ｖ２が反映されない。しかし、（ΣＶｉ）ｒ−ΣＶｉｒ＝２・Ｖ１・Ｖ２となり、Ｖ１≫Ｖ２でも、Ｖ２が反映される。
【０００６】
ちなみに、Ｐ＝ΣＶｉ／ｎとした場合、全てのＶｉを含んでいるため、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来るように見えるが、実際は一つの大きな値のデータが支配的になってしまう。例えばデータがＶ１，Ｖ２の２個とすると、ΣＶｉ＝Ｖ１＋Ｖ２となるが、Ｖ１≫Ｖ２なら、ΣＶｉ≒Ｖ１になってしまい、Ｖ２が反映されない。
【０００７】
第２の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第２の観点による高次モーメントを用いた画像投影方法は、前記第１の観点による高次モーメントを用いた画像投影方法におけるＰをそのまま画素値Ｇとするものである。
【０００８】
第３の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第３の観点による高次モーメントを用いた画像投影方法は、前記第１の観点による高次モーメントを用いた画像投影方法におけるＰの指数関数の値を画素値Ｇとするものである。
【０００９】
第４の観点では、本発明は、３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第４の観点による高次モーメントを用いた画像投影方法では、投影面の画素値を上記Ｐに基づいて決めるが、このＰは全てのＶｉを含んでいる。よって、投影画像に投影軸上の全てのデータＶｉの値を反映することが出来る。
【００１０】
第５の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第５の観点による高次モーメントを用いた画像投影方法は、前記第４の観点による高次モーメントを用いた画像投影方法におけるＰをそのまま画素値Ｇとするものである。
【００１１】
第６の観点では、本発明は、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第６の観点による高次モーメントを用いた画像投影方法は、前記第４の観点による高次モーメントを用いた画像投影方法におけるＰの指数関数の値を画素値Ｇとするものである。
【００１２】
第７の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、２≦ｒ≦１２８であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第７の観点による高次モーメントを用いた画像投影方法では、次数ｒ＝２，３，４，…とすることで、投影画像のコントラストが変わってくるため、目的に応じたコントラストとなる次数ｒを選べばよい。ただし、ｒ＝１２８，１２９，１３０，…としても、投影画像のコントラストはほとんど変わらないため、実用上ｒ＝１２８までで足る。
【００１３】
第８の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、操作者がｒを変更可能であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第８の観点による高次モーメントを用いた画像投影方法では、操作者が次数ｒを変更可能であるため、操作者の所望のコントラストとなる次数ｒを選択できる。
【００１４】
第８の観点では、本発明は、上記構成の高次モーメントを用いた画像投影方法において、２≦ｒ≦１２８であることを特徴とする高次モーメントを用いた画像投影方法を提供する。
上記第８の観点による高次モーメントを用いた画像投影方法では、前記第２の観点による高次モーメントを用いた画像投影方法において説明した理由により、実用上十分な次数ｒの範囲で投影画像のコントラストを選択できる。
【００１５】
第９の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第９の観点による画像処理装置では、前記第１の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１６】
第１０の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１０の観点による画像処理装置では、前記第２の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１７】
第１１の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１１の観点による画像処理装置では、前記第３の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１８】
第１２の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１２の観点による画像処理装置では、前記第４の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００１９】
第１３の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１３の観点による画像処理装置では、前記第５の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２０】
第１４の観点では、本発明は、３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置を提供する。
上記第１４の観点による画像処理装置では、前記第６の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２１】
第１５の観点では、本発明は、上記構成の画像処理装置において、２≦ｒ≦１２８であることを特徴とする画像処理装置を提供する。
上記第１５の観点による画像処理装置では、前記第７の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２２】
第１６の観点では、本発明は、上記構成の画像処理装置において、操作者がｒを指定するための次数指定手段を具備したことを特徴とする画像処理装置を提供する。
上記第１６の観点による画像処理装置では、前記第８の観点による高次モーメントを用いた画像投影方法を好適に実施できる。
【００２３】
【発明の実施の形態】
以下、本発明の実施の形態を説明する。なお、これにより本発明が限定されるものではない。
【００２４】
−第１の実施形態−
図１は、第１の実施形態にかかる医用画像診断装置の構成図である。
この医用画像診断装置１００は、撮像装置１と、画像処理装置２とを具備している。
【００２５】
撮像装置１は、Ｘ線ＣＴ装置やＭＲＩ装置や超音波診断装置などであり、被検体Ｋを撮像して得たデータを画像処理装置２に渡す。
【００２６】
画像処理装置２は、撮像装置１から渡されたデータを基に３次元データを構築し記憶する３次元データ構築・記憶部２ａと、操作者が投影方向を指定するための投影方向指定部２ｂと、操作者が次数ｒを指定するための次数指定部２ｃと、高次モーメントを用いた画像投影演算を行なう投影演算部２ｄと，投影画像をディスプレイ画面に表示する投影画像表示部２ｅとを具備している。
【００２７】
図２は、次数ｒを変更するために操作者が操作するスライド・レバーの外観図である。
スライド・レバーを動かすことにより、次数ｒを「２」から「１２８」まで変更できる。
【００２８】
図３は、画像処理装置２における高次モーメントを用いた画像投影処理を示すフローチャートである。
ステップＳＴ１では、３次元データ構築・記憶部２ａは、撮像装置１から渡されたデータを基に３次元データを構築し、記憶する。
【００２９】
ステップＳＴ２では、投影方向指定部２ｂは、投影方向を指定するために操作者が操作するデバイス（例えば、トラックボール）から投影方向を読み込む。
ステップＳＴ３では、投影演算部２ｄは、投影方向に垂直な投影面を設定する。
【００３０】
ステップＳＴ４では、投影演算部２ｄは、投影面上の一つの画素に着目する。
ステップＳＴ５では、着目した画素に対応する投影軸上のｎ個のデータＶｉを３次元データより取り出す。
ステップＳＴ６では、次数指定部２ｃは、次数ｒを指定するために操作者が操作するデバイス（例えば、図２のスライド・レバー）から次数ｒを読み込む。
【００３１】
ステップＳＴ７では、投影演算部２ｄは、画素値Ｇを次式により算出する。
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
【００３２】
ステップＳＴ８では、投影演算部２ｄは、全画素の画素値Ｇを得るまで、上記ステップＳＴ４〜ＳＴ７を繰り返す。
ステップＳＴ９では、投影画像表示部２ｅは、得られた投影画像をディスプレイ画面に表示する。
ステップＳＴ１０では、操作者が「終了」を指示すれば処理を終了し、「終了」を指示しなければ上記ステップＳＴ２に戻る。
【００３３】
図４，図５は、高次モーメントを用いた画像投影演算の数値例を２次元的に示した説明図である。次数ｒ＝２としている。
図４に、高次モーメントを用いた画像投影により、３次元データＴＤ１を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。投影軸ａ，ｂ，ｃ，ｄは、それぞれ画素Ａ，Ｂ，Ｃ，Ｄに対応する投影軸である。
図５に、高次モーメントを用いた画像投影により、３次元データＴＤ２を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
【００３４】
図６，図７は、最大値投影方法を用いた画像投影演算の数値例を２次元的に示した説明図である。
図６に、最大値投影方法を用いた画像投影により、３次元データＴＤ１を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
図７に、最大値投影方法を用いた画像投影により、３次元データＴＤ２を投影して得られた投影画像の画素Ａ，Ｂ，Ｃ，Ｄの画素値を示す。
【００３５】
高次モーメントを用いた図４と図５とを比較すれば判るように、投影軸上の最大値以外のデータの値も投影画像に反映されている（例えば、投影軸上の最小値が「１０」か「０」かによって画素値が異なっている）。一方、最大値投影方法を用いた図６と図７とを比較すれば判るように、投影軸上の最大値以外のデータの値は投影画像に全く反映されていない（例えば、投影軸上の最小値が「１０」か「０」かによらず画素値は最大値「７０」になっている）。
このことは、投影方向に例えば骨だけが存在しているのか骨と血管が重なって存在しているのかを、最大値投影方法では投影画像上で判別できないが、本発明では投影画像上で判別できることを意味している。
【００３６】
また、図４と図６とを比較すれば判るように、高次モーメントを用いた図４では、最大値を持つデータがただ１点だけなのか複数点あるのかの情報が投影画像に反映されている（投影軸上に最大値「７０」が存在する個数によって画素値が異なっている）。一方、最大値投影方法を用いた図６では、最大値を持つデータがただ１点だけなのか複数点あるのかの情報が投影画像に全く反映されていない（投影軸上に最大値「７０」が存在する個数によらず画素値は「７０」になっている）。
図５と図７とを比較しても、同じことが判る。
このことは、投影方向に例えば骨が１本あるのか複数本重なって存在しているのかを、最大値投影方法では投影画像上で判別できないが、本発明では投影画像上で判別できることを意味している。
【００３７】
−第２の実施形態−
図３のステップＳＴ７において、投影演算部２ｄは、画素値Ｇを次式により算出してもよい。
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映できる。
【００３８】
−第３の実施形態−
図８は、画像処理装置２における高次モーメントを用いた画像投影処理を示すフロー図である。
このフロー図は、図３のステップＳＴ７がステップＳＴ７’になっている以外は図３のフロー図と同じである。そこで、ステップＳＴ７’のみを説明する。
ステップＳＴ７’において、投影演算部２ｄは、画素値Ｇを次式により算出する。
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映でき、第１の実施形態よりも結果が良い場合があった。
【００３９】
−第４の実施形態−
図８のステップＳＴ７’において、投影演算部２ｄは、画素値Ｇを次式により算出してもよい。
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
この画素値Ｇでも、投影軸上の全てのデータの値を投影画像に反映でき、第２の実施形態よりも結果が良い場合があった。
【００４０】
なお、第１の実施形態〜第４の実施形態のいずれの画素値Ｇとするか、その他の関数Ｇ（Ｐ）とするかは、画像を作成する目的や好みに応じて適宜選べばよい。
【００４１】
【発明の効果】
本発明の高次モーメントを用いた画像投影方法および画像処理装置によれば、３次元データから作成した投影画像に、投影軸上の全てのデータの値を反映することが出来る。
【図面の簡単な説明】
【図１】本発明にかかる医用画像診断装置の構成図である。
【図２】次数ｒを操作者が変更するために操作するスライド・レバーの外観図である。
【図３】第１の実施形態にかかる高次モーメントを用いた画像投影処理を示すフローチャートである。
【図４】高次モーメントを用いた画像投影演算の数値例を２次元的に示した説明図である。
【図５】高次モーメントを用いた画像投影演算の数値例を２次元的に示した別の説明図である。
【図６】最大値投影方法を用いた画像投影演算の数値例を２次元的に示した説明図である。
【図７】最大値投影方法を用いた画像投影演算の数値例を２次元的に示した別の説明図である。
【図８】第３の実施形態にかかる高次モーメントを用いた画像投影処理を示すフローチャートである。
【符号の説明】
１ 撮像装置
２ 画像処理装置
２ａ ３次元データ構築・記憶部
２ｂ 投影方向指定部
２ｃ 次数指定部
２ｄ 投影演算部
２ｅ 投影画像表示部
１００ 医用画像診断装置
Claims (16)
- ３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法。 - ３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法。 - ３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法。 - ３次元データを投影面に向かって投影するとき、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決めることを特徴とする高次モーメントを用いた画像投影方法。 - 投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とすることを特徴とする高次モーメントを用いた画像投影方法。 - 投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とすることを特徴とする高次モーメントを用いた画像投影方法。 - 請求項１から請求項６のいずれかに記載の高次モーメントを用いた画像投影方法において、２≦ｒ≦１２８であることを特徴とする高次モーメントを用いた画像投影方法。
- 請求項１から請求項７のいずれかに記載の高次モーメントを用いた画像投影方法において、操作者がｒを変更可能であることを特徴とする高次モーメントを用いた画像投影方法。
- ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉ／ｎ）ｒ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値を、
Ｐ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
に基づいて決める高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - ３次元データを記憶する３次元データ記憶手段と、操作者が投影方向を指定するための投影方向指定手段と、投影軸上での３次元データの個数をｎとし、データ値をＶｉとし、１より大きい実数をｒとするとき、前記投影軸と投影面の交点の画素値Ｇを、
Ｇ＝ｅｘｐ｛｜（ｉ＝１ΣｎＶｉ／ｎ）ｒ−ｉ＝１Σｎ（Ｖｉｒ）／ｎ｜１／ｒ｝
とする高次モーメントを用いた画像投影手段と、投影画像を表示する投影画像表示手段とを具備したことを特徴とする画像処理装置。 - 請求項９から請求項１４のいずれかに記載の画像処理装置において、２≦ｒ≦１２８であることを特徴とする画像処理装置。
- 請求項９から請求項１５のいずれかに記載の画像処理装置において、操作者がｒを指定するための次数指定手段を具備したことを特徴とする画像処理装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002344974A JP3968003B2 (ja) | 2002-11-28 | 2002-11-28 | 画像処理装置 |
| CNB2003101164369A CN1297233C (zh) | 2002-11-28 | 2003-11-21 | 高阶基于矩的图像投影方法和图像处理装置 |
| US10/723,972 US7218765B2 (en) | 2002-11-28 | 2003-11-26 | Higher-order moment-based image projection method and image processing apparatus |
| EP03257442A EP1434173B1 (en) | 2002-11-28 | 2003-11-26 | Higher-order moment-based image projection method and image processing apparatus |
| DE60311633T DE60311633T2 (de) | 2002-11-28 | 2003-11-26 | Verfahren und Vorrichtung zur Bildprojektion, basierend auf statistischen Momenten höherer Ordnung |
| KR1020030084803A KR20040047681A (ko) | 2002-11-28 | 2003-11-27 | 고차 모멘트 기반 화상 투영 방법 및 화상 처리 장치 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002344974A JP3968003B2 (ja) | 2002-11-28 | 2002-11-28 | 画像処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004178344A true JP2004178344A (ja) | 2004-06-24 |
| JP3968003B2 JP3968003B2 (ja) | 2007-08-29 |
Family
ID=32462819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002344974A Expired - Fee Related JP3968003B2 (ja) | 2002-11-28 | 2002-11-28 | 画像処理装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US7218765B2 (ja) |
| EP (1) | EP1434173B1 (ja) |
| JP (1) | JP3968003B2 (ja) |
| KR (1) | KR20040047681A (ja) |
| CN (1) | CN1297233C (ja) |
| DE (1) | DE60311633T2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8781214B2 (en) * | 2009-10-29 | 2014-07-15 | Optovue, Inc. | Enhanced imaging for optical coherence tomography |
| GB201117807D0 (en) * | 2011-10-14 | 2011-11-30 | Siemens Medical Solutions | Identifying hotspots hidden on mip |
| JP6293619B2 (ja) * | 2014-08-28 | 2018-03-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像処理方法および装置並びにプログラム |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2662525B1 (fr) | 1990-05-25 | 1992-08-28 | Gen Electric Cgr | Procede de visualisation d'une partie de l'image d'une structure physique. |
| US5233299A (en) * | 1991-03-25 | 1993-08-03 | General Electric Company | Projection methods for producing two-dimensional images from three-dimensional data |
| US5297551A (en) * | 1992-08-06 | 1994-03-29 | Picker International, Inc. | Weighted ray projection imaging for MR angiography |
| US6123733A (en) * | 1996-11-27 | 2000-09-26 | Voxel, Inc. | Method and apparatus for rapidly evaluating digital data processing parameters |
| US6421413B1 (en) * | 2000-11-27 | 2002-07-16 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for interactively displaying curved reformation images |
| US6643533B2 (en) | 2000-11-28 | 2003-11-04 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for displaying images of tubular structures |
| US6865246B2 (en) * | 2001-09-26 | 2005-03-08 | Massachusetts Institute Of Technology | True 3D cone-beam imaging method and apparatus |
| US6574297B2 (en) * | 2001-10-30 | 2003-06-03 | Siemens Corporate Research, Inc. | System and method for image reconstruction in a cone beam imaging system |
| US6445762B1 (en) * | 2001-11-26 | 2002-09-03 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for defining regions of interest |
-
2002
- 2002-11-28 JP JP2002344974A patent/JP3968003B2/ja not_active Expired - Fee Related
-
2003
- 2003-11-21 CN CNB2003101164369A patent/CN1297233C/zh not_active Expired - Fee Related
- 2003-11-26 EP EP03257442A patent/EP1434173B1/en not_active Expired - Lifetime
- 2003-11-26 DE DE60311633T patent/DE60311633T2/de not_active Expired - Lifetime
- 2003-11-26 US US10/723,972 patent/US7218765B2/en not_active Expired - Fee Related
- 2003-11-27 KR KR1020030084803A patent/KR20040047681A/ko not_active Application Discontinuation
Also Published As
| Publication number | Publication date |
|---|---|
| EP1434173A1 (en) | 2004-06-30 |
| KR20040047681A (ko) | 2004-06-05 |
| CN1504168A (zh) | 2004-06-16 |
| US7218765B2 (en) | 2007-05-15 |
| DE60311633D1 (de) | 2007-03-22 |
| US20040175025A1 (en) | 2004-09-09 |
| DE60311633T2 (de) | 2007-10-25 |
| CN1297233C (zh) | 2007-01-31 |
| JP3968003B2 (ja) | 2007-08-29 |
| EP1434173B1 (en) | 2007-02-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7518619B2 (en) | Method and apparatus for integrating three-dimensional and two-dimensional monitors with medical diagnostic imaging workstations | |
| US7154985B2 (en) | Method and system for simulating X-ray images | |
| JP4864397B2 (ja) | 対象領域情報を用いて連続的に多重解像度３次元画像を再構成する方法およびシステム | |
| US7397475B2 (en) | Interactive atlas extracted from volume data | |
| US7463262B2 (en) | Image processing apparatus and method | |
| US20080167551A1 (en) | Feature emphasis and contextual cutaways for image visualization | |
| US8244018B2 (en) | Visualizing a 3D volume dataset of an image at any position or orientation from within or outside | |
| JP2005182831A (ja) | 三次元データを視覚化する方法及びシステム | |
| JP2003242257A (ja) | 患者データの処理システム及び処理方法 | |
| KR20210146283A (ko) | 부분 깊이 맵으로부터의 합성 삼차원 이미징의 생성 | |
| JP6739897B2 (ja) | レンダリングを行う方法及びシステム | |
| JP6114266B2 (ja) | 画像をズームするシステム及び方法 | |
| JP2006122663A (ja) | 画像処理装置及び画像処理方法 | |
| JP2006075216A (ja) | 医用画像処理システム、医用画像処理プログラムおよび医用画像処理方法 | |
| US9237849B1 (en) | Relative maximum intensity projection | |
| JP4327171B2 (ja) | 最大値投影方法および装置 | |
| JP3968003B2 (ja) | 画像処理装置 | |
| US8957891B2 (en) | Anatomy-defined automated image generation | |
| US20060103678A1 (en) | Method and system for interactive visualization of locally oriented structures | |
| CA2566635A1 (en) | Method and system for simulating x-ray images | |
| Young et al. | Three-dimensional computed tomography for planning urologic surgery | |
| Solberg et al. | Integration of a real-time video grabber component with the open source image-guided surgery toolkit IGSTK | |
| WO2024110640A1 (en) | Method for generating a simulated bidimensional image of a portion of a patient's body | |
| JP2023132107A (ja) | 画像処理装置、画像処理方法及びプログラム | |
| JP2005058327A (ja) | 投影画像作成表示方法及びこの方法を用いる手術支援装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20040623 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070206 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070426 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070522 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070601 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100608 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100608 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110608 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110608 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120608 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120608 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130608 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130608 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130608 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |