CN108898582A - Cardiac image method for reconstructing, device and computer equipment - Google Patents
Cardiac image method for reconstructing, device and computer equipment Download PDFInfo
- Publication number
- CN108898582A CN108898582A CN201810597965.1A CN201810597965A CN108898582A CN 108898582 A CN108898582 A CN 108898582A CN 201810597965 A CN201810597965 A CN 201810597965A CN 108898582 A CN108898582 A CN 108898582A
- Authority
- CN
- China
- Prior art keywords
- thoracic cavity
- contour images
- image
- cavity contour
- 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
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 39
- 210000000115 thoracic cavity Anatomy 0.000 claims abstract description 325
- 238000010586 diagram Methods 0.000 claims description 33
- 238000004590 computer program Methods 0.000 claims description 21
- 239000002872 contrast media Substances 0.000 description 41
- 210000000988 bone and bone Anatomy 0.000 description 26
- 230000000007 visual effect Effects 0.000 description 16
- 239000011159 matrix material Substances 0.000 description 9
- 210000000038 chest Anatomy 0.000 description 8
- 239000003795 chemical substances by application Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 239000000835 fiber Substances 0.000 description 6
- 238000003709 image segmentation Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000000717 retained effect Effects 0.000 description 4
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/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/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- 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/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
This application involves a kind of cardiac image method for reconstructing, device and computer equipments.The method includes：Using multiple preview images as input, thoracic cavity contour images are obtained according to preview image described in multiple and maximum intensity projection；The boundary position of the thoracic cavity contour images is calculated according to thoracic cavity contour images；Reconstructed center is calculated according to the boundary position of the thoracic cavity contour images；Leggy reconstruction is carried out according to the reconstructed center and default field range of rebuilding, obtains cardiac image.By first choosing cardiac position, leggy image reconstruction is carried out according to cardiac position, reduces the data volume of input, further reduces the time of data input, improves operation efficiency.
Description
Technical field
This application involves Image Reconstruction Technology fields, more particularly to a kind of cardiac image method for reconstructing, device and meter
Calculate machine equipment.
Background technique
With the development of Medical Devices technology, the technology for carrying out patient part detection using CT scanner is more and more mature.
When using CT scanner, CT scanner obtains the data of leggy, and carries out image reconstruction according to the data of multiple phases, obtains
The image of patient part.In cardiac CT scan, due to the movement of heart, when image reconstruction, needs to select optimal phase data
Carry out image reconstruction.
In current traditional technology, according to the data of multiple phases, image reconstruction is carried out to the data of each phase respectively,
The reconstruction image of all phases is obtained, is compared to obtain optimum phase reconstruction image by the reconstruction image to all phases.
Region is reduced in optimum phase reconstruction image, finds cardiac position image.Due to needing to carry out figure to the data of all phases
As rebuilding, therefore the input of big data quantity can take a substantial amount of time, and the operation of big data quantity also results in operation efficiency
The problems such as low.
Summary of the invention
Based on this, it is necessary to which in view of the above technical problems, smaller reconstruction matrix, smaller reconstruction can be used by providing one kind
The visual field rebuilds cardiac image method for reconstructing, device, computer that speed obtains apparent more accurate multiphase bit image faster
Equipment and storage medium.
A kind of cardiac image method for reconstructing, the method includes：Using multiple preview images as input, according to multiple
Preview image and maximum intensity projection obtain thoracic cavity contour images；The thoracic cavity profile is calculated according to thoracic cavity contour images
The boundary position of image；Reconstructed center is calculated according to the boundary position of the thoracic cavity contour images；According in the reconstruction
The heart and default field range of rebuilding carry out leggy reconstruction, obtain cardiac image.
In one of the embodiments, it is described using multiple preview images as inputting, according to preview image described in multiple with
And maximum intensity projection obtains thoracic cavity contour images, including：Using multiple preview images as input, by multiple preview images
It is split, obtains the first image sequence；It carries out out operation respectively to the first image sequence, obtains the second image sequence；
Maximum intensity projection is carried out to the first image sequence and the second image sequence respectively and obtains thoracic cavity contour images.
It is described in one of the embodiments, that maximal density is carried out to the first image sequence and the second image sequence respectively
Projection obtains thoracic cavity contour images, including：Maximum intensity projection is axially carried out in thoracic cavity to the first image sequence, obtains the
The maximum-density projection image of one image；Maximum intensity projection is axially carried out in thoracic cavity to second image sequence, obtains the
The maximum-density projection image of two images；By the maximal density of the maximum-density projection image of the first image and the second image
Projected image makes the difference, and obtains thoracic cavity contour images.
The boundary position of the thoracic cavity contour images includes in one of the embodiments,：Leftmost border position, rightmost
Boundary position and the top circle position.
The boundary that the thoracic cavity contour images are calculated according to thoracic cavity contour images in one of the embodiments,
Position, including：Thoracic cavity contour images are split, left thoracic cavity contour images and right thoracic cavity contour images are obtained；According to left chest
Leftmost border position is calculated in chamber contour images；Rightmost circle position is calculated according to right thoracic cavity contour images；According to chest
The top circle position is calculated in chamber contour images.
It is described in one of the embodiments, that leftmost border position is calculated according to left thoracic cavity contour images, including：Root
Largest connected domain is calculated according to left thoracic cavity contour images, and chooses the left thoracic cavity contour images in the largest connected domain；Choose institute
The upper parts of images in the left thoracic cavity contour images in largest connected domain where right hat is stated, and leftmost border position is calculated.
It is described in one of the embodiments, that rightmost circle position is calculated according to right thoracic cavity contour images, including：Root
Largest connected domain is calculated according to right thoracic cavity contour images, and chooses the right thoracic cavity contour images in largest connected domain；Choose most Dalian
Upper parts of images in right thoracic cavity contour images in logical domain where right hat, and rightmost circle position is calculated.
It is described in one of the embodiments, that the top circle position is calculated according to thoracic cavity contour images, including：According to
It sets the lowermost position of the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images and the second image sequence
Region to be analyzed is established in the contour images of thoracic cavity；Largest connected domain is calculated according to the region to be analyzed, and is chosen largest connected
Thoracic cavity contour images in domain；According to the thoracic cavity contour images in largest connected domain, the top circle position is calculated.
A kind of cardiac image reconstructing device, described device include：Maximum intensity projection module, for multiple preview images
As input, thoracic cavity contour images are obtained according to preview image described in multiple and maximum intensity projection；Boundary obtains module, uses
In the boundary position that the thoracic cavity contour images are calculated according to thoracic cavity contour images；Center rebuilds module, for according to institute
Reconstructed center is calculated in the boundary position for stating thoracic cavity contour images；Image reconstruction module, for according to reconstructed center and in advance
If rebuilding field range carries out leggy reconstruction, cardiac image is obtained.
A kind of computer equipment, including memory and processor, the memory are stored with computer program, the processing
Device realizes following steps when executing the computer program：
Using multiple preview images as input, thoracic cavity wheel is obtained according to preview image described in multiple and maximum intensity projection
Wide image；
The boundary position of the thoracic cavity contour images is calculated according to thoracic cavity contour images；
Reconstructed center is calculated according to the boundary position of the thoracic cavity contour images；
Leggy reconstruction is carried out according to the reconstructed center and default field range of rebuilding, obtains cardiac image.
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor
Following steps are realized when row：
Using multiple preview images as input, thoracic cavity wheel is obtained according to preview image described in multiple and maximum intensity projection
Wide image；
The boundary position of the thoracic cavity contour images is calculated according to thoracic cavity contour images；
Reconstructed center is calculated according to the boundary position of the thoracic cavity contour images；
Leggy reconstruction is carried out according to the reconstructed center and default field range of rebuilding, obtains cardiac image.
Above-mentioned cardiac image method for reconstructing, device, computer equipment and storage medium, by using multiple preview images as
Input, and thoracic cavity contour images are obtained according to multiple preview images and maximum intensity projection algorithm, further according to thoracic cavity profile diagram
As calculating leftmost border position, rightmost circle position and the top circle position, reconstructed center is obtained.According to reconstructed center and
Default field range of rebuilding carries out leggy reconstruction, obtains cardiac image.By first choosing cardiac position, according to cardiac position into
Row leggy image reconstruction, reduces the data volume of input, further reduces the time of data input, improves operation effect
Rate.
Detailed description of the invention
Fig. 1 is the flow diagram of one embodiment cardiac image rebuilding method；
Fig. 2 is the method flow schematic diagram that thoracic cavity contour images are determined in one embodiment；
Fig. 3 is the method flow schematic diagram that the boundary position of thoracic cavity contour images is calculated in one embodiment；
Fig. 4 is the flow diagram of another embodiment cardiac image rebuilding method；
Fig. 5 is the maximum-density projection image of the image of bone and contrast agent in one embodiment；
Fig. 6 is the maximum-density projection image of the image of contrast agent in one embodiment；
Fig. 7 is thoracic cavity contour images in one embodiment；
Fig. 8 is left thoracic cavity contour images in one embodiment；
Fig. 9 is right thoracic cavity contour images in one embodiment；
Figure 10 is area image to be analyzed in one embodiment；
Figure 11 is the structural block diagram of one embodiment cardiac equipment for reconstructing image；
Figure 12 is the structural block diagram of maximum intensity projection module in one embodiment；
Figure 13 is the structural block diagram of maximum intensity projection unit in one embodiment；
Figure 14 is the structural block diagram that boundary obtains module in one embodiment；
Figure 15 is the structural block diagram of left margin computing unit in one embodiment；
Figure 16 is the structural block diagram of right margin computing unit in one embodiment；
Figure 17 is the structural block diagram of coboundary computing unit in one embodiment；
Figure 18 is the internal structure chart of computer equipment in one embodiment.
Appended drawing reference：100 be maximum intensity projection module, 110 be the first image segmentation unit, 120 for open operating unit,
130 it is maximum intensity projection unit, 131 be the first maximum intensity projection subelement, 132 is that the second maximum intensity projection is single
Member, 133 be difference computation subunit, 200 be boundary acquisition module, 210 be the second image segmentation unit, 220 be left margin meter
Calculate unit, 221 for the first connected domain computation subunit, 222 be left margin computation subunit, 230 be right margin computing unit,
231 be the second connected domain computation subunit, 232 be right margin computation subunit, 240 be coboundary computing unit, 241 for point
Analysis region establish subelement, 242 be third connecting domain computation subunit, 243 be to weigh centered on coboundary computation subunit, 300
Modeling block, 400 are image reconstruction module.
Specific embodiment
It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood
The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not
For limiting the application.
Ct apparatus (CT) generally includes rack, scanning bed and for doctor's operation console.Rack
Side be provided with bulb, the side opposite with bulb is provided with detector.Console is that control bulb and detector carry out
The computer equipment of scanning, computer equipment is also used to the collected data of pick-up probe, and carries out processing reconstructed to data,
Ultimately form CT image.When being scanned using CT, patient is lain on scanning bed, by the scanning bed hole that patient is sent into rack
In diameter, the bulb being arranged in rack issues X-ray, and X-ray passes through patient and is received by a detector to form data, and data are transmitted
To computer equipment, computer equipment carries out preliminary treatment to data and image reconstruction obtains CT image.
In one embodiment, as shown in Figure 1, providing a kind of cardiac image method for reconstructing, include the following steps：
Step S102 is obtained using multiple preview images as input according to multiple preview images and maximum intensity projection
Thoracic cavity contour images.
Specifically, preview image is that computer equipment carries out the original graph obtained after preliminary treatment to the data received
Picture.Wherein preview image is the obtained original image in the big reconstruction visual field in at least visual field 500mm, can reflect the entire thoracic cavity of target
Structure.In medical imaging technology field, in order to enhance the imaging effect at patients target position, it will usually be injected to target site
Or take contrast agent.Preview image is injection or takes the laggard data got of contrast agent, and carries out preliminary treatment to data
The original image in obtained thoracic cavity.Multiple preview images are chosen, multiple preview images are split respectively, obtain reflection thoracic cavity
First image sequence of bone site and contrast locations.It carries out out operation respectively to the first image sequence again, obtains second
Image sequence.Wherein opening operation makes the profile of object become smooth, and can disconnect narrow interruption and eliminate thin protrusion.
Thinner edge can be got rid of by carrying out out operation respectively to the first image sequence, obtain the second figure of reflection contrast locations
As sequence.Maximum intensity projection is carried out to obtained the first image sequence and the second image sequence respectively and obtains thoracic cavity profile diagram
Picture.
The boundary position of thoracic cavity contour images is calculated according to thoracic cavity contour images by step S104.
Specifically, the boundary position of thoracic cavity contour images includes：Leftmost border position, rightmost circle position and the top
Boundary position.Wherein, leftmost border position is the physical location on the leftmost boundary of thoracic cavity contour images；Rightmost circle position is chest
The physical location on the boundary of chamber contour images rightmost；The top circle position is the physics on the boundary of thoracic cavity contour images the top
Position.According to obtained thoracic cavity contour images, thoracic cavity contour images are split, obtain left thoracic cavity contour images and right thoracic cavity
Contour images.Leftmost border is calculated further according to left thoracic cavity contour images, right thoracic cavity contour images and thoracic cavity contour images
Position, rightmost circle position and the top circle position.
Reconstructed center is calculated according to the boundary position of thoracic cavity contour images in step S106.
Specifically, according to the leftmost border position, rightmost circle position and the top circle position being calculated.It can be true
Determine the center of thoracic cavity profile, and is determined reconstructed center on thoracic cavity profile center is to the left in intrathoracic position according to heart
Position at.
Step S108 carries out leggy reconstruction according to reconstructed center and default field range of rebuilding, obtains cardiac image.
Specifically, leggy reconstruction is carried out according to determining reconstructed center and default field range of rebuilding, due to coronal
Artery is not constant, but the shape of bending change in the position of thoracic cavity axial direction.Therefore the default reconstruction visual field cannot be set
That sets is too small, and it is 80mm that the visual field is preferably rebuild in the present embodiment.Wherein rebuilding the visual field is the visual field rebuild for leggy.
Above-mentioned cardiac image method for reconstructing, by using multiple preview images as input, and according to multiple preview images with
And maximum intensity projection algorithm obtains thoracic cavity contour images, calculates leftmost border position, rightmost further according to thoracic cavity contour images
Boundary position and the top circle position, obtain reconstructed center.Multiphase is carried out according to reconstructed center and default field range of rebuilding
Position is rebuild, and obtains cardiac image.It is comparing for 200mm with traditional reconstruction visual field, the pixel resolution of the present embodiment improves
Nearly three times, and optimize by the calculating process of global Pixel calcualting and local Pixel calcualting, improve operation efficiency.
In one embodiment, as shown in Fig. 2 and Fig. 5-7, a kind of method of determining thoracic cavity contour images is provided, is wrapped
Include following steps：
Multiple preview images are split using multiple preview images as input, obtain the first image sequence by step S202
Column.
Specifically, preview image is that computer equipment carries out the original graph obtained after preliminary treatment to the data received
Picture.Wherein preview image is the obtained original image in the big reconstruction visual field in at least visual field 500mm, can reflect the entire thoracic cavity of target
Structure.In medical imaging technology field, in order to enhance the imaging effect at patients target position, it will usually be injected to target site
Or take contrast agent.Preview image is injection or takes the laggard data got of contrast agent, and carries out preliminary treatment to data
The original image in obtained thoracic cavity.Multiple preview images are chosen, multiple preview images are split respectively, obtain reflection thoracic cavity
First image sequence of bone site and contrast locations.
In one embodiment, can be according to the equation for being split to obtain the first image sequence to multiple preview images：
Wherein IPOz(i, j) represents the pixel of z preview images；The threshold value of TB1 expression high density structures；IPB1 is pair
Preview image is split to obtain the image of bone and contrast agent；IPB1z(i, j) represents image z of bone and contrast agent
The pixel of image.
Step S204 carries out out operation respectively to the first image sequence, obtains the second image sequence.
Specifically, according to the first obtained image sequence, operation is carried out out to the first image sequence respectively, obtains the second figure
As sequence.It opens operation wherein to make the profile of object become smooth, and disconnects narrow interruption and eliminate thin protrusion.To
One image sequence carries out out operation respectively to get rid of thinner edge, obtain the second image sequence of reflection contrast locations
Column.
In one embodiment, the first image sequence is carried out out respectively operation obtain the equation of the second image sequence can
For：
Wherein, IPB1 is to be split to obtain the image of bone and contrast agent to preview image；IPB2 is to bone and to make
The image of shadow agent carries out out the image of the contrast agent obtained after operation；Se_b1 is the structural element in morphology basic operation.
Step S206 carries out maximum intensity projection to the first image sequence and the second image sequence respectively and obtains thoracic cavity wheel
Wide image.
Specifically, maximum intensity projection is carried out to obtained the first image sequence and the second image sequence respectively and obtains thoracic cavity
Contour images.Wherein maximum intensity projection is by calculating the maximal density encountered on the ray of every, patients target position
Pixel and generate.I.e. when fiber optic bundle passes through the original image at patients target position, the maximum pixel of density is protected in image
It stays, and is projected on a two-dimensional surface, to form the maximum-density projection image at patients target position.
More specifically, maximum intensity projection is axially carried out in thoracic cavity to the first image sequence, obtains the maximum of the first image
Intensity Projection image.When fiber optic bundle is when thoracic cavity is axially through the image of bone and contrast agent, the maximum pixel of density in image
It is retained, and is projected on a two-dimensional surface, to form the maximum-density projection image of the first image.
In one embodiment, maximum intensity projection is axially carried out in thoracic cavity to the first image sequence and obtains the first image
The equation of maximum-density projection image can be：
IMB1(i, j)=max (IPB1z(i, j)), z=1,2 ... N
Wherein, IPB1z(i, j) is the pixel of image sequence z images of bone and contrast agent；N is bone and radiography
The number of the image of agent；IMB1For the maximum-density projection image of bone and the image of contrast agent；IMB1(i, j) represent bone and
The pixel of the maximum-density projection image of the image of contrast agent.
Maximum intensity projection is axially carried out in thoracic cavity to the second image sequence, obtains the maximum intensity projection figure of the second image
Picture.When fiber optic bundle is in image of the thoracic cavity axially through contrast agent, the maximum pixel of density is retained in image, and is projected to
On one two-dimensional surface, to form the maximum-density projection image of the second image.
In one embodiment, maximum intensity projection is axially carried out in thoracic cavity to the second image sequence and obtains the second image
The equation of maximum-density projection image can be：
IMB2(i, j)=max (IPB2z(i, j)), z=1,2 ... N
Wherein, IPB2z(i, j) is the pixel of image sequence z images of contrast agent；N is of the image of contrast agent
Number；IMB2For the maximum-density projection image of the image of contrast agent；IMB2The maximal density that (i, j) represents the image of contrast agent is thrown
The pixel of shadow image.
The maximum-density projection image of first image and the maximum-density projection image of the second image are made the difference, thoracic cavity is obtained
Contour images.
In one embodiment, by the maximum intensity projection figure of the maximum-density projection image of the first image and the second image
As the equation for making the difference to obtain thoracic cavity contour images can be：
Wherein, IMB1For the maximum-density projection image of bone and the image of contrast agent；IMB1(i, j) represents bone and makes
The pixel of the maximum-density projection image of the image of shadow agent；IMB2For the maximum-density projection image of the image of contrast agent；IMB2
(i, j) represents the pixel of the maximum-density projection image of the image of contrast agent；IMst is thoracic cavity contour images；IMst(i, j) is chest
The pixel of chamber contour images.
The method of above-mentioned determining thoracic cavity contour images obtains the first image sequence by being split to multiple preview images
Column.It carries out out operation respectively to the first image sequence, the second image sequence is obtained, respectively to the first image sequence and the second figure
As sequence carries out maximum intensity projection, and by the maximum intensity projection of the maximum-density projection image of the first image and the second image
Image makes the difference, and obtains thoracic cavity contour images.By maximum intensity projection can be more accurate obtain thoracic cavity contour images, can
It is significantly more efficient to remove thinner interference.
In one embodiment, as shown in Fig. 3 and Fig. 8-10, a kind of boundary bit for calculating thoracic cavity contour images is provided
The method set, includes the following steps：
Step S302 is split thoracic cavity contour images, obtains left thoracic cavity contour images and right thoracic cavity contour images.
Specifically, the boundary position of thoracic cavity contour images includes：Leftmost border position, rightmost circle position and the top
Boundary position.Wherein, leftmost border position is the physical location on the leftmost boundary of thoracic cavity contour images；Rightmost circle position is chest
The physical location on the boundary of chamber contour images rightmost；The top circle position is the physics on the boundary of thoracic cavity contour images the top
Position.According to obtained thoracic cavity contour images, thoracic cavity contour images are split, obtain left thoracic cavity contour images and right thoracic cavity
Contour images.
In one embodiment, thoracic cavity contour images are split to obtain left thoracic cavity contour images and right thoracic cavity profile diagram
The equation of picture can be：
Wherein, IMst is thoracic cavity contour images；IMst(i, j) is the pixel of thoracic cavity contour images；IMlFor left thoracic cavity profile
Image；IMl(i, j) is the pixel of left thoracic cavity contour images；IMrFor right thoracic cavity contour images；IMr(i, j) is right thoracic cavity profile diagram
The pixel of picture；M is the picture element matrix size of image.
Leftmost border position is calculated according to left thoracic cavity contour images in step S304.
Specifically, calculating largest connected domain according to left thoracic cavity contour images, and choose the left thoracic cavity wheel in largest connected domain
Wide image.Wherein, connected domain is a region on complex plane, does a simple closed curve if appointed wherein, and closed curve
Inside always belong to this region, just this region is referred to as connected domain.It chooses right in the left thoracic cavity contour images in largest connected domain
Upper parts of images where being preced with, and leftmost border position is calculated.
In one embodiment, the equation that leftmost border position is calculated can be：
Pixl=min (j | (i, j) ∈ IMl(i, j))
Wherein, PixlFor the most left location of pixels of left thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PoslFor
Thoracic cavity contour images leftmost border position；M is the picture element matrix size of image；IMlFor left thoracic cavity contour images.
Rightmost circle position is calculated according to right thoracic cavity contour images in step S306.
Specifically, calculating largest connected domain according to right thoracic cavity contour images, and choose the right thoracic cavity wheel in largest connected domain
Wide image.Wherein, connected domain is a region on complex plane, does a simple closed curve if appointed wherein, and closed curve
Inside always belong to this region, just this region is referred to as connected domain.It chooses right in the right thoracic cavity contour images in largest connected domain
Upper parts of images where being preced with, and rightmost circle position is calculated.
In one embodiment, the equation that rightmost circle position is calculated can be：
Pixr=max (j | (i, j) ∈ IMr(i, j))
Wherein, PixrFor the most right location of pixels of right thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PosrFor
Thoracic cavity contour images rightmost circle position；M is the picture element matrix size of image；IMrFor right thoracic cavity contour images.
The top circle position is calculated according to thoracic cavity contour images in step S308.
Specifically, according to the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images and second
The lowermost position of image sequence is set establishes region to be analyzed in the contour images of thoracic cavity；Wherein, the second image sequence is reflection radiography
The image sequence of agent position.Largest connected domain is calculated according to region to be analyzed, and chooses the thoracic cavity profile diagram in largest connected domain
Picture；According to the thoracic cavity contour images in largest connected domain, the top circle position is calculated.
In one embodiment, according to the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images
And second the lowermost position of image sequence set and establish the equation in region to be analyzed in the contour images of thoracic cavity and can be：
Irib=IMst(1：Id, id1：id2)；
Wherein, id is the minimum y location of the image sequence of contrast agent, and id1 is the most right x position of left thoracic cavity contour images,
Id2 is the most left x position of right thoracic cavity contour images；Irib is region to be analyzed；IMst is thoracic cavity contour images.
The side of the top circle position is calculated according to the thoracic cavity contour images in largest connected domain in one embodiment
Cheng Kewei：
Wherein, PixupFor the most upper location of pixels of thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PosupFor
Thoracic cavity contour images the top circle position；M is the picture element matrix size of image.
The method of the boundary position of above-mentioned calculating thoracic cavity contour images, is split thoracic cavity contour images, obtains left chest
Chamber contour images and right thoracic cavity contour images, further according to left thoracic cavity contour images, right thoracic cavity contour images and thoracic cavity profile diagram
As leftmost border position, rightmost circle position and the top circle position is calculated.By being acquisition extremely to thoracic cavity boundary,
Heart can be further determined more accurately in thoracic cavity to the center of thoracic cavity profile in calculating that can be more accurate
Location keeps the selection of cardiac position more accurate.
In one embodiment, as shown in figs. 4 through 10, a kind of cardiac image method for reconstructing is provided, is included the following steps：
Multiple preview images are split using multiple preview images as input, obtain the first image sequence by step S402
Column.
Specifically, preview image is that computer equipment carries out the original graph obtained after preliminary treatment to the data received
Picture.Wherein preview image is the obtained original image in the big reconstruction visual field in at least visual field 500mm, can reflect the entire thoracic cavity of target
Structure.In medical imaging technology field, in order to enhance the imaging effect at patients target position, it will usually be injected to target site
Or take contrast agent.Preview image is injection or takes the laggard data got of contrast agent, and carries out preliminary treatment to data
The original image in obtained thoracic cavity.Multiple preview images are chosen, multiple preview images are split respectively, obtain reflection thoracic cavity
First image sequence of bone site and contrast locations.
In one embodiment, can be according to the equation for being split to obtain the first image sequence to multiple preview images：
Wherein IPOz(i, j) represents the pixel of z preview images；The threshold value of TB1 expression high density structures；IPB1 is pair
Preview image is split to obtain the image of bone and contrast agent；IPB1z(i, j) represents image z of bone and contrast agent
The pixel of image.
Step S404 carries out out operation respectively to the first image sequence, obtains the second image sequence.
Specifically, according to the first obtained image sequence, operation is carried out out to the first image sequence respectively, obtains the second figure
As sequence.It opens operation wherein to make the profile of object become smooth, and disconnects narrow interruption and eliminate thin protrusion.To
One image sequence carries out out operation respectively to get rid of thinner edge, obtain the second image sequence of reflection contrast locations
Column.
In one embodiment, the first image sequence is carried out out respectively operation obtain the equation of the second image sequence can
For：
Wherein, IPB1 is to be split to obtain the image of bone and contrast agent to preview image；IPB2 is to bone and to make
The image of shadow agent carries out out the image of the contrast agent obtained after operation；Se_b1 is the structural element in morphology basic operation.
Step S406 carries out maximum intensity projection to the first image sequence and the second image sequence respectively and obtains thoracic cavity wheel
Wide image.
Specifically, maximum intensity projection is carried out to obtained the first image sequence and the second image sequence respectively and obtains thoracic cavity
Contour images.Wherein maximum intensity projection is by calculating the maximal density encountered on the ray of every, patients target position
Pixel and generate.I.e. when fiber optic bundle passes through the original image at patients target position, the maximum pixel of density is protected in image
It stays, and is projected on a two-dimensional surface, to form the maximum-density projection image at patients target position.
More specifically, maximum intensity projection is axially carried out in thoracic cavity to the first image sequence, obtains the maximum of the first image
Intensity Projection image.When fiber optic bundle is when thoracic cavity is axially through the image of bone and contrast agent, the maximum pixel of density in image
It is retained, and is projected on a two-dimensional surface, to form the maximum-density projection image of the first image.
In one embodiment, maximum intensity projection is axially carried out in thoracic cavity to the first image sequence and obtains the first image
The equation of maximum-density projection image can be：
IMB1(i, j)=max (IPB1z(i, j)), z=1,2 ... N
Wherein, IPB1z(i, j) is the pixel of image sequence z images of bone and contrast agent；N is bone and radiography
The number of the image of agent；IMB1For the maximum-density projection image of bone and the image of contrast agent；IMB1(i, j) represent bone and
The pixel of the maximum-density projection image of the image of contrast agent.
Maximum intensity projection is axially carried out in thoracic cavity to the second image sequence, obtains the maximum intensity projection figure of the second image
Picture.When fiber optic bundle is in image of the thoracic cavity axially through contrast agent, the maximum pixel of density is retained in image, and is projected to
On one two-dimensional surface, to form the maximum-density projection image of the second image.
In one embodiment, maximum intensity projection is axially carried out in thoracic cavity to the second image sequence and obtains the second image
The equation of maximum-density projection image can be：
IMB2(i, j)=max (IPB2z(i, j)), z=1,2 ... N
Wherein, IPB2z(i, j) is the pixel of image sequence z images of contrast agent；N is of the image of contrast agent
Number；IMB2For the maximum-density projection image of the image of contrast agent；IMB2The maximal density that (i, j) represents the image of contrast agent is thrown
The pixel of shadow image.
The maximum-density projection image of first image and the maximum-density projection image of the second image are made the difference, thoracic cavity is obtained
Contour images.
In one embodiment, by the maximum intensity projection figure of the maximum-density projection image of the first image and the second image
As the equation for making the difference to obtain thoracic cavity contour images can be：
Wherein, IMB1For the maximum-density projection image of bone and the image of contrast agent；IMB1(i, j) represents bone and makes
The pixel of the maximum-density projection image of the image of shadow agent；IMB2For the maximum-density projection image of the image of contrast agent；IMB2
(i, j) represents the pixel of the maximum-density projection image of the image of contrast agent；IMst is thoracic cavity contour images；IMst(i, j) is chest
The pixel of chamber contour images.
Step S408 is split thoracic cavity contour images, obtains left thoracic cavity contour images and right thoracic cavity contour images.
Specifically, the boundary position of thoracic cavity contour images includes：Leftmost border position, rightmost circle position and the top
Boundary position.Wherein, leftmost border position is the physical location on the leftmost boundary of thoracic cavity contour images；Rightmost circle position is chest
The physical location on the boundary of chamber contour images rightmost；The top circle position is the physics on the boundary of thoracic cavity contour images the top
Position.According to obtained thoracic cavity contour images, thoracic cavity contour images are split, obtain left thoracic cavity contour images and right thoracic cavity
Contour images.
In one embodiment, thoracic cavity contour images are split to obtain left thoracic cavity contour images and right thoracic cavity profile diagram
The equation of picture can be：
Wherein, IMst is thoracic cavity contour images；IMst(i, j) is the pixel of thoracic cavity contour images；IMlFor left thoracic cavity profile
Image；IMl(i, j) is the pixel of left thoracic cavity contour images；IMrFor right thoracic cavity contour images；IMr(i, j) is right thoracic cavity profile diagram
The pixel of picture；M is the picture element matrix size of image.
Leftmost border position is calculated according to left thoracic cavity contour images in step S410.
Specifically, calculating largest connected domain according to left thoracic cavity contour images, and choose the left thoracic cavity wheel in largest connected domain
Wide image.Wherein, connected domain is a region on complex plane, does a simple closed curve if appointed wherein, and closed curve
Inside always belong to this region, just this region is referred to as connected domain.It chooses right in the left thoracic cavity contour images in largest connected domain
Upper parts of images where being preced with, and leftmost border position is calculated.
In one embodiment, the equation that leftmost border position is calculated can be：
Pixl=min (j | (i, j) ∈ IMl(i, j))
Wherein, PixlFor the most left location of pixels of left thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PoslFor
Thoracic cavity contour images leftmost border position；M is the picture element matrix size of image；IMlFor left thoracic cavity contour images.
Rightmost circle position is calculated according to right thoracic cavity contour images in step S412.
Specifically, calculating largest connected domain according to right thoracic cavity contour images, and choose the right thoracic cavity wheel in largest connected domain
Wide image.Wherein, connected domain is a region on complex plane, does a simple closed curve if appointed wherein, and closed curve
Inside always belong to this region, just this region is referred to as connected domain.It chooses right in the right thoracic cavity contour images in largest connected domain
Upper parts of images where being preced with, and rightmost circle position is calculated.
In one embodiment, the equation that rightmost circle position is calculated can be：
Pixr=max (j | (i, j) ∈ IMr (i, j))
Wherein, PixrFor the most right location of pixels of right thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PosrFor
Thoracic cavity contour images rightmost circle position；M is the picture element matrix size of image；IMrFor right thoracic cavity contour images.
The top circle position is calculated according to thoracic cavity contour images in step S414.
Specifically, according to the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images and second
The lowermost position of image sequence is set establishes region to be analyzed in the contour images of thoracic cavity；It is calculated according to region to be analyzed largest connected
Domain, and choose the thoracic cavity contour images in largest connected domain；According to the thoracic cavity contour images in largest connected domain, it is calculated most
Coboundary position.
In one embodiment, according to the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images
And second the lowermost position of image sequence set and establish the equation in region to be analyzed in the contour images of thoracic cavity and can be：
Irib=IMst(1：Id, id1：id2)；
Wherein, id is the minimum y location of the second image sequence, and id1 is the most right x position of left thoracic cavity contour images, and id2 is
The most left x position of right thoracic cavity contour images；Irib is region to be analyzed；IMst is thoracic cavity contour images.
In one embodiment, the side of the top circle position is calculated according to the thoracic cavity contour images in largest connected domain
Cheng Kewei：
Wherein, PixupFor the most upper location of pixels of thoracic cavity contour images；SpacingpreFor the resolution ratio of pixel；PosupFor
Thoracic cavity contour images the top circle position；M is the picture element matrix size of image.
Reconstructed center is calculated according to the boundary position of thoracic cavity contour images in step S416.
Specifically, according to the leftmost border position, rightmost circle position and the top circle position being calculated.It can be true
Determine the center of thoracic cavity profile, and is determined reconstructed center on thoracic cavity profile center is to the left in intrathoracic position according to heart
Position at.
In one embodiment, can be according to the equation that reconstructed center is calculated in the boundary position of thoracic cavity contour images：
Wherein, OffSetX is used to adjust the distance that reconstructed center is deviated to the left profile center；OffseY is used to adjust reconstruction
Center is upwardly deviated from the distance at profile center；PoslFor thoracic cavity contour images leftmost border position；PosrMost for thoracic cavity contour images
Right margin position；PosupFor thoracic cavity contour images the top circle position.
Step S418 carries out leggy reconstruction according to reconstructed center and default field range of rebuilding, obtains cardiac image.
Specifically, leggy reconstruction is carried out according to determining reconstructed center and default field range of rebuilding, due to coronal
Artery is not constant, but the shape of bending change in the position of thoracic cavity axial direction.Therefore the default reconstruction visual field cannot be set
That sets is too small, and it is 80mm that the visual field is preferably rebuild in the present embodiment.Wherein rebuilding the visual field is the visual field rebuild for leggy.
Above-mentioned cardiac image method for reconstructing, device, computer equipment and storage medium, by using multiple preview images as
Input, and thoracic cavity contour images are obtained according to multiple preview images and maximum intensity projection algorithm, further according to thoracic cavity profile diagram
As calculating leftmost border position, rightmost circle position and the top circle position, reconstructed center is obtained.According to reconstructed center and
Default field range of rebuilding carries out leggy reconstruction, obtains cardiac image.By first choosing cardiac position, according to cardiac position into
Row leggy image reconstruction, reduces the data volume of input, further reduces the time of data input, improves operation effect
Rate.
It should be understood that although each step in the flow chart of Fig. 1-4 is successively shown according to the instruction of arrow,
These steps are not that the inevitable sequence according to arrow instruction successively executes.Unless expressly stating otherwise herein, these steps
Execution there is no stringent sequences to limit, these steps can execute in other order.Moreover, at least one in Fig. 1-4
Part steps may include that perhaps these sub-steps of multiple stages or stage are not necessarily in synchronization to multiple sub-steps
Completion is executed, but can be executed at different times, the execution sequence in these sub-steps or stage is also not necessarily successively
It carries out, but can be at least part of the sub-step or stage of other steps or other steps in turn or alternately
It executes.
In one embodiment, as shown in figure 11, a kind of cardiac image reconstructing device is provided, including：Maximal density is thrown
Shadow module 100, boundary obtain module 200, center rebuilds module 300 and image reconstruction module 400, wherein：
Maximum intensity projection module 100 is used for using multiple preview images as input, according to multiple preview images and most
Big Intensity Projection obtains thoracic cavity contour images；
Boundary obtains module 200, for the boundary position of thoracic cavity contour images to be calculated according to thoracic cavity contour images；
Center rebuilds module 300, for reconstructed center to be calculated according to the boundary position of thoracic cavity contour images；
Image reconstruction module 400 obtains the heart for carrying out leggy reconstruction according to reconstructed center and default field range
Dirty image.
As shown in figure 12, a kind of structural block diagram of maximum intensity projection module 100 is provided, wherein maximum intensity projection mould
Block 100 includes：First image segmentation unit 110 opens operating unit 120, maximum intensity projection unit 130.
First image segmentation unit 110, for using multiple preview images as input, multiple preview images to be divided
It cuts, obtains the first image sequence；
Operating unit 120 is opened, for carrying out out operation respectively to the first image sequence, obtains the second image sequence；
Maximum intensity projection unit 130 carries out the first image sequence and the second image sequence for respectively maximum close
Degree projection obtains thoracic cavity contour images.
As shown in figure 13, a kind of structural block diagram of maximum intensity projection unit 130 is provided, wherein maximum intensity projection list
First 130 include：First maximum intensity projection subelement 131, the second maximum intensity projection subelement 132, difference computation subunit
133。
First maximum intensity projection subelement 131, for axially carrying out maximal density throwing in thoracic cavity to the first image sequence
Shadow obtains the maximum-density projection image of the first image；
Second maximum intensity projection subelement 132, for axially carrying out maximal density throwing in thoracic cavity to the second image sequence
Shadow obtains the maximum-density projection image of the second image；
Difference computation subunit 133, for the maximum-density projection image of the first image and the maximum of the second image is close
Degree projected image makes the difference, and obtains thoracic cavity contour images.
As shown in figure 14, a kind of structural block diagram of boundary acquisition module 200 is provided, wherein boundary obtains module 200 and wraps
It includes：Second image segmentation unit 210, left margin computing unit 220, right margin computing unit 230, coboundary computing unit 240.
Second image segmentation unit 210 obtains left thoracic cavity contour images and the right side for being split to thoracic cavity contour images
Thoracic cavity contour images.
Left margin computing unit 220, for leftmost border position to be calculated according to left thoracic cavity contour images.
Right margin computing unit 230, for rightmost circle position to be calculated according to right thoracic cavity contour images.
Coboundary computing unit 240, for the top circle position to be calculated according to thoracic cavity contour images.
As shown in figure 15, a kind of structural block diagram of left margin computing unit 220 is provided, wherein left margin computing unit
220 include：First connected domain computation subunit 221, left margin computation subunit 222.
First connected domain computation subunit 221 for calculating largest connected domain according to left thoracic cavity contour images, and is chosen most
Left thoracic cavity contour images in big connected domain；
Left margin computation subunit 222, for choosing in the left thoracic cavity contour images in largest connected domain where right hat
Upper parts of images, and leftmost border position is calculated.
As shown in figure 16, a kind of structural block diagram of right margin computing unit 230 is provided, wherein right margin computing unit
230 include：Second connected domain computation subunit 231, right margin computation subunit 232.
Second connected domain computation subunit 231 for calculating largest connected domain according to right thoracic cavity contour images, and is chosen most
Right thoracic cavity contour images in big connected domain.
Right margin computation subunit 232, for choosing in the right thoracic cavity contour images in largest connected domain where right hat
Upper parts of images, and rightmost circle position is calculated.
As shown in figure 17, a kind of structural block diagram of coboundary computing unit 240 is provided, wherein coboundary computing unit
240 include：Subelement 241, third connecting domain computation subunit 242, coboundary computation subunit 243 are established in region to be analyzed.
Subelement 241 is established in region to be analyzed, for the least significant according to left thoracic cavity contour images, right thoracic cavity profile diagram
The lowermost position of the leftmost position of picture and the second image sequence, which is set, establishes region to be analyzed in the contour images of thoracic cavity.
Third connecting domain computation subunit 242 for calculating largest connected domain according to region to be analyzed, and chooses most Dalian
Thoracic cavity contour images in logical domain.
Coboundary computation subunit 243, for being calculated topmost according to the thoracic cavity contour images in largest connected domain
Boundary position.
Specific about cardiac image reconstructing device limits the limit that may refer to above for cardiac image method for reconstructing
Fixed, details are not described herein.Modules in above-mentioned cardiac image reconstructing device can fully or partially through software, hardware and its
Combination is to realize.Above-mentioned each module can be embedded in the form of hardware or independently of in the processor in computer equipment, can also be with
It is stored in the memory in computer equipment in a software form, in order to which processor calls the above modules of execution corresponding
Operation.
In one embodiment, a kind of computer equipment is provided, which can be terminal, internal structure
Figure can be as shown in figure 18.The computer equipment includes the processor connected by system bus, memory, network interface, shows
Display screen and input unit.Wherein, the processor of the computer equipment is for providing calculating and control ability.The computer equipment
Memory includes non-volatile memory medium, built-in storage.The non-volatile memory medium is stored with operating system and computer
Program.The built-in storage provides environment for the operation of operating system and computer program in non-volatile memory medium.The meter
The network interface for calculating machine equipment is used to communicate with external terminal by network connection.When the computer program is executed by processor
To realize a kind of cardiac image method for reconstructing.The display screen of the computer equipment can be liquid crystal display or electric ink is aobvious
Display screen, the input unit of the computer equipment can be the touch layer covered on display screen, be also possible to computer equipment shell
Key, trace ball or the Trackpad of upper setting can also be external keyboard, Trackpad or mouse etc..
It will be understood by those skilled in the art that structure shown in Figure 18, only part relevant to application scheme
The block diagram of structure, does not constitute the restriction for the computer equipment being applied thereon to application scheme, and specific computer is set
Standby may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
In one embodiment, a kind of computer equipment, including memory and processor are provided, is stored in memory
Computer program, the processor realize following steps when executing computer program：
Using multiple preview images as input, thoracic cavity profile diagram is obtained according to multiple preview images and maximum intensity projection
Picture；The boundary position of thoracic cavity contour images is calculated according to thoracic cavity contour images；According to the boundary position of thoracic cavity contour images
Reconstructed center is calculated；Leggy reconstruction is carried out according to reconstructed center and default field range of rebuilding, obtains cardiac image.
In one embodiment, following steps are also realized when processor executes computer program：
Using multiple preview images as input, multiple preview images are split, the first image sequence is obtained；To first
Image sequence carries out out operation respectively, obtains the second image sequence；Respectively to the first image sequence and the second image sequence into
Row maximum intensity projection obtains thoracic cavity contour images.
In one embodiment, following steps are also realized when processor executes computer program：
Thoracic cavity contour images are split, left thoracic cavity contour images and right thoracic cavity contour images are obtained；According to left thoracic cavity
Leftmost border position is calculated in contour images；Rightmost circle position is calculated according to right thoracic cavity contour images；According to thoracic cavity
The top circle position is calculated in contour images.
In one embodiment, following steps are also realized when processor executes computer program：
Using multiple preview images as input, multiple preview images are split, the first image sequence is obtained；To first
Image sequence carries out out operation respectively, obtains the second image sequence；Respectively to the first image sequence and the second image sequence into
Row maximum intensity projection obtains thoracic cavity contour images；Thoracic cavity contour images are split, left thoracic cavity contour images and the right side are obtained
Thoracic cavity contour images；Leftmost border position is calculated according to left thoracic cavity contour images；It is calculated according to right thoracic cavity contour images
To rightmost circle position；The top circle position is calculated according to thoracic cavity contour images；According to the boundary bit of thoracic cavity contour images
It sets and reconstructed center is calculated；Leggy reconstruction is carried out according to reconstructed center and default field range of rebuilding, obtains cardiod diagram
Picture.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated
Machine program realizes following steps when being executed by processor：
Using multiple preview images as input, thoracic cavity profile diagram is obtained according to multiple preview images and maximum intensity projection
Picture；The boundary position of thoracic cavity contour images is calculated according to thoracic cavity contour images；According to the boundary position of thoracic cavity contour images
Reconstructed center is calculated；Leggy reconstruction is carried out according to reconstructed center and default field range of rebuilding, obtains cardiac image.
In one embodiment, following steps are also realized when computer program is executed by processor：
Using multiple preview images as input, multiple preview images are split, the first image sequence is obtained；To first
Image sequence carries out out operation respectively, obtains the second image sequence；Respectively to the first image sequence and the second image sequence into
Row maximum intensity projection obtains thoracic cavity contour images.
In one embodiment, following steps are also realized when computer program is executed by processor：
Thoracic cavity contour images are split, left thoracic cavity contour images and right thoracic cavity contour images are obtained；According to left thoracic cavity
Leftmost border position is calculated in contour images；Rightmost circle position is calculated according to right thoracic cavity contour images；According to thoracic cavity
The top circle position is calculated in contour images.
In one embodiment, following steps are also realized when computer program is executed by processor：
Using multiple preview images as input, multiple preview images are split, the first image sequence is obtained；To first
Image sequence carries out out operation respectively, obtains the second image sequence；Respectively to the first image sequence and the second image sequence into
Row maximum intensity projection obtains thoracic cavity contour images；Thoracic cavity contour images are split, left thoracic cavity contour images and the right side are obtained
Thoracic cavity contour images；Leftmost border position is calculated according to left thoracic cavity contour images；It is calculated according to right thoracic cavity contour images
To rightmost circle position；The top circle position is calculated according to thoracic cavity contour images；According to the boundary bit of thoracic cavity contour images
It sets and reconstructed center is calculated；Leggy reconstruction is carried out according to reconstructed center and default field range of rebuilding, obtains cardiod diagram
Picture.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer
In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein,
To any reference of memory, storage, database or other media used in each embodiment provided herein,
Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM
(PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include
Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms,
Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing
Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM
(RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment
In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance
Shield all should be considered as described in this specification.
The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application
Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.
Claims (10)
1. a kind of cardiac image method for reconstructing, which is characterized in that the method includes：
Using multiple preview images as input, thoracic cavity profile diagram is obtained according to preview image described in multiple and maximum intensity projection
Picture；
The boundary position of the thoracic cavity contour images is calculated according to thoracic cavity contour images；
Reconstructed center is calculated according to the boundary position of the thoracic cavity contour images；
Leggy reconstruction is carried out according to the reconstructed center and default field range of rebuilding, obtains cardiac image.
2. the method according to claim 1, wherein it is described using multiple preview images as input, according to multiple
The preview image and maximum intensity projection obtain thoracic cavity contour images, including：
Using multiple preview images as input, multiple described preview images are split, the first image sequence is obtained；
It carries out out operation respectively to the first image sequence, obtains the second image sequence；
Maximum intensity projection is carried out to the first image sequence and the second image sequence respectively and obtains thoracic cavity contour images.
3. according to the method described in claim 2, it is characterized in that, described respectively to the first image sequence and the second image sequence
Column carry out maximum intensity projection and obtain thoracic cavity contour images, including：
Maximum intensity projection is axially carried out in thoracic cavity to the first image sequence, obtains the maximum intensity projection figure of the first image
Picture；
Maximum intensity projection is axially carried out in thoracic cavity to second image sequence, obtains the maximum intensity projection figure of the second image
Picture；
The maximum-density projection image of the first image and the maximum-density projection image of the second image are made the difference, thoracic cavity is obtained
Contour images.
4. according to the method described in claim 2, it is characterized in that,
The boundary position of the thoracic cavity contour images includes：Leftmost border position, rightmost circle position and the top circle position.
5. according to the method described in claim 4, it is characterized in that, described be calculated the thoracic cavity according to thoracic cavity contour images
The boundary position of contour images, including：
Thoracic cavity contour images are split, left thoracic cavity contour images and right thoracic cavity contour images are obtained；
Leftmost border position is calculated according to left thoracic cavity contour images；
Rightmost circle position is calculated according to right thoracic cavity contour images；
The top circle position is calculated according to thoracic cavity contour images.
6. according to the method described in claim 5, it is characterized in that, described be calculated Far Left according to left thoracic cavity contour images
Boundary position, including：
Largest connected domain is calculated according to left thoracic cavity contour images, and chooses the left thoracic cavity contour images in the largest connected domain；
The upper parts of images in the left thoracic cavity contour images in the largest connected domain where right hat is chosen, and is calculated most left
Boundary position.
7. according to the method described in claim 5, it is characterized in that, described be calculated rightmost according to right thoracic cavity contour images
Boundary position, including：
Largest connected domain is calculated according to right thoracic cavity contour images, and chooses the right thoracic cavity contour images in largest connected domain；
The upper parts of images in the right thoracic cavity contour images in largest connected domain where right hat is chosen, and rightmost circle is calculated
Position.
8. according to the method described in claim 5, it is characterized in that, described be calculated the top circle according to thoracic cavity contour images
Position, including：
Most according to the least significant of left thoracic cavity contour images, the leftmost position of right thoracic cavity contour images and the second image sequence
Region to be analyzed is established in lower position in the contour images of thoracic cavity；
Largest connected domain is calculated according to the region to be analyzed, and chooses the thoracic cavity contour images in largest connected domain；
According to the thoracic cavity contour images in largest connected domain, the top circle position is calculated.
9. a kind of cardiac image reconstructing device, which is characterized in that described device includes：
Maximum intensity projection module is used for using multiple preview images as input, according to preview image described in multiple and maximum
Intensity Projection obtains thoracic cavity contour images；
Boundary obtains module, for the boundary position of the thoracic cavity contour images to be calculated according to thoracic cavity contour images；
Center rebuilds module, for reconstructed center to be calculated according to the boundary position of the thoracic cavity contour images；
Image reconstruction module obtains heart for carrying out leggy reconstruction according to reconstructed center and default field range of rebuilding
Image.
10. a kind of computer equipment, including memory and processor, the memory are stored with computer program, feature exists
In the step of processor realizes any one of claims 1 to 8 the method when executing the computer program.
Priority Applications (10)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810597965.1A CN108898582B (en) | 2018-06-11 | 2018-06-11 | Heart image reconstruction method and device and computer equipment |
| US16/437,003 US10950016B2 (en) | 2018-06-11 | 2019-06-11 | Systems and methods for reconstructing cardiac images |
| US16/437,006 US11024062B2 (en) | 2018-06-11 | 2019-06-11 | Systems and methods for evaluating image quality |
| US17/202,308 US11450038B2 (en) | 2018-06-11 | 2021-03-15 | Systems and methods for reconstructing cardiac images |
| US17/243,522 US11288849B2 (en) | 2018-06-11 | 2021-04-28 | Systems and methods for evaluating image quality |
| US17/326,374 US11367228B2 (en) | 2018-06-11 | 2021-05-21 | Systems and methods for evaluating image quality based on regularity degrees and sharpness degrees of images |
| US17/651,895 US11688110B2 (en) | 2018-06-11 | 2022-02-21 | Systems and methods for evaluating image quality |
| US17/933,469 US11915347B2 (en) | 2018-06-11 | 2022-09-19 | Systems and methods for reconstructing cardiac images |
| US18/341,737 US20230351649A1 (en) | 2018-06-11 | 2023-06-26 | Systems and methods for evaluating image quality |
| US18/587,896 US20240202995A1 (en) | 2018-06-11 | 2024-02-26 | Systems and methods for reconstructing cardiac images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810597965.1A CN108898582B (en) | 2018-06-11 | 2018-06-11 | Heart image reconstruction method and device and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108898582A true CN108898582A (en) | 2018-11-27 |
| CN108898582B CN108898582B (en) | 2021-08-17 |
Family
ID=64344697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810597965.1A Active CN108898582B (en) | 2018-06-11 | 2018-06-11 | Heart image reconstruction method and device and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108898582B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109697740A (en) * | 2018-12-26 | 2019-04-30 | 上海联影智能医疗科技有限公司 | Image rebuilding method, device and computer equipment |
| CN111062998A (en) * | 2019-12-11 | 2020-04-24 | 上海联影医疗科技有限公司 | Image reconstruction method, image reconstruction device, CT system and storage medium |
| US11288849B2 (en) | 2018-06-11 | 2022-03-29 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for evaluating image quality |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101727666A (en) * | 2008-11-03 | 2010-06-09 | 深圳迈瑞生物医疗电子股份有限公司 | Image segmentation method and device, and method for judging image inversion and distinguishing front side and back side of sternum |
| CN102622775A (en) * | 2012-01-12 | 2012-08-01 | 北京理工大学 | Heart real-time dynamic rebuilding technology based on model interpolation compensation |
| US20140219524A1 (en) * | 2011-05-11 | 2014-08-07 | Toshiba Medical Systems Corporation | Medical image processing apparatus and method thereof |
| CN104200500A (en) * | 2014-07-29 | 2014-12-10 | 沈阳东软医疗***有限公司 | Method and device for reconstructing heart image |
| CN104240180A (en) * | 2014-08-08 | 2014-12-24 | 沈阳东软医疗***有限公司 | Method and device for achieving automatic adjusting of images |
| CN104331914A (en) * | 2014-09-18 | 2015-02-04 | 东莞南方医大松山湖科技园有限公司 | Method and device for reconstructing heart phase image |
| CN106296764A (en) * | 2016-08-02 | 2017-01-04 | 上海联影医疗科技有限公司 | Image rebuilding method and system |
-
2018
- 2018-06-11 CN CN201810597965.1A patent/CN108898582B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101727666A (en) * | 2008-11-03 | 2010-06-09 | 深圳迈瑞生物医疗电子股份有限公司 | Image segmentation method and device, and method for judging image inversion and distinguishing front side and back side of sternum |
| US20140219524A1 (en) * | 2011-05-11 | 2014-08-07 | Toshiba Medical Systems Corporation | Medical image processing apparatus and method thereof |
| CN102622775A (en) * | 2012-01-12 | 2012-08-01 | 北京理工大学 | Heart real-time dynamic rebuilding technology based on model interpolation compensation |
| CN104200500A (en) * | 2014-07-29 | 2014-12-10 | 沈阳东软医疗***有限公司 | Method and device for reconstructing heart image |
| CN104240180A (en) * | 2014-08-08 | 2014-12-24 | 沈阳东软医疗***有限公司 | Method and device for achieving automatic adjusting of images |
| CN104331914A (en) * | 2014-09-18 | 2015-02-04 | 东莞南方医大松山湖科技园有限公司 | Method and device for reconstructing heart phase image |
| CN106296764A (en) * | 2016-08-02 | 2017-01-04 | 上海联影医疗科技有限公司 | Image rebuilding method and system |
Non-Patent Citations (2)
| Title |
|---|
| PRATOMO ADHI NUGROHO 等: "3D heart image reconstruction and visualization with marching cubes algorithm", 《2016 INTERNATIONAL CONFERENCE ON KNOWLEDGE CREATION AND INTELLIGENT COMPUTING (KCIC)》 * |
| 蒋峻 等: "基于中国数字人数据的心脏三维模型的重建", 《组织工程与重建外科杂志》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11288849B2 (en) | 2018-06-11 | 2022-03-29 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for evaluating image quality |
| US11367228B2 (en) | 2018-06-11 | 2022-06-21 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for evaluating image quality based on regularity degrees and sharpness degrees of images |
| US11688110B2 (en) | 2018-06-11 | 2023-06-27 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for evaluating image quality |
| CN109697740A (en) * | 2018-12-26 | 2019-04-30 | 上海联影智能医疗科技有限公司 | Image rebuilding method, device and computer equipment |
| CN109697740B (en) * | 2018-12-26 | 2020-09-22 | 上海联影智能医疗科技有限公司 | Image reconstruction method and device and computer equipment |
| US11494877B2 (en) | 2018-12-26 | 2022-11-08 | Shanghai United Imaging Intelligence Co., Ltd. | Systems and methods for image reconstruction |
| CN111062998A (en) * | 2019-12-11 | 2020-04-24 | 上海联影医疗科技有限公司 | Image reconstruction method, image reconstruction device, CT system and storage medium |
| CN111062998B (en) * | 2019-12-11 | 2023-12-12 | 上海联影医疗科技股份有限公司 | Image reconstruction method, device, CT system and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108898582B (en) | 2021-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11132792B2 (en) | Cross domain medical image segmentation | |
| JP7046240B2 (en) | Efficient sharing of medical imaging and medical imaging information | |
| US11568534B2 (en) | System and method of mitral valve quantification | |
| CN109544566B (en) | Coronary image segmentation method, device, computer equipment and storage medium | |
| Yuan et al. | SIPID: A deep learning framework for sinogram interpolation and image denoising in low-dose CT reconstruction | |
| Zhou et al. | MDPET: a unified motion correction and denoising adversarial network for low-dose gated PET | |
| CN108898582A (en) | Cardiac image method for reconstructing, device and computer equipment | |
| CN109199422A (en) | CT preview image rebuilds optimization method, device, computer equipment and storage medium | |
| US20180064409A1 (en) | Simultaneously displaying medical images | |
| CN109544657A (en) | Medical image iterative reconstruction approach, device, computer equipment and storage medium | |
| CN109712163A (en) | Coronary artery extracting method, device, image processing workstations and readable storage medium storing program for executing | |
| CN111476791A (en) | Image processing method, image processing apparatus, and non-transitory computer-readable medium | |
| CN109727242A (en) | Vessel centerline extracting method, device, computer equipment and storage medium | |
| CN113506331A (en) | Method, apparatus, computer device and storage medium for registering tissue and organ | |
| CN110084868A (en) | Method for correcting image, device, computer equipment and readable storage medium storing program for executing | |
| CN111028212A (en) | Key point detection method and device, computer equipment and storage medium | |
| Wu et al. | Dapc-net: Deformable alignment and pyramid context completion networks for video inpainting | |
| US9965875B2 (en) | Virtual projection image method | |
| Fu et al. | Oif-net: An optical flow registration-based pet/mr cross-modal interactive fusion network for low-count brain pet image denoising | |
| CN111161369B (en) | Image reconstruction storage method, device, computer equipment and storage medium | |
| CN113129418A (en) | Target surface reconstruction method, device, equipment and medium based on three-dimensional image | |
| CN110473241B (en) | Image registration method, storage medium and computer device | |
| Li et al. | Quad-Net: Quad-domain network for CT metal artifact reduction | |
| EP3272288B1 (en) | Apparatus and method for ct data reconstruction based on motion compensation | |
| CN108154541A (en) | Medical imaging devices and its operating method |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |