CN106651985A - Method and device for reconstructing CT (Computed Tomography) image - Google Patents
Method and device for reconstructing CT (Computed Tomography) image Download PDFInfo
- Publication number
- CN106651985A CN106651985A CN201611245662.0A CN201611245662A CN106651985A CN 106651985 A CN106651985 A CN 106651985A CN 201611245662 A CN201611245662 A CN 201611245662A CN 106651985 A CN106651985 A CN 106651985A
- Authority
- CN
- China
- Prior art keywords
- image
- phase place
- phase
- back projection
- images
- 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.)
- Granted
Links
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/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
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention provides a method and device for reconstructing a CT (Computed Tomography) image. The method comprises the following steps of: carrying out preprocessing process before back projection on an image set corresponding to a plurality of phases together; and for an image subset corresponding to each phase, respectively carrying out back reflection, wherein the image subset corresponding to each phase is one part of the image set, and the image subsets corresponding to at least two phases have an intersection. According to the method for reconstructing the CT image, which is provided by the invention, a reconstruction speed can be improved.
Description
Technical field
The invention mainly relates to computed tomography (Computed Tomography, CT) equipment, more particularly to it is a kind of
The method for reconstructing and device of CT images.
Background technology
Computed tomography (Computed Tomography, CT) technology is a kind of different material that is based on for ray
With different attenuation properties, testee is irradiated from all directions with radioactivity, measure the transmitted intensity through object, and passed through
Certain algorithm for reconstructing calculates linear attenuation coefficient of the interior of articles each point material for ray, so as to obtain testee
The radiodiagnostic techniques of faultage image.
During CT scan, the motion (such as heartbeat, breathing etc.) of patient body or organ can be to scan image quality
Impact, cause motion artifacts.For with cardiac CT scan, in order to eliminate impact of the heart movement to image, reduce and rebuild
Image in motion artifacts, can using gate curve be weighted for scan data after carry out image reconstruction again, to obtain
Relatively clearly cardiac image.Specifically, electrocardio (ECG) signal is monitored while heart is scanned, and by one and the heart
The related gate curve of electric signal (being weighting function curve) is weighted to scan data.Adding temporary, to heart movement width
The scan data that gathers participates in image reconstruction using larger weight when degree is less, and gather when larger to heart movement amplitude
Scan data participates in image reconstruction using less weight.
Even so, iff the reconstruction image that single phase place is obtained according to electrocardiosignal, this still cannot be guaranteed
Reconstruction image is best.Therefore use behind CT device scans heart or other positions, doctor is desirable to obtain multiple phases
The reconstruction image of position (phase) such that it is able to therefrom select the reconstruction image of that best phase place of image effect.For this purpose, the phase
Hope that the phase place rebuild is enough, so as to select the image of optimum phase.
The method for reconstructing of current CT equipment, rebuilds be only capable of being rebuild for single phase place every time.Therefore, in order to rebuild
Enough the image of leggy, certainly will need to expend long time.
The content of the invention
The technical problem to be solved in the present invention is to provide the method for reconstructing and device of CT images, can quickly carry out multiple
The image reconstruction of phase place.
To solve above-mentioned technical problem, the invention provides a kind of method for reconstructing of CT images, comprises the following steps：To be many
Image collection corresponding to individual phase place carries out in the lump the preprocessing process before back projection；For image corresponding to each phase place
Collection, carries out respectively back projection, and the wherein image subset corresponding to each phase place is a part for the image collection, and at least two
Exist between image subset corresponding to phase place and occur simultaneously.
An embodiment of the invention, it is each in the image subset corresponding to each phase place when carrying out the back projection
Weight of the individual image according to it in corresponding phase carries out back projection.
An embodiment of the invention, it is each in the image subset corresponding to each phase place when carrying out the back projection
In individual image, the calculating of each pixel value is performed by single thread.
An embodiment of the invention, the faultage image of the advance reconstructed quantity of each phase place.
An embodiment of the invention, said method also includes the selected phase place of basis, rebuilds the institute of the selected phase place
There is faultage image.
An embodiment of the invention, the quantity of the phase place is 2-20.
An embodiment of the invention, above-mentioned CT images are the cardiac tomogram image scanned in multiple phase places
The present invention also provides a kind of reconstructing device of CT images, including memory and processor, and the memory storage has weight
Program is built, the computing device reconstruction algorithm is implementing following step：Image collection one corresponding to multiple phase places is gone forward side by side
Preprocessing process before row back projection；For the image subset corresponding to each phase place, back projection, wherein each phase are carried out respectively
Image subset corresponding to position is a part for the image collection, and is existed between the image subset corresponding at least two phase places
Occur simultaneously.
The present invention also proposes a kind of reconstructing device of CT images, including：For by the image collection corresponding to multiple phase places
The module of the preprocessing process before back projection is carried out in the lump；For the image subset corresponding to each phase place, carry out respectively
The module of back projection, the wherein image subset corresponding to each phase place are a parts for the image collection, and at least two phase places
Exist between corresponding image subset and occur simultaneously.
The present invention also proposes a kind of medical image system, the reconstructing device comprising CT images as above.
Compared with prior art, above-mentioned technical proposal of the invention is by the multiple phase places that there is very big overlapping interval
Scan image merges process so that amount of calculation is greatly reduced, so as to improve the speed of reconstruction.
Description of the drawings
Fig. 1 is the schematic diagram of the image processing system shown in some embodiments of the invention.
Fig. 2 is the schematic diagram of the CT imaging systems shown in some embodiments of the invention.
Fig. 3 is that the framework of the computer collocations of the image processing system shown in some embodiments of the invention shows
It is intended to.
Fig. 4 is the reconstruction flow process of CT images.
Fig. 5 is the method for reconstructing flow chart of CT images according to an embodiment of the invention.
Fig. 6 is 4 Phase Build Out schematic diagrames of CT images according to an embodiment of the invention.
Fig. 7 is the schematic diagram of electrocardiosignal and corresponding weighting function curve.
Fig. 8 is that the reconstruction image of the single phase place of cardiac CT image according to an embodiment of the invention chooses schematic diagram.
Specific embodiment
It is that the above objects, features and advantages of the present invention can be become apparent, below in conjunction with tool of the accompanying drawing to the present invention
Body embodiment elaborates.
Many details are elaborated in the following description in order to fully understand the present invention, but the present invention can be with
It is different from alternate manner described here to implement using other, therefore the present invention is not limited by following public specific embodiment
System.
Fig. 1 is the schematic diagram of the image processing system shown in some embodiments of the invention.The image processing system
100 can include 110, image processing system 120 of an imaging system and a network 130.In certain embodiments,
Imaging system 110 can be single mode imaging device, or multi-mode imaging system.In certain embodiments, image processing system
120 can be that the view data to obtaining is processed to obtain image and/or relevant information.
Imaging system 110 can be single imaging system, or the combination of multiple different imaging systems.Imaging system can be with
It is imaged by scanning a target, in certain embodiments, imaging system can be a medical image system.Medical science into
As system can gather the image information of partes corporis humani position.Medical image system can be the C-arm system of X-ray, combined type
Medical image system etc..Imaging system 110 can include one or more scanners.Scanner can be computed tomography
Scanner (CT Scanner).
Image processing system 120 can process the data message of acquisition.In certain embodiments, data message can include
The combination of one or more such as text message, image information, acoustic information.In certain embodiments, image processing system 120 can
With including the combination of one or more in a processor, a process cores, one or more memories etc..For example, image
Processing system 120 can include central processing unit (Central Processing Unit, CPU), special IC
(Application-Specific Integrated Circuit, ASIC), dedicated instruction processor (Application-
Specific Instruction-Set Processor, ASIP), graphic process unit (Graphics Processing Unit,
GPU), physical manipulations processor (Physics Processing Unit, PPU), digital signal processor (Digital
Signal Processor, DSP), field programmable gate array (Field Programmable Gate Array,
FPGA), programmable logic device (Programmable Logic Device, PLD), controller (Controller), microcontroller
Unit (Microcontroller unit), processor (Processor), microprocessor (Microprocessor), ARM process
The combinations of one or more such as device (Advanced RISC Machines).In certain embodiments, image processing system 120 can
To process the image information obtained from imaging system 110.
Network 130 can be single network, or the combination of multiple heterogeneous networks.For example, network 130 is probably a local
Net (local area network (LAN)), wide area network (wide area network (WAN)), common network, private network,
It is proprietary network, public switch telephone network (public switched telephone network (PSTN)), internet, wireless
The combination of one or more in network, virtual network, MAN, telephone network etc..Network 130 can include multiple networks
Access point, for example, super joint, WAP, base station, internet exchange point etc. in interior wired or wireless access point.
By these access points, data source can send data message with access network 130 and by network 130.In certain embodiments,
Network 130 can be used for the communication of image processing system 120, receive the information inside or outside image processing system 120, Xiang Tu
As the inside other parts of processing system 120 or outside transmission information.
It should be noted that above-mentioned image processing system 120 can be actually existed in imaging system 110, or by cloud
Calculating platform completes corresponding function.The cloud computing platform can be including the storage-type cloud platform based on data storage, with process
Calculation type cloud platform based on data and take into account data storage and the comprehensive cloud computing platform for processing.Imaging system 110 is made
Cloud platform can be public cloud, private clound, community cloud or mixed cloud etc..For example, according to actual needs, imaging system 110
Some image informations and/or data message of output, can be calculated and/or be stored by user's cloud platform.Other figures
As information and/or data message, can be calculated and/or be stored by local image processing system 120.
Fig. 2 is the schematic diagram of the CT imaging systems 200 shown in some embodiments of the invention.Imaging system 200 can
Being a specific embodiment of imaging system 110.The imaging system 200 can include a frame 210 and an examination couch
250。
In certain embodiments, frame 210 can include the rotatable portion for having around the rotation of the axis of imaging system 200
220.The space structure of rotatable portion 220 can be the combination of one or more such as cylinder, ellipsoid, cuboid.One
In a little embodiments, rotatable portion 220 can include x-ray source 230, X-ray detector 240 and scanning cavity 270.It is rotatable
Part 220 can be rotated with the axis 260 of imaging system 200 as axle.X-ray source 230 and X-ray detector 240 can be with
Rotated as axle with axis 260 with rotatable portion 220.
When being checked, an object (for example, patient, die body etc.) can be placed on examination couch 250.Examination couch
250 can be pushed in scanning cavity 270 along Z-direction.When being rotated around axis 260, x-ray source 230 and X-ray
Detector 240 can gather the scan data of patient.Scan data can be used to rebuild, for example, image to be corrected, to be corrected
Reference picture picture of image etc..
In certain embodiments, imaging system 200 can carry out helical scanning.In helical scanning, scanned object can
With movable along axle 260, while x-ray source can be rotated around axle 260.X-ray source can produce spiral shell relative to object
Rotation track.
Fig. 3 is the frame of the computer collocations of the image processing system 120 shown in some embodiments of the invention
Structure schematic diagram.Computer 300 can be used the particular system for realizing implementing to be disclosed in the application.Specific system in the present embodiment
Unite and explain a hardware platform comprising user interface using functional block diagram.Computer 300 can implement present image process
The one or more assemblies of system 120, module, unit, subelement.In addition, image processing system 120 can be by computer 300
Realized by its hardware device, software program, firmware and combinations thereof.This computer can be a general purpose
Computer, an or computer for having a specific purpose.Two kinds of computers can be used for realizing specific in the present embodiment
System.For convenience's sake, a computer is only depicted in Fig. 3, but the offer image procossing institute described by the present embodiment
The correlation computer function of needing information can be to be implemented in a distributed fashion, by one group of similar platform, decentralized system
Process load.
As shown in figure 3, computer 300 may include internal communication bus 310, processor (processor) 320 is read-only to deposit
Reservoir (ROM) 330, random access memory (RAM) 340, COM1 350, input output assembly 360, hard disk 370, and use
Family interface 380.Internal communication bus 310 are capable of achieving the data communication of the inter-module of computer 300.Processor 320 can perform journey
Sequence instruction completes any function, component, module, unit, the subelement of the image processing system 120 described in here application.
Processor 320 can be made up of one or more processors.COM1 350 can realize computer 300 and system 100 other
Data communication between part (such as imaging system 110) (such as by network 130).Computer 300 can also include not similar shape
The program storage unit and data storage element of formula, such as hard disk 370, read-only storage (ROM) 330, random access memory
Device (RAM) 340, stores computer disposal and/or the various data files that use of communication, and performed by processor 320 can
The programmed instruction of energy.Input output assembly 360 support computer 300 and the other assemblies of system 100 (such as imaging system 110) it
Between input/output data stream.Computer 300 can also send and receive information and number by COM1 350 from network 130
According to.
It will be appreciated by those skilled in the art that can there is various modifications and improvement in content disclosed herein.For example,
Different system component described above is all realized by hardware device, but is likely to only by the solution party of software
Case is achieved.For example：The installation system on existing server.Additionally, the offer of positional information described herein may
It is to be achieved by a firmware, the combination of firmware/software, the combination of firmware/hardware or the combination of hardware/firmware/software.
The method for reconstructing and device of the CT images of embodiments of the invention description, can be real in the environment shown in Fig. 1-3
Apply.For example, method for reconstructing and device are implemented in image processing system 120.Method for reconstructing and device can be by once rebuilding
The image of multiple phase places is obtained, the purpose of fast run-up picture is reached.Here, CT images can be cardiac CT image, lung CT image
Or abdominal CT images etc..
Fig. 4 is the reconstruction flow process of cardiac CT image.With reference to shown in Fig. 4, once rebuild and usually contained twenty or thirty algorithm
Step, each step will carry out single treatment to data.Whole process of reconstruction is divided into by species：Pretreatment part
410th, build as part 420 and post processing part 430.Pretreatment part include AIR, SliceNormal, ClossTalk ... HU
Deng, build as part 420 include AziRbin, RadialRebin, filtering (Filter) ..., back projection (Back
Projection, BP) etc., post processing part 130 include RingOff, TV ..., Saver etc..
By known reconstruction mode, reconstruction flow process as shown in Figure 5 can individually be carried out to the image of each phase place, this is needed
Take a substantial amount of time.But the analysis found that there is substantial amounts of overlapped data between the image that each phase place is included, this
A little data have all carried out multiple process and have wasted amount of calculation in pretreatment part algorithm and reconstruction some algorithm.
Fig. 5 is the method for reconstructing flow chart of CT images according to an embodiment of the invention.With reference to shown in Fig. 5, the present embodiment
Method comprise the following steps：
In step 501, the image collection corresponding to multiple phase places is carried out in the lump the preprocessing process before back projection；
In step 502, for the image subset corresponding to each phase place, back projection is carried out respectively, wherein each phase place institute
Corresponding image subset is exist to hand between a part for the image collection, and the image subset corresponding at least two phase places
Collection.
This method for reconstructing may be embodied as computer program, in the ROM 330 being stored in shown in Fig. 3 or hard disk 370,
Processor 320 is loaded into and performs this reconstruction algorithm with execution step 501-503.
In step 501, preprocessing process includes steps necessary and optional step before back projection.By taking Fig. 4 as an example, front place
The step of managing part 410, build as before back projection in the part of part 420, can partly or entirely be included in as needed pre- place
During reason.
In step 502, the image collection corresponding to multiple phase places can be divided into multiple images subset according to phase place, each
Image subset one phase place of correspondence.It is appreciated that exist between the image subset of out of phase occuring simultaneously.That is, not homophase
In the image subset of position, parts of images is to overlap.But in preprocessing process, image collection is not divided according to phase place,
Therefore will not be repeatedly processed through vaporization, condensation and absorption.Name the example explanation of a reality.
Fig. 6 is 4 Phase Build Out schematic diagrames of CT images according to an embodiment of the invention.With reference to shown in Fig. 6,4 phase places
Phase point be respectively：20%, 30%, 40%, 50%.By taking 4 phase places in the 1st cycle as an example, sweeping corresponding to them
Tracing is respectively 900,899,899 and 899 as (rebuilding required scan image) collective number.Corresponding to each phase place
The quantity of scan image is relevant with visual field size is rebuild, usually, when rebuilding the visual field and being less, when the pixel of its image is corresponding
Between resolution ratio it is higher；And when rebuilding the visual field and being larger, in order to the pixel for ensureing big field-of-view image edge has enough projection numbers
According to it is required that weighting function width increase at its corresponding phase place, i.e., can use the longer time (more) at the phase place
Data for projection participates in image reconstruction.From fig. 6 it can be seen that 4 phase places correspond to have part to be weight between scan image respectively
Folded, their image collection only has 1700 scan images altogether.
According to the embodiment of Fig. 5, the image collection corresponding to multiple phase places can together be carried out in step 501 back projection
Front preprocessing process.Specifically, the image collection corresponding to multiple phase places can experience together pretreatment part 110 in Fig. 1,
The step of building as before back projection in part 120.The image set corresponding to 4 phase places still by taking Fig. 3 as an example, in the 1st cycle
Close, i.e., 1700 scan images can in the lump experience pretreatment part 110 in Fig. 1, build as the step before back projection in part 120
Suddenly.Likewise, each follow-up cycle can be also carried out similar step, until the scan image of all visual angles (View) is all processed
Finish.
In back projection's step, then the scan image corresponding to each phase place (hereinafter indicated as image subset) can respectively be entered
Row back projection.The weight in each phase place for considering each scan image is needed during back projection.
Still by taking cardiac CT scan as an example, Fig. 7 is the schematic diagram of electrocardiosignal and corresponding weighting function curve, Fig. 7
Top half is the curve of electrocardio (ECG) signal, and its abscissa is time t, and ordinate is heartbeat amplitude, which represent heartbeat width
Degree is over time；The latter half of Fig. 7 is weighting function curve corresponding with electrocardiosignal.
As a rule, the scan data gathered when less to heart movement amplitude adopts larger weight, so that
The scan data that heart movement amplitude is gathered when less more participates in image reconstruction；And gather when larger to heart movement amplitude
Scan data adopt less weight so that the less participation of the scan data gathered when heart movement amplitude is larger
Image reconstruction.Above method can reduce image of the heart movement for reconstruction image, reduce the motion artifacts in image.As schemed
Shown in 7, in the present embodiment, when larger for heart movement amplitude corresponding weight coefficient be 0, for heart movement amplitude compared with
Corresponding weight coefficient is 1 when big.
Fig. 8 is that the reconstruction image of the single phase place of cardiac CT image according to an embodiment of the invention chooses schematic diagram, such as
Shown in Fig. 8, in the scan image at multiple visual angles (view) of each phase place, each scan image has a weight
Gweight, the weight is the weight of the reconstruction image that the scan image is affected on it.The weight of all of scan image is being regarded
Angular direction forms a weighting curve, and the position of reconstruction image is ImagePosition.
When multiple phase places, each phase place has a weighting curve, and in the position of each reconstruction image
There are multiple images (its quantity is the quantity of phase place), can so reduce the amount of calculation of algorithm before back projection.
The formula of back projection is as follows：
Wherein：
T=yjsinθl+xicosθl
V=-xisinθl+yjcosθl (2)
Wherein：Refer to that the pixel of i row on image, in the position in x directions, refers to that the pixel of j rows on image, in the position in y directions, refers to
K-th image in the position in z directions, i, j, k=1,2,3,4 ....
CDirection is the direction (enter bed or go out bed) of sick bed.The detector that R refers to refers to the distance of pivot
It is the anglec of rotation of bulb, l refers to index of the current visual angle in visual sequence.Pitch refers to pitch, z'kRefer to figure
As rebuilding the z directions position for starting, Couchpos is the position of sick bed, and certain point that t is referred on image is projected in visual angle
Point arrives distance of the visual angle centrally along detector direction (channel).Parameter t be subpoint in channel directions to through reconstruction
The distance of that ray at center, v refers to distance of the subpoint along directions of rays to reconstructed center.
P～CPOriginal (Raw) data of parallel beam being input into are referred to, (θ q) weights weight (Aperture to W for aperture
Weighting), calculation is as follows：
Wherein：Q is the factor for controlling this function flatness, and (m refers in image sequence, the figure of angle effects for m width that dq refers to
The index of picture, scope is 1 to picture number) in image (i, j) individual pixel cut into slices in visual angle (slice) direction subpoint away from
With a distance from the center of visual angle.
According to above-mentioned calculation, back projection's formula of single phase place is just revised as：
Wherein：W (P) refers to the corresponding weight in visual angle.
Back projection's formula of multiple phase places is revised as：
Wherein：L refers to l-th phase place, W (P)lRefer to the weight of visual angle l-th phase place of correspondence.
Back projection that can be in execution step 502 according to above-mentioned formula (5).
In calculating process, pixel value can be calculated using GPU (GPU).It is preferred that each phase place institute is right
In each image in the image subset answered, to each pixel a thread is distributed, and initial data is tied to into texture
Internal memory, to lift the reading speed of data.Such advantage is, therefore the picture number that obtains of n phase place is the n of single phase place
Times, and the time spent in almost as the time of single phase place.
Once build in known reconstruction mode as a phase place can only be rebuild.The data processing of so all algorithm nodes
Amount is：Data1 (corresponding data when 20%)+Data2 (corresponding data when 30%)+Data3 (corresponding data when 40%)+
Data4 (corresponding data when 50%).By contrast, in embodiments of the invention the Data1 that there is very big overlapping interval～
Data4 data merge process so that amount of calculation is greatly reduced, so as to improve the speed of reconstruction.
Carrying out quick leggy reconstruction clinically has significant meaning.For cardiac CT image reconstruction, doctor is not
Know image is best at which phase place, thereby tend to be rebuild at multiple phase places.In general phase place
Reconstruction image quantity is 300 width, and typically patient can be carried out at 10 phase places to build picture.So it is directed to a client need weight
Build 3000 (300x10) width images, then the 300 width images for selecting optimum phase.Apparent this process of reconstruction is a large amount of
Amount of calculation is wasted, this is also the reason for reconstruction speed is not high.More seriously, optimum phase may be not included in 10 for rebuilding
At individual phase place.Using embodiments of the invention, more phase points (tens to hundreds of, such as 100) can be once rebuild
Image (phase point is respectively 1% to 100%), but each phase place only rebuilds the image of very limited amount of setting quantity.For example
1 image that each phase place only can be rebuild at coronary artery.So, it is allowed to which doctor is judged most by evaluating all phase place dot images
Good phase point.CT imaging systems carry out complete image reconstruction according to selected phase point, that is, rebuild the institute of the phase point
There is faultage image.The amount of calculation of this mode only has 400 width images (1x100+300x1), while greatling save amount of calculation
Also improve the accuracy of selection.
From another viewpoint, the present invention also proposes a kind of reconstructing device of CT images, including：For multiple phase places institute is right
The image collection answered carries out in the lump the module of the preprocessing process before back projection；For image corresponding to each phase place
Collection, carries out respectively the module of back projection, and the wherein image subset corresponding to each phase place is a part for the image collection, and extremely
Exist between image subset corresponding to few two phase places and occur simultaneously.
Although the present invention is described with reference to current specific embodiment, those of ordinary skill in the art
It should be appreciated that the embodiment of the above is intended merely to the explanation present invention, can also make in the case of without departing from spirit of the invention
Go out various equivalent change or replacement, therefore, as long as the change, change in the spirit of the present invention to above-described embodiment
Type all will fall in the range of following claims.
Claims (10)
1. a kind of method for reconstructing of CT images, comprises the following steps：
Image collection corresponding to multiple phase places is carried out in the lump the preprocessing process before back projection；
For the image subset corresponding to each phase place, back projection, the wherein image subset corresponding to each phase place are carried out respectively
It is exist to occur simultaneously between a part for the image collection, and the image subset corresponding at least two phase places.
2. method according to claim 1, it is characterised in that figure when carrying out the back projection, corresponding to each phase place
As weight of each image in subset according to it in corresponding phase carries out back projection.
3. method according to claim 1, it is characterised in that figure when carrying out the back projection, corresponding to each phase place
In as each image in subset, the calculating of each pixel value is performed by single thread.
4. method according to claim 1, it is characterised in that each phase place only rebuilds the faultage image of predetermined number.
5. method according to claim 4, it is characterised in that also include, according to selected phase place, rebuilding the selected phase place
All faultage images.
6. method according to claim 1, it is characterised in that the quantity of the phase place is 2-20.
7. method according to claim 1, it is characterised in that the CT images are that the heart scanned in multiple phase places breaks
Tomographic image.
8. a kind of reconstructing device of CT images, including memory and processor, the memory storage has reconstruction algorithm, the processor
Perform the reconstruction algorithm to implement following step：
Image collection corresponding to multiple phase places is carried out in the lump the preprocessing process before back projection；
For the image subset corresponding to each phase place, back projection, the wherein image subset corresponding to each phase place are carried out respectively
It is exist to occur simultaneously between a part for the image collection, and the image subset corresponding at least two phase places.
9. a kind of reconstructing device of CT images, including：
For the module that the image collection corresponding to multiple phase places is carried out in the lump into the preprocessing process before back projection；
For the image subset corresponding to each phase place, the module of back projection is carried out respectively, wherein corresponding to each phase place
Image subset be the image collection a part, and between the image subset corresponding at least two phase places exist occur simultaneously.
10. a kind of medical image system, the reconstructing device comprising CT images as claimed in claim 8 or 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611245662.0A CN106651985B (en) | 2016-12-29 | 2016-12-29 | Reconstruction method and device of CT image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611245662.0A CN106651985B (en) | 2016-12-29 | 2016-12-29 | Reconstruction method and device of CT image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106651985A true CN106651985A (en) | 2017-05-10 |
| CN106651985B CN106651985B (en) | 2020-10-16 |
Family
ID=58835766
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611245662.0A Active CN106651985B (en) | 2016-12-29 | 2016-12-29 | Reconstruction method and device of CT image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106651985B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107221014A (en) * | 2017-06-26 | 2017-09-29 | 上海联影医疗科技有限公司 | The dispatching method and device and medical image system of medical image reconstruction tasks |
| CN111443318A (en) * | 2019-01-16 | 2020-07-24 | 上海联影医疗科技有限公司 | Magnetic resonance image processing method, magnetic resonance image processing device, storage medium and magnetic resonance imaging system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1615104A (en) * | 2002-01-10 | 2005-05-11 | 株式会社日立医药 | X-ray CT imaging method and X-ray ct device |
| CN1708256A (en) * | 2002-10-25 | 2005-12-14 | 皇家飞利浦电子股份有限公司 | Four-dimensional helical tomographic scanner |
| EP1861825B1 (en) * | 2005-03-17 | 2008-10-08 | Philips Intellectual Property & Standards GmbH | Method and device for the iterative reconstruction of cardiac images |
| CN104684482A (en) * | 2012-10-01 | 2015-06-03 | 株式会社东芝 | Image processing device and X-ray CT device |
-
2016
- 2016-12-29 CN CN201611245662.0A patent/CN106651985B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1615104A (en) * | 2002-01-10 | 2005-05-11 | 株式会社日立医药 | X-ray CT imaging method and X-ray ct device |
| CN1708256A (en) * | 2002-10-25 | 2005-12-14 | 皇家飞利浦电子股份有限公司 | Four-dimensional helical tomographic scanner |
| EP1861825B1 (en) * | 2005-03-17 | 2008-10-08 | Philips Intellectual Property & Standards GmbH | Method and device for the iterative reconstruction of cardiac images |
| CN104684482A (en) * | 2012-10-01 | 2015-06-03 | 株式会社东芝 | Image processing device and X-ray CT device |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107221014A (en) * | 2017-06-26 | 2017-09-29 | 上海联影医疗科技有限公司 | The dispatching method and device and medical image system of medical image reconstruction tasks |
| CN111443318A (en) * | 2019-01-16 | 2020-07-24 | 上海联影医疗科技有限公司 | Magnetic resonance image processing method, magnetic resonance image processing device, storage medium and magnetic resonance imaging system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106651985B (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105593905B (en) | The partly-adjusting method to regularization parameter is used for for the image quality optimization in complete 3D iteration CT reconstruction | |
| JP5248648B2 (en) | Computer tomography system and method | |
| CN104644200B (en) | The method and apparatus for reducing pseudomorphism in computed tomography images reconstruct | |
| US20110044559A1 (en) | Image artifact reduction | |
| JP6050029B2 (en) | Method and system for reconstruction of tomographic images | |
| EP4105888A1 (en) | Systems and methods for computed tomography image reconstruction | |
| US12205253B2 (en) | Systems and methods for adaptive blending in computed tomography imaging | |
| US20140016847A1 (en) | Multi-phase computed tomography image reconstruction | |
| CN111150419B (en) | Method and device for reconstructing image by spiral CT scanning | |
| Taguchi et al. | Direct cone‐beam cardiac reconstruction algorithm with cardiac banding artifact correction | |
| US20240404166A1 (en) | Methods, systems, and storage mediums for image generation | |
| CN106651985A (en) | Method and device for reconstructing CT (Computed Tomography) image | |
| US8364244B2 (en) | Methods and systems to facilitate reducing banding artifacts in images | |
| US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
| Leng et al. | Motion artifact reduction in fan-beam and cone-beam computed tomography via the fan-beam data consistency condition (FDCC) | |
| CN109620273A (en) | A kind of quick CBCT algorithm for reconstructing calculating short scanning weight in real time | |
| US20250095239A1 (en) | Methods and systems for generating dual-energy images from a single-energy imaging system | |
| US20230013905A1 (en) | Apparatus for computer tomography x-ray data acquired at high relative pitch | |
| WO2017091864A1 (en) | Visualising three-dimensional images of variable objects | |
| CN119670544A (en) | Method for generating low-resolution and high-resolution data pairs and medical imaging system | |
| CN108366772B (en) | Image processing apparatus and method | |
| WO2024008764A1 (en) | Cone beam artifact reduction | |
| JP6777407B2 (en) | Helical CT device, medical image processing method, and medical image processing program | |
| US20060243914A1 (en) | Attenuation map generation from pet scans | |
| JP2022155027A (en) | Image processing device, learning device, radiation imaging system, image processing method, learning method, image processing program, and learning program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 no.2258, Chengbei Road, Jiading District, Shanghai Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |