TWI524198B - A morbidity assessment device, a morbidity assessment method, and a morbidity assessment program - Google Patents

Description

Susceptibility assessment device, attack rate assessment method and attack rate assessment program

The present invention relates to an attack rate evaluation device, an attack rate evaluation method, and an attack rate evaluation program.

Recently, there have been many techniques for analyzing medical images by a computer, performing or supporting the determination of whether or not there is an abnormality, and the specificity of a hypothetical disease name.

The image display device of Patent Document 1 acquires a feature amount indicating that the candidate region is more abnormal than the candidate region of the abnormal image in the diagnostic image, and displays the feature amount similarly among the medical images stored in advance, and displays them in a comparable manner. Diagnose portraits and matching medical portraits.

The diagnostic support device of Patent Document 2 stores a reference image that has been diagnosed in the database. Then, the feature amount on the image of the lesion position of the diagnostic image is extracted, and the extracted feature amount is compared with the feature amount of each of the reference images to calculate the similarity. In addition, a plural reference image with a high degree of similarity is selected, and the average value of the similarity is displayed for each disease name, that is, the probability of the disease name.

The diagnostic support device of Patent Document 3 is memory containing the detected The inspection history information of the information of the lesion part. Then, in accordance with the inspection history information, when processing the data obtained by the subject, the size of the lesion to be treated is changed. In addition, according to the examination history information, the probability of occurrence of the lesion is calculated, and the treatment method of the lesion is changed according to the probability of occurrence.

The image collection device of Patent Document 4 collects an image and detects the type of the abnormality and the degree of progress. Next, when the degree of progress is within the predetermined range, the output must be additionally checked, and the inspection condition for the additional inspection is determined.

The outpatient selection device of Patent Document 5 acquires the virtual patient information of the patient candidate based on the disease information including the type of the disease of the patient candidate, the time of occurrence, and the place of occurrence, and determines whether or not the patient candidate needs to perform the outpatient service.

The study of Non-Patent Document 1 compares MRI (Magnetic Resonance Imaging) images of brains of patients with mild cognitive impairment in a time series. Next, it was reported that brain atrophy actually occurred in all patients who were surveyed.

[Previous Technical Literature] [Patent Literature]

[Patent Document 1] Japanese Laid-Open Patent Publication No. 2006-34585 (Request Item 1 and paragraph 0031)

[Patent Document 2] Japanese Patent No. 4021179 (Request Item 7, Paragraph 0026)

[Patent Document 3] Japanese Patent No. 5,519, 924 (Requests 1, 3) Wait)

[Patent Document 4] Japanese Patent No. 3495327 (Request Item 1 and the like)

[Patent Document 5] Japanese Laid-Open Patent Publication No. 2010-113477 (Request Item 1 and the like)

[Non-patent literature]

[Non-Patent Document 1] Nakamura Kenji, and the other three, "Study on Time-Differential Portraits for Dementia Support", [online], Heisei 21 (Earth 2009), the first data worker Learning and Information Management Related Society, Publication No. E3-4, [Search on April 11th, 2013], Internet (URL: http://db-event.jpn.org/deim2009/proceedings/ Files/E3-4.pdf)

In the devices of Patent Documents 1, 2, and 3, when the two medical images are compared, the "features" of the respective images are calculated, and the "similarity" is calculated from the two "features" (see, for example, the paragraph of Patent Document 2). 0026). In this case, the "feature amount" is a multi-dimensional vector in which the complex physical quantity obtained by the medical image is used as a component. In terms of the plural component, it is assumed that there are many people including volume, brightness, and circularity. In order to obtain all the physical quantities shown above, the image quality of the medical image is quite high, which is a necessary condition. However, in general, group health diagnoses, health checkups, and short-term health diagnoses that have not been completed for one year are conducted for a large number of people. many. In the case of the above, there are restrictions on the technical side and the cost side, and it is difficult to obtain a high-quality image with a large amount of information.

The image collection device of Patent Document 4 detects the degree of progress of an abnormality, but this detection is based on the premise that there is an example of precision diagnosis in which an abnormality diagnosis has occurred in the past. Therefore, the device of Patent Document 4 is not applicable even when it is not expected to obtain a high-quality image.

The outpatient selection device of Patent Document 5 targets a patient who has some abnormality but actually has an outpatient history, and is not a certain diagnosis for a patient who is newly diagnosed. The study of Non-Patent Document 1 utilizes an image of a patient who has been diagnosed as a mild cognitive disorder in the past. Neither of them is a healthy person who is included in a group of health care, health checkups, and short-term health checkups that are less than one year old.

Therefore, when it is necessary to image a large number of medical images in a short period of time, such as group health diagnosis, etc., it is determined that there is an abnormality based on a high-quality image, and each disease name is used. Calculate the possibility of suffering.

The attack rate evaluation device of the present invention is characterized in that: the memory unit and the control unit are provided, and the memory unit stores a difference between two images at different imaging points when the image is captured by the memory in association with the body part. And the diagnosis information of the disease name, the control unit accepts the type of the feature quantity input by the user and the part to be diagnosed; In the type of the image to be diagnosed, the difference between the two images at the time of the imaging of the part to be diagnosed is calculated, and the diagnosis target is searched for as the search key, and the diagnostic information is retrieved. The disease name that matches the record is calculated based on the ratio of the number of coincident records to the number of records having the part, and the calculated attack rate is displayed for each disease name obtained.

Other means are described in the form for carrying out the invention.

According to the present invention, when it is necessary to image a large number of medical images in a short period of time, such as a group health diagnosis, it is possible to determine whether or not there is an abnormality based on a high-quality image, and calculate the problem for each disease name. The possibility.

1‧‧‧ attack rate assessment device

2‧‧‧ Terminal devices

3‧‧‧ camera

4‧‧‧Network

11‧‧‧Central Control Unit (Control Department)

12‧‧‧ Input device

13‧‧‧Output device

14‧‧‧Main memory device (memory department)

15‧‧‧Auxiliary memory device (memory department)

16‧‧‧Communication device

21‧‧‧ Relative Position Correction Department

22‧‧‧Differential image calculation unit

23‧‧‧Change Rate Calculation Department

24‧‧‧About the presence or absence of the judgment department

25‧‧‧Occurrence rate calculation department

26‧‧‧Import and Control Department

31‧‧‧Portrait Management Information

32‧‧‧Abnormality information

33‧‧‧Physician Diagnostic Information

34‧‧‧ medical portrait

35‧‧‧Type list

51a, 51b‧‧‧Diagnosis results display screen

101‧‧‧ medical portrait

102‧‧‧ medical portrait

103‧‧‧ Medical Image ID

104‧‧‧ Corresponding information

111‧‧‧The medical portrait of the past

111a‧‧‧The heart of the skull

111b‧‧‧ the apex of the skull

112‧‧‧ Current medical portrait

112a, 112b‧‧ points

113‧‧‧ Medical portrait

114‧‧‧ Past diagnostic examples

121‧‧‧The medical portrait of the past

122‧‧‧ Medical portraits of the present and the past

123‧‧‧Extracting the difference

124‧‧‧ Calculate the rate of change of the difference from the baseline

125‧‧‧ Time series chart for change rate

126‧‧‧Applicable time series chart of rate of change to plural "types"

127‧‧‧Complete information on memory with 1 or multiple diseases and the number of diagnoses

131‧‧‧ Patient ID column

132‧‧‧ Patient Name Column

133‧‧‧ part name column

134a‧‧‧"Previous" column

134b‧‧‧Last time bar

135a‧‧‧ "This time" column

135b‧‧‧Time point bar

136‧‧‧"Differential" column

137‧‧‧Results table

137a‧‧‧Characteristics column

137b‧‧‧change rate column

137c‧‧‧Changes column

137d‧‧‧Results column

138‧‧‧ nuisance rate table

138a‧‧‧ disease name column

138b‧‧‧ Suffering rate column

138c‧‧‧determination column

139‧‧‧"Determining Results Login" button

140a‧‧‧"Baseline" column

140b‧‧‧"Previous previous time" column

140c‧‧‧"Previous" column

140d‧‧‧ "This time" column

141a‧‧‧ benchmark time point

141b‧‧‧ ago time ago

141c‧‧‧Last time bar

141d‧‧‧Time point bar

142a, 142b, 142c, 143‧‧‧ columns

201‧‧‧ Patient ID column

202‧‧‧ patient name column

203‧‧‧Portrait ID column

204‧‧‧Camera ID column

205‧‧‧Photographing time bar

206‧‧‧ part name column

207‧‧‧Portrait file name column

211‧‧‧ part name column

212‧‧‧Characteristics column

213‧‧‧Priority column

214‧‧‧Normal rate of change range

215‧‧‧Normal variation range column

216‧‧‧Unit column

217‧‧‧Change period column

221‧‧‧ Patient ID column

222‧‧‧Physician ID column

223‧‧‧Diagnosis Image ID column

224‧‧‧ part name column

225‧‧‧Change rate column

226‧‧‧Changes column

227‧‧‧ disease name column

228‧‧‧Date of diagnosis

231, 232, 233, 234‧‧‧ time series chart

241‧‧‧Type column

247‧‧‧ symbol shift bar

243‧‧‧ increase and decrease tendency column

244‧‧‧Example column of time series chart

Figure 1 is a diagram illustrating the prior art.

Fig. 2 is a view for explaining the features of the embodiment.

Fig. 3 is a view for explaining the features of the embodiment.

Fig. 4 is a view for explaining the features of the embodiment.

Fig. 5 is a view for explaining the features of the embodiment.

Fig. 6 is a configuration diagram of the attack rate evaluation device.

Fig. 7 is a view showing an example of portrait management information.

Fig. 8 is a view showing an example of abnormality determination information.

Fig. 9 is a view showing an example of physician diagnosis information.

Figure 10 is a diagram illustrating the type of time series chart.

Fig. 11 is a flow chart showing the first processing procedure.

Fig. 12 is a flow chart showing the second processing procedure.

Fig. 13 is a view showing an example of a diagnosis result display screen.

Fig. 14 is a view showing an example of a diagnosis result display screen.

Hereinafter, a mode for carrying out the present invention (referred to as "this embodiment") will be described with reference to the drawings and the like.

(existing technology)

An example of the prior art will be described with reference to FIG. This example is among the aforementioned documents and is closest to "Patent Document 2" (paragraph 0026, etc.). The current medical image 102 of a certain part (in the example of FIG. 1) of the patient to be diagnosed is acquired. On the other hand, the past medical image 101 of the part of the other patient is accumulated. The medical image 101 is a medical image ID 103 in which a part of a majority of patients is photographed and a medical image of a specific medical image is attached. On the other hand, there is correspondence information 104 in which the name of the disease diagnosed by the doctor is stored in association with the medical image ID.

In the first device of the prior art, the feature amount of the medical image 102 is extracted for each of n types such as "area" and "lightness". In the second aspect of the apparatus, a "feature amount vector Q" having a plurality of feature quantities q _i (i = 1, 2, ..., n) as described above is created. The number attached to the components of each vector is, for example, "1" indicates the feature amount of "area", and "2" indicates the feature amount of "lightness". This device is the third, and the "feature amount vector P" is created in the same manner for all of the past medical images 101. This device is the fourth, and the similarity S is calculated. The degree of similarity S is obtained by inputting the feature quantity vector Q and any one of the feature quantity vectors P, and the scalar quantity calculated by the following formula 1.

S=W ^t (E-(QP))/|W| (Expression 1)

Here, "E" is a vector in which all n components are "1". "W" is a vector having weights of each feature amount as n components. "|W|" is the sum (quantity) of the components of the vector "W". "W ^t " is the "transpose matrix" of the vector "W". If the vector "E-(QP)" is multiplied by "W ^t ", the scalar quantity is calculated. In general, the components of the vectors "P" and "Q" are normalized in the range of 0 to 1.

As mentioned above, when the medical image 101 and the medical image 102 are completely identical, P = Q, that is, for all i (i = 1, 2, ... n), p _i = q _i , so it becomes S=1. When the medical image 101 is completely different from the medical image 102, since QP=1, that is, for all i(i=1, 2, . . . , n), q _i -p _i =1, so S =0. The S system is calculated only by the number of past medical images 101.

This device is the fifth, and the similarity S as shown above is associated with the past medical image in a one-to-one correspondence, and the corresponding information 104 is completed. In the sixth aspect of the present invention, the disease name of the past medical image corresponding to the maximum degree of similarity (close to "1") is specified as the disease name of the patient who is currently diagnosed.

In the example of Fig. 1, a bad situation as shown below occurs.

. Originally, it is not possible to generally expect a high-quality medical image in which a large variety of feature quantities can be obtained.

. When medical images that are not sufficiently high in quality are compared with each other, a general calculation can be performed, but a meaningless result is output. For example, the result of the difference cannot be detected, and the similarity of the same level is calculated for all past medical images.

. The processing takes a long time.

. Without the idea of changing according to the time axis, the level of the feature quantity itself is compared. Therefore, there is also a case where the characteristic amount of the medical portrait of the patient with a healthy subject is accidentally similar to the characteristic amount of the medical portrait of the rickety patient, and is misdiagnosed as a situation in which certain diseases are afflicted.

An outline of information processing in the present embodiment will be described with reference to Figs. 2 to 5 . The details are described in the following flow chart.

First, Fig. 2 will be explained. A past medical image 111 is acquired for a certain part of the patient to be diagnosed (the brain is also in the example of FIG. 2). On the other hand, the current medical image 112 is also obtained for the part of the patient. In the attack rate evaluation device 1 (Fig. 6) of the present embodiment, the past diagnosis example 114 can be used. The diagnostic example 114 can be assumed to be various. Here, for example, for each disease, a plurality of graphs showing the relationship between the rate of change and the number of diagnoses are prepared. For example, in the past 100 diagnostic cases, it is assumed that there are "the result of a 5% reduction in the area of the brain, five cases of being diagnosed as dementia" and "the result of a 8% reduction in the area of the brain, which is diagnosed as dementia. Case There are 10 pieces of information. Then, by plotting the points "-5%, 5 pieces", "-8%, 10 pieces", ..., a chart about "dementia" is created. Similarly, a chart of "mild cognitive impairment", "brain tumor", ... is created.

The attack rate evaluation device 1 of the present embodiment is the first, and determines a certain feature amount (category). It is assumed that the determined feature amount is, for example, "area". The attack rate evaluation device 1 is the second, and the difference in the feature amount is extracted. The difference in the "area" is only the difference between the sectional area of the portion in which the medical image 111 is expressed as the brain and the sectional area of the portion in which the medical image 112 is expressed as the brain. The difference can also be expressed as a numerical value (unit: cm ² ) by converting the number of pixels into an area. Further, the portion (contracted portion) having a specific pixel value which is present in the medical image 111 and not in the medical image 112 may be expressed as the medical image 113.

The attack rate evaluation device 1 is the third, and the rate of change R (percentage) of the difference is calculated by the following formula 2.

R=(q ₁ -p ₁ )/p ₁ ×100 (Expression 2)

Here, the area of the brain of the p ₁ -based medical image 111. The area of the brain of the medical image 112 of the q ₁ system.

The attack rate evaluation device 1 is the fourth, and the difference is applied to all the charts, and the number of diagnoses is obtained for each disease. The attack rate evaluation device 1 is the fifth, and the disease name is displayed in association with the attack rate in the order of large and small. The attack rate is, for example, a value (percentage) obtained by dividing the number of diagnoses obtained by the number of full diagnoses for the site (through all diseases).

Next, Fig. 3 will be described. Compared with the example of Fig. 2, in Fig. 3, only the feature quantity is "lightness" is different. The difference in "lightness" is only the brightness of the portion of the medical image 111 that is expressed as the brain (for example, the average value of the gray scale scale of 0 to 255 of each pixel), and the portion of the medical image 112 that is expressed as the brain. The difference in brightness. The differential system can be represented as a unitless value. Further, a color or the like determined in advance according to the difference size may be attached, such as a medical image 113.

In the disease system of the brain, there are various types of people whose incidence is "area", those whose incidence is "lightness", and those whose incidence is other characteristic quantities. The user can also retrieve the disease name by using the rate of change of a certain feature amount as an index key. Further, the combination of the rate of change of the complex feature amount (for example, "area" is reduced and "lightness" is increased) can be used as an index key to search for the disease name.

Next, Fig. 4 will be described. Figure 4 is similar to Figure 2 ("area" is a feature quantity). In Fig. 2, the rate of change of the area is a single value. On the other hand, in FIG. 4, the rate of change of the area indicates a complex value which changes in accordance with the change of the time series, which is greatly different. A past medical image 121 of a certain part (brain) of the patient to be diagnosed is acquired. On the other hand, a medical image 122 of the current and recent past of the part of the patient is also obtained.

The affliction rate evaluation device 1 is the first, and the medical image 121 of the past (std) which is the longest reference (photo) of the medical image 122 (now now, the previous now-1, the previous and the previous now-2) is extracted. The difference (symbol 123) calculates the rate of change 124 of the difference from the reference. Differential rate of change The number of medical images 122 is calculated. The attack rate evaluation device 1 is the second, and a time series chart 125 of the change rate is created. The attack rate evaluation device 1 is the third, and the time series chart 125 of the change rate is applied to the plural "type" (symbol 126). "Type" refers to a pattern of a time-series chart of the rate of change made based on past diagnostic examples (detailed later).

The attack rate evaluation device 1 is the fourth, and specifies the "type" of the time series chart 125 that is most approximated. The attack rate evaluation device 1 is the fifth, and displays a specific type of disease name in association with the percentage of the number of diagnoses of the disease divided by the number of full diagnostics of the type, that is, the attack rate. Among them, the attack rate evaluation device 1 is configured to associate one or more diseases and the number of pieces of diagnostic information 127 in association with the type according to the past diagnosis example.

In addition, FIG. 5 will be explained. Compared with the example of Fig. 4, in Fig. 5, only the feature quantity is "lightness" is different.

In the disease of the brain, there are various types such as the change of the time series of the change in the "area", the change of the time series of the change of the "lightness", and the change of the time series of the change of the other characteristic quantity. Everybody. The user can also search for the disease name by using the transition of the time series of the change in one feature amount as an index key. In addition, the combination of the time series transition of the change in the complex feature amount (for example, the "area" is repeatedly reduced and enlarged, and the "lightness" is continuously increased) is used as an index key to search for the disease name.

(machine composition)

The configuration of the attack rate evaluation device 1 will be described with reference to Fig. 6 . The attack rate evaluation device 1 is a general computer. The attack rate evaluation device 1 includes a central control device 11, an input device 12 such as a keyboard, a mouse, and a touch screen, an output device 13 such as a display, a main memory device 14, an auxiliary memory device 15, and a communication device 16. These are interconnected by a system bus.

The relative position correction unit 21, the difference image calculation unit 22, the change rate calculation unit 23, the abnormality presence determination unit 24, the attack rate calculation unit 25, and the input/output control unit 26 in the main memory device 14 are programs. In the following, when the main body is referred to as "○○ Department", the central control unit 11 reads each program from the auxiliary storage device 15 and loads it into the main memory device 14, and realizes the function of each program (details will be described later). The auxiliary memory device 15 stores image management information 31, abnormality determination information 32, physician diagnostic information 33, and medical image 34. Details of these are detailed later.

The input/output control unit 26 corresponds to an "input unit" and an "output unit".

The attack rate evaluation device 1 is connected to the terminal device 2 via the network 4 in a manner that allows communication with the terminal device 2. The terminal device 2 is also a general computer, and has a central control device connected to each other by a bus bar, an input device such as a keyboard, a mouse, a touch screen, an output device such as a display, a main memory device, an auxiliary memory device, and a communication device (not shown). ).

The imaging device 3 is a device that images a medical image by the body. There are many types of medical images captured by the imaging device 3, and for example, CT (Computed Tomography) images, MRI images, PET (Positron Emission Tomography) paintings can be used. Image, X-ray image, ultrasonic image, endoscope image, etc. Next, the imaging device 3 can output the captured image as a computer-readable digital image (consisting of pixels and pixel values) including the attack rate evaluation device 1 and the terminal device 2.

In general, the attack rate evaluation device 1, the terminal device 2, and the imaging device 3 are often placed in a hospital or the like. The doctor or the like uses the imaging device 3 to acquire a medical image of the body of the patient (including the healthy person) based on an opportunity such as a health diagnosis. Next, the doctor or the like inputs the medical image to the terminal device 2, and displays the medical image on the output device (display), for example. Next, the medical image is transmitted to the attack rate evaluation device 1 via the network 4. In this way, the attack rate evaluation device 1 responds to the terminal device 2 by using the received medical image and creating a "diagnosis result".

In FIG. 6, it is assumed that the plurality of terminal devices 2 are accessed by one attack rate evaluation device 1. However, it is also possible to form a configuration in which the attack rate evaluation device 1 and the terminal device 2 are housed in one housing. In addition, the attack rate evaluation device 1 may be formed into, for example, a plurality of frames for each body part (brain, lung, ...), each medical department (brain neurology, circulatory department, ...). Body composition.

The user who uses the terminal device 2 or the like may be other than the doctor, for example, a staff member of an insurance company, a patient himself who includes a healthy person, or the like.

Reflecting the health upsurge, collecting medical images captured by themselves or doctors to manage their own health status has recently been a lot. Next, the diagnosis of a large number of physicians is not intended to be published. However, it can also be assumed that the future In the case of the re-evaluation of the decree or the signing of a contract, etc., the case of the diagnosis is published on the condition of maintaining anonymity. In this way, the person who manages his/her own state of health can also obtain the diagnostic example and operate the attack rate evaluation device 1 of the present embodiment. Further, the health-related company may obtain a diagnosis example by a doctor or the like, receive a medical image of the customer from a general customer, and operate the attack rate evaluation device 1 to provide a diagnosis service to the customer.

(Portrait Management Information)

The image management information 31 will be described with reference to Fig. 7 . In the portrait management information 31, in association with the patient ID stored in the patient ID column 201, the patient name is stored in the patient name column 202, the portrait ID is stored in the portrait ID field 203, and the image is stored in the imaging device ID column 204. In the imaging device ID, the imaging point is stored in the point bar 205 at the time of recording, the part name is stored in the part name column 206, and the portrait file name is stored in the image file name field 207.

The patient ID of the patient ID column 201 is the identification code of the unique patient. In the present embodiment, "patient" means a general subject who obtains an image of the body. That is to say, "patients" also include the concept of "healthy people" who are not ill.

The patient name in the patient name field 202 is the patient's name.

The image ID of the portrait ID column 203 is an identification code of the unique medical image.

The imaging device ID of the imaging device ID column 204 is an identification code of the unique imaging device 3.

The imaging time point of the imaging time point column 205 is the date of the time when the medical image is captured.

The part name of the part name column 206 is the name of a part (part) of the patient's body. The name of the part may also be a bone such as a brain, a sputum, a nerve center, a heart, a lung, or the like, a pelvis, a thigh bone, and the like.

The image file name of the portrait file name column 207 is the name of the medical image as the digital portrait information.

(abnormal determination information)

The abnormality determination information 32 will be described with reference to Fig. 8 . In the abnormality determination information 32, in association with the part name stored in the part name column 211, the feature amount is stored in the feature amount column 212, the priority is stored in the priority column 213, and the memory is stored in the normal change rate range column 214. There is a range of normal change rate, and a range of normal change amounts is stored in the normal change amount range column 215, a unit is stored in the unit column 216, and a change period is stored in the change period column 217.

The part name of the part name column 211 is the same as the part name of FIG.

The feature quantity of the feature amount column 212 is a type of physical quantity of the part that can be acquired by the medical image. The characteristic quantity is the shape, the trait, and the like of the main performance part, and the specific examples of these are described later. In the following, the attack rate evaluation device 1 acquires the difference in the feature amount at two different points, and estimates the attack rate of each disease name (candidate) from the past diagnosis example based on the difference (details will be described later).

The priority of the priority column 213 is if there are a plurality of parts and one part name The corresponding feature quantity is a priority order between the feature quantities. The smaller the number, the higher the priority.

The normal rate of change range of the normal change rate range column 214 is a combination of the upper limit and the lower limit of the rate of change diagnosed by the physician or the like. The upper and lower limits do not have to be equally spaced around "0". For example, the upper limit is "+3%" and the lower limit is "-2%". The rate of change is the ratio (percentage) of the change in the characteristic quantity. The rate of change is calculated by the following formula 3.

Rate of change = (feature amount before change - feature amount after change) / feature amount before change × 100 (Expression 3)

The normal variation amount range of the normal variation amount range column 215 is a combination of the upper limit and the lower limit of the amount of change that the physician or the like diagnoses as normal. The amount of change is the absolute amount of change in the feature quantity. The amount of change is calculated by the following formula 4.

Change amount = feature quantity before change - characteristic quantity after change (Expression 4)

The unit of the unit column 216 is a unit of feature quantity. There are also cases where there is no unit. In the present embodiment, as in the case of "lightness", for example, a feature quantity expressed by a 2-bit binary number (in the case of a decimal number, 0 to 255) is formed as "no unit".

The change period of the change period 217 indicates that the rate of change of the feature amount and the amount of change are those that are converted into the rate of change and the amount of change over the period of time. For example, the rate of change of "+1%" during the change period is equal to a change of "+12%" in a 12-month period and "+6%" in a 6-month period. "+0.5%",... Similarly, the "+1" change in the "one month" period of the change period is equal to, for example, that a certain feature quantity changes "+12" in 12 months, "+6" in 6 months, and changes in half a month. "+0.5",...

However, in FIG. 8, "..." omits the value represented by each column, and does not mean that there is no numerical value in the field.

(Specific example of feature amount)

I assume a medical image of "brain". This medical image is assumed to be, for example, an MRI image shown by reference numeral 111 in Fig. 2 .

"Area": The attack rate evaluation device 1 can count, for example, the number of pixels having a pixel value of a specific range among the pixels of the medical image located in the region enclosed by the outline, and the number The area of the brain (correctly the sectional area) is calculated by converting it into an area or the like according to a predetermined rule.

"Brightness": The attack rate evaluation device 1 can calculate, for example, an average value of a pixel value (for example, a gray scale scale of 0 to 255) of a pixel located in a region enclosed by a contour in a medical image. In order to calculate the brightness of the brain. Further, the region may be subdivided, and the dispersion of the average value of the pixel values in the respective subdivided regions may be calculated to obtain the "lightness dispersion".

"Curve degree": The attack rate evaluation device 1 calculates the aspect ratio of the region enclosed by the contour in the medical image, and applies the contour shape or the like to a template or the like to calculate the brain. The circularity (correctly the circularity of the section).

In other examples, an endoscope image (color image) of the "stomach" is assumed.

"Color": The attack rate evaluation device 1 can calculate the stomach value of the pixel value of the pixel in the imaged region (for example, 0 to 255 for each of R, G, and B), thereby calculating the stomach. The color of the inner wall (hue).

"Texture": The attack rate evaluation device 1 can perform spatial frequency analysis, Fourier transform, etc., for pixels located in the imaged region, thereby calculating the texture statistics (concavity, Texture, etc.).

In other examples, an X-ray image of a hypothetical "lung" is assumed.

"Number of nodules": The attack rate evaluation device 1 can count the number of shadows having a predetermined brightness and circularity between straight contours (ribs), thereby calculating the nodules existing in the lungs (mostly The number of lesions such as tumors.

Although the other specific examples are not listed, the attack rate evaluation device 1 can acquire other characteristic amounts such as "blood flow velocity", "temperature", and "chroma" of a specific portion using an existing technique.

As shown in FIG. 8, for one part, there is a complex feature quantity. Of course, the specific feature quantity is individual case that is easy to obtain for a certain part, and is difficult to obtain for other parts (or cannot be assumed originally). However, in Fig. 8, in order to simplify the description, the individual cases as shown above are discarded.

(Physician Diagnostic Information)

The physician diagnostic information 33 is illustrated in accordance with FIG. Physician diagnostic information 33, in association with the patient ID stored in the patient ID column 221, the physician ID is stored in the physician ID column 222, the diagnostic image ID is stored in the diagnostic image ID column 223, and the site name is stored in the part name column 224. The rate of change is stored in the rate of change column 225, the amount of change is stored in the amount of change column 226, the name of the disease is stored in the disease name column 227, and the point of diagnosis is stored in the point bar 228 of the diagnosis.

The patient ID of the patient ID column 221 is the same as the patient ID of FIG.

The physician ID of the physician ID column 222 is the identification code of the physician who uniquely diagnoses the patient.

The diagnostic image ID of the diagnostic image ID column 223 is the same as the image ID of FIG. However, the image ID of the medical image (the medical image having the feature amount "p" in FIGS. 2 to 5) which is the reference is compared with the feature amount, and the doctor who compares the feature amounts is targeted. The combination of the image IDs of the portraits (the medical images having the feature amount "q" in FIGS. 2 to 5).

The part name of the part name column 224 is the same as the part name of FIG. However, the name of the part here may be a part name based on the patient's self-informed symptom, and is not limited to the presence of a disease in the part by the doctor's eye.

The rate of change of the rate of change column 225 is a vector in which the component has a rate of change of each feature amount. In Fig. 8, the feature quantity for the part name "brain" is "area", "lightness", "temperature", "color", "circularity", and "texture" in order of priority from high to low. "." The composition of the vector of the rate of change of FIG. 9 ranks the rate of change of the feature quantities in this order. By. That is, for example, in the first line of recording of Fig. 9, it is formed (the rate of change of area, the rate of change of brightness, ...) = (+0%, +0%, ...). The vector of the rate of change does not necessarily have a component for the rate of change of all feature quantities. For example, when the brain is only formed with "area" and "lightness" as a problem, other feature quantities are not memorized. At this time, the rate of change becomes (change rate of area, rate of change of brightness, rate of change of temperature, rate of change of color, ...) = (+0%, +0%, *, *, ...) . "*" indicates that the value does not exist.

The amount of change in the amount of change column 226 is a vector in which the component has an absolute amount of change in each feature amount. Similarly to the rate of change, for example, in the first line of recording of FIG. 9, the amount of change in area (the amount of change in brightness, ...) is (=0, +0, ...). The "*" is also the same as the rate of change.

The disease name of the disease name column 227 is the name of the disease to be diagnosed based on the medical image, the rate of change, and the amount of change specified by the diagnostic image ID. When viewed by the doctor's eye, if there is no disease at the site, the memory is "(no match)".

The diagnosis time point of the diagnosis time point column 228 is the date of the time point diagnosed by the physician.

Note that the records in the third to sixth rows of Fig. 9 and Fig. 7 show the following.

. The patient with the patient ID "P002" "Koizumi Erlang" is aware that the brain is abnormal.

. "Koizumi Erlang" received treatment from the doctor "D002" four times. The four calendar days are "20120415", "20120515" and "20120615". And "20120715".

. The doctor "D002" compares the medical images "I101" and "I102" at the initial diagnosis point "20120415". The medical image "I101" is obtained from the time point "20120315" before the first diagnosis. Koizumi Erlang brought the medical image "I101" to the doctor "D002" at the time of the initial diagnosis. Then, the medical image "I102" is the doctor "D002" who uses the imaging device "M004" to obtain the first diagnosis.

. The doctor "D002" uses the imaging device "M004" every time to obtain the medical image, and the acquired images "1103", "I104", and "I105" and the medical image " I101" for comparison.

. Regarding the "area" of the brain, the rate of change is "-1%" and the amount of change is "-5". These values are converted into the rate of change (amount) of the "change period of one month". As a result, if the rate of change is "-1%" for one month and "-5cm ² " for each change, the area of the brain continues to shrink.

. Regarding the "lightness" of the brain, the rate of change is as "+4%, -4%, +6%, -2%", and the amount of change is as "+8, -8, +12%, -4%". . As a result, the brightness system tends to increase while gradually increasing and decreasing the amplitude.

. The physician "D002" has consistently made a diagnosis of "mild cognitive impairment" after the initial diagnosis.

(Types of)

The type of time series chart is illustrated in accordance with FIG. The records in the third to sixth rows in Fig. 9 are for the record of the patient "P002". These records correspond to time series charts 231 and 232 of FIG. Note that the first (leftmost) component of each vector of the change rate column 225 of the record in the third row to the sixth row of FIG. 9 is memorized. This component indicates the rate of change of the area of the brain. The above records are sequentially, that is, the time series are "-1%", "-1%", "-1%", and "-1%". The transition of this time series is graphically represented by the time series chart 231 of FIG. Similarly, the second component (by left) of the respective vectors of the change rate column 225 of the record in the third row to the sixth row of Fig. 9 is memorized. This component indicates the rate of change of the brightness of the brain. The above records are sequentially, that is, the time series are "+4%", "-4%", "+6%", and "-2%". The time chart of this time series is the time series chart 232 of FIG.

The records in the 7th to 10th lines in Fig. 9 are for the record of the patient "P003". These records correspond to time series charts 233 and 234 of FIG. The first (most left) component of each vector of the change rate column 225 of the record in the seventh to the tenth rows of Fig. 9 is memorized. This component indicates the rate of change in the number of nodules in the lungs. The above records are sequentially, that is, the time series are "+3%", "+3%", "+4%", and "+4%". The transition of this time series is graphed as time series chart 233 of FIG. Similarly, the second component (by left) of the respective vectors of the change rate column 225 of the record in the 7th to 10th rows of Fig. 9 is memorized. This component is the rate of change in the lightness of the lungs, from the above record In order, the time series is "+2%", "+2%", "+3%", and "+3%". The transition of this time series is graphed as time series chart 234 of FIG.

The above is an example of charting the rate of change. However, similarly, the amount of change can also be graphed. Next, it is considered that the time series charts are divided into groups (patterning) according to their shapes, positions, and the like. An example of the group distinction is displayed as "type list" 35.

(type list)

In the type list 35 of FIG. 10, in association with the type stored in the type column 241, the sign transition is stored in the symbol change column 242, and the increase/decrease tendency column 243 stores the tendency of increase and decrease, and is shown in the time series chart. Example column 244 stores an example of a time series chart.

The type of the type column 241 is an identification code of the unique meaning type. The type here refers to the combination of "the transition of the symbol" and the "increasing or decreasing tendency".

The change of the sign of the symbol shift column 242 indicates how the sign of the change rate changes.

The tendency of increase or decrease in the increase/decrease tendency column 243 indicates how the level tendency of the change rate changes.

An example of a time series chart of the time series chart example column 244 is a symbol of a time series chart that matches the type. It can be seen that the time series chart 231 belongs to the type "a", the time series chart 232 belongs to the type "b", and the time series charts 233 and 234 all belong to the type "c".

The combination of the above-mentioned symbols and the tendency to increase or decrease are only As an example. In addition to these, the concave and convex (symbol of the second derivative function), the number of inflection points, the time series chart and the horizontal axis (rate of change = 0) of the time series chart are crossed, and the distance from the horizontal axis is imaginable. (amplitude), integral value (area between time series chart and horizontal axis), frequency, wavelength, presence or absence of convergence value (asymptote), level of convergence value, time series of non-continuous complex curve The combination of the transition, etc. as a type. That is, in general, the type means that if the time series chart is not repeatedly classified according to the result of mathematically and physics processing the time series chart, it may be any.

(processing order)

The processing procedure of this embodiment will be described below. There are two processing orders, and these are the basic first processing order and the second processing order belonging to the development type. The attack rate evaluation device 1 can calculate the attack rate by executing either one.

(first processing order)

The first processing sequence will be described with reference to Fig. 11 . In the case where the first processing sequence is started, the image management information 31, the abnormality determination information 32, and the physician diagnostic information 33 are stored in the auxiliary storage device 15 in the completed state. Further, all the medical images 34 in which the portrait ID is stored in the portrait management information 31 are stored in the auxiliary storage device 15. Next, for simplification of description, the attack rate evaluation device 1 and the terminal device 2 are formed as a component that is formed in one frame. That is, doctors, etc. It is assumed to be the input device 12 or the like that operates the attack rate evaluation device 1.

In step S301, the input/output control unit 26 acquires the medical image to be compared. Specifically, the input/output control unit 26 receives a patient ID, a part name, and a medical image through the input device 12 by a user such as a doctor. The medical image input here is the current medical image (symbol 112 of FIG. 2), and may be referred to as "comparison medical image" hereinafter. In addition, the patient ID and the part name which are received may be referred to as "target patient ID" and "target part name", respectively.

In step S302, the input/output control unit 26 acquires the medical image of the comparison standard. Specifically, the input/output control unit 26 is the first, and the target patient ID and the target part name are used as search keys, and the image management information 31 (FIG. 7) is searched for, and the record at the time of imaging is the latest record in the coincidence record.

The input/output control unit 26 is the second, and the auxiliary storage device 15 acquires the medical image having the image ID of the record acquired in the "first" of step S302. It is assumed that the portrait ID is attached to all medical images stored in the auxiliary memory device 15. The medical image obtained here may be referred to as a "comparative medical image" hereinafter.

In step S303, the relative position correcting unit 21 performs alignment of the medical image. Here, it is temporarily returned to FIG. 2, and the reference medical image 111 and the comparison medical image 112 are compared. The patient and the part to which the two images are to be applied are the same, but the imaging devices that obtain the images are different. In addition, if it is assumed that the two images are taken by the same imaging device, there is also a setting bar of the imaging device in the respective imaging time points. Parts, camera environment, etc. are different situations. As a result, the two portraits cannot be compared in the same state. Therefore, in order to superimpose the two images to obtain the difference (symbol 113), it is necessary to perform the alignment of the two images.

Specifically, the relative position correcting unit 21 is the first, and a point in which two bones are changed with little change in time is selected from among the comparative standard medical images. How to choose the 2 points shown above depends on the location. In the case of a medical image of the brain, the two points shown above are, for example, the center of gravity 111a of the skull (Fig. 2) and the apex 111b of the skull (Fig. 2). The following is an example of the brain to continue the explanation.

The relative position correction unit 21 is the second point, and the points corresponding to the two points selected in the "first" in step S303 are selected from the comparison target medical images (points 112a and 112b in Fig. 2).

The relative position correction unit 21 is the third, and the comparison target medical image and the comparison standard medical image are placed on the same plane. Next, the point 112a overlaps with the point 111a, and the point 112b overlaps with the point 111b to reduce or enlarge the entire comparison medical image. The medical image of the comparison target that has been reduced or enlarged as described above may be referred to as "the medical image for comparison after the position correction".

In step S304, the difference image calculation unit 22 determines the feature amount to be compared. Specifically, the difference image calculation unit 22 first searches the abnormality determination information 32 ( FIG. 8 ) using the target part name as a search key, and displays all the matching records on the output device 13 .

The difference image calculation unit 22 is the second, and the user receives the transmission input from among the feature quantities of the record displayed in the "first" of step S304. Set 12 to select 1 or a complex feature quantity (category). At this time, the difference image calculation unit 22 can also accept the number of categories in which the user inputs the feature amount. For example, when the "2" is input, the difference image calculation unit 22 is regarded as the feature quantity ("area") of the record having the target part name and having the priority of "2" or less in the abnormality determination information 32. "Brightness"). Here, the following description will be continued as the "area only" selected. The feature amount selected here is sometimes referred to as "selection feature amount" hereinafter.

The difference image calculation unit 22 is the third, and the user is allowed to select any one of the "change rate" and the "change amount" through the input device 12. Here, as the person who has selected the "change rate", the following description will be continued.

In step S305, the difference image calculation unit 22 creates a difference image. Specifically, the difference image calculation unit 22 is the first, and the comparison target medical image and the position correction comparison target medical image are superimposed, and the pixel value of the comparison medical image is determined for any pixel in the collection region. The pixel value of the comparison medical image is subtracted. Here, the "and collection area" refers to the in-contour area of the brain in the comparison medical image, and the sum of the in-contour areas of the brain in the medical image after comparison with the position correction (any one is consistent) section). The difference image calculation unit 22 can also perform the processing on the pixels located at all positions in the collection area. In addition, in the sum and region, from top to bottom, from top to bottom, such as the first, second, third, ... specific pixel position, for example, only for the "10" The integer multiples of the pixels in the position are processed (extraction processing).

The difference image calculation unit 22 is the second, and the pixel of the "first" subtraction result of the step S305 is exceeded by a certain threshold. Usually in the brain There is a significant difference between the pixel value of the pixel being imaged and the pixel value of the pixel (the space inside the brain).

The difference image calculation unit 22 is third, and sets the pixel value of the pixel specified in the "second" in step S305 to a predetermined pixel value (for example, a pixel value indicating "white"), and sets the pixel value of the other pixel. A difference image of other predetermined pixel values (for example, pixel values indicating "black"). With the difference image, the user can easily visualize which part of the part is reduced or enlarged.

In step S306, the change rate calculation unit 23 calculates a change rate and the like. Specifically, the change rate calculation unit 23 first counts the number of pixels belonging to the region within the outline of the brain among the pixels of the comparison standard medical image, and multiplies the number by a predetermined conversion. The result after the coefficient is set to the area "p".

The change rate calculation unit 23 counts the number of pixels belonging to the region within the contour of the brain among the pixels of the comparison medical image, and multiplies the number by the conversion factor. It is the area "q".

The change rate calculation unit 23 is the third, and calculates the change rate ((q-p)/p×100). When the "change amount" is received in "3rd" of step S304, the change rate calculation unit 23 calculates the amount of change (q-p).

In the fourth, the change rate calculation unit 23 converts the change rate calculated in the "third" of the step S306 into a change using the difference between the imaging time point of the comparison medical image and the number of days of the imaging time of the comparison medical image. The rate of change during the period. For example, if the difference in days is "50 days", When the "change period" corresponding to the selected feature amount is one month, the rate of change calculated in "third" of step S306 is multiplied by "30/50".

In step S307, the abnormality presence/absence determination unit 24 determines whether or not there is an abnormality. Specifically, the abnormality presence/absence determination unit 24 is the first, and the target portion name and the selected feature amount are used as search keys, and the abnormality determination information 32 (FIG. 8) is searched for, and the normal change rate range of the coincidence record is obtained. However, when the "change amount" is received in "3rd" of step S304, the abnormality presence/absence determination unit 24 acquires the normal change amount range.

The abnormality determination unit 24 is the second, and it is determined whether or not the change rate converted in the "fourth" of the step S306 is within the range of the normal change rate obtained in the "first" of the step S307, and the determination is made within the range. The result is "normal". If it is not within the range, the judgment result "abnormal" is generated.

In step S308, the attack rate calculating unit 25 calculates the attack rate. Specifically, the attack rate calculation unit 25 is the first, and the physician diagnosis information 33 (FIG. 9) is obtained, and the record that satisfies the following conditions 1 and 2 is obtained.

(Condition 1) The part name (column 224) coincides with the target part name.

(Condition 2) The rate of change corresponding to the selected feature amount (column 225) coincides with the rate of change converted in "4th" of step S306. Or, if not completely consistent, it is included within a predetermined tolerance.

Now, it is assumed that the target portion is named "brain", and the feature amount is selected as "area", and the rate of change converted in "4th" in step S306 is "-1%". In this way, the records of the third row to the sixth row and the eleventh row of FIG. 9 are obtained.

The attack rate calculation unit 25 is the second, and acquires the disease names of all the records acquired in the "first" of step S308. Next, the number of records obtained is kept for each disease name obtained. Here, "(mild cognitive impairment, 4), (dementia, 1)" is maintained.

The attack rate calculation unit 25 is the third, and counts the number of records (the number of target part records) that satisfy only "Condition 1". Here, the number of recorded parts is "7".

The attack rate calculation unit 25 is the fourth, and keeps the value of the record of the disease name divided by the number of the target part records for each of the recorded disease names acquired in the "first" of the step S308. rate). Here, "(mild cognitive impairment, 57%), (dementia, 14%)" is maintained. "57%" is a calculation result of 4/7×100, and "14%" is a calculation result of 1/7×100.

In step S309, the input/output control unit 26 displays the diagnosis result. Specifically, the input/output control unit 26 displays the diagnosis result display screen 51a (FIG. 13) on the output device 13, and displays the following information on the following fields of the diagnosis result display screen 51a.

. Patient ID column 131: Subject patient ID

. Patient Name Column 132: Patient Name Specific to Subject Patient ID

. Part name column 133: object part name

. "Previous" column 134a: Comparative standard medical portrait

. The previous time point column 134b: the imaging time point of the comparison medical image

. "This time" column 135a: Comparison of medical images after position correction

. At this time point column 135b: comparing the imaging time points of the medical image of the object

. Difference section 136: Difference portrait

The determination result table 137 includes a feature amount column 137a, a change rate column 137b, a change amount column 137c, and a determination result column 137d. In the feature amount column 137a, all the feature amounts corresponding to the target part name are displayed. Next, the input/output control unit 26 displays the converted rate of change (or the amount of change) in the "fourth" of step S306 in the change rate column 137b (or the change amount column 137c) of the row having the selected feature amount. Further, the determination result generated in "2nd" of step S307 is displayed in the determination result column 137d. Where "*" indicates that no processing has been performed for this field.

The attack rate table 138 has a disease name column 138a, an attack rate column 138b, and a determination column 138c. The input/output control unit 26 displays a combination of the disease name and the attack rate ("mild cognitive impairment, 57%", (dementia, 14%)" held in "4th" of step S308. In the disease name column 138a and the attack rate column 138b.

In step S310, the input/output control unit 26 registers the diagnosis result. Specifically, the input/output control unit 26 is the first, and the user such as the doctor accepts the check box input check mark displayed on the determination column 138c of any line of the attack rate table 138, and presses the "determination result registration" button. 139. The user such as the doctor confirms the comparison medical image 134a, the position-corrected comparison target image 135a, the difference image 136, and the determination result table 137, and inputs the check mark based on his/her own knowledge. It is not necessary to enter a check mark in the line where the attack rate of the attack rate table 138 is the highest. In addition, there is a case where the disease name which can be judged to be appropriate among the disease names showing the attack rate is not displayed. At this time, the user is doubled by the disease name column 138a. In the area below the line, select the name of the disease that is considered appropriate. Next, the check box under the double line enters the check mark, and the "OK result registration" button 139 is pressed. If the user does not input the check mark in any of the check boxes, and the "determination result registration" button 139 is pressed, the input/output control unit 26 regards that the disease name "(no match)" is selected.

The input/output control unit 26 is the second to create a new record of the physician diagnosis information 33 (Fig. 9), and the following information is stored in the following fields of the new record.

. Patient ID column 221: Subject patient ID

. Physician ID column 222: User's own physician ID

. Diagnostic Image ID column 223: Image ID of the comparison medical image and the image ID of the re-numbered medical image of the comparison

. Part name column 224: object part name

. The position of the component corresponding to the selected feature amount of the change rate column 225 (or the change amount column 226): the converted rate of change (or the amount of change) in "4th" of step S306

. The position of the other components of the change rate column 225 (or the change amount column 226): "*"

. Disease name column 227: Name of the disease corresponding to the entered check mark or "(no match)"

. Diagnostic time point 228: date of the current point

After that, the first processing sequence is ended.

(second processing order)

The second processing sequence will be described with reference to Fig. 12 . However, when the second processing sequence is started, it is assumed that the premise is the same as the first processing order.

In step S401, the input/output control unit 26 acquires the medical image to be compared. The processing of step S401 is the same as step S301 of the first processing sequence.

In step S402, the input/output control unit 26 acquires the medical image of the comparison standard. Specifically, the input/output control unit 26 is the first, and the target patient ID and the target part name are used as search keys, and the image management information 31 (FIG. 7) is searched, and the coincidence record is in the order of new and old at the time of imaging. A predetermined number of records is obtained. Here, the predetermined number is set to "3", and the following description will be continued.

The input/output control unit 26 is the second, and the auxiliary storage device 15 acquires the medical image having the image ID of the record acquired in the "first" of step S402. It is assumed that all the medical images stored in the auxiliary memory device 15 are attached with the image ID. In the case of the medical image obtained here, the image is the longest, and may be referred to as a "comparative medical image" hereinafter, and the other (two) are referred to as "intermediate medical image".

In step S403, the relative position correcting unit 21 performs alignment of the medical image. The processing of step S403 is the same as step S303 of the first processing procedure. However, the relative position correction unit 21 is configured to perform alignment with the "comparative medical image" for the "intermediate medical image" and to create a "post-correction intermediate medical image".

In step S404, the difference image calculation unit 22 determines The amount of features that should be compared. The processing of step S404 is the same as step S304 of the first processing sequence.

In step S405, the difference image calculation unit 22 creates a difference image. The processing of step S405 is the same as step S305 of the first processing sequence. However, the difference image calculation unit 22 is configured to create all of the following difference imagers.

. The difference between the "basic medical image" and the "image of the intermediate medical image after the position correction is earlier"

. "Differential image of the middle of the intermediate medical image after the position correction" and "the later of the intermediate medical image after the position correction"

. "Differential image of the image of the medical image in the middle of the intermediate medical image after the position correction" and the medical image of the comparison after the position correction

In step S406, the change rate calculation unit 23 calculates a change rate and the like. The processing of step S406 is the same as step S306 of the first processing sequence. However, the change rate calculation unit 23 is configured to calculate all of the following change rates (or change amounts).

. The rate of change (or the amount of change) of the "comparative medical image" in the "middle medical image"

. The rate of change (or the amount of change) of the "comparative medical image" in the "middle medical image"

. The rate of change (or amount of change) of the "comparison target medical image" relative to the "comparative medical image"

Here, the change rate calculation unit 23 is for three different imaging points. In the medical image, the rate of change of the medical image that is the standard is obtained. Of course, it is also possible to obtain the rate of change of the medical image to be used as a reference for four or more medical images having different imaging points.

In step S407, the abnormality presence/absence determination unit 24 determines whether or not there is an abnormality. Specifically, the abnormality presence/absence determination unit 24 is the first, and the target portion name and the selected feature amount are used as search keys, and the abnormality determination information 32 (FIG. 8) is searched for, and the normal change rate range of the coincidence record is obtained. However, when the "change amount" is received in "3rd" of step S404 (S304), the abnormality presence/absence determination unit 24 acquires the normal change amount range.

The abnormality presence/absence determination unit 24 is the second, and determines whether or not all of the converted rate of change in the "fourth" of the step S406 (S306) is within the range of the normal change rate obtained in the "first" of the step S407. When the rate of change is within the range, the determination result is "normal", and in other cases, the determination result "abnormal" is generated.

In step S408, the attack rate calculation unit 25 classifies the past diagnosis examples. Specifically, the attack rate calculation unit 25 is the first, and the physician diagnosis information 33 (FIG. 9) is obtained, and all the records satisfying the following conditions 3 to 5 are acquired.

(Condition 3) The part name (column 224) coincides with the target part name.

(Condition 4) The rate of change corresponding to the selected feature amount (not the "*") is stored in the change rate column 225.

(Condition 5) The patient ID is not the subject patient ID.

The attack rate calculation unit 25 is the second, and the records acquired in the "first" of the step S408 are sorted for each patient ID, and created. A record group per patient ID.

The attack rate calculation unit 25 is the third, and deletes the number of records belonging to the record group among the record groups created in the "second" in step S408, which does not satisfy the predetermined number. The predetermined number at this time is, for example, "4". In other words, if there is a record of "4" or more, the time series of the change rate of the records can be handled as an intentional person.

The attack rate calculating unit 25 is the fourth time series, and creates a time-series chart of the selected feature amount for the patient based on the recorded group that has not been deleted in the "third" of step S408. In other words, the attack rate calculating unit 25 creates a time-series chart for the rate of change of the target portion "brain" and the selected feature amount "area" for each patient ID (see, for example, reference numeral 231 in FIG. 10). Here, the method of creating the time series chart is based on the method described in step S406.

The attack rate calculating unit 25 is the fifth type, creates a plural type, and classifies the time series chart created in the "fourth" of step S408 by the type of creation. Next, the patient ID, type, and disease name of the time series chart are temporarily held in association with each other. The type means a pattern in which the time series chart is processed mathematically and physics in an arbitrary manner, and the time series chart is not repeatedly classified (as described above with reference to FIG. 10). The information temporarily held at this time is, for example, "(P001, type d, mild cognitive impairment), (P002, type e, dementia), (P003, type f, brain tumor), (P004, type g, meningitis) ), (P006, type d, mild cognitive impairment), (P007, type f, brain tumor), (P009, type g, meningitis),...". The following system will have "ID, type, disease name)" is called "classification completion information". One classified completion information corresponds to one time series chart.

In step S409, the attack rate calculation unit 25 classifies the transition of the change in the image of the medical image to be compared. Specifically, the attack rate calculation unit 25 is the first, and based on the complex change rate calculated in step S406, the comparison target medical image is used to obtain the medical image to be compared via the intermediate medical image. Time series chart of the rate of change of patients with patient ID.

The attack rate calculation unit 25 is the second, and determines whether or not the time series chart created in the "first" of the step S409 matches any of the types created in the "5th" of the step S408. The type determined at this time is sometimes referred to as "patient type" below.

In step S410, the attack rate calculation unit 25 calculates the attack rate. Specifically, the attack rate calculating unit 25 is the first, and counts the number of patients having the patient type for each disease name among the classified information held in "5th" of step S408. Next, the counted amount is temporarily remembered in association with the disease name. The information temporarily held at this time is, for example, "(mild cognitive impairment, 8), (dementia, 4), (meningitis, 2), ...".

The attack rate calculation unit 25 is the second, and divides the number (8, 4, 2, ...) of each disease name counted in "1st" of step S410 by the total number of classified information having the patient type. The result, that is, the percentage, is temporarily maintained in association with the disease name. Now suppose that the total number of classified information with patient types is "20". As a result, in the foregoing In the example, the attack rate calculation unit 25 maintains ((mild cognitive impairment, 40%), (dementia, 20%), (meningitis, 10%), ...". "40%" is the calculation result of 8/20×100, "20%" is the calculation result of 4/20×100, and "10%" is the calculation result of 2/20×100.

In step S411, the input/output control unit 26 displays the diagnosis result. Specifically, the input/output control unit 26 displays the diagnosis result display screen 51b (FIG. 14) on the output device 13, and displays the following information on the following fields of the diagnosis result display screen 51b.

. Patient ID column 131: Subject patient ID

. Patient Name Column 132: Patient Name Specific to Subject Patient ID

. Part name column 133: object part name

. "Base" column 140a: comparison standard medical portrait

. "Previous and previous time" column 140b: The image is taken earlier in the middle medical image after the position correction

. "Previous" column 140c: After the position correction, the intermediate medical image is later than the camera

. "This time" column 140d: Comparison of medical images after correction of position

. Reference time point column 141a: comparison imaging time point of the reference medical image

. The front and the next time point column 141b: the imaging time point of the earlier imaging point in the middle medical image

. The previous time point column 141c: the imaging time point of the person who is late in the middle medical image

. This time point column 141d: comparing the imaging time points of the medical image of the object

. Column 142a: comparing the reference medical image with the position correction A difference image of an earlier image taken during a medical image

. Column 142b: a difference image of the middle of the intermediate medical image after the position correction, and the later of the middle medical image after the position correction

. Column 142c: a difference image of the medical image that is compared with the position correction after the position correction in the intermediate medical image after the position correction

. Column 143: Time series chart created in "1st" of step S409

Regarding the determination result table 137, the description about FIG. 13 is directly applied. However, the rate of change (or amount of change) is a combination of the lower limit and the upper limit (indicating the upper and lower amplitudes of the time series chart). Regarding the attack rate table 138, the description about FIG. 13 is also directly applied.

In step S412, the input/output control unit 26 registers the diagnosis result. Specifically, the input/output control unit 26 performs the same processing as the processing of "first" and "second" in step S310 of the first processing procedure.

The input/output control unit 26 is the second, and the record having the target patient ID is acquired by the physician diagnosis information 33 (FIG. 9). Next, the acquired disease name of the record is updated to the disease name or "(no match)" corresponding to the input check mark, and is overwritten. If there is no record with the target patient ID in the physician diagnosis information 33, nothing is done.

After that, the second processing sequence is ended.

(Modification 1)

In the "second" of the step S304 of the first processing procedure, the processing when the plural "selected feature amount" (category) is selected is as follows.

(1) The difference image calculation unit 22 or the like repeats the processing of steps S305 to S308 for each selection feature amount.

(2) The input/output control unit 26 performs the same processing as that of the above-described step S309, except for the processing of step S309.

. In the "Difference" column 136 (Fig. 13), a difference portrait regarding the complex selection feature amount is displayed. It can also be displayed in a slideshow format such as a difference image (5 seconds) of "area" → a difference image (5 seconds) of "lightness".

. The rate of change (or the amount of change) and the determination result are displayed on all the lines of the selected feature amount having the determination result table 137.

. In the attack rate column 138b, the attack rate, such as "(area, brightness, temperature, color, ...) = (20%, 15%, *, 10%, ...)", and the selected feature amount Produce a related display.

(3) The input/output control unit 26 performs the same processing as that of the above-described step S310 with respect to the processing of step S310.

(Modification 2)

In the "second" of the step S404 (S304) of the second processing procedure, the processing when the plural "selected feature amount" (category) is selected is as follows.

(1) The difference image calculation unit 22 or the like repeats the processing of steps S405 to S410 for each selection feature amount.

(2) The input/output control unit 26 performs the same processing as that of the above-described step S411 except for the processing of step S411.

. In the "difference" columns 142a, 142b, and 142c (Fig. 14), a difference image for the complex selection feature amount is displayed. It can also be displayed in a slideshow format.

. In column 143, a time series chart showing the number of feature selections is displayed. It can also be displayed in a slideshow format.

. The combination of the upper limit and the lower limit of the change rate (or change amount) and the determination result are displayed on all the lines of the selected feature amount having the determination result table 137.

. In the attack rate column 138b, the attack rate, such as "(area, brightness, temperature, color, ...) = (20%, 15%, *, 10%, ...)", and the selected feature amount Produce a related display.

(3) The input/output control unit 26 performs the same processing as that of the above-described step S412 with respect to the processing of step S412.

(Effects of the embodiment)

The attack rate evaluation device 1 of the present embodiment achieves the following effects.

(1) The attack rate evaluation device 1 searches for a past diagnosis example using the difference of the received feature amounts. At this time, if the number of categories of the feature amount is limited, the medical image of the past diagnosis example or the medical image received is not required to be of high image quality. Therefore, it is possible to effectively use a large amount of data that does not necessarily require high image quality such as regular health diagnosis. Further, the attack rate evaluation device 1 calculates the difference. Therefore, if the medical image captured by the same camera 3 is continuously used, the measurement error (deformation of the image, etc.) inherent to the camera can be eliminated.

(2) The attack rate evaluation device 1 uses either or both of the "rate" and "quantity" of the difference. Therefore, the multi-faceted attack rate can be evaluated according to the category of the feature amount.

(3) The attack rate evaluation device 1 feeds back the diagnosis result as the physician diagnosis information 33. Therefore, the more the number of uses is increased, the smaller the influence of the random measurement error occurring in the imaging of the medical image is, the higher the diagnostic accuracy is.

(4) The attack rate evaluation device 1 determines "abnormal" or "normal" according to the type of the feature amount. Therefore, it is easy to know which type of feature quantity has undergone a large change.

(5) The attack rate evaluation device 1 searches for a diagnosis example using transitions at different time points of three or more different times. Therefore, the diagnosis results are formed to reflect the longer-term perspective, and have a good impact on the patient's psychology.

(6) The attack rate evaluation device 1 uses the feature amount of the plural category. Therefore, it is possible to flexibly and flexibly use the feature amount that conforms to the characteristics of the site or disease.

(7) The attack rate evaluation device 1 displays a difference image. Therefore, the variation of the difference can be easily recognized.

(8) The attack rate evaluation device 1 converts the difference into a difference with respect to a period of a predetermined length. Therefore, the difference system is formed to be easy to understand and easy to compare according to the type of the part, the disease, and the feature amount.

(9) The attack rate evaluation device 1 performs alignment of medical images. Therefore, it is possible to compare medical images having different setting conditions or imaging environments.

However, the present invention is not limited to the foregoing embodiments, and includes various modifications. For example, the above-described embodiments are described in detail for easy understanding of the present invention, and are not necessarily limited to those having all of the constituents described. Further, a part of the configuration of a certain embodiment may be replaced with a configuration of another embodiment, and a configuration of another embodiment may be added to the configuration of a certain embodiment. Further, addition, deletion, and replacement of other configurations may be performed for a part of the configuration of each embodiment.

Further, each of the above-described configurations, functions, processing units, processing means, and the like may be implemented by a part of or all of the components, for example, in an integrated circuit design or the like. Further, each of the above-described configurations, functions, and the like may be implemented in software by explaining and executing a program in which the processor realizes respective functions. The programs, tables, files, and other information for realizing each function can be placed in a recording device such as a memory or a hard disk, an SSD (Solid State Drive), or a recording medium such as an IC card, an SD card, or a DVD.

In addition, the control line or the information line is indicated as indicating that it is convenient to be considered as required, and all control lines or information lines are not necessarily displayed on the product. In fact, it can be considered that almost all the components are connected to each other.

51a‧‧‧Diagnosis results display screen

131‧‧‧ Patient ID column

132‧‧‧ Patient Name Column

133‧‧‧ part name column

134a‧‧‧"Previous" column

134b‧‧‧Last time bar

135a‧‧‧ "This time" column

135b‧‧‧Time point bar

136‧‧‧"Differential" column

137‧‧‧Results table

137a‧‧‧Characteristics column

137b‧‧‧change rate column

137c‧‧‧Changes column

137d‧‧‧Results column

138‧‧‧ nuisance rate table

138a‧‧‧ disease name column

138b‧‧‧ Suffering rate column

138c‧‧‧determination column

139‧‧‧"Determining Results Login" button

Claims (11)

An attack rate evaluation device is provided with a memory unit and a control unit that stores a difference between two images at different imaging points when the image is captured in association with a part of the body. And the diagnosis information of the disease name, the control unit includes: an input unit that receives a portion to be diagnosed by the user, and a difference image calculation unit that receives the type of the feature amount of the user input image, and In the above-mentioned part of the diagnosis target, the difference in the type of the feature quantity of the image having the different imaging point is extracted, and the change rate calculation unit calculates the diagnosis target for the same diagnosis target for the type of the feature quantity received as described above. In the above-mentioned part, the difference between the two images at the time of imaging is performed, and the attack rate calculation unit searches for the diagnostic information and acquires the coincidence record by using the portion to be diagnosed and the calculated difference as a search key. The name of the disease, based on the number of the aforementioned coincident records relative to the portion having the aforementioned Recording the percentage of the number of calculated attack rate; and an output unit, which were based diseases acquired by each of the foregoing, the display of the calculated prevalence. The attack rate evaluation device according to the first aspect of the patent application, wherein the difference is caused by two different imaging points when imaging the portion The amount of change and/or absolute value of the value obtained by the image. The prevalence rate evaluation device according to the second aspect of the invention, wherein the control unit adds a record having the received portion, the calculated difference, and the acquired disease name as a new record of the diagnostic information. . The prevalence rate evaluation device according to the third aspect of the invention, wherein the memory unit stores the determination information that is diagnosed as a numerical range in which the difference is normal in association with the part and the type, and the control unit The portion to be diagnosed, the accepted type, and the calculated difference are used as search keys, and the determination information is searched to obtain a numerical range of the coincidence record, and whether the calculated value range includes the calculated value The difference is used to determine whether or not the portion to be diagnosed is normal. The attack rate evaluation device according to the fourth aspect of the invention, wherein the difference is a time-series display of a change in the difference between the three or more images having different imaging points. The prevalence rate evaluation device according to claim 5, wherein the accepted category is a plural type, and the control unit calculates the attack rate for the plural category. In the attack rate evaluation device according to the sixth aspect of the invention, the control unit displays an image in which only the difference between the two images at different imaging points when the imaging target is captured is displayed. The attack rate evaluation device according to claim 7, wherein the control unit converts the calculated difference into a predetermined length The difference during the period. The attack rate evaluation device according to the eighth aspect of the invention, wherein the control unit performs two images having different imaging points when imaging the portion to be diagnosed by superimposing two sets of corresponding points on each other. The opposite. An attack rate evaluation method according to the present invention, wherein the control unit is configured to receive a type of feature quantity of an image input by a user through an input device, and to In the part of the body to be diagnosed, the difference in the type of the feature amount of the image at the time of imaging is extracted for the part to be diagnosed, and the type of the feature amount received is calculated as the same diagnosis target. The difference between the two images at the time of imaging at the imaging target is the difference between the two images that are the target of the diagnosis, and the calculated difference is used as a search key, and the search is performed in association with the part. When the imaging is performed, the difference between the two portraits and the diagnosis information of the disease name are obtained, and the matching disease name is obtained, and the attack rate is calculated based on the ratio of the number of the matching records to the number of records having the aforementioned portion. Each of the aforementioned disease names shows the aforementioned calculated disease . An attack rate assessment program that is used to assess the prevalence rate The attack rate evaluation program that functions as a function is characterized in that the memory unit of the attack rate evaluation device performs the following processing: storing two different imaging points at which the image is captured in association with the body part. In the diagnosis of the disease rate, the control unit of the attack rate evaluation device performs a process of accepting the type of the feature amount of the user input image and the part to be diagnosed, and the diagnosis is performed. In the above-mentioned part of the object, the difference in the type of the feature amount of the image at the time of imaging is extracted, and the imaging time at which the image of the same diagnosis target is imaged is calculated for the type of the received feature amount. In the difference between the two images, the portion to be diagnosed and the calculated difference are used as search keys, and the diagnostic information is searched for, and the disease name of the coincidence record is obtained, and the number of the matching records is relative to the portion having the aforementioned portion. The ratio of the number of records, calculate the attack rate, per Disease name of the aforementioned achievements, the aforementioned display the calculated prevalence.