CN113808232A - Grid artifact suppression method and device, electronic device and storage medium - Google Patents
Grid artifact suppression method and device, electronic device and storage medium Download PDFInfo
- Publication number
- CN113808232A CN113808232A CN202111105620.8A CN202111105620A CN113808232A CN 113808232 A CN113808232 A CN 113808232A CN 202111105620 A CN202111105620 A CN 202111105620A CN 113808232 A CN113808232 A CN 113808232A
- Authority
- CN
- China
- Prior art keywords
- image
- grid
- radiographic image
- frequency domain
- region
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration by non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Abstract
The application provides a grid artifact suppression method, a grid artifact suppression device, an electronic device and a storage medium, wherein the method comprises the following steps: acquiring a radiographic image, which is an image of the X-ray after passing through the grid and the region of the object to be inspected; fourier transform is carried out on the ray image to obtain a transformed ray image; filtering and inhibiting the grid artifacts in the transformed radiographic image by using an adaptive filter to obtain an inhibited radiographic image; and performing inverse Fourier transform on the suppressed radiographic image to obtain a result image. The method comprises the steps of carrying out Fourier transform on a ray image to obtain a transformed ray image, filtering grid artifacts in the transformed ray image by using an adaptive filter, and carrying out inverse Fourier transform on the filtered ray image to effectively filter and inhibit the grid artifacts in the ray image.
Description
Technical Field
The present application relates to the technical field of medical equipment, image processing, and grid artifact suppression, and in particular, to a grid artifact suppression method, an apparatus, an electronic device, and a storage medium.
Background
The grid is a medical instrument part used between the front of a detector and the back of a human body and is mainly used for filtering scattered rays, reducing gray haze and improving image contrast.
At present, in consideration of economical efficiency, a static low-density grid is generally used in a Digital Radiography (DR) system, and when X-rays pass through the static low-density grid, streak artifacts of the grid with alternate light and shade are generated on an image.
Disclosure of Invention
An object of the embodiments of the present application is to provide a method, an apparatus, an electronic device, and a storage medium for suppressing artifacts caused by light and dark interleaving of a low-density grid.
The embodiment of the application provides a grid artifact suppression method, which comprises the following steps: acquiring a radiographic image, which is an image of the X-ray after passing through the grid and the region of the object to be inspected; fourier transform is carried out on the ray image to obtain a transformed ray image; filtering and inhibiting the grid artifacts in the transformed radiographic image by using an adaptive filter to obtain an inhibited radiographic image; and performing inverse Fourier transform on the suppressed radiographic image to obtain a result image. In the implementation process, the Fourier transform is performed on the radiographic image to obtain the transformed radiographic image, the adaptive filter is used for filtering and suppressing the grid artifacts in the transformed radiographic image, and the inverse Fourier transform is performed on the suppressed radiographic image to effectively filter and suppress the grid artifacts in the radiographic image, so that the quality of the obtained radiographic image is improved.
Optionally, in an embodiment of the present application, after acquiring the radiographic image and before performing fourier transform on the radiographic image, the method further includes: judging whether a grid artifact exists in the radiographic image or not; and if so, performing Fourier transform on the ray image. In the implementation process, whether the grid artifact exists in the radiographic image or not is judged firstly, and if the grid artifact exists, the grid artifact is further removed, so that the situation that the grid artifact is still removed due to the fact that the grid is removed and moved away, and a vibrating grid or a high-density grid is used, and the grid artifact does not exist is avoided, and image information is effectively protected and computing resources are saved.
Optionally, in an embodiment of the present application, before filtering the grid artifact in the transformed radiographic image by using the adaptive filter, the method further includes: acquiring a region of interest with grid artifacts in radiographic images; carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest; calculating the frequency domain characteristics of the grid artifacts in the frequency domain image; an adaptive filter is constructed from the frequency domain features in the frequency domain image. In the implementation process, a region of interest with grid artifacts in radiographic images is acquired; carrying out Fourier transform on the region of interest, and calculating the frequency domain characteristics of the grid artifact in a frequency domain image; and then, an adaptive filter is constructed according to the frequency domain features in the frequency domain image, so that the filter constructed according to the frequency domain features can be used for effectively filtering and inhibiting grid artifacts in the radiographic image, and the quality of the acquired radiographic image is improved.
Optionally, in an embodiment of the present application, acquiring a region of interest in which grid artifacts are present in radiographic images includes: and selecting the region of interest from the radiographic image according to a preset rule, wherein the preset rule comprises the following steps: and (4) dividing lines identified by shutter software and/or coordinate parameters sent by shutter hardware. In the implementation process, by acquiring the region of interest in which the grid artifact exists in the radiographic image and constructing the adaptive filter according to the spectral characteristics of the region of interest, the filter constructed by using the region of interest can effectively filter and suppress the grid artifact in the radiographic image, and the quality of the acquired radiographic image is improved.
Optionally, in this embodiment of the present application, before performing fourier transform on the region of interest, the method further includes: and performing boundary clipping, mirror image continuation, spatial filtering, spatial superposition, multi-scale decomposition and/or windowing on the region of interest and the like. In the implementation process, the region of interest is preprocessed, so that the frequency spectrum leakage is reduced, the processing speed and the frequency point positioning precision are improved, and the details of the image are effectively protected.
Optionally, in an embodiment of the present application, the frequency domain features include: frequency range, curve shape and image direction; calculating a frequency domain characteristic of the grid artifact in the frequency domain image, comprising: counting the frequency range of the grid artifact in the frequency domain image; performing curve fitting on the frequency domain image to obtain the curve shape of the grid artifact in the frequency domain image; the image orientation of the grid with respect to the detector is determined from the position and the characteristics of the grid artifact in the frequency domain image.
Optionally, in this embodiment of the present application, before performing fourier transform on the ray image, the method further includes: and performing boundary clipping, mirror image continuation, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing on the ray image. In the implementation process, the ray image is preprocessed, so that the frequency spectrum leakage is reduced, the processing speed and the frequency point positioning precision are improved, and the details of the image are effectively protected.
Optionally, in an embodiment of the present application, the adaptive filter includes: a gaussian filter or a butterworth filter.
An embodiment of the present application further provides a grid artifact suppression apparatus, including: a radiographic image acquisition module for acquiring a radiographic image which is an image of the X-ray after passing through the grid and the region of the object to be inspected; the radiographic image transformation module is used for carrying out Fourier transformation on the radiographic image to obtain a transformed radiographic image; a radiographic filtering module, configured to perform filtering suppression on the grid artifacts in the transformed radiographic image by using an adaptive filter, so as to obtain a suppressed radiographic image; and the radiographic image inverse transformation module is used for performing inverse Fourier transform on the suppressed radiographic image to obtain a result image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes: the radiographic image judging module is used for judging whether a grid artifact exists in the radiographic image or not; and the image transformation executing module is used for executing Fourier transformation on the ray image if the grid artifact exists in the ray image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes: the interest region acquisition module is used for acquiring an interest region with grid artifacts in the radiographic image; the frequency domain image obtaining module is used for carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest; the frequency domain characteristic calculation module is used for calculating the frequency domain characteristics of the grid artifacts in the frequency domain image; and the adaptive filter module is used for constructing an adaptive filter according to the frequency domain characteristics in the frequency domain image.
Optionally, in an embodiment of the present application, the interest region obtaining module includes: the region of interest selecting module is used for selecting the region of interest from the radiographic image according to preset rules, and the preset rules comprise: and (4) dividing lines identified by shutter software and/or coordinate parameters sent by shutter hardware.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes: and the interest region preprocessing module is used for performing boundary clipping, mirror extension, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing and the like on the interest region.
Optionally, in an embodiment of the present application, the frequency domain features include: frequency range, curve shape and image direction; a frequency domain feature computation module comprising: the frequency range counting module is used for counting the frequency range of the grid artifact in the frequency domain image; the curve shape obtaining module is used for performing curve fitting on the frequency domain image to obtain the curve shape of the grid artifact in the frequency domain image; and the image direction determining module is used for determining the image direction of the grid relative to the detector according to the position and the characteristics of the grid artifact in the frequency domain image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes: a radiographic image preprocessing module for performing boundary clipping, mirror extension, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing, etc. on the radiographic image.
Optionally, in an embodiment of the present application, the adaptive filter includes: a gaussian filter or a butterworth filter.
An embodiment of the present application further provides an electronic device, including: a processor and a memory, the memory storing processor-executable machine-readable instructions, the machine-readable instructions when executed by the processor performing the method as described above.
Embodiments of the present application also provide a computer-readable storage medium having a computer program stored thereon, where the computer program is executed by a processor to perform the method as described above.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and that those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a schematic flow chart of a grid artifact reduction method provided by an embodiment of the present application;
FIG. 2 is a schematic diagram of a flow chart for constructing an adaptive filter according to an embodiment of the present application;
fig. 3 is a schematic structural diagram of a grid artifact reduction apparatus provided in an embodiment of the present application;
fig. 4 is a schematic structural diagram of an electronic device provided in an embodiment of the present application.
Detailed Description
The technical solution 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.
Before introducing the grid artifact reduction method provided by the embodiments of the present application, some concepts involved in the embodiments of the present application are introduced:
it should be noted that the grid artifact reduction method provided in the embodiments of the present application may be executed by an electronic device, where the electronic device refers to a device terminal having a function of executing a computer program or the server described above, and the device terminal includes, for example: a smart phone, a personal computer, a tablet computer, a personal digital assistant, or a mobile internet device, etc. The server is, for example: x86 server and non-x 86 server, non-x 86 server includes: mainframe, minicomputer, and UNIX server. Of course, in the specific implementation process, a Graphics Processing Unit (GPU) may be used to accelerate the execution speed of the above-mentioned grid artifact reduction method, or the GPU and a Central Processing Unit (CPU) may be used to execute the grid artifact reduction method together, or only the CPU is used to execute the grid artifact reduction method.
Application scenarios to which the grid artifact reduction method is applicable are described below, where the application scenarios include, but are not limited to: it can be noted that when the X-ray is in human body interaction, part of the X-ray can penetrate the human body, another part of the energy can generate photoelectric effect, and another part of the X-ray can generate Compton effect when passing through the grid, so that the original propagation direction is changed, and the scattered ray is formed; the use of a grid can effectively reduce scattered radiation and improve the contrast of the image, but in the process of using a low-density grid, grid artifacts between light and dark can be introduced into the X-ray image. In this case, the grid artifact or the like in the X-ray image can be effectively removed by the grid artifact reduction method, thereby improving the quality of the obtained radiographic image. The grid artifact mitigation method may also be used to enhance the functionality of medical software systems, such as in particular: the grid artifact suppression method is used for enhancing the function of a Digital Radiography (DR) system, so that the DR system can effectively acquire an X-ray image without grid artifacts, the bright-dark staggered artifacts generated by a low-density grid are effectively suppressed, and the quality of the X-ray image is improved.
It can be appreciated that the grid artifact reduction method has wide applicability to different detectors and grid specifications, i.e. the method can be used for grid artifact reduction of X-ray images generated by different types of detectors and grids of different specifications, and can well reduce grid artifacts while maintaining detail and balance ringing.
Please refer to fig. 1, which is a schematic flow chart of a grid artifact reduction method provided in the embodiment of the present application; the main idea of the grid artifact suppression method is that Fourier transformation is carried out on a ray image to obtain a transformed ray image, filtering suppression is carried out on grid artifacts in the transformed ray image by using an adaptive filter, and then inverse Fourier transformation is carried out on the ray image after filtering suppression, so that the grid artifacts in the ray image are effectively filtered and suppressed, and the quality of the obtained ray image is improved. The grid artifact reduction method described above may include:
step S110: a radiographic image is acquired, which is an image of the X-rays after they have passed through the grid and the area of the object to be examined.
The method for acquiring the radiographic image comprises the following steps: a first acquisition mode in which an X-ray image of an inspection object is acquired by a transmitting and receiving device such as a bulb or an X-ray detector; then the terminal device sends the radiographic image to the electronic device, then the electronic device receives the radiographic image sent by the terminal device, and the electronic device can store the radiographic image into a file system, a database or a mobile storage device; a second acquisition mode, which acquires a pre-stored radiographic image, specifically includes: obtaining a radiographic image from a file system, or from a database, or from a mobile storage device; in the third mode, the radiographic image is acquired by using software such as a browser or the like, or by accessing the internet by using other application programs.
Optionally, after acquiring the radiographic image, the grid artifact may be removed only when the grid artifact is determined to be present in the radiographic image, for example: judging whether a grid artifact exists in the radiographic image or not; if it is determined that the grid artifact exists in the radiological image, the following step S120 is performed.
Step S120: fourier transform is performed on the radiographic image to obtain a transformed radiographic image.
Optionally, before fourier transforming the ray image, preprocessing is further required to be performed on the ray image, and an embodiment of the preprocessing may include: and (3) performing boundary clipping, mirror image continuation, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing on the ray image. And when the ROI area is larger than the empirical threshold, performing boundary clipping, mirror image continuation or copy continuation.
The embodiment of step S120 described above is, for example: one-dimensional fourier transform equations may be used
Or a two-dimensional Fourier transform formula
Fourier transform is carried out on the ray image to obtain a transformed ray image; where f (u) represents the result of one-dimensional Fourier transform, and f (u, v) represents the result of two-dimensional Fourier transform.
Step S130: and filtering and suppressing the grid artifact in the transformed radiographic image by using an adaptive filter to obtain a suppressed radiographic image.
The above-mentioned step S130 can be implemented in two ways: in a first embodiment, an adaptive filter constructed in advance is obtained, and filtering suppression is performed by using the adaptive filter, and the first embodiment may include: a pre-constructed adaptive filter is acquired, and the adaptive filter is used to filter and suppress grid artifacts in the transformed radiograph, thereby obtaining a suppressed radiograph. Wherein the adaptive filter includes: a gaussian filter or a butterworth filter. In the second embodiment, each time the adaptive filter is constructed and used for filtering, the process of constructing the adaptive filter will be described in detail below.
Step S140: and performing inverse Fourier transform on the suppressed radiographic image to obtain a result image.
The embodiment of step S140 described above is, for example: the filtered and suppressed radiographic image is subjected to Inverse fourier transform (Inverse fourier transform), and the filtered and suppressed radiographic image is converted from the frequency domain to the spatial domain, thereby obtaining a resultant image.
In the implementation process, the Fourier transform is performed on the radiographic image to obtain the transformed radiographic image, the adaptive filter is used for filtering the grid artifacts in the transformed radiographic image, and the inverse Fourier transform is performed on the filtered radiographic image to effectively filter and suppress the grid artifacts in the radiographic image, so that the quality of the obtained radiographic image is improved.
Please refer to fig. 2, which is a schematic flow chart of the adaptive filter construction provided in the embodiment of the present application; it is understood that before the adaptive filter is used, the adaptive filter is required to be constructed, and the process of constructing the adaptive filter may include:
step S210: a region of interest in which grid artifacts are present in the radiographs is acquired.
The Region of interest (ROI) is a Region that needs to be processed and is delineated from a processed image in a form of a box, a circle, an ellipse, an irregular polygon, or the like in machine vision and image processing.
The embodiment of step S210 described above is, for example: selecting the region of interest from the radiographic image according to a preset rule, wherein the preset rule can be obtained through configuration items in a configuration file, and the method specifically comprises the following steps: the cut line identified by the shutter software and/or the coordinate parameters sent by the shutter hardware, etc. Here, the ROI is taken as a rectangular region for example, and a rectangular target region, Rect (x, y), is obtained in a coordinate system with the upper left corner of the radiographic image as the origin, where point (x, y) is the coordinate of the upper left corner of the rectangle, x is the width of the rectangle, and y is the height of the rectangle. The ROI can be transmitted by a shutter hardware coordinate parameter, or can be used for carrying out target identification on a ray image through software on the shutter, so that a foreground region and a background region are segmented.
Step S220: and carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest.
Optionally, before performing fourier transform on the region of interest, the region of interest may be further preprocessed, and an embodiment of the preprocessing may include: performing boundary clipping, mirror extension, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing on the region of interest; and when the ROI area is larger than the empirical threshold, performing boundary clipping, mirror image continuation or copy continuation.
The embodiment of step S220 described above is, for example: one-dimensional fourier transform equations may be used
Or a two-dimensional Fourier transform formula
Carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest; where f (u) represents the result of one-dimensional Fourier transform, and f (u, v) represents the result of two-dimensional Fourier transform.
Step S230: the frequency domain characteristics of the grid artifacts in the frequency domain image are calculated.
Wherein the frequency domain features include: the frequency range, the curve shape and the image direction are described below in these three aspects, respectively.
In a first aspect, the frequency range of grid artifacts in a frequency domain image is calculated; specific examples thereof include: assuming that a pixel point x in the frequency domain image is represented by u (x) ═ x-a, x + a, the domain of the fourier transform of the pixel point x can be represented by f (x), the frequency range of the grid in the frequency domain image (referred to as the range) can be counted by searching from the high frequency to the low frequency according to the formula f (x) ═ max (u (x + a))) & & avg (u (x))) > avg (u (x-a))) × weight & & max (u (x)) & max (u (x-a &) and the frequency point (x)) & max (u (x) + max (u (x + a)) and the position and the characteristic of the grid in the frequency domain image can be obtained. Wherein max represents the maximum value on the calculation bit, avg represents the bit mean calculation, and weight is the weight value.
In a second aspect, the curve shape of the grid artifact in the frequency domain image is calculated; specific examples thereof include: and performing curve fitting on the frequency domain image to obtain the curve shape of the grid artifact in the frequency domain image.
In a third aspect, the image direction of the grid artifact in the frequency domain image is calculated; specific examples thereof include: the image orientation of the grid with respect to the detector is determined from the position and the characteristics of the grid artifact in the frequency domain image. The judgment of the image direction is divided into judgment during one-dimensional Fourier transform and judgment during two-dimensional Fourier transform, the judgment modes of the two are different, and the judgment needs to be carried out according to a series of relevant rules of frequency domain position and frequency domain characteristics.
Step S240: an adaptive filter is constructed from the frequency domain features in the frequency domain image.
The embodiment of step S240 described above is, for example: constructing an adaptive filter according to the frequency domain characteristics in the frequency domain image; among these, adaptive filters include but are not limited to: gaussian filters and butterworth filters, etc. Here, a gaussian filter is taken as an example, and the gaussian filter may be expressed as f (x) ═ a (exp (- (x-b)2/2c2) Where a, b and c are all coefficients of the adaptive filter, assuming that the resulting coefficient of the adaptive filter is denoted as w, the resulting adaptive filter can be denoted as filter (w). The process of filtering above may be expressed using the formula y (w) ═ x (w) x filter (w); wherein w is the coefficient of the adaptive filter, y (w) is the filtering result, x (w) is the result of performing fast fourier transform on the transformed radiographic image, and filter (w) is the adaptive filter.
In the implementation process, a region of interest with grid artifacts in radiographic images is acquired; carrying out Fourier transform on the region of interest, and calculating the frequency domain characteristics of the grid artifact in a frequency domain image; and then constructing an adaptive filter according to the frequency domain features in the frequency domain image, so that the filter constructed by using the targeted frequency domain features can effectively filter and inhibit grid artifacts in the radiographic image, thereby improving the quality of the acquired radiographic image.
Please refer to fig. 3, which is a schematic structural diagram of a grid artifact reduction apparatus provided in an embodiment of the present application; the embodiment of the present application provides a grid artifact suppression apparatus 300, including:
a radiograph acquisition module 310 for acquiring a radiograph, which is an image after X-rays pass through the grid and the region of the object to be examined.
And a radiation image transforming module 320 for performing fourier transform on the radiation image to obtain a transformed radiation image.
A radiograph filtering module 330 for performing filter suppression on the grid artifacts in the transformed radiograph using an adaptive filter to obtain a suppressed radiograph.
And the inverse radiographic image transformation module 340 is configured to perform inverse fourier transformation on the suppressed radiographic image to obtain a result image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes:
and the radiographic image judging module is used for judging whether the radiographic image has grid artifacts.
And the image transformation executing module is used for executing Fourier transformation on the ray image if the grid artifact exists in the ray image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus further includes:
and the interest region acquisition module is used for acquiring the interest region with the grid artifact in the radiographic image.
And the frequency domain image obtaining module is used for carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest.
And the frequency domain characteristic calculation module is used for calculating the frequency domain characteristics of the grid artifacts in the frequency domain image.
And the adaptive filter module is used for constructing an adaptive filter according to the frequency domain characteristics in the frequency domain image.
Optionally, in an embodiment of the present application, the interest region obtaining module includes:
the region of interest selecting module is used for selecting the region of interest from the radiographic image according to preset rules, and the preset rules comprise: and (4) dividing lines identified by shutter software and/or coordinate parameters sent by shutter hardware.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus may further include:
and the interest region preprocessing module is used for performing boundary clipping, mirror extension, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing and the like on the interest region.
Optionally, in an embodiment of the present application, the frequency domain features include: frequency range, curve shape and image direction; a frequency domain feature computation module comprising:
and the frequency range counting module is used for counting the frequency range of the grid artifact in the frequency domain image.
And the curve shape obtaining module is used for performing curve fitting on the frequency domain image to obtain the curve shape of the grid artifact in the frequency domain image.
And the image direction determining module is used for determining the image direction of the grid relative to the detector according to the position and the characteristics of the grid artifact in the frequency domain image.
Optionally, in an embodiment of the present application, the grid artifact reduction apparatus may further include:
a radiographic image preprocessing module for performing boundary clipping, mirror extension, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing, etc. on the radiographic image.
Optionally, in an embodiment of the present application, the adaptive filter includes: a gaussian filter or a butterworth filter.
It should be understood that the apparatus corresponds to the above-mentioned embodiment of the grid artifact reduction method, and can perform the steps related to the above-mentioned embodiment of the method, and the specific functions of the apparatus can be referred to the above description, and the detailed description is appropriately omitted here to avoid redundancy. The device includes at least one software function that can be stored in memory in the form of software or firmware (firmware) or solidified in the Operating System (OS) of the device.
Please refer to fig. 4 for a schematic structural diagram of an electronic device according to an embodiment of the present application. An electronic device 400 provided in an embodiment of the present application includes: a processor 410 and a memory 420, the memory 420 storing machine-readable instructions executable by the processor 410, the machine-readable instructions when executed by the processor 410 performing the method as above.
Embodiments of the present application also provide a computer-readable storage medium 430, where the computer-readable storage medium 430 stores a computer program, and the computer program is executed by the processor 410 to perform the above method. The computer-readable storage medium 430 may be implemented by any type of volatile or nonvolatile storage device or combination thereof, such as a Static Random Access Memory (SRAM), an Electrically Erasable Programmable Read-Only Memory (EEPROM), an Erasable Programmable Read-Only Memory (EPROM), a Programmable Read-Only Memory (PROM), a Read-Only Memory (ROM), a magnetic Memory, a flash Memory, a magnetic disk, or an optical disk.
Embodiments of the present application further provide a Graphics Processing Unit (GPU)440, where the memory 420 stores machine-readable instructions executable by the GPU 440, and the machine-readable instructions, when executed by the GPU 440, perform the method as above.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved.
In addition, functional modules of the embodiments in the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
The above description is only an alternative embodiment of the embodiments of the present application, but the scope of the embodiments of the present application is not limited thereto, and any person skilled in the art can easily conceive of changes or substitutions within the technical scope of the embodiments of the present application, and all the changes or substitutions should be covered by the scope of the embodiments of the present application.
Claims (10)
1. A grid artifact suppression method, comprising:
acquiring the radiographic image, wherein the radiographic image is an image of X-rays after the X-rays pass through the grid and the region of the object to be inspected;
performing Fourier transform on the radiographic image to obtain a transformed radiographic image;
filtering and suppressing the grid artifacts in the transformed radiographic image by using an adaptive filter to obtain a suppressed radiographic image;
and carrying out inverse Fourier transform on the suppressed radiographic image to obtain a result image.
2. The method of claim 1, after said acquiring the radiograph and before said fourier transforming the radiograph, further comprising:
judging whether grid artifacts exist in the radiographic image or not;
if yes, performing Fourier transform on the radiographic image.
3. The method of claim 1, further comprising, prior to said filtering grid artifacts in the transformed radiographs using an adaptive filter:
acquiring a region of interest in which grid artifacts exist in the radiographic image;
carrying out Fourier transform on the region of interest to obtain a frequency domain image of the region of interest;
calculating the frequency domain characteristics of the grid artifact in the frequency domain image;
and constructing an adaptive filter according to the frequency domain features in the frequency domain image.
4. The method of claim 3, wherein acquiring a region of interest in the radiograph in which grid artifacts are present comprises:
selecting the region of interest from the radiographic image according to preset rules, wherein the preset rules comprise: the cut line identified by the shutter software and/or the coordinate parameters sent by the shutter hardware.
5. The method of claim 3, further comprising, prior to said Fourier transforming said region of interest:
and performing boundary clipping, mirror image continuation, spatial filtering, spatial stacking, multi-scale decomposition and/or windowing on the region of interest.
6. The method of claim 3, wherein the frequency domain features comprise: frequency range, curve shape and image direction; the calculating the frequency domain features of the grid artifacts in the frequency domain image comprises:
counting a frequency range of the grid artifact in the frequency domain image;
performing curve fitting on the frequency domain image to obtain the curve shape of the grid artifact in the frequency domain image;
and determining the image direction of the grid relative to the detector according to the position and the characteristics of the grid artifact in the frequency domain image.
7. The method of any of claims 1-6, further comprising, prior to said Fourier transforming said radiograph:
and performing boundary clipping, mirror image continuation, spatial filtering, spatial superposition, multi-scale decomposition and/or windowing on the radiographic image.
8. A grid artifact suppression device, comprising:
a radiographic image acquisition module for acquiring the radiographic image, which is an image of the X-ray after passing through the grid and the region of the object to be inspected;
the radiographic image transformation module is used for carrying out Fourier transformation on the radiographic image to obtain a transformed radiographic image;
a radiographic filtering module, configured to perform filtering suppression on the grid artifacts in the transformed radiographic image by using an adaptive filter, so as to obtain a suppressed radiographic image;
and the inverse transformation module of the radiographic image is used for carrying out inverse Fourier transform on the suppressed radiographic image to obtain a result image.
9. An electronic device, comprising: a processor and a memory, the memory storing machine-readable instructions executable by the processor, the machine-readable instructions, when executed by the processor, performing the method of any of claims 1 to 7.
10. A computer-readable storage medium, having stored thereon a computer program which, when executed by a processor, performs the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111105620.8A CN113808232A (en) | 2021-09-22 | 2021-09-22 | Grid artifact suppression method and device, electronic device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111105620.8A CN113808232A (en) | 2021-09-22 | 2021-09-22 | Grid artifact suppression method and device, electronic device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113808232A true CN113808232A (en) | 2021-12-17 |
Family
ID=78939944
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111105620.8A Pending CN113808232A (en) | 2021-09-22 | 2021-09-22 | Grid artifact suppression method and device, electronic device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113808232A (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222318A (en) * | 2010-04-16 | 2011-10-19 | 深圳迈瑞生物医疗电子股份有限公司 | Grid artifact detection and suppression method and device |
| CN104077743A (en) * | 2013-03-25 | 2014-10-01 | 深圳市蓝韵实业有限公司 | Suppression method and device for grid artifact in X-ray image |
| CN104574288A (en) * | 2013-10-22 | 2015-04-29 | 辽宁开普医疗系统有限公司 | Self-adaptive grid artifact suppression method and device |
| WO2017045618A1 (en) * | 2015-09-15 | 2017-03-23 | Shanghai United Imaging Healthcare Co., Ltd. | Image reconstruction system and method |
| GB201704783D0 (en) * | 2016-04-07 | 2017-05-10 | Gen Electric | Method and system for reducing ringing artifacts in x-ray image |
| CN106651814A (en) * | 2017-01-12 | 2017-05-10 | 沈阳东软医疗系统有限公司 | Image processing method and apparatus |
| CN107203983A (en) * | 2016-03-17 | 2017-09-26 | 通用电气公司 | Method and system for reducing the grid line artifact in radioscopic image |
| CN111110259A (en) * | 2018-10-30 | 2020-05-08 | 辽宁开普医疗系统有限公司 | Method for identifying installation state of grid and inhibiting grid shadow |
| CN111383186A (en) * | 2018-12-29 | 2020-07-07 | Tcl集团股份有限公司 | Image processing method and device and terminal equipment |
-
2021
- 2021-09-22 CN CN202111105620.8A patent/CN113808232A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222318A (en) * | 2010-04-16 | 2011-10-19 | 深圳迈瑞生物医疗电子股份有限公司 | Grid artifact detection and suppression method and device |
| CN104077743A (en) * | 2013-03-25 | 2014-10-01 | 深圳市蓝韵实业有限公司 | Suppression method and device for grid artifact in X-ray image |
| CN104574288A (en) * | 2013-10-22 | 2015-04-29 | 辽宁开普医疗系统有限公司 | Self-adaptive grid artifact suppression method and device |
| WO2017045618A1 (en) * | 2015-09-15 | 2017-03-23 | Shanghai United Imaging Healthcare Co., Ltd. | Image reconstruction system and method |
| CN107203983A (en) * | 2016-03-17 | 2017-09-26 | 通用电气公司 | Method and system for reducing the grid line artifact in radioscopic image |
| GB201704783D0 (en) * | 2016-04-07 | 2017-05-10 | Gen Electric | Method and system for reducing ringing artifacts in x-ray image |
| CN106651814A (en) * | 2017-01-12 | 2017-05-10 | 沈阳东软医疗系统有限公司 | Image processing method and apparatus |
| CN111110259A (en) * | 2018-10-30 | 2020-05-08 | 辽宁开普医疗系统有限公司 | Method for identifying installation state of grid and inhibiting grid shadow |
| CN111383186A (en) * | 2018-12-29 | 2020-07-07 | Tcl集团股份有限公司 | Image processing method and device and terminal equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10304217B2 (en) | Method and system for generating image using filtered backprojection with noise weighting and or prior in | |
| US7831097B2 (en) | System and method for image reconstruction | |
| RU2737018C1 (en) | Image processing device, image processing method and data medium | |
| US9953441B2 (en) | Processing metal artifact in CT image | |
| KR20130038794A (en) | Method of noise reduction in digital x-ray frames series | |
| US7929746B2 (en) | System and method for processing imaging data | |
| US9418417B2 (en) | Image processing apparatus, tomography apparatus, image processing method, and storage medium | |
| CN110163857B (en) | Image background area detection method and device, storage medium and X-ray system | |
| WO2005110232A1 (en) | Image processing device and method thereof | |
| CN109146803B (en) | SAR image radiation resolution improving method and device based on multi-angle image | |
| CN111353958A (en) | Image processing method, device and system | |
| US8718348B2 (en) | Grid suppression in imaging | |
| CN106033598B (en) | Digital grid imaging method for eliminating scattered radiation influence | |
| Toledo et al. | Image denoising based on genetic algorithm | |
| US10755389B2 (en) | Image processing apparatus, image processing method, and medium | |
| CN113808232A (en) | Grid artifact suppression method and device, electronic device and storage medium | |
| CN111080534B (en) | Image filtering method and device and electronic equipment | |
| US10891764B2 (en) | Image processing apparatus, image processing method, and image processing program | |
| JP2018503436A (en) | Anti-correlation noise filter | |
| US20190221010A1 (en) | Image processing apparatus, image processing method, and image processing program | |
| US20220172461A1 (en) | Machine learning apparatus, machine learning method, and non-transitory computer-readable storage medium | |
| Boby et al. | Medical Image Denoising Techniques against Hazardous Noises: An IQA Metrics Based Comparative Analysis | |
| KR102438068B1 (en) | Method and apparatus for enhancing the image quality using non-local total variation denoiser based on mutual information | |
| WO2021157403A1 (en) | Image processing device, method, and program | |
| CN109813259B (en) | High dynamic X-ray imaging method, storage medium and apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |