CN105007818A - X-ray diagnostic device, and image processing device - Google Patents

Abstract

An X-ray diagnostic device (1) according to an embodiment of the present invention is provided with: a profile measurement unit (25); a correction-factor identification unit (26); a DSA image correction unit (27); and a controller (34). The profile measurement unit (25) respectively measures profiles in 2 difference images of a head of a subject imaged from substantially identical directions at different imaging times, said profiles being related to contrast-medium concentration in blood-vessel-containing regions of interest respectively set in substantially identical positions. The correction-factor identification unit (26) identifies a correction factor such that two profiles respectively measured by the profile measurement unit (25) substantially correspond. The DSA image correction unit (27) corrects at least one of the 2 difference images on the basis of the identified correction factor. The controller (34) implements control such that information based upon the 2 difference images, at least one of which having been corrected by the DSA image correction unit (27), is displayed on a display (40).

Description

Radiographic apparatus and image processing apparatus

Technical field

Embodiments of the present invention relate 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) analyze be known.In brain perfusion analysis, as by injecting contrast agent and the layer image of the brain having carried out photographing and obtained to be used to the diagnosis of cerebral infarction to the brain perfusion image of the image representing the circulation state of blood.

In addition, in recent years, to the head profile image obtained of being made a video recording by X ray CT device, manually or automatically set and the two cerebral hemispheres of head is divided into two-part boundary line, make an image inversion being partitioned into according to boundary line and overlap on another image, the analytical method generating the difference image of the image overlapped is known.According to this analytical method, owing to appearing lesion region in one's mind in difference image, therefore, the shadow of relatively reading of left and right becomes easy, can check blood flow anomalies such as cerebral infarction.But, in above-mentioned prior art, sometimes can not image more preoperative exactly and postoperative image.

Patent documentation 1: Japanese Unexamined Patent Publication 2009-153870 publication

Summary of the invention

The problem to be solved in the present invention is, provide a kind of can the radiographic apparatus of image more preoperative exactly and postoperative image and image processing apparatus.

The radiographic apparatus of embodiment possesses measurement section, determination portion, correction unit and display control unit.Measurement section measures the profile relevant to the contrast concentration in Region Of Interest respectively, and above-mentioned Region Of Interest is individually set the roughly same position comprising blood vessel in photography period of photographing from the head of roughly the same direction to subject different two difference images.Determination portion in measured respectively by above-mentioned measurement section two modes that profile is roughly consistent to determine correction coefficient.Correction unit, according to the correction coefficient determined by above-mentioned determination portion, corrects at least one in above-mentioned two difference images.Display control unit controls, to show the information based on two difference images being corrected at least one by above-mentioned correction unit by the display part specified.

Detailed description of the invention

Below, with reference to the accompanying drawings, radiographic apparatus involved in the present invention and the embodiment of image processing apparatus is explained.Wherein, embodiment shown below does not limit the present invention.

(the 1st embodiment)

Fig. 1 is the figure of an example of the structure of the radiographic apparatus 1 represented involved by the 1st embodiment.As shown in Figure 1, the radiographic apparatus 1 involved by the 1st embodiment has X-ray mechanism 10 and image processing apparatus 100.X-ray 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 that is arranged on base (omit diagram) as propeller around subject P high speed rotating.

As shown in Figure 1, image processing apparatus 100 has A/D (Analog/Digital) converter section 21, position skew determination portion 22, offset correction portion, position 23, ROI (Region ofInterest: Region Of Interest) configuration part 24, profile (profile) measurement section 25, correction coefficient (factor) determination portion 26, DSA (Digital Subtraction Angiography) image correcting section 27, 1st subtracts shadow (subtraction) portion 28, 2nd subtracts shadow portion 29, filtering part 30, affine converter section 31, LUT (Look Up Table) 32, image storage 33, control part 34, and display part 40.In addition, although do not illustrate, image processing apparatus 100 such as has mouse, keyboard, trace ball, positioning equipment etc. and accepts the input part of operator to the various operations that radiographic apparatus 1 carries out.

Display part 40 shows the various information such as the various images after being processed by image processing apparatus 100, GUI (Graphical User Interface).Such as, display part 40 is CRT (Cathode Ray Tube) display or liquid crystal display etc.A/D converter section 21 is connected with detector 12, converts the analogue signal inputted from detector 12 to digital signal, and the digital signal after conversion is collected Image Saving in image storage 33 as X-ray.Image storage 33 storing X collection of rays image.

Position skew determination portion 22 determines the position skew of two the DSA images obtained with different time-lapse photographys.Offset correction portion, position 23 offsets according to the position determined by position skew determination portion 22, corrects the position skew of two DSA images.ROI configuration part 24 sets ROI on DSA image.Specifically, the blood vessel diameter that ROI configuration part 24 is extracted the blood vessel comprised in the scope of regulation from the lower end of DSA image is maximum blood vessel, and in the blood vessel extracted, by the peak value of the partial product score value of the region calculating concentration value of each regulation, and be that minimum region is set as ROI by the peak value calculated.

Fig. 2 is the figure of an example of the structure of the ROI configuration part 24 represented involved by the 1st embodiment.As shown in Figure 2, ROI configuration part 24 has blood vessel exploration portion 24a, maximum blood vessel determination portion 24b and most ebb region determination portion 24c.Blood vessel exploration portion 24a explores blood vessel away from the transverse axis of certain limit according to DSA image in the lower end from DSA image.Maximum blood vessel determination portion 24b determines according to the blood vessel explored by blood vessel exploration portion 24a the blood vessel that blood vessel diameter is maximum.Most ebb region determination portion 24c follows the trail of the maximum blood vessel determined by maximum blood vessel determination portion 24b downwards, determines the position that the peak value of the partial product score value of concentration value is minimum.

Turn back to Fig. 1, profile is measured on the ROI being set in DSA image by profile measurement portion 25.Specifically, the profile relevant to the contrast concentration in ROI is measured respectively by profile measurement portion 25, and above-mentioned ROI is set at respectively and photographs from the head of roughly the same direction to subject and the roughly the same position comprising blood vessel different two the DSA images of the phase (photography period) that obtains.Such as, profile measurement portion 25 is measured and the large artery trunks comprised as blood vessel or the profile relevant to the contrast concentration for the treatment of in the Region Of Interest not having influential blood capillary respectively.Regional correction coefficient determination portion 26 determines correction coefficient in the mode that two profiles are roughly consistent.Specifically, correction coefficient determination portion 26 in measured respectively by profile measurement portion 25 two modes that profile is roughly consistent to determine correction coefficient.Such as, correction coefficient determination portion 26, in the roughly consistent mode of two profiles in two difference images, is determined until at least one in the skew of time of arriving of contrast agent, the gain relevant to the stroke volume of subject and blood flow rate.

DSA image correcting section 27 corrects at least one the DSA image in two DSA images according to the correction coefficient determined by correction coefficient determination portion 26.1st subtract the 28 pairs of contrast agent in shadow portion throw in before image and contrast agent throw in after image subtract shadow.2nd subtracts shadow portion 29 subtracts shadow to the DSA image after being corrected by DSA image correcting section 27 and another DSA image.Filtering part 30 carries out high frequency emphasis filtering etc.Affine converter section 31 carry out image amplification, to reduce or mobile etc.LUT32 carries out tone conversion.

It is overall that control part 34 controls radiographic apparatus 1.Specifically, control part 34 controls X-ray and collects the involved various process such as the display of the display image in the collection of image, the display generation of image, display part 40.Such as, control part 34 controls, to show the information (such as, difference information etc.) based on two the DSA images being corrected at least one by DSA image correcting section 27 by display part 40.

Radiographic apparatus 1 involved by 1st embodiment, according to above-mentioned structure, can be difficult to the preoperative image that compares and postoperative image exactly more in the prior art.

Such as, in the treatment of head employing radiographic apparatus, exist and conduit is inserted into narrow, make the airbag inflation of the surrounding being arranged on conduit to expand the treatment of narrow.This treatment is called as intervention.In intervention, when making airbag inflation, the little sector-meeting of a part for narrow is flowed to erase, because these small pieces block sometimes in the blood capillary of brain.

Given this, postoperative what get involved, carry out by brain perfusion image the diagnosis that whether there occurs infraction.In addition, in the prior art, such as compared by the blood flow state of the brain to left and right and read shadow, carry out the inspection of blood flow anomalies.But, in above-mentioned prior art, owing to not reflecting preoperative state, therefore, the tendency etc. of the infraction that occurred before the treatment of mistaken diagnosis or some dyeing sometimes.Such as, even chronically infraction occur and do not need to block to it patient treated, also cause the infraction occurred to be diagnosed as the de novo infraction because for the treatment of, implement insignificant thrombolytic therapy.

Given this, the radiographic apparatus 1 involved by present embodiment, can exactly to get involved or the preoperative image of the operation such as thrombolytic therapy and postoperative image compare by the process of following detailed description.Fig. 3 is the flow chart of the step of the process represented based on the radiographic apparatus 1 involved by the 1st embodiment.

As shown in Figure 3, in radiographic apparatus 1, first before the operation such as intervention or thrombolytic therapy starts, collect preoperative DSA image (step S101).Specifically, in radiographic apparatus 1, C-arm 13 is arranged on arbitrary direction, before radiography, takes the frisket image of number frame amount, afterwards, during contrast agent flows in the blood vessel, take contrast (contrast) image continuously.Frisket image and multiple contrast image of the several frame amount before radiography are converted into digital signal by A/D converter section 21, and are stored in image storage 33.

1st subtracts the frisket image that shadow portion 28 reads the several frame amount stored by image storage 33, carries out summation averaging to read-out frisket image, generates the average frisket image that noise is few.And the 1st subtracts shadow portion 28 by subtracting shadow (Log addition) according to multiple contrast image respectively to average frisket image, generates preoperative DSA image.

If subtract shadow portion 28 by the 1st to generate preoperative DSA image, then control part 34 makes the preoperative DSA image of generation in the enterprising Mobile state image display of display part 40.Wherein, the preoperative DSA image being comprised of image memorizer 33 generated stores.

Afterwards, such as carry out getting involved or the operation such as thrombolytic therapy, when the treatment is completed, radiographic apparatus 1 carries out and preoperative identical data collection with preoperative identical data collection program by specifying, and collects postoperative DSA image (step S102).

And the postoperative DSA image generated shows at the enterprising Mobile state image of display part 40 by control part 34.Wherein, the postoperative DSA image being comprised of image memorizer 33 generated stores.

Then, if receive from user the instruction that the blood flow based on preoperative DSA image and postoperative DSA image compares via not shown GUI, then skew determination portion 22 in position reads preoperative DSA image and postoperative DSA image average frisket image separately from image storage 33, determines position offset (step S103).Such as, position skew determination portion 22, by formula (1) shown below, determines the frisket image (M of preoperative DSA image _pre(i, j)) and the frisket image (M of postoperative DSA image _post(i, j)) position offset.At this, the CR (Δ i, Δ j) in formula (1) represents (M _pre(i, j)) and (M _post(i, j)) position offset.In addition, the N in formula (1) represents picture size.

[formula 1]

$CR(\mathrm{\Δ}i,\mathrm{\Δ}j)=\underset{i=1}{\overset{N}{\Σ}}\underset{H}{\overset{N}{\Σ}}{\{M_{post}(i,j)-M_{pre}(i+\mathrm{\Δ}i,j+\mathrm{\Δ}j)\}}^2\mathrm{...}\left(1\right)$

Position skew determination portion 22, while change (Δ i, Δ j) with successive approximation type (algorithm), explores CR (Δ i, Δ j) for minimum (Δ i, Δ j).And skew determination portion 22 in position sends exploring (Δ i, the Δ j) that arrive to offset correction portion, position 23.Wherein, in this case simplified illustration, is described for detecting situation about offset to the position in the direction of two dimension, but also wishes that the position detecting direction of rotation offsets.That is, desired location skew determination portion 22 is determined (Δ i, Δ j, Δ θ).

Then, if receive position offset data from position skew determination portion 22, then position offset correction portion 23 uses the position offset data " (Δ i; Δ j) or (Δ i; Δ j; Δ θ) " received, and carrys out (step S104) in the DSA image before revision procedure or postoperative DSA image.

Then, 24 pairs, ROI configuration part DSA image setting ROI (step S105).Specifically, ROI configuration part 24 uses the preoperative DSA image after correcting position skew and a DSA image in postoperative DSA image, the roughly the same position setting ROI of DSA image in the preoperative and postoperative DSA image.Such as, ROI configuration part 24 calculates the aortic position at least one image in two difference images, and is set as ROI by comprising the aortic region calculated.Fig. 4 is the figure of the example representing the ROI that the ROI configuration part 24 involved by the 1st embodiment sets.

As shown in Figure 4, the roughly the same position of ROI configuration part 24 DSA image and postoperative DSA image in the preoperative sets ROI50 and ROI51 respectively.At this, the details below for the process based on ROI configuration part 24 is described.Fig. 5 is the flow chart of the step of the process represented based on the ROI configuration part 24 involved by the 1st embodiment.Wherein, the process shown in Fig. 5 is equivalent to the process in the step S105 of Fig. 3.

As shown in Figure 5, in ROI configuration part 24, blood vessel exploration portion 24a uses a DSA image in the preoperative DSA image and postoperative DSA image correcting position skew, explores from the lower end of DSA image the blood vessel (step S201) comprised in the scope of regulation.Fig. 6 is the figure of an example for illustration of the process based on the blood vessel exploration portion 24a involved by the 1st embodiment.At this, in (A) of Fig. 6, situation about exploring the blood vessel in the DSA image of picture size " N × N " is shown.

Such as, as shown in (A) of Fig. 6, the position of the lower end N/4 of distance DSA image is set to the initial position that blood vessel is explored by blood vessel exploration portion 24a.And blood vessel exploration portion 24a, as shown in (B) of Fig. 6, measures the profile of the DSA value on the transverse axis of the position of the N/4 of DSA image.Afterwards, blood vessel exploration portion 24a is in the profile be measured to, and peak DSA value being exceeded the threshold value " Th " of regulation is extracted as blood vessel.Such as, as shown in Figure 6, blood vessel exploration portion 24a extracts the blood vessel 60 corresponding with the peak exceeding " Th " and blood vessel 61.And the positional information of profile and blood vessel 60 and blood vessel 61 sends to maximum blood vessel determination portion 24b by blood vessel exploration portion 24a.

Maximum blood vessel determination portion 24b determines that at the blood vessel medium vessels diameter extracted by blood vessel exploration portion 24a be maximum blood vessel (step S202).Such as, the result of the profile of maximum blood vessel determination portion 24b according to (B) of Fig. 6 measures blood vessel diameter (width at peak), is that maximum blood vessel 61 is defined as maximum blood vessel by blood vessel diameter.And the positional information of the information of the blood vessel diameter measured and blood vessel 61 sends to most ebb region determination portion 24c by maximum blood vessel determination portion 24b.As described above, the initial position that the position of the lower end N/4 of distance DSA image is explored as blood vessel is explored blood vessel by blood vessel exploration portion 24a, maximum blood vessel determination portion 24b, by determining that blood vessel diameter is maximum blood vessel, can extract the blood vessel of large artery trunks as setting ROI.

Most ebb region determination portion 24c, by each subregion of the blood vessel determined by maximum blood vessel determination portion 24b, calculates the integrated value (step S203) of the concentration value of each time.Specifically, most ebb region determination portion 24c follows the trail of blood vessel center downwards from the position of the N/4 of the longitudinal axis of DSA image, while the integrated value of concentration value by each time computing section region.Fig. 7 A ~ Fig. 7 C is the figure of the computing of the integrated value by each subregion for illustration of most qualification bonding part, the ebb region 24c involved by the 1st embodiment.Wherein, in figure 7b, the Fig. 1 be exaggerated the region 53 of Fig. 7 A is shown.

Such as, as shown in Figure 7 A, most qualification portion, ebb region 24c follows the trail of the blood vessel center of blood vessel 61 downwards from the position of the lower end N/4 of the longitudinal axis of distance DSA image, while calculate the integrated value of the concentration value of each subregion.At this, such as shown in Figure 7 B, the rectangle of horizontal " 1.5L (=1.5 × blood vessel diameter) " × vertical " 3 pixel " is set as subregion, by each calculating integral value in each several part region by most ebb region determination portion 24c.In addition, the size of subregion at random can be set by user.

DSA image during most ebb region determination portion 24c flows for contrast agent respectively, by the integrated value of each calculating concentration value in each several part region.That is, most ebb region determination portion 24c by subregion each calculation chart 7C shown in the time variations of such integrated value.

Then, most ebb region determination portion 24c extracts the peak value (step S204) of the integrated value of each time calculated by each of subregion.Such as, as seen in figure 7 c, most ebb region determination portion 24c extracts the representative value that integrated value becomes peak by subregion each.

Afterwards, most ebb region determination portion 24c represents the subregion (step S205) of minimum peak to ROI setting.Specifically, most ebb region determination portion 24c compares the representative value in each several part region, will represent that the subregion of minimum is set as ROI.Fig. 8 is the figure of the extraction process of most ebb for illustration of the most ebb region determination portion 24c involved by the 1st embodiment.In fig. 8, transverse axis represents that " 0 " shown in Fig. 7 A is to " N/4 ", and the longitudinal axis represents integrated value.That is, in fig. 8, the curve chart of the representative value depicting each several part region is shown.

Such as, as shown in Figure 8, the representative value in each several part region is drawn on the graph by most ebb region determination portion 24c, and determines most ebb.And the subregion corresponding with the most ebb determined is set as ROI by most ebb region determination portion 24c.

As mentioned above, most ebb region determination portion 24c will represent in the representative value in each several part region that the subregion of minimum is set as ROI.That is, subregion minimum for the time variations of integrated value is set as ROI by most ebb region determination portion 24c.Although the large artery trunks blood vessel diameter of setting ROI has almost no change, the concentration value in the large artery trunks of reality significantly changes.This is mainly because the traveling of blood vessel and beam hardening cause.

Such as, when the traveling of blood vessel becomes 90 degree relative to the direction of illumination of X-ray, the amount of the contrast agent of illuminated X-ray is only the amount of thickness of blood vessel diameter, and concentration value is minimum.But, when the traveling of blood vessel and the direction of illumination of X-ray parallel time, the amount of the contrast agent of illuminated X-ray becomes the depth amount of blood vessel traveling, and concentration value is maximum.

Therefore, due to when the traveling of blood vessel and the irradiation of X-ray parallel time, can not assess blood flow exactly, so wish the traveling of blood vessel to become the region of the angle close to 90 degree to be set as ROI relative to the irradiation of X-ray.Given this, most ebb region determination portion 24c by being that minimum subregion is set as ROI by the integrated value of concentration value, the traveling of extracting blood vessel relative to the irradiation of X-ray closest to the region of 90 degree.

In addition, because of beam hardening, concentration value reduces such as overlapping with thick skeleton angiosomes.Given this, the concentration value of the frisket image that ROI configuration part 24 also can use according to the generation of DSA image, calculates the X-ray transparent degree of each subregion, using region low for the X-ray transparent degree that calculates as outside the object of ROI.

Above, be illustrated for the process based on ROI configuration part 24.Turn back to Fig. 3, if as described above to DSA image setting ROI, then the information of ROI is sent by profile measurement portion 25.The profile (step S106) in the ROI that set by ROI configuration part 24 is measured in profile measurement portion 25.

Specifically, the meansigma methods of DSA pixel value in ROI or the profile of total value, from the DSA image after image storage 33 read-out position offset correction and another DSA image, are measured respectively for the DSA image read by profile measurement portion 25.Fig. 9 is the figure of an example for illustration of the process based on the profile measurement portion 25 involved by the 1st embodiment.

Such as, as shown in Figure 9, the profile " f (t) " of DSA image and the profile " g (t) " of postoperative DSA image before profile measurement portion 25 surveying.At this, " the Δ T " shown in Fig. 9 represents until contrast agent arrives the time delay of ROI.The information of the profile measured by profile measurement portion 25 is sent by correction coefficient determination portion 26.

If obtain the information of the profile in ROI, then correction coefficient determination portion 26 determines the correction coefficient (step S107) that profile in the ROI of preoperative DSA image and postoperative DSA image is roughly consistent.Such as, correction coefficient determination portion 26, by formula (2) shown below, determines to make the correction coefficient that the profile " f (t) " of preoperative DSA image is roughly consistent with the profile " g (t) " of postoperative DSA image.

[formula 2]

E＝||f(t)-αg{T(t-Δt)}|| ²...(2)

Correction coefficient determination portion 26 is by successive approximation type (algorithm), while change gain " α " until contrast agent arrive time delay " Δ t " and blood flow rate " T ", exploration " E " become minimum " α ", " Δ t " and " T ".Namely, correction coefficient determination portion 26 use the different gain " α " of the stroke volume of the heart being used for patient before revision procedure, postoperative, for the position skew corrected because of conduit during preoperative, postoperative input contrast agent cause until contrast agent " the Δ t " of time delay that arrive and for before revision procedure, the different blood flow rate " T " of the beats of postoperative patient, make two profiles roughly consistent.And " α ", " Δ t " and " T " that explore sends to DSA image correcting section 27 by correction coefficient determination portion 26.

Wherein, in above-mentioned example, for obtaining " α ", the situation of all correction coefficients of " Δ t " and " T " is illustrated.But embodiment is not limited thereto.Such as, when the heart beating of supposition patient is stablized, a blood flow rate roughly timing also can be the situation being set to " T=1 ".Thereby, it is possible to perform the determination of correction coefficient rapidly.

DSA image correcting section 27 use receive from correction coefficient determination portion 26 correction coefficient " α ", " Δ t " and " T ", come the DSA image before revision procedure or postoperative DSA image (step S108).At this, the DSA image before the general revision procedure of correction of DSA image, but in the present embodiment, any DSA image can be corrected.Image storage 33 is stored in by the DSA image after DSA image correcting section 27 corrects.

Wherein, the DSA image of preoperative, the postoperative blood flow state of inspection is used to determine correction coefficient at this.But, also can use roughly the photography of phase same time obtain to the DSA images that can not have an impact such as treatments to determine correction coefficient.Such as, specifically, when carrying out the endovascular treatment of RICA, at preoperative, the postoperative DSA image of the front and back of the timing of photographing to preoperative, the postoperative DSA image of RICA shooting left internal carotid.Owing to thinking the impact that there is not treatment etc. in the images, therefore, there is the advantage can stablized and determine correction coefficient.

If corrected by DSA image correcting section 27 pairs of DSA images, then the 2nd subtracts shadow portion 29 from the DSA image after image storage 33 reading correction and another DSA image, performs and subtracts shadow (step S109).Such as, the 2nd subtracts shadow portion 29 by formula (3) shown below, subtracts shadow to the DSA image after correction and another DSA image.

[formula 3]

C(i，j)＝αg _(i，j){T(t-Δt)}-f _(i，j)(t) ...(3)

At this, in formula (3), C (i, j) represents blood flow check image.In addition, in formula (3), f _{(i, j)}t () represents preoperative DSA image, g _{(i, j)}t () represents postoperative DSA image.In addition, in formula (3), represent formula when postoperative DSA image is corrected.Shown in (3), the 2nd subtracts shadow portion 29 passes through the DSA image that difference (difference) is preoperative from the postoperative DSA image after correction, generates blood flow check image.Subtract by the 2nd the blood flow check image (subtraction image) that shadow portion 29 generates to be stored in image storage 33.

And, subtract by the 2nd the blood flow check image that shadow portion 29 generates and be stored in image storage 33, send to control part 34 simultaneously, carried out showing (step S110) with colour by control part 34 on display part 40.Such as, display part 40 shows the perfusion image based on carrying out the result of difference to two the DSA images correcting at least one.Figure 10 is the figure of an example of the blood flow check image representing display on the display part 40 involved by the 1st embodiment.

Such as, as shown in Figure 10, control part 34 converts pixel value to colour (such as, the region (the region R1 etc. of Figure 10) fully increased by blood flow converts redness to, the region (the region R2 etc. of Figure 10) fully reduced by blood flow converts blueness to, convert the region (the region R3 etc. of Figure 10) of centre to yellow), and the image after conversion is presented on display part 40.In addition, shown blood flow check image also can be the image of black and white.

As described above, radiographic apparatus 1 involved by 1st embodiment is by showing along time series the dynamic image that blood flow check image comes display color or black and white continuously, and above-mentioned blood flow check image carries out difference and the image obtained to correcting at least one two preoperative and postoperative DSA images.At this, the radiographic apparatus 1 involved by the 1st embodiment, except above-mentioned dynamic image, can also show various analysis result.Such as, radiographic apparatus 1 can also show the curve chart of the change of the state by the preoperative and postoperative blood flow of the region representation of each regulation of DSA image.

Figure 11 A ~ Figure 11 C is the figure of the example representing the analysis result shown on the display part involved by the 1st embodiment.In Figure 11 A ~ Figure 11 C, (A) in figure represents the preoperative and postoperative profile in the region of regulation, and (B) in figure represents the difference information of preoperative profile and postoperative profile.Wherein, in Figure 11 A ~ Figure 11 C, difference information time preoperative profile " f (t) " is shown from postoperative profile " g (t) " difference.

Such as, control part 34 treats for by based on comprising aortic region or comprising the region that the correction coefficient not having the region of influential blood capillary to determine corrects the regulation of at least one two preoperative and postoperative DSA images, the profile of preoperative and postoperative DSA image and difference information is presented on display part 40.Such as, as shown in (A) of Figure 11 A, control part 34 makes the profile " f (t) " in the region of the regulation in preoperative DSA image be presented on display part 40 with the profile " g (t) " in the identical region in postoperative DSA image.And, as shown in (B) of Figure 11 A, control part 34 has made from postoperative profile " g (the t) " difference shown in (A) of Figure 11 A preoperative profile " f (t) " and the curve Figure 71 obtained is presented at display part 40.

At this, for the region of regulation of the illustrated profile of display and curve chart, at random can be specified by operator.In addition, also can be the situation being covered display integral image by the region of the size specified.Thus, operator can observe the region demonstrating profile or curve chart in detail and how to change with postoperative in the preoperative.Such as, when showing such profile or curve chart shown in Figure 11 A, operator can be judged as that, in the region of correspondence, the state of blood flow is improved after surgery.

Equally, operator, by Figure 11 B and the profile shown in Figure 11 C and curve chart, easily can judge how corresponding region changes with postoperative in the preoperative.Such as, when showing the profile shown in Figure 11 B or curve chart on display part 40, operator can be judged as there occurs some delays in blood flow after surgery, can suspect may have occurred slight narrow.In addition, such as, when showing the profile shown in Figure 11 C or curve chart on display part 40, operator can be judged as that postoperative blood flow is deteriorated, and can suspect and may have occurred infraction.

In addition, for shown in Figure 11 A ~ Figure 11 C, from postoperative profile " g (t) " difference, preoperative profile " f (t) " and the situation of the curve chart obtained are illustrated, but embodiment is not limited thereto.That is, postoperative profile " g (t) " also can be shown from preoperative profile " f (t) " difference and the curve chart obtained.

In addition, the radiographic apparatus 1 involved by the 1st embodiment can also show blood flow check image and the profile etc. corresponding with the region of blood flow check image.Figure 12 is the figure of the example representing the display information shown on the display part involved by the 1st embodiment.As shown in figure 12, if operator specifies a little or region in blood flow check image via input part, then control part 34 makes the preoperative and postoperative profile in specified point or region show.In addition, only profile is shown in fig. 12, but the curve chart of difference can also be shown.

Like this, in the radiographic apparatus 1 involved by the 1st embodiment, except blood flow check image, the profile in the region of each regulation or the curve chart etc. of difference can also be shown.Further, in the radiographic apparatus 1 involved by the 1st embodiment, the blood flow of difference that can also show the blood flow in the region representing regulation is poor, the different diversity factor that represents preoperative from postoperative blood flow.Below, formula (4) ~ formula (7) is used to be described for blood flow difference and diversity factor.Wherein, in formula (4) ~ formula (7), enumerating the situation carrying out correcting according to " α " and " Δ t " in above-mentioned correction coefficient " α ", " Δ t " and " T " is that an example is described.

That is, correction coefficient determination portion 26 is determined profile " f (t) " in the ROI of the preoperative DSA image measured by profile measurement portion 25, is become formula (4) such " α " and " Δ t " with the profile " g (t) " in the ROI of postoperative DSA image.Wherein, " α " be for before revision procedure, the different gain of the stroke volume of the heart of postoperative patient, " Δ t " for before revision procedure, the postoperative position skew because of conduit when throwing in contrast agent cause until contrast agent time delay of arriving.

[formula 4]

And the 2nd subtracts shadow portion 29 by formula (5) shown below, shadow is subtracted to the DSA image after correction and another DSA image.At this, in formula (5), COMP (t) represents the difference of DSA image.In addition, in formula (5), f _{(i, j)}t () represents preoperative DSA image, g _{(i, j)}t () represents postoperative DSA image.In addition, in formula (5), illustrate and the formula of timing has been carried out to postoperative DSA image.

[formula 5]

COMP(t)＝αg _(i，j)(t-Δt)-f _(i，j)(t)...(5)

Shown in (5), the 2nd subtracts shadow portion 29 is undertaken subtracting shadow by each pixel by the postoperative DSA image after correction and preoperative DSA image.And control part 34 uses and subtracts by the 2nd the COMP (t) that shadow portion 29 calculates, and by formula (6) shown below, calculates blood flow difference.That is, control part 34 is according to the profile at each position calculated based on the difference image after correction, calculates blood flow difference and diversity factor.At this, in formula (6), DifV represents that blood flow is poor.In addition, in formula (6), N represents picture size.In addition, in formula (6), COMP (n) represents the difference of each pixel.

[formula 6]

$DifV=\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}COMP\left(n\right)\mathrm{...}\left(6\right)$

Shown in (6), control part 34 calculates the average difference as blood flow of the difference of each pixel, and is presented on display part 40.In addition, control part 34 uses and subtracts by the 2nd the COMP (t) that shadow portion 29 calculates, and by formula (7) shown below, carrys out calculated difference degree.At this, 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]

$DifE=\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}COMP{\left(n\right)}^2\mathrm{...}\left(7\right)$

Shown in (7), control part 34 calculate the difference of each pixel square the average diversity factor as blood flow, and to be presented on display part 40.Like this, the radiographic apparatus 1 involved by the 1st embodiment, except the curve chart of blood flow check image, profile and difference, can also calculate blood flow difference and diversity factor showing.Wherein, for blood flow difference and diversity factor, according to circumstances use is distinguished by operator.

At this, above-mentioned blood flow difference and a formula only example of diversity factor, also can be calculated by other formula.Such as, also standardization can be carried out according to preoperative, postoperative profile.Below, carry out standardized situation for according to preoperative, postoperative profile, use formula (8) ~ formula (13) to be described.Wherein, in formula (8) ~ formula (13), the right of blood flow difference DifV and diversity factor DifE is only shown.

Such as, control part 34, by formula (8) shown below, calculates that to have carried out standardized blood flow according to preoperative profile poor, by formula (9) shown below, calculates and has carried out standardized diversity factor according to preoperative profile.At this, in formula (8) and (9), " f (n) " represents preoperative profile.At this, " β " in formula (8) and (9) for avoid denominator be 0 situation, be set to the value of the value of the angiosomes much smaller than " f (n) ".

[formula 8]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\frac{COMP\left(n\right)}{f\left(n\right)+\mathrm{\β}}\mathrm{...}\left(8\right)$

[formula 9]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}{\left\{\frac{COMP\left(n\right)}{f\left(n\right)+\mathrm{\β}}\right\}}^2\mathrm{...}\left(9\right)$

In addition, such as, control part 34, by formula (10) shown below, calculates that to have carried out standardized blood flow according to postoperative profile poor, by formula (11) shown below, calculate and carried out standardized diversity factor according to postoperative profile.At this, in formula (10) and (11), " g (n) " represents preoperative profile, illustrates that postoperative image carries out situation about correcting according to correction coefficient.At this, " β " in formula (10) and (11) is the situation of 0 in order to avoid denominator, is set to the value of the value of the angiosomes much smaller than " f (n) ".

[formula 10]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\frac{COMP\left(n\right)}{\mathrm{\α}g(n-\mathrm{\Δ}t)+\mathrm{\β}}\mathrm{...}\left(10\right)$

[formula 11]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}{\left\{\frac{COMP\left(n\right)}{\mathrm{\α}g(n-\mathrm{\Δ}t)+\mathrm{\β}}\right\}}^2\mathrm{...}\left(11\right)$

In addition, such as, control part 34 is by formula (12) shown below, calculate based on preoperative and postoperative profile on average to have carried out standardized blood flow poor, by formula (13) shown below, calculate and on average carried out standardized diversity factor based on preoperative and postoperative profile.At this, in formula (12) and (13), " f (n) " represents preoperative profile, and " g (n) " represents preoperative profile, illustrates that postoperative image carries out situation about correcting according to correction coefficient.At this, " β " in formula (12) and (13), for the situation avoiding denominator to become 0, is set to the value of the value of the angiosomes much smaller than " f (n) ".

[formula 12]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\frac{COMP\left(n\right)}{\frac{f\left(n\right)+\mathrm{\α}g(n-\mathrm{\Δ}t)}{2}+\mathrm{\β}}\mathrm{...}\left(12\right)$

[formula 13]

$\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}{\left\{\frac{COMP\left(n\right)}{\frac{f\left(n\right)+\mathrm{\α}g(n-\mathrm{\Δ}t)}{2}+\mathrm{\β}}\right\}}^2\mathrm{...}\left(13\right)$

In above-mentioned example, enumerating the situation carrying out correcting according to correction coefficient " α " and " Δ t " is that an example is illustrated, but embodiment is not limited thereto.That is, also can use correction coefficient " α ", " Δ t " and " T " calculate blood flow difference and diversity factor.Now, the profile " g (t) " that correction coefficient determination portion 26 is determined in the ROI of profile " f (t) " in the ROI of the preoperative DSA image measured by profile measurement portion 25 and postoperative DSA image becomes formula (14) shown below such " α ", " Δ t " and " T ".Wherein, " α " be for before revision procedure, the different gain of the stroke volume of the heart of postoperative patient, " Δ t " for before revision procedure, postoperative because of the position skew of conduit when throwing in contrast agent cause until contrast agent time delay of arriving.In addition, " T " represents the different blood flow rate of the beats being used for patient before revision procedure, postoperative.

[formula 14]

Shown in (14), the radiographic apparatus 1 involved by the 1st embodiment can correct for before revision procedure, the different blood flow rate of the beats of postoperative patient carry out logistic before, postoperative blood flow difference and diversity factor.Thus, though in the preoperative, postoperative beats significantly different when, also can analyze on the basis that they are corrected.

As mentioned above, according to the 1st embodiment, the profile relevant to the contrast concentration in ROI is measured respectively by profile measurement portion 25, and above-mentioned ROI is set in respectively and photographs from the head of roughly the same direction to subject and the roughly the same position comprising blood vessel different two the DSA images of the phase (photography period) that obtains.And, correction coefficient determination portion 26 in measured respectively by profile measurement portion 25 two modes that profile is roughly consistent to determine correction coefficient.DSA image correcting section 27, according to the correction coefficient determined by correction coefficient determination portion 26, corrects at least one in two DSA images.And control part 34 controls, with make by based on the information displaying of two DSA images being corrected at least one by DSA image correcting section 27 on display part 40.Thus, radiographic apparatus 1 involved by 1st embodiment can generate the difference image after making the Color based on contrast agent in preoperative DSA image and postoperative DSA image roughly consistent and show, can image more preoperative exactly and postoperative image.

In addition, according to the 1st embodiment, 1st subtracts shadow portion 28 to photograph different two the X-ray dynamic images of the phase (photographing period) that obtains on one side according to injecting contrast agent to the head of subject while from roughly the same direction, by there is multiple image of the impact of contrast agent and there is the image subtraction of impact of contrast agent hardly, calculate DSA image.According to the image that there is the impact of contrast agent hardly in two X-ray dynamic images, position skew determination portion 22 determines that position offsets.Offset correction portion, position 23, according to the position offset information determined, carries out position offset correction at least one in two DSA images.The profile of ROI is measured in profile measurement portion 25 respectively, and above-mentioned ROI is set in the roughly the same position comprising blood vessel in two the DSA images correcting at least one respectively.Correction coefficient determination portion 26 in measured respectively by profile measurement portion 25 two modes that profile is roughly consistent to determine correction coefficient.DSA image correcting section 27, according to the correction coefficient determined by correction coefficient determination portion 26, corrects at least one in two DSA images.Control part 34 controls, with make by based on the information displaying of two DSA images being corrected at least one by DSA image correcting section 27 on display part 40.Thus, radiographic apparatus 1 involved by 1st embodiment can generate the difference image after making the Color based on contrast agent in preoperative DSA image and postoperative DSA image roughly consistent and show, can image more preoperative exactly and postoperative image.

In addition, according to the 1st embodiment, ROI comprises as the large artery trunks of blood vessel or to treating the region not having influential blood capillary.Thus the radiographic apparatus 1 involved by the 1st embodiment can use various regions in image to determine correction coefficient.

In addition, according to the 1st embodiment, ROI configuration part 24 calculates the aortic position at least one image in two DSA images, and is set as ROI by comprising the aortic region calculated.Thus, the radiographic apparatus 1 involved by the 1st embodiment can use comprise to treatment can not have an impact and the change of contrast concentration by the aortic region significantly represented, determine correction coefficient more accurately.

In addition, according to the 1st embodiment, display part 40 shows the perfusion image based on carrying out the result of difference to two the DSA images correcting at least one.Thus the radiographic apparatus 1 involved by the 1st embodiment can show analysis result more accurately.

In addition, according to the 1st embodiment, correction coefficient determination portion 26, in the roughly consistent mode of two profiles in two DSA images, is determined until at least one in the skew of time of arriving of contrast agent, gain and blood flow rate.Thus the correction coefficient that the radiographic apparatus 1 involved by the 1st embodiment can use the effect of contrast agent significantly to change, to correct DSA image, can carry out comparing of preoperative DSA image and postoperative DSA image more accurately.

In addition, according to the 1st embodiment, the blood vessel diameter that ROI configuration part 24 is extracted the blood vessel comprised in the scope of regulation from the lower end of DSA image is maximum blood vessel, in the blood vessel extracted, by the peak value of the integrated value of each subregion calculating concentration value, and be that minimum subregion is set as ROI by the peak value calculated.Thus region high for the effect of the contrast agent on large artery trunks can be set as ROI by the radiographic apparatus 1 involved by the 1st embodiment, can compare accurately and read shadow.

In addition, according to the 1st embodiment, the concentration value of the frisket image that ROI configuration part 24 uses according to the generation of DSA image, calculates the X-ray transparent degree of each subregion, and using subregion low for the X-ray transparent degree that calculates as outside the object of ROI.Thus the radiographic apparatus 1 involved by the 1st embodiment can make the operation high speed involved by setting of ROI.

In addition, according to the 1st embodiment, blood flow check image (subtraction image) is formed into dynamic image, but also within dynamic image, by each pixel determination maximum error, maximum error image can be shown as still image.

(the 2nd embodiment)

Be illustrated for the 1st embodiment above, except the 1st above-mentioned embodiment, also can implement in a variety of ways.

(variation)

In the above-described first embodiment, be illustrated for the situation using the data collection program identical with during the DSA image that collection is preoperative to collect postoperative DSA image.But embodiment is not limited thereto, the coefficient such as also can at random selecting the effect of the dyeing making contrast agent to change, makes selected coefficient consistent with postoperative in the preoperative.

Now, such as control part 34 controls, to make the position of X-ray condition, camera angle, FOV (Field Of View), SID (Source Image Distance), collimator, the position of compensating filter, x-ray focus size, radiation quality adjustment wave filter until at least one of contrast agent time of injecting and Photo condition is consistent with during preoperative DSA image photography, postoperative DSA image is photographed.

In the above-described first embodiment, automatically the situation of DSA image setting ROI is illustrated for ROI configuration part 24.But embodiment is not limited thereto, such as, also can be set by user.

Now, such as, ROI configuration part 24 asks user to set ROI.Row are given one example, and ROI configuration part 24 makes display part 40 show the information of supervising setting ROI.User sets ROI via not shown GUI.Now, ROI configuration part 24 with relative to user, close to throw in contrast agent conduit, away from affected part, and the blood vessel mode parallel with detector face sets ROI.

In the above-described first embodiment, the situation generating blood flow check image for radiographic apparatus 1 is illustrated, but above-mentioned process also can be performed by image processing apparatus such as work stations.Now, the work station be such as connected with radiographic apparatus or image archive apparatus etc. via network obtains view data from radiographic apparatus or image archive apparatus etc.And the view data acquired by work station uses performs above-mentioned process.

As described above, according to 1st ~ 2 embodiments, the radiographic apparatus of present embodiment and image processing apparatus can image more preoperative exactly and postoperative images.

Although the description of several embodiment of the present invention, but these embodiments are pointed out as an example, is not intended to limit scope of the present invention.These embodiments can be implemented in other various modes, within a range not departing from the gist of the invention, can carry out various omission, displacement, change.These embodiments or its distortion be contained in scope of invention or purport the same, be contained in claims record invention and equivalent scope in.

Accompanying drawing explanation

Fig. 1 is the figure of an example of the structure of the radiographic apparatus represented involved by the 1st embodiment.

Fig. 2 is the figure of an example of the structure of the ROI configuration part represented involved by the 1st embodiment.

Fig. 3 is the flow chart of the step of the process represented based on the radiographic apparatus involved by the 1st embodiment.

Fig. 4 is the figure of an example of the ROI of the ROI configuration part setting represented involved by the 1st embodiment.

Fig. 5 is the flow chart of the step of the process represented based on the ROI configuration part involved by the 1st embodiment.

Fig. 6 is the figure of an example for illustration of the process based on the blood vessel exploration portion involved by the 1st embodiment.

Fig. 7 A is the figure of the computing of the integrated value by each subregion for illustration of the most ebb region determination portion involved by the 1st embodiment.

Fig. 7 B is the figure of the computing of the integrated value by each subregion for illustration of the most ebb region determination portion involved by the 1st embodiment.

Fig. 7 C is the figure of the computing of the integrated value by each subregion for illustration of the most ebb region determination portion involved by the 1st embodiment.

Fig. 8 is the figure of the extraction process of most ebb for illustration of the most ebb region determination portion involved by the 1st embodiment.

Fig. 9 is the figure of an example for illustration of the process based on the profile measurement portion involved by the 1st embodiment.

Figure 10 is the figure of the example representing the blood flow check image shown on the display part involved by the 1st embodiment.

Figure 11 A is the figure of the example representing the analysis result shown on the display part involved by the 1st embodiment.

Figure 11 B is the figure of the example representing the analysis result shown on the display part involved by the 1st embodiment.

Figure 11 C is the figure of the example representing the analysis result shown on the display part involved by the 1st embodiment.

Figure 12 is the figure of the example representing the display information shown on the display part involved by the 1st embodiment.

Claims (13)

1. a radiographic apparatus, possesses:

Measurement section, measure the profile relevant to the contrast concentration in Region Of Interest respectively, above-mentioned Region Of Interest is set in respectively photographs from the head of roughly the same direction to subject and the roughly the same position comprising blood vessel photography period of obtaining different two difference images;

Determination portion, in measured respectively by above-mentioned measurement section two modes that profile is roughly consistent to determine correction coefficient;

Correction unit, according to the correction coefficient determined by above-mentioned determination portion, corrects at least one in above-mentioned two difference images; And

Control part, controls, to show the information based on two difference images being corrected at least one by above-mentioned correction unit by the display part specified.

2. radiographic apparatus according to claim 1, wherein,

Above-mentioned Region Of Interest comprises as the large artery trunks of above-mentioned blood vessel or to treating the region not having influential blood capillary.

3. radiographic apparatus according to claim 2, wherein,

Above-mentioned radiographic apparatus also has calculating part, and this calculating part calculates the above-mentioned aortic position at least one image in above-mentioned two difference images,

Above-mentioned Region Of Interest comprises the above-mentioned aortic region calculated by above-mentioned calculating part.

4. radiodiagnosis dress according to claim 1, wherein,

Above-mentioned display part display is to the perfusion image carrying out the result of difference based at least one two difference image be corrected above-mentioned.

5. radiographic apparatus according to claim 1, wherein,

Above-mentioned determination portion, in the roughly consistent mode of two profiles in above-mentioned two difference images, is determined until at least one in the skew of time of arriving of contrast agent, the gain relevant to the stroke volume of above-mentioned subject and blood flow rate.

6. radiographic apparatus according to claim 1, wherein,

Above-mentioned two difference images are period the image obtained before treatment and after treating in photography.

7. radiographic apparatus according to claim 1, wherein,

Also possesses configuration part, the blood vessel diameter that this configuration part is extracted the blood vessel comprised in the scope of regulation from the lower end of above-mentioned difference image is maximum blood vessel, in the blood vessel extracted, by the peak value of the partial product score value of the region calculating concentration value of each regulation, and be that minimum region is set as above-mentioned Region Of Interest by the peak value calculated.

8. radiographic apparatus according to claim 7, wherein,

The concentration value of the frisket image that above-mentioned configuration part uses according to the generation of above-mentioned difference image, calculates the X-ray transparent degree in the region of above-mentioned each regulation, and using region low for the X-ray transparent degree that calculates as outside the object of above-mentioned Region Of Interest.

9. radiographic apparatus according to claim 1, wherein,

Above-mentioned control part, according to the profile at each position calculated based on the difference image after correction, calculates blood flow difference and diversity factor.

10. radiographic apparatus according to claim 1, wherein,

Above-mentioned control part controls, to make the position of X-ray condition, camera angle, FOV, SID, collimator, the position of compensating filter, x-ray focus size, radiation quality adjustment wave filter until at least one of contrast agent time of injecting and Photo condition is consistent, two difference images different to above-mentioned photography period are photographed.

11. 1 kinds of radiographic apparatus, possess:

Difference image calculating part, according to while inject from roughly the same direction to the head of subject two X-ray dynamic images that contrast agent is photographed and photography period of obtaining is different on one side, by there is multiple image of the impact of contrast agent and there is the image subtraction of impact of contrast agent hardly, calculate difference image;

According to the image that there is the impact of contrast agent hardly in above-mentioned two X-ray dynamic images, position skew determination portion, determines that position offsets;

Offset correction portion, position, according to the above-mentioned position offset information determined, carries out position offset correction at least one in two X-ray difference images;

Measurement section, measure the profile relevant to the contrast concentration in Region Of Interest respectively, above-mentioned Region Of Interest is set in the roughly the same position comprising blood vessel that at least one two X-ray difference image be corrected above-mentioned comprise respectively;

Determination portion, in measured respectively by above-mentioned measurement section two modes that profile is roughly consistent, determines correction coefficient;

Correction unit, according to the correction coefficient determined by above-mentioned determination portion, corrects at least one in above-mentioned two difference images; And

Control part, controls, to show the information based on two difference images being corrected at least one by above-mentioned correction unit by the display part specified.

12. 1 kinds of image processing apparatus, possess:

Measurement section, measure the profile relevant to the contrast concentration in Region Of Interest respectively, above-mentioned Region Of Interest is set in the roughly the same position comprising blood vessel different two difference images in photography period of photographing from the head of roughly the same direction to subject respectively;

Determination portion, in measured respectively by above-mentioned measurement section two modes that profile is roughly consistent to determine correction coefficient;

Correction unit, according to the correction coefficient determined by above-mentioned determination portion, corrects at least one in above-mentioned two difference images; And

Control part, controls, to show the information based on two difference images being corrected at least one by above-mentioned correction unit by the display part specified.

13. 1 kinds of image processing apparatus, possess:

Difference image calculating part, according to while inject from roughly the same direction to the head of subject two X-ray dynamic images that contrast agent is photographed and photography period of obtaining is different on one side, by there is multiple image of the impact of contrast agent and there is the image subtraction of impact of contrast agent hardly, calculate difference image;

According to the image that there is the impact of contrast agent hardly in above-mentioned two X-ray dynamic images, position skew determination portion, determines that position offsets;

Offset correction portion, position, according to the above-mentioned position offset information determined, carries out position offset correction at least one in two X-ray difference images;

Measurement section, measure the profile relevant to the contrast concentration in Region Of Interest respectively, above-mentioned Region Of Interest is set in the roughly the same position comprising blood vessel that at least one two X-ray difference image be corrected above-mentioned comprise respectively;

Determination portion, in measured respectively by above-mentioned measurement section two modes that profile is roughly consistent to determine correction coefficient;

Correction unit, according to the correction coefficient determined by above-mentioned determination portion, corrects at least one in above-mentioned two difference images; And

Control part, controls, to show the information based on two difference images being corrected at least one by above-mentioned correction unit by the display part specified.