CN113345043A - Method, device, medium and electronic equipment for eliminating metal artifacts of CT image - Google Patents

Method, device, medium and electronic equipment for eliminating metal artifacts of CT image Download PDF

Info

Publication number
CN113345043A
CN113345043A CN202110610275.7A CN202110610275A CN113345043A CN 113345043 A CN113345043 A CN 113345043A CN 202110610275 A CN202110610275 A CN 202110610275A CN 113345043 A CN113345043 A CN 113345043A
Authority
CN
China
Prior art keywords
curve
sampling
bulb
radar
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110610275.7A
Other languages
Chinese (zh)
Other versions
CN113345043B (en
Inventor
任毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Shengnuo Medical Technology Co ltd
Original Assignee
Suzhou Shengnuo Medical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Suzhou Shengnuo Medical Technology Co ltd filed Critical Suzhou Shengnuo Medical Technology Co ltd
Priority to CN202110610275.7A priority Critical patent/CN113345043B/en
Publication of CN113345043A publication Critical patent/CN113345043A/en
Application granted granted Critical
Publication of CN113345043B publication Critical patent/CN113345043B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10044Radar image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed 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 an electronic device for correcting a bulb sampling value. The method comprises the following steps: determining a location of a metal in the subject based on the pre-scan information; performing tomography on a metal-containing part of a subject, wherein the tomography comprises scanning by using a bulb X-ray and at least one laser radar; acquiring corresponding number of bulb tube sampling values and radar sampling values at the same position of a detected part; generating a first curve according to a radar sampling value, and generating a second curve according to the bulb sampling value; and correcting the second curve based on the first curve, eliminating the metal artifact in the second curve according to the characteristics of the radar sensor, simultaneously sampling the radar sensor and the bulb tube sensor, completely registering the sampling positions of the radar sensor and the bulb tube sensor, correcting the bulb tube sampling value comprising the metal artifact information by means of the first curve information generated by the radar sampling value, and eliminating the metal artifact in the medical image generated based on the corrected bulb tube sampling value.

Description

Method, device, medium and electronic equipment for eliminating metal artifacts of CT image
Technical Field
The present disclosure relates to the field of medical imaging, and in particular, to a method, an apparatus, a medium, and an electronic device for eliminating a metal artifact in a CT image.
Background
When a patient is subjected to Computed Tomography (CT) scanning, if metal (such as a prosthesis or a stent) exists in the patient, metal artifacts exist in a CT image, which affects diagnosis of a doctor.
At present, for a correction algorithm of 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 through front and back projection, and then data in the metal region is replaced with calibrated data.
However, this method of eliminating metal artifacts results in loss of resolution in the image. Therefore, the present disclosure provides a method for correcting a sampling value of a bulb, so as to solve one of the above technical problems.
Disclosure of Invention
The present disclosure is directed to a method, an apparatus, a medium, and an electronic device for eliminating a metal artifact in a CT image, which can solve at least one of the above-mentioned technical problems. The specific scheme is as follows:
according to a specific embodiment of the present disclosure, in a first aspect, the present disclosure provides a method for eliminating a metal artifact in a CT image, including:
determining a location of a metal in the subject based on the pre-scan information;
performing a tomography scan of the subject site containing the metal, the tomography scan including an X-ray scan using a bulb and at least one lidar scan;
acquiring bulb sampling values and radar sampling values of the detected part, wherein the number of the bulb sampling values corresponds to the same position of the detected part, and the bulb sampling positions corresponding to the bulb sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence relationship;
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 value includes:
normalizing the radar sampling value of each sampling position of the radar to obtain a normalized value of each sampling position of the radar;
and fitting the normalized value based on each sampling position of the radar to generate 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 sampling value of each sampling position of the bulb to obtain a normalized value of each sampling position of the bulb;
and fitting the normalized value based on each sampling position of the bulb to generate a second curve, wherein the second curve comprises a depression formed by metal artifacts.
Optionally, the correcting the second curve based on the first curve to eliminate metal artifacts in the second curve includes:
correcting the second curve equation with the first curve equation at the metal artifact location in the second curve.
Optionally, correcting the second curve equation with the first curve equation at the metal artifact location in the second curve includes:
when the number of the sampling quantity values of the bulb tubes is M and the number of the radar sampling values is N, M, N is a natural number, M is greater 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 artifacts, and determining second ordinate values corresponding to all the second abscissa values;
and fitting all the second ordinate values to determine a second curve equation corresponding to the metal artifact.
Optionally, the number of the laser radars is multiple, and the laser radars are distributed at positions corresponding to the scanning of the bulb tube.
According to a specific embodiment of the present disclosure, in a second aspect, the present disclosure provides an apparatus for eliminating metal artifacts in CT images, including:
a determination unit configured to determine a position of a metal in the subject based on the pre-scan information;
a scanning unit for performing tomography scanning on the part of the subject containing the metal, wherein the tomography scanning comprises scanning by using a bulb X-ray and at least one laser radar;
the acquisition unit is used for acquiring the bulb sampling values and the radar sampling values of the detected part, wherein the bulb sampling positions corresponding to the bulb sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence relationship;
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 eliminating unit is used for correcting the second curve based on the first curve and eliminating the metal artifacts in the second curve.
According to a third aspect, the present disclosure provides a computer-readable storage medium, on which a computer program is stored, which when executed by a processor implements the method for eliminating metal artifacts in CT images as described in any one of the above.
According to a fourth aspect thereof, the present disclosure provides an electronic device, comprising: one or more processors; a storage device 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 for eliminating metal artifacts in CT images as described in any one of the above.
Compared with the prior art, the scheme of the embodiment of the disclosure at least has the following beneficial effects:
the disclosure provides a method, a device, a medium and an electronic device for eliminating metal artifacts of CT images. According to the characteristics of the radar sensor, the radar sensor and the bulb tube sensor are simultaneously sampled, the sampling positions of the radar sensor and the bulb tube sensor are completely registered, the bulb tube sampling value containing metal artifact information is corrected by means of first curve information generated by the radar sampling value, and the metal artifact is eliminated on the basis of a medical image generated by the corrected bulb tube sampling value.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and together with the description, serve to explain the principles of the disclosure. It is to be understood that the drawings in the following description are merely exemplary of the disclosure, and that other drawings may be derived from those drawings by one of ordinary skill in the art without the exercise of inventive faculty. In the drawings:
FIG. 1 shows a flowchart of a method for eliminating metal artifacts in CT images according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram illustrating a correspondence relationship between a bulb and a radar in the CT machine according to the embodiment of the disclosure;
FIG. 3 shows a schematic diagram of a bulb and radar data acquisition correspondence in an embodiment of the disclosure;
FIG. 4 illustrates a first graphical schematic of an embodiment of the present disclosure;
FIG. 5 illustrates a second graphical schematic of an embodiment of the present disclosure;
FIG. 6 shows a block diagram of a device for eliminating metal artifacts in CT images according to an embodiment of the present disclosure;
fig. 7 shows an electronic device connection structure schematic according to an embodiment of the present disclosure.
Detailed Description
To make the objects, technical solutions and advantages of the present disclosure clearer, the present disclosure will be described in further detail with reference to the accompanying drawings, and it is apparent that the described embodiments are only a part of the embodiments of the present disclosure, rather than all embodiments. All other embodiments, which can be derived by one of ordinary skill in the art from the embodiments disclosed herein without making any creative effort, shall fall within the scope of protection of the present disclosure.
The terminology used in the embodiments of the present disclosure is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. As used in the disclosed 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, and "a plurality" typically includes at least two.
It should be understood that the term "and/or" as used herein is merely one type of association that describes an associated object, meaning that three relationships may exist, e.g., a and/or B may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
It should be understood that although the terms first, second, third, etc. may be used in the embodiments of the present disclosure, these descriptions should not be limited to these terms. These terms are only used to distinguish one description from another. For example, a first could also be termed a second, and, similarly, a second could also be termed 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 phrases "if determined" or "if detected (a stated condition or event)" may be interpreted as "when determined" or "in response to a determination" or "when detected (a stated condition or event)" or "in response to a detection (a stated condition or event)", depending on the context.
It is also noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that an article 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 article or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in the article or device in which the element is included.
Alternative embodiments of the present disclosure are described in detail below with reference to the accompanying drawings.
The embodiment of the present disclosure provides a method for eliminating a metal artifact of a CT image, as shown in fig. 1, specifically including the following method steps:
step S102: determining a location of a metal in the subject based on the pre-scan information;
the conventional scanning is carried out by a CT machine, the position of metal in the scanned object is determined according to scanning data acquired by the conventional scanning, for example, the metal on the leg of the scanned object can be obtained by X-ray scanning, and the approximate position of the metal can be determined.
Step S104: and carrying out tomography scanning on the part of the object containing the metal, wherein the tomography scanning comprises scanning by using a bulb X-ray and at least one laser radar.
As shown in fig. 2, a CT machine is used to perform tomography on a part of a subject containing the metal, wherein the subject may be a human body or an animal body, and the part of the subject may be a part of the human body or the animal body, and the part includes the metal determined in the above steps.
As shown in fig. 2, the CT machine includes a bulb X-ray emitter for emitting X-rays to scan a subject, an X-ray array sensor is disposed on a side opposite to the bulb for receiving the X-rays after passing through the subject, and the X-ray array sensor uploads the received data to a data processing device to form a contour image of a scanned portion. 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 tomography, the scanning position of the lidar corresponds to the scanning position of the X-ray one-to-one, for example, the X-ray starts from a 90-degree position, and the lidar starts from a 0-degree position and scans once every one-degree rotation, so that after the CT machine rotates for 360 degrees, the X-ray and the lidar acquire scanning data once at each of 0-degree, 1-degree, 2-degree, … …, and 359 degrees, and form one-to-one corresponding scanning images, and the scanning data of 360 scanning positions are acquired in total.
The radar sensor adopts the working principle that laser is used as a signal source, the radar sensor emits pulse laser, 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 the target object to obtain the sampling values of all target points on the surface of the target object.
Step S106: and acquiring the bulb sampling values and the radar sampling values of the detected part, wherein the bulb sampling positions corresponding to the bulb sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence relationship.
The scanning process is tomography, namely the detected object is not moved, and the CT machine rotates, and simultaneously, the detected part is scanned through the bulb tube and the laser radar, and scanning data is obtained. After the CT machine rotates for tomography, the scanning position of the lidar corresponds to the scanning position of the X-ray one-to-one, for example, the X-ray starts from a 90-degree position, and the lidar starts from a 0-degree position and scans once every one-degree rotation, so that after the CT machine rotates for 360 degrees, the X-ray and the lidar acquire scanning data once at each of 0-degree, 1-degree, 2-degree, … …, and 359 degrees, and form one-to-one corresponding scanning images, and the scanning data of 360 scanning positions are acquired in total. For the scanning of each position, because the principle of a hardware structure, the number of the X-ray sensor arrays is far greater than that of the laser radar sensors, therefore, the scanning data obtained by the X-ray sensors at each position is far greater than that of the laser radar sensors, as an example, the X-ray sensor arrays are provided with 800 sensor units which can receive 800 scanning signals at one time, the laser radar sensors are provided with 80 sensor units which can receive 80 scanning signals at one time, the X-ray sensors at each scanning time are used for obtaining 800 data, and the laser radar obtains 80 data for example for explanation. To further illustrate, when data is collected at the same position, for example, for a position with 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, as will be understood with reference to 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 position as an example, for example, at the 0 degree position, the bulb obtains 800 scan data of X-ray scan, the lidar obtains 80 scan data, the 80 scan data obtained based on the radar can form a first curve, and the 800 scan data obtained based on the bulb can form a second curve. 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: normalizing the radar sampling value of each sampling position of the radar to obtain a normalized value of each sampling position of the radar; and fitting the normalized value based on each sampling position of the radar to generate 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 with the abscissa x1 of the first curve being 1, 2, … …, 80 are respectively between 1 ten thousand and 5 ten thousand, the maximum value is selected as a divisor (the sum of all numbers can also be selected as a divisor), the radar sampling value is used as a dividend, and the obtained quotient is the normalized value of the radar sampling position; for example, the radar samples are: 21000. 32000, 13600, 24800 and 42100, wherein the maximum value is 42100, the normalized values obtained are: 0.499, 0.76, 0.323, 0.589 and 1.
Based on the normalized y coordinate values, a normalized fit first curve is generated, e.g., from 80 data points, to arrive at a first curve equation, e.g., y1=a1x1 2+b1x1+c1And can obtain accurate a1、b1、c1The numerical values of (a) are shown in FIG. 4.
As an optional implementation, the generating a second curve according to the bulb sampled value includes: normalizing the bulb sampling value of each sampling position of the bulb to obtain a normalized value of each sampling position of the bulb; and fitting the normalized value based on each sampling position of the bulb to generate a second curve, wherein the second curve comprises a depression formed by metal artifacts.
For the same reason, at the 0-degree position, performing normalization processing on the bulb sampling value of each sampling position of the bulb to obtain a normalization value of each sampling position of the bulb, for example, the bulb sampling value of each sampling position of the bulb is respectively between 0 and 255, selecting the maximum value as a divisor, the bulb sampling value as a dividend, and obtaining a quotient which is the normalization value of the sampling position of the bulb; for example, the bulb sampling value is: 210. 62, 53, 78 and 190, where 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, the artifacts exist between the 300 th and 500 th 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, comprising: correcting the second curve equation with the first curve equation at the metal artifact location in the second curve.
Correcting the second curve equation with the first curve equation at metal artifact locations in the second curve, comprising: when the number of the sampling quantity values of the bulb tubes is M and the number of the radar sampling values is N, M, N is a natural number, M is greater 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 artifacts, and determining second ordinate values corresponding to all the second abscissa values; and fitting all the second ordinate values to determine a second curve equation corresponding to the metal artifact.
As described above, in order to eliminate the artifacts between the data points 300-. As an example, the following is illustrated:
describing by correcting data of 400 positions, for data of 400X 2 in the second curve equation of the tube X-ray, corresponding to the position where the radar scan curve X1 is 40, the data of 40X 1 is substituted into the first curve equation to obtain the value of y1, so as to obtain the value of y2 when 400X 2 is obtained, and sequentially traverse each value of 300-.
In an alternative embodiment, the number of the lidar is multiple, and the plurality of the lidar is distributed at positions corresponding to the scanning of the bulb tube. Because the scanning data of one laser radar is far lower than that of one bulb tube, 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 plurality of lidar need to be distributed on the scanning position of the bulb, for example, the plurality of lidar are distributed at positions of 0 degree, 1 degree, 2 degrees, … … degrees, and 359 degrees, and optionally, 2 to 4 lidar obtained data are used for correction.
The present disclosure provides a method for eliminating metal artifacts in CT images. According to the characteristics of the radar sensor, the radar sensor and the bulb tube sensor are simultaneously sampled, the sampling positions of the radar sensor and the bulb tube sensor are completely registered, the bulb tube sampling value containing metal artifact information is corrected by means of first curve information generated by the radar sampling value, and the metal artifact is eliminated on the basis of a medical image generated by the corrected bulb tube sampling value.
As an optional implementation, determining the artifact value based on the second curve information and the first curve information includes the following steps:
firstly, comparing the normalization value of each bulb tube sampling position in the second curve information with the normalization value of the corresponding radar sampling position in the first curve information to obtain a comparison result.
In the step, radar sampling positions and bulb tube sampling positions have 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, because the radar sampling position B is the same as the position of the bulb sampling position A, and the normalization value corresponding to the radar sampling position B is B; therefore, a comparison result is obtained by comparing a with b, such as calculating the difference between a and b.
And secondly, when the comparison result meets a preset artifact condition, determining a bulb tube sampling position corresponding to the comparison result.
The comparison is performed to check the variation rule of the normalized value of the sampling position of the bulb by taking the normalized value of the sampling position of the radar as a reference. For example, continuing with 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 satisfying a preset artifact condition), indicating that the curve represented by the value a in the second curve information and the curve represented by the value b in the first curve information have a large change, the change can determine that the tube sample value at the tube sample position is a value that generates an image artifact.
And thirdly, determining a bulb sampling value corresponding to the bulb sampling position as the artifact value.
The embodiment of the disclosure uses the first curve information as a reference for identifying information associated with the artifact value in the second curve information, so as to find the artifact value in the bulb sampling value.
As an optional implementation, the step of correcting the artifact value in the tube sample value based on the first curve information includes the following specific steps:
firstly, a corresponding radar sampling position is determined based on a bulb tube sampling position corresponding to the artifact value.
Because the bulb sampling position and the radar sampling position have a one-to-one correspondence relationship, the corresponding radar sampling position can be determined according to the bulb sampling position corresponding to the artifact value. It can be understood that the position information related to the artifact value is found in all radar sample positions. 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 normalization value from the first curve information based on the radar sampling position.
Since the first curve information includes the correspondence relationship of the radar sampling position and the normalized value of the radar sampling position. Therefore, the corresponding normalized value can be acquired from the first curve information by the radar sampling position.
Thirdly, correcting the normalization value of the bulb tube sampling position corresponding to the current radar sampling position based on the normalization value, and generating a first normalization value of the bulb tube sampling position.
The normalization value of the bulb tube sampling position corresponding to the current radar sampling position is the normalization value between two virtual lines in the bulb tube acquisition image, and the normalization value corresponds to the artifact value.
Optionally, the step of correcting the normalization value of the tube sampling position corresponding to the current radar sampling position based on the normalization value to generate a first normalization value of the tube sampling position includes the following specific steps:
first, it is determined that the first normalization value is equal to the normalization value of the current radar sampling position.
It is understood that the normalized value of the radar sample position replaces the normalized value of the radar sample position associated with the artifact value. Third curve information generated after the replacement of the second curve information is made 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 the inverse processing of normalization processing on the sampling value of the bulb, and the processing steps are opposite to those of the normalization processing. The second curve information comprises the corresponding relation between the bulb sampling position and the normalization value of the bulb sampling position. The normalization value of the bulb tube sampling position irrelevant to the artifact value in the second curve information is subjected to inverse normalization processing, and the obtained value is the same as the original bulb tube sampling value; and after the normalization value of the bulb sampling position related to the artifact value in the second curve information is subjected to inverse normalization processing, the obtained value is different from the original bulb sampling value. For example, the bulb samples are: 210. 62, 53, 78 and 190, and the normalized value of the bulb sampling position in the second curve information is: 1. 0.295, 0.252, 0.371, and 0.905; normalization values of the bulb sampling positions irrelevant to the artifact values in the second curve information are 1 and 0.905, and after inverse normalization processing is carried out, obtained values are 210 and 190 and are the same as the original bulb sampling values; the normalized values of the tube sample positions associated with the artifact values in the second curve information are 0.295, 0.252, and 0.371, the first normalized values generated after correction are 0.973, 0.926, and 0.918, and after the inverse normalization processing, the obtained values are 204.33, 194.46, and 192.78, which are different from the original tube sample values 62, 53, and 78.
The medical image without the artifact can be generated through the correction value of the artifact value and the bulb sampling value irrelevant to the artifact value, and a doctor can be effectively assisted to diagnose a patient.
The present disclosure provides a method for eliminating metal artifacts in CT images. According to the characteristics of the radar sensor, the radar sensor and the bulb tube sensor are simultaneously sampled, the sampling positions of the radar sensor and the bulb tube sensor are completely registered, the bulb tube sampling value containing metal artifact information is corrected by means of first curve information generated by the radar sampling value, and the metal artifact is eliminated on the basis of a medical image generated by the corrected bulb tube sampling value.
The present embodiment provides a device for eliminating metal artifacts in CT images, as shown in fig. 6, which is used to implement the method described in the above embodiments, and the same features have the same technical effects, which are not described herein again, and the method includes:
a determination 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 tomography scan of the subject part containing the metal, the tomography scan including an X-ray scan using a bulb and at least one lidar scan;
an obtaining unit 606, configured to obtain a number of bulb sampling values and radar sampling values corresponding to the same position of the detected part, where a one-to-one correspondence relationship exists between a bulb sampling position corresponding to the bulb sampling value and a radar sampling position corresponding to the radar sampling value;
a generating unit 608, configured to generate a first curve according to the radar sampling value, and generate a second curve according to the bulb sampling value;
a removing unit 610, configured to correct the second curve based on the first curve, and remove metal artifacts in the second curve.
As an optional implementation, the generating unit 608 is configured to: normalizing the radar sampling value of each sampling position of the radar to obtain a normalized value of each sampling position of the radar; and fitting the normalized value based on each sampling position of the radar to generate a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
As an optional implementation, the generating unit 608 is configured to: normalizing the bulb sampling value of each sampling position of the bulb to obtain a normalized value of each sampling position of the bulb; and fitting the normalized value based on each sampling position of the bulb to generate a second curve, wherein the second curve comprises a depression formed by metal artifacts.
As an optional implementation, the eliminating unit 610 is configured to: correcting the second curve equation with the first curve equation at the metal artifact location in the second curve.
Correcting the second curve equation with the first curve equation at metal artifact locations in the second curve, comprising: when the number of the sampling quantity values of the bulb tubes is M and the number of the radar sampling values is N, M, N is a natural number, M is greater 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 artifacts, and determining second ordinate values corresponding to all the second abscissa values; and fitting all the second ordinate values to determine a second curve equation corresponding to the metal artifact.
In an alternative embodiment, the number of the lidar is multiple, and the plurality of the lidar is distributed at positions corresponding to the scanning of the bulb tube.
The present disclosure provides a device for eliminating metal artifacts in CT images. According to the characteristics of the radar sensor, the radar sensor and the bulb tube sensor are simultaneously sampled, the sampling positions of the radar sensor and the bulb tube sensor are completely registered, the bulb tube sampling value containing metal artifact information is corrected by means of first curve information generated by the radar sampling value, and the metal artifact is eliminated on the basis of a medical image generated by the corrected bulb tube sampling value.
This embodiment provides an electronic device, electronic device includes: 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 cause the at least one processor to perform the method steps of the above embodiments.
The disclosed embodiments provide a non-volatile computer storage medium having stored thereon computer-executable instructions that may perform the method steps as described in the embodiments above.
Referring now to FIG. 7, shown is a schematic diagram of an electronic device suitable for use in implementing embodiments of the present disclosure. The terminal device in the embodiments of the present disclosure may include, but is not limited to, a mobile terminal such as a mobile phone, a notebook computer, a digital broadcast receiver, a PDA (personal digital assistant), a PAD (tablet computer), a PMP (portable multimedia player), a vehicle terminal (e.g., a car navigation terminal), and the like, and a stationary terminal such as a digital TV, a desktop computer, and the like. The electronic device shown in fig. 6 is only an example, and should not bring any limitation to the functions and the scope of use of the embodiments of the present disclosure.
As shown in fig. 7, the electronic device may include a processing device (e.g., central processing unit, 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 device 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 by a bus 704. An input/output (I/O) interface 705 is also connected to bus 704.
Generally, the following devices may be connected to the I/O interface 705: input devices 706 including, for example, a touch screen, touch pad, keyboard, mouse, camera, microphone, accelerometer, gyroscope, etc.; 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 device 709 may allow the electronic device to communicate wirelessly or by wire with other devices to exchange data. While fig. 7 illustrates an electronic device having various means, it is to be understood that not all illustrated means are required to be implemented or provided. More or fewer devices may alternatively be implemented or provided.
In particular, according to an embodiment of the present disclosure, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, embodiments of the present disclosure include a computer program product comprising a computer program carried on a non-transitory computer readable medium, the computer program containing program code for performing the method illustrated by the flow chart. In such embodiments, the computer program may be downloaded and installed from a network via the communication means 709, or may be installed from the storage means 708, or may be installed from the ROM 702. The computer program, when executed by the processing device 701, performs the above-described functions defined in the methods of the embodiments of the present disclosure.
It should be noted that the computer readable medium in the present disclosure can be a computer readable signal medium or a computer readable storage medium or any combination of the two. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination 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 present 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 contrast, in the present disclosure, a computer readable signal medium may comprise a propagated data signal with computer readable program code embodied therein, either in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. 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, optical cables, RF (radio frequency), etc., or any suitable combination of the foregoing.
The computer readable medium may be embodied in the electronic device; or may exist separately without being assembled into the electronic device.
The flowchart 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 described in the embodiments of the present disclosure may be implemented by software or hardware. Where the name of an element does not in some cases constitute a limitation on the element itself.
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. A 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 described in the embodiments of the present disclosure may be implemented by software or hardware. Where the name of an element does not in some cases constitute a limitation on the element itself.

Claims (9)

1. A method for eliminating metal artifacts in CT images, comprising:
determining a location of a metal in the subject based on the pre-scan information;
performing a tomography scan of the subject site containing the metal, the tomography scan including an X-ray scan using a bulb and at least one lidar scan;
acquiring bulb sampling values and radar sampling values of the detected part, wherein the number of the bulb sampling values corresponds to the same position of the detected part, and the bulb sampling positions corresponding to the bulb sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence relationship;
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.
2. The method of claim 1, wherein generating a first curve from the radar sample values comprises:
normalizing the radar sampling value of each sampling position of the radar to obtain a normalized value of each sampling position of the radar;
and fitting the normalized value based on each sampling position of the radar to generate a first curve, wherein the first curve comprises a first coefficient meeting a first curve equation.
3. The method of claim 2, wherein generating a second curve from the bulb sample values comprises:
normalizing the bulb sampling value of each sampling position of the bulb to obtain a normalized value of each sampling position of the bulb;
and fitting the normalized value based on each sampling position of the bulb to generate a second curve, wherein the second curve comprises a depression formed by metal artifacts.
4. The method of claim 3, wherein said correcting said second curve based on said first curve, eliminating metal artifacts in said second curve, comprises:
correcting the second curve equation with the first curve equation at the metal artifact location in the second curve.
5. The method of claim 4, wherein correcting the second curve equation with the first curve equation at metal artifact locations in the second curve comprises:
when the number of the sampling quantity values of the bulb tubes is M and the number of the radar sampling values is N, M, N is a natural number, M is greater 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 artifacts, and determining second ordinate values corresponding to all the second abscissa values;
and fitting all the second ordinate values to determine a second curve equation corresponding to the metal artifact.
6. The method of claim 1, wherein the number of lidar is multiple, and wherein the multiple lidar is distributed at locations corresponding to a scanning of a bulb.
7. An apparatus for eliminating metal artifacts in CT images, comprising:
a determination unit configured to determine a position of a metal in the subject based on the pre-scan information;
a scanning unit for performing tomography scanning on the part of the subject containing the metal, wherein the tomography scanning comprises scanning by using a bulb X-ray and at least one laser radar;
the acquisition unit is used for acquiring the bulb sampling values and the radar sampling values of the detected part, wherein the bulb sampling positions corresponding to the bulb sampling values and the radar sampling positions corresponding to the radar sampling values have a one-to-one correspondence relationship;
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 eliminating unit is used for correcting the second curve based on the first curve and eliminating the metal artifacts in the second curve.
8. A computer-readable storage medium, on which a computer program is stored, which program, when being executed by a processor, carries out the method according to any one of claims 1 to 6.
9. 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 carry out the method of any one of claims 1 to 6.
CN202110610275.7A 2021-06-01 2021-06-01 Method, device, medium and electronic equipment for eliminating CT image metal artifact Active CN113345043B (en)

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 true CN113345043A (en) 2021-09-03
CN113345043B 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)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080247502A1 (en) * 2007-04-05 2008-10-09 Liao Hstau Y System and methods for tomography image reconstruction
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
US20190206096A1 (en) * 2017-11-30 2019-07-04 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining at least one artifact calibration coefficient
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080247502A1 (en) * 2007-04-05 2008-10-09 Liao Hstau Y System and methods for tomography image reconstruction
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
US20190206096A1 (en) * 2017-11-30 2019-07-04 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining at least one artifact calibration coefficient
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)

* Cited by examiner, † Cited by third party
Title
李斌: "能谱CT成像算法及其在质子治疗中的应用研究", 中国博士学位论文全文数据库 医药卫生科技辑 *
高留刚;孙鸿飞;林涛;眭建峰;倪昕晔;: "运用兆伏锥形束计算机断层成像校正千伏计算机断层成像中假牙的金属伪影", 生物医学工程学杂志, no. 05 *

Also Published As

Publication number Publication date
CN113345043B (en) 2024-05-17

Similar Documents

Publication Publication Date Title
KR102560911B1 (en) Image processing apparatus, image processing method, and storage medium
US20170301066A1 (en) System and method for image correction
JP4792463B2 (en) System and method for correction of temporal artifacts in tomographic images
US10607099B2 (en) Methods and systems for imaging performance analysis
CN108292428A (en) The system and method for image reconstruction
EP2433261A2 (en) Continuous time-of-flight scatter simulation method
CN112037170B (en) Method and device for detecting vascular stenosis and computer storage medium
JP2015167613A (en) Apparatus, method and program for processing radiation image
CN112581554B (en) CT imaging method, device, storage equipment and medical imaging system
EP3349655B1 (en) Tomography apparatus and controlling method for the same
US10231688B2 (en) Radiation image analysis device, method, and program
CN113345043B (en) Method, device, medium and electronic equipment for eliminating CT image metal artifact
US10417795B2 (en) Iterative reconstruction with system optics modeling using filters
CN112006713B (en) Dynamic reconstruction, scan object positioning, count rate display method and PET scanning device
CN109493396B (en) CT image display method, device, equipment and medium
CN112052885B (en) Image processing method, device, equipment and PET-CT system
CN109613593B (en) Gamma photon action position sampling method, gamma photon action position sampling device, electronic equipment and storage medium
KR101958099B1 (en) Attenuation correction method using time-of-flight information in positron emission tomography
CN112562030A (en) Image reconstruction method and device and electronic equipment
CN112767347A (en) Image registration method and device, electronic equipment and storage medium
JPH1048341A (en) Imaging method and imaging system
CN113379862A (en) Method, apparatus, medium, and electronic device for obtaining scan parameter value
CN112598760B (en) Image truncation artifact correction method and system
CN107736895B (en) Protocol parameter configuration method, device and terminal
CN109727297B (en) Medical image reconstruction method, system, readable storage medium and apparatus

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