CN117788680A - Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium - Google Patents

Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium Download PDF

Info

Publication number
CN117788680A
CN117788680A CN202211158589.9A CN202211158589A CN117788680A CN 117788680 A CN117788680 A CN 117788680A CN 202211158589 A CN202211158589 A CN 202211158589A CN 117788680 A CN117788680 A CN 117788680A
Authority
CN
China
Prior art keywords
blood vessel
pixel
body layer
view angle
projection view
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.)
Pending
Application number
CN202211158589.9A
Other languages
Chinese (zh)
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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN202211158589.9A priority Critical patent/CN117788680A/en
Publication of CN117788680A publication Critical patent/CN117788680A/en
Pending legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The embodiment of the invention discloses a blood vessel image processing method, a device, a computer device and a storage medium, wherein the method comprises the following steps: determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer aiming at a projection view angle of each projection moment in a projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer; determining intersections of the target virtual rays and the current body layer, and determining pixel values of corresponding vascular pixel points of the current projection view angle in the current body layer according to types of the intersections and corresponding relations of types and pixel determination rules; and determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the corresponding blood vessel pixel points of each projection view angle in each body layer. Solves the problem that the time resolution of CT blood vessel images can not be improved in the prior art.

Description

Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium
Technical Field
The embodiment of the invention relates to the field of medical image processing, in particular to a blood vessel image processing method, a blood vessel image processing device, computer equipment and a storage medium.
Background
The current method for improving the time resolution and the dynamic imaging quantitative accuracy of CT (Computed Tomography, namely, electronic computed tomography) images mainly comprises the following steps: 1) A time decay curve for each pixel is fitted using a time basis function. Such methods first decompose the corresponding time decay curve for each pixel onto a known basis function, such as Gaussian, gamma-variate, linear or higher order spline, with the decomposition coefficients being the only unknowns. The number of unknowns to be estimated can be greatly reduced by introducing the known basis function, so that the pathological problem can be well restrained. The fitted time-curve function has a high time sampling rate. 2) Deconvolution in the time domain. The time decay curve corresponding to each pixel in the method can be modeled as a convolution of the true value curve with a known convolution kernel. The estimated truth curve can be obtained by solving a deconvolution problem. The truth curve has a high temporal sampling rate. 3) The local time decay curve is calibrated using the projection data. The method first models a time decay curve of a region of interest and then locally calibrates the time decay curve.
In summary, for various CT images, such as CT blood vessel images, the prior art can only increase the time sampling rate and cannot increase the time resolution.
Disclosure of Invention
The embodiment of the invention provides a blood vessel image processing method, a blood vessel image processing device, computer equipment and a storage medium, which solve the problem that the time resolution of CT blood vessel images cannot be improved in the prior art.
In a first aspect, an embodiment of the present invention provides a blood vessel image processing method, including:
acquiring a three-dimensional blood vessel image of a scanning part, deduction projection data for creating the three-dimensional blood vessel image, and a projection view angle sequence corresponding to the deduction projection data;
determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer according to a projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
determining intersection sets of the target virtual rays and the current body layer, and determining pixel values of all blood vessel pixel points corresponding to the current projection view angle in the current body layer according to the type of each intersection set and the corresponding relation of the type and the pixel determining rule, wherein the number of uncorrected blood vessel pixel points included in the intersection sets of different types is different;
And determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the corresponding blood vessel pixel points of each projection view angle in each body layer.
In a second aspect, an embodiment of the present invention further provides a blood vessel image processing apparatus, including:
the acquisition module is used for acquiring a three-dimensional blood vessel image of a scanning part, deduction projection data used for creating the three-dimensional blood vessel image and a projection view angle sequence corresponding to the deduction projection data;
the virtual light module is used for determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual light corresponding to the current body layer according to the projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual light is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
the pixel value determining module is used for determining intersections of the target virtual rays and the current body layer, and determining pixel values of blood vessel pixel points corresponding to the current projection view angle in the current body layer according to the types of the intersections and the corresponding relation between the types and pixel determining rules, wherein the number of uncorrected blood vessel pixel points included in the intersections of different types is different;
And the four-dimensional image determining module is used for determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the blood vessel pixel points of each body layer under each projection view angle.
In a third aspect, an embodiment of the present invention further provides a computer apparatus, including:
one or more processors;
a storage means for storing one or more programs;
the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the vascular image processing method as described in any of the embodiments.
In a fourth aspect, embodiments of the present invention also provide a storage medium containing computer-executable instructions, which when executed by a computer processor, are configured to perform the vascular image processing method of any of the embodiments.
Compared with the prior art, the technical scheme of the blood vessel image processing method provided by the embodiment is that for each projection view angle of the current body layer, the pixel value of each blood vessel pixel point corresponding to each projection view angle in the current body layer is determined through the type of intersection of the target virtual light and the current body layer and the corresponding relation of the type and the pixel determining rule, and the pixel value of each blood vessel pixel point corresponding to each projection view angle in each body layer is determined by analogy, so that a four-dimensional blood vessel image is determined; when the four-dimensional blood vessel image is displayed, the pixel value of each blood vessel pixel point in each body layer can be changed along with the switching of the projection view angle, so that the technical aim of improving the time resolution of the four-dimensional blood vessel image to the time resolution of the projection view angle is fulfilled, and the four-dimensional blood vessel image with the improved time resolution can provide more detailed blood vessel information for clinical diagnosis.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a blood vessel image processing method provided by an embodiment of the present invention;
FIG. 2 is a flow chart of a three-dimensional blood vessel image determination method provided by an embodiment of the invention;
FIG. 3 is a schematic diagram of a set image processing sequence according to an embodiment of the present invention;
FIG. 4 is a flowchart of yet another method for processing a blood vessel image provided by an embodiment of the present invention;
fig. 5 is a block diagram of a blood vessel image processing apparatus according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a CT system according to an embodiment of the present invention;
fig. 7 is a schematic structural diagram of a computer device in a CT system according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described by means of implementation examples with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Examples
Fig. 1 is a flowchart of a blood vessel image processing method according to an embodiment of the present invention. The technical scheme of the embodiment is suitable for improving the time resolution of CT blood vessel images. The method can be executed by the vascular image processing device provided by the embodiment of the invention, and the device can be realized in a software and/or hardware mode and is configured to be applied in a processor of electronic computer equipment. The method specifically comprises the following steps:
s101, acquiring a three-dimensional blood vessel image of a scanning part, subtracting projection data used for creating the three-dimensional blood vessel image, and a projection view angle sequence corresponding to the subtracting projection data.
Wherein, as shown in fig. 2, the three-dimensional blood vessel image is reconstructed by subtracting projection data and comprises at least two body layers. The subtracted projection data is the difference between the contrast enhanced projection data and the swept projection data. The flat scan projection data is CT projection data of a scanning part acquired before the patient is injected with the reinforcing agent, and the contrast enhancement projection data is CT projection data acquired after the patient is injected with the reinforcing agent. Since the enhancer is used to enhance angiography, the contrast enhanced projection data is subtracted from the swept projection data, and the projection data of non-vascular tissue in the contrast enhanced projection data can be removed, resulting in subtracted projection data comprising only blood vessels.
The flat scan projection data, the contrast enhancement projection data and the subtraction projection data all carry projection view angle sequences. Since the data amount of the subtracted projection data is smaller than the swept projection data and the subtracted projection data, the present embodiment preferably acquires the projection view angle sequence from the subtracted projection data to reduce the data calculation amount and increase the blood vessel image processing speed. The sequence of projection view angles includes a projection time instant and one or more projection view angles corresponding to each projection time instant.
It will be appreciated that for a single-source CT system, each projection instant corresponds to a projection view angle; for a multi-ray source CT system, each projection instant corresponds to at least two projection views.
S102, determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer according to the projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer.
The three-dimensional blood vessel image comprises at least one segment, and the set image processing sequence comprises a segment processing sequence and a body layer processing sequence in each segment, wherein the number of blood vessels of a starting body layer in each segment is smaller than the number of blood vessels of an ending body layer in each segment. Specifically, the three-dimensional blood vessel image is a brain three-dimensional blood vessel image, and the brain three-dimensional blood vessel image is divided into an upper section and a lower section by the center of the cone beam; setting the image processing sequence includes going from the base of the lower segment to the center of the cone beam, and then from the top of the upper segment to the center of the cone beam. Illustratively, as shown in fig. 3, the cone-beam center divides the three-dimensional brain blood vessel image into an upper segment and a lower segment, sequentially from a volume layer identified as 1 in the lower segment to a volume layer identified as S1, and then from a volume layer identified as s1+1 in the upper segment to a volume layer of s1+s2. Wherein S1 and S2 are natural numbers, and can be the same or different.
Wherein the light source is a radiation source of an imaging device, such as a radiation source of a CT imaging device.
The radiation emitted by the light source (radiation source) is referred to as light in this embodiment, and it can be understood that the light starts from the light source and ends in a certain detector element of the detector.
In one embodiment, the determination of the target virtual ray may be selected from: and connecting the light source with each blood vessel pixel point of the current body layer to obtain target virtual light corresponding to each blood vessel pixel point. It can be understood that the target virtual light reaches a certain probe element of the detector after being prolonged, and the projection data detected by the probe element is the projection data corresponding to the target virtual light. The target virtual light determination mode is simple, accurate and quick.
In one embodiment, the determination of the target virtual ray may be selected from: determining virtual light rays between the light source and each probe element of the detector; and taking the virtual light passing through the vascular pixel point of the current body layer as a target virtual light. Specifically, this embodiment takes as the target virtual ray the virtual ray that passes through one or more vessel pixels of the current precursor layer, and takes as the obsolete virtual ray the virtual ray that does not pass through any vessel pixel of the current precursor layer.
It can be appreciated that for a target virtual ray perpendicular to the current precursor layer, there is only one pixel point at its intersection with the current precursor layer; for a target virtual ray oblique to the current layer, the intersection of the target virtual ray and the current layer may have only one pixel point, or may have a plurality of pixel points; for a target virtual ray parallel to the current precursor layer, there are multiple pixels at its intersection with the current precursor layer.
S103, determining intersections of the target virtual rays and the current layer, and determining pixel values of all blood vessel pixel points corresponding to the current projection view angle in the current layer according to the types of the intersections and the corresponding relation of the types and pixel determination rules, wherein the number of uncorrected blood vessel pixel points included in the intersections of different types is different.
The uncorrected vascular pixel points are vascular pixel points with unknown pixel values in the current body layer at the current projection view angle. For the target virtual ray, its intersection with the current precursor layer includes at least one uncorrected vessel pixel.
The number of uncorrected vessel pixels in the one or more pixels included in each intersection is determined. The type of each intersection is determined based on the number of uncorrected vessel pixel points that each intersection includes. In this embodiment, a pixel determination rule corresponding to each type is created in advance, and then, according to the type of intersection and the correspondence between the type and the pixel determination rule, the pixel value of the vascular pixel point corresponding to the current projection view angle in the current body layer is determined.
And S140, determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the corresponding blood vessel pixel points in each body layer of each projection view angle.
Because each projection view angle corresponds to one projection moment, after the pixel values of the corresponding vascular pixel points of each projection view angle in all body layers are determined, a four-dimensional vascular image sequence can be determined according to the three-dimensional vascular image and the pixel values of the corresponding vascular pixel points of each projection view angle in each body layer.
When the four-dimensional blood vessel image sequence is displayed, the pixel value of each blood vessel pixel point of each body layer changes along with the switching of the projection view angles, for example, at the time T1 corresponding to the first projection view angle, the pixel value of each blood vessel pixel point corresponding to the first projection view angle in each body layer is displayed, at the time T2 corresponding to the second projection view angle, the pixel values of each blood vessel pixel point corresponding to the first projection view angle and the second projection view angle in each body layer are displayed, at the time TN corresponding to the Nth projection view angle, the pixel values of each blood vessel pixel point corresponding to the first projection view angle to the Nth projection view angle in each body layer are displayed, and the technical effect of improving the time resolution of the four-dimensional blood vessel image sequence to the sampling resolution of the projection view angle is achieved. Of course, the number of projection view angles displayed each time may be set according to the requirement, for example, one or M projection view angles are displayed each time, and when the projection time corresponding to the projection view angle is reached, only the pixel value of each vascular pixel point corresponding to the projection view angle in each body layer, or the pixel values of each vascular pixel point corresponding to the projection view angle and M-1 projection view angles before the projection view angle in each body layer are displayed, where M is greater than or equal to 1.
The serial data processing is performed on each volume layer in the three-dimensional blood vessel image according to the set image processing sequence, and the pixel value of each blood vessel pixel point corresponding to each view angle in the current volume layer is determined by adopting a parallel processing mode with respect to the current volume layer.
Compared with the prior art, the technical scheme of the blood vessel image processing method provided by the embodiment determines the pixel value of each blood vessel pixel point corresponding to each projection view angle in the current body layer through the type of the intersection of the target virtual light and the current body layer and the corresponding relation between the type and the pixel of the rule, and so on, and determines the pixel value of each blood vessel pixel point corresponding to each projection view angle in each body layer, thereby determining a four-dimensional blood vessel image, having low data calculation complexity and higher image processing speed and meeting the clinical requirement on image processing time; when the four-dimensional blood vessel image is displayed, the pixel value of each blood vessel pixel point in each body layer can be changed along with the switching of the projection view angle, so that the technical aim of improving the time resolution of the four-dimensional blood vessel image to the time resolution of the projection view angle is fulfilled, and the four-dimensional blood vessel image with the improved time resolution can provide more detailed blood vessel information for clinical diagnosis.
Fig. 4 is a flowchart of a blood vessel image processing method according to still another embodiment of the present invention. The embodiment of the invention is used for further optimizing the pixel value of each vascular pixel point corresponding to the current projection view angle in the current body layer according to the type of each intersection and the corresponding relation between the type and the pixel determining rule in the previous embodiment, and specifically comprises the following steps:
s210, acquiring a three-dimensional blood vessel image of a scanning part, subtracting projection data used for creating the three-dimensional blood vessel image, and a projection view angle sequence corresponding to the subtracting projection data.
S220, determining a current body layer in the three-dimensional blood vessel image based on the set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer according to the projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer.
S2301, if the intersection of all the target virtual rays and the current layer includes an uncorrected blood vessel pixel, determining a pixel value of the uncorrected blood vessel pixel based on the first pixel rule.
The intersection comprising an uncorrected blood vessel pixel point specifically comprises two cases, specifically: 1) Only one uncorrected vascular pixel is included; 2) Comprises an uncorrected blood vessel pixel point and at least one corrected blood vessel pixel point. The corrected vascular pixel point is a vascular pixel point with the determined pixel value corresponding to the current projection view angle in the current body layer.
Accordingly, this step comprises the sub-steps of:
a1, if the intersection of the target virtual light and the current body layer only comprises an uncorrected blood vessel pixel point, acquiring projection data corresponding to the target virtual light from the subtracted projection data, and taking the projection data as the pixel value of the uncorrected blood vessel pixel point.
If the intersection of the target virtual ray and the current layer only comprises an uncorrected blood vessel pixel point, projection data corresponding to the target virtual ray is generated by the uncorrected blood vessel pixel point, so that the projection data corresponding to the target virtual ray is obtained from the subtracted projection data, and the projection data is used as the pixel value of the uncorrected blood vessel pixel point. The method for determining the projection data corresponding to the target virtual light in the subtracted projection data comprises the following steps: determining a corresponding probe element identifier of the target virtual light in the detector; determining partial projection data corresponding to the current projection view angle identification in the current projection data; and taking the projection data corresponding to the probe element identification in the part of projection data as projection data of the target virtual light beam for drinking in the deduction projection data.
a2, if the intersection of the target virtual light and the current layer comprises an uncorrected blood vessel pixel point and at least one corrected blood vessel pixel point, acquiring projection data corresponding to the target virtual light from subtracted projection data, and taking the difference value of the projection data and the pixel value of the at least one corrected blood vessel pixel point as the pixel value of the uncorrected blood vessel pixel point.
If the intersection of the target virtual ray and the current layer comprises an uncorrected blood vessel pixel point and at least one corrected blood vessel pixel point, the projection data corresponding to the target virtual ray is obtained from the uncorrected blood vessel pixel point and the at least one corrected blood vessel pixel point, so that the projection data corresponding to the target virtual ray is obtained from the subtracted projection data, and then the pixel value of the at least one corrected blood vessel pixel point is subtracted from the projection data to obtain a corresponding difference value, and the difference value is taken as the pixel value of the uncorrected blood vessel pixel point.
It can be understood that when the pixel values of uncorrected blood vessel pixels in the intersection of all the target virtual rays and the current body layer are determined, the pixel values of each blood vessel pixel corresponding to the current projection view angle in the current body layer are determined.
S2302, if the intersection of any target virtual ray and the current layer includes at least two uncorrected blood vessel pixel points, determining the pixel value of each blood vessel pixel corresponding to the current projection view angle in the current layer based on the second pixel rule.
It will be appreciated that if the intersection of any of the target virtual lines with the current precursor layer includes at least two uncorrected vessel pixel points, then the projection data corresponding to the target virtual ray in the subtracted projection data needs to be assigned to the at least two uncorrected vessel pixel points. Since the projection data distribution ratio between the at least two uncorrected blood vessel pixel points cannot be known, the pixel value of each uncorrected blood vessel pixel point cannot be directly determined. To this end, in one embodiment, the pixel value of each vessel pixel point corresponding to the current projection view angle in the current body layer is determined by the following steps:
b1, if the intersection of any target virtual ray and the current precursor layer comprises at least two uncorrected blood vessel pixel points, marking the projection view angle corresponding to the target virtual ray as a pending state.
And b2, determining the pixel value of each blood vessel pixel point corresponding to the projection view angle in the pending state in the current body layer according to the pixel value of each blood vessel pixel point corresponding to the projection view angle in the first setting neighborhood of the projection view angle in the pending state.
The first setting neighborhood may be a projection view angle adjacent to the left side and the right side of the projection view angle in the undetermined state or two projection view angles nearest to the left side and the right side of the projection view angle in the undetermined state. In one embodiment, two undetermined state projection view angles closest to the undetermined state projection view angle are determined, and pixel values of each vascular pixel point corresponding to the undetermined state projection view angle in the current body layer are determined according to pixel values of each vascular pixel point corresponding to the two undetermined state projection view angles in the current body layer.
The projection view angle in the undetermined state refers to the projection view angle in which the pixel value of each vascular pixel point in the current body layer is already determined.
For example, for projection view sequence t= [ A1, A2, A3 … AM ], where M is a natural number. The states of all projection view angles in the projection view angle sequence are determined, wherein the projection view angle A3 is set as a projection view angle in a pending state, other projection view angles are projection view angles in a non-pending state, and two projection view angles in a non-pending state, which correspond to the projection view angle A3 and have the nearest distance, are respectively a projection view angle A2 and a projection view angle A4, so that the pixel value of each vascular pixel point corresponding to the projection view angle A3 in the current layer is obtained by linear interpolation according to the pixel values of each vascular pixel point corresponding to the projection view angle A2 and the projection view angle A4 in the current layer.
For example, for projection view sequence t= [ A1, A2, A3 … AM ], where M is a natural number. Setting a projection view angle A2 and a projection view angle A3 as projection view angles in a undetermined state, wherein other projection view angles are projection view angles in undetermined states, the projection view angles in the undetermined states, which correspond to the projection view angle A2, are respectively the projection view angle A1 and the projection view angle A4, so that the pixel value of each vascular pixel point corresponding to the projection view angle A2 in a current body layer is obtained through linear interpolation according to the pixel value of each vascular pixel point corresponding to the projection view angle A1 and the projection view angle A4 in the current body layer, and then the pixel value of each vascular pixel point corresponding to the projection view angle A3 in the current body layer is obtained through linear interpolation according to the pixel value of each vascular pixel point corresponding to the projection view angle A1 and the projection view angle A4 in the current body layer.
It should be noted that, the technical scheme of determining the pixel value of each blood vessel pixel point corresponding to the projection view angle in the pending state in the current body layer based on linear interpolation can be used for limb blood vessel image processing.
In one embodiment, the pixel value of each vascular pixel point corresponding to the current projection view angle in the current body layer is determined by the following steps:
And c1, if the intersection of any target virtual ray and the current precursor layer comprises at least two uncorrected blood vessel pixel points, marking the projection view angle corresponding to the target virtual ray as a pending state.
And c2, solving undetermined parameters in the initial gamma fitting model according to the pixel values corresponding to the undetermined state projection view angles in the current body layer in a set number to obtain a gamma fitting model, and determining the pixel values of all blood vessel pixel points corresponding to the undetermined state projection view angles in the current body layer according to the gamma fitting model.
Among them, the Gamma-fitted model (Gamma-variable model) is also called a time intensity quantitative analysis model. The initial gamma fitting model formula in this embodiment is as follows:
wherein,
wherein f (t) is a gamma function, t 0 Is the delay time. A. t is t 0 Alpha, beta and gamma are all undetermined parameters. The integration operation in the model may be approximately solved based on the "trapezoidal rule".
Because the undetermined parameters of the initial gamma fitting model are determined based on the corresponding pixel values of the undetermined state projection view angles in the current body layer, and each projection view angle corresponds to one projection moment, the pixel values of the corresponding vascular pixel points of the undetermined state projection view angles in the current body layer can be quantified based on the gamma fitting model. Because the undetermined parameters need to be determined through the pixel values of all the vascular pixel points corresponding to the current body layer through at least 5 undetermined state head projection visual angles, the set number in the step is more than or equal to 5.
c3, if the number of the projection view angles in the undetermined state in the current precursor layer is smaller than the set number, determining the same type of blood vessel pixel points in a second set neighborhood of each blood vessel pixel point in the current precursor layer, and obtaining the pixel value of each corresponding blood vessel pixel point of each undetermined state projection view angle in the current precursor layer by linear interpolation according to the pixel value of the same type of blood vessel pixel point in the second set neighborhood under each undetermined state projection view angle; the same type of blood vessel pixel points in the second setting neighborhood are positioned in the setting adjacent body layer of the current body layer, and each blood vessel pixel point and the same type of blood vessel pixel point in the corresponding second setting neighborhood belong to the same blood vessel.
The second set neighborhood is the intersection of the sphere neighborhood with the set radius corresponding to the vascular pixel point to be solved and the set adjacent body layer, or the intersection of the cuboid neighborhood with the set side length and the set adjacent body layer. The adjacent body layer is preferably a front-back adjacent body layer, and may be a front adjacent body layer or a back adjacent body layer, wherein the direction from the head to the foot of the three-dimensional blood vessel image is a front-back direction. For example, the precursor layer is a bulk layer N, N is greater than 1 and less than N M Wherein N is M For the maximum number of layers, it is set that the adjacent bulk layers are preferably bulk layer N-1 and bulk layer N+1, but also only bulk layer N-1 or only bulk layer N+1.
In one embodiment, the same type of vessel pixel points refer to vessel pixel points of the same vessel, and for example, the vessel pixel points to be solved in the current layer correspond to vessel walls, and the vessel pixel points of the same type in the corresponding second setting neighborhood are vessel pixel points which are located in the second setting neighborhood range of the vessel pixel points to be solved and belong to the same vessel as the vessel pixel points to be solved in the setting adjacent body layer of the current layer.
In one embodiment, the same type of vessel pixel points refer to vessel pixel points of the same portion in the same vessel. The first set neighborhood of the current layer corresponds to a first set neighborhood of the to-be-obtained vascular pixel, and the second set neighborhood of the to-be-obtained vascular pixel corresponds to a second set neighborhood of the to-be-obtained vascular pixel.
The method can be applied to a cone beam central layer or a body layer close to the cone beam central layer, and solves the problem that the pixel values of all blood vessel pixel points corresponding to all the projection view angles in the undetermined state cannot be determined by adopting a gamma fitting model because the projection view angles in the undetermined state are less.
It should be noted that, the technical scheme of determining the pixel value of each blood vessel pixel point corresponding to the projection view angle in the undetermined state in the current body layer based on the gamma fitting model can be used for various blood vessel image processing, such as cerebral blood vessel image processing.
S240, determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel value of each blood vessel pixel point corresponding to each projection view angle in each body layer.
According to the type of intersection of the target virtual light and the current layer, the pixel value of each pixel point corresponding to each projection view angle in the current layer is determined by adopting a corresponding pixel determination rule, so that the efficiency, accuracy and flexibility of determining the pixel value of each pixel point corresponding to each projection view angle in the current layer are improved, and the time resolution of the four-dimensional vascular image is improved to the time resolution of the projection view angle.
Fig. 5 is a block diagram of a blood vessel image processing apparatus according to a third embodiment of the present invention. The apparatus is used for executing the method of the vascular image processing apparatus provided in any of the above embodiments, and the apparatus may be implemented in software or hardware. The device comprises:
an acquisition module 11, configured to acquire a three-dimensional blood vessel image of a scan site, subtracted projection data for creating the three-dimensional blood vessel image, and a projection view sequence corresponding to the subtracted projection data;
A virtual light module 12, configured to determine a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determine, for each projection view angle of the sequence of projection views, at least one target virtual light corresponding to the current body layer, where the target virtual light is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
a pixel value determining module 13, configured to determine intersections of the target virtual rays and the current layer, and determine pixel values of each vessel pixel point corresponding to the current projection view angle in the current layer according to types of the intersections and a correspondence between types and pixel determining rules, where the number of uncorrected vessel pixel points included in the intersections of different types is different;
the four-dimensional image determining module 14 is configured to determine a four-dimensional vascular image sequence according to the three-dimensional vascular image and pixel values of vascular pixels of each body layer under each projection view angle.
Optionally, the three-dimensional blood vessel image comprises at least one segment; the set image processing sequence comprises a segmentation processing sequence and a body layer processing sequence in each segment, wherein the number of blood vessels of a starting body layer in each segment is smaller than the number of blood vessels of an ending body layer in each segment.
Optionally, the three-dimensional blood vessel image is a brain three-dimensional blood vessel image, and the brain three-dimensional blood vessel image is divided into an upper section and a lower section by the center of the cone beam; the set image processing sequence includes a lower segment from the skull base to the cone beam center and an upper segment from the skull top to the cone beam center.
Optionally, the number of projection views at each projection instant is greater than or equal to 1.
Optionally, the virtual light module 12 is configured to connect the light source with each blood vessel pixel of the current body layer, so as to obtain a target virtual light corresponding to each blood vessel pixel; or determining virtual light rays between the light source and each probe element of the detector, and taking the virtual light rays passing through the vascular pixel points of the current body layer as target virtual light rays.
Optionally, the pixel value determining module 13 includes:
a first unit configured to determine, if all intersections of the target virtual rays with the current layer include an uncorrected blood vessel pixel, a pixel value of the uncorrected blood vessel pixel based on a first pixel rule;
and the second unit is used for determining the pixel value of each vascular pixel corresponding to the current projection view angle in the current body layer based on a second pixel rule if the intersection of any one of the target virtual rays and the current body layer comprises at least two uncorrected vascular pixel points.
Optionally, the first unit is configured to obtain, if the intersection of the target virtual ray and the current layer includes only one uncorrected vascular pixel, projection data corresponding to the target virtual ray from the subtracted projection data, and take the projection data as a pixel value of the uncorrected vascular pixel; if the intersection of the target virtual ray and the current layer comprises an uncorrected blood vessel pixel point and at least one corrected blood vessel pixel point, acquiring projection data corresponding to the target virtual ray from the subtracted projection data, and taking the difference value of the projection data and the pixel value of the at least one corrected blood vessel pixel point as the pixel value of the uncorrected blood vessel pixel point.
Optionally, the second unit is configured to mark the projection view angle corresponding to the target virtual ray as a pending state if the intersection of any one of the target virtual rays and the current body layer includes at least two uncorrected vascular pixels; and setting pixel values of the corresponding vascular pixel points of the projection view angles in the neighborhood in the current body layer according to the projection view angles in the undetermined state, and determining the pixel values of the corresponding vascular pixel points of the projection view angles in the undetermined state in the current body layer.
Optionally, the second unit is specifically configured to mark the projection view angle corresponding to the target virtual ray as a pending state if the intersection of any one of the target virtual rays and the current body layer includes at least two uncorrected vascular pixel points; according to the pixel values corresponding to the undetermined state projection view angles in the current body layer, undetermined parameters in an initial gamma fitting model are solved to obtain a gamma fitting model, and according to the gamma fitting model, the pixel values of all blood vessel pixel points corresponding to the undetermined state projection view angles in the current body layer are determined.
Optionally, the second unit is specifically configured to determine the same type of vascular pixel point in the second set neighborhood of each vascular pixel point of the current layer if the number of projection views in the undetermined state in the current layer is less than the set number; according to the pixel values of the same type of blood vessel pixel points in the second setting neighborhood under each pending state projection view angle, linear interpolation is carried out to obtain the pixel values of the corresponding blood vessel pixel points of each pending state projection view angle in the current body layer, wherein the same type of blood vessel pixel points in the second setting neighborhood are positioned in the setting adjacent body layer of the current body layer, and the blood vessel pixel points belong to the same blood vessel with the same type of blood vessel pixel points in the corresponding second setting neighborhood.
Compared with the prior art, the technical scheme of the blood vessel image processing device provided by the embodiment determines the pixel value of each blood vessel pixel point corresponding to each projection view angle in the current body layer through the type of the intersection of the target virtual light and the current body layer and the corresponding relation between the type and the pixel of the rule, and so on, and determines the pixel value of each blood vessel pixel point corresponding to each projection view angle in each body layer, thereby determining a four-dimensional blood vessel image, having low data calculation complexity and higher image processing speed and meeting the clinical requirement on image processing time; when the four-dimensional blood vessel image is displayed, the pixel value of each blood vessel pixel point in each body layer can be changed along with the switching of the projection view angle, so that the technical aim of improving the time resolution of the four-dimensional blood vessel image to the time resolution of the projection view angle is fulfilled, and the four-dimensional blood vessel image with the improved time resolution can provide more detailed blood vessel information for clinical diagnosis.
The blood vessel image processing device provided by the embodiment of the invention can execute the blood vessel image processing method provided by any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the execution method.
Fig. 6 is a schematic structural diagram of a CT system according to another embodiment of the present invention. The system is for acquiring projection data for creating a three-dimensional vessel image of a scan site and subtracted projection data of the three-dimensional vessel image. The system includes a gantry 1211, a detector 1212, a couch 1214, an X-ray tube 1215, a C-arm drive shaft 1216, a spindle 1217, and a base 1219. The X-ray tube 1215 and the detector 1212 are provided at both ends of the C-shaped gantry 1211 with their center line perpendicular to the center axis 1218 of rotation. The C-shaped gantry 1211 rotates about a central axis of rotation 1218 to capture image data of a patient 1213 on a couch at different projection views. The X-ray tube 1215 is controlled by the X-ray generator 123 for current, voltage, exposure time, etc., and projection data acquired by the detector 1212 is transmitted by the communication system 126 to a computer device, and the gantry 1211 is coupled to the C-arm drive shaft 1216, the power of which is provided by the shaft 1217. The base 1219 is responsible for load bearing. The C-arm control unit 121 controls the rotational speed, angle, position, etc. of the gantry 1211. The spindle control unit 122 is connected to the base 1219 and provides power support for the entire C-arm system. The X-ray generator 123 controls the current, voltage and exposure time of the X-ray tube 1215. The data acquisition system 124 coordinates the gantry 1211, the detector 1212, and the X-ray generator 1215 and collects the acquired data. The couch control system 125 controls the position and movement speed of the couch 1214 to achieve different scan trajectories for the patient 1213. The communication system 126 connects the C-arm control unit 121, the spindle control unit 122, the X-ray generator 124, the data acquisition system 124, and the couch board control system 125, and transmits the acquired projection data to the memory of the computer apparatus 2.
Fig. 7 is a schematic structural diagram of a computer device according to another embodiment of the present invention, and as shown in fig. 7, the computer device 2 includes a processor 201, a memory 202, an input device 203, and an output device 204; the number of processors 201 in the device may be one or more, one processor 201 being taken as an example in fig. 7; the processor 201, memory 202, input devices 203, and output devices 204 in the apparatus may be connected by a bus or other means, for example in fig. 7.
The memory 202 serves as a computer-readable storage medium for storing a software program, a computer-executable program, and modules such as program instructions/modules (e.g., the acquisition module 11 and the virtual light module 12, the pixel value determination module 13, and the four-dimensional image determination module 14) corresponding to the blood vessel image processing method in the embodiment of the present invention. The processor 201 executes various functional applications of the apparatus and data processing, that is, implements the above-described blood vessel image processing method, by running software programs, instructions, and modules stored in the memory 202.
The memory 202 mainly includes a storage program area and a storage data area, wherein the storage program area stores an operating system, at least one application program required for functions; the storage data area stores data created according to the use of the terminal, and the like. In addition, memory 202 may include high-speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid-state storage device. In some examples, memory 202 may further include memory located remotely from processor 201, which may be connected to the device via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The input means 203 may be used to receive entered numeric or character information and to generate key signal inputs related to user settings and function control of the device. The input device may be configured at an operator workstation through which the operator controls the operation of the C-arm CT system.
The output device 204 may include a display device such as a display screen, for example, a display screen of an operating workstation.
Still another embodiment of the present invention provides a storage medium containing computer-executable instructions, which when executed by a computer processor, are for performing a vascular image processing method, the method comprising:
acquiring a three-dimensional blood vessel image of a scanning part, deduction projection data for creating the three-dimensional blood vessel image, and a projection view angle sequence corresponding to the deduction projection data;
determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer according to a projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
Determining intersection sets of the target virtual rays and the current body layer, and determining pixel values of all blood vessel pixel points corresponding to the current projection view angle in the current body layer according to the type of each intersection set and the corresponding relation of the type and the pixel determining rule, wherein the number of uncorrected blood vessel pixel points included in the intersection sets of different types is different;
and determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the corresponding blood vessel pixel points of each projection view angle in each body layer.
Of course, the storage medium containing the computer executable instructions provided in the embodiments of the present invention is not limited to the method operations described above, and may also perform the related operations in the blood vessel image processing method provided in any embodiment of the present invention.
From the above description of embodiments, it will be clear to a person skilled in the art that the present invention may be implemented by means of software and necessary general purpose hardware, but of course also by means of hardware, although in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a computer readable storage medium, such as a floppy disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a FLASH Memory (FLASH), a hard disk or an optical disk of a computer, and the like, and include several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the vascular image processing method according to the embodiments of the present invention.
It should be noted that, in the above-described embodiment of the blood vessel image processing apparatus, each unit and module included is divided according to the functional logic only, but is not limited to the above-described division, as long as the corresponding functions can be realized; in addition, the specific names of the functional units are also only for distinguishing from each other, and are not used to limit the protection scope of the present invention.
Note that the above is only a preferred embodiment of the present invention and the technical principle applied. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments described herein, but is capable of various obvious changes, rearrangements and substitutions as will now become apparent to those skilled in the art without departing from the scope of the invention. Therefore, while the invention has been described in connection with the above embodiments, the invention is not limited to the embodiments, but may be embodied in many other equivalent forms without departing from the spirit or scope of the invention, which is set forth in the following claims.

Claims (13)

1. A blood vessel image processing method, comprising:
acquiring a three-dimensional blood vessel image of a scanning part, deduction projection data for creating the three-dimensional blood vessel image, and a projection view angle sequence corresponding to the deduction projection data;
Determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual ray corresponding to the current body layer according to a projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual ray is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
determining intersection sets of the target virtual rays and the current body layer, and determining pixel values of all blood vessel pixel points corresponding to the current projection view angle in the current body layer according to the type of each intersection set and the corresponding relation of the type and the pixel determining rule, wherein the number of uncorrected blood vessel pixel points included in the intersection sets of different types is different;
and determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the corresponding blood vessel pixel points of each projection view angle in each body layer.
2. The method of claim 1, wherein the step of determining the position of the substrate comprises,
the three-dimensional vessel image comprises at least one segment;
the set image processing sequence comprises a segmentation processing sequence and a body layer processing sequence in each segment, wherein the number of blood vessels of a starting body layer in each segment is smaller than the number of blood vessels of an ending body layer in each segment.
3. The method of claim 1, wherein the step of determining the position of the substrate comprises,
the three-dimensional blood vessel image is a brain three-dimensional blood vessel image, and the brain three-dimensional blood vessel image is divided into an upper section and a lower section by the center of a cone beam;
the set image processing sequence includes from the base of the lower segment to the center of the cone beam and from the top of the upper segment to the center of the cone beam.
4. The method of claim 1, wherein the number of projection views per projection instant is greater than or equal to 1.
5. The method of claim 1, wherein determining the target virtual ray comprises:
connecting a light source with each blood vessel pixel point of the current precursor layer to obtain a target virtual ray corresponding to each blood vessel pixel point; or alternatively
Virtual light rays between the light source and each probe element of the detector are determined, and the virtual light rays passing through the vascular pixel points of the current body layer are taken as target virtual light rays.
6. The method according to claim 1, wherein determining the pixel value of each vessel pixel point corresponding to the current projection view angle in the current body layer according to the type of each intersection and the correspondence between the type and the pixel determination rule includes:
If all the intersections of the target virtual rays and the current layer comprise an uncorrected blood vessel pixel point, determining a pixel value of the uncorrected blood vessel pixel point based on a first pixel rule;
if the intersection of any one of the target virtual rays and the current body layer comprises at least two uncorrected blood vessel pixel points, determining the pixel value of each blood vessel pixel corresponding to the current projection view angle in the current body layer based on a second pixel rule.
7. The method of claim 6, wherein if all of the intersections of the target virtual rays with the current layer include an uncorrected vessel pixel, determining a pixel value for the uncorrected vessel pixel based on a first pixel rule comprises:
if the intersection of the target virtual light and the current layer only comprises an uncorrected blood vessel pixel point, acquiring projection data corresponding to the target virtual light from the subtracted projection data, and taking the projection data as a pixel value of the uncorrected blood vessel pixel point;
if the intersection of the target virtual ray and the current layer comprises an uncorrected blood vessel pixel point and at least one corrected blood vessel pixel point, acquiring projection data corresponding to the target virtual ray from the subtracted projection data, and taking the difference value of the projection data and the pixel value of the at least one corrected blood vessel pixel point as the pixel value of the uncorrected blood vessel pixel point.
8. The method of claim 6, wherein if the intersection of any of the target virtual rays with the current body layer includes at least two uncorrected vessel pixel points, determining, based on a second pixel rule, a pixel value for each vessel pixel corresponding to the current projection view angle in the current body layer, comprises:
if the intersection of any one of the target virtual light rays and the current body layer comprises at least two uncorrected blood vessel pixel points, marking the corresponding projection view angle of the target virtual light rays as a pending state;
and determining the pixel value of each blood vessel pixel point corresponding to the projection view angle in the current body layer according to the pixel value of each blood vessel pixel point corresponding to the projection view angle in the first setting neighborhood of the projection view angle in the undetermined state.
9. The method of claim 6, wherein if the intersection of any of the target virtual rays with the current body layer includes at least two uncorrected vessel pixel points, determining, based on a second pixel rule, a pixel value for each vessel pixel corresponding to the current projection view angle in the current body layer, comprises:
if the intersection of any one of the target virtual light rays and the current body layer comprises at least two uncorrected blood vessel pixel points, marking the corresponding projection view angle of the target virtual light rays as a pending state;
According to the pixel values corresponding to the undetermined state projection view angles in the current body layer, undetermined parameters in an initial gamma fitting model are solved to obtain a gamma fitting model, and according to the gamma fitting model, the pixel values of all blood vessel pixel points corresponding to the undetermined state projection view angles in the current body layer are determined.
10. The method of claim 9, wherein if the intersection of any of the target virtual rays with the current body layer includes at least two uncorrected vessel pixel points, then marking the corresponding projection view angle of the target virtual ray as pending further comprises:
if the number of projection view angles in the undetermined state in the current precursor layer is smaller than the set number, determining the same type of blood vessel pixel points in a second set neighborhood of each blood vessel pixel point in the current precursor layer, wherein the same type of blood vessel pixel points in the second set neighborhood are positioned in a set adjacent body layer of the current precursor layer, and each blood vessel pixel point and the same type of blood vessel pixel points in the corresponding second set neighborhood belong to the same blood vessel;
and linearly interpolating according to the pixel values of the same type of blood vessel pixel points in the second set neighborhood under each pending state projection view angle to obtain the pixel values of the corresponding blood vessel pixel points of each pending state projection view angle in the current body layer.
11. A blood vessel image processing apparatus, comprising:
the acquisition module is used for acquiring a three-dimensional blood vessel image of a scanning part, deduction projection data used for creating the three-dimensional blood vessel image and a projection view angle sequence corresponding to the deduction projection data;
the virtual light module is used for determining a current body layer in the three-dimensional blood vessel image based on a set image processing sequence, and determining at least one target virtual light corresponding to the current body layer according to the projection view angle of each projection moment in the projection view angle sequence, wherein the target virtual light is emitted by a light source and passes through at least one blood vessel pixel point of the current body layer;
the pixel value determining module is used for determining intersections of the target virtual rays and the current body layer, and determining pixel values of blood vessel pixel points corresponding to the current projection view angle in the current body layer according to the types of the intersections and the corresponding relation between the types and pixel determining rules, wherein the number of uncorrected blood vessel pixel points included in the intersections of different types is different;
and the four-dimensional image determining module is used for determining a four-dimensional blood vessel image sequence according to the three-dimensional blood vessel image and the pixel values of the blood vessel pixel points of each body layer under each projection view angle.
12. A computer device, the computer device comprising:
one or more processors;
a storage means for storing one or more programs;
the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the vascular image processing method of any of claims 1-10.
13. A storage medium containing computer executable instructions which, when executed by a computer processor, are for performing the vascular image processing method of any of claims 1-10.
CN202211158589.9A 2022-09-22 2022-09-22 Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium Pending CN117788680A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211158589.9A CN117788680A (en) 2022-09-22 2022-09-22 Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211158589.9A CN117788680A (en) 2022-09-22 2022-09-22 Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium

Publications (1)

Publication Number Publication Date
CN117788680A true CN117788680A (en) 2024-03-29

Family

ID=90378801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211158589.9A Pending CN117788680A (en) 2022-09-22 2022-09-22 Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium

Country Status (1)

Country Link
CN (1) CN117788680A (en)

Similar Documents

Publication Publication Date Title
US9754390B2 (en) Reconstruction of time-varying data
EP1869643B1 (en) Image processing device and method for blood flow imaging
US8131041B2 (en) System and method for selective blending of 2D x-ray images and 3D ultrasound images
EP1522045B1 (en) Motion artifact correction of tomographical images
US8111889B2 (en) Method and apparatus for efficient calculation and use of reconstructed pixel variance in tomography images
US20100201786A1 (en) Method and apparatus for reconstructing an image
US9968324B2 (en) Generating a 2D projection image of a vascular system
US10083511B2 (en) Angiographic roadmapping mask
JP2004320771A (en) Method for performing digital subtraction angiography
CN110473269B (en) Image reconstruction method, system, equipment and storage medium
WO2009147605A1 (en) Reconstruction of dynamical cardiac spect for measuring tracer uptake and redistribution
US20150279084A1 (en) Guided Noise Reduction with Streak Removal for High Speed C-Arm CT
US11179122B2 (en) Bolus imaging
US10043268B2 (en) Medical image processing apparatus and method to generate and display third parameters based on first and second images
Stiel et al. Digital flashing tomosynthesis: a promising technique for angiocardiographic screening
US11317875B2 (en) Reconstruction of flow data
EP3649957B1 (en) Device and method for editing a panoramic radiography image
Manhart et al. Fast dynamic reconstruction algorithm with joint bilateral filtering for perfusion C-arm CT
CN116051463A (en) Medical image processing method, medical image processing device, computer equipment and storage medium
CN117788680A (en) Vascular image processing method, vascular image processing device, vascular image processing computer equipment and storage medium
CN113313649B (en) Image reconstruction method and device
US20220414832A1 (en) X-ray imaging restoration using deep learning algorithms
CN111062998B (en) Image reconstruction method, device, CT system and storage medium
JP2847539B2 (en) Digital angiography imaging system
CN110840479A (en) Bone component determination method, device, equipment and storage medium

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