CN111184523B - Three-dimensional image reconstruction method and system based on DR equipment - Google Patents
Three-dimensional image reconstruction method and system based on DR equipment Download PDFInfo
- Publication number
- CN111184523B CN111184523B CN202010052376.2A CN202010052376A CN111184523B CN 111184523 B CN111184523 B CN 111184523B CN 202010052376 A CN202010052376 A CN 202010052376A CN 111184523 B CN111184523 B CN 111184523B
- Authority
- CN
- China
- Prior art keywords
- scanning data
- correction
- ray
- pixel
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 230000006870 function Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 7
- 210000000056 organ Anatomy 0.000 abstract description 5
- 238000002591 computed tomography Methods 0.000 abstract description 3
- 230000002159 abnormal effect Effects 0.000 description 8
- 230000003068 static effect Effects 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000002496 gastric effect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 238000001739 density measurement Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
-
- 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]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- 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
- G06T2207/10124—Digitally reconstructed radiograph [DRR]
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Computer Graphics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a three-dimensional image reconstruction method and a three-dimensional image reconstruction system based on DR equipment, wherein the method comprises the following steps: acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current; carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data; and carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data in sequence to obtain a three-dimensional reconstruction image. The invention realizes the extension of the DR equipment image from two dimensions to three dimensions, the obtained three-dimensional image has no tissue and organ aliasing, and the condition of the body to be detected can be reflected more clearly and accurately. A two-dimensional DR image and a three-dimensional CT image can be obtained through one-time DR scanning, and CT scanning is not needed after DR detection, so that the X-ray irradiation dose of a patient is greatly reduced, and the detection cost can be reduced.
Description
Technical Field
The invention relates to the technical field of image processing, in particular to a three-dimensional image reconstruction method and a three-dimensional image reconstruction system based on DR equipment.
Background
DR (Digital Radiography) refers to a technology of directly performing Digital Radiography under the control of a computer, that is, an X-ray detector converts X-ray information penetrating a human body into a Digital signal, and a computer performs post-processing and display on an image. The DR system mainly comprises an X-ray generating device, an X-ray detector, a system controller, an image monitor, an image processing workstation and the like. The DR adopts a digital technology, so that various image post-processing can be performed according to clinical requirements, and abundant functions such as automatic image processing, edge enhancement, magnification roaming, image splicing, window width and level adjustment of interest area windows, distance, area and density measurement and the like can be realized. The DR technology has wide dynamic range, high X-ray photon detection efficiency (DQE), wide exposure latitude and good image quality even under slightly poor exposure conditions. In addition, DR has the characteristics of low dose, high spatial resolution, short scanning time, low cost and the like, is widely applied to the fields of physical examination and medical image diagnosis, and is one of main equipment for medical image diagnosis.
In general, DR devices fall into two broad categories: static DR and dynamic DR. Static DR, as the name implies, can only take a static image. And the doctor guides the patient to complete the positioning and equipment setting and then adopts a single exposure form to acquire a static two-dimensional image of the patient. The shooting mode is similar to the shooting mode in daily life, so that the static DR image can only feed back two-dimensional image information of the exposure time of the patient. Failure to take a picture may result if the setup is poor or if the patient moves during the exposure. The dynamic DR changes the traditional 'photographing' into 'shooting', and introduces a continuous exposure concept into the DR shooting process, so that not only can a DR image of a patient be shot, but also clear and accurate images can be successfully shot at one time in the shooting of a plurality of complex body positions by observing the condition of the patient, and the condition that the patient needs to repeatedly shoot for many times in the static DR can be effectively avoided. Because of the addition of the perspective function, doctors can observe the lesion part at a plurality of different angles, thereby further avoiding missed diagnosis, and the dynamic DR can replace the traditional gastrointestinal machine to carry out gastrointestinal radiography.
Dynamic DR converts traditional "photographic" static acquisition to "photographic" perspective acquisition, providing more information to the physician to assist in diagnosis. However, the images acquired by the static DR device and the dynamic DR device are both two-dimensional images, that is, aliasing phenomena exist in all organs on the X-ray penetration path, and although the DR device has a high spatial resolution, the actual positions and related information of various organs and tissues cannot be accurately displayed, an accurate three-dimensional spatial image cannot be provided, and clinical application of the DR device is limited.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: a three-dimensional image reconstruction method and system based on DR equipment are provided, which can obtain a three-dimensional image without aliasing of tissues and organs.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a three-dimensional image reconstruction method based on DR equipment comprises the following steps:
acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
and sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain a three-dimensional reconstruction image.
The invention adopts another technical scheme that:
a three-dimensional image reconstruction system based on DR equipment comprises a high-voltage generator, an X-ray emitter, a rotating body with an angle measuring and returning function, an X-ray detector, a rack and a three-dimensional image reconstruction terminal, wherein the high-voltage generator is electrically connected with the X-ray emitter, the rotating body with the angle measuring and returning function is used for placing an object to be detected, the rotating body with the angle measuring and returning function is arranged between the X-ray emitter and the X-ray detector, the X-ray emitter and the X-ray detector are respectively and fixedly arranged on the rack, and the three-dimensional image reconstruction terminal is electrically connected with the X-ray detector;
the three-dimensional image reconstruction terminal comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the following steps:
acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
and sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain a three-dimensional reconstruction image.
The invention has the beneficial effects that: the X-ray scanning data of the body to be detected at different angles are collected, then the scanning data are corrected and three-dimensionally reconstructed, the DR equipment image is expanded from two dimensions to three dimensions, the obtained three-dimensional image has no tissues and viscera aliasing, and the condition of the body to be detected can be more clearly and accurately reflected. A two-dimensional DR image and a three-dimensional CT image can be obtained through one-time DR scanning, and CT scanning is not needed after DR detection, so that the X-ray irradiation dose of a patient is greatly reduced, and the detection cost can be reduced.
Drawings
Fig. 1 is a flowchart of a three-dimensional image reconstruction method based on a DR device according to a first embodiment of the present invention;
fig. 2 is a schematic structural diagram of a three-dimensional image reconstruction system according to a second embodiment of the present invention;
fig. 3 is another schematic structural diagram of a three-dimensional image reconstruction system according to a second embodiment of the present invention;
fig. 4 is another schematic structural diagram of a three-dimensional image reconstruction system according to a second embodiment of the present invention;
fig. 5 is a schematic diagram of a three-dimensional image reconstruction terminal according to a second embodiment of the present invention.
Description of reference numerals:
100. a three-dimensional image reconstruction terminal; 1. a memory; 2. a processor.
200. A high voltage generator; 300. an X-ray emitter; 400. a rotating body with angle measuring and returning functions; 500. an X-ray detector; 600. a frame; 700. the body is to be detected.
Detailed Description
In order to explain technical contents, achieved objects, and effects of the present invention in detail, the following description is made with reference to the accompanying drawings in combination with the embodiments.
The most key concept of the invention is as follows: the DR equipment is adopted to collect X-ray scanning data of the object to be detected at different angles, and then the scanning data is corrected and three-dimensionally reconstructed, so that the DR equipment image is expanded from two dimensions to three dimensions.
Referring to fig. 1, a three-dimensional image reconstruction method based on a DR apparatus includes:
acquiring X-ray scanning data of objects to be detected at different angles on DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
and sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain a three-dimensional reconstruction image.
From the above description, the beneficial effects of the present invention are: the X-ray scanning data of the body to be detected at different angles are collected, then the scanning data are corrected and three-dimensionally reconstructed, the DR equipment image is expanded from two dimensions to three dimensions, the obtained three-dimensional image has no tissues and viscera aliasing, and the condition of the body to be detected can be more clearly and accurately reflected. A two-dimensional DR image and a three-dimensional image can be obtained through one-time DR scanning, and CT scanning is not needed after DR detection, so that the X-ray irradiation dose of a patient is greatly reduced, and the detection cost can be reduced.
Further, performing dark correction on the acquired X-ray scanning data specifically includes:
collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data;
and subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction.
As can be seen from the above description, the purpose of dark correction is to remove the problem of the non-uniformity of the data collected by the X-ray detector of the DR equipment.
Further, the step of correcting the dead pixel and the dead line specifically comprises the following steps:
presetting a pixel range;
judging whether a pixel point of the scanning data after dark correction is in the preset pixel range;
if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point.
As can be seen from the above description, the purpose of performing the dead pixel and dead line correction is to remove the dead pixel and the dead line existing in the X-ray detector. The general DR equipment X-ray detector has different numbers of dead pixels and dead lines, the surface of the problem pixel points in the image is 0 or other abnormal values, the abnormal values are certain values and do not change along with the scanning conditions, if the abnormal values are not corrected, the filtering is abnormal, and the problem of serious ring artifacts in the image or image does not exist.
Further, the air correction specifically includes:
collecting air scanning data under different voltage and current conditions to obtain an air correction table;
calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current;
performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel and dead pixel line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data;
and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data.
As can be seen from the above description, the purpose of performing the air correction is to correct the problems of the X-ray detector such as non-uniformity and ghost of the X-ray reception, by calculating the incident ray intensity of the X-ray and performing ln operation on the scan data to obtain the scan data that can be finally used for three-dimensional reconstruction.
Further, cone angle weighting on the corrected scan data specifically includes:
and performing cone angle cosine weighting processing on the corrected scanning data.
As can be seen from the above description, the X-ray detector of the DR apparatus has a large X-ray cone angle, which can cause the scanned data to be mismatched from the actual data to cause artifacts. Therefore, it is desirable to reduce the effect of cone angle on the scan data, with higher cone angles being weighted lower and lower cone angles being weighted higher.
Further, the filtering process specifically includes:
and performing filtering processing on the scanning data after cone angle weighting processing by adopting a filter.
As can be seen from the above description, since the back projection operation needs to model the data at different projection angles back into the three-dimensional reconstructed image, this process enhances the low frequency components in the reconstructed image, which causes the image to become blurred and uneven. To suppress this phenomenon, a filter is required to suppress the low frequency region of the image while enhancing the high frequency information. And filtering the scanning data by adopting a filter to enhance the boundary information of the scanning data to obtain a final reconstructed image.
Further, the sequentially performing cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain a three-dimensional reconstructed image further includes:
removing noise points in the three-dimensional reconstruction image and smoothing the boundary of the three-dimensional reconstruction image.
From the above description, it can be known that removing the noise point and the smooth boundary can obtain a three-dimensional reconstructed image with better quality, and the removing of the noise point and the smooth boundary can be performed by using the existing known means.
Referring to fig. 2, another technical solution related to the present invention is:
a three-dimensional image reconstruction system based on DR equipment comprises a high voltage generator 200, an X-ray emitter 300, a rotating body 400 with an angle measuring and returning function, an X-ray detector 500, a rack 600 and a three-dimensional image reconstruction terminal 100, wherein the high voltage generator 200 is electrically connected with the X-ray emitter 300, the rotating body 400 with the angle measuring and returning function is used for placing a body 700 to be detected, the rotating body 400 with the angle measuring and returning function is arranged between the X-ray emitter 300 and the X-ray detector 500, the X-ray emitter 300 and the X-ray detector 500 are respectively and fixedly arranged on the rack 600, and the three-dimensional image reconstruction terminal 100 is electrically connected with the X-ray detector 500;
the three-dimensional image reconstruction terminal 100 comprises a memory 1, a processor 2 and a computer program stored on the memory 1 and executable on the processor 2, the processor 2 implementing the following steps when executing the computer program:
acquiring X-ray scanning data of objects to be detected at different angles on DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
and sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain a three-dimensional reconstruction image.
Further, performing dark correction on the acquired X-ray scanning data specifically includes:
collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data;
and subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction.
Further, the step of correcting the dead pixel and the dead line specifically comprises the following steps:
presetting a pixel range;
judging whether a pixel point of the scanning data after dark correction is in the preset pixel range;
if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point.
Further, the air correction specifically includes:
collecting air scanning data under different voltage and current conditions to obtain an air correction meter;
calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current;
performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel and dead pixel line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data;
and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data.
Further, cone angle weighting on the corrected scan data specifically includes:
and performing cone angle cosine weighting processing on the corrected scanning data.
Further, the filtering process specifically includes:
and performing filtering processing on the scanning data after cone angle weighting processing by adopting a filter.
Further, the processor 2, when executing the computer program, further implements the following steps:
the method for sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain the three-dimensional reconstruction image further comprises the following steps:
removing noise points in the three-dimensional reconstruction image and smoothing the boundary of the three-dimensional reconstruction image.
Example one
Referring to fig. 1, a first embodiment of the present invention is:
a three-dimensional image reconstruction method based on DR equipment comprises the following steps:
s1, acquiring X-ray scanning data of a to-be-detected body on DR equipment at different angles under the conditions of preset voltage and preset current.
In this embodiment, X-ray scan data is acquired by a DR device. The body to be detected can be a person or an animal, the body to be detected can stand or be fixed on a rotatable rotating body with the angle measuring and returning functions, the rotating body with the angle measuring and returning functions drives the body to be detected to rotate, so that X-ray scanning data at different angles can be obtained, the rotating speed can be adjusted according to needs, and the body to be detected needs to be kept still in the rotating process. The return function means that the rotating body can record and transmit corresponding angle information. The preset voltage and the preset current can also be adjusted according to actual conditions. The X-ray emitter can be fixed through a rack, and the rack can be a fixed rack or an adjustable rack with a slide rail.
And S2, carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data.
In this embodiment, the dark correction on the acquired X-ray scanning data specifically includes: collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data; and subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction.
The purpose of dark correction is to remove the problem of data non-uniformity of X-ray detector of DR equipment. And performing dark correction, namely subtracting the pixel value of the corresponding pixel point in the corresponding dark correction image from the pixel value of each pixel point in the obtained scanning data.
In this embodiment, the specific step of correcting the dead pixel and the dead line is as follows: presetting a pixel range; judging whether a pixel point of the scanning data after dark correction is in the preset pixel range or not; if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point; if yes, the pixel value of the pixel point is unchanged.
The purpose of the dead pixel dead line correction is to remove the dead pixel and the dead line existing in the X-ray detector. The general DR equipment X-ray detector has different numbers of dead pixels and dead lines, the surface of the problem pixel points in the image is 0 or other abnormal values, the abnormal values are certain values and do not change along with the scanning conditions, if the abnormal values are not corrected, the filtering is abnormal, and no image or serious ring artifact problem occurs in the image.
In this embodiment, the air correction specifically includes: collecting air scanning data under different voltage and current conditions to obtain an air correction table; calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current; performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel and dead pixel line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data; and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data. In this embodiment, the average value of all the pixel points in the air correction table may be calculated as the incident light intensity, or the average value of any one or more pixel points may be selected as the incident light intensity.
The purpose of the air correction is to correct the problems of the X-ray detector such as the receiving unevenness and the ghost of the X-ray, and the scanning data can be finally used for three-dimensional reconstruction by calculating the incident ray intensity of the X-ray and performing ln operation on the scanning data. The air scan data, i.e., the data obtained by X-rays that penetrate only the air and are received by the detector.
And S3, cone angle weighting, filtering processing and back projection are sequentially carried out on the corrected scanning data, and a three-dimensional reconstruction image is obtained.
In this embodiment, the cone angle weighting on the corrected scan data specifically includes: and performing cone angle cosine weighting processing on the corrected scanning data. The X-ray detector of a DR apparatus has a large X-ray cone angle, which can cause scan data to mismatch with reality and cause artifacts. Therefore, it is desirable to reduce the effect of cone angle on the scan data, with higher cone angles being weighted lower and lower cone angles being weighted higher. In addition, the scanning data can be weighted by designing linear, gaussian and other weights.
In this embodiment, the filtering process specifically includes: and performing filtering processing on the scanning data after cone angle weighting processing by adopting a filter. Since the back projection operation needs to model the data at different projection angles back into the three-dimensional reconstructed image, this process relatively enhances the low-frequency components in the reconstructed image, which causes the image to become blurred and uneven. To suppress this phenomenon, a filter is required to suppress the low frequency region of the image while enhancing the high frequency information. And filtering the scanning data by adopting a filter to enhance the boundary information of the scanning data to obtain a final reconstructed image.
The back projection process is a process of calculating scanning data corresponding to any three-dimensional reconstruction pixel point under any projection angle and accumulating the data to the position of a reconstruction point to obtain final reconstruction data.
In this embodiment, after step S3, the method further includes:
removing noise points in the three-dimensional reconstruction image and smoothing the boundary of the three-dimensional reconstruction image. And finally, displaying the processed three-dimensional reconstruction image.
Example two
Referring to fig. 2 to 5, a second embodiment of the present invention is:
a three-dimensional image reconstruction system based on DR equipment, as shown in fig. 2 to 4, includes a high voltage generator 200, an X-ray emitter 300, a rotating body 400 with an angle measuring and returning function, an X-ray detector 500, a rack 600 and a three-dimensional image reconstruction terminal 100, wherein the high voltage generator 200 is electrically connected to the X-ray emitter 300, the rotating body 400 with an angle measuring and returning function is used for placing an object to be detected 700, the rotating body 400 with an angle measuring and returning function is disposed between the X-ray emitter 300 and the X-ray detector 500, the X-ray emitter 300 and the X-ray detector 500 are respectively and fixedly disposed on the rack 600, and the three-dimensional image reconstruction terminal 100 is electrically connected to the X-ray detector 500. In this embodiment, the return function means that the rotating body can record and transmit the corresponding angle information to the three-dimensional image reconstruction terminal 100.
In fig. 2, a gantry 600 is an arc structure, an X-ray emitter 300 and an X-ray detector 500 are respectively and fixedly disposed at two ends of the gantry 600, and a rotating body 400 with angle measuring and returning functions is a device with an angle controller, which can adjust the angle of an object 700 to be detected according to needs.
In fig. 3, the gantry 600 is a slide rail type, the gantry 600 includes two slide rails, the X-ray emitter 300 is fixed on one of the slide rails, and the X-ray detector 500 is fixed on the other slide rail, so that the positions of the X-ray emitter 300 and the X-ray detector 500 on the gantry 600 can be adjusted as required.
In fig. 4, the gantry 600 is a column, and includes two columns, the X-ray emitter 300 is fixed on one of the columns, and the X-ray detector 500 is fixed on the other column, so that the positions of the two columns can be adjusted as required.
As shown in fig. 5, the three-dimensional image reconstruction terminal 100 includes a memory 1, a processor 2 and a computer program stored on the memory 1 and executable on the processor 2, and the processor 2 implements the following steps when executing the computer program:
acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
and carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data in sequence to obtain a three-dimensional reconstruction image.
Further, performing dark correction on the acquired X-ray scanning data specifically includes:
collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data;
and subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction.
Further, the step of correcting the dead pixel dead line specifically comprises the following steps:
presetting a pixel range;
judging whether a pixel point of the scanning data after dark correction is in the preset pixel range or not;
if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point.
Further, the air correction specifically includes:
collecting air scanning data under different voltage and current conditions to obtain an air correction meter;
calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current;
performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel dead line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data;
and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data.
Further, cone angle weighting on the corrected scan data specifically includes:
and performing cone angle cosine weighting processing on the corrected scanning data.
Further, the filtering specifically includes:
and performing filtering processing on the scanning data after cone angle weighting processing by adopting a filter.
Further, the processor 2, when executing the computer program, further implements the following steps:
the method for sequentially carrying out cone angle weighting, filtering processing and back projection on the corrected scanning data to obtain the three-dimensional reconstruction image further comprises the following steps:
removing noise points in the three-dimensional reconstruction image and smoothing the boundary of the three-dimensional reconstruction image.
In conclusion, the three-dimensional image reconstruction method and the three-dimensional image reconstruction system based on the DR equipment provided by the invention expand the DR equipment image from two dimensions to three dimensions, and the obtained three-dimensional image has no tissue and organ aliasing, so that the condition of the object to be detected can be reflected more clearly and accurately.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all equivalent changes made by using the contents of the present specification and the drawings, or applied directly or indirectly to the related technical fields, are included in the scope of the present invention.
Claims (5)
1. A three-dimensional image reconstruction method based on DR equipment is characterized by comprising the following steps:
acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
cone angle weighting, filtering and back projection are sequentially carried out on the corrected scanning data to obtain a three-dimensional reconstruction image;
the dark correction of the acquired X-ray scanning data specifically includes:
collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data;
subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction;
the dead pixel and dead line correction specifically comprises the following steps:
presetting a pixel range;
judging whether a pixel point of the scanning data after dark correction is in the preset pixel range;
if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point;
the air correction is specifically carried out as follows:
collecting air scanning data under different voltage and current conditions to obtain an air correction meter;
calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current;
performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel and dead pixel line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data;
and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data.
2. The DR device-based three-dimensional image reconstruction method of claim 1, wherein the cone angle weighting of the corrected scan data is specifically:
and performing cone angle cosine weighting processing on the corrected scanning data.
3. The DR device-based three-dimensional image reconstruction method of claim 2, wherein the filtering process is specifically:
and performing filtering processing on the scanning data after cone angle weighting processing by adopting a filter.
4. The DR equipment-based three-dimensional image reconstruction method of claim 1, wherein the sequentially performing cone angle weighting, filtering and back-projection on the corrected scan data to obtain a three-dimensional reconstructed image further comprises:
removing noise points in the three-dimensional reconstruction image and smoothing the boundary of the three-dimensional reconstruction image.
5. A three-dimensional image reconstruction system based on DR equipment is characterized by comprising a high-voltage generator, an X-ray emitter, a rotating body with an angle measuring and returning function, an X-ray detector, a rack and a three-dimensional image reconstruction terminal, wherein the high-voltage generator is electrically connected with the X-ray emitter, the rotating body with the angle measuring and returning function is used for placing a to-be-detected body, the rotating body with the angle measuring and returning function is arranged between the X-ray emitter and the X-ray detector, the X-ray emitter and the X-ray detector are respectively and fixedly arranged on the rack, and the three-dimensional image reconstruction terminal is electrically connected with the X-ray detector;
the three-dimensional image reconstruction terminal comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the following steps:
acquiring X-ray scanning data of the to-be-detected body at different angles on the DR equipment under the conditions of preset voltage and preset current;
carrying out dark correction, dead pixel dead line correction and air correction on the acquired X-ray scanning data in sequence to obtain corrected scanning data;
cone angle weighting, filtering and back projection are sequentially carried out on the corrected scanning data to obtain a three-dimensional reconstruction image;
the dark correction of the acquired X-ray scanning data specifically includes:
collecting data of an X-ray detector under the condition of not emitting X-rays to obtain dark correction map data;
subtracting the dark correction image data from the acquired X-ray scanning data to obtain scanning data after dark correction;
the dead pixel and dead line correction specifically comprises the following steps:
presetting a pixel range;
judging whether a pixel point of the scanning data after dark correction is in the preset pixel range;
if not, taking the pixel value of the pixel point with the adjacent pixel value in the preset pixel range as the pixel value of the pixel point;
the air correction is specifically carried out as follows:
collecting air scanning data under different voltage and current conditions to obtain an air correction table;
calculating according to the air correction table to obtain the incident light intensity under the conditions of the preset voltage and the preset current;
performing ln logarithmic operation on the incident light intensity, the scanning data after dead pixel and dead pixel line correction and the air scanning data respectively to obtain logarithmic incident light intensity, logarithmic scanning data and logarithmic air scanning data;
and subtracting the difference value of the logarithmic scanning data and the logarithmic air scanning data from the logarithmic incident light intensity to obtain corrected scanning data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010052376.2A CN111184523B (en) | 2020-01-17 | 2020-01-17 | Three-dimensional image reconstruction method and system based on DR equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010052376.2A CN111184523B (en) | 2020-01-17 | 2020-01-17 | Three-dimensional image reconstruction method and system based on DR equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111184523A CN111184523A (en) | 2020-05-22 |
| CN111184523B true CN111184523B (en) | 2023-03-10 |
Family
ID=70685030
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010052376.2A Active CN111184523B (en) | 2020-01-17 | 2020-01-17 | Three-dimensional image reconstruction method and system based on DR equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111184523B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111419255A (en) * | 2020-03-24 | 2020-07-17 | 深圳锐探科技有限公司 | CT imaging method and device |
| CN113223109B (en) * | 2021-05-27 | 2022-08-12 | 明峰医疗系统股份有限公司 | Dead pixel correction method based on CT axis scanning |
| CN113570705B (en) * | 2021-07-28 | 2024-04-30 | 广州瑞多思医疗科技有限公司 | Three-dimensional dose reconstruction method, device, computer equipment and storage medium |
| CN113570708A (en) * | 2021-07-30 | 2021-10-29 | 重庆市特种设备检测研究院 | Defect three-dimensional modeling method and device and computer readable storage medium |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005077627A (en) * | 2003-08-29 | 2005-03-24 | Fuji Photo Film Co Ltd | Dark correction method and device for radiation image information reader |
| JP2005296340A (en) * | 2004-04-12 | 2005-10-27 | Canon Inc | Cone-beam x-ray ct apparatus and method of obtaining image using the same |
| CN1881075A (en) * | 2005-06-15 | 2006-12-20 | 上海医疗器械厂有限公司 | X-ray digital imaging correction method |
| CN101371786A (en) * | 2007-08-24 | 2009-02-25 | 北京师范大学珠海分校 | Method and system of X ray image three-dimensional reconstruction |
| CN101795349A (en) * | 2010-03-24 | 2010-08-04 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing faulty wire of x-ray flat panel detector image and device thereof |
| CN102762151A (en) * | 2010-01-21 | 2012-10-31 | 卡尔斯特里姆保健公司 | Four-dimensional volume imaging system |
| CN103366389A (en) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | CT (computed tomography) image reconstruction method |
| CN103767723A (en) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | Cone beam CT three-dimensional digital imaging method based on C-shaped arm |
| WO2019000601A1 (en) * | 2017-06-27 | 2019-01-03 | 西安立人医学科技有限公司 | Cone beam ct multi-directional scanning apparatus |
| CN109730713A (en) * | 2019-01-28 | 2019-05-10 | 深圳市纳诺艾医疗科技有限公司 | A kind of second-order fluorescence wide-angle Multifunctional imaging equipment |
| WO2019149809A1 (en) * | 2018-02-02 | 2019-08-08 | Gom Gmbh | Method for calibrating a measurement device, and computer program for same |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4172201B2 (en) * | 2002-04-30 | 2008-10-29 | コニカミノルタホールディングス株式会社 | Radiation imaging apparatus and radiation image forming apparatus |
| KR102104534B1 (en) * | 2013-06-12 | 2020-04-27 | 삼성전자주식회사 | X-ray imaging apparatus and x-ray imaging apparatus control method |
-
2020
- 2020-01-17 CN CN202010052376.2A patent/CN111184523B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005077627A (en) * | 2003-08-29 | 2005-03-24 | Fuji Photo Film Co Ltd | Dark correction method and device for radiation image information reader |
| JP2005296340A (en) * | 2004-04-12 | 2005-10-27 | Canon Inc | Cone-beam x-ray ct apparatus and method of obtaining image using the same |
| CN1881075A (en) * | 2005-06-15 | 2006-12-20 | 上海医疗器械厂有限公司 | X-ray digital imaging correction method |
| CN101371786A (en) * | 2007-08-24 | 2009-02-25 | 北京师范大学珠海分校 | Method and system of X ray image three-dimensional reconstruction |
| CN102762151A (en) * | 2010-01-21 | 2012-10-31 | 卡尔斯特里姆保健公司 | Four-dimensional volume imaging system |
| CN101795349A (en) * | 2010-03-24 | 2010-08-04 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing faulty wire of x-ray flat panel detector image and device thereof |
| CN103767723A (en) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | Cone beam CT three-dimensional digital imaging method based on C-shaped arm |
| CN103366389A (en) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | CT (computed tomography) image reconstruction method |
| WO2019000601A1 (en) * | 2017-06-27 | 2019-01-03 | 西安立人医学科技有限公司 | Cone beam ct multi-directional scanning apparatus |
| WO2019149809A1 (en) * | 2018-02-02 | 2019-08-08 | Gom Gmbh | Method for calibrating a measurement device, and computer program for same |
| CN109730713A (en) * | 2019-01-28 | 2019-05-10 | 深圳市纳诺艾医疗科技有限公司 | A kind of second-order fluorescence wide-angle Multifunctional imaging equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111184523A (en) | 2020-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111184523B (en) | Three-dimensional image reconstruction method and system based on DR equipment | |
| JP3257630B2 (en) | Display device for identifying coronary arteries | |
| US9044190B2 (en) | C-arm computerized tomography system | |
| US20090060121A1 (en) | Computed tomography data acquisition apparatus and method | |
| CN105030266A (en) | X-ray computer tomographic apparatus and scan plan setting supporting apparatus | |
| US8249216B2 (en) | X-ray moving image radiographing apparatus | |
| JP4537129B2 (en) | System for scanning objects in tomosynthesis applications | |
| JP6640527B2 (en) | X-ray CT system | |
| JP2007089674A (en) | Shape of appearance measuring apparatus and x-ray ct apparatus | |
| US9194965B2 (en) | System and method for X-ray image acquisition and processing | |
| CN103860185B (en) | Determining of multipotency spirogram picture | |
| JPWO2005011502A1 (en) | Radiation tomography equipment | |
| JP3897925B2 (en) | Cone beam CT system | |
| EP3349655B1 (en) | Tomography apparatus and controlling method for the same | |
| CN101721222B (en) | Method for correcting effect of bed board and positioning auxiliary device on image quality | |
| KR101768520B1 (en) | A method of integrated operation of chest X-ray digital radiography and chest digital tomosynthesis | |
| US20220071578A1 (en) | Improved method of acquiring a radiographic scan of a region-of-interest in a metal containing object | |
| US20230145920A1 (en) | Systems and methods for motion detection in medical images | |
| US20220414832A1 (en) | X-ray imaging restoration using deep learning algorithms | |
| CN110084753B (en) | Dynamic DR image splicing method and terminal based on multi-sensor fusion | |
| CN109243585B (en) | Medical image generation method, medical image processing system and interaction method thereof | |
| CN110730977B (en) | Low dose imaging method and device | |
| US11622741B2 (en) | System and method for imaging a subject | |
| JP2014195492A (en) | X-ray ct scanner | |
| KR20190002960A (en) | Method for task based CT reconstruction protocol optimization using body-recognition |
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 |