X-ray scattering correction method and device, cone beam computed tomography system and product

Abstract

The embodiment of the invention discloses an X-ray scattering correction method and device, a cone beam computed tomography system and a product, and the method comprises the steps: obtaining first projection information of a first detection region of a first detector and second projection information of a second detection region; wherein the second detection area is a shielding area corresponding to the collimators in the first detector under the condition that the collimators are arranged on the two axial sides of the X-ray machine; the first detection area is an area except the second detection area in the first detector; generating a first scatter correction template with the same dimension as the first detection area based on the second projection information; and performing scattering correction on the first projection information according to the first scattering correction template to obtain third projection information. According to the technical scheme of the embodiment of the invention, the problem of poor scattering correction effect is solved, and the quality and uniformity of the cone beam computed tomography medical image can be improved.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of medical imaging technology, and in particular to an X-ray scatter correction method, device, cone-beam computed tomography system, and product. Background Art

[0002] In cone-beam computed tomography (CBCT) systems with offset-mounted flat panels, commonly used scatter correction methods are ineffective due to the limited projection area. The larger the scanned object, the more scattering there is, which significantly impacts the uniformity of the scanned image. Summary of the Invention

[0003] Embodiments of the present invention provide an X-ray scatter correction method, apparatus, cone-beam computed tomography system, and product, which can improve the correction effect of scatter correction of projection information, thereby improving the quality and uniformity of cone-beam computed tomography medical images.

[0004] In a first aspect, an embodiment of the present invention provides an X-ray scatter correction method, the method comprising:

[0005] Obtaining first projection information of a first detection area of ​​the first detector and second projection information of a second detection area; wherein the second detection area is a shielded area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area;

[0006] generating a first scatter correction template having the same dimensions as the first detection area based on the second projection information;

[0007] Scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information.

[0008] In a second aspect, an embodiment of the present invention provides an X-ray scatter correction device, the device comprising:

[0009] A projection information acquisition module is configured to acquire first projection information of a first detection area of ​​the first detector and second projection information of a second detection area; wherein the second detection area is an occluded area of ​​the first detector corresponding to the collimator when collimators are provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area;

[0010] a correction template generating module, configured to generate a first scatter correction template having the same dimension as the first detection area based on the second projection information;

[0011] The scatter correction module is configured to perform scatter correction on the first projection information according to the first scatter correction template to obtain third projection information.

[0012] In a third aspect, an embodiment of the present invention provides a cone-beam computed tomography system, the system comprising:

[0013] X-ray machines, collimators, flat panel detectors, and computer equipment;

[0014] Wherein, the collimators are symmetrically arranged on both sides of the axis of the X-ray machine;

[0015] Computer equipment includes:

[0016] one or more processors;

[0017] a memory for storing one or more programs;

[0018] When the one or more programs are executed by one or more processors, the one or more processors implement the X-ray scatter correction method provided by any embodiment of the present invention.

[0019] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the X-ray scatter correction method provided by any embodiment of the present invention.

[0020] In a fifth aspect, an embodiment of the present invention further provides a computer program product, comprising a computer program, which, when executed by a processor, implements the X-ray scatter correction method provided by any embodiment of the present invention.

[0021] The embodiments of the above invention have the following advantages or beneficial effects:

[0022] In an embodiment of the present invention, first projection information of a first detection area of ​​a first detector and second projection information of a second detection area are obtained; wherein the second detection area is the occluded area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; the first detection area is the area of ​​the first detector excluding the second detection area; based on the second projection information, a first scatter correction template with the same dimensions as the first detection area is generated; and scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information. The technical solution of the embodiment of the present invention solves the problem of poor scatter correction effect and can improve the effect of scatter correction on projection information, thereby improving the quality and uniformity of cone-beam computed tomography medical images. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 is a flow chart of an X-ray scatter correction method provided by an embodiment of the present invention;

[0024] Figure 2 This is a flow chart of another X-ray scatter correction method provided by an embodiment of the present invention;

[0025] Figure 3A This is a flow chart of another X-ray scatter correction method provided by an embodiment of the present invention;

[0026] Figure 3B 1 is a schematic diagram of an image reconstruction result without scatter correction provided by an embodiment of the present invention;

[0027] Figure 3C 1 is a schematic diagram of a scatter-corrected image reconstruction result provided by an embodiment of the present invention;

[0028] Figure 3D This is a schematic diagram comparing horizontal reconstructed image analysis diagrams provided by an embodiment of the present invention;

[0029] Figure 3E This is a schematic diagram comparing vertical reconstructed image analysis diagrams provided by an embodiment of the present invention;

[0030] Figure 4 1 is a schematic structural diagram of an X-ray scatter correction device provided by an embodiment of the present invention;

[0031] Figure 5A 1 is a schematic structural diagram of a cone-beam computed tomography system provided by an embodiment of the present invention;

[0032] Figure 5B This is a front view of a detector in a cone-beam computed tomography system provided by an embodiment of the present invention;

[0033] Figure 6 It is a structural diagram of a computer device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0034] The present invention will be further described in detail below with reference to the accompanying drawings and examples. It will be understood that the specific embodiments described herein are intended only to illustrate the present invention and are not intended to limit the present invention. It should also be noted that, for ease of description, the accompanying drawings only illustrate portions relevant to the present invention, not all structures.

[0035] Figure 1 This is a flowchart of an X-ray scatter correction method provided in an embodiment of the present invention. This embodiment is applicable to cone-beam computed tomography imaging scenarios. This method can be performed by an X-ray scatter correction device, which can be implemented using software and / or hardware and integrated into a computer device with application development capabilities, such as a medical image acquisition and processing device.

[0036] like Figure 1 As shown, the X-ray scatter correction method of this embodiment includes the following steps:

[0037] S110 : Acquire first projection information of a first detection area and second projection information of a second detection area of ​​a first detector.

[0038] The first detector may be a detector that acquires X-ray projection information during medical image acquisition, and may be a flat panel detector, an arc detector, or other different types of detectors.

[0039] Projection information is the modulation of the incident X-ray signal by the internal structure of the object being scanned during X-ray imaging. It reflects the object's properties, such as density, thickness, and chemical composition. As X-rays pass through an object, different tissues (such as bone and soft tissue) attenuate them to varying degrees. The intensity differences recorded by the detector are the projection information.

[0040] The first projection information of the first detection area and the second projection information of the second detection area are projection information collected from different areas of the first detector. The second detection area is the area of ​​the first detector blocked by the collimator when collimators are installed on both sides of the X-ray machine axis. The first detection area is the area of ​​the first detector excluding the second detection area.

[0041] For example, the positional relationship between the first detection area and the second detection area can be referred to Figure 5A The detectors in a cone-beam computed tomography system are shown. Figure 5A In the X-ray machine, two collimators are placed symmetrically on the axis of the A-ray machine at the exit of the X-ray machine, corresponding to Figure 5A The collimator changes the X-ray exit angle, which in turn limits the amount of projection information that passes through the scanned object and reaches the detector, forming two symmetrical projection shielding areas on the detector. Figure 5A The detector in the middle is a side view of the detector. The upper and lower projected occlusion areas correspond to the second detection area, and the middle non-projected occlusion area corresponds to the first detection area.

[0042] S120 : Generate a first scatter correction template having the same dimension as the first detection area based on the second projection information.

[0043] exist Figure 5A In the case of the collimator setting shown, it can be considered that the Del projection information detected in the second projection area is all scattered information, which can reflect the scattering distribution law.

[0044] Specifically, when generating a first scatter correction template with the same dimensions as the first detection area based on the second projection information, the second projection data can be preprocessed, such as by dark current correction, gain correction, and / or bad pixel repair. Then, if the second projection information and the first detection area have different detector physical dimensions, the second projection information needs to be interpolated or resampled to unify the pixel dimensions, thereby generating projection information with the same dimensions as the first detection area as the first scatter correction template.

[0045] Furthermore, template calibration can be performed on the first scatter correction template. For example, background noise can be subtracted by performing flat-field correction to remove residual noise in uniform areas (e.g., calculating the template mean and normalizing it); and range truncation can be performed to limit scatter values ​​to a reasonable range (e.g., 0-1) to avoid negative values ​​or overestimation.

[0046] S130 : Perform scatter correction on the first projection information according to the first scatter correction template to obtain third projection information.

[0047] The first scatter correction template has the same dimensions as the first detection area, so each pixel in the first scatter correction template can correspond one-to-one with each pixel in the first projection information. Scatter correction is completed by subtracting the information value of the pixel corresponding to the first scatter correction template from the information value of each pixel in the first projection information.

[0048] The technical solution of the present invention obtains first projection information from a first detection area of ​​a first detector and second projection information from a second detection area; wherein the second detection area is the occluded area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; the first detection area is the area of ​​the first detector excluding the second detection area; based on the second projection information, a first scatter correction template with the same dimensions as the first detection area is generated; and scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information. The technical solution of the embodiment of the present invention solves the problem of poor scatter correction effect, can optimize the effect of scatter correction, and improve the quality and uniformity of cone-beam computed tomography medical images.

[0049] Figure 2 This is a flowchart of an X-ray scatter correction method provided in an embodiment of the present invention. This embodiment shares the same inventive concept as the data processing method described in the previous embodiment and further describes the process of obtaining a first scatter correction template. This method can be performed by an X-ray scatter correction device, which can be implemented using software and / or hardware and integrated into a computer device with application development capabilities, such as a medical image acquisition and processing device.

[0050] like Figure 2As shown, the X-ray scatter correction method of this embodiment includes the following steps:

[0051] S210 : Acquire first projection information of a first detection area and second projection information of a second detection area of ​​a first detector.

[0052] Among them, the second detection area is the shielding area corresponding to the collimator in the first detector when collimators are set on both sides of the X-ray machine axis; the first detection area is the area in the first detector except the second detection area.

[0053] S220 , calculating an average value of the first sub-region projection information and the second sub-region projection information in the second projection information to obtain average projection information of the second detection area.

[0054] In the embodiment, the first sub-region projection information and the second sub-region projection information respectively correspond to the projection information of the shielding region corresponding to the collimator. For example, the first sub-region projection information and the second sub-region projection information correspond to the projection information of the upper and lower sub-regions in the second detection region.

[0055] For example, Figure 5A The two occluded areas above and below the central axis can be averaged, and then filtered. This can be considered the horizontal distribution of scattered information within the projection. The projection information of the first and second sub-areas can be mapped pixel by pixel, and the projection information values ​​of the corresponding pixels can be averaged to obtain the average projection information.

[0056] In another optional embodiment, since the number of columns corresponding to the first sub-area projection information and the second sub-area projection information is the same, the average information value of the first sub-area projection information and the second sub-area projection information in the column dimension can be calculated to obtain a one-dimensional array in the horizontal direction as the average projection information. That is, the projection information of the pixels in the columns corresponding to the first sub-area projection information and the second sub-area projection information is averaged, with each column corresponding to a projection information mean value. This can result in a one-dimensional array in the horizontal direction as the average projection information, thereby obtaining the average projection information. In some scenarios, the average information value of the first sub-area projection information and the second sub-area projection information in the row dimension can also be calculated to obtain a one-dimensional array in the vertical direction as the average projection information.

[0057] After obtaining the average projection information, the average projection information can also be filtered to remove random noise in the projection data; perform edge enhancement, highlight the object boundary, and improve contrast; reduce artifacts caused by factors such as scattering and motion, make the data smooth, and optimize data continuity to facilitate subsequent reconstruction.

[0058] S230 : Copy the average projection information in the row or column dimension according to the dimensional information of the first detection area to obtain a first scatter correction template.

[0059] The average projection information of the one-dimensional array in the horizontal direction may be copied according to the corresponding number of rows in the dimensional information of the first detection area, thereby obtaining a first scatter correction template with the same dimension as the first detection area.

[0060] The average projection information of the one-dimensional array in the numerical direction can be copied according to the number of columns corresponding to the dimensional information of the first detection area. If the number of rows in the copied result does not match the dimensionality of the first detection area, a first scatter correction template with the same dimensionality as the first detection area can be obtained by copying only the row dimensions or by interpolating the values.

[0061] S240: Perform scatter correction on the first projection information according to the first scatter correction template to obtain third projection information.

[0062] The third projection information is the projection information that has undergone projection correction.

[0063] The technical solution of the present invention obtains first projection information of a first detection area of ​​a first detector and second projection information of a second detection area; wherein the second detection area is the occlusion area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; the first detection area is the area of ​​the first detector excluding the second detection area; the first sub-area projection information and the second sub-area projection information in the second projection information are averaged to obtain the average projection information of the second detection area; the average projection information is copied in the row or column dimension according to the dimensional information of the first detection area to obtain a first scatter correction template; the first projection information is scatter corrected according to the first scatter correction template to obtain third projection information. The technical solution of the embodiment of the present invention solves the problem of poor scatter correction effect, can optimize the effect of scatter correction, and improve the quality and uniformity of cone-beam computed tomography medical images.

[0064] Figure 3A This is a flowchart of an X-ray scatter correction method provided by an embodiment of the present invention. This embodiment, which shares the same inventive concept as the data processing method described in the previous embodiment, further describes the process. This method can be performed by an X-ray scatter correction device, which can be implemented using software and / or hardware and integrated into a computer device with application development capabilities, such as a medical image acquisition and processing device.

[0065] like Figure 3A As shown, the X-ray scatter correction method of this embodiment includes the following steps:

[0066] S310 : Acquire first projection information of a first detection area and second projection information of a second detection area of ​​a first detector.

[0067] Among them, the second detection area is the shielding area corresponding to the collimator in the first detector when collimators are set on both sides of the X-ray machine axis; the first detection area is the area in the first detector except the second detection area.

[0068] S320 : Generate a first scatter correction template having the same dimension as the first detection area based on the second projection information.

[0069] S330, determining the pixel points corresponding to the first projection information and the first scattering template; and for each pixel point, subtracting the projection information in the first scattering template from the first projection information corresponding to the pixel point to obtain third projection information, thereby completing scattering correction of the first projection information.

[0070] S340: Perform image reconstruction according to the third projection information to obtain a first image.

[0071] A filtered back projection algorithm or a neural network-based image reconstruction algorithm may be used to perform image reconstruction to obtain the first image.

[0072] In an example of image reconstruction, the image reconstruction result without scatter correction is as follows Figure 3B As shown, the image reconstruction result after scatter correction by the scatter correction method in this embodiment is as follows: Figure 3C shown.

[0073] against Figure 3B and Figure 3C For analysis, the corresponding reconstructed image analysis diagram can be referred to Figure 3D and Figure 3E .in, Figure 3D This is a horizontal cross-sectional view of the middle 10 layers. Figure 3E The following is a profile of the middle 10 vertical layers. The horizontal axis of the profile represents pixel position, and the vertical axis represents the pixel value corresponding to each pixel. The blue curve represents the profile data of the original reconstructed image (the image reconstructed without the X-ray scatter correction method of the present embodiment), and the red curve represents the profile data of the reconstructed image after scatter correction using the X-ray scatter correction method of the present embodiment. The performance of the reconstructed image corrected by the X-ray scatter correction method provided in this embodiment is significantly better than that of the conventional reconstructed image without the X-ray scatter correction method provided in this embodiment.

[0074] The technical solution of the present invention obtains first projection information of a first detection area of ​​a first detector and second projection information of a second detection area; wherein the second detection area is the occlusion area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; the first detection area is the area of ​​the first detector excluding the second detection area; based on the second projection information, a first scatter correction template with the same dimension as the first detection area is generated; the pixel points corresponding to the first projection information and the first scatter template are determined; and for each pixel point, the first projection information corresponding to the pixel point is subtracted from the projection information in the corresponding first scatter template to obtain third projection information, thereby completing the scatter correction of the first projection information; and image reconstruction is performed based on the third projection information to obtain the first image. The technical solution of the embodiment of the present invention solves the problem of poor scatter correction effect, can optimize the effect of scatter correction, and improve the quality and uniformity of cone-beam computed tomography medical images.

[0075] Figure 4 This is a schematic diagram of the structure of an X-ray scatter correction device provided in an embodiment of the present invention. This embodiment is applicable to X-ray scatter correction scenarios in cone-beam computed tomography imaging. The X-ray scatter correction device can be implemented using software and / or hardware and integrated into a computer terminal device with application development capabilities.

[0076] like Figure 4 As shown, the X-ray scatter correction device includes: a projection information acquisition module 410 , a correction template generation module 420 and a scatter correction module 430 .

[0077] Among them, the projection information acquisition module 410 is used to obtain the first projection information of the first detection area of ​​the first detector and the second projection information of the second detection area; wherein, the second detection area is the occlusion area corresponding to the collimator in the first detector when a collimator is arranged on both sides of the X-ray machine axis; the first detection area is the area of ​​the first detector other than the second detection area; the correction template generation module 420 is used to generate a first scattering correction template with the same dimension as the first detection area based on the second projection information; the scattering correction module 430 is used to perform scattering correction on the first projection information according to the first scattering correction template to obtain third projection information.

[0078] The technical solution of this embodiment obtains first projection information from a first detection area of ​​a first detector and second projection information from a second detection area; wherein the second detection area is the occluded area of ​​the first detector corresponding to the collimator when collimators are provided on both sides of the X-ray machine axis; and the first detection area is the area of ​​the first detector excluding the second detection area; based on the second projection information, a first scatter correction template having the same dimensions as the first detection area is generated; and scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information. The technical solution of this embodiment of the present invention solves the problem of poor scatter correction effect and can improve the effect of scatter correction on projection information, thereby improving the quality and uniformity of cone-beam computed tomography medical images.

[0079] In an optional implementation, the correction template generation module 420 is specifically configured to:

[0080] Calculating an average value of the first sub-region projection information and the second sub-region projection information in the second projection information to obtain average projection information of the second detection area;

[0081] Copying the average projection information in the row or column dimension according to the dimensional information of the first detection area to obtain a first scatter correction template;

[0082] The first sub-region projection information and the second sub-region projection information respectively correspond to projection information of the shielding region corresponding to the collimator.

[0083] In an optional implementation, the correction template generation module 420 is specifically configured to:

[0084] The average information value of the first sub-region projection information and the second sub-region projection information in the column dimension is calculated to obtain a one-dimensional array in the horizontal direction as the average projection information.

[0085] In an optional implementation, the correction template generation module 420 may also be used to:

[0086] After obtaining the average projection information, filtering is performed on the average projection information.

[0087] In an optional embodiment, the X-ray scatter correction device further includes an image reconstruction module, configured to:

[0088] Image reconstruction is performed according to the third projection information to obtain a first image.

[0089] In an optional embodiment, the scatter correction module 430 is specifically configured to:

[0090] Determining a pixel point corresponding to the first projection information and the first scattering template;

[0091] For each pixel point, the projection information in the corresponding first scattering template is subtracted from the first projection information corresponding to the pixel point to complete the scattering correction of the first projection information.

[0092] The X-ray scatter correction device provided in the embodiment of the present invention can execute the X-ray scatter correction method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0093] Figure 5A This is a structural diagram of a cone-beam computed tomography system provided by an embodiment of the present invention. This embodiment is applicable to scenarios of cone-beam-based medical imaging, especially scanning imaging scenarios with a large field of view.

[0094] like Figure 5A As shown, the cone-beam computed tomography system provided by the embodiment of the present invention specifically includes:

[0095] X-ray machines, collimators, flat panel detectors and computer equipment.

[0096] The collimators are symmetrically arranged on both sides of the X-ray machine axis; the collimators are arranged so that the ratio of the area of ​​the regions blocked at both ends of the flat panel detector to the overall area of ​​the flat panel detector is a first preset ratio value.

[0097] For example, a collimator can be added to the axial direction of the X-ray machine's light outlet to adjust the light output angle so that approximately 10% of the entire flat-panel detector is blocked at both ends of the axial direction.

[0098] For further information, please refer to Figure 5B Front view of the flat-panel detector, with the upper and lower ends blocked.

[0099] The computer device includes: one or more processors;

[0100] a memory for storing one or more programs;

[0101] When one or more programs are executed by one or more processors, the one or more processors implement the X-ray scatter correction method of any of the above embodiments.

[0102] Figure 6 A schematic structural diagram of a computer device provided in an embodiment of the present invention. Figure 6 A block diagram of an exemplary computer device 12 suitable for use in implementing embodiments of the present invention is shown. Figure 6 The computer device 12 shown is only an example and should not limit the functionality and scope of use of the embodiments of the present invention. The computer device 12 can be any terminal device with computing capabilities, such as an intelligent controller, a server, a mobile phone, or other terminal devices.

[0103] like Figure 6 As shown, computer device 12 is implemented as a general-purpose computing device. Components of computer device 12 may include, but are not limited to, one or more processors or processing units 16, system memory 28, and a bus 18 that connects various system components (including system memory 28 and processing unit 16).

[0104] Bus 18 represents one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus architectures. Examples of these architectures include, but are not limited to, an Industry Standard Architecture (ISA) bus, a Micro Channel Architecture (MAC) bus, an Enhanced ISA bus, a Video Electronics Standards Association (VESA) local bus, and a Peripheral Component Interconnect (PCI) bus.

[0105] The computer device 12 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the computer device 12, including volatile and non-volatile media, removable and non-removable media.

[0106] System memory 28 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 30 and / or cache 32. Computer device 12 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, storage system 34 may be configured to read and write non-removable, non-volatile magnetic media ( Figure 6 Not shown, often called a "hard drive"). Although Figure 6 Not shown, a magnetic disk drive for reading and writing to a removable non-volatile magnetic disk (e.g., a "floppy disk"), and an optical disk drive for reading and writing to a removable non-volatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical media) may be provided. In these cases, each drive may be connected to bus 18 via one or more data media interfaces. System memory 28 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of various embodiments of the present invention.

[0107] A program / utility 40 having a set (at least one) of program modules 42 may be stored, for example, in system memory 28. Such program modules 42 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data, each of which, or some combination thereof, may include an implementation of a network environment. Program modules 42 generally perform the functions and / or methods of the embodiments described herein.

[0108] The computer device 12 may also communicate with one or more external devices 14 (e.g., a keyboard, a pointing device, a display 24, etc.), one or more devices that enable a user to interact with the computer device 12, and / or any device that enables the computer device 12 to communicate with one or more other computing devices (e.g., a network card, a modem, etc.). Such communication may be performed via an input / output (I / O) interface 22. Furthermore, the computer device 12 may also communicate with one or more networks (e.g., a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) via a network adapter 20. As shown, the network adapter 20 communicates with the other modules of the computer device 12 via the bus 18. It should be understood that although Figure 6 Not shown, other hardware and / or software modules may be used in conjunction with the computer device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RFID systems, tape drives, and data backup storage systems.

[0109] The processing unit 16 executes various functional applications and data processing by running programs stored in the system memory 28, such as implementing the X-ray scatter correction method provided in the embodiment of the present invention, which includes:

[0110] Obtaining first projection information of a first detection area of ​​the first detector and second projection information of a second detection area; wherein the second detection area is a shielded area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area;

[0111] generating a first scatter correction template having the same dimensions as the first detection area based on the second projection information;

[0112] Scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information.

[0113] An embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the X-ray scatter correction method provided in any embodiment of the present invention is implemented. The method includes:

[0114] Obtaining first projection information of a first detection area of ​​the first detector and second projection information of a second detection area; wherein the second detection area is a shielded area corresponding to the collimator in the first detector when collimators are provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area;

[0115] generating a first scatter correction template having the same dimensions as the first detection area based on the second projection information;

[0116] Scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information.

[0117] The computer storage medium of the embodiment of the present invention can adopt any combination of one or more computer-readable media. The computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. The computer-readable storage medium can be, for example, but not limited to: an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or component, or any combination of the above. More specific examples (non-exhaustive list) of computer-readable storage media include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this document, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, device or device.

[0118] A computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, which carries computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0119] Program code embodied on a computer-readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

[0120] Computer program code for performing the operations of the present invention can be written in one or more programming languages, or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, Python, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0121] An embodiment of the present invention further provides a computer program product, including a computer program, which, when executed by a processor, implements the X-ray scatter correction method provided in any embodiment of the present application.

[0122] The computer program product, during implementation, may be written in one or more programming languages, or a combination thereof, for performing the operations of the present invention. The programming languages ​​include object-oriented programming languages ​​such as Java, Smalltalk, Python, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0123] Those skilled in the art will appreciate that the modules or steps of the present invention described above can be implemented using a general-purpose computing device. They can be centralized on a single computing device or distributed across a network of multiple computing devices. Alternatively, they can be implemented using program code executable by a computer device, which can then be stored in a storage device and executed by the computing device. Alternatively, they can be fabricated into separate integrated circuit modules, or multiple modules or steps can be fabricated into a single integrated circuit module. Thus, the present invention is not limited to any specific combination of hardware and software.

[0124] Note that the above are only preferred embodiments of the present invention and the technical principles employed. Those skilled in the art will understand that the present invention is not limited to the specific embodiments described herein, and that various obvious changes, readjustments, and substitutions can be made by those skilled in the art without departing from the scope of protection of the present invention. Therefore, although the present invention has been described in detail through the above embodiments, the present invention is not limited to the above embodiments and may include many other equivalent embodiments without departing from the concept of the present invention. The scope of the present invention is determined by the scope of the appended claims.

Claims

1. An X-ray scattering correction method, characterized in that: include: Obtaining first projection information of a first detection area and second projection information of a second detection area of ​​a first detector; wherein the second detection area is a shielded area of ​​the first detector corresponding to a collimator provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area; generating a first scatter correction template having the same dimensions as the first detection area based on the second projection information; Scatter correction is performed on the first projection information according to the first scatter correction template to obtain third projection information.

2. The method according to claim 1, characterized in that Generating a first scatter correction template based on the second projection information includes: Calculating an average value of the first sub-region projection information and the second sub-region projection information in the second projection information to obtain average projection information of the second detection area; Copying the average projection information in a row or column dimension according to the dimensional information of the first detection area to obtain the first scatter correction template; The first sub-region projection information and the second sub-region projection information respectively correspond to projection information of the shielding region corresponding to the collimator.

3. The method according to claim 2, characterized in that The calculating an average value of the first sub-region projection information and the second sub-region projection information in the second projection information to obtain the average projection information of the second detection area includes: The average information value of the first sub-region projection information and the second sub-region projection information in the column dimension is calculated to obtain a one-dimensional array in the horizontal direction as the average projection information.

4. The method according to claim 2, characterized in that After obtaining the average projection information, the method further includes: Perform filtering processing on the average projection information.

5. The method according to claim 1, wherein The method further comprises: Image reconstruction is performed according to the third projection information to obtain a first image.

6. The method according to claim 1, characterized in that The performing scatter correction on the first projection information according to the first scatter correction template includes: Determining a pixel point corresponding to the first projection information and the first scattering template; For each pixel point, the first projection information corresponding to the pixel point is subtracted from the corresponding projection information in the first scattering template to complete the scattering correction of the first projection information.

7. An X-ray scatter correction device, characterized in that: include: a projection information acquisition module, configured to acquire first projection information of a first detection area of ​​a first detector and second projection information of a second detection area; wherein the second detection area is an occlusion area of ​​the first detector corresponding to a collimator provided on both sides of the X-ray machine axis; and the first detection area is an area of ​​the first detector excluding the second detection area; a correction template generating module, configured to generate a first scatter correction template having the same dimension as the first detection area based on the second projection information; The scatter correction module is configured to perform scatter correction on the first projection information according to the first scatter correction template to obtain third projection information.

8. A cone-beam computed tomography system, characterized in that: include: X-ray machines, collimators, flat panel detectors, and computer equipment; Wherein, the collimators are symmetrically arranged on both sides of the axial direction of the X-ray machine; The computer device includes: one or more processors; a memory for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the X-ray scatter correction method according to any one of claims 1 to 6.

9. The system according to claim 8, characterized in that The collimator is arranged at a position such that the ratio of the area of ​​the regions respectively shielded at the two ends of the flat panel detector to the overall area of ​​the flat panel detector is a first preset ratio value.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the computer program implements the X-ray scatter correction method according to any one of claims 1 to 6.

Images