CN113345043B - Method, device, medium and electronic equipment for eliminating CT image metal artifact - Google Patents
Method, device, medium and electronic equipment for eliminating CT image metal artifact Download PDFInfo
- Publication number
- CN113345043B CN113345043B CN202110610275.7A CN202110610275A CN113345043B CN 113345043 B CN113345043 B CN 113345043B CN 202110610275 A CN202110610275 A CN 202110610275A CN 113345043 B CN113345043 B CN 113345043B
- Authority
- CN
- China
- Prior art keywords
- curve
- sampling
- radar
- values
- bulb
- 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
- 239000002184 metal Substances 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005070 sampling Methods 0.000 claims abstract description 184
- 238000010606 normalization Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 230000003287 optical effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- 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/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- 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/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The disclosure provides a method, a device, a medium and electronic equipment for correcting a bulb tube sampling value. The method comprises the following steps: determining a position of a metal in the subject based on the pre-scan information; performing a tomographic scan of a subject site comprising a metal, the tomographic scan including a bulb X-ray scan and at least one lidar scan; obtaining a corresponding number of bulb tube sampling values and radar sampling values at the same position of the detected part; generating a first curve according to the radar sampling value, and generating a second curve according to the bulb sampling value; correcting the second curve based on the first curve, and eliminating metal artifacts in the second curve according to the characteristics of the radar sensor, the method and the device enable the radar sensor and the bulb sensor to be sampled simultaneously, enable sampling positions of the radar sensor and the bulb sensor to be completely registered, correct the bulb sampling value comprising metal artifact information by means of first curve information generated by radar sampling values, and enable medical images generated based on the corrected bulb sampling values to eliminate metal artifacts.
Description
Technical Field
The disclosure relates to the field of medical imaging, in particular to a method, a device, a medium and electronic equipment for eliminating metal artifacts of CT images.
Background
In the case of an electronic computed tomography (CT scan, english full Computed Tomography) of a patient, if metal (e.g., a prosthesis or a stent) is present in the patient, metal artifacts may be present in the CT image, thereby affecting the doctor's diagnosis.
At present, in a correction algorithm for metal artifacts in a CT image, generally, a metal region is segmented according to a CT value obtained in the CT image, data correction is performed by back-projection, and then the data after correction is used to replace the data in the metal region.
But this method of eliminating metal artifacts results in image loss resolution. Therefore, the disclosure proposes a method for correcting a sampling value of a bulb tube, so as to solve one of the above technical problems.
Disclosure of Invention
The disclosure aims to provide a method, a device, a medium and an electronic device for eliminating metal artifacts of a CT image, which can solve at least one technical problem. The specific scheme is as follows:
according to a first aspect of the present disclosure, there is provided a method for eliminating metal artifacts of CT images, comprising:
determining a position of a metal in the subject based on the pre-scan information;
Performing a tomographic scan of a subject site containing the metal, the tomographic scan including a bulb X-ray scan and at least one lidar scan;
obtaining a ball tube sampling value and a radar sampling value which are corresponding to the same position of the detected part, wherein the ball tube sampling position corresponding to the ball tube sampling value has a one-to-one correspondence with the radar sampling position corresponding to the radar sampling value;
generating a first curve according to the radar sampling value, and generating a second curve according to the bulb sampling value;
And correcting the second curve based on the first curve, and eliminating metal artifacts in the second curve.
Optionally, the generating a first curve according to the radar sampling values includes:
Carrying out normalization processing on the radar sampling value of each sampling position of the radar to obtain a normalization value of each sampling position of the radar;
And carrying out fitting processing on the normalized values based on each sampling position of the radar, and generating a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
Optionally, the generating a second curve according to the bulb sampling value includes:
Normalizing the bulb tube sampling value of each sampling position of the bulb tube to obtain a normalized value of each sampling position of the bulb tube;
and fitting the normalized values based on each sampling position of the bulb tube to generate a second curve, wherein the second curve comprises a recess formed by metal artifacts.
Optionally, the correcting the second curve based on the first curve, and eliminating metal artifacts in the second curve, includes:
and correcting the second curve equation by the first curve equation at the metal artifact position in the second curve.
Optionally, correcting the second curve equation with the first curve equation at the metal artifact position in the second curve includes:
When the number of the bulb tube sampling number values is M and the number of the radar sampling number values is N, M, N is a natural number and M is larger than N, a first ordinate value determined by a first abscissa under a first curve equation is used as a second ordinate value corresponding to a second abscissa under a second curve equation, wherein the second abscissa is M/N times of the first abscissa;
Traversing all second abscissa values corresponding to the metal artifact, and determining second ordinate values corresponding to all second abscissa values;
And fitting all the second ordinate values, and determining a second curve equation corresponding to the metal artifact.
Optionally, the number of the laser radars is multiple, and the multiple laser radars are distributed at positions corresponding to the bulb tube scanning.
According to a second aspect of the present disclosure, an apparatus for eliminating metal artifacts of CT images includes:
a determining unit configured to determine a position of a metal in the subject based on the pre-scan information;
a scanning unit for performing a tomographic scan of a subject site including the metal, the tomographic scan including a bulb tube X-ray scan and at least one lidar scan;
The device comprises an acquisition unit, a detection unit and a detection unit, wherein the acquisition unit is used for acquiring a ball tube sampling value and a radar sampling value which are corresponding to the same position of the detected part, and the ball tube sampling position corresponding to the ball tube sampling value has a one-to-one correspondence with the radar sampling position corresponding to the radar sampling value;
the generating unit is used for generating a first curve according to the radar sampling value and generating a second curve according to the bulb sampling value;
and the elimination unit is used for correcting the second curve based on the first curve and eliminating metal artifacts in the second curve.
According to a third aspect of the present disclosure, there is provided a computer readable storage medium having stored thereon a computer program which when executed by a processor implements a method of eliminating CT image metal artifacts as described in any of the above.
According to a fourth aspect of the present disclosure, there is provided an electronic device comprising: one or more processors; storage means for storing one or more programs that when executed by the one or more processors cause the one or more processors to implement a method of eliminating CT image metal artifacts as described in any preceding claim.
Compared with the prior art, the scheme of the embodiment of the disclosure has at least the following beneficial effects:
The present disclosure provides a method, apparatus, medium and electronic device for eliminating CT image metal artifacts. According to the characteristics of the radar sensor, the radar sensor and the bulb sensor are sampled simultaneously, the sampling positions of the radar sensor and the bulb sensor are completely registered, and the bulb sampling value comprising metal artifact information is corrected by means of first curve information generated by the radar sampling value, so that the metal artifact is eliminated on the basis of medical images generated by the corrected bulb sampling value.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description, serve to explain the principles of the disclosure. It will be apparent to those of ordinary skill in the art that the drawings in the following description are merely examples of the disclosure and that other drawings may be derived from them without undue effort. In the drawings:
FIG. 1 illustrates a flow chart of a method of eliminating CT image metal artifacts in an embodiment of the present disclosure;
fig. 2 is a schematic diagram illustrating a correspondence between a bulb and a radar in a CT machine according to an embodiment of the disclosure;
fig. 3 illustrates a schematic diagram of a bulb and radar data acquisition correspondence of an embodiment of the present disclosure;
FIG. 4 shows a first curvilinear schematic of an embodiment of the present disclosure;
FIG. 5 shows a second curvilinear schematic of an embodiment of the present disclosure;
FIG. 6 shows a block diagram of a unit of an apparatus for eliminating CT image metal artifacts in an embodiment of the present disclosure;
Fig. 7 shows a schematic diagram of an electronic device connection structure according to an embodiment of the present disclosure.
Detailed Description
For the purpose of promoting an understanding of the principles and advantages of the disclosure, reference will now be made in detail to the drawings, in which it is apparent that the embodiments described are only some, but not all embodiments of the disclosure. Based on the embodiments in this disclosure, all other embodiments that a person of ordinary skill in the art would obtain without making any inventive effort are within the scope of protection of this disclosure.
The terminology used in the embodiments of the disclosure is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. As used in this disclosure of embodiments and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise, the "plurality" generally includes at least two.
It should be understood that the term "and/or" as used herein is merely one relationship describing the association of the associated objects, meaning that there may be three relationships, e.g., a and/or B, may represent: a exists alone, A and B exist together, and B exists alone. In addition, the character "/" herein generally indicates that the front and rear associated objects are an "or" relationship.
It should be understood that although the terms first, second, third, etc. may be used in embodiments of the present disclosure, these descriptions should not be limited to these terms. These terms are only used to distinguish one from another. For example, a first may also be referred to as a second, and similarly, a second may also be referred to as a first, without departing from the scope of embodiments of the present disclosure.
The words "if", as used herein, may be interpreted as "at … …" or "at … …" or "in response to a determination" or "in response to a detection", depending on the context. Similarly, the phrase "if determined" or "if detected (stated condition or event)" may be interpreted as "when determined" or "in response to determination" or "when detected (stated condition or event)" or "in response to detection (stated condition or event), depending on the context.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a product or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such product or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a commodity or device comprising such elements.
Alternative embodiments of the present disclosure are described in detail below with reference to the drawings.
The embodiment of the disclosure provides a method for eliminating metal artifacts of a CT image, as shown in fig. 1, specifically comprising the following method steps:
step S102: determining a position of a metal in the subject based on the pre-scan information;
the position of the metal in the scanned object is determined according to the scanning data acquired by the conventional scanning by the CT machine, for example, the metal on the leg of the scanned object can be obtained by X-ray scanning, and then the approximate position of the metal can be determined.
Step S104: a tomographic scan is performed on a subject site containing the metal, the tomographic scan including a bulb X-ray scan and at least one lidar scan.
As shown in fig. 2, a tomographic scan is performed by a CT machine on a subject site including the metal, wherein the subject may be a human or animal body, and the subject site may be a part of the human or animal body including the metal determined in the above steps.
As shown in fig. 2, the CT machine includes a bulb X-ray emitter for scanning X-rays emitted from a subject, and an X-ray array sensor for receiving the X-rays after passing through the subject is provided on a side opposite to the bulb, and the X-ray array sensor uploads the received data to a data processing device to form a contour image of a scanned region. The CT machine also comprises at least one laser radar scanning device which is used for scanning the same detected part, and according to the laser radar reflection type scanning principle, the laser radar scanning image can not generate the artifact of the scanning image due to the existence of metal. After the CT machine rotates for tomographic scanning, the scanning positions of the laser radar and the scanning positions of the X-rays are in one-to-one correspondence, for example, the X-rays start from a 90-degree position and scan once every rotation, the laser radar starts from a 0-degree position and scans once every rotation, so that after the CT machine rotates 360 degrees, the X-rays and the laser radar acquire scanning data once at each position of 0-degree, 1-degree, 2-degree, … … -degree and 359-degree, a one-to-one corresponding scanning image is formed, and the scanning data of 360 scanning positions are acquired.
The working principle of the radar sensor is that the radar sensor emits pulse laser by taking laser as a signal source, the pulse laser strikes the surface of an object to cause scattering, a part of light waves can be reflected to a receiver of the radar sensor, the distance from the radar sensor to a target point can be obtained according to the laser ranging principle, and the radar sensor continuously scans the surface of a target object, so that sampling values of all target points on the surface of the target object can be obtained.
Step S106: and obtaining a ball tube sampling value and a radar sampling value which are corresponding to the same position of the detected part, wherein the ball tube sampling positions corresponding to the ball tube sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence.
The scanning process is tomographic scanning, i.e. the detected body is stationary, and the CT machine rotates and scans the detected part through the bulb tube and the laser radar to obtain scanning data. After the CT machine rotates for tomographic scanning, the scanning positions of the laser radar and the scanning positions of the X-rays are in one-to-one correspondence, for example, the X-rays start from a 90-degree position and scan once every rotation, the laser radar starts from a 0-degree position and scans once every rotation, so that after the CT machine rotates 360 degrees, the X-rays and the laser radar acquire scanning data once at each position of 0-degree, 1-degree, 2-degree, … … -degree and 359-degree, a one-to-one corresponding scanning image is formed, and the scanning data of 360 scanning positions are acquired. For scanning of each position, due to the principle of a hardware structure, the number of the X-ray sensor arrays is far greater than that of the laser radar sensors, so that scanning data obtained by the X-ray sensors are far greater than that of the laser radar sensors at each position, and as an example, the X-ray sensor arrays have 800 sensor units, can receive 800 scanning signals at one time, the laser radar sensors have 80 sensor units, can receive 80 scanning signals at one time, and then take the X-ray sensors obtain 800 data at one time, and the laser radar obtains 80 data for illustration. To illustrate this further, when data is acquired at the same location, for example, for a location of the same length L, the X-ray sensor equally divides L into 800 points to obtain 800 data, and the lidar equally divides L into 80 points to obtain 80 data, which is understood with reference to the schematic diagram of fig. 3.
Step S108: and generating a first curve according to the radar sampling value, and generating a second curve according to the bulb sampling value.
As described above, taking one of the positions as an example, for example, at the 0 degree position, the bulb tube obtains 800 scan data of the X-ray scan, the laser radar obtains 80 scan data, a first curve can be formed based on the 80 scan data obtained by the radar, and a second curve can be formed based on the 800 scan data obtained by the bulb tube. Wherein, the abscissa x1 of the first curve is 1, 2, … …,80, the ordinate y1 of the first curve is the measured value corresponding to each point, the abscissa x2 of the second curve is 1, 2, … …, 800, and the ordinate y2 of the second curve is the measured value corresponding to each point.
As an alternative embodiment, generating a first curve from the radar sample values comprises: carrying out normalization processing on the radar sampling value of each sampling position of the radar to obtain a normalization value of each sampling position of the radar; and carrying out fitting processing on the normalized values based on each sampling position of the radar, and generating a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
For example, at the 0 degree position, the corresponding radar sampling values of which the abscissa x1 of the first curve is 1, 2, … … and 80 are respectively between 1 ten thousand and 5 ten thousand, the maximum value is selected as a divisor (the sum of all the numbers can be selected as the divisor), the radar sampling value is taken as a dividend, and the quotient is the normalized value of the radar sampling position; for example, the radar sampling value is: 21000. 32000, 13600, 24800 and 42100, wherein the maximum value is 42100, the normalized value is: 0.499, 0.76, 0.323, 0.589, and 1.
Based on the normalized y-coordinate values described above, a normalized fit first curve is generated, for example, from 80 data points, a first curve equation, for example, y 1＝a1x1 2+b1x1+c1, can be derived, and an accurate value of a 1、b1、c1 can be derived, as shown in FIG. 4.
As an alternative embodiment, the generating a second curve according to the tube sampled values includes: normalizing the bulb tube sampling value of each sampling position of the bulb tube to obtain a normalized value of each sampling position of the bulb tube; and fitting the normalized values based on each sampling position of the bulb tube to generate a second curve, wherein the second curve comprises a recess formed by metal artifacts.
For the same reason, normalizing the ball tube sampling value of each sampling position of the ball tube at the position of 0 degrees to obtain a normalized value of each sampling position of the ball tube, for example, the ball tube sampling value of each ball tube sampling position is respectively between 0 and 255, the maximum value is selected as a divisor, the ball tube sampling value is selected as a divisor, and the quotient is the normalized value of the ball tube sampling position; for example, the bulb sampling value is: 210. 62, 53, 78 and 190, wherein the maximum value is 210, the normalized value obtained is: 1. 0.295, 0.252, 0.371, and 0.905.
And fitting the normalized value based on each sampling position of the bulb tube to generate a second curve. However, due to the presence of metal artifacts, the second curve has an irregular concave curve, as shown in fig. 5, for example, artifacts exist between 300-500 data points.
Step S110: and correcting the second curve based on the first curve, and eliminating metal artifacts in the second curve.
The correcting the second curve based on the first curve, eliminating metal artifacts in the second curve, includes: and correcting the second curve equation by the first curve equation at the metal artifact position in the second curve.
Correcting the second curve equation with the first curve equation at the metal artifact location in the second curve, comprising: when the number of the bulb tube sampling number values is M and the number of the radar sampling number values is N, M, N is a natural number and M is larger than N, a first ordinate value determined by a first abscissa under a first curve equation is used as a second ordinate value corresponding to a second abscissa under a second curve equation, wherein the second abscissa is M/N times of the first abscissa; traversing all second abscissa values corresponding to the metal artifact, and determining second ordinate values corresponding to all second abscissa values; and fitting all the second ordinate values, and determining a second curve equation corresponding to the metal artifact.
As described above, to eliminate artifacts between 300-500 data points, values of 300-500 data points may be traversed, and corresponding y values of 300-500 data points may be repaired, correcting the artifacts. As an example, the following is illustrated:
For the data of 400 positions corrected, for the data of X2 in the second curve equation of the bulb X-ray, corresponding to the position of X1 of the radar scanning curve being 40, substituting the data of x1=40 into the first curve equation to obtain the value of y1, traversing each value of 300-500 in turn corresponding to the value of y2 when x2=400 is obtained, and repairing the artifact value of 300-500 can be achieved corresponding to the value of y2 when x2=300-500.
As an alternative embodiment, the number of the laser radars is plural, and the plurality of the laser radars are distributed at positions corresponding to the bulb scan. Because the scanning data of one laser radar is far lower than the scanning data of one bulb, a mode of simultaneously scanning a plurality of laser radars can be adopted to increase the scanning data and obtain accurate comparison data. At this time, the positions of the multiple lidars need to be distributed on the scanning position of the bulb, for example, the multiple lidars are distributed on the positions of 0 degree, 1 degree, 2 degrees, … … degrees and 359 degrees, and optionally, 2-4 lidars are used to obtain data for correction.
The present disclosure provides a method of eliminating CT image metal artifacts. According to the characteristics of the radar sensor, the radar sensor and the bulb sensor are sampled simultaneously, the sampling positions of the radar sensor and the bulb sensor are completely registered, and the bulb sampling value comprising metal artifact information is corrected by means of first curve information generated by the radar sampling value, so that the metal artifact is eliminated on the basis of medical images generated by the corrected bulb sampling value.
As an alternative embodiment, determining the artifact value based on the second curve information and the first curve information comprises the steps of:
and comparing the normalized value of each bulb tube sampling position in the second curve information with the normalized value of the corresponding radar sampling position in the first curve information to obtain a comparison result.
In the step, the radar sampling positions and the bulb sampling positions have a one-to-one correspondence. For example, in the second curve information, the normalized value corresponding to the bulb sampling position a is a; in the first curve information, the radar sampling position B is the same as the bulb sampling position A in position, and the normalization value corresponding to the radar sampling position B is B; thus, a comparison result is obtained by comparing a with b, for example, calculating the difference between a and b.
Secondly, when the comparison result meets the preset artifact condition, determining a bulb tube sampling position corresponding to the comparison result.
The purpose of the comparison is to check the change rule of the normalized value of the bulb tube sampling position by taking the normalized value of the radar sampling position as a reference. For example, continuing the above example, comparing a with b, for example, calculating a difference between a and b, and when the difference is greater than a preset threshold (corresponding to the comparison result meeting a preset artifact condition), indicating that a curve represented by the a value in the second curve information and a curve represented by the b value in the first curve information change greatly, it can be determined that the bulb sampling value at the bulb sampling position is a value generating an image artifact.
Thirdly, determining a bulb tube sampling value corresponding to the bulb tube sampling position as the artifact value.
The embodiment of the disclosure uses the first curve information as a reference to identify the information associated with the artifact value in the second curve information, thereby finding the artifact value in the bulb tube sampling value.
As an alternative embodiment, correcting the artifact value in the bulb sample value based on the first curve information includes the following specific steps:
and firstly, determining a corresponding radar sampling position based on the bulb sampling position corresponding to the artifact value.
Because the bulb tube sampling positions and the radar sampling positions have a one-to-one correspondence, the corresponding radar sampling positions can be determined through the bulb tube sampling positions corresponding to the artifact values. It is understood that the position information related to the artifact values is found in all radar sampling locations. The radar sampling value corresponding to the position information can be used as a reference value for correcting the artifact value.
And secondly, acquiring a corresponding normalized value from the first curve information based on the radar sampling position.
Since the first curve information includes a correspondence of radar sampling positions to normalized values of radar sampling positions. Therefore, the corresponding normalized value can be obtained from the first curve information by the radar sampling position.
Thirdly, correcting the normalized value of the bulb tube sampling position corresponding to the current radar sampling position based on the normalized value, and generating a first normalized value of the bulb tube sampling position.
The normalized value of the bulb tube sampling position corresponding to the current radar sampling position is the normalized value between two virtual lines in the bulb tube acquisition image, and the normalized value corresponds to the artifact value.
Optionally, the correcting the normalized value of the bulb tube sampling position corresponding to the current radar sampling position based on the normalized value generates a first normalized value of the bulb tube sampling position, which includes the following specific steps:
First, determining that the first normalized value is equal to a normalized value of a current radar sampling position.
It is understood that the normalized value of the radar sampling location is replaced with the normalized value of the radar sampling location associated with the artifact value. Third curve information generated after the replacement of the second curve information becomes similar to the first curve information.
And secondly, carrying out inverse normalization processing on the first normalization value to generate a correction value of the artifact value.
The inverse normalization processing is inverse processing of normalization processing on the bulb tube sampling value, and the processing steps are opposite to those of the normalization processing. Since the second curve information includes the correspondence of the bulb sampling position and the normalized value of the bulb sampling position. The normalized value of the bulb tube sampling position irrelevant to the artifact value in the second curve information is the same as the original bulb tube sampling value after the inverse normalization processing; and the normalized value of the bulb tube sampling position related to the artifact value in the second curve information is different from the original bulb tube sampling value after the inverse normalization processing. For example, the bulb sample value is: 210. 62, 53, 78 and 190, the normalized values of the bulb sampling position in the second curve information are: 1. 0.295, 0.252, 0.371, and 0.905; normalized values of the bulb tube sampling positions irrelevant to the artifact value in the second curve information are 1 and 0.905, and after the inverse normalization processing is carried out, the obtained values are 210 and 190, which are the same as the original bulb tube sampling values; normalized values of the bulb sampling positions related to the artifact values in the second curve information are 0.295, 0.252 and 0.371, and first normalized values generated after correction are 0.973, 0.926 and 0.918, and values 204.33, 194.46 and 192.78 are obtained after inverse normalization processing, which are different from the original bulb sampling values 62, 53 and 78.
The medical image after artifact removal can be generated through the correction value of the artifact value and the bulb tube sampling value which is irrelevant to the artifact value, and diagnosis of a patient can be effectively assisted by a doctor.
The present disclosure provides a method of eliminating CT image metal artifacts. According to the characteristics of the radar sensor, the radar sensor and the bulb sensor are sampled simultaneously, the sampling positions of the radar sensor and the bulb sensor are completely registered, and the bulb sampling value comprising metal artifact information is corrected by means of first curve information generated by the radar sampling value, so that the metal artifact is eliminated on the basis of medical images generated by the corrected bulb sampling value.
The present embodiment provides a device for eliminating metal artifacts in CT images, as shown in fig. 6, which is configured to implement the method described in the foregoing embodiment, where the same features have the same technical effects, and are not described in detail, and includes:
A determining unit 602 for determining a position of a metal in the subject based on the pre-scan information;
A scanning unit 604 for performing a tomographic scan of a subject site including the metal, the tomographic scan including a bulb X-ray scan and at least one lidar scan;
An obtaining unit 606, configured to obtain a number of bulb tube sampling values and radar sampling values corresponding to the same position of the detected part, where the bulb tube sampling positions corresponding to the bulb tube sampling values have a one-to-one correspondence with the radar sampling positions corresponding to the radar sampling values;
The generating unit 608 is configured to generate a first curve according to the radar sampling value, and generate a second curve according to the bulb sampling value;
And an elimination unit 610, configured to correct the second curve based on the first curve, and eliminate metal artifacts in the second curve.
As an alternative embodiment, the generating unit 608 is configured to: carrying out normalization processing on the radar sampling value of each sampling position of the radar to obtain a normalization value of each sampling position of the radar; and carrying out fitting processing on the normalized values based on each sampling position of the radar, and generating a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
As an alternative embodiment, the generating unit 608 is configured to: normalizing the bulb tube sampling value of each sampling position of the bulb tube to obtain a normalized value of each sampling position of the bulb tube; and fitting the normalized values based on each sampling position of the bulb tube to generate a second curve, wherein the second curve comprises a recess formed by metal artifacts.
As an alternative embodiment, the eliminating unit 610 is configured to: and correcting the second curve equation by the first curve equation at the metal artifact position in the second curve.
Correcting the second curve equation with the first curve equation at the metal artifact location in the second curve, comprising: when the number of the bulb tube sampling number values is M and the number of the radar sampling number values is N, M, N is a natural number and M is larger than N, a first ordinate value determined by a first abscissa under a first curve equation is used as a second ordinate value corresponding to a second abscissa under a second curve equation, wherein the second abscissa is M/N times of the first abscissa; traversing all second abscissa values corresponding to the metal artifact, and determining second ordinate values corresponding to all second abscissa values; and fitting all the second ordinate values, and determining a second curve equation corresponding to the metal artifact.
As an alternative embodiment, the number of the laser radars is plural, and the plurality of the laser radars are distributed at positions corresponding to the bulb scan.
The present disclosure provides an apparatus for eliminating CT image metal artifacts. According to the characteristics of the radar sensor, the radar sensor and the bulb sensor are sampled simultaneously, the sampling positions of the radar sensor and the bulb sensor are completely registered, and the bulb sampling value comprising metal artifact information is corrected by means of first curve information generated by the radar sampling value, so that the metal artifact is eliminated on the basis of medical images generated by the corrected bulb sampling value.
The present embodiment provides an electronic device, including: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the one processor to enable the at least one processor to perform the method steps described in the embodiments above.
The disclosed embodiments provide a non-transitory computer storage medium storing computer executable instructions that perform the method steps described in the embodiments above.
Referring now to fig. 7, a schematic diagram of an electronic device suitable for use in implementing embodiments of the present disclosure is shown. The terminal devices in the embodiments of the present disclosure may include, but are not limited to, mobile terminals such as mobile phones, notebook computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), in-vehicle terminals (e.g., in-vehicle navigation terminals), and the like, and stationary terminals such as digital TVs, desktop computers, and the like. The electronic device shown in fig. 6 is merely an example and should not be construed to limit the functionality and scope of use of the disclosed embodiments.
As shown in fig. 7, the electronic device may include a processing means (e.g., a central processing unit, a graphics processor, etc.) 701, which may perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 702 or a program loaded from a storage means 708 into a Random Access Memory (RAM) 703. In the RAM 703, various programs and data necessary for the operation of the electronic apparatus are also stored. The processing device 701, the ROM 702, and the RAM 703 are connected to each other through a bus 704. An input/output (I/O) interface 705 is also connected to bus 704.
In general, the following devices may be connected to the I/O interface 705: input devices 706 including, for example, a touch screen, touchpad, keyboard, mouse, camera, microphone, accelerometer, gyroscope, and the like; an output device 707 including, for example, a Liquid Crystal Display (LCD), a speaker, a vibrator, and the like; storage 708 including, for example, magnetic tape, hard disk, etc.; and a communication device 709. The communication means 709 may allow the electronic device to communicate with other devices wirelessly or by wire to exchange data. While fig. 7 shows an electronic device having various means, it is to be understood that not all of the illustrated means are required to be implemented or provided. More or fewer devices may be implemented or provided instead.
In particular, according to embodiments of the present disclosure, the processes described above with reference to flowcharts may be implemented as computer software programs. For example, embodiments of the present disclosure include a computer program product comprising a computer program embodied on a non-transitory computer readable medium, the computer program comprising program code for performing the method shown in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network via communication device 709, or installed from storage 708, or installed from ROM 702. The above-described functions defined in the methods of the embodiments of the present disclosure are performed when the computer program is executed by the processing device 701.
It should be noted that the computer readable medium described in the present disclosure may be a computer readable signal medium or a computer readable storage medium, or any combination of the two. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples of the computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this disclosure, a computer-readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present disclosure, however, the computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, with the computer-readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: electrical wires, fiber optic cables, RF (radio frequency), and the like, or any suitable combination of the foregoing.
The computer readable medium may be contained in the electronic device; or may exist alone without being incorporated into the electronic device.
The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units involved in the embodiments of the present disclosure may be implemented by means of software, or may be implemented by means of hardware. Wherein the names of the units do not constitute a limitation of the units themselves in some cases.
In the context of this disclosure, a machine-readable medium may be a tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. The machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.
The units involved in the embodiments of the present disclosure may be implemented by means of software, or may be implemented by means of hardware. Wherein the names of the units do not constitute a limitation of the units themselves in some cases.
Claims (7)
1. A method of eliminating CT image metal artifacts, comprising:
determining a position of a metal in the subject based on the pre-scan information;
Performing a tomographic scan of a subject site containing the metal, the tomographic scan including a bulb X-ray scan and at least one lidar scan;
Obtaining a corresponding number of bulb tube sampling values and radar sampling values of the same position of the detected body part, wherein the bulb tube sampling positions corresponding to the bulb tube sampling values have a one-to-one correspondence with the radar sampling positions corresponding to the radar sampling values;
generating a first curve according to the radar sampling value, and generating a second curve according to the bulb sampling value;
correcting the second curve based on the first curve, and eliminating metal artifacts in the second curve;
the correcting the second curve based on the first curve, eliminating metal artifacts in the second curve, includes:
correcting a second curve equation with the first curve equation at the metal artifact position in the second curve;
correcting a second curve equation with a first curve equation at a metal artifact location in the second curve, comprising:
When the number of the bulb tube sampling number values is M and the number of the radar sampling number values is N, M, N is a natural number and M is larger than N, a first ordinate value determined by a first abscissa under a first curve equation is used as a second ordinate value corresponding to a second abscissa under a second curve equation, wherein the second abscissa is M/N times of the first abscissa;
traversing all second abscissa values corresponding to the metal artifact, and determining all second ordinate values corresponding to the second abscissa values;
And fitting all the second ordinate values, and determining a second curve equation corresponding to the metal artifact.
2. The method of claim 1, wherein the generating a first curve from the radar sample values comprises:
Carrying out normalization processing on the radar sampling value of each sampling position of the radar to obtain a normalization value of each sampling position of the radar;
And carrying out fitting processing on the normalized values based on each sampling position of the radar, and generating a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
3. The method of claim 2, wherein the generating a second curve from the bulb sample values comprises:
Normalizing the bulb tube sampling value of each sampling position of the bulb tube to obtain a normalized value of each sampling position of the bulb tube;
and fitting the normalized values based on each sampling position of the bulb tube to generate a second curve, wherein the second curve comprises a recess formed by metal artifacts.
4. The method of claim 1, wherein the number of lidars is a plurality, the plurality of lidars being distributed at locations corresponding to the bulb scan.
5. An apparatus for eliminating metal artifacts in CT images, comprising:
a determining unit configured to determine a position of a metal in the subject based on the pre-scan information;
a scanning unit for performing a tomographic scan of a subject site including the metal, the tomographic scan including a bulb tube X-ray scan and at least one lidar scan;
The device comprises an acquisition unit, a detection unit and a detection unit, wherein the acquisition unit is used for acquiring a corresponding number of bulb tube sampling values and radar sampling values at the same position of the detected body part, and the bulb tube sampling positions corresponding to the bulb tube sampling values have a one-to-one correspondence with the radar sampling positions corresponding to the radar sampling values;
the generating unit is used for generating a first curve according to the radar sampling value and generating a second curve according to the bulb sampling value;
An elimination unit configured to correct the second curve based on the first curve, and eliminate metal artifacts in the second curve; the correcting the second curve based on the first curve, eliminating metal artifacts in the second curve, includes: correcting a second curve equation with the first curve equation at the metal artifact position in the second curve;
Correcting a second curve equation with a first curve equation at a metal artifact location in the second curve, comprising: when the number of the bulb tube sampling number values is M and the number of the radar sampling number values is N, M, N is a natural number and M is larger than N, a first ordinate value determined by a first abscissa under a first curve equation is used as a second ordinate value corresponding to a second abscissa under a second curve equation, wherein the second abscissa is M/N times of the first abscissa; traversing all second abscissa values corresponding to the metal artifact, and determining all second ordinate values corresponding to the second abscissa values; and fitting all the second ordinate values, and determining a second curve equation corresponding to the metal artifact.
6. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the method according to any one of claims 1 to 4.
7. An electronic device, comprising:
one or more processors;
storage means for storing one or more programs which when executed by the one or more processors cause the one or more processors to implement the method of any of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110610275.7A CN113345043B (en) | 2021-06-01 | 2021-06-01 | Method, device, medium and electronic equipment for eliminating CT image metal artifact |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110610275.7A CN113345043B (en) | 2021-06-01 | 2021-06-01 | Method, device, medium and electronic equipment for eliminating CT image metal artifact |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113345043A CN113345043A (en) | 2021-09-03 |
| CN113345043B true CN113345043B (en) | 2024-05-17 |
Family
ID=77474515
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110610275.7A Active CN113345043B (en) | 2021-06-01 | 2021-06-01 | Method, device, medium and electronic equipment for eliminating CT image metal artifact |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113345043B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2024103412A1 (en) * | 2022-11-18 | 2024-05-23 | 上海联影医疗科技股份有限公司 | Metal artifact correction method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105528773A (en) * | 2015-12-31 | 2016-04-27 | 深圳市一体医疗科技有限公司 | Multi-modal microwave imaging method and system based on labeling method |
| CN105943071A (en) * | 2016-05-25 | 2016-09-21 | 厦门大学 | X-ray ct imaging system |
| CN107789001A (en) * | 2017-10-31 | 2018-03-13 | 上海联影医疗科技有限公司 | A kind of pendulum position method and system for image scanning |
| CN111084635A (en) * | 2020-01-02 | 2020-05-01 | 沈阳先进医疗设备技术孵化中心有限公司 | CT scanning method and device, CT equipment and CT system |
| CN111297386A (en) * | 2020-02-18 | 2020-06-19 | 苏州晟诺医疗科技有限公司 | CT data positioning method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7840053B2 (en) * | 2007-04-05 | 2010-11-23 | Liao Hstau Y | System and methods for tomography image reconstruction |
| US10922855B2 (en) * | 2017-11-30 | 2021-02-16 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining at least one artifact calibration coefficient |
-
2021
- 2021-06-01 CN CN202110610275.7A patent/CN113345043B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105528773A (en) * | 2015-12-31 | 2016-04-27 | 深圳市一体医疗科技有限公司 | Multi-modal microwave imaging method and system based on labeling method |
| CN105943071A (en) * | 2016-05-25 | 2016-09-21 | 厦门大学 | X-ray ct imaging system |
| CN107789001A (en) * | 2017-10-31 | 2018-03-13 | 上海联影医疗科技有限公司 | A kind of pendulum position method and system for image scanning |
| CN111084635A (en) * | 2020-01-02 | 2020-05-01 | 沈阳先进医疗设备技术孵化中心有限公司 | CT scanning method and device, CT equipment and CT system |
| CN111297386A (en) * | 2020-02-18 | 2020-06-19 | 苏州晟诺医疗科技有限公司 | CT data positioning method |
Non-Patent Citations (2)
| Title |
|---|
| 能谱CT成像算法及其在质子治疗中的应用研究;李斌;中国博士学位论文全文数据库 医药卫生科技辑;全文 * |
| 运用兆伏锥形束计算机断层成像校正千伏计算机断层成像中假牙的金属伪影;高留刚;孙鸿飞;林涛;眭建峰;倪昕晔;;生物医学工程学杂志(第05期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113345043A (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11756164B2 (en) | System and method for image correction | |
| CN109171793B (en) | Angle detection and correction method, device, equipment and medium | |
| US10339634B2 (en) | System and method for image reconstruction | |
| US9437017B2 (en) | Method and apparatus for metal artifact elimination in a medical image | |
| US9706972B1 (en) | Systems and methods for reconstruction of emission activity image | |
| US7801264B2 (en) | Method for calibrating a dual -spectral computed tomography (CT) system | |
| JP2001351091A (en) | Image processor, photographing device, system and method for image processing, and storage medium | |
| EP2662835B1 (en) | Method of reconstructing image and apparatus and tomography apparatus for performing the same | |
| CN110176010A (en) | A kind of image detecting method, device, equipment and storage medium | |
| CN113345043B (en) | Method, device, medium and electronic equipment for eliminating CT image metal artifact | |
| CN112037170A (en) | Method and device for detecting vascular stenosis and computer storage medium | |
| GB2552412A (en) | Systems and methods for detecting luggage in an imaging system | |
| CN112581554B (en) | CT imaging method, device, storage equipment and medical imaging system | |
| US20060066911A1 (en) | Edge detection and correcting system and method | |
| US20200240934A1 (en) | Tomography apparatus and controlling method for the same | |
| US10417795B2 (en) | Iterative reconstruction with system optics modeling using filters | |
| CN109493396B (en) | CT image display method, device, equipment and medium | |
| KR101958099B1 (en) | Attenuation correction method using time-of-flight information in positron emission tomography | |
| JP2006230904A (en) | Image processor and image processing method | |
| CN113129239B (en) | Image correction method, device, electronic equipment and storage medium | |
| EP4242965A1 (en) | Optimal path finding based spinal center line extraction | |
| TW201417769A (en) | A method for improving image quality and imaging system using the same | |
| CN107736895B (en) | Protocol parameter configuration method, device and terminal | |
| CN113781593B (en) | Four-dimensional CT image generation method, device, terminal equipment and storage medium | |
| JP7233211B2 (en) | Image processing device, image processing method and program |
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 |