CN105007818B - Radiographic apparatus and image processing apparatus - Google Patents
Radiographic apparatus and image processing apparatus Download PDFInfo
- Publication number
- CN105007818B CN105007818B CN201480010340.0A CN201480010340A CN105007818B CN 105007818 B CN105007818 B CN 105007818B CN 201480010340 A CN201480010340 A CN 201480010340A CN 105007818 B CN105007818 B CN 105007818B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- image
- region
- determining section
- blood vessel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The radiographic apparatus of embodiment possesses profile measurement portion, correction coefficient determining section, DSA image correcting sections and control unit.Profile measurement portion measure respectively with the relevant profile of contrast medium concentration in Region Of Interest, above-mentioned Region Of Interest is set in respectively photograph from substantially common direction to the head of subject obtained from the roughly the same position for including blood vessel photographed in period different two difference images.Correction coefficient determining section determines correction coefficient in a manner that two profiles measured respectively by profile measurement portion are substantially uniform.DSA image correcting sections according to determine after correction coefficient, correct two difference images in it is at least one.Control unit is controlled, to be shown by display unit based on the information that two at least one difference images are corrected by DSA image correcting sections.
Description
Technical field
Embodiments of the present invention are related to radiographic apparatus and image processing apparatus.
Background technology
In the past, based on X ray CT (Computed Tomography) device in the diagnosis of cerebral infarction, brain perfusion
(Brain Perfusion) analysis is known.In brain perfusion analysis, as by injection contrast agent and carried out photography obtain
The layer image of brain represent that the brain perfusion image of the image of the circulation state of blood is used in the diagnosis of cerebral infarction.
In addition, it in recent years, to the head profile image imaged by X ray CT device, manually or automatically sets
It is fixed that the two cerebral hemispheres on head are divided into two-part boundary line, make an image inversion being partitioned into according to boundary line and be overlapped
On another image, the analysis method for generating the difference image of the image of coincidence is known.According to the analysis method, due to
Emerge lesion region in difference image, therefore, left and right comparison read shadow become easy, can to the blood flow anomalies such as cerebral infarction into
Row checks.However, it is above-mentioned in the prior art, sometimes cannot more preoperative image and postoperative image exactly.
Patent document 1：Japanese Unexamined Patent Publication 2009-153870 publications
The content of the invention
The problem to be solved in the present invention is, provide it is a kind of being capable of more preoperative image and postoperative image exactly
Radiographic apparatus and image processing apparatus.
The radiographic apparatus of embodiment possesses measurement portion, determining section, correction unit and display control unit.Measurement portion
Measurement and the relevant profile of contrast medium concentration in Region Of Interest, above-mentioned Region Of Interest are individually set from roughly the same respectively
It is substantially same comprising blood vessel in photography period that direction photographs to the head of subject different two difference images
Position.Determining section determines correction coefficient in a manner that two profiles measured respectively by above-mentioned measurement portion are substantially uniform.School
Positive portion is at least one in above-mentioned two difference image to correct according to the correction coefficient determined by above-mentioned determining section.Display
Control unit is controlled, to be shown by defined display unit based on being corrected two at least one difference by above-mentioned correction unit
The information of image.Radiographic apparatus according to the above configuration, being capable of more preoperative image and postoperative image exactly.
Description of the drawings
Fig. 1 is the figure of an example of the structure for representing the radiographic apparatus involved by the 1st embodiment.
Fig. 2 is the figure of an example of the structure for representing the ROI configuration parts involved by the 1st embodiment.
Fig. 3 is flow chart the step of representing the processing based on the radiographic apparatus involved by the 1st embodiment.
Fig. 4 is the figure for an example for representing the ROI set as the ROI configuration parts involved by the 1st embodiment.
Fig. 5 is flow chart the step of representing the processing based on the ROI configuration parts involved by the 1st embodiment.
Fig. 6 is for illustrating the figure of an example of the processing based on the blood vessel exploration portion involved by the 1st embodiment.
Fig. 7 A are for illustrating based on the most ebb region determining section involved by the 1st embodiment by each subregion
Integrated value calculating processing figure.
Fig. 7 B are for illustrating based on the most ebb region determining section involved by the 1st embodiment by each subregion
Integrated value calculating processing figure.
Fig. 7 C are for illustrating based on the most ebb region determining section involved by the 1st embodiment by each subregion
Integrated value calculating processing figure.
Fig. 8 is for illustrating at the extraction of the most ebb based on the most ebb region determining section involved by the 1st embodiment
The figure of reason.
Fig. 9 is for illustrating the figure of an example of the processing based on the profile measurement portion involved by the 1st embodiment.
Figure 10 is an example for representing the blood flow check image shown on the display unit involved by the 1st embodiment
Figure.
Figure 11 A are the figures for an example for representing the analysis result shown on the display unit involved by the 1st embodiment.
Figure 11 B are the figures for an example for representing the analysis result shown on the display unit involved by the 1st embodiment.
Figure 11 C are the figures for an example for representing the analysis result shown on the display unit involved by the 1st embodiment.
Figure 12 is the figure for an example for representing the display information shown on the display unit involved by the 1st embodiment.
Specific embodiment
Hereinafter, with reference to the accompanying drawings, radiographic apparatus and image processing apparatus according to the present invention are explained in detail
Embodiment.Wherein, embodiment as shown below does not limit the present invention.
(the 1st embodiment)
Fig. 1 is the figure of an example of the structure for representing the radiographic apparatus 1 involved by the 1st embodiment.Such as Fig. 1
Shown, the radiographic apparatus 1 involved by the 1st embodiment has X-ray mechanism 10 and image processing apparatus 100.X
Radiography mechanism 10 has X-ray tube 11, detector (FPD (Flat Panel Detector)) 12, C-arm 13 and bed
14.C-arm 13 supports X-ray tube 11 and detector 12, by the motor such as propeller for being arranged on pedestal (illustration omitted)
It is rotated at a high speed around subject P like that.
As shown in Figure 1, there is image processing apparatus 100 A/D (Analog/Digital) converter section 21, position offset to determine
Portion 22, position offset correction unit 23, ROI (Region of Interest：Region Of Interest) configuration part 24, profile (profile)
Measurement portion 25, correction coefficient (factor) determining section 26, DSA (Digital Subtraction Angiography) image calibration
Positive portion the 27, the 1st subtracts shadow (subtraction) portion the 28, the 2nd and subtracts shadow portion 29, filtering part 30, affine transform portion 31, LUT (Look Up
Table) 32, video memory 33, control unit 34 and display unit 40.In addition, though do not illustrate, but image processing apparatus
100 receive the various behaviour that operator carries out radiographic apparatus 1 such as with mouse, keyboard, trace ball, positioning device
The input unit of work.
Display unit 40 is shown by treated various images, GUI (the Graphical User of image processing apparatus 100
The various information such as Interface).For example, display unit 40 is CRT (Cathode Ray Tube) displays or liquid crystal display
Deng.A/D converter sections 21 are connected with detector 12, and the analog signal inputted from detector 12 is converted into digital signal, and will be turned
Digital signal after changing is collected image as X-ray and is stored in video memory 33.Video memory 33 stores X-ray and collects
Image.
Position offset determining section 22 determines the position offset of the two DSA images obtained with different time-lapse photographies.Position
Offset correction portion 23 is inclined to correct the position of two DSA images according to the position offset determined by position offset determining section 22
It moves.ROI configuration parts 24 set ROI on DSA images.Specifically, ROI configuration parts 24 are extracted from the lower end of DSA images and are provided
In the range of blood vessel diameter in the blood vessel that is included be maximum blood vessel, and in the blood vessel extracted, by each defined
Region calculates the peak value of the partial product score value of concentration value, and the peak value calculated is set as ROI for minimum region.
Fig. 2 is the figure of an example of the structure for representing the ROI configuration parts 24 involved by the 1st embodiment.Such as Fig. 2 institutes
Show, ROI configuration parts 24 have blood vessel exploration portion 24a, maximum blood vessel determining section 24b and most ebb region determining section 24c.Blood
Pipe exploration portion 24a from the lower end of DSA images away from exploring blood vessel according to DSA images on a certain range of transverse axis.It is maximum
Blood vessel determining section 24b determines the blood vessel of blood vessel diameter maximum according to the blood vessel explored by blood vessel exploration portion 24a.Most ebb area
Domain determining section 24c tracks downwards the maximum blood vessel determined by maximum blood vessel determining section 24b, determines the partial integration of concentration value
The minimum position of the peak value of value.
Back to Fig. 1, profile measurement portion 25 measures profile on the ROI of DSA images is set in.Specifically, profile is surveyed
Amount portion 25 measures respectively to be set at respectively with the relevant profile of contrast medium concentration in ROI, above-mentioned ROI from roughly the same direction
It is big comprising blood vessel in two different DSA images of phase obtained from photographing to the head of subject (photography period)
Cause same position.For example, profile measurement portion 25 measure respectively with comprising the main artery as blood vessel or on treatment do not have it is influential
The relevant profile of contrast medium concentration in the Region Of Interest of capillary.Regional correction coefficient determining section 26 with two profiles substantially
Consistent mode determines correction coefficient.Specifically, correction coefficient determining section 26 to be measured respectively by profile measurement portion 25
Two profiles substantially uniform modes determine correction coefficient.For example, correction coefficient determining section 26 is in two difference images
The substantially uniform mode of two profiles, the offset for determining time until contrast agent reaches, the output with subject
It is at least one in relevant gain and blood flow velocity.
DSA image correcting sections 27 correct two DSA figures according to the correction coefficient determined by correction coefficient determining section 26
At least one DSA images as in.1st subtracts the image before contrast agent is launched in shadow portion 28 and the image progress after contrast agent dispensing
Subtract shadow.2nd subtracts shadow portion 29 carries out subtracting shadow to the DSA images after being corrected by DSA image correcting sections 27 and another DSA image.Filter
Ripple portion 30 carries out high frequency emphasis filtering etc..Affine transform portion 31 carries out amplification, diminution or movement of image etc..LUT32 carries out color
It turns and changes.
Control unit 34 controls radiographic apparatus 1 whole.Specifically, control unit 34 controls the receipts of X-ray collection image
Various processing involved by display of display image in collection, the generation of display image, display unit 40 etc..For example, control unit 34
It is controlled, to be shown by display unit 40 based on correcting two at least one DSA images by DSA image correcting sections 27
Information (for example, difference information etc.).
Radiographic apparatus 1 involved by 1st embodiment can exactly compare existing according to above-mentioned structure
The preoperative image compared and postoperative image are difficult in technology.
For example, in the treatment on head of radiographic apparatus has been used, exist and insert a catheter into narrow, make to set
The air bag around conduit is put to expand to expand the treatment of narrow.The treatment is referred to as intervening.In intervention, when making air bag
During expansion, the small sector-meeting of a part for narrow is flowed with erasing, since the small pieces are sent out in the capillary of brain sometimes
Raw infraction.
In consideration of it, in the postoperative of intervention, the diagnosis blocked is made whether by brain perfusion image.In addition,
In the prior art, such as by the blood flow state of the brain to left and right reading shadow is compared, to carry out the inspection of blood flow anomalies.So
And it is above-mentioned in the prior art, due to not reflecting preoperative state, mistaken diagnosis sometimes has been sent out before the treatment
Raw infraction or the tendency of some dyeing etc..It is treated for example, even infraction chronically occurs without blocking to it
Patient, the infraction that will be had occurred and that also is caused to be diagnosed as the infraction of kainogenesis due to treatment, implements insignificant thrombolytic therapy.
In consideration of it, the radiographic apparatus 1 involved by present embodiment passes through processing described further below, Neng Gouzhun
Really the preoperative image and postoperative image of the operations such as intervention or thrombolytic therapy are compared.Fig. 3 is to represent real based on the 1st
The flow chart for the step of applying the processing of the radiographic apparatus 1 involved by mode.
As shown in figure 3, in radiographic apparatus 1, collected first before the operations such as intervention or thrombolytic therapy start preoperative
DSA images (step S101).Specifically, in radiographic apparatus 1, C-arm 13 is arranged on arbitrary direction,
The frisket image of number frame amount is shot before radiography, afterwards, contrast is continuously shot during contrast agent flows in the blood vessel
(contrast) image.The frisket image of number frame amount before radiography and multiple contrast images are converted by A/D converter sections 21
Into digital signal, and it is stored in video memory 33.
1st subtracts the frisket image that shadow portion 28 reads the number frame amount stored by video memory 33, to read-out frisket figure
As carrying out summation averaging, the few average frisket image of generation noise.Moreover, the 1st subtracts shadow portion 28 and passes through according to multiple contrast figures
As carrying out subtracting shadow (Log additions) respectively to average frisket image, to generate preoperative DSA images.
If subtracting shadow portion 28 by the 1st generates preoperative DSA images, control unit 34 makes the preoperative DSA images of generation exist
40 enterprising Mobile state image display of display unit.Wherein, the preoperative DSA images generated are stored by video memory 33.
Afterwards, such as intervened or the operations such as thrombolytic therapy, when the treatment is completed, radiographic apparatus 1, which passes through, specifies
With preoperative identical data collection program come carry out with preoperative identical data collection, collect postoperative DSA image (steps
S102)。
Moreover, control unit 34 by the postoperative DSA images of generation in 40 enterprising Mobile state image display of display unit.Wherein,
The postoperative DSA images generated are stored by video memory 33.
Then, if received through not shown GUI from user based on preoperative DSA images and postoperative DSA images
The instruction compared of blood flow, then position offset determining section 22 preoperative DSA images and postoperative are read from video memory 33
The respective average frisket image of DSA images, determine position offset (step S103).For example, position offset determining section 22 is logical
Formula as shown below (1) is crossed, to determine the frisket image (M of preoperative DSA imagespre(i, j)) and postoperative DSA images illiteracy
Picture (Mpost(i, j)) position offset.Here, (Δ i, Δ j) represent (M to the CR in formula (1)pre(i, j)) and (Mpost
(i, j)) position offset.In addition, the N in formula (1) represents picture size.
【Formula 1】
22 one side of position offset determining section changes (Δ i, Δ j), while exploring CR with successive approximation type (algorithm)
(Δ i, Δ j) are minimum (Δ i, Δ j).Moreover, position offset determining section 22 by explore (Δ i, Δ j) are to position offset
Correction unit 23 is sent.Wherein, herein for the purpose of simplifying the description, said for the situation of detection to the position offset in the direction of two dimension
It is bright, but also want to the position offset of detection direction of rotation.That is, it is desirable to position offset determining section 22 determines (Δ i, Δ j, Δ θ).
Then, if receiving position offset data from position offset determining section 22, position offset correction unit 23 uses
The position offset data " (Δ i, Δ j) or (Δ i, Δ j, Δ θ) " received, come the DSA images before revision procedure or postoperative DSA
One (step S104) in image.
Then, ROI configuration parts 24 are to DSA image setting ROI (step S105).Specifically, ROI configuration parts 24 use school
The just preoperative DSA images after position offset and a DSA image in postoperative DSA images, DSA images in the preoperative
And the roughly the same position setting ROI of postoperative DSA images.For example, ROI configuration parts 24 are calculated in two difference images extremely
Aortic position in a few image, and will be set as ROI comprising the aortic region calculated.Fig. 4 be represent by
The figure of the example for the ROI that ROI configuration parts 24 involved by 1st embodiment are set.
As shown in figure 4, the roughly the same position of DSA images and postoperative DSA images is set respectively in the preoperative for ROI configuration parts 24
Determine ROI50 and ROI51.Here, it is illustrated below for the details of the processing based on ROI configuration parts 24.Fig. 5 is to represent base
Flow chart in the processing of the ROI configuration parts 24 involved by the 1st embodiment the step of.Wherein, the processing shown in Fig. 5 is equivalent to
Processing in the step S105 of Fig. 3.
As shown in figure 5, in ROI configuration parts 24, blood vessel exploration portion 24a uses the preoperative DSA for correcting position offset
A DSA image in image and postoperative DSA images, from what is included as defined in the exploration of the lower end of DSA images in scope
Blood vessel (step S201).Fig. 6 is for illustrating one of the processing based on the blood vessel exploration portion 24a involved by the 1st embodiment
The figure of example.Here, in (A) of Fig. 6, the feelings that the blood vessel in the DSA images to picture size " N × N " is explored are shown
Condition.
For example, as shown in (A) of Fig. 6, the position of the lower end N/4 of distance DSA images is set to blood vessel by blood vessel exploration portion 24a
The initial position of exploration.Moreover, shown in (B) of blood vessel exploration portion 24a such as Fig. 6, the transverse axis of the position of the N/4 of DSA images is measured
On DSA values profile.Afterwards, DSA values are more than defined threshold value " Th " in the profile being measured to by blood vessel exploration portion 24a
Peak be extracted as blood vessel.For example, as shown in fig. 6, blood vessel exploration portion 24a extractions be more than " Th " the corresponding blood vessel 60 in peak and
Blood vessel 61.Moreover, blood vessel exploration portion 24a by profile and the location information of blood vessel 60 and blood vessel 61 to maximum blood vessel determining section
24b is sent.
Maximum blood vessel determining section 24b is determined in a diameter of maximum of blood vessel medium vessels extracted by blood vessel exploration portion 24a
Blood vessel (step S202).For example, maximum blood vessel determining section 24b is straight to measure blood vessel according to the result of the profile shown in (B) of Fig. 6
Blood vessel diameter is determined as maximum blood vessel by footpath (width at peak) for maximum blood vessel 61.Moreover, maximum blood vessel determining section 24b will
The information of the blood vessel diameter measured and the location information of blood vessel 61 are sent to most ebb region determining section 24c.As it is above-mentioned that
Sample, blood vessel exploration portion 24a explore blood vessel using the position of the lower end N/4 of distance DSA images as the initial position that blood vessel is explored,
Maximum blood vessel determining section 24b is by determining that blood vessel diameter for maximum blood vessel, can extract blood of the main artery as setting ROI
Pipe.
Most ebb region determining section 24c is by each subregion for the blood vessel determined by maximum blood vessel determining section 24b, meter
Calculate the integrated value (step S203) of the concentration value of each time.Specifically, most ebb region determining section 24c from DSA while scheme
Blood vessel center is tracked downwards in the position of the N/4 of the longitudinal axis of picture, while the integration of the concentration value by each time computing section region
Value.Fig. 7 A~Fig. 7 C be for illustrate based on the most ebb region involved by the 1st embodiment determine portion 24c by each portion
The figure of the calculating processing of subregional integrated value.Wherein, in figure 7b, the Fig. 1 being exaggerated to the region 53 of Fig. 7 A is shown.
For example, as shown in Figure 7 A, most ebb region identification portion 24c is while from the lower end N/4's of the longitudinal axis of distance DSA images
Position is downwards tracked the blood vessel center of blood vessel 61, while calculating the integrated value of the concentration value of each subregion.Here,
Such as shown in Figure 7 B, most ebb region determining section 24c is by the square of horizontal " 1.5L (=1.5 × blood vessel diameter) " × vertical " 3 pixel "
Shape is set as subregion, by each calculating integral value in each several part region.In addition, the size of subregion can be by user
Arbitrarily set.
Most ebb region determining section 24c is directed to the DSA images during contrast agent is flowing respectively, by each several part region
Each calculate concentration value integrated value.That is, most ebb region determining section 24c calculates Fig. 7 C institutes by each of subregion
Show the time change of such integrated value.
Then, the integration for each time that most ebb region determining section 24c is calculated by each extraction of subregion
The peak value (step S204) of value.For example, as seen in figure 7 c, most ebb region determining section 24c is extracted by each of subregion
Integrated value becomes the typical value at peak.
Afterwards, most ebb region determining section 24c sets ROI the subregion (step S205) for representing minimum peak.Tool
For body, most ebb region determining section 24c compares the typical value in each several part region, and the subregion for representing minimum is set as
ROI.Fig. 8 is for illustrating the extraction process of the most ebb based on the most ebb region determining section 24c involved by the 1st embodiment
Figure.In fig. 8, transverse axis represents " 0 " to " N/4 " shown in Fig. 7 A, and the longitudinal axis represents integrated value.That is, in fig. 8 it is shown that drawing
The graph of the typical value in each several part region.
For example, as shown in figure 8, most ebb region determining section 24c draws the typical value in each several part region on the graph,
And determine most ebb.Moreover, subregion corresponding with the most ebb determined is set as by most ebb region determining section 24c
ROI。
As described above, most ebb region determining section 24c will represent the part area of minimum in the typical value in each several part region
Domain is set as ROI.That is, the minimum subregion of the time change of integrated value is set as ROI by most ebb region determining section 24c.
Although the main artery blood vessel diameter of setting ROI has almost no change, the concentration value significantly change in actual main artery.This master
Traveling and beam hardening that will be because of blood vessel cause.
For example, when the traveling of blood vessel compared with the direction of illumination of X-ray into 90 degree when, the contrast agent of illuminated X-ray
Amount is only the amount of thickness of blood vessel diameter, and concentration value is minimum.However, when the traveling of blood vessel and the parallel direction of illumination of X-ray,
The amount of the contrast agent of illuminated X-ray becomes the depth amount of blood vessel traveling, and concentration value is maximum.
Therefore, because when the traveling of blood vessel and the parallel irradiation of X-ray, it is impossible to blood flow is assessed exactly, it is desirable that
The traveling of blood vessel is set as ROI compared with the region of irradiation into the angle close to 90 degree of X-ray.In consideration of it, most ebb area
Domain determining section 24c is opposite to extract the traveling of blood vessel by the way that the integrated value of concentration value is set as ROI for minimum subregion
In X-ray irradiation closest to 90 degree of region.
In addition, such as angiosomes Chong Die with thick bone concentration value due to beam hardening reduces.In consideration of it, ROI is set
Determine portion 24 can also according to used in the generation of DSA images frisket image concentration value, the X to calculate each subregion penetrates
Line permeability, using the low region of the X-ray transparent degree calculated as the object of ROI outside.
More than, it is illustrated for the processing based on ROI configuration parts 24.Back to Fig. 3, if right as described above
DSA image setting ROI, then the information of ROI sent to profile measurement portion 25.Profile measurement portion 25 is measured by ROI configuration parts 24
Profile (step S106) in the ROI of setting.
Specifically, profile measurement portion 25 is from DSA images after 33 read-out position offset correction of video memory and another
A DSA images measure the average value of DSA pixel values or the profile of total value in ROI respectively for the DSA images of reading.Fig. 9
It is for illustrating the figure of an example of the processing based on the profile measurement portion 25 involved by the 1st embodiment.
For example, as shown in figure 9, the profile " f (t) " of DSA images and postoperative DSA images before 25 surveying of profile measurement portion
Profile " g (t) ".Here, " Δ T " shown in Fig. 9 represents the time delay until contrast agent reaches ROI.By profile measurement
The information for the profile that portion 25 measures is sent to correction coefficient determining section 26.
If obtaining the information of the profile in ROI, correction coefficient determining section 26 determines preoperative DSA images and postoperative
DSA images ROI in the substantially uniform correction coefficient (step S107) of profile.For example, correction coefficient determining section 26 passes through
Formula (2) as shown below, to determine the profile " g (t) " for profile " f (t) " and the postoperative DSA images for making preoperative DSA images substantially
Consistent correction coefficient.
【Formula 2】
E=| | f (t)-α g T (t- Δs t) } | |2 ...(2)
Correction coefficient determining section 26 is by successive approximation type (algori thm), while changing gain " α ", until radiography
Agent reach until time delay " Δ t " and blood flow velocity " T ", on one side explore " E " become minimum " α ", " Δ t " and
“T”.That is, correction coefficient determining section 26 using for before revision procedure, the different gains of the output of the heart of postoperative patient
" α ", for correct because it is preoperative, it is postoperative launch contrast agent when conduit position offset caused by until contrast agent reach until prolonging
" the Δ t " of slow time and for before revision procedure, the different blood flow velocity " T " of the beats of postoperative patient, to make two
Profile is substantially uniform.Moreover, correction coefficient determining section 26 by " α " that explores, " Δ t " and " T " to DSA image correcting sections 27
It sends.
Wherein, in the above example, carried out for the situation for all correction coefficient that " α ", " Δ t " and " T " is obtained
Explanation.However, it's not limited to that for embodiment.For example, when the heartbeat for assuming patient is stablized, blood flow velocity is substantially certain
When or be set to the situation of " T=1 ".Thereby, it is possible to determining for correction coefficient is performed quickly.
DSA image correcting sections 27 using receive from correction coefficient determining section 26 correction coefficient " α ", " Δ t " and
" T " comes the DSA images before revision procedure or postoperative DSA images (step S108).Here, the general revision procedure of the correction of DSA images
Preceding DSA images, but in the present embodiment, arbitrary DSA images can be corrected.DSA after being corrected by DSA image correcting sections 27
Image is stored in video memory 33.
Wherein, correction coefficient is determined using the DSA images for checking preoperative, postoperative blood flow state herein.However, it is possible to
The DSA images that will not be had an impact to treatment etc. substantially photographed in the phase same time to be used to determine correction coefficient.Tool
For body, such as when carrying out the endovascular treatment of right internal carotid artery, in preoperative, the postoperative DSA images to right internal carotid artery
Preoperative, the postoperative DSA images of the front and rear shooting left internal carotid for the timing photographed.Due to thinking not deposit in the images
In the influence for the treatment of etc., therefore, have the advantages that definite correction coefficient can be stablized.
If corrected by DSA image correcting sections 27 to DSA images, the 2nd subtracts shadow portion 29 from video memory 33
DSA images and another DSA image after correction are read, execution subtracts shadow (step S109).For example, the 2nd subtract shadow portion 29 by with
Lower shown formula (3), carries out subtracting shadow to the DSA images after correction and another DSA image.
【Formula 3】
C(i, j)=α g(i, j){T(t-Δt)}-f(i, j) ...(3)
Here, in formula (3), C (i, j) represents blood flow check image.In addition, in formula (3), f(i,j)(t) represent preoperative
DSA images, g(i,j)(t) postoperative DSA images are represented.In addition, in formula (3), represent when postoperative DSA images are corrected
Formula.As shown in formula (3), the 2nd subtract shadow portion 29 by from the postoperative DSA images after correction difference (difference) it is preoperative
DSA images, to generate blood flow check image.Subtract the blood flow check image (subtraction image) that shadow portion 29 generates by the 2nd to be stored
In video memory 33.
Moreover, subtracting the blood flow check image that shadow portion 29 generates by the 2nd is stored in video memory 33, while to control unit
34 send, and are shown (step S110) with colour on display unit 40 by control unit 34.For example, display unit 40 display based on pair
Correct the perfusion image that two at least one DSA images carry out the result of difference.Figure 10 is represented in the 1st embodiment institute
The figure of one example of the blood flow check image shown on the display unit 40 being related to.
For example, as shown in Figure 10, pixel value is converted into colored (for example, by the fully increased region of blood flow by control unit 34
(region R1 of Figure 10 etc.) is converted into red, region (region R2 of Figure 10 etc.) the conversion au bleu that blood flow is fully reduced, will
Intermediate region (region R3 of Figure 10 etc.) conversion yellowly), and by transformed image display on display unit 40.In addition,
Shown blood flow check image can also be the image of black and white.
As described above, the radiographic apparatus 1 involved by the 1st embodiment along time series by continuously displaying
Blood flow check image shows the dynamic image of colored or black and white, and above-mentioned blood flow check image is preoperative and postoperative to correcting
At least one two DSA images carry out difference obtained from image.Here, the radiodiagnosis involved by the 1st embodiment
Device 1 is in addition to above-mentioned dynamic image, additionally it is possible to show various analysis results.For example, radiographic apparatus 1 can also
Display is represented the graph of the variation of the state of preoperative and postoperative blood flow by each defined region of DSA images.
Figure 11 A~Figure 11 C are an examples for representing the analysis result shown on the display unit involved by the 1st embodiment
The figure of son.In Figure 11 A~Figure 11 C, (A) in figure represent as defined in preoperative and postoperative profile in region, in figure
(B) difference information of preoperative profile and postoperative profile is represented.Wherein, in Figure 11 A~Figure 11 C, show from postoperative wheel
Difference information during wide " g (t) " difference preoperative profile " f (t) ".
For example, control unit 34 be directed to by based on comprising aortic region or comprising treatment do not have influential capillary
The correction coefficient that region is determined corrects the defined region of preoperative and postoperative at least one two DSA images, will
The profile and difference information of preoperative and postoperative DSA images are shown on display unit 40.For example, as shown in (A) of Figure 11 A,
Control unit 34 makes the identical area in the profile " f (t) " and postoperative DSA images in the defined region in preoperative DSA images
The profile " g (t) " in domain is shown on display unit 40.Also, as shown in (B) of Figure 11 A, control unit 34 makes (A) institute from Figure 11 A
Postoperative profile " g (t) " difference for showing graph 71 obtained from preoperative profile " f (t) " is shown on display unit 40.
Here, the defined region of the profile and graph for display diagram, can arbitrarily be specified by operator.
Alternatively, it is also possible to be by the region of defined size cover display image entirety situation.Operator can see in detail as a result,
It examines and shows how the region of profile or graph changes in the preoperative with postoperative.For example, when such profile shown in display Figure 11 A
Or during graph, operator can interpolate that as in corresponding region, the state of blood flow is improving after surgery.
Equally, operator can easily judge to correspond to by the profile and graph shown in Figure 11 B and Figure 11 C
Region how to change with postoperative in the preoperative.For example, work as the profile or graph shown on display unit 40 shown in Figure 11 B
When, operator can interpolate that as several delays has occurred in blood flow after surgery, can suspect may have occurred it is slight narrow
It is narrow.In addition, for example when showing the profile or graph shown in Figure 11 C on display unit 40, operator can interpolate that as art
Blood flow afterwards is deteriorated, and can suspect and may have occurred infraction.
In addition, for having shown from postoperative profile " g (t) " difference preoperative profile " f (t) " in Figure 11 A~Figure 11 C
Obtained from the situation of graph be illustrated, but it's not limited to that for embodiment.That is, can also show from preoperative wheel
Wide " f (t) " difference graph obtained from postoperative profile " g (t) ".
In addition, the radiographic apparatus 1 involved by the 1st embodiment can also show blood flow check image and and blood flow
Corresponding profile in region of check image etc..Figure 12 is to represent the display shown on the display unit involved by the 1st embodiment
The figure of one example of information.As shown in figure 12, if operator specified via input unit in blood flow check image a little or
Region, then control unit 34 shown the preoperative and postoperative profile in specified point or region.In addition, in fig. 12 only
Profile is shown, but can also show the graph of difference.
In this way, in the radiographic apparatus 1 involved by the 1st embodiment, in addition to blood flow check image, moreover it is possible to
Graph of the profile in region or difference as defined in enough displays are each etc..Also, it is examined in the X-ray involved by the 1st embodiment
In disconnected device 1, additionally it is possible to which the blood flow of the difference of the blood flow in region is poor as defined in display expression, represents preoperative and postoperative blood flow
Different diversity factoies.Hereinafter, illustrated using formula (4)~formula (7) for blood flow difference and diversity factor.Wherein, in formula
(4) in~formula (7), " α " and " Δ t " that enumerates in above-mentioned correction coefficient " α ", " Δ t " and " T " is corrected
Situation illustrated for example.
That is, correction coefficient determining section 26 determines the wheel in the ROI of the preoperative DSA images measured by profile measurement portion 25
It is wide " f (t) ", such " α " and " Δ t " as formula (4) with the profile " g (t) " in the ROI of postoperative DSA images.Wherein,
" α " be for before revision procedure, the different gains of the output of the heart of postoperative patient, " Δ t " be used for revision procedure before, it is postoperative
Conduit when launching contrast agent position offset caused by until contrast agent reach until time delay.
【Formula 4】
Moreover, the 2nd subtracts shadow portion 29 by formula (5) as shown below, to the DSA images after correction and another DSA image
It carries out subtracting shadow.Here, in formula (5), COMP (t) represents the difference of DSA images.In addition, in formula (5), f(i,j)(t) art is represented
Preceding DSA images, g(i,j)(t) postoperative DSA images are represented.In addition, in formula (5), show to carry out postoperative DSA images
The formula of timing.
【Formula 5】
COMP (t)=α g(i, j)(t-Δt)-f(i, j)(t) ...(5)
As shown in formula (5), the 2nd subtracts shadow portion 29 by the postoperative DSA images after correction and preoperative DSA images by each picture
Element carries out subtracting shadow.Moreover, control unit 34 subtracts the COMP (t) that shadow portion 29 calculates using by the 2nd, by formula (6) as shown below,
It is poor to calculate blood flow.That is, control unit 34 is according to the profile at each position calculated based on the difference image after correction, meter
Calculate blood flow difference and diversity factor.Here, in formula (6), DifV represents that blood flow is poor.In addition, in formula (6), N represents image
Size.In addition, in formula (6), COMP (n) represents the difference of each pixel.
【Formula 6】
As shown in formula (6), control unit 34 calculates the average difference as blood flow of the difference of each pixel, and is shown in
On display unit 40.In addition, control unit 34 subtracts the COMP (t) that shadow portion 29 calculates using by the 2nd, by formula (7) as shown below,
To calculate diversity factor.Here, in formula (7), DifE represents diversity factor.In addition, in formula (7), N represents picture size.In addition,
In formula (7), COMP (n) represents the difference of each pixel.
【Formula 7】
As shown in formula (7), control unit 34 calculate the difference of each pixel square the average diversity factor as blood flow,
And it is shown on display unit 40.In this way, the radiographic apparatus 1 involved by the 1st embodiment is except blood flow check image, wheel
Outside wide and difference graph, additionally it is possible to calculate blood flow difference and diversity factor and show.Wherein, for blood flow difference with
And diversity factor, use is according to circumstances distinguished by operator.
Here, the formula of an above-mentioned blood flow difference and diversity factor only example, can also pass through others
Formula calculates.For example, it is also possible to it is standardized according to preoperative, postoperative profile.Hereinafter, for according to preoperative, postoperative
The situation that profile is standardized is illustrated using formula (8)~formula (13).Wherein, in formula (8)~formula (13), only show
The right of blood flow difference DifV and diversity factor DifE.
For example, control unit 34 calculates the blood standardized according to preoperative profile by formula (8) as shown below
Difference in flow by formula (9) as shown below, calculates the diversity factor standardized according to preoperative profile.Here, in formula
(8) and in (9), " f (n) " represents preoperative profile.Here, " β " in formula (8) and (9) is used for the feelings that denominator is avoided to be 0
Condition is set to be much smaller than the value of the value of the angiosomes of " f (n) ".
【Formula 8】
【Formula 9】
In addition, for example, control unit 34 calculates by formula (10) as shown below and has carried out standard according to postoperative profile
The blood flow of change is poor, by formula (11) as shown below, calculates the diversity factor standardized according to postoperative profile.
This, in formula (10) and (11), " g (n) " represents preoperative profile, shows that postoperative image is corrected according to correction coefficient
Situation.Here, " β " in formula (10) and (11) is set in order to avoid the situation that denominator is 0 much smaller than " f (n) "
The value of the value of angiosomes.
【Formula 10】
【Formula 11】
In addition, for example, control unit 34 is calculated by formula (12) as shown below based on preoperative and postoperative profile
The blood flow averagely standardized is poor, by formula (13) as shown below, to calculate based on preoperative and postoperative profile
The average diversity factor standardized.Here, in formula (12) and (13), " f (n) " represents preoperative profile, " g
(n) " represent preoperative profile, the situation that postoperative image is corrected according to correction coefficient is shown.Here, formula (12) and
(13) " β " in is used for the situation that denominator is avoided to become 0, is set to be much smaller than the value of the value of the angiosomes of " f (n) ".
【Formula 12】
【Formula 13】
In the above example, enumerate situation about being corrected according to correction coefficient " α " and " Δ t " for an example into
Explanation is gone, but it's not limited to that for embodiment.That is, blood can also be calculated using correction coefficient " α ", " Δ t " and " T "
Difference in flow and diversity factor.At this point, correction coefficient determining section 26 determines the preoperative DSA images measured by profile measurement portion 25
ROI in profile " f (t) " and postoperative DSA images ROI in profile " g (t) " become formula (14) as shown below that
" α " of sample, " Δ t " and " T ".Wherein, " α " be for before revision procedure, the output of the heart of postoperative patient it is different
Gain, " Δ t " be used for revision procedure before, caused by the position offset of postoperative conduit when launching contrast agent until contrast agent arrival
Until time delay.In addition, " T " represent for before revision procedure, the different blood flow velocity of the beats of postoperative patient.
【Formula 14】
As shown in formula (14), the radiographic apparatus 1 involved by the 1st embodiment can correct for before revision procedure, art
The different blood flow velocity of the beats of patient afterwards carry out logistic before, postoperative blood flow difference and diversity factor.Even if as a result,
In the preoperative, in the case that postoperative beats are significantly different, can also be analyzed on the basis of being corrected to them.
As described above, according to the 1st embodiment, profile measurement portion 25 measures related to the contrast medium concentration in ROI respectively
Profile, above-mentioned ROI be set in respectively from substantially common direction photograph to the head of subject obtained from phase (take the photograph
Shadow period) the roughly the same position for including blood vessel in different two DSA images.Moreover, correction coefficient determining section 26 with by
Two profiles that profile measurement portion 25 measures respectively substantially uniform modes determines correction coefficient.DSA image correcting sections 27
According to the correction coefficient determined by correction coefficient determining section 26, correct at least one in two DSA images.Moreover, control
Portion 34 is controlled, so that will be shown based on the information that two at least one DSA images are corrected by DSA image correcting sections 27
Show on display unit 40.So as to which the radiographic apparatus 1 involved by the 1st embodiment, which can generate, makes preoperative DSA images
And the coloring based on contrast agent in postoperative DSA images it is substantially uniform after difference image and shown, can
More preoperative image and postoperative image exactly.
Subtract shadow portion 28 according to the head to subject while from roughly the same side in addition, according to the 1st embodiment, the 1st
Two different X-ray dynamic images of the phase (photograph period) that is obtained to the photography on one side of injection contrast agent, there will be contrast agent
Influence multiple images and there's almost no contrast agent influence image subtraction, to calculate DSA images.Position offset determines
The image of the influence that there's almost no contrast agent of the portion 22 in two X-ray dynamic images determines position offset.Position
Offset correction portion 23 is according to the position offset information determined, at least one carry out position offset school in two DSA images
Just.Profile measurement portion 25 measures the profile of ROI respectively, and above-mentioned ROI is set in respectively corrects two at least one DSA images
In the roughly the same position for including blood vessel.Two wheels of the correction coefficient determining section 26 to be measured respectively by profile measurement portion 25
Wide substantially uniform mode determines correction coefficient.DSA image correcting sections 27 are according to being determined by correction coefficient determining section 26
Correction coefficient, it is at least one in two DSA images to correct.Control unit 34 is controlled, so that will be based on by DSA images
Correction unit 27 corrects the presentation of information of two at least one DSA images on display unit 40.So as to the 1st embodiment institute
The radiographic apparatus 1 being related to can generate make in preoperative DSA images and postoperative DSA images based on contrast agent
Difference image after coloring is substantially uniform simultaneously shown, being capable of more preoperative image and postoperative image exactly.
In addition, according to the 1st embodiment, ROI is comprising the main artery as blood vessel or does not have influential capillary to treatment
The region of blood vessel.So as to which the radiographic apparatus 1 involved by the 1st embodiment can be come true using the various regions in image
Determine correction coefficient.
In addition, according to the 1st embodiment, ROI configuration parts 24 calculate big at least one image in two DSA images
The position of artery, and will be set as ROI comprising the aortic region calculated.So as to which the X involved by the 1st embodiment is penetrated
Ray diagnosis apparatus 1 can be used comprising the main artery that will not be had an impact to treatment and the variation of contrast medium concentration is significantly represented
Region, determine more accurately correction coefficient.
In addition, according to the 1st embodiment, display unit 40 show based on to correct two at least one DSA images into
The perfusion image of the result of row difference.So as to which the radiographic apparatus 1 involved by the 1st embodiment can be shown more accurately
Analysis result.
In addition, according to the 1st embodiment, correction coefficient determining section 26 is with two profiles substantially one in two DSA images
The mode of cause, to determine until at least one in the offset, gain and blood flow velocity of the time until contrast agent reaches.From
And the radiographic apparatus 1 involved by the 1st embodiment can be come using the correction coefficient of the effect significantly change of contrast agent
DSA images are corrected, can more accurately carry out the comparison of preoperative DSA images and postoperative DSA images.
In addition, according to the 1st embodiment, included as defined in ROI configuration parts 24 from the lower end of DSA images extraction in scope
Blood vessel in blood vessel diameter be maximum blood vessel, in the blood vessel extracted, the product of concentration value is calculated by each subregion
The peak value of score value, and the peak value calculated is set as ROI for minimum subregion.So as to involved by the 1st embodiment
The high region of the effect of contrast agent on main artery can be set as ROI by radiographic apparatus 1, can accurately be compared
Compared with reading shadow.
In addition, according to the 1st embodiment, ROI configuration parts 24 according to used in the generation of DSA images frisket image it is dense
Angle value calculates the X-ray transparent degree of each subregion, and using the low subregion of the X-ray transparent degree calculated as ROI
Object outside.So as to which the radiographic apparatus 1 involved by the 1st embodiment can make the process involved by the setting of ROI high
Speedization.
In addition, according to the 1st embodiment, blood flow check image (subtraction image) is formed into dynamic image, but can also
Within dynamic image, worst error is determined by each pixel, is still image by worst error image display.
(the 2nd embodiment)
It is illustrated above in relation to the 1st embodiment, it, can also be with various in addition to the 1st above-mentioned embodiment
Different mode is implemented.
(variation)
In the above-described first embodiment, for data collection program identical when using with collecting preoperative DSA images
Situation to collect postoperative DSA images is illustrated.However, it's not limited to that for embodiment, such as can also be arbitrary
Ground selection make contrast agent dyeing effect variation coefficient, make selected coefficient in the preoperative with it is postoperative consistent.
At this point, for example control unit 34 is controlled, so that X-ray condition, camera angle, FOV (Field Of
View), SID (Source Image Distance), the position of collimator, the position of compensating filter, x-ray focus size,
Radiation quality adjusts wave filter, until at least one and preoperative in the time until contrast agent injects and Photo condition
It is consistent during DSA image photographies, to photograph to postoperative DSA images.
In the above-described first embodiment, automatically the situation of DSA image settings ROI is carried out for ROI configuration parts 24
Explanation.However, it's not limited to that for embodiment, for example, it is also possible to be set by user.
At this point, for example, ROI configuration parts 24 ask user to set ROI.Row are given one example, and ROI configuration parts 24 make display
Portion 40 shows the information for supervising setting ROI.User sets ROI through not shown GUI.At this point, ROI configuration parts 24 are with opposite
In user, close to launch the conduit of contrast agent, away from affected part, and the blood vessel mode parallel with detector face sets ROI.
In the above-described first embodiment, the situation that blood flow check image is generated for radiographic apparatus 1 carries out
Illustrate, but above-mentioned processing can also be performed by image processing apparatus such as work stations.At this point, for example via network and X-ray
The work station of the connections such as diagnostic device or image archive apparatus obtains picture number from radiographic apparatus or image archive apparatus etc.
According to.Moreover, image data acquired by work station use performs above-mentioned processing.
As described above, according to the 1st~2 embodiment, at the radiographic apparatus and image of present embodiment
Managing device being capable of more preoperative image and postoperative image exactly.
Although the description of several embodiments of the invention, but these embodiments are to prompt as an example, not
Intended limitation the scope of the present invention.These embodiments can be implemented in a manner that others are various, not depart from invention master
In the range of purport, various omissions, displacement, change can be carried out.The scope of these embodiments or its deformation with being contained in invention
Or it is the same in purport, it is contained in invention and its equivalent scope of claims record.
Claims (13)
1. a kind of radiographic apparatus, possesses：
Measurement portion, respectively measurement change over time relevant profile, above-mentioned care area with the contrast medium concentration in Region Of Interest
Domain is set in respectively photograph from substantially common direction to the head of subject obtained from photograph period different two
The roughly the same position for including blood vessel in difference image；
Determining section determines correction coefficient in a manner that two profiles measured respectively by above-mentioned measurement portion are substantially uniform；
Correction unit according to the correction coefficient determined by above-mentioned determining section, corrects at least one in above-mentioned two difference image；
And
Control unit is controlled, to be shown by defined display unit based on correcting at least one two by above-mentioned correction unit
The information of a difference image.
2. radiographic apparatus according to claim 1, wherein,
Above-mentioned Region Of Interest is the region for including as the main artery of above-mentioned blood vessel or not having on treatment influential capillary.
3. radiographic apparatus according to claim 2, wherein,
Above-mentioned radiographic apparatus also has calculating part, which calculates at least one figure in above-mentioned two difference image
Above-mentioned aortic position as in,
Above-mentioned Region Of Interest is comprising the above-mentioned aortic region calculated by above-mentioned calculating part.
4. radiographic apparatus according to claim 1, wherein,
Above-mentioned display unit shows the filling of the result to carrying out difference based on above-mentioned at least one two difference images being corrected by
Note image.
5. radiographic apparatus according to claim 1, wherein,
Above-mentioned determining section is in a manner that two profiles in above-mentioned two difference image are substantially uniform, to determine to arrive until contrast agent
It is at least one in the offset of time until reaching, the relevant gain of output with above-mentioned subject and blood flow velocity.
6. radiographic apparatus according to claim 1, wherein,
Above-mentioned two difference image is the image obtained before photography period is treatment and after treatment.
7. radiographic apparatus according to claim 1, wherein,
Configuration part is also equipped with, in the blood vessel included as defined in the configuration part from the lower end of above-mentioned difference image extraction in scope
Blood vessel diameter is maximum blood vessel, and in the blood vessel extracted, the partial product score value of concentration value is calculated by each defined region
Peak value, and the peak value calculated is set as above-mentioned Region Of Interest for minimum region.
8. radiographic apparatus according to claim 7, wherein,
The concentration value of above-mentioned configuration part frisket image according to used in the generation of above-mentioned difference image calculates above-mentioned each regulation
Region X-ray transparent degree, and using the low region of the X-ray transparent degree calculated as the object of above-mentioned Region Of Interest outside.
9. radiographic apparatus according to claim 1, wherein,
Above-mentioned control unit is poor to calculate blood flow according to the profile at each position calculated based on the difference image after correction
And diversity factor.
10. radiographic apparatus according to claim 1, wherein,
Above-mentioned control unit is controlled, so that the position of X-ray condition, camera angle, FOV, SID, collimator, compensation filter
Position, x-ray focus size, radiation quality adjustment wave filter, the time until contrast agent injects and the radiography of ripple device
It is at least one consistent in condition, to photograph to two different difference images of above-mentioned photography period.
11. a kind of radiographic apparatus, possesses：
Difference image calculating part, photographed on one side to the head of subject injection contrast agent from substantially common direction according to one side and
Obtained photography period different two groups of X-ray dynamic images, there will be the influence of contrast agent multiple images and hardly deposit
In the image subtraction of the influence of contrast agent, to calculate difference image；
Position offset determining section, the image of the influence that there's almost no contrast agent in above-mentioned two groups of X-ray dynamic images
To determine the position offset between above-mentioned two groups of X-ray dynamic images；
Position offset correction unit, according to the above-mentioned position offset information determined, to distinguishing from above-mentioned two groups of X-ray dynamic images
At least one carry out position offset correction in two difference images extracted；
Measurement portion, respectively measurement change over time relevant profile, above-mentioned care area with the contrast medium concentration in Region Of Interest
Domain is set in the roughly the same position for including blood vessel that above-mentioned at least one two difference images being corrected are included respectively；
Determining section, in a manner that two profiles measured respectively by above-mentioned measurement portion are substantially uniform, to determine correction coefficient；
Correction unit according to the correction coefficient determined by above-mentioned determining section, corrects at least one in above-mentioned two difference image；
And
Control unit is controlled, to be shown by defined display unit based on correcting at least one two by above-mentioned correction unit
The information of a difference image.
12. a kind of image processing apparatus, possesses：
Measurement portion, respectively measurement change over time relevant profile, above-mentioned care area with the contrast medium concentration in Region Of Interest
Domain is set in two different difference diagrams of the photography period photographed to the head of subject from substantially common direction respectively
The roughly the same position for including blood vessel as in；
Determining section determines correction coefficient in a manner that two profiles measured respectively by above-mentioned measurement portion are substantially uniform；
Correction unit according to the correction coefficient determined by above-mentioned determining section, corrects at least one in above-mentioned two difference image；
And
Control unit is controlled, to be shown by defined display unit based on correcting at least one two by above-mentioned correction unit
The information of a difference image.
13. a kind of image processing apparatus, possesses：
Difference image calculating part, photographed on one side to the head of subject injection contrast agent from substantially common direction according to one side and
Obtained photography period different two groups of X-ray dynamic images, there will be the influence of contrast agent multiple images and hardly deposit
In the image subtraction of the influence of contrast agent, to calculate difference image；
Position offset determining section, the image of the influence that there's almost no contrast agent in above-mentioned two groups of X-ray dynamic images
To determine the position offset between above-mentioned two groups of X-ray dynamic images；
Position offset correction unit, according to the above-mentioned position offset information determined, to distinguishing from above-mentioned two groups of X-ray dynamic images
At least one carry out position offset correction in two difference images extracted；
Measurement portion, respectively measurement change over time relevant profile, above-mentioned care area with the contrast medium concentration in Region Of Interest
Domain is set in the roughly the same position for including blood vessel that above-mentioned at least one two difference images being corrected are included respectively；
Determining section determines correction coefficient in a manner that two profiles measured respectively by above-mentioned measurement portion are substantially uniform；
Correction unit according to the correction coefficient determined by above-mentioned determining section, corrects at least one in above-mentioned two difference image；
And
Control unit is controlled, to be shown by defined display unit based on correcting at least one two by above-mentioned correction unit
The information of a difference image.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013-037984 | 2013-02-27 | ||
| JP2013037984 | 2013-02-27 | ||
| PCT/JP2014/054957 WO2014133104A1 (en) | 2013-02-27 | 2014-02-27 | X-ray diagnostic device, and image processing device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105007818A CN105007818A (en) | 2015-10-28 |
| CN105007818B true CN105007818B (en) | 2018-05-18 |
Family
ID=54380169
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201480010340.0A Active CN105007818B (en) | 2013-02-27 | 2014-02-27 | Radiographic apparatus and image processing apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105007818B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114723844B (en) * | 2022-06-10 | 2022-08-19 | 中国医学科学院北京协和医院 | Method, system and equipment for reconstructing pulsatile perfusion image corrected by saline contrast |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003210456A (en) * | 2002-01-21 | 2003-07-29 | Toshiba Corp | Processor for time series image |
| CN1882950A (en) * | 2003-09-25 | 2006-12-20 | 派昂公司 | System and method for three-dimensional reconstruction of a tubular organ |
| CN1895170A (en) * | 2005-07-15 | 2007-01-17 | 株式会社东芝 | X-ray diagnostic apparatus, imaging angle determination device, program storage medium, and method |
| CN101028194A (en) * | 2006-02-28 | 2007-09-05 | 东芝医疗系统株式会社 | X-ray diagnostic apparatus |
| CN102918563A (en) * | 2010-06-02 | 2013-02-06 | 皇家飞利浦电子股份有限公司 | Automated quantification of intravascular embolization success |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060079746A1 (en) * | 2004-10-11 | 2006-04-13 | Perret Florence M | Apparatus and method for analysis of tissue classes along tubular structures |
-
2014
- 2014-02-27 CN CN201480010340.0A patent/CN105007818B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003210456A (en) * | 2002-01-21 | 2003-07-29 | Toshiba Corp | Processor for time series image |
| CN1882950A (en) * | 2003-09-25 | 2006-12-20 | 派昂公司 | System and method for three-dimensional reconstruction of a tubular organ |
| CN1895170A (en) * | 2005-07-15 | 2007-01-17 | 株式会社东芝 | X-ray diagnostic apparatus, imaging angle determination device, program storage medium, and method |
| CN101028194A (en) * | 2006-02-28 | 2007-09-05 | 东芝医疗系统株式会社 | X-ray diagnostic apparatus |
| CN102918563A (en) * | 2010-06-02 | 2013-02-06 | 皇家飞利浦电子股份有限公司 | Automated quantification of intravascular embolization success |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105007818A (en) | 2015-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11694339B2 (en) | Method and apparatus for determining blood velocity in X-ray angiography images | |
| EP1302163B1 (en) | Apparatus for calculating an index of local blood flows | |
| CN102085100B (en) | Medical image processing apparatus and medical image processing method | |
| KR101774095B1 (en) | Control method and control system | |
| US9569839B2 (en) | Image processing apparatus, method and medical image device | |
| US8929633B2 (en) | Diagnostic X-ray system and method | |
| JP4013128B2 (en) | Blood flow dynamic analysis apparatus, method, and diagnostic imaging apparatus | |
| JP4804039B2 (en) | Blood flow dynamic analysis apparatus, X-ray CT apparatus, MRI apparatus, and blood flow dynamic analysis program | |
| US20090238434A1 (en) | Method for reproducing the spatial orientation of an immobilized subject in a multi-modal imaging system | |
| US20120308106A1 (en) | Method and apparatus for measuring activity of a tracer | |
| US20160051218A1 (en) | X-ray diagnostic apparatus and image processing apparatus | |
| US10043268B2 (en) | Medical image processing apparatus and method to generate and display third parameters based on first and second images | |
| CN105007818B (en) | Radiographic apparatus and image processing apparatus | |
| WO2012153231A9 (en) | Myocardial clusterification and orientation | |
| JP5358855B2 (en) | Apparatus and method for analyzing examination image by nuclear cardiology examination method | |
| EP2501287B1 (en) | Functional imaging | |
| EP2285282B1 (en) | Spatial orientation method of an immobilized subject | |
| JP2013510621A5 (en) | ||
| US11145094B2 (en) | Image reconstruction apparatus and image reconstruction method | |
| JP2010022667A (en) | Cerebral blood flow analyzer | |
| van der Geest et al. | Advances in PET Imaging of Large Vessel Vasculitis: An Update and Future Trends | |
| JP2023176825A (en) | Dynamic image analysis device and program | |
| Anoosha | Approval Sheet |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20160711 Address after: Japan Tochigi Tian Yuan City Applicant after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Applicant before: Toshiba Corp Applicant before: Toshiba Medical System Co., Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |