JP6495945B2 - 微分位相コントラストイメージングにおけるスペクトル位相アンラッピングのためのバイアスフリー正則化 - Google Patents
微分位相コントラストイメージングにおけるスペクトル位相アンラッピングのためのバイアスフリー正則化 Download PDFInfo
- Publication number
- JP6495945B2 JP6495945B2 JP2016570099A JP2016570099A JP6495945B2 JP 6495945 B2 JP6495945 B2 JP 6495945B2 JP 2016570099 A JP2016570099 A JP 2016570099A JP 2016570099 A JP2016570099 A JP 2016570099A JP 6495945 B2 JP6495945 B2 JP 6495945B2
- Authority
- JP
- Japan
- Prior art keywords
- objective function
- phase shift
- critical point
- regularization
- shift data
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title description 17
- 230000003595 spectral effect Effects 0.000 title 1
- 230000010363 phase shift Effects 0.000 claims description 74
- 238000000034 method Methods 0.000 claims description 51
- 238000012545 processing Methods 0.000 claims description 31
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 138
- 238000005457 optimization Methods 0.000 description 24
- 238000004422 calculation algorithm Methods 0.000 description 21
- 238000013459 approach Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 6
- 238000007476 Maximum Likelihood Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 3
- 230000002074 deregulated effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000003936 working memory Effects 0.000 description 2
- PEQMJVGRHNZPAM-UHFFFAOYSA-N 1,4-dichloro-2-isocyanatobenzene Chemical compound ClC1=CC=C(Cl)C(N=C=O)=C1 PEQMJVGRHNZPAM-UHFFFAOYSA-N 0.000 description 1
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000003831 deregulation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/48—Diagnostic techniques
- A61B6/484—Diagnostic techniques involving phase contrast X-ray imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
−検出器が測定した位相シフトデータを受け取るステップと、
−正則化項を有しない非凸状又は非凹状の目的関数に、パラメータとして、前記位相シフトデータを含めるステップと、
−前記目的関数の臨界点を見つけるステップと、
−前記臨界点のみの中から、正則化目的関数のターゲット臨界点を見つけるステップであって、前記正則化目的関数は前記目的関数と正則化項とにより校正されるステップと、
−アンラッピングされた位相シフトデータの推定値として、前記ターゲット臨界点を出力するステップとを含む。
−検出器が測定した位相シフトデータを受け取るステップと、
−非凸状又は非凹状の目的関数と正則化項とにより構成された正則化目的関数に、パラメータとして、前記位相シフトデータを含めるステップと、
−前記正則化目的関数の臨界点を見つけるステップであって、前記臨界点はリード臨界点を含むステップと、
−前記正則化項を有しない前記目的関数の非正則化臨界点を見つけるステップと、
−見つかった前記非正則化臨界点のうちから、少なくとも1つの他の非正則化臨界点のうち少なくとも1つより前記リード臨界点に近いターゲットポイントを選択するステップと、
−アンラッピングされた位相シフトデータの推定値として、選択された前記ターゲット点を出力するステップとを含む。
wは異なるエネルギービンを示し、
φmは測定された、潜在的にラッピングされた位相シフトであり、
Mは推定されたアンラッピングされた位相シフト
(外1)
であり、
Ewはエネルギービンwの有効エネルギーであり、
(外2)
は測定された位相シフトφmの標準偏差であり、
Rは正則化関数であり、
λは正則化関数Rの強さを制御する正則化パラメータである。
(外4)
を推定するのに使えるミーゼス分布(von Mises distribution)の変形例を説明している。好適な密度のその他のオプションは、H Gudbjartssonその他により、「The Rician Distribution of Noisy MRI Data」(Magn. Reson. Med., 1995, Dec, 34(6), pp 910−914)の式(5)に記載されている。
Claims (16)
- 位相シフトデータをアンラッピングする方法であって、
検出器が測定した位相シフトデータを受け取るステップと、
正則化項を有しない非凸状又は非凹状の目的関数に、パラメータとして、前記位相シフトデータを含めるステップと、
前記目的関数の臨界点を見つけるステップと、
前記臨界点のみの中から、正則化目的関数のターゲット臨界点を見つけるステップであって、前記正則化目的関数は前記目的関数と正則化項とにより校正されるステップと、
アンラッピングされた位相シフトデータの推定値として、前記ターゲット臨界点を出力するステップとを含む、方法。 - 位相シフトデータをアンラッピングする方法であって、
検出器が測定した位相シフトデータを受け取るステップと、
非凸状又は非凹状の目的関数と正則化項とにより構成された正則化目的関数に、パラメータとして、前記位相シフトデータを含めるステップと、
前記正則化目的関数の臨界点を見つけるステップであって、前記臨界点はリード臨界点を含むステップと、
前記正則化項を有しない前記目的関数の非正則化臨界点を見つけるステップと、
見つかった前記非正則化臨界点のうちから、少なくとも1つの他の非正則化臨界点のうち少なくとも1つより前記リード臨界点に近いターゲット臨界点を選択するステップと、
アンラッピングされた位相シフトデータの推定値として、選択された前記ターゲット臨界点を出力するステップとを含む、方法。 - 前記目的関数は、測定された前記位相シフトデータに関する尤度関数から得られる、
請求項1または2に記載の方法。 - 前記尤度関数の土台となる確率密度はラッピングされたガウシアン密度である、
請求項3に記載の方法。 - 前記尤度関数の土台となる確率密度はミーゼス型密度である、
請求項3に記載の方法。 - 前記正則化項はチホノフ正則化である、請求項1乃至5いずれか一項に記載の方法。
- 前記方法において、あるピクセルにおいて測定された位相シフトデータは、隣接するピクセルにおいて測定された位相シフトデータとは独立して処理される、
請求項1乃至6いずれか一項に記載の方法。 - 位相シフトデータをアンラッピングする信号処理装置であって、
検出器が測定した位相シフトデータを受け取る手段と、
正則化項を有しない非凸状又は非凹状の目的関数に、パラメータとして、前記位相シフトデータを含める手段と、
前記目的関数の臨界点を見つける手段と、
前記臨界点のみの中から、正則化目的関数のターゲット臨界点を見つけるステップであって、前記正則化目的関数は前記目的関数と正則化項とにより校正される手段と、
アンラッピングされた位相シフトデータの推定値として、前記ターゲット臨界点を出力する手段とを有する、
信号処理装置。 - 位相シフトデータをアンラッピングする信号処理装置であって、
検出器が測定した位相シフトデータを受け取る手段と、
非凸状又は非凹状の目的関数と正則化項とにより構成された正則化目的関数に、パラメータとして、前記位相シフトデータを含める手段と、
前記正則化目的関数の臨界点を見つける手段であって、前記臨界点はリード臨界点を含む手段と、
前記正則化項を有しない前記目的関数の非正則化臨界点を見つける手段と、
見つかった前記非正則化臨界点のうちから、少なくとも1つの他の非正則化臨界点のうち少なくとも1つより前記リード臨界点に近いターゲット臨界点を選択する手段と、
アンラッピングされた位相シフトデータの推定値として、選択された前記ターゲット臨界点を出力する手段とを有する、
信号処理装置。 - 前記目的関数を、測定された前記位相シフトデータに関する尤度関数から求めるように構成された、請求項8または9に記載の信号処理装置。
- 前記尤度関数の土台となる確率密度はラッピングされたガウシアン密度である、請求項10に記載の信号処理装置。
- 前記正則化項としてチホノフ正則化を用いるように構成された、
請求項8乃至11いずれか一項に記載の信号処理装置。 - あるピクセルにおいて測定された位相シフトデータを、隣接するピクセルで測定された位相シフトデータとは独立して処理するように構成された、
請求項8乃至12いずれか一項に記載の信号処理装置。 - 前記検出器はエネルギー分解的である、
請求項8乃至13いずれか一項に記載の信号処理装置。 - 請求項8乃至14いずれか一項に記載の信号処理装置を制御するコンピュータプログラム要素であって、データ処理ユニットにより実行されると、請求項1乃至7いずれか一項に記載の方法の各ステップを実行させるコンピュータプログラム要素。
- 請求項15に記載のプログラム要素を記憶したコンピュータ読み取り可能媒体。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14170740.6 | 2014-06-02 | ||
EP14170740 | 2014-06-02 | ||
PCT/EP2015/058566 WO2015185259A1 (en) | 2014-06-02 | 2015-04-21 | Biais-free regularization for spectral phase-unwrapping in differential phase contrast imaging. |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017520294A JP2017520294A (ja) | 2017-07-27 |
JP6495945B2 true JP6495945B2 (ja) | 2019-04-03 |
Family
ID=50943073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016570099A Expired - Fee Related JP6495945B2 (ja) | 2014-06-02 | 2015-04-21 | 微分位相コントラストイメージングにおけるスペクトル位相アンラッピングのためのバイアスフリー正則化 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10037600B2 (ja) |
EP (1) | EP2953097B1 (ja) |
JP (1) | JP6495945B2 (ja) |
CN (1) | CN106462984A (ja) |
WO (1) | WO2015185259A1 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140111818A (ko) * | 2013-03-12 | 2014-09-22 | 삼성전자주식회사 | 엑스선 영상 장치 및 그 제어 방법 |
JP6838531B2 (ja) * | 2017-09-06 | 2021-03-03 | 株式会社島津製作所 | 放射線位相差撮影装置 |
WO2019056309A1 (en) | 2017-09-22 | 2019-03-28 | Shenzhen United Imaging Healthcare Co., Ltd. | METHOD AND SYSTEM FOR GENERATING A PHASE CONTRAST IMAGE |
CN108537842B (zh) * | 2017-12-29 | 2022-01-07 | 南京理工大学 | 差分相衬显微成像中背景非均匀性的校正与补偿方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6741357B2 (en) * | 2001-08-14 | 2004-05-25 | Seagate Technology Llc | Quadrature phase shift interferometer with unwrapping of phase |
US7372393B2 (en) * | 2006-07-07 | 2008-05-13 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images |
US8433739B2 (en) * | 2007-06-27 | 2013-04-30 | L-3 Communications Integrated Systems, L.P. | Methods and systems for detecting repetitive synchronized signal events |
US8855265B2 (en) | 2009-06-16 | 2014-10-07 | Koninklijke Philips N.V. | Correction method for differential phase contrast imaging |
WO2012038857A1 (en) | 2010-09-20 | 2012-03-29 | Koninklijke Philips Electronics N.V. | Phase gradient unwrapping in differential phase contrast imaging |
KR101795837B1 (ko) * | 2011-01-18 | 2017-11-08 | 사반트 시스템즈 엘엘씨 | 헤드-업 동작 및 비주얼 피드백을 제공하는 원격 제어 인터페이스 |
CN102217934A (zh) * | 2011-04-08 | 2011-10-19 | 中国科学院深圳先进技术研究院 | 磁共振成像方法及*** |
RU2598310C2 (ru) * | 2011-08-31 | 2016-09-20 | Конинклейке Филипс Н.В. | Формирование дифференциальных фазовых контрастных изображений с помощью чувствительного к энергии обнаружения |
FR2983464B1 (fr) * | 2011-12-05 | 2014-06-27 | Solystic | Dispositif d'empilage pour objets plats empiles sur chant et machine de tri postal |
WO2013091078A1 (en) * | 2011-12-23 | 2013-06-27 | Liu Junmin | Method for phase unwrapping |
JP6177800B2 (ja) * | 2012-01-12 | 2017-08-09 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | X線システムにおける減衰画像データおよび位相画像データの生成 |
CN103632345B (zh) * | 2013-11-27 | 2016-03-30 | 中国科学技术大学 | 一种基于正则化的mri图像非均匀性校正方法 |
-
2015
- 2015-04-21 CN CN201580029369.8A patent/CN106462984A/zh active Pending
- 2015-04-21 EP EP15164391.3A patent/EP2953097B1/en not_active Not-in-force
- 2015-04-21 US US15/309,223 patent/US10037600B2/en active Active
- 2015-04-21 JP JP2016570099A patent/JP6495945B2/ja not_active Expired - Fee Related
- 2015-04-21 WO PCT/EP2015/058566 patent/WO2015185259A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
EP2953097B1 (en) | 2016-10-26 |
JP2017520294A (ja) | 2017-07-27 |
WO2015185259A1 (en) | 2015-12-10 |
EP2953097A1 (en) | 2015-12-09 |
US20170091933A1 (en) | 2017-03-30 |
US10037600B2 (en) | 2018-07-31 |
CN106462984A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6150940B2 (ja) | 位相コントラストctを使うことによる単色減衰コントラスト画像生成 | |
CN107810524B (zh) | 针对暗场和相衬ct的鲁棒重建 | |
JP5997180B2 (ja) | デュアルエネルギーct画像再構成のための方法及びシステム | |
JP6214819B1 (ja) | 微分位相コントラストx線撮像における暗視野信号の最適なエネルギ加重 | |
JP6495945B2 (ja) | 微分位相コントラストイメージングにおけるスペクトル位相アンラッピングのためのバイアスフリー正則化 | |
US9613441B2 (en) | Joint reconstruction of electron density images | |
JP6105586B2 (ja) | エネルギー高感度検出の微分位相コントラストイメージング | |
KR101892321B1 (ko) | 타겟 대상의 영역의 이미지를 생성하기 위한 이미지 데이터 제공 방법 및 장치 | |
Köhler et al. | Noise properties of grating‐based x‐ray phase contrast computed tomography | |
JP6805173B2 (ja) | 暗視野・位相コントラストイメージングをスキャンするためのビーム硬化補正 | |
JP5739976B2 (ja) | 位相画像を再構成する方法、コンピュータプログラム、位相画像を再構成する演算装置 | |
EP3170148B1 (en) | Iterative reconstruction method for spectral, phase-contrast imaging | |
US20140169524A1 (en) | Image reconstruction method for differential phase contrast x-ray imaging | |
WO2012029039A1 (en) | Beam hardening correction for phase-contrast imaging | |
WO2013073453A1 (en) | Imaging apparatus and image processing method | |
EP3510563A2 (en) | Improved phase-contrast and dark-field ct reconstruction algorithm | |
WO2016091594A1 (en) | Statistically weighted regularization in multi-contrast imaging | |
US9330456B2 (en) | Systems and methods for regularized Fourier analysis in x-ray phase contrast imaging | |
Wang et al. | Graphics processing unit-based dispersion encoded full-range frequency-domain optical coherence tomography | |
US20230153968A1 (en) | Bias correction for dark-field imaging based on sliding window phase retrieval | |
TWI500412B (zh) | 一種影像品質的改善處理方法及其造影系統 | |
KR20200017188A (ko) | 기계 학습을 통한 격자간섭계에서의 민감도 및 해상도 향상을 위한 장치 | |
Simo et al. | Non-uniform projection angle processing in computed tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180419 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20180419 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20180531 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180619 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180910 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181023 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190115 |
|
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: 20190212 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190307 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6495945 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |