CN110070589B - CT image correction method and device, storage medium and electronic equipment - Google Patents
CT image correction method and device, storage medium and electronic equipment Download PDFInfo
- Publication number
- CN110070589B CN110070589B CN201910349470.1A CN201910349470A CN110070589B CN 110070589 B CN110070589 B CN 110070589B CN 201910349470 A CN201910349470 A CN 201910349470A CN 110070589 B CN110070589 B CN 110070589B
- Authority
- CN
- China
- Prior art keywords
- image
- edge layer
- frequency information
- information
- layer image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000003702 image correction Methods 0.000 title description 3
- 238000012937 correction Methods 0.000 claims abstract description 19
- 238000002591 computed tomography Methods 0.000 claims description 66
- 238000001914 filtration Methods 0.000 claims description 55
- 230000002146 bilateral effect Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 12
- 238000000605 extraction Methods 0.000 claims description 10
- 230000002123 temporal effect Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 description 10
- 210000001519 tissue Anatomy 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 210000000988 bone and bone Anatomy 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Abstract
The application provides a correction method, a device, a storage medium and electronic equipment of a CT image, which are used for eliminating the integral difference between an edge layer image and other layer images so as to solve the problem of integral brightness or darkness of the edge layer image, wherein the correction method of the CT image comprises the following steps: obtaining a CT image obtained by scanning a subject; extracting non-low frequency information from an edge layer image of the CT image, and obtaining an edge layer image only containing low frequency information; extracting non-artifact information from the non-low frequency information; and compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image.
Description
Technical Field
The present disclosure relates to the field of image processing technologies, and in particular, to a method and apparatus for correcting a CT image, a storage medium, and an electronic device.
Background
CT (Computed Tomography) images are currently a common medical image, CT images are image sequences, comprising edge layer images and other layer images, wherein the edge layer images are images of the edge layer during tomography, and the other layer images are images of layers other than the edge layer. In the tomographic scanning, the edge layer is influenced by factors such as upper slice shielding, so that the signals received by the edge layer have deviation, and the formed edge layer image is brighter or darker than other layer images, so that obvious streak artifacts exist in the MPR image when the MPR (Multiplanar Reformation, multi-plane reconstruction) of the CT image is carried out, and the image quality is seriously influenced.
Disclosure of Invention
In view of this, the present application provides a method, apparatus, storage medium and electronic device for correcting CT images, which are used to eliminate the overall difference between the edge layer image and other layer images, so as to solve the problem of overall brightness or darkness of the edge layer image.
In a first aspect, embodiments of the present application provide a method for correcting a CT image, the method including:
obtaining a CT image obtained by scanning a subject;
extracting non-low frequency information from an edge layer image of the CT image, and obtaining an edge layer image only containing low frequency information;
extracting non-artifact information from the non-low frequency information;
and compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image.
According to the method, the non-low frequency information is firstly extracted from the edge layer image of the CT image, the edge layer image only containing the low frequency information is obtained, then the non-artifact information is extracted from the non-low frequency information, the edge layer image only containing the low frequency information is compensated through the non-artifact information, and the corrected edge layer image is obtained, and because the corrected edge layer image does not contain the artifact information, the integral difference between the edge layer image and other layer images can be eliminated, and therefore the problem that the edge layer image is shiny or dark integrally is solved.
In one possible implementation manner, the extracting non-low frequency information from the edge layer image of the CT image and obtaining the edge layer image only including the low frequency information includes:
performing time domain filtering on the edge layer image of the CT image to obtain a first filtering image which is used as the edge layer image only containing low-frequency information;
and subtracting the edge layer image of the CT image from the first filtering image to obtain a first difference image containing non-low frequency information.
In a possible implementation manner, the extracting non-artifact information from the non-low frequency information includes:
performing spatial domain filtering on the first difference image to obtain a second filtered image only containing artifact information;
subtracting the first difference image from the second filtered image to obtain a second difference image containing non-artifact information;
the compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image comprises the following steps:
and adding the first filtered image and the second difference image to obtain a corrected edge layer image.
In a possible implementation manner, the performing temporal filtering on the edge layer image of the CT image to obtain a first filtered image includes:
and performing time domain filtering on an edge layer image of the CT image by adopting a one-dimensional bilateral filter in the CT scanning direction to obtain a first filtering image.
In the method, the bilateral filter is adopted to carry out time domain filtering on the edge layer image of the CT image, so that the edge detail information of the image can be well reserved, and the image boundary blurring is prevented.
In a possible implementation manner, the weight of the one-dimensional bilateral filter is calculated according to a first specified formula, where the first specified formula is:
wherein n is the coordinate in a one-dimensional convolution window of the CT scanning direction, z is the central coordinate of the one-dimensional convolution window, g (z) is the CT value corresponding to the z point, and sigma s1 Sum sigma r1 Is a filter parameter.
In a possible implementation manner, the performing spatial domain filtering on the first difference image to obtain a second filtered image containing non-artifact information includes:
and performing spatial domain filtering on the first difference image by adopting a two-dimensional bilateral filter to obtain a second filtered image only containing artifact information.
In the method, the bilateral filter is adopted to carry out spatial domain filtering on the first difference image, so that the edge detail information of the image can be well reserved, and the image boundary blurring is prevented.
In a possible implementation manner, the weight of the two-dimensional bilateral filter is calculated according to a second specified formula, where the second specified formula is:
where (i, j) is the seat in the two-dimensional convolution windowThe label (x, y) is the central coordinate of the two-dimensional convolution window, g (x, y) is the CT value corresponding to the (x, y) point, sigma s2 Sum sigma r2 Is a filter parameter.
In a second aspect, embodiments of the present application further provide a device for correcting a CT image, including a module for performing the method for correcting a CT image in the first aspect or any possible implementation manner of the first aspect.
In a third aspect, embodiments of the present application further provide a storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the method for correcting CT images in the first aspect or any possible implementation manner of the first aspect.
In a fourth aspect, embodiments of the present application further provide an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method for correcting CT images in the first aspect or any possible implementation manner of the first aspect when the program is executed.
Drawings
Fig. 1 is a flowchart of a method for correcting a CT image according to an embodiment of the present application;
FIG. 2 is a pre-correction MPR image provided in an embodiment of the present application;
FIG. 3 is a corrected MPR image provided by an embodiment of the present application;
fig. 4 is a schematic structural diagram of a CT image correction apparatus according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples are not representative of all implementations consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with some aspects of the present application as detailed in the accompanying claims.
The terminology used in the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application. As used in this application and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.
It should be understood that although the terms first, second, third, etc. may be used herein to describe various information, these information should not be limited by these terms. These terms are only used to distinguish one type of information from another. For example, a first message may also be referred to as a second message, and similarly, a second message may also be referred to as a first message, without departing from the scope of the present application. The word "if" as used herein may be interpreted as "at … …" or "at … …" or "responsive to a determination", depending on the context.
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
Referring to fig. 1, an embodiment of the present application provides a method for correcting a CT image, which may include the steps of:
s101, obtaining a CT image obtained by scanning a detected body;
wherein the CT image includes an edge layer image and other layer images. For example: in the CT scanning process, if the detected body needs to sequentially perform 10 section scans (or 10 circles of scans), each section scan can obtain 32 tomographic images, and the tomographic images are respectively 1-32 tomographic images along the CT scanning direction, the 1 st tomographic image and the 32 st tomographic image are edge layer images, the 2-31 st tomographic images are other layer images, and the obtained CT image is an image sequence containing 10×32 images. Of course, the 1, 2 and 31, 32 tomograms may be edge layer images, and the 3-30 tomograms may be other layer images, which is not limited in this embodiment, but the number of edge layer images on one side of each section generally does not exceed 3.
S102, extracting non-low frequency information from an edge layer image of a CT image, and obtaining an edge layer image only containing the low frequency information;
wherein the non-low frequency information includes artifact information and non-artifact information.
In some embodiments, the extracting non-low frequency information from the edge layer image of the CT image to obtain the edge layer image only including the low frequency information may include:
performing time domain filtering on an edge layer image of the CT image to obtain a first filtering image which is used as an edge layer image only containing low-frequency information;
and subtracting the edge layer image of the CT image from the first filtering image to obtain a first difference image containing non-low frequency information.
S103, extracting non-artifact information from the non-low frequency information;
in some embodiments, the extracting the non-artifact information from the non-low frequency information may include:
performing spatial domain filtering on the first difference image to obtain a second filtered image only containing artifact information;
and subtracting the first difference image from the second filtered image to obtain a second difference image containing non-artifact information.
And S104, compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image.
In some embodiments, the compensating the edge layer image including only the low frequency information by the non-artifact information to obtain the corrected edge layer image may include:
and adding the first filtered image and the second difference image to obtain a corrected edge layer image.
In a possible implementation manner, the performing temporal filtering on the edge layer image of the CT image to obtain a first filtered image may include:
and performing time domain filtering on the edge layer image of the CT image by adopting a one-dimensional bilateral filter in the CT scanning direction to obtain a first filtered image.
The weight of the one-dimensional bilateral filter can be calculated according to a first specified formula (1), and the first specified formula is as follows:
wherein n is the coordinate in a one-dimensional convolution window of the CT scanning direction, z is the central coordinate of the one-dimensional convolution window, g (z) is the CT value corresponding to the z point, and sigma s1 Sum sigma r1 For the filter parameters, by adjusting the filter parameters sigma s1 、σ r1 The filtered edge layer image can be made to contain only low frequency information. And carrying out normalization processing on the calculated weight.
The above-mentioned temporal filtering of the edge layer image of the CT image may be understood as filtering the edge layer image using other layer images (for example, images of the number of half-convolution window widths above and below the edge layer image) adjacent to the edge layer image in the CT scanning direction.
For example: time domain filtering is performed on the No. 32 tomogram obtained by the first circle of scanning (the No. 32 tomogram is an edge layer image), and assuming that a one-dimensional convolution window contains 11 tomograms, the No. 27-32 tomograms obtained by the first circle of scanning and the No. 1-5 tomograms obtained by the second circle of scanning can be used for filtering the No. 32 tomogram obtained by the first circle of scanning.
In a possible implementation manner, the performing spatial domain filtering on the first difference image to obtain a second filtered image that only includes artifact information may include:
and performing spatial domain filtering on the first difference image by adopting a two-dimensional bilateral filter to obtain a second filtered image only containing artifact information.
The weight of the two-dimensional bilateral filter can be calculated according to a second specified formula (2), and the second specified formula is:
wherein (i, j) is the coordinate in the two-dimensional convolution window, (x, y) is the central coordinate of the two-dimensional convolution window, g (x, y) is the CT value corresponding to the (x, y) point, sigma s2 Sum sigma r2 For the filter parameters, by adjusting the filter parameters sigma s2 、σ r2 The second filtered image obtained by filtering the first difference image may be made to contain only artifact information. And carrying out normalization processing on the calculated weight.
It should be noted that, the above spatial domain filtering of the first difference image may be understood as that the pixel point in the image is filtered by the pixel point adjacent to the pixel point.
For a more visual understanding of the embodiments of the present application, the following is a specific example:
assuming that an edge layer image of a CT image contains basic tissue information, bone medium information, detail tissue information, high-frequency noise information and artifact information, after time domain filtering, the obtained first filter image only contains the basic tissue information and the bone medium information, then subtracting the edge layer image of the CT image from the first filter image, the obtained first difference image contains the detail tissue information, the high-frequency noise information and the artifact information, then carrying out space domain filtering on the first difference image, the obtained second filter image only contains the artifact information, then subtracting the first difference image from the second filter image, the obtained second difference image contains the detail tissue information and the high-frequency noise information, finally adding the first filter image and the second difference image (namely, recombining the first filter image only containing the basic tissue information and the bone medium information with the second difference image containing the detail tissue information and the high-frequency noise information) to obtain a corrected edge layer image, and the corrected edge layer image does not contain artifact information, so that the integral difference between the edge layer image and other layer images can be eliminated, and the problem of bright edge layer image or the MPR can be obviously improved when the MPR is carried out.
In addition, it should be noted that, through a large number of experimental verification, the correction method of the CT image provided in the embodiment of the present application is proved to be effective. The following shows a set of experimental results as shown in fig. 2 and 3, where fig. 2 is an MPR image before correction formed by MPR based on an original edge layer image, and as can be seen from fig. 2, there is a distinct streak artifact in the MPR image before correction (as shown by a dashed box in fig. 2), and fig. 3 is an MPR image after correction formed by MPR based on a corrected edge layer image obtained by a correction method for CT image provided by an embodiment of the present application, and as can be seen from fig. 3, the streak artifact in the MPR image after correction is significantly improved, and detailed tissue information in the image is also well preserved, so that the correction method for CT image provided by an embodiment of the present application can effectively remove the streak artifact in the MPR image.
Based on the same inventive concept, referring to fig. 4, an embodiment of the present application further provides a correction device for a CT image, including: a CT image acquisition module 11, a first extraction module 12, a second extraction module 13, and a correction module 14.
Wherein, the CT image acquisition module 11 is used for acquiring a CT image obtained by scanning the object;
a first extraction module 12 for extracting non-low frequency information from an edge layer image of the CT image and obtaining an edge layer image containing only low frequency information;
a second extraction module 13 for extracting non-artifact information from the non-low frequency information;
the correction module 14 is configured to compensate the edge layer image only including the low frequency information through the non-artifact information, so as to obtain a corrected edge layer image.
In one possible implementation, the first extraction module 12 may include:
the time domain filtering sub-module is used for performing time domain filtering on the edge layer image of the CT image to obtain a first filtering image which is used as the edge layer image only containing low-frequency information;
and the first image subtraction sub-module is used for subtracting the edge layer image of the CT image from the first filtering image to obtain a first difference image containing non-low frequency information.
In a possible implementation, the second extraction module 13 may include:
the spatial domain filtering sub-module is used for performing spatial domain filtering on the first difference image to obtain a second filtered image only containing artifact information;
the second image subtraction sub-module is used for subtracting the first difference image from the second filter image to obtain a second difference image containing non-artifact information;
the correction module 14 may be specifically configured to:
and adding the first filtered image and the second difference image to obtain a corrected edge layer image.
In one possible implementation, the temporal filtering submodule may be specifically configured to:
and performing time domain filtering on the edge layer image of the CT image by adopting a one-dimensional bilateral filter in the CT scanning direction to obtain a first filtered image.
The weight of the one-dimensional bilateral filter can be calculated according to the first specified formula (1).
In one possible implementation, the spatial domain filtering submodule may be specifically configured to:
and performing spatial domain filtering on the first difference image by adopting a two-dimensional bilateral filter to obtain a second filtered image only containing artifact information.
The weight of the two-dimensional bilateral filter can be calculated according to the second specified formula (2).
The implementation process of the functions and roles of each unit in the above device is specifically shown in the implementation process of the corresponding steps in the above method, and will not be described herein again.
For the device embodiments, reference is made to the description of the method embodiments for the relevant points, since they essentially correspond to the method embodiments. The apparatus embodiments described above are merely illustrative, wherein the elements illustrated as separate elements may or may not be physically separate, and the elements shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purposes of the present application. Those of ordinary skill in the art will understand and implement the present invention without undue burden.
Based on the same inventive concept, the embodiments of the present application also provide a storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the method for correcting a CT image in any of the possible implementations described above.
Alternatively, the storage medium may be a memory.
Based on the same inventive concept, referring to fig. 5, the embodiment of the present application further provides an electronic device, including a memory 71 (e.g. a non-volatile memory), a processor 72, and a computer program stored on the memory 71 and executable on the processor 72, wherein the processor 72 implements the steps of the method for correcting a CT image in any of the possible implementations described above when executing the program. The electronic device may be, for example, a PC.
As shown in fig. 5, the electronic device may generally further include: memory 73, network interface 74, and internal bus 75. In addition to these components, other hardware may be included, which is not described in detail.
It should be noted that the above-mentioned correction device for CT images may be implemented by software, which is a device in a logic sense, and is formed by reading, by the processor 72 of the electronic device in which the correction device is located, the computer program instructions stored in the nonvolatile memory into the memory 73 for execution.
Embodiments of the subject matter and the functional operations described in this specification can be implemented in: digital electronic circuitry, tangibly embodied computer software or firmware, computer hardware including the structures disclosed in this specification and structural equivalents thereof, or a combination of one or more of them. Embodiments of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on a tangible, non-transitory program carrier for execution by, or to control the operation of, data processing apparatus. Alternatively or additionally, the program instructions may be encoded on a manually-generated propagated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode and transmit information to suitable receiver apparatus for execution by data processing apparatus. The computer storage medium may be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them.
The processes and logic flows described in this specification can be performed by one or more programmable computers executing one or more computer programs to perform corresponding functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application-specific integrated circuit).
Computers suitable for executing computer programs include, for example, general purpose and/or special purpose microprocessors, or any other type of central processing unit. Typically, the central processing unit will receive instructions and data from a read only memory and/or a random access memory. The essential elements of a computer include a central processing unit for carrying out or executing instructions and one or more memory devices for storing instructions and data. Typically, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto-optical disks, or optical disks, etc. However, a computer does not have to have such a device. Furthermore, the computer may be embedded in another device, such as a mobile phone, a Personal Digital Assistant (PDA), a mobile audio or video player, a game console, a Global Positioning System (GPS) receiver, or a portable storage device such as a Universal Serial Bus (USB) flash drive, to name a few.
Computer readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media and memory devices including, for example, semiconductor memory devices (e.g., EPROM, EEPROM, and flash memory devices), magnetic disks (e.g., internal hard disk or removable disks), magneto-optical disks, and CD-ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
While this specification contains many specific implementation details, these should not be construed as limitations on the scope of any invention or of what may be claimed, but rather as descriptions of features of specific embodiments of particular inventions. Certain features that are described in this specification in the context of separate embodiments can also be implemented in combination in a single embodiment. On the other hand, the various features described in the individual embodiments may also be implemented separately in the various embodiments or in any suitable subcombination. Furthermore, although features may be acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a subcombination or variation of a subcombination.
Similarly, although operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results. In some cases, multitasking and parallel processing may be advantageous. Moreover, the separation of various system modules and components in the embodiments described above should not be understood as requiring such separation in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
Thus, particular embodiments of the subject matter have been described. Other embodiments are within the scope of the following claims. In some cases, the actions recited in the claims can be performed in a different order and still achieve desirable results. Furthermore, the processes depicted in the accompanying drawings are not necessarily required to be in the particular order shown, or sequential order, to achieve desirable results. In some implementations, multitasking and parallel processing may be advantageous.
The foregoing description of the preferred embodiments of the present invention is not intended to limit the invention to the precise form disclosed, and any modifications, equivalents, improvements and alternatives falling within the spirit and principles of the present invention are intended to be included within the scope of the present invention.
Claims (8)
1. A method of correcting a CT image, the method comprising:
obtaining a CT image obtained by scanning a subject;
extracting non-low frequency information from an edge layer image of the CT image, and obtaining an edge layer image only containing low frequency information;
extracting non-artifact information from the non-low frequency information;
compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image;
the extracting non-low frequency information from the edge layer image of the CT image and obtaining the edge layer image only containing the low frequency information comprises the following steps:
performing time domain filtering on the edge layer image of the CT image to obtain a first filtering image which is used as the edge layer image only containing low-frequency information;
subtracting the edge layer image of the CT image from the first filtering image to obtain a first difference image containing non-low frequency information;
the extracting non-artifact information from the non-low frequency information includes:
performing spatial domain filtering on the first difference image to obtain a second filtered image only containing artifact information;
subtracting the first difference image from the second filtered image to obtain a second difference image containing non-artifact information;
the compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image comprises the following steps:
adding the first filtered image and the second difference image to obtain a corrected edge layer image;
the spatial domain filtering is performed on the first difference image to obtain a second filtered image only containing artifact information, including:
performing spatial domain filtering on the first difference image by adopting a two-dimensional bilateral filter to obtain a second filtered image only containing artifact information;
the performing temporal filtering on the edge layer image of the CT image to obtain a first filtered image includes:
and filtering the edge layer image by adopting other layer images adjacent to the edge layer image in the CT scanning direction to obtain the first filtered image.
2. The method of claim 1, wherein the weights of the one-dimensional bilateral filter are calculated according to a first specified formula:
wherein n is the coordinate in a one-dimensional convolution window of the CT scanning direction, z is the central coordinate of the one-dimensional convolution window, g (z) is the CT value corresponding to the z point, and sigma s1 Sum sigma r1 Is a filter parameter.
3. The method of claim 1, wherein the weights of the two-dimensional bilateral filter are calculated according to a second specified formula, the second specified formula being:
wherein (i, j) is the coordinate in the two-dimensional convolution window, (x, y) is the central coordinate of the two-dimensional convolution window, g (x, y) is the CT value corresponding to the (x, y) point, sigma s2 Sum sigma r2 Is a filter parameter.
4. A correction device for CT images, said device comprising:
the CT image acquisition module is used for acquiring a CT image obtained by scanning the object;
a first extraction module, configured to extract non-low frequency information from an edge layer image of the CT image, and obtain an edge layer image only containing low frequency information;
the second extraction module is used for extracting non-artifact information from the non-low frequency information;
the correction module is used for compensating the edge layer image only containing the low-frequency information through the non-artifact information to obtain a corrected edge layer image;
the first extraction module includes:
the time domain filtering sub-module is used for performing time domain filtering on the edge layer image of the CT image to obtain a first filtering image which is used as the edge layer image only containing low-frequency information;
a first image subtraction sub-module, configured to subtract the edge layer image of the CT image from the first filtered image, to obtain a first difference image that includes non-low frequency information;
the second extraction module includes:
the spatial domain filtering sub-module is used for performing spatial domain filtering on the first difference image to obtain a second filtering image only containing artifact information;
a second image subtraction sub-module, configured to subtract the first difference image from the second filtered image to obtain a second difference image that includes non-artifact information;
the correction module is specifically configured to:
adding the first filtered image and the second difference image to obtain a corrected edge layer image;
the spatial domain filtering submodule is specifically used for:
performing spatial domain filtering on the first difference image by adopting a two-dimensional bilateral filter to obtain a second filtered image only containing artifact information;
the time domain filtering sub-module is specifically configured to:
and filtering the edge layer image by adopting other layer images adjacent to the edge layer image in the CT scanning direction to obtain the first filtered image.
5. The apparatus of claim 4, wherein the weights of the one-dimensional bilateral filter are calculated according to a first specified formula:
wherein n is the coordinate in a one-dimensional convolution window of the CT scanning direction, z is the central coordinate of the one-dimensional convolution window, g (z) is the CT value corresponding to the z point, and sigma s1 Sum sigma r1 Is a filter parameter.
6. The apparatus of claim 4, wherein the weights of the two-dimensional bilateral filter are calculated according to a second specified formula, the second specified formula being:
wherein (i, j) is the coordinate in the two-dimensional convolution window, (x, y) is the central coordinate of the two-dimensional convolution window, g (x, y) is the CT value corresponding to the (x, y) point, sigma s2 Sum sigma r2 Is a filter parameter.
7. A storage medium having stored thereon a computer program, which when executed by a processor performs the steps of the method of any of claims 1-3.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the method of any of claims 1-3 when the program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910349470.1A CN110070589B (en) | 2019-04-28 | 2019-04-28 | CT image correction method and device, storage medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910349470.1A CN110070589B (en) | 2019-04-28 | 2019-04-28 | CT image correction method and device, storage medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110070589A CN110070589A (en) | 2019-07-30 |
| CN110070589B true CN110070589B (en) | 2024-03-22 |
Family
ID=67369267
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910349470.1A Active CN110070589B (en) | 2019-04-28 | 2019-04-28 | CT image correction method and device, storage medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110070589B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111798391A (en) * | 2020-06-29 | 2020-10-20 | 东软医疗系统股份有限公司 | Image processing method and device, medical imaging equipment and system |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20160012559A (en) * | 2014-07-24 | 2016-02-03 | 삼성전자주식회사 | Magnetic resonance imaging apparatus and imaging method for magnetic resonance image thereof |
| CN104504655A (en) * | 2014-12-04 | 2015-04-08 | 沈阳东软医疗系统有限公司 | CT (computer tomography) metal artifact processing method and device |
| WO2016158138A1 (en) * | 2015-04-01 | 2016-10-06 | 株式会社日立製作所 | X-ray ct apparatus, reconfiguration arithmetic apparatus, and x-ray ct image generation method |
| CN107203983B (en) * | 2016-03-17 | 2024-03-22 | 通用电气公司 | Method and system for reducing grid line artifacts in X-ray images |
-
2019
- 2019-04-28 CN CN201910349470.1A patent/CN110070589B/en active Active
Non-Patent Citations (2)
| Title |
|---|
| 基于双边滤波迭代修正的CT欠投影ART重建;齐宏亮 等;《生物医学工程学杂志》;20130430;第30卷(第2期);第1.2节 * |
| 改进Canny的核磁共振图像伪影校正算法;张琳梅，赵莉;《南京理工大学学报》;20160228;第40卷(第1期);第1.3.1节，第2节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110070589A (en) | 2019-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8891885B2 (en) | Method, computing unit, CT system and C-arm system for reducing metal artifacts in CT image datasets | |
| Zheng et al. | Magnetic resonance–based automatic air segmentation for generation of synthetic computed tomography scans in the head region | |
| US20200349449A1 (en) | 3-d convolutional autoencoder for low-dose ct via transfer learning from a 2-d trained network | |
| US10964072B2 (en) | Methods, systems, and media for noise reduction in computed tomography images | |
| CN105528766B (en) | CT metal artifacts treating method and apparatus | |
| US8139835B2 (en) | Method for noise reduction in digital images with locally different and directional noise | |
| CN109171727B (en) | Magnetic resonance imaging method and device | |
| CN106530236B (en) | Medical image processing method and system | |
| Manjón et al. | Multicomponent MR image denoising | |
| CN108287324B (en) | Reconstruction method and device of magnetic resonance multi-contrast image | |
| US9646366B2 (en) | Method and apparatus for enhancing medical images | |
| CN107622481B (en) | Method and device for reducing CT image noise and computer equipment | |
| US11360180B2 (en) | Methods, systems, and computer readable media for using a trained adversarial network for performing retrospective magnetic resonance imaging (MRI) artifact correction | |
| US20170172534A1 (en) | Thoracic imaging for cone beam computed tomography | |
| CN110070589B (en) | CT image correction method and device, storage medium and electronic equipment | |
| WO2017190968A1 (en) | Device and method for denoising a vector-valued image | |
| CN111091517A (en) | Residual weighted imaging method and device | |
| CN108351396B (en) | Method, computer program product and magnetic resonance imaging system for tissue classification | |
| WO2018108848A1 (en) | Information-adaptive regularization for iterative image reconstruction | |
| Tian et al. | Aliasing artifact reduction in spiral real‐time MRI | |
| Pham | Supervised restoration of degraded medical images using multiple-point geostatistics | |
| US9436978B2 (en) | Information processing apparatus, information processing method, and storage medium | |
| CN102890817A (en) | Method and system for processing medical images | |
| Anas et al. | High-quality 3D correction of ring and radiant artifacts in flat panel detector-based cone beam volume CT imaging | |
| Yim et al. | A deep convolutional neural network for simultaneous denoising and deblurring in computed tomography |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |