Image reconstruction method, device, storage medium and electronic device

Abstract

The present application discloses an image reconstruction method, device, storage medium, and electronic device, wherein the method includes: acquiring an original scanned image; the original scanned image is a three-dimensional stereoscopic image; performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image; the first plane is a sagittal plane or a coronal plane; extracting a cross-sectional high-frequency component from the original scanned image to obtain a target high-frequency component; the target high-frequency component includes at least a first high-frequency component corresponding to the first plane; compensating the first image based on the target high-frequency component to obtain a first target image. The present application obtains a processed image by performing a sagittal or coronal high-frequency component removal process on the original scanned image, and then compensating the target high-frequency information extracted from the original scanned image to the processed image, thereby obtaining an image after removing windmill artifacts, making the reconstructed image clearer.

Description

Technical Field

[0001] The present application relates to the field of tomography technology, and in particular to an image reconstruction method, device, storage medium and electronic device. Background Art

[0002] With the rapid development of CT technology, multi-slice CT is becoming increasingly common. Windmill artifact is a common artifact in multi-slice CT. This artifact appears as alternating white and black bars, with the number of white and black bars equal to the number of detector slices. Currently, there are two main methods for removing windmill artifacts. One method primarily uses hardware to increase the axial sampling rate to remove windmill artifacts, such as the Z-flying focus technique. This method has the disadvantage of requiring expensive hardware support equipment. The other method primarily uses software to remove windmill artifacts through algorithms. The current main algorithm removes windmill artifacts through image post-processing, but the drawback is that the removal effect is not very ideal.

[0003] Therefore, an image reconstruction method is urgently needed to obtain an image after removing windmill artifacts. Summary of the Invention

[0004] In view of this, the present invention provides an image reconstruction method, device, storage medium and electronic device, the main purpose of which is to solve the problem of inaccurate windmill artifact removal in the current image reconstruction process.

[0005] To solve the above problems, the present application provides an image reconstruction method, comprising:

[0006] Acquire an original scanned image, where the original scanned image is a three-dimensional stereo image;

[0007] performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, wherein the first plane is a sagittal plane or a coronal plane;

[0008] Extracting high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components at least include a first high-frequency component corresponding to the first surface;

[0009] The first image is compensated based on the target high-frequency component to obtain a first target image.

[0010] Optionally, the target high-frequency component further includes a second high-frequency component corresponding to the second surface; after obtaining the reconstructed first target image, the method further includes:

[0011] performing a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane;

[0012] The second image is compensated based on the second high-frequency component in the target high-frequency component to obtain a second target image.

[0013] Optionally, the method further includes:

[0014] Determine the target pitch by doubling the sampling in the z direction;

[0015] Scanning the target object based on the target pitch to obtain scanning data;

[0016] Image reconstruction is performed based on the scan data to obtain the original scan image.

[0017] Optionally, the determining the target pitch by doubling the sampling in the z direction specifically includes:

[0018] Get the total number of detection units in the detector;

[0019] Determine the target moving layer number; the target moving layer number is the number of layers that the scanning bed moves relative to the detector during one scanning circle;

[0020] The target pitch is calculated using a predetermined formula based on the total number of detection units in the detector and the number of target moving layers;

[0021] The predetermined calculation formula is:

[0022] Where P represents the target pitch; m represents the number of target moving layers; and N represents the total number of detection units in the detector.

[0023] Optionally, the first high-frequency component removal processing method includes a Z-direction filtering method or a TV noise reduction method; the second high-frequency component removal processing method includes a Z-direction filtering method or a TV noise reduction method.

[0024] To solve the above technical problems, the present application provides an image reconstruction device, comprising:

[0025] An acquisition module, configured to acquire an original scanned image, wherein the original scanned image is a three-dimensional stereoscopic image;

[0026] a removal module, configured to perform a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, wherein the first plane is a sagittal plane or a coronal plane;

[0027] an extraction module, configured to extract high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components include at least a first high-frequency component corresponding to the first surface;

[0028] A compensation module is used to compensate the first image based on the target high-frequency component to obtain a first target image.

[0029] Optionally, the target high-frequency component further includes a second high-frequency component corresponding to the second surface;

[0030] The removal module is further configured to: perform a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane;

[0031] The compensation module is further configured to compensate the second image based on the second high-frequency component in the target high-frequency component to obtain a second target image.

[0032] Optionally, the image reconstruction device also includes a determination module, which is used to: determine the target pitch by doubling the sampling in the z direction; the acquisition module is used to: scan the target object based on the target pitch to obtain scanning data; and perform image reconstruction based on the scanning data to obtain the original scanning image.

[0033] To solve the above problems, the present application provides a storage medium, wherein the storage medium stores a computer program, and when the computer program is executed by a processor, the steps of any of the above-mentioned image reconstruction methods are implemented.

[0034] To solve the above problems, the present application provides an electronic device, which includes at least a memory and a processor, wherein a computer program is stored on the memory, and the processor implements the steps of any of the above-mentioned image reconstruction methods when executing the computer program on the memory.

[0035] The image reconstruction method, device, storage medium and electronic device in the present application obtain a processed image by processing the high-frequency components of the sagittal plane or coronal plane of the original three-dimensional scan image, and then compensate the high-frequency information of the axial image corresponding to the sagittal plane or coronal plane extracted from the original scan image to the processed image, thereby obtaining a complete image after removing the windmill artifacts, improving the accuracy of windmill artifact removal, and making the reconstructed image clearer.

[0036] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are specifically listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Various other advantages and benefits will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiment below. The accompanying drawings are for illustration purposes only and are not to be considered as limiting the present invention. The same reference symbols are used throughout the drawings to represent the same components. In the drawings:

[0038] Figure 1 This is a flowchart of an image reconstruction method according to an embodiment of the present application;

[0039] Figure 2 This is a flowchart of an image reconstruction method according to another embodiment of the present application;

[0040] Figure 3 This is a flowchart of an image reconstruction method according to another embodiment of the present application;

[0041] Figure 4 This is a schematic diagram of the sampling distribution in the z direction according to another embodiment of the present application;

[0042] Figure 5 This is a flowchart of an image reconstruction device according to another embodiment of the present application.

[0043] Figure 6 a is the image corresponding to the cross section in the original scan image without removing the windmill artifacts;

[0044] Figure 6 b is an image corresponding to the cross section of the second target image after the windmill artifacts are removed using the method of this application;

[0045] Figure 7 a is the image corresponding to the cross section in the original scan image without removing the windmill artifacts;

[0046] Figure 7 b is an image corresponding to the cross section of the second target image after the windmill artifacts are removed using the method of this application;

[0047] Figure 8 a is the image corresponding to the sagittal plane in the original scan image without removing the windmill artifact;

[0048] Figure 8 b is the image corresponding to the sagittal plane in the second target image after the windmill artifacts are removed using the method in this application;

[0049] Figure 9 a is the image corresponding to the coronal plane in the original scan image without removing the windmill artifact;

[0050] Figure 9 b is the image corresponding to the coronal plane in the second target image after the windmill artifacts are removed using the method in this application. DETAILED DESCRIPTION

[0051] Various aspects and features of the present application are described herein with reference to the accompanying drawings.

[0052] It should be understood that various modifications may be made to the embodiments of the present application. Therefore, the above description should not be considered as limiting, but merely as an example of an embodiment. Other modifications within the scope and spirit of the present application will occur to those skilled in the art.

[0053] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the present application and, together with the general description of the present application given above and the detailed description of the embodiments given below, serve to explain the principles of the present application.

[0054] These and other characteristics of the present application will become apparent from the following description of a preferred form of embodiment given as a non-limiting example with reference to the accompanying drawings.

[0055] It should also be understood that although the present application has been described with reference to certain specific examples, those skilled in the art will readily be able to implement many other equivalent forms of the present application.

[0056] The above and other aspects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings.

[0057] Specific embodiments of the present application will be described hereinafter with reference to the accompanying drawings; however, it should be understood that the embodiments described are merely examples of the present application and may be implemented in a variety of ways. Familiar and / or repetitive functions and structures are not described in detail to avoid obscuring the present application with unnecessary or redundant details. Therefore, the specific structural and functional details described herein are not intended to be limiting, but rather serve merely as a basis and representative basis for the claims to teach those skilled in the art to variously utilize the present application with substantially any suitable detailed structure.

[0058] This specification may use the phrases "in one embodiment," "in another embodiment," "in yet another embodiment," or "in other embodiments," which may all refer to one or more of the same or different embodiments according to the present application.

[0059] The present application embodiment provides an image reconstruction method, such as Figure 1 As shown, it includes the following steps:

[0060] Step S101, obtaining an original scanned image, wherein the original scanned image is a three-dimensional stereo image;

[0061] In the specific implementation of this step, the target object can be scanned according to predetermined scanning parameters to obtain scan data, and then the scan data can be reconstructed to obtain the original scan image. When setting the scanning parameters, the scanning parameters can be set with the goal of doubling the sampling in the Z direction (along the scanning bed).

[0062] Step S102, performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, where the first plane is a sagittal plane or a coronal plane;

[0063] After obtaining the original scanned image, high-frequency components in the sagittal plane or coronal plane can be removed to obtain a first image, which is also a three-dimensional image. Specifically, Z-direction filtering or TV noise reduction can be used to remove high-frequency components in the sagittal or coronal planes.

[0064] Step S103, extracting high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components at least include a first high-frequency component corresponding to the first shape surface;

[0065] In this step, when the high-frequency component removal operation is performed on the sagittal plane in step S102, the cross-sectional high-frequency component extraction is performed on the original scanned image to obtain the target high-frequency component, which includes the first high-frequency component corresponding to the sagittal plane; when the high-frequency component removal operation is performed on the coronal plane in step S102, the cross-sectional high-frequency component extraction is performed on the original scanned image to obtain the target high-frequency component, which includes the first high-frequency component corresponding to the coronal plane.

[0066] Step S104 : compensating the first image based on the target high-frequency component to obtain a first target image.

[0067] During the specific implementation of this step, the target high-frequency information can be compensated into the processed image by image fusion, thereby obtaining the target image after removing the windmill artifacts.

[0068] The image reconstruction method in the present application obtains a processed image by removing the high-frequency components of the sagittal or coronal plane of the original three-dimensional scanned image, and then compensates the high-frequency information of the axial image corresponding to the sagittal or coronal plane extracted from the original scanned image to the processed image, thereby obtaining a complete image after removing the windmill artifacts, improving the accuracy of windmill artifact removal, and making the reconstructed image clearer and more complete.

[0069] Another embodiment of the present application provides an image reconstruction method, such as Figure 2 As shown, it includes the following steps:

[0070] Step S201, obtaining an original scanned image, wherein the original scanned image is a three-dimensional image;

[0071] Step S202, performing high-frequency component removal processing on a first plane of the original scanned image to obtain a processed first image, where the first plane is a sagittal plane or a coronal plane;

[0072] Step S203, extracting high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components include first high-frequency components corresponding to the first shape surface;

[0073] During the specific implementation of this step, for example, a high-pass filtering method may be used to extract the target high-frequency component, that is, to extract the first high-frequency component corresponding to the first surface.

[0074] Step S204, compensating the first image based on the target high-frequency component to obtain a first target image;

[0075] Step S205: performing a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane;

[0076] Step S206 : Compensating the second image based on the second high-frequency component in the target high-frequency component to obtain a second target image.

[0077] During the specific implementation of this embodiment, the first high-frequency information and the second high-frequency information may be compensated into the processed image by image fusion, thereby obtaining the second target image after removing the windmill artifacts.

[0078] During the specific implementation of this embodiment, after obtaining the original scanned image, the high-frequency component removal processing of the sagittal plane and the high-frequency component removal processing of the coronal plane can also be directly performed on the original scanned image in sequence to obtain a processed image. For example, the high-frequency component removal operation of the sagittal plane is performed on the original scanned image to obtain a first image. Then, the high-frequency component removal operation of the coronal plane is directly performed on the first image to obtain a processed second image, and both the first image and the second image are three-dimensional stereo images. Subsequently, the first high-frequency component and the second high-frequency component in the target high-frequency component are used to perform compensation operations in sequence, thereby obtaining a second target image. In this step, the high-frequency components of the sagittal plane and the coronal plane can be removed by Z-direction filtering or TV noise reduction.

[0079] In this embodiment, the Z-direction filtering method is used to remove the high-frequency components of the sagittal and coronal planes as an example for specific description. The specific process of obtaining the target image is as follows:

[0080] (1) Let I(x,y,z) be the original scanned image, I 矢 (x, y, z) is the first image after removing high-frequency components from the sagittal plane, then I 矢 (x,y,z) can be obtained using the following formula:

[0081]

[0082] Where x' indicates that the dimension x is fixed, u is the image variable, ▽u refers to the image difference, and λ refers to the regularization parameter of the TV method.

[0083] (2)In I 矢 (x, y, z) based on the coronal plane to remove the high-frequency components to obtain the second image I 矢冠 (x,y,z), can be obtained using the following formula:

[0084]

[0085] (3) Extract the target high-frequency component of the axial (coronal) image from the original image I (x, y, z) to obtain the target high-frequency component H including the first high-frequency component and the second high-frequency component 轴 (x,y,z), which can be obtained specifically in the following way:

[0086] H 轴 (x,y,z)=I(x,y,z)-I 轴 (x,y,z)

[0087]

[0088] Among them, I 轴 represents the axial image; u represents the image variable; ▽u represents the difference of the image; λ represents the canonical parameter of the TV method.

[0089] (4) Compensate the axial target high-frequency component to I 矢冠 Get the target image I on (x,y,z) new (x, y, z), the following formula can be used for compensation.

[0090] I new (x,y,z)=I 矢冠 (x,y,z)+H 轴 (x,y,z)

[0091] In the present application, the characteristic of windmill artifacts in the image domain that they appear as low-frequency components in the axial (xy direction) image and as high-frequency components in the sagittal plane (yz direction) and the coronal plane (xz direction) are utilized. The high-frequency components of the sagittal and coronal planes of the original scanned image can be processed to obtain a processed image, and then the high-frequency information of the axial image extracted from the original scanned image is compensated to the processed image, thereby obtaining a complete image after removing the windmill artifacts, thereby improving the accuracy of windmill artifact removal.

[0092] On the basis of the above embodiments, in order to improve the removal effect of windmill artifacts, another embodiment of the present application provides an image reconstruction method, which can be specifically as follows: Figure 3 As shown, the following steps are included:

[0093] Step S301 : determining a target pitch by doubling sampling in the z direction; scanning a target object based on the target pitch to obtain scan data; and reconstructing an image based on the scan data to obtain an original scan image.

[0094] In this step, the z direction refers to the direction along the length of the scanning bed. When determining the pitch, the following method can be used to obtain it: obtain the total number of detection units in the detector; determine the target number of movement layers; the target number of movement layers is the number of layers that the scanning bed moves relative to the detector during one scan; based on the total number of detection units in the detector and the target number of movement layers, use a predetermined calculation formula to calculate the target pitch. The predetermined calculation formula is: Where P represents the target pitch; m represents the number of target moving layers; and N represents the total number of detection units in the detector.

[0095] In this step, the z-direction sampling is doubled to achieve the z-direction sampling of the current scan and the opposite scan in the spiral scan, which is exactly crossed / alternated (i.e., the scanning boundaries of the current scan and the opposite scan overlap in the z-direction sampling. This overlap can be expressed as actual overlap, or the distance difference is within 200 microns). This achieves the effect of doubling the z-direction sampling. θ for:

[0096]

[0097] Therefore, in order to achieve detector sampling doubling through the opposite side, the following relationship needs to be satisfied:

[0098]

[0099] Where m = 0, 1, …, (N-1) / 2; m represents the number of layers the scanning bed moves relative to the detector during one scan; θ represents the rotation angle; p represents the pitch; T represents the distance between two adjacent detection units; N represents the total number of detection units; and R represents the distance from the light source to the rotation center.

[0100] Therefore, when the pitch satisfies the following relationship, the z-direction sampling can be doubled.

[0101]

[0102] In this step, N=16, m=4, and then the target pitch is calculated. After scanning with the target pitch, the following can be obtained: Figure 4 Schematic diagram of the z-direction sampling distribution shown. Figure 3 In the figure, the dots “●”, “ο” and “*” represent the sampling in the z direction (θ-π, θ, θ+π) respectively.

[0103] Step S302, performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image; the first plane is a sagittal plane or a coronal plane;

[0104] Step S303, extracting high-frequency components of the cross section of the original scanned image to obtain target high-frequency components; the target high-frequency components at least include a first high-frequency component corresponding to the first surface;

[0105] Step S304, compensating the first image based on the target high-frequency component to obtain a first target image;

[0106] Step S305: performing a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane;

[0107] Step S306, compensating the second image based on the second high-frequency component in the target high-frequency component to obtain a second target image;

[0108] In this step, the second high-frequency component is a high-frequency component corresponding to the second shape surface.

[0109] In this application, windmill artifacts are removed by combining a data acquisition method with an image domain. The data acquisition method mainly refers to doubling the sampling in the z direction by selecting appropriate scanning parameters, thereby reducing some windmill artifacts. The image domain method mainly uses the characteristics of windmill artifacts that appear as low-frequency components in the axial (xy direction) image and as high-frequency components in the sagittal plane (yz direction) and coronal plane (xz direction) to remove them. That is, by combining the two, windmill artifacts can be removed more accurately, while maintaining the image resolution unchanged, making the effect of removing windmill artifacts more ideal, and thus making the reconstructed image clearer. For details, please refer to Figure 6 、 Figure 7 、 Figure 8 as well as Figure 9 The comparison between a and b is before and after removing the windmill artifact. Figure 6 The positions pointed by arrows / circled in a, 7a, 8a, and 9a represent windmill artifacts that have not been removed.

[0110] Another embodiment of the present application provides an image reconstruction device, such as Figure 5 As shown, including:

[0111] An acquisition module 1 is used to acquire an original scanned image, where the original scanned image is a three-dimensional image;

[0112] a removal module 2, configured to perform a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, wherein the first plane is a sagittal plane or a coronal plane;

[0113] Extraction module 3, configured to extract high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components include at least a first high-frequency component corresponding to the first surface;

[0114] The compensation module 4 is configured to compensate the first image based on the target high-frequency component to obtain a first target image.

[0115] During the specific implementation of this embodiment, the target high-frequency component also includes a second high-frequency component corresponding to the second shape plane; the second shape plane of the first target image is processed to remove the second high-frequency component to obtain a processed second image; the second shape plane is a sagittal plane or a coronal plane, and the second shape plane is different from the first shape plane; the compensation module is also used to: compensate the second image based on the second high-frequency component in the target high-frequency component to obtain a second target image.

[0116] Specifically, the image reconstruction device in this embodiment also includes a determination module, which is used to determine the target pitch by doubling the sampling in the z direction; the acquisition module is specifically used to: scan the target object based on the target pitch to obtain scanning data; and perform image reconstruction based on the scanning data to obtain the original scanning image.

[0117] Specifically, the determination module is specifically used to: obtain the total number of detection units in the detector; determine the target movement layer number; the target movement layer number is the number of layers that the scanning bed moves relative to the detector during one scan; based on the total number of detection units in the detector and the target movement layer number, calculate the target pitch using a predetermined formula. The predetermined calculation formula is: Where P represents the target pitch; m represents the number of target moving layers; and N represents the total number of detection units in the detector.

[0118] Specifically, the removal module is specifically used to: perform sagittal plane high-frequency component removal processing and coronal plane high-frequency component removal processing on the original scan image in sequence based on the Z-direction filtering method; or use the TV noise reduction method to perform sagittal plane high-frequency component removal processing and coronal plane high-frequency component removal processing on the original scan image in sequence.

[0119] In this application, the high-frequency components of the sagittal or coronal plane of the three-dimensional original scanned image are removed to obtain a processed image, and then the target high-frequency information of the axial image extracted from the original scanned image is compensated to the processed image, thereby obtaining a complete image after removing the windmill artifacts, thereby improving the accuracy of windmill artifact removal.

[0120] Another embodiment of the present application provides a storage medium storing a computer program. When the computer program is executed by a processor, the steps of the following method are implemented:

[0121] Step 1: Obtain an original scanned image; the original scanned image is a three-dimensional image;

[0122] Step 2: performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image; the first plane is a sagittal plane or a coronal plane;

[0123] Step 3: extracting high-frequency components of the cross section of the original scanned image to obtain target high-frequency components; the target high-frequency components at least include a first high-frequency component corresponding to the first surface;

[0124] Step 4: Compensate the processed first image based on the target high-frequency component to obtain a first target image.

[0125] The specific implementation process of the above method steps can be found in any of the above-mentioned embodiments of the image reconstruction method, and will not be repeated in this embodiment.

[0126] In this application, the high-frequency components of the sagittal or coronal plane of the original three-dimensional scanned image are removed to obtain a processed image, and then the target high-frequency information of the axial image extracted from the original scanned image is compensated to the processed image, thereby obtaining a complete image after removing the windmill artifacts, improving the accuracy of windmill artifact removal, and making the reconstructed image clearer and more complete.

[0127] Another embodiment of the present application provides an electronic device including at least a memory and a processor, wherein the memory stores a computer program, and the processor implements the following method steps when executing the computer program in the memory:

[0128] Step 1: Obtain an original scanned image; the original scanned image is a three-dimensional image;

[0129] Step 2: performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image; the first plane is a sagittal plane or a coronal plane;

[0130] Step 3: extracting high-frequency components of the cross section of the original scanned image to obtain target high-frequency components; the target high-frequency components at least include a first high-frequency component corresponding to the first surface;

[0131] Step 4: Compensate the processed first image based on the target high-frequency component to obtain a first target image.

[0132] The specific implementation process of the above method steps can be found in any of the above-mentioned embodiments of the image reconstruction method, and will not be repeated in this embodiment.

[0133] In this application, the high-frequency components of the sagittal or coronal plane of the original three-dimensional scanned image are removed to obtain a processed image, and then the target high-frequency information of the axial image extracted from the original scanned image is compensated to the processed image, thereby obtaining a complete image after removing the windmill artifacts, improving the accuracy of windmill artifact removal, and making the reconstructed image clearer and more complete.

[0134] The above embodiments are merely exemplary embodiments of the present application and are not intended to limit the scope of the present application. The scope of protection of the present application is defined by the claims. Those skilled in the art may make various modifications or equivalent substitutions to the present application within the essence and scope of protection of the present application, and such modifications or equivalent substitutions shall also be deemed to fall within the scope of protection of the present application.

Claims

1. An image reconstruction method, characterized in that: The method is used to remove windmill artifacts, where the windmill artifacts have the characteristics of appearing as low-frequency components in the axial image and appearing as high-frequency components in the sagittal and coronal planes in the image domain. The method includes: Acquire an original scanned image, where the original scanned image is a three-dimensional stereo image; performing a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, wherein the first plane is a sagittal plane or a coronal plane; Extracting high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components at least include a first high-frequency component corresponding to the first surface; Compensating the first image based on the target high-frequency component to obtain a first target image; The method further comprises: Determine the target pitch by doubling the sampling in the z direction; Scanning the target object based on the target pitch to obtain scanning data; Performing image reconstruction based on the scan data to obtain the original scan image; Determining the target pitch by doubling the sampling in the z direction specifically includes: Get the total number of detection units in the detector; Determine the target moving layer number; the target moving layer number is the number of layers that the scanning bed moves relative to the detector during one scanning circle; The target pitch is calculated using a predetermined formula based on the total number of detection units in the detector and the number of target moving layers; The predetermined calculation formula is: Where P represents the target pitch; m represents the number of target moving layers; and N represents the total number of detection units in the detector.

2. The method according to claim 1, wherein The target high-frequency component further includes a second high-frequency component corresponding to the second surface; after obtaining the reconstructed first target image, the method further includes: performing a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane; The second image is compensated based on the second high-frequency component in the target high-frequency component to obtain a second target image.

3. The method according to claim 2, wherein The first high-frequency component removal processing method includes a Z-direction filtering method or a TV noise reduction method; the second high-frequency component removal processing method includes a Z-direction filtering method or a TV noise reduction method.

4. An image reconstruction device, characterized in that The device is used to remove windmill artifacts, wherein the windmill artifacts have the characteristics of appearing as low-frequency components in the axial image and appearing as high-frequency components in the sagittal and coronal planes in the image domain, and comprises: An acquisition module, configured to acquire an original scanned image, wherein the original scanned image is a three-dimensional stereoscopic image; a removal module, configured to perform a first high-frequency component removal process on a first plane of the original scanned image to obtain a processed first image, wherein the first plane is a sagittal plane or a coronal plane; an extraction module, configured to extract high-frequency components of a cross section of the original scanned image to obtain target high-frequency components, wherein the target high-frequency components include at least a first high-frequency component corresponding to the first surface; a compensation module, configured to compensate the first image based on the target high-frequency component to obtain a first target image; Wherein, it also includes a determination module, the determination module is used to: determine the target pitch by doubling the sampling in the z direction; The acquisition module is configured to: scan the target object based on the target pitch to obtain scan data, and reconstruct an image based on the scan data to obtain the original scan image; The determination module is further configured to: obtain the total number of detection units in the detector; Determine the target moving layer number; the target moving layer number is the number of layers that the scanning bed moves relative to the detector during one scanning circle; The target pitch is calculated using a predetermined formula based on the total number of detection units in the detector and the number of target moving layers; The predetermined calculation formula is: Where P represents the target pitch; m represents the number of target moving layers; and N represents the total number of detection units in the detector.

5. The device according to claim 4, characterized in that The target high-frequency component also includes a second high-frequency component corresponding to the second shape surface; The removal module is further configured to: perform a second high-frequency component removal process on a second plane of the first target image to obtain a processed second image, wherein the second plane is a sagittal plane or a coronal plane, and the second plane is different from the first plane; The compensation module is further configured to compensate the second image based on the second high-frequency component in the target high-frequency component to obtain a second target image.

6. A storage medium, characterized in that The storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the image reconstruction method according to any one of claims 1 to 3 are implemented.

7. An electronic device, characterized in that: The apparatus comprises at least a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps of the image reconstruction method according to any one of claims 1 to 3 when executing the computer program in the memory.

Images