[go: up one dir, main page]

CN111643100B - Information Characterization Method and System for Phase Contrast Imaging System - Google Patents

Information Characterization Method and System for Phase Contrast Imaging System Download PDF

Info

Publication number
CN111643100B
CN111643100B CN201911147099.7A CN201911147099A CN111643100B CN 111643100 B CN111643100 B CN 111643100B CN 201911147099 A CN201911147099 A CN 201911147099A CN 111643100 B CN111643100 B CN 111643100B
Authority
CN
China
Prior art keywords
information
phase contrast
dark field
phase
contrast
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.)
Active
Application number
CN201911147099.7A
Other languages
Chinese (zh)
Other versions
CN111643100A (en
Inventor
张丽
吴承鹏
高河伟
邢宇翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201911147099.7A priority Critical patent/CN111643100B/en
Publication of CN111643100A publication Critical patent/CN111643100A/en
Application granted granted Critical
Publication of CN111643100B publication Critical patent/CN111643100B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/484Diagnostic techniques involving phase contrast X-ray imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种相衬成像信息表征方法,包括:获取相衬成像的原始数据;根据所述原始数据对相衬信息和暗场信息进行提取;对提取的所述相衬信息和暗场信息通过公式:H=D′·sin(P′)进行表征,其中H为相衬信息和暗场信息的混合场对比度,P'为对所述相衬信息的对应数据调整到[﹣π,π]内的相位偏移,D'对所述暗场信息的对应数据进行指数化操作后的图像对比度降低比值。通过对提取的相衬信息和暗场信息进行处理,能够将原有的相衬信息和暗场信息表征在一张混合场图像中,该混合场图像与原有的相衬图像和暗场图像相比,图像表征清楚,相比于普通的图像融合方法具有更加明确的物理含义。

Figure 201911147099

The invention discloses a method for characterizing phase contrast imaging information, comprising: acquiring original data of phase contrast imaging; extracting phase contrast information and dark field information according to the original data; and extracting the extracted phase contrast information and dark field information The information is characterized by the formula: H=D'·sin(P'), where H is the mixed field contrast of the phase contrast information and the dark field information, and P' is the corresponding data of the phase contrast information adjusted to [﹣π, The phase shift within π], D', the image contrast reduction ratio after performing the indexing operation on the corresponding data of the dark field information. By processing the extracted phase contrast information and dark field information, the original phase contrast information and dark field information can be represented in a mixed field image, which is different from the original phase contrast image and dark field image. In contrast, the image representation is clear, and it has a clearer physical meaning than ordinary image fusion methods.

Figure 201911147099

Description

相衬成像系统信息表征方法及系统Information Characterization Method and System for Phase Contrast Imaging System

技术领域technical field

本发明涉及辐射成像领域,特别涉及一种相衬成像信息表征方法及系统。The invention relates to the field of radiation imaging, in particular to a method and system for characterizing phase contrast imaging information.

背景技术Background technique

相衬成像技术主要采用X射线对待测物进行照射,能够实现微米甚至亚微米级上的局部结构分辨,是对传统X射线衰减成像技术的很好补充。该技术能够同时表征吸收、相衬和暗场三种对比度信息,适用于低原子序数、低密度物质,特别是针对包括乳腺在内的生物软组织结构。目前X射线相衬成像技术主要是表征出三种对比度信息,包括吸收信息、相衬信息和暗场信息。这三种对比度信息分别可以表征出X射线在物体中的衰减量、相位偏移量以及小角散射量,这几种表征量目前均已有比较明确的物理含义,如衰减量对应于物体的线性衰减系数、相位偏移量对应于物体的复折射率实部、小角散射量对应于物体的广义散射参数。Phase contrast imaging technology mainly uses X-rays to irradiate the object to be tested, which can achieve local structure resolution at the micron or even sub-micron level, which is a good complement to the traditional X-ray attenuation imaging technology. This technique can simultaneously characterize three kinds of contrast information of absorption, phase contrast and dark field, and is suitable for low atomic number, low density substances, especially for biological soft tissue structures including breast. At present, X-ray phase contrast imaging technology mainly characterizes three kinds of contrast information, including absorption information, phase contrast information and dark field information. These three kinds of contrast information can respectively characterize the attenuation, phase shift and small-angle scattering of X-rays in the object. These kinds of characterization quantities have clear physical meanings. For example, the attenuation corresponds to the linearity of the object. The attenuation coefficient and the phase offset correspond to the real part of the complex refractive index of the object, and the small-angle scattering amount corresponds to the generalized scattering parameter of the object.

在实际应用中(例如医学成像应用),对于观察图像的用户(如放射物理师)而言,如果能够将所有最相关的信息尽可能精简地呈现出来将具有非常重要的实际价值。因此对于X射线相衬成像技术而言,这也就正是其中的一个重要研究方向,即选取一种合适的表征量来尽可能整合以上所述的三种对比度信息,尤其是整合相衬成像特有的相衬信息和暗场信息。在现有的技术和研究中,传统技术手段是通过有两种方式向用户反馈吸收、相衬、暗场三种信息,第一是将三种信息通过三个不同的吸收图像、相衬图像、暗场图像给用户,比较准确,但用户需要同时查看三个图像,非常不方便;第二是得到方式一种的三个不同的图像,对其中两个图像进行图像处理操作之后得到一个图像,可以向用户同时反映三种信息,即通过表征出的几种对比度信息图像进行后处理的图像融合,这些融合方法得到的图像虽然也能反映出多种信息,但是由于是通过图像处理获得的,并没有明确的物理含义,而且图像处理算法中的一些缺陷,并不能解释得到的单一图像中与原有信息表征的三种对比度图像之间的联系,即会造成最终获得的单个图像的三种信息无法明确对应实际的物质信息。In practical applications (such as medical imaging applications), it is of great practical value for users viewing images (such as radiophysicists) to be able to present all the most relevant information as condensed as possible. Therefore, for X-ray phase contrast imaging technology, this is an important research direction, that is, to select an appropriate characterization quantity to integrate the above three contrast information as much as possible, especially the integration of phase contrast imaging Unique phase contrast information and dark field information. In the existing technology and research, the traditional technical means is to feed back three kinds of information of absorption, phase contrast and dark field to the user in two ways. The first is to pass the three kinds of information through three different absorption images and phase contrast images. , Dark field images are more accurate for users, but users need to view three images at the same time, which is very inconvenient; the second is to obtain three different images in one way, and one image is obtained after image processing operation on two of the images. , which can reflect three kinds of information to the user at the same time, that is, image fusion through post-processing of several contrast information images represented. Although the images obtained by these fusion methods can also reflect a variety of information, they are obtained through image processing. , there is no clear physical meaning, and some defects in the image processing algorithm cannot explain the connection between the obtained single image and the three contrast images represented by the original information, that is, it will cause the three contrast images of the single image finally obtained. Such information cannot clearly correspond to actual substance information.

为了解决以上问题,本发明提出了一种相衬成像信息表征方法及系统,能够实现将原有的相衬信息和暗场信息直接表征在一张图像中,并且具有明确的物理定义。In order to solve the above problems, the present invention proposes a phase contrast imaging information representation method and system, which can realize the direct representation of the original phase contrast information and dark field information in an image, and has a clear physical definition.

发明内容SUMMARY OF THE INVENTION

(一)要解决的技术问题(1) Technical problems to be solved

为使得用户可以采用同一图像同时查看相衬和暗场两种信息,同时可以保证图像信息清晰准确,本发明提出了一种相衬成像信息表征方法及系统、可读存储介质,具体如下:In order to enable users to use the same image to view two kinds of information of phase contrast and dark field at the same time, and at the same time to ensure that the image information is clear and accurate, the present invention proposes a method and system for characterizing phase contrast imaging information, and a readable storage medium, as follows:

(二)技术方案(2) Technical solutions

本发明的一个方面提出了一种相衬成像信息表征方法,方法包括:获取相衬成像的原始数据;根据原始数据对相衬信息和暗场信息进行提取;对提取的相衬信息和暗场信息通过以下公式进行表征:One aspect of the present invention provides a method for characterizing phase contrast imaging information, the method includes: acquiring original data of phase contrast imaging; extracting phase contrast information and dark field information according to the original data; extracting phase contrast information and dark field information Information is characterized by the following formula:

H=D′·sin(P′)H=D′·sin(P′)

其中,H为相衬信息和暗场信息的混合场对比度,P′为对相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Among them, H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the corresponding data of the dark field information The ratio of image contrast reduction after exponential operation.

可选地,根据原始数据对相衬信息和暗场信息进行提取,是通过包括FCA算法、SAXS算法或ASAXS算法实现的。Optionally, extracting the phase contrast information and the dark field information according to the original data is realized by including the FCA algorithm, the SAXS algorithm or the AAXS algorithm.

可选地,根据原始数据对相衬信息和暗场信息进行提取之后,还包括:对相衬信息的对应数据进行归一化操作获取调整到[-π,π]内的相位偏移P′;对暗场信息的对应数据进行指数化操作获取的图像对比度降低比值D′。Optionally, after extracting the phase contrast information and the dark field information according to the original data, the method further includes: performing a normalization operation on the corresponding data of the phase contrast information to obtain a phase offset P' adjusted to [-π, π] ; The image contrast reduction ratio D' obtained by performing the indexing operation on the corresponding data of the dark field information.

本发明的另一个方面提出了一种相衬成像信息表征系统,应用上述的方法,系统包括:光学成像装置,用于获取相衬成像的原始数据;信息提取装置,用于根据光学成像装置获得的原始数据对相衬信息和暗场信息进行提取;信息表征装置,用于根据信息提取装置提取的相衬信息和暗场信息通过公式:Another aspect of the present invention provides a phase contrast imaging information characterization system, applying the above method, the system includes: an optical imaging device for acquiring original data of phase contrast imaging; an information extraction device for obtaining according to the optical imaging device Extract the phase contrast information and dark field information from the original data of

H=D′·sin(P′)H=D′·sin(P′)

进行表征,其中H为相衬信息和暗场信息的混合场对比度,P′为对相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Characterization, where H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the contrast of the dark field information The ratio of image contrast reduction after indexing the corresponding data.

可选地,光学成像装置为分析成像装置、边缘照射成像装置、光斑成像装置或光栅成像装置。Optionally, the optical imaging device is an analytical imaging device, an edge illumination imaging device, a spot imaging device or a grating imaging device.

可选地,光栅成像装置为Talbot-Lau类型、几何投影类型或双相位光栅类型。Optionally, the grating imaging device is a Talbot-Lau type, a geometric projection type or a dual phase grating type.

可选地,光栅成像装置的照射待测位置位于光源与第一光栅之间或位于第一光栅和第二光栅之间。Optionally, the position to be measured for illumination of the grating imaging device is located between the light source and the first grating or between the first grating and the second grating.

可选地,光学成像装置的光源为常规X射线光源、微焦点X射线光源或同步辐射X射线光源,光学成像装置的探测器为能量积分型探测器或光子计数型探测器。Optionally, the light source of the optical imaging device is a conventional X-ray light source, a microfocus X-ray light source or a synchrotron radiation X-ray light source, and the detector of the optical imaging device is an energy integration detector or a photon counting detector.

本发明的另一个方面提出了一种相衬成像信息表征系统,系统包括:处理器;存储器,其上存储有可在处理器上运行的计算机程序,其中计算机程序被处理器执行时,实现上述的方法。Another aspect of the present invention provides a system for characterizing phase contrast imaging information. The system includes: a processor; Methods.

本发明的另一个方面提出了一种计算机可读存储介质,计算机可读存储介质上存储有数据处理程序,数据处理程序被处理器执行时,实现上述的方法。Another aspect of the present invention provides a computer-readable storage medium. The computer-readable storage medium stores a data processing program. When the data processing program is executed by a processor, the above method is implemented.

(三)有益效果(3) Beneficial effects

本发明提出的一种相衬成像信息表征方法及系统,通过对提取的相衬信息和暗场信息进行处理,能够将原有的相衬信息和暗场信息表征在一张混合场图像中,该混合场图像与原有的相衬图像和暗场图像相比,图像表征清楚,相比于普通的图像融合方法具有更加明确的物理含义。The method and system for characterizing phase contrast imaging information proposed by the present invention can characterize the original phase contrast information and dark field information in a mixed field image by processing the extracted phase contrast information and dark field information. Compared with the original phase-contrast image and dark-field image, the mixed-field image has a clear image representation, and has a clearer physical meaning than the common image fusion method.

附图说明Description of drawings

图1是本发明实施例中相衬成像信息表征方法流程示意图;1 is a schematic flowchart of a method for characterizing phase contrast imaging information in an embodiment of the present invention;

图2是本发明实施例中分别采用FCA、SAXS、ASAXS三种算法对应提取的吸收、相衬、暗场三种信息的成像对比图;2 is an imaging comparison diagram of three kinds of information of absorption, phase contrast, and dark field respectively extracted by three algorithms of FCA, SAXS, and AAXS in the embodiment of the present invention;

图3是本发明实施例中对ASAXS算法进行混合表征与FCA、SAXS算法分别提取的暗场信息的成像对比图。FIG. 3 is an imaging comparison diagram of the dark field information extracted by the hybrid characterization of the AAXS algorithm and the FCA and SAXS algorithms, respectively, in an embodiment of the present invention.

具体实施方式Detailed ways

为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to specific embodiments and accompanying drawings.

为使得用户可以采用同一图像同时查看相衬和暗场两种信息,同时可以保证图像信息清晰准确,本发明提出了一种相衬成像信息表征方法及系统、可读存储介质,具体如下:In order to enable users to use the same image to view two kinds of information of phase contrast and dark field at the same time, and at the same time to ensure that the image information is clear and accurate, the present invention proposes a method and system for characterizing phase contrast imaging information, and a readable storage medium, as follows:

本发明的一个方面提出了一种相衬成像信息表征方法,如图1是本发明实施例中相衬成像信息表征方法流程示意图所示,方法包括:One aspect of the present invention proposes a method for characterizing phase contrast imaging information. FIG. 1 is a schematic flowchart of the method for characterizing phase contrast imaging information in an embodiment of the present invention. The method includes:

S110、获取相衬成像的原始数据;S110, acquiring raw data of phase contrast imaging;

S120、根据原始数据对相衬信息和暗场信息进行提取;S120, extracting phase contrast information and dark field information according to the original data;

S130、对提取的相衬信息和暗场信息通过如下公式进行表征:S130, characterize the extracted phase contrast information and dark field information by the following formula:

H=D′·sin(P′)H=D′·sin(P′)

其中,H为相衬信息和暗场信息的混合场对比度,P′为对相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Among them, H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the corresponding data of the dark field information The ratio of image contrast reduction after exponential operation.

具体地,在本发明实施例中,采用X射线相衬光栅成像装置作为相衬成像系统,X射线相衬光栅成像装置包括X射线光源,源光栅、相位光栅、分析光栅源光栅及探测器,具体地,可以采用相位光栅间距为1020mm,源光栅到分析光栅间距为1280mm,源光栅、相位光栅和分析光栅周期分别为16.8μm、4.2μm和2.4μm,X射线球管参数为35kVp和35mA,探测器像素为74.8μm×74.8μm,相位步进数为8。Specifically, in the embodiment of the present invention, an X-ray phase contrast grating imaging device is used as the phase contrast imaging system, and the X-ray phase contrast grating imaging device includes an X-ray light source, a source grating, a phase grating, an analysis grating source grating and a detector, Specifically, the phase grating spacing can be 1020mm, the source grating to analysis grating spacing is 1280mm, the source grating, phase grating and analysis grating periods are 16.8μm, 4.2μm and 2.4μm respectively, and the X-ray tube parameters are 35kVp and 35mA, The detector pixel is 74.8 μm × 74.8 μm, and the number of phase steps is 8.

在上述X射线相衬成像系统采集到原始数据后,需要通过信息提取算法对原始数据进行吸收(Absorption)信息、相衬(Phase-contrast)信息、暗场(Dark-field)信息三种信息的提取,以下分别用A、P、D表示。本发明以相衬信息和暗场信息的表征为实施例对本发明进行详细说明。相衬信息可以理解为由于物体受到X射线照射后,X射线经过物体发生折射造成的X射线相位偏移;暗场信息可以理解为由于物体受到X射线照射后,X射线经过物体发生散射造成的X射线图像对比度改变。本发明定义了一种全新的X射线相衬成像物理表征量H,用于将暗场信息和相衬信息表征在同一混合场下,该混合场信息可以单一图像的形式向用户进行展现。其中,该X射线相衬成像物理表征量即“混合场”(Hybrid-field)对比度H,满足如下公式(1):After the above-mentioned X-ray phase contrast imaging system collects the original data, it is necessary to use the information extraction algorithm to extract three kinds of information: absorption information, phase-contrast information, and dark-field information. Extraction, hereinafter represented by A, P, and D, respectively. The present invention will be described in detail by taking the characterization of phase contrast information and dark field information as an embodiment. Phase contrast information can be understood as the X-ray phase shift caused by the refraction of X-rays through the object after the object is irradiated by X-rays; dark field information can be understood as the X-rays scattered by the object after the object is irradiated by X-rays. X-ray image contrast changes. The present invention defines a new X-ray phase contrast imaging physical characterization quantity H, which is used to characterize the dark field information and the phase contrast information in the same mixed field, and the mixed field information can be displayed to the user in the form of a single image. Among them, the X-ray phase contrast imaging physical characterization quantity, that is, the "Hybrid-field" contrast H, satisfies the following formula (1):

H=D′·sin(P′)H=D′·sin(P′)

其中,H为相衬信息和暗场信息的混合场对比度,P′为对相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Among them, H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the corresponding data of the dark field information The ratio of image contrast reduction after exponential operation.

可选地,根据原始数据对相衬信息和暗场信息进行提取,是通过包括FCA算法、SAXS算法或ASAXS算法实现的。Optionally, extracting the phase contrast information and the dark field information according to the original data is realized by including the FCA algorithm, the SAXS algorithm or the AAXS algorithm.

具体地,在本发明实施例中,在基于光栅的X射线相衬成像系统中,根据原始数据对相衬信息和暗场信息进行提取在本发明中至少存在FCA(Fourier Component Analysis)算法、SAXS(Small Angle X-ray Scattering)算法或ASAXS算法三种提取算法,通过这三种提取算法都可以获得原始数据对应的相衬信息和暗场信息,以下对三种提取算法分别进行详细介绍:Specifically, in the embodiment of the present invention, in the grating-based X-ray phase contrast imaging system, the phase contrast information and the dark field information are extracted according to the original data. In the present invention, there are at least FCA (Fourier Component Analysis) algorithm, SAXS algorithm There are three extraction algorithms (Small Angle X-ray Scattering) algorithm or AAXS algorithm. Through these three extraction algorithms, the phase contrast information and dark field information corresponding to the original data can be obtained. The three extraction algorithms are described in detail below:

在FCA算法中,基于光栅的X射线相衬成像系统采集到的背景位移曲线f(φ)和物体位移曲线s(φ)均假设为余弦模型,即公式(2)和公式(3):In the FCA algorithm, both the background displacement curve f(φ) and the object displacement curve s(φ) collected by the grating-based X-ray phase contrast imaging system are assumed to be cosine models, namely formulas (2) and (3):

Figure GDA0003206699890000051
Figure GDA0003206699890000051

Figure GDA0003206699890000052
Figure GDA0003206699890000052

其中,

Figure GDA0003206699890000053
Figure GDA0003206699890000054
分别为背景位移取线和物体位移曲线的零阶、一阶Fourier成分幅度,
Figure GDA0003206699890000055
为两者的一阶Fourier成分幅角。in,
Figure GDA0003206699890000053
and
Figure GDA0003206699890000054
are the zero-order and first-order Fourier component amplitudes of the background displacement line and the object displacement curve, respectively,
Figure GDA0003206699890000055
is the first-order Fourier component argument of both.

FCA算法通过快速傅里叶变换可以计算背景位移曲线f(φ)和物体位移曲线s(φ)的参数,进而获得提取到的吸收信息AFCA、相衬信息PFCA和暗场信息DFCA,具体可以表达为如下形式,即公式(4)、(5)和(6):The FCA algorithm can calculate the parameters of the background displacement curve f(φ) and the object displacement curve s(φ) through the fast Fourier transform, and then obtain the extracted absorption information A FCA , phase contrast information P FCA and dark field information D FCA , Specifically, it can be expressed in the following forms, namely formulas (4), (5) and (6):

Figure GDA0003206699890000056
Figure GDA0003206699890000056

Figure GDA0003206699890000061
Figure GDA0003206699890000061

Figure GDA0003206699890000062
Figure GDA0003206699890000062

其中,A、P、D分别代表吸收信息、相衬信息、暗场信息,p2为光栅相衬成像中靠近探测器的最后一块光栅周期,d为最后两块光栅之间的距离。

Figure GDA0003206699890000063
为两位移曲线的相位差,
Figure GDA0003206699890000064
分别为物体位移曲线和背景位移曲线的对比度。Among them, A, P, and D represent absorption information, phase contrast information, and dark field information, respectively, p 2 is the period of the last grating near the detector in grating phase contrast imaging, and d is the distance between the last two gratings.
Figure GDA0003206699890000063
is the phase difference between the two displacement curves,
Figure GDA0003206699890000064
are the contrast of the object displacement curve and the background displacement curve, respectively.

因此,可见FCA算法基于余弦模型假设,通过快速傅里叶变换即可提取出三种信息,方便高效,但是存在固有的相位卷绕问题(Phase-wrapping),在实际应用中可能对图像质量有较大影响,容易造成图像质量较差。Therefore, it can be seen that the FCA algorithm is based on the cosine model assumption, and can extract three kinds of information through fast Fourier transform, which is convenient and efficient, but there is an inherent phase-wrapping problem (Phase-wrapping), which may affect the image quality in practical applications. Larger impact, it is easy to cause poor image quality.

在SAXS算法中,基于光栅的X射线相衬成像系统基于光栅的X射线相衬成像系统采集到的背景位移曲线和物体位移曲线。其中,基本假设为物体位移曲线s(φ)可以表示为背景位移曲线f(φ)和小角散射分布g(φ)的卷积,其表达式如下,即公式(7):In the SAXS algorithm, the background displacement curve and the object displacement curve collected by the grating-based X-ray phase contrast imaging system are collected. Among them, the basic assumption is that the object displacement curve s(φ) can be expressed as the convolution of the background displacement curve f(φ) and the small-angle scattering distribution g(φ), and its expression is as follows, namely formula (7):

Figure GDA0003206699890000065
Figure GDA0003206699890000065

基于上述(7)式,先通过Lucy-Richardson迭代的方法,从采集到的每个探测器像素上的背景位移曲线f(φ)和物体位移曲线s(φ)中解卷积获得小角散射分布,解卷积过程中的第k次迭代计算公式为公式(8):Based on the above formula (7), the small-angle scattering distribution is obtained by deconvolution from the background displacement curve f(φ) and the object displacement curve s(φ) on each detector pixel collected by the Lucy-Richardson iteration method. , the k-th iteration calculation formula in the deconvolution process is formula (8):

Figure GDA0003206699890000066
Figure GDA0003206699890000066

其中,

Figure GDA0003206699890000067
表示f关于原点的镜像对称,通常选择g0=s作为初始值。然后可以将传统FCA方法提取出的吸收(A)、相衬(P)和暗场(D)三种对比度信息表示为小角散射分布的零阶矩M0(g)、一阶矩M1(g)和二阶矩M2(g),即公式(9)、(10)和(11):in,
Figure GDA0003206699890000067
Indicates the mirror symmetry of f about the origin, and g 0 =s is usually chosen as the initial value. Then the three contrast information of absorption (A), phase contrast (P) and dark field (D) extracted by the traditional FCA method can be expressed as the zero-order moment M 0 (g) and the first-order moment M 1 ( g) and the second moment M 2 (g), namely formulas (9), (10) and (11):

A→M0(g)=∫g(φ)dφA→M 0 (g)=∫g(φ)dφ

P→M1(g)=∫φg(φ)dφ/M0(g)P→M 1 (g)=∫φg(φ)dφ/M 0 (g)

D→M2(g)=∫(φ-M1(g))2g(φ)dφ/M0(g))D→M 2 (g)=∫(φ-M 1 (g)) 2 g(φ)dφ/M 0 (g))

因此,可见SAXS算法基于卷积假设,通过解卷积的方法得到小角散射分布,然后计算其多阶矩即可获得三种信息,避免了FCA的相位卷绕问题,但是由于需要进行解卷积迭代,运算效率较低,同时需要预先设定迭代方法和停止条件,这些差异对迭代结果都会有较大影响。Therefore, it can be seen that the SAXS algorithm is based on the convolution assumption, obtains the small-angle scattering distribution by deconvolution, and then calculates its multi-order moment to obtain three kinds of information, avoiding the phase warping problem of FCA, but due to the need for deconvolution Iteration, the operation efficiency is low, and the iterative method and stopping conditions need to be set in advance, and these differences will have a greater impact on the iterative results.

在ASAXS信息提取算法与SAXS方法一样不依赖于余弦模型假设,其可以直接基于光栅的X射线相衬成像系统采集到的原始数据背景位移曲线f(φ)和物体位移曲线s(φ),对吸收(M0)、相衬(M1)和暗场(M2)三种信息进行提取,具体可以表示为三种多阶矩,即公式(12)、(13)和(14):Like the SAXS method, the AAXS information extraction algorithm does not depend on the assumption of the cosine model. It can be directly based on the background displacement curve f(φ) and object displacement curve s(φ) of the original data collected by the grating X-ray phase contrast imaging system. Absorption (M 0 ), phase contrast (M 1 ) and dark field (M 2 ) three kinds of information are extracted, which can be expressed as three kinds of multi-order moments, namely formulas (12), (13) and (14):

Figure GDA0003206699890000071
Figure GDA0003206699890000071

Figure GDA0003206699890000072
Figure GDA0003206699890000072

Figure GDA0003206699890000073
Figure GDA0003206699890000073

其中,SinM1(·)和CosM1(·)的定义如下,即公式(15)和(16):Among them, SinM 1 (·) and CosM 1 (·) are defined as follows, namely formulas (15) and (16):

SinM1(y)=∫sin(φ)y(φ)dφ,SinM 1 (y)=∫sin(φ)y(φ)dφ,

CosM1(y)=∫cos(φ)y(φ)dφ.CosM 1 (y)=∫cos(φ)y(φ)dφ.

因此,可见ASAXS算法避免了FCA算法的相位卷绕问题,以及SAX算法的解卷积迭代问题,运算效率更高。Therefore, it can be seen that the ASAXS algorithm avoids the phase winding problem of the FCA algorithm and the deconvolution iteration problem of the SAX algorithm, and the operation efficiency is higher.

需要说明的是,以上三种信息提取算法,均可以采用公式(1)的形式,通过对相衬信息和暗场信息进行处理,得到对应的P′、D′,进而得到相衬信息和暗场信息的混合场信息H,即相衬信息和暗场信息的混合场对比度,进而可以将暗场信息和相衬信息表征在同一混合场下,该混合场信息可以单一图像的形式向用户进行展现。It should be noted that the above three information extraction algorithms can all be in the form of formula (1). By processing the phase contrast information and the dark field information, the corresponding P' and D' can be obtained, and then the phase contrast information and dark field information can be obtained. The mixed field information H of the field information is the mixed field contrast of the phase contrast information and the dark field information, so that the dark field information and the phase contrast information can be represented in the same mixed field, and the mixed field information can be displayed to the user in the form of a single image. show.

可选地,根据原始数据对相衬信息和暗场信息进行提取,包括:根据原始数据获取吸收信息;从原始数据和/或吸收信息提取相衬信息和暗场信息。在本发明实施例中,如公式(12)、(13)、(14)所示,根据原始数据获取吸收信息M0,相衬信息M1则根据原始数据和吸收信息M0获取,同时暗场信息M2则根据原始数据和吸收信息M0、相衬信息M1获取。Optionally, extracting the phase contrast information and the dark field information according to the original data includes: acquiring absorption information according to the original data; and extracting the phase contrast information and the dark field information from the original data and/or the absorption information. In the embodiment of the present invention, as shown in formulas (12), (13) and (14), the absorption information M 0 is obtained according to the original data, and the phase contrast information M 1 is obtained according to the original data and the absorption information M 0 . The field information M 2 is obtained according to the original data, absorption information M 0 , and phase contrast information M 1 .

可选地,根据原始数据对相衬信息和暗场信息进行提取之后,还包括:对相衬信息的对应数据进行归一化操作获取调整到[-π,π]内的相位偏移P′;对暗场信息的对应数据进行指数化操作获取的图像对比度降低比值D′。Optionally, after extracting the phase contrast information and the dark field information according to the original data, the method further includes: performing a normalization operation on the corresponding data of the phase contrast information to obtain a phase offset P' adjusted to [-π, π] ; The image contrast reduction ratio D' obtained by performing the indexing operation on the corresponding data of the dark field information.

在本发明实施例中,以ASAXS算法为例,将公式(2)(3)中的余弦模型代入公式(12)(13)(14)中,建立FCA算法与ASAXS算法之间的联系,即可以将ASAXS算法提取出的吸收(M0)、相衬(M1)和暗场(M2)三种信息重新表示为吸收信息(AASAXS)、相衬信息(PASAXS)和暗场信息(DASAXS),即公式(17)、(18)和(19):In the embodiment of the present invention, taking the AAXS algorithm as an example, the cosine model in the formula (2) (3) is substituted into the formula (12) (13) (14) to establish the connection between the FCA algorithm and the AAXS algorithm, that is, The absorption (M 0 ), phase contrast (M 1 ) and dark field (M 2 ) information extracted by the ASAXS algorithm can be re-expressed as absorption information (A ASAXS ), phase contrast information (P ASAXS ) and dark field information (D ASAXS ), i.e. equations (17), (18) and (19):

Figure GDA0003206699890000081
Figure GDA0003206699890000081

Figure GDA0003206699890000082
Figure GDA0003206699890000082

Figure GDA0003206699890000083
Figure GDA0003206699890000083

公式(18)获得的一阶矩信息实质上是原有的FCA算法提取出的相衬信息与暗场信息的混合形式,这种混合场信息,从物理上明确地融合了相衬信息和暗场信息,并且在某些情况下可以经过简单处理便突出其中的某一种信息。例如,在对轻材料均匀物体进行成像时,相位偏移量φc往往接近于0,因此这时候可以对公式(18)取对数,便可以获得类似公式(6)中的暗场信息,即暗场信息D′,即公式(20):The first-order moment information obtained by formula (18) is essentially a mixed form of phase contrast information and dark field information extracted by the original FCA algorithm. This mixed field information physically combines phase contrast information and dark field information. Field information, and in some cases can be simply processed to highlight one of them. For example, when imaging a uniform object with light materials, the phase offset φ c is often close to 0, so the logarithm of formula (18) can be taken at this time, and the dark field information similar to formula (6) can be obtained, That is, the dark field information D', that is, formula (20):

Figure GDA0003206699890000084
Figure GDA0003206699890000084

其中,第一项便是(6)中的暗场信息,第二项则是小量,第三项为常数量。Among them, the first item is the dark field information in (6), the second item is a small quantity, and the third item is a constant quantity.

另外,在对散射较少的物体进行成像时,对比度降低的程度不明显,则公式(18)则可以近似为公式(5)中的相衬信息,相应地,可以得到相衬信息P′。因此,公式(1)中的混合场对比度H可以据此进行信息定义。基于上述表达,无论是依据例如ASAXS算法对原始数据的直接提取方法,抑或是依据例如FCA算法对原始数据的间接提取方法,都可以通过本发明提供的混合场对比度H对相衬信息和暗场信息进行表达,因此,本发明所提出的混合场对比度H具有一定的普适性。In addition, when imaging an object with less scattering, the degree of contrast reduction is not obvious, then formula (18) can be approximated as the phase contrast information in formula (5), correspondingly, the phase contrast information P' can be obtained. Therefore, the mixed field contrast H in the formula (1) can be defined according to the information. Based on the above expression, whether it is based on the direct extraction method of the original data such as the AAXS algorithm, or the indirect extraction method of the original data according to, for example, the FCA algorithm, the mixed field contrast H provided by the present invention can be used to compare the phase contrast information and the dark field. Therefore, the hybrid field contrast H proposed by the present invention has certain universality.

在本发明实施例中,光栅相衬成像系统中,实验样本为活体小鼠,源光栅到相位光栅间距为1020mm,源光栅到分析光栅间距为1280mm,源光栅、相位光栅和分析光栅周期分别为16.8μm、4.2μm和2.4μm,X射线球管参数为35kVp和35mA,探测器像素为74.8μm×74.8μm,相位步进数为8。如图2本发明实施例中分别采用FCA、SAXS、ASAXS三种算法对应提取的吸收、相衬、暗场三种信息的成像对比图,以及图3本发明实施例中对ASAXS算法进行混合表征与FCA、SAXS算法分别提取的暗场信息的成像对比图所示,从图2可以看出本发明提出的混合场图像(即ASAXS的相衬图像)大部分区域与原有的相衬图像十分一致,而从图3可以看出经过公式(15)的取绝对值对数的简单操作后,本发明提出的混合场图像与图1中的暗场信息十分相似,以上结果表明提出的混合场物理表征量能够较好地反映相衬和暗场两种信息的特征。In the embodiment of the present invention, in the grating phase contrast imaging system, the experimental sample is a living mouse, the distance from the source grating to the phase grating is 1020 mm, the distance from the source grating to the analysis grating is 1280 mm, and the periods of the source grating, the phase grating and the analysis grating are respectively 16.8μm, 4.2μm and 2.4μm, the X-ray tube parameters are 35kVp and 35mA, the detector pixel is 74.8μm×74.8μm, and the number of phase steps is 8. In the embodiment of the present invention, as shown in FIG. 2, the imaging comparison diagram of the absorption, phase contrast, and dark field information corresponding to the extraction of the three algorithms of FCA, SAXS, and AAXS, and FIG. 3, the hybrid characterization of the AAXS algorithm in the embodiment of the present invention As shown in the imaging comparison diagram of the dark field information extracted by the FCA and SAXS algorithms respectively, it can be seen from Figure 2 that the mixed field image proposed by the present invention (that is, the phase contrast image of AAXS) is very similar to the original phase contrast image in most areas. It can be seen from Figure 3 that after the simple operation of taking the logarithm of the absolute value of formula (15), the mixed field image proposed by the present invention is very similar to the dark field information in Figure 1. The above results show that the proposed mixed field Physical characterization quantities can better reflect the characteristics of phase contrast and dark field information.

以上为基于光栅的X射线相衬成像中的一种实施例,对于其它X射线相衬成像实施例而言,可以按照公式(1)计算出本发明提出的混合场物理表征量,即可以先用传统方法计算出相衬信息和暗场信息,然后按照公式(1)进行计算得到混合场信息。The above is an embodiment of grating-based X-ray phase contrast imaging. For other X-ray phase contrast imaging embodiments, the physical characterization quantity of the mixed field proposed by the present invention can be calculated according to formula (1), that is, the first The phase contrast information and dark field information are calculated by the traditional method, and then the mixed field information is obtained by calculating according to formula (1).

本发明的另一个方面提出了一种相衬成像信息表征系统,应用上述的方法,系统包括:光学成像装置,用于获取相衬成像的原始数据;信息提取装置,用于根据光学成像装置获得的原始数据对相衬信息和暗场信息进行提取;信息表征装置,用于根据信息提取装置提取的相衬信息和暗场信息通过公式:Another aspect of the present invention provides a phase contrast imaging information characterization system, applying the above method, the system includes: an optical imaging device for acquiring original data of phase contrast imaging; an information extraction device for obtaining according to the optical imaging device Extract the phase contrast information and dark field information from the original data of

H=D′·sin(P′)H=D′·sin(P′)

进行表征,其中H为相衬信息和暗场信息的混合场对比度,P′为对相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Characterization, where H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the contrast of the dark field information The ratio of image contrast reduction after indexing the corresponding data.

可选地,光学成像装置为分析成像装置、边缘照射成像装置、光斑成像装置或光栅成像装置。光学成像装置作为原始数据的获取装置,其不限于基于分析晶体的分析成像装置(Analyzer-based Imaging,ABI)、基于边缘照射的边缘照射成像装置(Edge-Illumination,EI)、基于斑点的光斑成像装置(Speckle-based Imaging)以及基于光栅的光栅成像装置(Grating-based Imaging,GI)等各类相衬成像系统,以可以获取吸收信息、相衬信息、暗场信息三种信息的原始数据。Optionally, the optical imaging device is an analytical imaging device, an edge illumination imaging device, a spot imaging device or a grating imaging device. The optical imaging device is used as the acquisition device of raw data, which is not limited to the analysis imaging device based on analysis crystal (Analyzer-based Imaging, ABI), the edge illumination imaging device based on edge illumination (Edge-Illumination, EI), the spot imaging based on speckle Various phase contrast imaging systems such as Specckle-based Imaging and Grating-based Imaging (GI) can be used to obtain the original data of absorption information, phase contrast information and dark field information.

可选地,光栅成像装置为Talbot-Lau类型、几何投影类型或双相位光栅类型。光栅成像装置类别不限,即可以包括Talbot-Lau类型、几何投影类型和双相位光栅类型等各种类型系统的光栅成像装置。以可以获取吸收信息、相衬信息、暗场信息三种信息的原始数据。Optionally, the grating imaging device is a Talbot-Lau type, a geometric projection type or a dual phase grating type. The types of grating imaging devices are not limited, that is, grating imaging devices of various types of systems such as Talbot-Lau type, geometric projection type, and dual-phase grating type can be included. In order to obtain the original data of three kinds of information: absorption information, phase contrast information and dark field information.

可选地,光栅成像装置的照射待测位置位于光源与第一光栅之间或位于第一光栅和第二光栅之间。此外物体位置既可以是放置在光源与第一块光栅之间,也可以是放置在两块光栅之间。光栅成像装置一般为光源、第一光栅以及探测器,或光源、第一光栅、第二光栅以及探测器,因此,待测物体的照射待测位置可以位于光源和第一光栅之间,也可以位于第一光栅和第二光栅之间。照射待测位置的变化对于成像信息的获取以及成像信息的准确度有关。Optionally, the position to be measured for illumination of the grating imaging device is located between the light source and the first grating or between the first grating and the second grating. In addition, the position of the object can be placed between the light source and the first grating, or between the two gratings. The grating imaging device is generally a light source, a first grating, and a detector, or a light source, a first grating, a second grating, and a detector. Therefore, the position of the object to be measured can be located between the light source and the first grating, or it can be between the first grating and the second grating. The change of the position to be measured is related to the acquisition of imaging information and the accuracy of the imaging information.

可选地,光学成像装置的光源为常规X射线光源、微焦点X射线光源或同步辐射X射线光源,光学成像装置的探测器为能量积分型探测器或光子计数型探测器。此外X射线源既可以是常规X光源,也可以是微焦点X光源或者同步辐射X光源。此外探测器既可以是能量积分型探测器也可以是光子计数型探测器。Optionally, the light source of the optical imaging device is a conventional X-ray light source, a microfocus X-ray light source or a synchrotron radiation X-ray light source, and the detector of the optical imaging device is an energy integration detector or a photon counting detector. In addition, the X-ray source can be either a conventional X-ray source, a microfocus X-ray source or a synchrotron radiation X-ray source. In addition, the detector can be either an energy-integrating detector or a photon-counting detector.

本发明提出的一种相衬成像信息表征方法及系统,通过对提取的相衬信息和暗场信息进行处理,能够将原有的相衬信息和暗场信息表征在一张混合场图像中,该混合场图像与原有的相衬图像和暗场图像相比,图像表征清楚,相比于普通的图像融合方法具有更加明确的物理含义。The method and system for characterizing phase contrast imaging information proposed by the present invention can characterize the original phase contrast information and dark field information in a mixed field image by processing the extracted phase contrast information and dark field information. Compared with the original phase-contrast image and dark-field image, the mixed-field image has a clear image representation, and has a clearer physical meaning than the common image fusion method.

本发明的另一个方面提出了一种相衬成像信息表征系统,系统包括:处理器;存储器,其上存储有可在处理器上运行的计算机程序,其中所述计算机程序被所述处理器执行时实现以上实施例所述的方法步骤。本发明实施例中的光栅成像系统信息提取装置包括处理器(例如,微处理器(μP)、数字信号处理器(DSP)等)。处理器可以是用于执行本发明实施例描述的方法流程的不同动作的单一处理单元或者是多个处理单元。此外,装置可以包括具有非易失性或易失性存储器形式的至少一个计算机可读存储介质,以下将具体介绍。Another aspect of the present invention provides a system for characterizing phase contrast imaging information, the system comprising: a processor; a memory on which a computer program executable on the processor is stored, wherein the computer program is executed by the processor At the same time, the method steps described in the above embodiments are implemented. The apparatus for extracting information of a raster imaging system in the embodiment of the present invention includes a processor (eg, a microprocessor (μP), a digital signal processor (DSP), etc.). The processor may be a single processing unit or multiple processing units for executing different actions of the method flow described in the embodiments of the present invention. In addition, the apparatus may include at least one computer-readable storage medium in the form of non-volatile or volatile memory, as described in detail below.

本发明的另一个方面提出了一种计算机可读存储介质,所述计算机可读存储介质上存储有数据处理程序,所述数据处理程序被处理器执行时实现以上实施例所述的方法。本发明实施例的信息提取方法可以采取存储有指令的计算机可读存储介质上的计算机程序产品的形式,该计算机程序产品可供指令执行系统(例如,一个或多个处理器)使用或者结合指令执行系统使用。在本发明实施例的上下文中,计算机可读存储介质可以是能够包含、存储、传送、传播或传输指令的任意介质。例如,计算机可读存储介质可以包括但不限于电、磁、光、电磁、红外或半导体系统、装置、器件或传播介质。计算机可读存储介质的具体示例包括:磁存储装置,如磁带或硬盘(HDD);光存储装置,如光盘(CD-ROM);存储器,如随机存取存储器(RAM)或闪存;和/或有线/无线通信链路。Another aspect of the present invention provides a computer-readable storage medium, where a data processing program is stored thereon, and when the data processing program is executed by a processor, the methods described in the above embodiments are implemented. The information extraction methods of embodiments of the present invention may take the form of a computer program product on a computer-readable storage medium storing instructions, the computer program product being usable by an instruction execution system (eg, one or more processors) or in combination with instructions Execute system usage. In the context of embodiments of the present invention, a computer-readable storage medium can be any medium that can contain, store, communicate, propagate, or transmit instructions. For example, a computer-readable storage medium may include, but is not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, device, or propagation medium. Specific examples of computer-readable storage media include: magnetic storage devices, such as magnetic tapes or hard disks (HDDs); optical storage devices, such as compact disks (CD-ROMs); memories, such as random access memory (RAM) or flash memory; and/or wired/wireless communication link.

以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further describe the purpose, technical solutions and beneficial effects of the present invention in further detail. It should be understood that the above descriptions are only specific embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims (10)

1.一种相衬成像信息表征方法,其特征在于,所述方法包括:1. A method for characterizing phase contrast imaging information, wherein the method comprises: 获取相衬成像的原始数据;Obtain the raw data of phase contrast imaging; 根据所述原始数据对相衬信息和暗场信息进行提取;extracting phase contrast information and dark field information according to the original data; 对提取的所述相衬信息和暗场信息通过以下公式进行表征:The extracted phase contrast information and dark field information are characterized by the following formula: H=D′·sin(P′)H=D′·sin(P′) 其中,H为相衬信息和暗场信息的混合场对比度,P′为对所述相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对所述暗场信息的对应数据进行指数化操作后的图像对比度降低比值。Among them, H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the phase shift of the dark field The ratio of the image contrast reduction after the indexing operation of the corresponding data of the information. 2.根据权利要求1所述的方法,其特征在于,根据所述原始数据对相衬信息和暗场信息进行提取,是通过包括FCA算法、SAXS算法或ASAXS算法实现的。2 . The method according to claim 1 , wherein the extraction of phase contrast information and dark field information according to the original data is realized by including FCA algorithm, SAXS algorithm or AAXS algorithm. 3 . 3.根据权利要求1所述的方法,其特征在于,所述根据原始数据对相衬信息和暗场信息进行提取之后,还包括:3. The method according to claim 1, wherein after the phase contrast information and the dark field information are extracted according to the original data, the method further comprises: 对所述相衬信息的对应数据进行归一化操作获取调整到[-π,π]内的相位偏移P′;Perform a normalization operation on the corresponding data of the phase contrast information to obtain a phase offset P' adjusted to [-π, π]; 对所述暗场信息的对应数据进行指数化操作获取的图像对比度降低比值D′。The image contrast reduction ratio D' obtained by performing the indexing operation on the corresponding data of the dark field information. 4.一种相衬成像信息表征系统,应用根据权利要求1-3任一项所述的方法,其特征在于,所述系统包括:4. A phase contrast imaging information characterization system, applying the method according to any one of claims 1-3, wherein the system comprises: 光学成像装置,用于获取相衬成像的原始数据;An optical imaging device for acquiring raw data of phase contrast imaging; 信息提取装置,用于根据所述光学成像装置获得的所述原始数据对相衬信息和暗场信息进行提取;an information extraction device, configured to extract phase contrast information and dark field information according to the original data obtained by the optical imaging device; 信息表征装置,用于根据所述信息提取装置提取的所述相衬信息和暗场信息通过公式:An information characterization device, configured to pass the formula according to the phase contrast information and the dark field information extracted by the information extraction device: H=D′·sin(P′)H=D′·sin(P′) 进行表征,其中H为相衬信息和暗场信息的混合场对比度,P′为对所述相衬信息的对应数据调整到[-π,π]内的相位偏移,D′为对所述暗场信息的对应数据进行指数化操作后的图像对比度降低比值。to characterize, where H is the mixed field contrast of the phase contrast information and dark field information, P' is the phase offset adjusted to [-π, π] for the corresponding data of the phase contrast information, D' is the phase shift for the phase contrast information The ratio of image contrast reduction after the indexing operation of the corresponding data of the dark field information. 5.根据权利要求4所述的系统,其特征在于,所述光学成像装置为分析成像装置、边缘照射成像装置、光斑成像装置或光栅成像装置。5. The system according to claim 4, wherein the optical imaging device is an analytical imaging device, an edge illumination imaging device, a spot imaging device or a grating imaging device. 6.根据权利要求5所述的系统,其特征在于,所述光栅成像装置为Talbot-Lau类型、几何投影类型或双相位光栅类型。6. The system according to claim 5, wherein the grating imaging device is a Talbot-Lau type, a geometric projection type or a dual-phase grating type. 7.根据权利要求5所述的系统,其特征在于,所述光栅成像装置的照射待测位置位于光源与第一光栅之间或位于第一光栅和第二光栅之间。7 . The system according to claim 5 , wherein the position to be measured for illumination of the grating imaging device is located between the light source and the first grating or between the first grating and the second grating. 8 . 8.根据权利要求4所述的系统,其特征在于,所述光学成像装置的光源为常规X射线光源、微焦点X射线光源或同步辐射X射线光源,所述光学成像装置的探测器为能量积分型探测器或光子计数型探测器。8. The system according to claim 4, wherein the light source of the optical imaging device is a conventional X-ray light source, a microfocus X-ray light source or a synchrotron radiation X-ray light source, and the detector of the optical imaging device is an energy source Integrating detectors or photon counting detectors. 9.一种相衬成像信息表征系统,其特征在于,所述系统包括:9. A phase contrast imaging information characterization system, wherein the system comprises: 处理器;processor; 存储器,其上存储有可在所述处理器上运行的计算机程序,其中所述计算机程序被所述处理器执行时,实现权利要求1-3任一项所述的方法。A memory having stored thereon a computer program executable on the processor, wherein the computer program, when executed by the processor, implements the method of any one of claims 1-3. 10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有数据处理程序,所述数据处理程序被处理搭执行时,实现权利要求1-3任一项所述的方法。10. A computer-readable storage medium, characterized in that, a data processing program is stored on the computer-readable storage medium, and when the data processing program is processed and executed, the data processing program of any one of claims 1-3 is realized. Methods.
CN201911147099.7A 2019-11-21 2019-11-21 Information Characterization Method and System for Phase Contrast Imaging System Active CN111643100B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911147099.7A CN111643100B (en) 2019-11-21 2019-11-21 Information Characterization Method and System for Phase Contrast Imaging System

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911147099.7A CN111643100B (en) 2019-11-21 2019-11-21 Information Characterization Method and System for Phase Contrast Imaging System

Publications (2)

Publication Number Publication Date
CN111643100A CN111643100A (en) 2020-09-11
CN111643100B true CN111643100B (en) 2021-10-15

Family

ID=72342532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911147099.7A Active CN111643100B (en) 2019-11-21 2019-11-21 Information Characterization Method and System for Phase Contrast Imaging System

Country Status (1)

Country Link
CN (1) CN111643100B (en)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2509504B1 (en) * 2009-12-10 2018-07-18 Koninklijke Philips N.V. Scanning system for differential phase contrast imaging
US9117296B2 (en) * 2011-07-28 2015-08-25 Paul Scherrer Institut Method for image fusion based on principal component analysis
US9357975B2 (en) * 2013-12-30 2016-06-07 Carestream Health, Inc. Large FOV phase contrast imaging based on detuned configuration including acquisition and reconstruction techniques
US10096098B2 (en) * 2013-12-30 2018-10-09 Carestream Health, Inc. Phase retrieval from differential phase contrast imaging
DE102013219553A1 (en) * 2013-09-27 2015-04-02 Friedrich-Alexander-Universität Erlangen-Nürnberg Examination method of an object with an X-ray system for phase contrast imaging with stochastic phase scanning
EP3349657B1 (en) * 2015-09-16 2019-06-05 Koninklijke Philips N.V. Dark-field enhanced virtual x-ray colonoscopy
EP3452982B1 (en) * 2016-05-04 2019-10-02 Koninklijke Philips N.V. Feature suppression in dark field or phase contrast x-ray imaging
EP3378397A1 (en) * 2017-03-24 2018-09-26 Koninklijke Philips N.V. Sensitivity optimized patient positioning system for dark-field x-ray imaging
CN109377533B (en) * 2018-09-21 2023-01-24 上海交通大学 X-ray Grating Phase Contrast Imaging Reconstruction Method and System
CN110197486B (en) * 2019-07-02 2021-04-06 清华大学 X-ray grating integrated imaging information extraction method, system and storage medium

Also Published As

Publication number Publication date
CN111643100A (en) 2020-09-11

Similar Documents

Publication Publication Date Title
Croton et al. In situ phase contrast X-ray brain CT
Enders et al. A computational framework for ptychographic reconstructions
Chapman et al. High-resolution ab initio three-dimensional x-ray diffraction microscopy
Almasian et al. Validation of quantitative attenuation and backscattering coefficient measurements by optical coherence tomography in the concentration-dependent and multiple scattering regime
Yi et al. Imaging a full set of optical scattering properties of biological tissue by inverse spectroscopic optical coherence tomography
Almasian et al. OCT amplitude and speckle statistics of discrete random media
Chen et al. Volumetric quantification of lung nodules in CT with iterative reconstruction (ASiR and MBIR)
Kandel et al. Label-free tissue scanner for colorectal cancer screening
Zanette et al. Interlaced phase stepping in phase-contrast x-ray tomography
Zdora et al. Simulations of x-ray speckle-based dark-field and phase-contrast imaging with a polychromatic beam
Viermetz et al. High resolution laboratory grating-based X-ray phase-contrast CT
Kuo et al. Airway tapering: an objective image biomarker for bronchiectasis
CN110310347B (en) A filter iterative reconstruction method and system for X-ray grating phase contrast imaging
Garcia-Vicente et al. Textural features and SUV-based variables assessed by dual time point 18F-FDG PET/CT in locally advanced breast cancer
CN104434150B (en) Two-dimensional modulation transfer function measuring method of digital X-ray imaging system
CN116563413B (en) Correction parameter determining method, image reconstructing method, device, equipment and medium
Piccolomini et al. A model-based optimization framework for iterative digital breast tomosynthesis image reconstruction
Patel et al. Energy-specific optimization of attenuation thresholds for low-energy virtual monoenergetic images in renal lesion evaluation
Ongie et al. Optimizing model observer performance in learning-based CT reconstruction
CN111643100B (en) Information Characterization Method and System for Phase Contrast Imaging System
Leatham et al. X-ray phase and dark-field computed tomography without optical elements
Golbus et al. Ultra-low dose chest CT with silver filter and deep learning reconstruction significantly reduces radiation dose and retains quantitative information in the investigation and monitoring of lymphangioleiomyomatosis (LAM)
Ji et al. Quantitative evaluation on 3D fetus morphology via X‐ray grating based imaging technique
US20130236923A1 (en) Characteristic Parameters of Cells and Tissue from Quantitative Phase Imaging
JP7080214B6 (en) Identification of adipose tissue type

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant