WO2015146520A1 - 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム - Google Patents

類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム Download PDF

Info

Publication number
WO2015146520A1
WO2015146520A1 PCT/JP2015/056370 JP2015056370W WO2015146520A1 WO 2015146520 A1 WO2015146520 A1 WO 2015146520A1 JP 2015056370 W JP2015056370 W JP 2015056370W WO 2015146520 A1 WO2015146520 A1 WO 2015146520A1
Authority
WO
WIPO (PCT)
Prior art keywords
case
similarity
similar
interest
individual
Prior art date
Application number
PCT/JP2015/056370
Other languages
English (en)
French (fr)
Inventor
大沢 哲
Original Assignee
富士フイルム株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士フイルム株式会社 filed Critical 富士フイルム株式会社
Priority to CN201580016673.9A priority Critical patent/CN106164903B/zh
Publication of WO2015146520A1 publication Critical patent/WO2015146520A1/ja
Priority to US15/275,851 priority patent/US10839511B2/en

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/56Details of data transmission or power supply, e.g. use of slip rings
    • A61B6/563Details of data transmission or power supply, e.g. use of slip rings involving image data transmission via a network
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/51Indexing; Data structures therefor; Storage structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/58Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
    • G06F16/583Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using metadata automatically derived from the content
    • G06F16/5838Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using metadata automatically derived from the content using colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/502Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion

Definitions

  • the present invention relates to a similar case search device, a similar case search method, and a similar case search program.
  • the examination image is an image taken with a modality such as a CT (Computed Tomography) apparatus that performs tomography or a general X-ray imaging apparatus that captures a simple fluoroscopic image, and diagnoses the patient such as identifying a patient's disease. Used for.
  • CT Computer Tomography
  • a general X-ray imaging apparatus only one inspection image may be captured or a plurality of inspection images may be captured.
  • a plurality of tomographic images are acquired. Therefore, one inspection data includes one or more inspection images. Since cases are often created by collecting past examination data, one case data also includes one or more case images.
  • a region of interest is designated in the examination image by a user such as a doctor.
  • the region of interest is a region in which the doctor is particularly interested in the examination image and includes a lesion to be diagnosed.
  • the similar case retrieval apparatus compares a feature quantity obtained by quantifying the feature of one region of interest designated in the examination image with a feature quantity obtained by quantifying the feature of one lesion existing in the case image. The degree of similarity is determined.
  • a lesion included in the region of interest of the examination image is referred to as a target lesion
  • a lesion included in the case image is referred to as a case lesion.
  • the similar case search device searches for a case including a case lesion similar to the region of interest from a case database storing a plurality of cases.
  • Japanese Patent Application Laid-Open No. 2010-237930 discloses that there is an individual difference for each user in the method of specifying a region of interest including a target lesion, and that the search variation in which the search result varies due to the individual difference is suppressed.
  • the target invention is disclosed. Specifically, even when an area including the same target lesion is designated as a region of interest, the shape and size of the designated region may differ depending on how the region of interest is designated due to individual differences among the designated users. Since it changes, the feature amount may change. When the feature amount changes, the similarity also changes, and thus a search variation in which a search result changes for each user occurs.
  • 2010-237930 calculates feature amounts for each region of interest for a plurality of regions of interest that are specified differently for one target lesion, for example, in order to suppress such search variations. Similarity image search is performed by calculating the similarity based on the average value of the feature values of a plurality of regions of interest. According to this, the search dispersion
  • Japanese Patent Application Laid-Open No. 2012-118583 (US Publication US2012 / 134555) relates to a technique for outputting a search result more suitable for the subjectivity that a user feels similar to. Specifically, when target lesions of the same type exist for a plurality of examination images, a region of interest including a plurality of target lesions of the same type that the user feels similar when designating the region of interest As one group. Then, in one examination data, a range of feature amounts including all the feature amounts of a plurality of target lesions belonging to the same type group is obtained, and similar case search is performed using the range of feature amounts as a search condition. Since the range of the feature amount of the same type group is considered to coincide with the range that the user feels to be subjectively similar, the search result can be made more suitable for the user's subjectivity.
  • the appearance of a plurality of target lesions in the examination image may be the basis for specifying the disease.
  • the disease may be identified by the presence of three types of target lesions in the examination image: cavity shadows, punctate shadows, and ground glass shadows.
  • a disease may be identified by the appearance of two types of target lesions, an abnormal bronchial shadow and a punctate shadow.
  • a similar search with a single target lesion was sufficient.
  • a similar case search with these multiple target lesions is necessary.
  • JP 2010-237930 A and JP 2012-118583 A both apply to one target lesion included in the examination image. Focusing on searching for similar cases based on the feature amount of the region of interest including one focused target lesion, and considering focusing on each of the plurality of targeted lesions included in the examination image is considered There wasn't.
  • Japanese Patent Application Laid-Open No. 2010-237930 calculates the feature amount for each region of interest, but the plurality of regions of interest are different in the way of designation, and the region of interest of the same target lesion is the same. There is no disclosure about searching for similar cases by paying attention to the feature amount of each region of interest of different target lesions.
  • Japanese Patent Laying-Open No. 2012-118583 discloses 1 for one homogeneous group in which the same type of target lesion is grouped with respect to a plurality of target lesions included in a plurality of examination images. One search condition is created, and similar cases are searched using the created search condition.
  • 2012-118583 (US Publication US2012 / 134555) only calculates a feature amount common to each region of interest including a plurality of target lesions of the same type so as to meet the user's preference. Thus, there is no disclosure about searching for similar cases by paying attention to the feature amount of each region of interest of a plurality of target lesions.
  • An object of the present invention is to provide a similar case search device, a similar case search method, and a similar case search program capable of searching for an appropriate similar case search in a short time even when there are a plurality of regions of interest. .
  • the similar case retrieval apparatus of the present invention is an apparatus for retrieving a similar case similar to an examination image used for diagnosis of a patient from a case database in which a plurality of cases including one or more case images are registered.
  • a quantity acquisition unit, an individual similarity calculation unit, an overall similarity calculation unit, and a similar case search unit are provided.
  • the feature amount acquisition unit is a region of interest that is designated in examination data including one or more examination images and is designated to include a target lesion that is a lesion existing in the examination image, and each of one or more different target lesions. With respect to a plurality of regions of interest that are included, feature quantities for each region of interest are acquired.
  • the individual similarity calculation unit associates the region of interest with a plurality of case lesions existing in the case image on a one-to-one basis, and compares the feature amount of the corresponding region of interest with the feature amount of the case lesion. Thus, the individual similarity for each region of interest is calculated.
  • the total similarity calculation unit is a total similarity calculation unit that calculates a total similarity based on a plurality of individual similarities calculated for a plurality of regions of interest, and a plurality of cases of different types in the same case Only for completely different combinations configured by combining a lesion and a plurality of regions of interest, the total similarity is calculated using the individual similarity.
  • the similar case search unit searches for similar cases based on the total similarity corresponding to completely different combinations.
  • case lesions present in a case image when there are a plurality of case lesions within one case image, for example, when there is one case lesion for each of two case images.
  • the total similarity calculation unit can create a number of combinations according to the number of regions of interest and the number of types of case lesions and the number of complete heterogeneous combinations, and calculate the total similarity for each completely different combination.
  • the total similarity calculation unit obtains individual similarities for case lesions included in each case for each case, and creates completely different combinations for the obtained individual similarities. You may ask for the degree. Alternatively, first, individual similarities for the case lesions of all cases may be obtained, and then a completely different combination of individual similarities may be created for each case. In any case, after calculating the individual similarity ISM, completely different combinations of the individual similarity ISM are created.
  • the overall similarity calculation unit distinguishes the plurality of case lesions of the same type, It is preferable to create a completely different combination for each case lesion of the same type that has been distinguished.
  • the similar case search unit creates a similar case list in which information on a plurality of similar cases is listed based on the total similarity.
  • the similar cases are preferably sorted in the order of the overall similarity.
  • the display items of the similar case list include the value of the total similarity and the breakdown information regarding the total similarity, and the breakdown information may include the correspondence between the region of interest and the case lesion for calculating the individual similarity. preferable.
  • the similar case list in addition to the value of the overall similarity, it is preferable to display a plurality of individual similarity values that are calculation elements of the overall similarity.
  • the similar case list preferably displays images of the region of interest and the case lesion.
  • the case database preferably stores information on the type of case lesion.
  • the total similarity is preferably the sum of a plurality of individual similarities included in a completely different combination.
  • the individual similarity calculation unit may create an individual similarity table in which a plurality of individual similarities calculated in association with a plurality of case lesions are recorded for each region of interest.
  • the total similarity calculation unit reads out the individual similarity one by one from each of the plurality of individual similarity tables for each region of interest, and creates a completely different combination using the read out individual similarity as an element. Good.
  • the individual similarity calculation unit excludes cases for which there are at least a plurality of case lesion types from among individual cases, and does not calculate the individual similarity for cases with one case lesion type. It is preferable.
  • the total similarity calculation unit may perform a weighting process on the total similarity according to the value of the individual similarity that is an element for calculating the total similarity.
  • the weighting process is preferably a process for increasing the overall similarity when the individual similarity is equal to or greater than a threshold value.
  • the similar case search unit may exclude cases where the number of types of case lesions is less than the number of regions of interest. In addition, the similar case search unit may search for cases where the number of types of case lesions is less than the number of regions of interest when there are a plurality of types of case lesions.
  • one representative value is selected from the plurality of total similarities. It is preferable to have a representative value determination unit for determining, and for the similar case search unit to search for a similar case based on the representative value.
  • the similar case retrieval method of the present invention is a method for retrieving a similar case similar to an examination image used for diagnosis of a patient from a case database in which a plurality of cases including one or more case images are registered.
  • the feature amount acquisition step is a region of interest that is designated in examination data including one or more examination images and is designated to include a target lesion that is a lesion existing in the examination image, and each of one or more different target lesions. With respect to a plurality of regions of interest that are included, feature quantities for each region of interest are acquired.
  • the individual similarity calculation step associates the region of interest with a plurality of case lesions existing in the case image on a one-to-one basis, and compares the feature amount of the corresponding region of interest with the feature amount of the case lesion. Thus, the individual similarity for each region of interest is calculated.
  • the total similarity calculation step is a total similarity calculation unit that calculates a total similarity based on a plurality of individual similarities calculated respectively for a plurality of regions of interest, and a plurality of cases of different types in the same case Only for completely different combinations configured by combining a lesion and a plurality of regions of interest, the total similarity is calculated using the individual similarity.
  • similar case search step similar cases are searched based on the total similarity corresponding to completely different combinations.
  • a completely different combination of case lesions may be created for each case, and the individual similarity corresponding to the created combination may be obtained. Moreover, you may make the completely different combination of individual similarity for every case. In any case, after calculating the individual similarity ISM, completely different combinations of the individual similarity ISM are created.
  • the similar case search program of the present invention causes a computer to execute a process of searching for a similar case similar to a test image used for diagnosis of a patient from a case database in which a plurality of cases including one or more case images are registered.
  • the program includes a feature amount acquisition step, an individual similarity calculation step, an overall similarity calculation step, and a similar case search step.
  • the feature amount acquisition step is a region of interest that is designated in examination data including one or more examination images and is designated to include a target lesion that is a lesion existing in the examination image, and each of one or more different target lesions. With respect to a plurality of regions of interest that are included, feature quantities for each region of interest are acquired.
  • the individual similarity calculation step associates the region of interest with a plurality of case lesions existing in the case image on a one-to-one basis, and compares the feature amount of the corresponding region of interest with the feature amount of the case lesion. Thus, the individual similarity for each region of interest is calculated.
  • the total similarity calculation step is a total similarity calculation unit that calculates a total similarity based on a plurality of individual similarities calculated respectively for a plurality of regions of interest, and a plurality of cases of different types in the same case Only for completely different combinations configured by combining a lesion and a plurality of regions of interest, the total similarity is calculated using the individual similarity.
  • similar case search step similar cases are searched based on the total similarity corresponding to completely different combinations.
  • the present invention compares the feature amount for each region of interest with the feature amount of the case lesion, calculates the individual similarity for each region of interest, and calculates a plurality of case lesions and a plurality of regions of interest within the same case. Even when there are multiple regions of interest, the total similarity is calculated using the individual similarity only for completely different combinations configured in combination, and similar cases are searched based on the calculated total similarity.
  • a similar case search device a similar case search method, and a similar case search program capable of searching for an appropriate similar case search in a short time.
  • FIG. 30 is an explanatory diagram of combinations that are not subject to calculation in the example of FIG. 29. It is explanatory drawing which shows the completely different combination of the example of FIG.
  • FIG. 32 is an explanatory diagram of combinations that are not subject to calculation in the example of FIG. 31. It is explanatory drawing of the case which becomes a search object. It is explanatory drawing of another case which is not a search object. It is a table
  • a medical information system 9 shown in FIG. 1 is constructed in a medical facility such as a hospital.
  • the medical information system 9 includes a medical department terminal 11 arranged in the medical department 10, a modality (medical image capturing apparatus) 13 installed in the examination department 12, an order management terminal 14, and an examination image database (hereinafter referred to as “DB”). ) Server 15, case DB server 16, and similar case search server 17. These are communicably connected through the network 18.
  • the network 18 is, for example, a LAN (Local Area Network) laid in the hospital.
  • the modality 13 is, for example, a CT (Computed Tomography) apparatus or MRI (Magnetic Resonance Imaging) apparatus that captures tomographic images, and a general X-ray imaging apparatus (DR: Digital Radiography, CR: Computed Radiography, etc.) that captures simple fluoroscopic images. ) Etc.
  • CT Computer Tomography
  • MRI Magnetic Resonance Imaging
  • DR Digital Radiography
  • CR Computed Radiography, etc.
  • the clinical department terminal 11 is operated by a doctor of the clinical department 10 (indicated by Dr in the figure), and in addition to inputting and viewing the electronic medical record, issuing an examination order for requesting the examination to the examination department 12. Used.
  • the medical department terminal 11 also displays an examination image 19 photographed by the examination department 12 and stored in the examination image DB server 15, and is also used as an image display terminal for a doctor to view the examination image 19.
  • the order management terminal 14 accepts the examination order from the clinical department 10, and manages the accepted examination order.
  • the technician of the examination department 12 takes a picture of the patient with the modality 13 according to the contents of the examination order.
  • One or a plurality of inspection images 19 are taken for one inspection order.
  • the modality 13 transmits the photographed inspection image 19 to the inspection image DB server 15.
  • the examination department 12 notifies the doctor of the medical department 10 that the examination is completed, and the storage destination of the examination image 19 in the examination image DB server 15 is also notified.
  • the doctor of the medical department 10 accesses the inspection image DB server 15 through the medical department terminal 11 and browses the inspection image 19 on the medical department terminal 11.
  • the inspection image DB server 15 includes an inspection image DB 20 in which a plurality of inspection images 19 are stored, and is a so-called PACS (Picture Archiving and Communication System) server.
  • the inspection image DB 20 is a database that can be searched by keywords, and distributes the inspection image 19 that matches the search condition or the designated inspection image 19 in response to a search request or distribution request from the medical department terminal 11 or the like.
  • one inspection data 21 including one or more inspection images 19 is stored in association with one inspection order.
  • the inspection image 19 photographed by the CT apparatus or the MRI apparatus is a tomographic image (also called a slice image), and one inspection data 21 includes a plurality of inspection images 19.
  • the inspection image 19 taken by a general X-ray imaging apparatus it is a simple fluoroscopic image, and one inspection data 21 may include only one inspection image 19. In some cases, there are multiple sheets.
  • the examination order includes requester information such as ID (Identification Data) and affiliation of the doctor of the medical department 10, patient information, examination type, and the like.
  • the image file of the inspection image 19 includes image data and incidental information such as DICOM (Digital Imaging and COmmunication in Medicine) headers. Information on the inspection order is stored as incidental information of the inspection image 19. Further, the incidental information includes an inspection ID and an image ID assigned to each inspection image 19. In the example of FIGS. 2 and 3, the inspection ID is “O901” or “O902”, the image IDs are “O901-03”, “O901-01”, and one inspection image 19 is assigned to the inspection ID. It is given in the form of adding a serial number for identification.
  • the inspection image DB server 15 can search using items included in such DICOM tags as search keys.
  • the similar case retrieval server 17 accepts the examination image 19 as a retrieval condition and retrieves a case including the case image 22 similar to the accepted examination image 19.
  • the case image 22 is an examination image used for diagnosis in the past.
  • the case DB server 16 has a case DB 23 that stores a plurality of cases in a searchable manner.
  • the similar case search server 17 accesses the case DB server 16, reads out one case at a time, compares the inspection image 19 received as a search condition with the case image 22 in the case, and displays the inspection image 19. Search for similar cases.
  • the doctor of the medical department 10 operates the medical department terminal 11 to access the inspection image DB server 15 and download the inspection data 21 including the inspection image 19.
  • the examination image 19 is displayed on the medical department terminal 11 and viewed by a doctor.
  • the examination image 19 of the patient includes a lesion (referred to as a target lesion OL) that shows a symptom of the disease.
  • the doctor of the medical department 10 selects the examination image 19 including the target lesion OL from the examination images 19 included in the examination data 21.
  • the examination image 19 is attached to a similar case search request issued to the similar case search server 17 by the medical department terminal 11 and transmitted to the similar case search server 17.
  • the similar case search server 17 searches the case DB server 16 for cases similar to the examination image 19 and distributes the search result to the requesting medical department terminal 11.
  • the doctor of the medical department 10 confirms the cases included in the test results.
  • the case includes an interpretation report accompanying the case image 22.
  • the doctor makes a definitive diagnosis such as specifying a disease in the examination image 19 with reference to the findings on the case image 22 described in the interpretation report.
  • the case DB 23 is provided with a case image DB 23A and a feature amount DB 23B.
  • the case image DB 23A is a database that stores the case images 22 so as to be searchable.
  • a case ID is assigned to each case.
  • the case ID corresponds to the examination ID in the examination image 19.
  • One case includes one or more case images 22.
  • each case image 22 is given an image ID obtained by adding a serial number to the case ID.
  • the case data 24 with the case ID “C101” includes, for example, 60 tomographic images.
  • the case image 22 includes a lesion representing a disease symptom (case lesion CL).
  • case lesion CL a lesion representing a disease symptom
  • one or more case lesions CL are registered.
  • three case lesions CL are registered in the case ID “C101”, three from No1 to No3, two in “C102”, and one in “C103”.
  • the case lesion CL is an area designated as a lesion by a doctor when the case image 22 has been used as an examination image in the past, and is registered as a case lesion CL after a definitive diagnosis by the doctor. .
  • the method for specifying the case lesion CL is, for example, the same as the region of interest ROI described later.
  • the feature amount DB 23B is a database that stores the feature amount CAC of the image of the case lesion CL.
  • the feature amount CAC is given an identification ID consisting of a case ID and a lesion number.
  • the case ID “C101” includes three case lesions CL, and each case lesion CL is given No. 1 to No. 3 which are serial numbers within one case.
  • the number after the feature amount CAC corresponds to the serial number in the case.
  • lesion patterns which are typical lesion image patterns, are classified into, for example, eight types A to H. That is, A: abnormal shadow in a low respiratory region (emphysema, pneumothorax, bra, etc.), B: cavity, C: abnormal shadow of bronchi (bronchial wall thickening, bronchodilation, traction bronchodilation, bronchial transillumination, etc.), D : Honeycomb, E: Frosted glass-like shadow, F: Point-like shadow (granular shadow, TIB, etc.), G: Abnormal shadow in high absorption area (consolidation, nodule, bronchial mucous gland, etc.), H: Linear / reticular shadow is there.
  • A abnormal shadow in a low respiratory region (emphysema, pneumothorax, bra, etc.)
  • B cavity
  • C abnormal shadow of bronchi (bronchial wall thickening, bronchodilation, traction bronchodilation, bronchi
  • the type of lesion of the case lesion CL is determined by which of the eight lesion patterns the image pattern of the case lesion CL is closest to.
  • the type of case lesion CL is determined based on the feature amount CAC.
  • the feature amount CAC is an eight-dimensional feature vector composed of output values of eight types of discriminators corresponding to eight types of lesion patterns, similarly to a region of interest ROI described later.
  • the type information of the case lesion CL is stored in the feature amount DB 23B together with the feature amount CAC.
  • a typical lesion pattern is classified into 8 types, and the corresponding lesion types are classified into 8 types. However, the number may be less than 8 types, or more than 8 types. Good.
  • the region including the target lesion OL in the examination image 19 is designated as the region of interest ROI by the doctor.
  • the similar case search request is attached with the examination image 19 including the target lesion OL and area information (for example, coordinate information in the examination image 19) in the examination image 19 corresponding to the designated region of interest ROI.
  • the similar case search server 17 specifies the region of interest ROI based on the image data of the examination image 19 and the region information. Then, the feature amount of the region of interest ROI is calculated.
  • the similar case search server 17 After calculating the feature amount, the similar case search server 17 reads cases one by one from the case DB server 16, compares the feature amounts of the region of interest ROI and the case lesion CL, and searches for similar cases. Then, a similar case list in which information on a plurality of similar cases is listed is created, and is distributed to the medical department terminal 11 as a search result.
  • the medical department terminal 11, the order management terminal 14, the examination image DB server 15, the case DB server 16, and the similar case search server 17 are based on a computer such as a personal computer, a server computer, and a workstation, and a control program such as an operating system.
  • An application program such as a client program or a server program is installed.
  • the computers constituting the DB servers 15 to 17 and the terminals 11 and 14 have the same basic configuration, and are a CPU (Central Processing Unit) 41, a memory 42, and a storage device 43, respectively. , A communication I / F 44, and an input / output unit 46. These are connected via a data bus 47.
  • the input / output unit 46 includes a display unit 48 and an input device 49 such as a keyboard and a mouse.
  • the storage device 43 is, for example, an HDD (Hard Disk Drive), and stores a control program and an application program (hereinafter referred to as AP) 50.
  • AP application program
  • a disk array in which a plurality of HDDs are connected is provided in the server in which the DB is constructed, as the DB storage device 43.
  • the disk array may be built in the server main body, or may be provided separately from the server main body and connected to the server main body through a cable or a network.
  • the memory 42 is a work memory for the CPU 41 to execute processing.
  • the CPU 41 loads the control program stored in the storage device 43 into the memory 42 and executes processing according to the program, thereby comprehensively controlling each part of the computer.
  • the communication I / F 44 is a network interface that performs transmission control with the network 18.
  • client programs such as electronic medical record software for browsing and editing electronic medical records and viewer software for browsing examination images and similar case lists are installed as AP50.
  • the viewer software may be dedicated software or a general-purpose WEB browser, for example.
  • the display unit 48 ⁇ / b> A of the medical department terminal 11 has an operation function based on GUI (Graphical User Interface).
  • the inspection image display screen 52 is displayed.
  • the CPU 41 ⁇ / b> A of the medical department terminal 11 functions as the GUI control unit 53 and the search request issue unit 54.
  • On the examination image display screen 52 it is possible to specify a region of interest ROI for the examination image 19 and issue a similar case search request.
  • the GUI control unit 53 receives an operation instruction from the input device 49A through the inspection image display screen 52, and performs screen control according to the received operation instruction.
  • the input issuance instruction is input from the GUI control unit 53 to the search request issuing unit 54.
  • the search request issuing unit 54 issues a similar case search request by attaching the specified examination image 19 and region information of the region of interest ROI to the similar case search request.
  • the inspection image display screen 52 includes an image display area 52A for displaying the inspection image 19 and various operation units.
  • the image display area 52A for example, three inspection images 19 are displayed side by side.
  • the inspection image 19 to be displayed can be switched by a scroll operation or a frame advance operation.
  • An input box 52B for inputting an examination ID is provided above the image display area 52A.
  • the examination ID is entered in the input box 52B, the examination data 21 of the entered examination ID is downloaded from the examination image DB server 15 and the examination image 19 is displayed in the image display area 52A.
  • an area designation button 52C, a clear button 52D, and a similar case search button 52E are provided below the image display area 52A.
  • the region designation button 52C is an operation button for designating the region of interest ROI in the inspection image 19.
  • an area designation operation for designating an arbitrary area in the inspection image 19 becomes possible.
  • the pointer 56 is operated, and the outer periphery of the region including the target lesion OL is designated by, for example, a spline.
  • the spline is a smooth curve passing through a plurality of designated control points, and is input by designating the control points with the pointer 56.
  • the region including the target lesion OL is designated as the region of interest ROI.
  • the clear button 52D is a button for canceling the designated region of interest ROI.
  • a plurality of regions of interest ROI can be specified.
  • No. 1 to No. 3 ROIs of interest are designated for each of the three inspection images 19 having the image IDs “O901-01” to “O901-03”. If there is no other region of interest ROI designation in the examination data 21 with the examination ID “901”, a total of three regions of interest ROI are designated in one examination data 21.
  • two regions of interest ROI (No1, No2) are designated in the inspection image 19 with the image ID “O906-01”, and the image IDs are “O906-02” and “O906-03”.
  • Regions of interest ROI (No 3 and No 4) are designated one by one in the inspection image 19.
  • the region of interest ROI of No3 includes two target lesions OL (No3, No4). As described above, a region including a plurality of target lesions OL may be designated as one region of interest ROI. If there is no other region of interest ROI designation in the examination data 21 with the examination ID “906”, a total of four regions of interest ROI are designated in one examination data 21. Each region of interest ROI designated in this way includes one or more different target lesions OL.
  • a similar case search server program is installed as an AP 50 in the similar case search server 17.
  • the CPU 41 ⁇ / b> B of the similar case search server 17 receives a request reception unit 61, a feature amount calculation.
  • the request reception unit 61 receives a similar case search request transmitted from the medical department terminal 11 and stores the received examination image 19 and region information of the region of interest ROI in the storage device 43 of the similar case search server 17 or the like.
  • the feature amount calculation unit 62 calculates the feature amount of the region of interest ROI based on the received inspection image 19 and region information.
  • the feature amount calculation unit 62 functions as a feature amount acquisition unit.
  • the feature amount of the region of interest ROI of No1 is “RAC1”
  • the feature amount of the region of interest ROI of No2 is “RAC2”
  • the region of interest ROI of No3 is calculated for each region of interest ROI.
  • the feature amount RAC is a feature vector composed of multi-dimensional numerical values (discriminator output values described later) corresponding to a plurality of types of lesion patterns preset as typical lesion image patterns. Typical lesion patterns are classified into, for example, eight types A to H as described with reference to FIG.
  • the feature amount calculation unit 62 includes discriminators 62A to 62H corresponding to eight typical types of lesion patterns.
  • Each of the discriminators 62A to 62H outputs a numerical value corresponding to each typical lesion pattern based on the image pattern of the region of interest ROI.
  • Each numerical value output from each discriminator 62A to 62H is a multi-dimensional numerical value constituting the feature vector, and here, each numerical value is referred to as a discriminator output value.
  • the discriminator output value represents the likelihood of a typical lesion pattern, and is a value representing the degree to which a typical lesion pattern exists in the region of interest ROI. Therefore, the larger the discriminator output value, the higher the degree of existence of a typical lesion pattern in the region of interest ROI, and the smaller the discriminator output value, the lower the degree of existence. More specifically, a discriminator output value of “+” indicates that a typical lesion pattern exists in the region of interest ROI, and a discriminator output value of “ ⁇ ” does not exist in the region of interest ROI. Indicates. An output value “+” is indicated, and the greater the discriminator output value, the higher the degree of presence.
  • the lesion pattern discriminator 62B of “B: cavity” and the discriminator 62G of “G: high absorption region” show an output value of “+”, and the output value of “B: cavity” is the highest. Since the region of interest ROI is large, the region of interest ROI includes lesion patterns of “B: cavity” and “G: high absorption region”. Among the eight types of lesion patterns, “B: cavity” has a dominant image pattern. It turns out that it is.
  • Each classifier for typical lesion patterns is, for example, “Document name: Computer Vision and Image Image Understanding Volume 88, pages 119 to 151, published December 2002, Using Human Perceptual Categories for content-based retrieval from a medical. Created by machine learning algorithms such as “Ada-boost” using well-known features such as Image Database Author Chi-Ren Shyu, Christina Pavlopoulou Avinash C.kak, and Cala E.Brodley be able to.
  • the feature amount calculation unit 62 calculates a feature amount RAC for each region of interest ROI for all the plurality of regions of interest ROI specified in the examination data 21 attached to the similar case search request.
  • the feature quantity CAC of each case lesion CL stored in the feature quantity DB 23B in the case DB 23 is also composed of feature vectors corresponding to the above eight types of lesion patterns.
  • Such a feature quantity CAC is also calculated by a configuration similar to that of the feature quantity calculation unit 62 described in FIG. 14, and is configured by an eight-dimensional feature vector including eight types of output values of the discriminators 62A to 62H. .
  • each case lesion CL stores the associated lesion type.
  • the individual similarity calculation unit 65 compares the feature amount RAC of the region of interest ROI with the feature amount CAC of the case lesion CL, and calculates the individual similarity ISM. Specifically, the individual similarity ISM is calculated by comparing the feature quantity RAC and the eight-dimensional feature vector included in the feature quantity CAC. The value of the individual similarity ISM is calculated by, for example, the least square distance or the correlation. In the former case, the smaller the value, the higher the similarity between the region of interest ROI and the case lesion CL. In the latter case, the greater the value, the higher the similarity between the region of interest ROI and the case lesion CL.
  • the individual similarity calculation unit 65 has a one-to-one correspondence between a plurality of regions of interest ROI included in one examination data 21 and a case lesion CL included in one case data 24.
  • each feature amount RAC is compared with the feature amount CAC, and the individual similarity ISM is calculated. Since the individual similarity ISM is the similarity for each region of interest ROI, it is the individual similarity for each region of interest, and is calculated in a one-to-one correspondence with the case lesion CL, so the individual similarity for each case lesion CL But there is.
  • the individual similarity calculation unit 65 calculates the individual similarity ISM by the number of case lesions CL for one region of interest ROI. Since the individual similarity ISM is calculated for all regions of interest ROI, the individual similarity calculator 65 calculates the individual similarity ISM by the number obtained by multiplying the number of regions of interest ROI and the number of case lesions CL.
  • the case data 24 with the case ID “C101” has three cases No1 to No3.
  • a lesion CL is registered.
  • the types of the three case lesions No. 1 to No. 3 are “B: cavity”, “F: dot-like shadow”, and “E: ground glass-like shadow”, respectively. Therefore, three types of three case lesions CL are registered.
  • the individual similarity calculation unit 65 of this example calculates the individual similarity ISM noting the type of the case lesion CL but focusing only on the number of the case lesions CL. Therefore, as shown in FIG. 18, a total of nine individual similarity ISMs of 3 ⁇ 3 are calculated between the examination data 21 of “O901” and the case data 24 of “C101”.
  • the identification code in parentheses attached to each individual similarity ISM is obtained by adding the serial numbers of the region of interest ROI and the case lesion CL to the case ID.
  • this indicates an individual similarity ISM between the region of interest ROI of No1 and the case lesion CL of No1 registered in the case data 24 with the case ID “C101”.
  • C101-12 indicates that the region of interest ROI of No1 and the individual similarity ISM of the case lesion CL of No2 with the case ID “C101”.
  • the individual similarity calculation unit 65 associates the five case lesions CL of “C102” with respect to each of the three regions of interest ROI regardless of the type, a total of 15 individual similarities of 3 ⁇ 5 The degree ISM is calculated.
  • the case ID “C103” shown in FIG. 20 is an example in which two types of four case lesions CL are registered.
  • the four case lesions CL of No. 1 to No. 4 of “C103” have two types, “B: cavity” and “F: dot”, and “B: cavity” has the same case lesion CL ( No. 1 to No. 3), and “F: dotted shadow” has one case lesion CL (No. 4).
  • the individual similarity calculation unit 65 associates four case lesions CL of “C103” with respect to each of the three regions of interest ROI regardless of the type, and a total of 12 individual similarities in 3 ⁇ 4.
  • the degree ISM is calculated.
  • the 21 is an example in which four types of case lesions CL of one type are registered.
  • the four case lesions CL of No. 1 to No. 4 of “C104” are all of the same type with the type “B: Cavity”.
  • the individual similarity calculation unit 65 associates four case lesions CL of “C103” with respect to each of the three regions of interest ROI regardless of the type, and a total of 12 individual similarities in 3 ⁇ 4.
  • the degree ISM is calculated.
  • the case with the case ID “C106” shown in FIG. 22 is an example in which the same three types and five case lesions CL as “C101” shown in FIG. 19 are registered, but the breakdown is different. That is, there are three types of lesions CL of No. 1 to No. 5 of “C106”: “B: cavity”, “F: dotted shadow”, and “E: ground glass shadow”. Is the same as “C101”. However, in “C106”, two case lesions CL (No1, No2) of the same type are registered for “B: cavity”, and the same case lesion CL (No3) for “F: punctiform shadow” is registered. , No4) are registered. One (No. 5) is registered for “E: ground glass shadow”.
  • the individual similarity calculation unit 65 associates each of the three regions of interest ROI with five case lesions CL of “C106” regardless of the type, and a total of 15 individual similarities in 3 ⁇ 5. The degree ISM is calculated.
  • the individual similarity calculation unit 65 creates an individual similarity table (hereinafter referred to as an ISM table) 71 in, for example, the memory 42B or the storage device 43B of the similar case search server 17 to generate an ISM.
  • the calculated individual similarity ISM is recorded in the table 71.
  • the ISM table 71 is created for each region of interest ROI.
  • the example of FIG. 23 shows the ISM table 71 for the region of interest ROI of No1.
  • the ISM table 71 is a table that stores a case ID, a lesion No., a lesion type, and a lesion image in association with each individual similarity ISM.
  • the lesion image is image data of the case lesion CL. That is, in the ISM table 71, one record is composed of five items of data of case ID, lesion No., lesion type, lesion image, and individual similarity ISM.
  • the individual similarity calculation unit 65 first records the individual similarity ISM in the order in which the individual similarity ISM is calculated in the ISM table 71. Each individual similarity ISM is recorded in order from the smallest case ID number, for example, “C101”, “C102”, and “C103”. Since the value of the individual similarity ISM is calculated by the correlation between the feature amount RAC of the region of interest ROI and the feature amount CAC of the case lesion CL, the larger the numerical value, the higher the similarity.
  • the individual similarity calculation unit 65 creates an ISM table 71 for each region of interest ROI. When there are three regions of interest ROI No1 to No3, three ISM tables 71 are created. At this stage, as shown in FIG. 23, in the ISM table 71, the records are arranged in the order of the case IDs. The individual similarity calculation unit 65 passes the ISM table 71 to the total similarity calculation unit 66 when the creation of the ISM table 71 is completed.
  • the total similarity calculation unit 66 creates a total similarity TSM table 72 (hereinafter referred to as a TSM table 72) based on a plurality of ISM tables 71 created for each region of interest ROI. Specifically, as shown in FIG. 26, individual similarity ISMs calculated by correspondence with case lesions CL in the same case are read one by one from a plurality of ISM tables 71 and read from each ISM table 71. Based on the plurality of individual similarities ISM, the total similarity TSM is calculated. More specifically, the total similarity calculation unit 66 sets the number of regions of interest ROI and the number of case lesions CL as the individual similarity ISM read from each ISM table 71 as an element of the total similarity TSM. A corresponding number of individual similarity ISM combinations are created, and an overall similarity TSM is calculated for each combination.
  • the individual similarity ISM is a one-to-one correspondence between a plurality of regions of interest ROI and a plurality of case lesions CL, and depends on which one of each region of interest ROI and each case lesion CL is associated with which. Since the value of the individual similarity ISM changes, the total similarity TSM changes depending on the combination pattern of the individual similarity ISM.
  • Each individual similarity ISM is the similarity when the region of interest ROI and the case lesion CL are made to correspond one-to-one. Therefore, the higher the total similarity TSM, the three regions of interest ROI and the three case lesions CL. This means that the average value of each individual similarity ISM is high. In this example, since the individual similarity ISM is represented by a correlation value, the larger the numerical value, the higher the similarity. Therefore, the overall similarity TSM also indicates that the similarity is higher as the numerical value is larger.
  • the total similarity TSM of “C101-2” is the highest at “2.04”.
  • the total similarity TSM of “C101-3” is the lowest at “1.32”.
  • the total similarity TSM “C101-5” includes the individual similarity ISM (C101-13) indicating the highest value “0.91” in the individual similarity ISM.
  • the total similarity TSM of “C101-2” that does not include the value has a higher numerical value because the individual similarity ISM is higher on average.
  • the total similarity TSM (C101-2) having the highest similarity among the six total similarity TSMs is the individual similarity ISM (C101-11) between the No. 1 region of interest ROI and the No. 1 case lesion CL, and No. 2.
  • the total similarity calculation unit 66 creates a completely different type of combination of the individual similarities ISM for calculating the total similarity TSM.
  • the completely heterogeneous combination is a combination of a plurality of case lesions CL in the same case, each of which is extracted for each type, and combined. Furthermore, it means all combinations when corresponding to each region of interest ROI. Then, the total similarity calculation unit 66 determines a completely different combination as a calculation target of the total similarity TSM, and determines a combination other than a completely different combination as a calculation target. That is, the total similarity calculation unit 66 calculates the total similarity TSM using the individual similarity ISM only for completely different combinations.
  • Combinations other than completely different combinations include completely the same kind of combinations and some of the same kinds of combinations.
  • a plurality of case lesions CL of the same type are extracted one by one from a plurality of case lesions CL in the same case, and a combination of the extracted plurality of case lesions CL of the same type is further interested. This means all combinations when corresponding to each of the regions ROI.
  • Some of the same types of combinations are different in type from a plurality of case lesions CL of the same type from among cases including at least one case lesion CL of different types and two or more case lesions CL of the same type.
  • each individual similarity ISM is symbols representing the types of case lesions CL corresponding to each individual similarity ISM.
  • Each individual similarity level ISM represents which type of case lesion CL is calculated corresponding to the region of interest ROI.
  • the number of permutations such as “C101”, is obtained when the three regions of interest ROI are extracted from the five case lesions CL of “C102”.
  • the number of completely different combinations is the same kind of combination “B1, B2, B3”, “B1, B3, B2” shown in FIG. 30, or “B1, B2, F”, “B1”. , B2, F ”are the numbers excluding the same kind of combination.
  • the total similarity calculation unit 66 determines that the total similarity combination shown in FIG. 29 is not subject to calculation of the total similarity TSM, and calculates a combination other than the complete heterogeneity combination shown in FIG. Judged out of scope.
  • the search time for similar cases can be shortened.
  • the main reason for searching for similar cases by paying attention to a plurality of regions of interest ROI is the presence of three types of target lesions: a cavity, a dotted shadow, and a ground glass-like shadow. This is because it is used for diagnosis of non-cancer diseases such as tuberculosis that is the basis for disease identification. Therefore, it is necessary to search for a case including a different case lesion CL as a similar case as a case similar to a plurality of regions of interest ROI including a different target lesion OL.
  • a plurality of regions of interest ROI are often designated to include different types of target lesions OL, and if so, the same type of case lesion CL is included in the first place. Since the necessity of extracting cases is low, even if the calculation target of the total similarity TSM is limited to completely different combinations, there is little influence on actual diagnosis.
  • the total similarity calculation unit 66 further includes the case of “C103” including two types and four case lesions CL shown in FIG. 33 and one type and four case lesions shown in FIG.
  • the total similarity TSM is not calculated.
  • the reason for excluding cases such as “C103” where the number of types of case lesions CL is less than the number of regions of interest ROI is not necessary to be searched. That is, in the diagnosis of a non-cancer disease, the doctor designates three regions of interest ROI because it desires to search for cases including three types of case lesions CL, and the number of types less than the number of regions of interest ROI. This is because even if the case lesion CL is extracted, the necessity is considered to be low in relation to the diagnosis. Note that in this example, the case of “C103” including two types of case lesions CL is excluded from the calculation target because this example assumes three regions of interest ROI. In the case of two, of course, “C103” is also a target of calculation of the total similarity TSM.
  • the total similarity calculation unit 66 indicates the case lesion CL as in the cases with the case IDs “C101”, “C102”, “C105”, and “C106”. Only for cases where the number of types is equal to or greater than the number of regions of interest ROI, the total similarity TSM is calculated and used as a search target for similar cases. The total similarity TSM is calculated only for completely different combinations. Conversely, for cases where the number of registered case lesions CL is less than the number of regions of interest ROI, such as “C103” and “C104” (in this example, there are less than three case lesions CL), the total similarity TSM is calculated. Not calculated and excluded from search.
  • the total similarity calculation unit 66 records the total similarity TSM calculated for a plurality of cases to be searched in the TSM table 72.
  • the TSM table 72 one record is composed of three items of data of a combination pattern of case ID, total similarity TSM, and individual similarity ISM.
  • the combination pattern is a combination of individual similarities ISM that is a calculation source of each total similarity TSM.
  • the overall similarity calculation unit 66 delivers the created TSM table 72 to the similar case search unit 67.
  • the similar case retrieval unit 67 is provided with a representative value determination unit 67B.
  • the representative value determining unit 67B determines a representative value for each case from a plurality of total similarity TSMs in the TSM table 72.
  • the total similarity TSM (C101-2) having the highest similarity (correlation value is the maximum) among the six total similarity TSMs is determined as the representative value.
  • the similarity is the highest (correlation value is maximum), and the total similarity TSM (C102-3) is determined to be the representative value. Is done.
  • the representative value determination unit 67B performs the representative value determination for all cases in the TSM table 72. Thereby, only the representative total similarity TSM determined for each case is extracted from the TSM table 72.
  • the similar case search unit 67 sorts the records in descending order of the total similarity TSM in the TSM table 72 from which representatives have been extracted. Since only one total similarity TSM is extracted from one case in the TSM table 72, sorting each record means sorting the cases in descending order of the total similarity TSM. Thereby, ranking is performed for each case, and a similar case having a high similarity is extracted at the top of the TSM table 72.
  • the similar case search unit 67 is provided with a list creation unit 67A (see FIG. 12). Based on the TSM table 72, the list creating unit 67A creates a similar case list 74 shown in FIG. The similar case list 74 is displayed in the search result display screen 76.
  • the similar case list 74 is a list of information on a plurality of similar cases.
  • the search result display screen 76 is an example of a screen in which the similar case search server 17 distributes the search result to the medical department terminal 11 that is the request source of the similar case search request.
  • the list creation unit 67A extracts the total similarity TSM from the TSM table 72, and creates a similar case list 74 in which similar cases are arranged in descending order of the total similarity TSM.
  • the display items of the similar case list 74 include, for each total similarity TSM, the value of the total similarity TSM, the rank based on the total similarity TSM, the case ID, and the breakdown information regarding the total similarity TSM.
  • the breakdown information is a correspondence relationship between the region of interest ROI and the case lesion CL for calculating each individual similarity ISM that is a calculation element of the total similarity TSM.
  • the breakdown information includes the type of the case lesion CL (lesion type).
  • the total similarity TSM with the case ID “C106” corresponds to the total similarity TSM of “C106-1” shown in FIG.
  • the combination pattern of the individual similarity ISM that is a calculation element of the total similarity TSM of “C106-1” is the individual similarity ISM (C106-11), (C106-22), and (C106-33).
  • the breakdown information displayed in the similar case list 74 is a correspondence relationship between the region of interest ROI and the case lesion CL for calculating each individual similarity ISM.
  • the individual similarity ISM (C106-11) is the region of interest ROI of No1 and the case lesion CL of No1 of “C106”, and the individual similarity ISM (C106-22) is the region of interest ROI of No2 and “C106 And the individual similarity ISM (C106-33) is the region of interest ROI of No3 and the case lesion CL of No3 of “C106”. Since the type is also displayed in the correspondence relationship, it is possible to confirm what type of combination the combination pattern of the individual similarity degree ISM is.
  • the display items of the similar case list 74 include lesion images of the case lesions CL of No1, No2, and No3.
  • the lesion image is read from the ISM table 71, for example.
  • the examination image 19 including the region of interest ROI of No1 to No3 is displayed.
  • the top six cases are displayed.
  • the output control unit 69 creates XML data for WEB distribution by using a markup language such as XML (Extensible Markup Language) with respect to the search result display screen 76 thus created.
  • XML Extensible Markup Language
  • the WEB browser reproduces the search result display screen 76 based on the XML data and displays it on the display unit 48A. Thereby, the search result display screen 76 including the similar case list 74 is provided for viewing by the doctor.
  • the doctor of the medical department 10 accesses the inspection image DB server 15 at the medical department terminal 11 and acquires the inspection data 21 of the inspection requested from the inspection department 12 (S1100).
  • the medical department terminal 11 displays the examination data 21 on the display unit 48A (S1200).
  • the inspection image 19 included in the acquired inspection data 21 is displayed on the inspection image display screen 52 shown in FIG.
  • the region of interest ROI is designated in the examination image 19 by the doctor through the examination image display screen 52.
  • the medical department terminal 11 accepts designation of a plurality of regions of interest ROI by a doctor designation operation (S1300).
  • the region of interest ROI when diagnosing a non-cancer disease such as tuberculosis or diffuse pancreatic schizophrenia, a plurality of regions of interest ROI are included so that different types of target lesions OL of different types of lesions are included. Is specified.
  • the similar case search button 52E is operated. Thereby, the medical department terminal 11 receives a search instruction (S1400).
  • the search request issuing unit 54 issues a similar case search request with the examination image 19 and the region information attached thereto, and this is transmitted to the similar case search server 17 (S1500).
  • the request reception unit 61 receives the request (S2100). Then, the feature amount calculation unit 62 calculates a feature amount for each region of interest ROI based on the inspection image 19 and the region information of the region of interest ROI (S2200). Thereafter, a similar case search is executed (S2300).
  • the individual similarity calculation unit 65 reads one case data 24 from the case DB server 16 (S2310).
  • the individual similarity calculation unit 65 calculates the individual similarity ISM by associating the plurality of regions of interest ROI in the examination data 21 with the case lesions CL included in one case data 24 on a one-to-one basis. (S2320).
  • the individual similarity ISM is calculated for each case lesion CL.
  • the individual similarity calculation unit 65 records the calculated individual similarity ISM in the ISM table 71, and creates the ISM table 71 for each region of interest ROI (S2330). After such processing is performed for one case data 24, the same processing is performed for the next case data 24. The same processing is repeated until calculation of the individual similarity ISM and creation of the ISM table 71 are completed for a plurality of case data 24, for example, all case data 24 in the case DB 23 (N in S2340).
  • the total similarity calculation unit 66 creates the TSM table 72 based on the plurality of ISM tables 71 created for each region of interest ROI (S2350). In creating the TSM table 72, the total similarity calculation unit 66 first determines a case to be searched based on the number of types of case lesions CL for each case (S2350). As a result of this determination, the number of types of case lesions CL is the region of interest ROI as shown in FIG. 34 where the case ID is “C104” and there is only one type of case lesion CL or “C105” shown in FIG. Cases less than the number of cases are determined not to be searched.
  • the total similarity calculation unit 66 creates a completely different combination of individual similarity degrees ISM for each case for the search target case (S2352), and the total similarity only for the created completely different combination. TSM is calculated (S2353).
  • the total similarity calculation unit 66 creates 18 completely different combinations shown in FIG. 29 for the case of “C102”, and determines 18 combinations as calculation targets.
  • 18 total similarity TSMs (C102-1 to C102-18) are calculated for each determined combination of calculation targets.
  • the total similarity calculation unit 66 determines that the completely same-type combination or the partially same-type combination shown in FIG. 30 is not subject to calculation, and does not calculate the total similarity TSM.
  • the total similarity calculation unit 66 creates 24 completely different combinations as shown in FIG. 31, and 24 total similarity TSMs (24 corresponding to the 24 combinations). C106-1 to C106-24) are calculated.
  • the total similarity calculation unit 66 determines that the completely same type combination or the partially same type combination shown in FIG. 32 is not subject to calculation, and does not calculate the total similarity TSM. Then, the total similarity calculation unit 66 records the calculated total similarity TSM in the TSM table 72 for each case (S2354).
  • the similar case search unit 67 creates a similar case list 74 based on the created TSM table 72 (S2360). In creating the list, first, the representative value determining unit 67B determines a representative value from a plurality of total similarity TSM for each case in the TSM table 72 as shown in FIG. 38 (S2361), A TSM table 72 is created in which only the total similarity TSM of case representatives is extracted. Thereafter, the similar case search unit 67 sorts the cases in descending order of the total similarity TSM in the TSM table 72 (S2362). As a result, a similar case having a high similarity is extracted at the top of the TSM table 72.
  • the list creation unit 67A extracts cases up to a predetermined rank based on the TSM table 72, and creates a similar case list 74 in which similar cases are arranged in descending order of the overall similarity TSM (S2363).
  • the output control unit 69 converts the search result display screen 76 including the similar case list 74 created as the search result by the list creation unit 67A into XML data for delivery, and delivers it to the clinical department terminal 11. (S2400).
  • the medical department terminal 11 receives the XML data including the similar case list 74 (S1600), reproduces the search result display screen 76 (FIG. 30) based on the XML data, and displays it on the display unit 48A (S1700). .
  • each of a plurality of regions of interest ROI including each target lesion OL and a plurality of case lesions CL is obtained.
  • an individual similarity ISM with a case lesion CL is calculated for each region of interest ROI, and a similar case search is performed by calculating an overall similarity TSM based on the calculated individual similarity ISM.
  • the total similarity TSM is an index for evaluating a case having a plurality of individual similarity ISMs having a high average as a case having a high similarity.
  • Tuberculosis that requires attention to three types of target lesions OL: hollow shadows, point-like shadows, and ground glass shadows, and diffuse generalization that requires attention to two types of target lesions OL, bronchial abnormal shadows and point-like shadows
  • the appearance of multiple target lesions OL may be the specific basis for the disease.
  • the present invention is useful in the diagnosis of non-cancer diseases in which attention must be paid to the feature quantities of a plurality of regions of interest ROI.
  • the processing time and processing load of the total similarity calculation unit 66 are reduced, and the search time for similar cases can be shortened.
  • the case where the total similarity calculation unit 66 determines that the total similarity TSM is not calculated is also the individual similarity ISM calculation target. However, the total similarity calculation unit 66 does not calculate the case. About the case to determine, it is good also as a calculation object of individual similarity degree ISM.
  • the case lesion CL such as the case of “C103” having two types and four case lesions CL and the case of “C104” having one type and four case lesions CL shown in FIG.
  • the number of types is less than the number of regions of interest ROI, as shown in FIG.
  • the individual similarity calculation unit 65 determines that the individual similarity ISM is not a search target when calculating the individual similarity ISM, and does not calculate the individual similarity ISM.
  • the individual similarity calculation unit 65 determines whether or not to make a search target based on the number of types of case lesions CL and the number of regions of interest ROI. In this way, useless processing time is reduced, and search time can be further shortened.
  • the total similarity TSM is the sum total of a plurality of individual similarities ISM, but may be a multiplied value.
  • the representative value determination unit 67B determines a representative value for each case from a plurality of total similarity TSM calculated for each case, and searches for similar cases based only on the representative value. Is going.
  • the representative value By determining the representative value, there is an effect that the processing time can be shortened because data handled in the search processing is reduced, for example, the number of total similarity TSMs recorded in the TSM table 72 is reduced.
  • Diagnosis results such as doctor's findings regarding the case lesion CL described in the interpretation report exist for each case.
  • diagnosis based on similar cases can be performed appropriately. It can also be performed efficiently.
  • a similar case search may be performed without determining the representative value.
  • a plurality of total similarity TSMs of the same case may be displayed in the similar case list 74. Since the combination pattern of the individual similarity degree ISM used as a calculation element differs in these some total similarity TSM, it can refer with respect to one case changing a viewpoint.
  • each individual similarity ISM which is a calculation element of the total similarity TSM, grasps which one of a plurality of regions of interest ROI and a plurality of case lesions CL are associated with each other and which is calculated. can do.
  • the examination image 19 and the lesion image of the case lesion CL are also displayed in the similar case list 74, it is easy to compare and refer to the image patterns, and it is easy to intuitively determine the similarity between the image patterns.
  • the value of the individual similarity ISM that is a calculation element of the total similarity TSM may be displayed. If the value of the individual similarity ISM as a breakdown is displayed in addition to the value of the total similarity TSM, the individual similarity ISM is high or low in the correspondence between each region of interest ROI and each case lesion CL It is convenient because you can check which one is. For example, when a doctor wants to emphasize one of a plurality of regions of interest ROI, an appropriate similar case can be found by looking at the value of the individual similarity ISM of the region of interest ROI to be emphasized.
  • the overall similarity TSM is an index for evaluating a case having a plurality of individual similarity ISMs having a high average as a case having a high similarity. Therefore, in the similar case list, the rank of cases having a high average value of a plurality of individual similarity degrees ISM is high, and for cases having a low average value, the rank is low even if one individual similarity degree ISM is prominently high. .
  • the doctor's subjective evaluation regarding the similarity may be greatly influenced by the doctor's impression of the specific case lesion CL and the specific region of interest ROI rather than the average value. There may be a discrepancy between the subjective evaluation and the objective evaluation (ranking) based on the total similarity TSM.
  • the doctor can confirm the individual similarity ISM, You can verify your subjective evaluation. If the value of the individual similarity ISM is displayed, the doctor considers the value of the individual similarity ISM so that the doctor can evaluate the objective evaluation based on the total similarity TSM based on the similar case list 75. While making corrections based on subjective evaluations, it is possible to find appropriate similar cases suitable for individual diagnoses.
  • the calculation target of the total similarity TSM is set to a completely different type of combination of a plurality of individual similarities ISM. Limited to combinations. Therefore, even if a plurality of regions of interest ROI each including the same type of target lesion OL are specified (for example, “B: two cavities”), the same type of case lesions are applied to the plurality of regions of interest ROI of the same type. Since the combination of the individual similarity ISM calculated in correspondence with CL is excluded from the calculation target of the total similarity TSM, it is not extracted as a similar case in the similar case list.
  • the total similarity TSM is calculated.
  • pathology CL contained in the case is calculated about one case instead of all cases.
  • a completely different combination of the calculated individual similarities ISM is created, and the total similarity TSM is calculated.
  • each procedure of calculating the individual similarity ISM, creating a completely different combination, and calculating the total similarity is executed for the next case. This is repeated for all cases.
  • the total similarity TSM is the sum of a plurality of types of individual similarity ISM corresponding to a plurality of regions of interest ROI, and therefore the total similarity TSM cannot be calculated when the number of types is “1”. .
  • the number of types of case lesions CL is preferably greater than or equal to the number of regions of interest ROI. This is because when the similar case search is used for diagnosis of a non-cancer disease, the number of regions of interest ROI corresponds to the number of types of target lesions OL. However, in the case DB 23, there are cases where the number of cases accumulated with a large number of types of case lesions CL is small. In this case, if all cases of the number of types less than the number of regions of interest ROI are excluded from the search target, the search target is too small and an appropriate similar case may not be searched. Even in cases where the number of types of case lesions CL is less than the number of regions of interest ROI, there are cases where the individual case lesions CL have valuable values for diagnosis.
  • the individual similarity calculation unit 65 calculates the individual similarity ISM by associating two types of case lesions CL with the regions of interest ROI No1 to No3.
  • the hatched column shown in FIG. 35 indicates a location where the individual similarity ISM is not calculated because there is no corresponding case lesion CL.
  • the total similarity calculation unit 66 creates completely different combinations of individual similarities ISM calculated by associating two types of four case lesions CL with respect to the three regions of interest ROI, and the combination is determined as the total similarity. It is determined as a TSM calculation target. In the case of “C103”, the completely different combinations are “B1, F”, “B2, F”, and the like.
  • the similar case search unit 67 creates a similar case list 81 shown in FIG.
  • the total similarity TSM (1.67) of the case “C103” is the sum of the two types of individual similarity ISM and the total of the three types of individual similarity ISM. It cannot be compared with the value of the total similarity TSM of other cases. Therefore, the similar case search unit 67 performs normalization so that the total similarity TSMs having different numbers of individual similarities ISM as calculation elements can be compared. Each case is ranked based on a normalized value that is a normalized similarity. Normalization is, for example, a process of dividing the total similarity TSM by the number of individual similarities ISM.
  • the similar case search server 17 may perform a search again according to a change in the designated number of regions of interest ROI.
  • the first similar case search request issues a designated number of regions of interest ROI of 2.
  • the similar case search server 17 performs a similar case search based on the designation, and distributes the similar case list 74 as a search result.
  • the doctor sees the search result, newly designates a region of interest ROI, and issues a second similar case search request.
  • the similar case search server 17 performs a similar case search based on the request including the added designation, and distributes the similar case list 74 as a search result. In this way, the doctor can change the search conditions as necessary while looking at the search results, and therefore, it is easy to find an appropriate similar case.
  • the designated number of regions of interest ROI may be decreased as well as increased.
  • the similar case search server 17 may store data created as intermediate processes and processing results of similar case searches such as the ISM table 71 and the TSM table 72 in preparation for re-search. If such data is used in the re-search, the search time can be shortened.
  • the total similarity TSM has been described as an example of a simple total value of a plurality of individual similarities ISM, but weighting is performed based on the value of the individual similarity ISM as a calculation element. Processing may be performed.
  • the weighting process for example, when the individual similarity ISM is equal to or greater than a threshold value, the total similarity TSM is increased by multiplying by a positive weighting coefficient, or conversely, the individual similarity ISM is less than the threshold value. In some cases, a negative weighting coefficient is multiplied to lower the overall similarity TSM. If there is even one case lesion CL that closely resembles the region of interest ROI, the case may be valuable in diagnosis. Positive weighting makes it easier to find such cases as similar cases.
  • the type of the case lesion CL is determined in advance, and the determined type information is stored in the case DB 23.
  • the feature amount at the time of the search is described.
  • Type determination based on CAC may be performed. Of course, since this takes time for the search process, it is preferable to store the type information in advance.
  • the similar case search is performed by calculating the individual similarity ISM in association with the case lesion CL without determining the type of the target lesion OL included in the region of interest ROI.
  • the lesion type is determined for the target lesion OL and the case lesion CL included in the region of interest ROI, and only the same type of lesions are associated with each other.
  • a similar case search may be performed by calculating the individual similarity ISM.
  • the lesion pattern is typically distinguished by the type of lesion. Therefore, the type of lesion can be determined based on the feature amount at the stage of calculating the feature amount.
  • the third embodiment is a form using such lesion type determination.
  • the similar case search server 17 includes a lesion type determination unit 86.
  • the lesion type determination unit 86 determines the type of the target lesion OL included in the region of interest ROI based on the feature amount RAC of the region of interest ROI calculated by the feature amount calculation unit 62.
  • the lesion type determination unit 86 includes, for example, a lesion type corresponding to a discriminator indicating the maximum discriminator output value among the discriminator output values output from the respective discriminators 62A to 62H, as a target included in the region of interest ROI. Determined as the type of lesion OL. In this example, since the discriminator output value of the discriminator 62B of “B: cavity” is the maximum, the type of the target lesion OL is determined to be “B: cavity”.
  • the individual similarity calculation unit 65 when calculating the individual similarity ISM between the region of interest ROI and each case lesion CL, the individual similarity calculation unit 65 only determines the individual similarity ISM for the same kind of lesions. And the individual similarity ISM is not calculated for those having different types of lesions.
  • the type of the region of interest ROI of No1 is “B: Cavity”
  • the individual similarities to the case lesion CL of No3 whose type is “B: Cavity” from the cases of “C101”. Only the degree ISM is calculated.
  • a plurality of case lesions CL of the same type as the region of interest ROI are registered in one case, a plurality of individual similarities ISM are calculated. If no case lesion CL of the same type as the region of interest ROI has been registered, the individual similarity ISM is not calculated for that case.
  • the calculation processing time of the individual similarity calculation unit 65 can be reduced.
  • the number of individual similarities ISM is reduced, as shown in FIG. 50, the number of completely different combinations to be calculated for the total similarity TSM is also reduced.
  • the total similarity calculation unit 66 calculates three total similarity TSMs C102-1 to C102-3 corresponding to this combination.
  • the calculation time of the individual similarity ISM and the total similarity TSM is shorter than that of the first embodiment in which the individual similarity ISM is calculated without distinguishing between the types of lesions. Further, since the size of the ISM table 71 and the TSM table 72 is reduced, the work area of the memory can be reduced. Therefore, the load on the CPU 41B of the similar case search server 17 is reduced, and the search time can be shortened.
  • the lesion type is determined in advance and the individual similarity degree ISM is performed only for the same type
  • the case lesion CL to be searched for as a similar case is omitted.
  • leakage There are also concerns about leakage.
  • a determination biased to one of the plurality of target lesions OL is made by the type determination. . Therefore, it is preferable to implement the fifth embodiment after determining the accuracy of determining the type of lesion.
  • the fourth embodiment shown in FIGS. 51 and 52 is a form in which the feature amount of the region of interest ROI is calculated not in the similar case search server 17 but in the medical department terminal 11.
  • the feature amount of the region of interest ROI may be calculated by the clinical department terminal 11 as in the third embodiment.
  • the similar case search server 17 is not provided with the feature amount calculation unit 62, and FIG.
  • the components other than the feature amount calculation unit 62 such as the individual similarity calculation unit 65, the total similarity calculation unit 66, and the similar case search unit 67 shown in FIG.
  • the clinical department terminal 11 is provided with a feature quantity calculation unit 88 similar to the feature quantity calculation unit 62.
  • the feature amount calculation unit 88 is realized, for example, when the CPU 41A executes software installed in the medical department terminal 11.
  • the feature amount calculation unit 88 calculates the feature amount RAC based on the inspection data 21 including the inspection image 19 and the region information of the region of interest ROI input through the GUI control unit 53.
  • the search request issuing unit 54 issues a similar case search request by attaching the image corresponding to the region of interest ROI and the calculated feature amount RAC.
  • the similar case search request is transmitted from the clinical department terminal 11 to the similar case search server 17.
  • the similar case search server 17 performs a similar search based on the received similar case search request and distributes the search result to the medical department terminal 11.
  • the request reception unit 61 of the similar case search server 17 functions as a feature amount acquisition unit.
  • the similar case search device of the present invention has been described in the form of the similar case search server 17 that performs a similar case search based on a request from the medical department terminal 11.
  • the clinical department terminal 11 may be provided with a similar case search function such that the clinical department terminal 11 accesses the case DB server 16 and searches for a similar case without using it.
  • the medical department terminal 11 serves as a similar case search device.
  • the similar case search server 17 and the case DB server 16 are configured as separate servers. However, they may be integrated into a single server. As described above, a plurality of functions may be integrated into one server, or the servers may be separated for each function.
  • the hardware configuration of the computer system can be variously modified.
  • the similar case search server 17 may be configured by a plurality of server computers separated as hardware for the purpose of improving processing capability and reliability.
  • the hardware configuration of the computer system can be appropriately changed according to required performance such as processing capability, safety, and reliability.
  • programs such as the case DB 23 and the AP 50 can of course be duplicated or distributed and stored in a plurality of storage devices for the purpose of ensuring safety and reliability. .
  • the similar case search server 17 has been described as being used in one medical facility, but may be in a form in which a plurality of medical facilities can be used.
  • the similar case search server 17 is configured such that a client terminal installed in one medical facility such as the medical department terminal 11 is connected to be communicable via a LAN, and a request from the client terminal
  • the application service related to similar case search is provided based on the above.
  • the similar case search server 17 communicates with client terminals installed in a plurality of medical facilities via a WAN (Wide (Area Network) such as the Internet or a public communication network, for example. Connected as possible.
  • the similar case search server 17 receives the request from the client terminal of a some medical facility, and provides the application service regarding a similar case search with respect to each client terminal.
  • WAN Wide (Area Network)
  • the installation location and operating entity of the similar case search server 17 may be, for example, a data center different from the medical facility or one of a plurality of medical facilities.
  • WAN Virtual Private Network
  • HTTPS Hypertext Transfer Protocol Secure
  • the present invention is not limited to the above-described embodiments, and various configurations can be adopted without departing from the gist of the present invention.
  • CT, MRI, and a simple X-ray image are taken as examples of inspection images, but the present invention can also be applied to inspection images taken with other modalities such as mammography and an endoscope.
  • the above-described various embodiments and various modifications can be appropriately combined.
  • the present invention extends to a storage medium for storing a program in addition to the program for realizing the present invention.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Data Mining & Analysis (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Databases & Information Systems (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • General Engineering & Computer Science (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Library & Information Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Software Systems (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Quality & Reliability (AREA)
  • Medical Treatment And Welfare Office Work (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Processing Or Creating Images (AREA)
  • Image Analysis (AREA)

Abstract

 複数の関心領域を基に適切な類似症例を迅速に検索する類似症例検索装置、方法及びプログラムを提供する。 類似症例検索サーバ(17)は、個別類似度算出部(65)、総合類似度算出部(66)、類似症例検索部(67)を備える。個別類似度算出部(65)は、関心領域と症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた関心領域の特徴量と症例病変の特徴量とを比較して、関心領域毎の個別類似度をそれぞれ算出する。総合類似度算出部(66)は、複数の関心領域に対してそれぞれ算出した複数の個別類似度に基づいて、総合類似度を算出するものであり、同一症例内の種類が異なる複数の症例病変と複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、個別類似度を用いて総合類似度を算出する。類似症例検索部(67)は、完全異種の組み合わせに対応する総合類似度に基づいて類似症例を検索する。

Description

類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム
 本発明は、類似症例検索装置、類似症例検索方法、及び類似症例検索プログラムに関する。
 医療分野において、検査画像に基づいて、検査画像に類似する過去の症例を検索する類似症例検索装置が知られている(例えば特開2010-237930号公報、特開2012-118583号公報(米国公開公報US2012/134555号)参照)。検査画像は、例えば、断層撮影を行うCT(Computed Tomography)装置や単純透視画像を撮影する一般X線撮影装置などのモダリティで撮影された画像であり、患者の疾患の特定など患者の診断を行うために用いられる。一般X線撮影装置による1回の検査では、検査画像が1枚だけ撮影される場合もあれば、複数枚撮影される場合もある。また、CT装置による1回の検査では、複数枚の断層画像(スライス画像)が取得される。そのため、1件の検査データには、1枚以上の検査画像が含まれる。症例は、過去の検査データの集積によって作成される場合が多いため、1件の症例のデータにも、1枚以上の症例画像が含まれる。
 類似症例検索を行う場合には、まず、医師などのユーザによって、検査画像内に関心領域が指定される。関心領域は、検査画像内において医師が特に関心を示す領域であり、診断対象となる病変を含む領域である。類似症例検索装置は、検査画像内に指定された1つの関心領域の特徴を数値化した特徴量と、症例画像内に存在する1つの病変の特徴を数値化した特徴量とを比較して両者の類似度を判定する。ここで、便宜上、検査画像の関心領域内に含まれる病変を対象病変といい、症例画像に含まれる病変を症例病変という。そして、類似症例検索装置は、複数の症例を記憶した症例データベースの中から、関心領域に類似する症例病変を含む症例を検索する。
 特開2010-237930号公報には、対象病変を含む関心領域の指定の仕方にはユーザ毎に個人差があり、個人差に起因して検索結果にばらつきが生じるという検索ばらつきを抑制することを目的する発明が開示されている。具体的には、同じ1つの対象病変を含む領域を関心領域として指定する場合でも、指定するユーザの個人差などが原因で、関心領域の指定の仕方によっては指定した領域の形状や大きさが変化するため、特徴量が変化してしまう場合がある。特徴量が変化すると、類似度も変化するため、ユーザ毎に検索結果が変わるという検索ばらつきが生じてしまう。特開2010-237930号公報は、こうした検索ばらつきを抑制するために、1例として、1つの対象病変について指定の仕方が異なる複数の関心領域について、関心領域毎に特徴量を算出し、算出した複数の関心領域の特徴量の平均値に基づいて類似度を算出して、類似画像検索を行っている。これによれば、ユーザ毎の個人差に起因する検索ばらつきを抑制することができる。
 特開2012-118583号公報(米国公開公報US2012/134555号)は、ユーザが類似していると感じる主観に対して、より適合的な検索結果を出力する技術に関するものである。具体的には、同種の対象病変が複数の検査画像に対して存在する場合に、関心領域の指定に際して、ユーザが類似していると感じる同種の複数の対象病変を含む関心領域を、同種グループとして1つのグループにまとめる。そして、1件の検査データにおいて、同種グループに属する複数の対象病変の特徴量をすべて包含する特徴量の範囲を求め、この特徴量の範囲を検索条件として、類似症例検索を行う。同種グループの特徴量の範囲は、ユーザが主観的に類似すると感じる範囲と一致すると考えられるため、検索結果もユーザの主観に対して、より適合的なものとすることができる。
 ところで、疾患によっては、検査画像において複数の対象病変が現れていることが、疾患特定の根拠となる場合がある。例えば、結核の場合には、検査画像において、空洞影、点状影、すりガラス影という3種の対象病変が現れていることを以て疾患が特定される場合があり、びまん性汎細気支炎では、気管支異常影と点状影の2種の対象病変が現れていることを以て疾患が特定される場合がある。がんの場合には、単一の対象病変での類似検索でよかったが、がん以外の非がん疾患である場合には、これら複数の対象病変での類似症例検索が必要である。
 特開2010-237930号公報及び特開2012-118583号公報(米国公開公報US2012/134555号)に記載の従来の類似症例検索装置は、いずれも、検査画像に含まれている1つの対象病変に着目して、着目した1つの対象病変を含む関心領域の特徴量に基づいて類似症例を検索するものであり、検査画像に含まれている複数の対象病変のそれぞれに着目することは考慮されていなかった。
 上述したように、特開2010-237930号公報は、関心領域毎に特徴量を算出しているが、複数の関心領域は、指定の仕方が異なるだけで、同じ1つの対象病変の関心領域であり、異なる対象病変のそれぞれの関心領域の特徴量に着目して類似症例を検索することについては開示が無い。また、特開2012-118583号公報(米国公開公報US2012/134555号)は、複数の検査画像に含まれる複数の対象病変に対して、同種の対象病変がグループ化された1つの同種グループについて1つの検索条件を作成し、作成した検索条件で類似症例を検索している。特開2012-118583号公報(米国公開公報US2012/134555号)は、いわば、同種の複数の対象病変を含むそれぞれの関心領域に共通する特徴量をユーザの嗜好に合うように算出しているのみで、複数の対象病変のそれぞれの関心領域の特徴量に着目して類似症例を検索することについては開示が無い。
 特開2010-237930号公報及び特開2012-118583号公報(米国公開公報US2012/134555号)のように、1つの関心領域の特徴量に着目する技術では、関心領域が複数有る場合には、適切な類似症例を検索することができない。そこで、発明者は、関心領域が複数有る場合においては複数の関心領域に着目して、複数の症例病変が登録されている症例の中から類似する症例を検索する技術を検討している。しかし、複数の関心領域に着目するためには、複数の関心領域と複数の症例病変のそれぞれを対応させて両者の特徴量を比較しなければならない。その場合には、対応させる組み合わせパターンが多くなるため、検索処理に時間が掛かるという新たな課題が生じることが分かってきた。関心領域の数が増えるほど、こうした問題は顕著となる。
 本発明は、関心領域が複数有る場合でも、適切な類似症例検索を短時間で検索することが可能な類似症例検索装置、類似症例検索方法、及び類似症例検索プログラムを提供することを目的とする。
 本発明の類似症例検索装置は、症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する装置であって、特徴量取得部、個別類似度算出部、総合類似度算出部、類似症例検索部を備える。特徴量取得部は、検査画像を1枚以上含む検査データにおいて指定され、検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、関心領域毎の特徴量を取得する。個別類似度算出部は、関心領域と、症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた関心領域の特徴量と症例病変の特徴量とを比較して、関心領域毎の個別類似度をそれぞれ算出する。総合類似度算出部は、複数の関心領域に対してそれぞれ算出した複数の個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の症例病変と複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、個別類似度を用いて総合類似度を算出する。類似症例検索部は、完全異種の組み合わせに対応する総合類似度に基づいて、類似症例を検索する。
 ここで、症例画像内に存在する複数件の症例病変は、1枚の症例画像内に複数存在する場合と、例えば2枚の症例画像のそれぞれに、症例病変が1つずつ存在する場合のように、複数枚の症例画像内に存在する症例病変の合計が複数件になる場合を含む。
 総合類似度算出部は、完全異種の組み合わせに関して、関心領域の個数と、症例病変の種類数及び個数とに応じた数の組み合わせを作り、完全異種の組み合わせ毎に総合類似度を算出することが好ましい。具体的には、総合類似度算出部は、1件の症例毎に、その症例に含まれる症例病変についての個別類似度を求め、求めた個別類似度についての完全異種の組み合わせを作って総合類似度を求めてもよい。また、最初に、症例全件の症例病変についての個別類似度をそれぞれ求め、その後に、症例毎に、個別類似度についての完全異種の組み合わせを作ってもよい。いずれにしても、個別類似度ISMを算出後に、個別類似度ISMについての完全異種の組み合わせを作る。
 同一症例内に含まれる複数の異種の症例病変のうち少なくとも1種に関して、同種の症例病変が複数個存在する場合には、総合類似度算出部は、同種の複数個の症例病変を区別して、区別した同種の症例病変毎に完全異種の組み合わせを作ることが好ましい。
 類似症例検索部は、総合類似度に基づいて、複数件の類似症例に関する情報をリスト化した類似症例リストを作成することが好ましい。類似症例リストにおいて、類似症例は、総合類似度順にソートされることが好ましい。
 類似症例リストの表示項目には、総合類似度の値と総合類似度に関する内訳情報とを含み、内訳情報は、個別類似度を算出するための関心領域と症例病変との対応関係を含むことが好ましい。類似症例リストには、総合類似度の値に加えて、総合類似度の算出要素となる複数の個別類似度の値が表示されることが好ましい。類似症例リストには、関心領域と症例病変のそれぞれの画像が表示されることが好ましい。
 症例データベースには、症例病変の種類の情報が格納されていることが好ましい。総合類似度は、完全異種の組み合わせに含まれる複数の個別類似度の総和であることが好ましい。
 個別類似度算出部は、関心領域毎に、複数の症例病変を対応させて算出した複数の個別類似度を記録した個別類似度テーブルを作成してもよい。総合類似度算出部は、関心領域毎の複数の個別類似度テーブルのそれぞれから、1つずつ個別類似度を読み出して、読み出した複数の個別類似度を要素として、完全異種の組み合わせを作ってもよい。
 個別類似度算出部は、症例のうち、少なくとも症例病変の種類が複数の症例については個別類似度の算出対象とし、症例病変の種類が1種類の症例については個別類似度の算出対象外とすることが好ましい。
 総合類似度算出部は、総合類似度を算出する要素となる個別類似度の値に応じて、総合類似度に対して重み付け処理を行ってもよい。重み付け処理は、個別類似度が閾値以上の場合に、総合類似度を高くする処理であることが好ましい。
 類似症例検索部は、症例病変の種類数が、関心領域の数未満の症例については、検索の対象外としてもよい。また、類似症例検索部は、症例病変の種類数が、関心領域の数未満の症例についても、症例病変の種類数が複数有る場合には、検索の対象としてもよい。
 1つの関心領域と、1件の症例に含まれる複数の症例病変とを対応させることにより、複数の総合類似度が算出された場合に、複数の総合類似度の中から、1つの代表値を判定する代表値判定部を有しており、類似症例検索部は、代表値に基づいて、類似症例を検索することが好ましい。
 本発明の類似症例検索方法は、症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する方法において、特徴量取得ステップ、個別類似度算出ステップ、総合類似度算出ステップ、類似症例検索ステップを含む。特徴量取得ステップは、検査画像を1枚以上含む検査データにおいて指定され、検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、関心領域毎の特徴量を取得する。個別類似度算出ステップは、関心領域と、症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた関心領域の特徴量と症例病変の特徴量とを比較して、関心領域毎の個別類似度をそれぞれ算出する。総合類似度算出ステップは、複数の関心領域に対してそれぞれ算出した複数の個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の症例病変と複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、個別類似度を用いて総合類似度を算出する。類似症例検索ステップは、完全異種の組み合わせに対応する総合類似度に基づいて、類似症例を検索する。
 総合類似度算出ステップは、症例毎に症例病変についての完全異種の組み合わせを作り、作った組み合わせに対応する個別類似度を求めてもよいし、最初に症例全件の個別類似度を求め、その後に症例毎に個別類似度の完全異種の組み合わせを作ってもよい。いずれにしても、個別類似度ISMを算出後に、個別類似度ISMについての完全異種の組み合わせを作る。
 本発明の類似症例検索プログラムは、症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する処理をコンピュータに実行させるプログラムにおいて、特徴量取得ステップ、個別類似度算出ステップ、総合類似度算出ステップ、類似症例検索ステップを含む。特徴量取得ステップは、検査画像を1枚以上含む検査データにおいて指定され、検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、関心領域毎の特徴量を取得する。個別類似度算出ステップは、関心領域と、症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた関心領域の特徴量と症例病変の特徴量とを比較して、関心領域毎の個別類似度をそれぞれ算出する。総合類似度算出ステップは、複数の関心領域に対してそれぞれ算出した複数の個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の症例病変と複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、個別類似度を用いて総合類似度を算出する。類似症例検索ステップは、完全異種の組み合わせに対応する総合類似度に基づいて、類似症例を検索する。
 本発明は、関心領域毎の特徴量と症例病変の特徴量とを比較して、関心領域毎の個別類似度を算出し、同一症例内の種類が異なる複数の症例病変と複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、個別類似度を用いて総合類似度を算出し、算出した総合類似度に基づいて類似症例を検索するため、関心領域が複数有る場合でも、適切な類似症例検索を短時間で検索することが可能な類似症例検索装置、類似症例検索方法、及び類似症例検索プログラムを提供することができる。
類似症例検索サーバを含む医療情報システムを示す構成図である。 複数枚の検査画像で構成される検査データの概略を示す説明図である。 1枚の検査画像で構成される検査データの概略を示す説明図である。 診療科、検査科、検査画像DBサーバ、及び症例DBサーバの機能を説明する説明図である。 症例DBの説明図である。 病変の種類の画像パターンを示す説明図である。 類似症例検索サーバの機能概要の説明図である。 各DBサーバや各端末を構成するコンピュータを示す構成図である。 診療科端末の概略を示す構成図である。 関心領域を指定する検査画像表示画面の説明図である。 関心領域の指定方法について、図9とは別の例を説明図である。 類似症例検索サーバの概略を示す構成図である。 関心領域の特徴量の説明図である。 特徴量算出部の構成図である。 症例病変の特徴量の説明図である。 個別類似度算出部の説明図である。 個別類似度算出方法の説明図である。 症例病変が3種3個の場合の個別類似度の説明図である。 症例病変が3種5個の場合の個別類似度の説明図である。 症例病変が2種4個の場合の個別類似度の説明図である。 症例病変が1種4個の場合の個別類似度の説明図である。 症例病変が3種5個の場合の個別類似度の説明図である。 個別類似度テーブルの概略を示す説明図である。 関心領域毎に作られる個別類似度テーブルの説明図である。 総合類似度テーブルの説明図である。 総合類似度の算出方法を示す説明図である。 図19の例の完全異種の組み合わせを示す説明図である。 図27の完全異種の組み合わせに対応する総合類似度の説明図である。 3種5個の症例病変の場合の完全異種の組み合わせを示す説明図である。 図29の例において算出対象外となる組み合わせの説明図である。 図22の例の完全異種の組み合わせを示す説明図である。 図31の例において算出対象外となる組み合わせの説明図である。 検索対象外となる症例の説明図である。 検索対象外となる別の症例の説明図である。 検索対象と対象外となる症例を説明する表である。 総合類似度を算出した総合類似度テーブルを示す説明図である。 症例毎に総合類似度の代表値を判定する方法を示す説明図である。 症例毎の代表値をテーブルに並べた説明図である。 症例毎の代表値をソートしたテーブルを示す説明図である。 類似症例リストを表示した画面を示す説明図である。 類似症例画像検索装置の手順を示すフローチャートである。 類似症例検索の手順を示すフローチャートである。 総合類似度と個別類似度との両方を表示する例を示す説明図である。 第2実施形態の説明図である。 第2実施形態における2種2個の症例の完全異種の組み合わせの説明図である。 第2実施形態の検索結果表示画面の説明図である。 第3実施形態の説明図である。 病変種類判定部の説明図である。 第3実施形態の個別類似度算出部の説明図である。 第3実施形態の完全異種の組み合わせテーブルの説明図である。 特徴量算出部を診療科端末に備えた別の例を示す説明図である。 図51で説明した診療科端末と類似症例検索サーバとの説明図である。
[第1実施形態]
 図1に示す医療情報システム9は、病院などの医療施設に構築される。医療情報システム9は、診療科10に配置される診療科端末11、検査科12に設置されるモダリティ(医用画像撮影装置)13、及びオーダ管理端末14、検査画像データベース(以下「DB」と称す)サーバ15、症例DBサーバ16、及び類似症例検索サーバ17を有する。これらは、ネットワーク18を通じて通信可能に接続される。ネットワーク18は、例えば、院内に敷設されたLAN(Local Area Network)である。モダリティ13は、例えば、断層画像を撮影するCT(Computed Tomography)装置やMRI(Magnetic Resonance Imaging)装置、及び単純透視画像を撮影する一般X線撮影装置(DR:Digital Radiographyや、CR:Computed Radiographyなど)等を含む。
 診療科端末11は、診療科10の医師(図においてDrの符号を付す)によって操作され、電子カルテの入力や閲覧の他、検査科12に対して検査を依頼するための検査オーダの発行に利用される。また、診療科端末11は、検査科12で撮影され検査画像DBサーバ15に保存された検査画像19を表示して、医師が検査画像19を閲覧するための画像表示端末としても利用される。
 検査科12において、オーダ管理端末14は、診療科10からの検査オーダを受け付けて、受け付けた検査オーダを管理する。検査科12の技師は、検査オーダの内容に従って、モダリティ13により患者の撮影を行う。検査画像19は、1件の検査オーダに対して1枚あるいは複数枚撮影される。モダリティ13は、撮影が終了すると、撮影した検査画像19を検査画像DBサーバ15に送信する。検査が終了すると、検査科12から診療科10の医師に検査終了が通知され、合わせて、検査画像DBサーバ15内の検査画像19の保存先が通知される。診療科10の医師は、診療科端末11を通じて、検査画像DBサーバ15にアクセスして、診療科端末11で検査画像19を閲覧する。
 検査画像DBサーバ15は、複数枚の検査画像19が格納される検査画像DB20を有しており、いわゆるPACS(Picture Archiving and Communication System)サーバである。検査画像DB20は、キーワードによる検索が可能なデータベースであり、診療科端末11などの検索要求や配信要求に応じて、検索条件に合致する検査画像19や指定された検査画像19を配信する。
 図2及び図3に示すように、検査画像DB20においては、1件の検査オーダに対して、1枚以上の検査画像19を含む1件の検査データ21が対応付けて格納される。図2に示すように、CT装置やMRI装置で撮影された検査画像19は、断層画像(スライス画像とも呼ばれる)であり、1件の検査データ21には、複数枚の検査画像19が含まれる。図3に示すように、一般X線撮影装置で撮影された検査画像19の場合には、単純透視画像であり、1件の検査データ21は、検査画像19が1枚だけの場合もあれば、複数枚の場合もある。
 検査オーダには、診療科10の医師のID(Identification Data)や所属などの依頼元情報、患者情報、検査種別などが含まれる。検査画像19の画像ファイルは、画像データとDICOM(Digital Imaging and COmmunication in Medicine)ヘッダなどの付帯情報とを有する。検査オーダの情報は、検査画像19の付帯情報として格納される。また、付帯情報には、検査IDや、1枚の検査画像19にそれぞれ付与される画像IDが含まれる。図2、図3の例では、検査IDは「O901」や「O902」で、画像IDは「O901-03」、「O901-01」というように、検査IDに、1枚の検査画像19を識別するためのシリアルナンバーを加えた形態で付与されている。検査画像DBサーバ15は、こうしたDICOMタグに含まれる項目を検索キーとして、検索が可能になっている。
 類似症例検索サーバ17は、検索条件として検査画像19を受け付けて、受け付けた検査画像19に類似する症例画像22を含む症例を検索する。症例画像22は、過去に診断に利用された検査画像である。症例DBサーバ16は、複数件の症例を検索可能に格納した症例DB23を有している。類似症例検索サーバ17は、症例DBサーバ16にアクセスして、症例を1件ずつ読み出して、検索条件として受け付けた検査画像19と症例内の症例画像22の比較照合を行って、検査画像19に類似する症例を検索する。
 図4に示すように、診療科10の医師は、診療科端末11を操作して、検査画像DBサーバ15にアクセスして検査画像19を含む検査データ21をダウンロードする。診療科端末11で検査画像19は表示され、医師によって閲覧される。患者が罹患している場合には、患者の検査画像19には、疾患の症状を現す病変(対象病変OLという)が含まれている。診療科10の医師は、検査データ21に含まれる検査画像19の中から、対象病変OLを含む検査画像19を選択する。検査画像19は、診療科端末11が類似症例検索サーバ17に対して発行する類似症例検索要求に添付されて、類似症例検索サーバ17に送信される。類似症例検索サーバ17は、類似症例検索要求を受け付けると、症例DBサーバ16から検査画像19に類似する症例を検索して、検索結果を依頼元の診療科端末11に配信する。
 診療科10の医師は、検査結果に含まれる症例を確認する。症例には、症例画像22に付随して読影レポートが含まれている。医師は、読影レポートに記載された、症例画像22に対する所見などを参照して、検査画像19の疾患の特定など確定診断を下す。
 図5に示すように、症例DB23には、症例画像DB23Aと、特徴量DB23Bが設けられている。症例画像DB23Aは、症例画像22を検索可能に格納するデータベースである。症例には、1件毎に症例IDが付与される。症例IDは、検査画像19における検査IDに相当するものである。1件の症例には、1枚以上の症例画像22が含まれる。また、検査画像19と同様に、各症例画像22には、症例IDにシリアルナンバーを加えた画像IDが付与される。図5において、症例IDが「C101」の症例データ24には、例えば、60枚の断層画像が含まれている。
 症例画像22には、疾患の症状を表す病変(症例病変CL)が含まれている。1件の症例には、1つ以上の症例病変CLが登録されている。本例では、症例ID「C101」には、No1~No3の3個、「C102」には2個、「C103」には1個、それぞれ症例病変CLが登録されている。症例病変CLは、過去に症例画像22が検査画像として診断に利用された際に、医師によって病変として指定された領域であり、医師による確定診断を経て、症例病変CLとして登録されたものである。症例病変CLの指定方法については、例えば、後述する関心領域ROIと同様である。
 特徴量DB23Bは、症例病変CLの画像の特徴量CACを格納するデータベースである。特徴量CACは、症例IDと病変Noからなる識別IDが付与されている。例えば、症例ID「C101」には、3個の症例病変CLがあり、各症例病変CLには、1件の症例内のシリアルナンバーであるNo1~No3がそれぞれ付与される。特徴量CACの後の数字は、症例内のシリアルナンバーに対応している。
 また、特徴量DB23Bには、各症例病変CLの病変の種類に関する情報(種類情報)が格納されている。図6に示すように、典型的な病変の画像パターンである病変パターンは、例えば、A~Hの8種類に分類される。すなわち、A:低呼吸域の異常影(気腫、気胸、ブラなど)、B:空洞、C:気管支の異常影(気管支壁肥厚、気管支拡張、牽引性気管支拡張、気管支透亮像など)、D:蜂窩、E:すりガラス状陰影、F:点状影(粒状影、TIBなど)、G:高吸収域の異常影(コンソリデーション、ノジュール、気管支粘液腺など)、H:線状・網状影である。
 症例病変CLの病変の種類は、症例病変CLの画像パターンが、これら8種類の病変パターンのどれに最も近いかによって判定される。症例病変CLの種類は、特徴量CACに基づいて判定される。特徴量CACは、後述する関心領域ROIと同様に、8種類の病変パターンに対応する8種類の判別器の出力値で構成される8次元の特徴ベクトルである。症例病変CLの種類情報は、特徴量CACとともに、特徴量DB23Bに格納される。なお、本例では、典型的な病変パターンを8種類に分類し、それに対応して病変の種類も8種類に分類する例で説明しているが、8種類未満でもよいし、8種類以上でもよい。
 図7に示すように、診療科端末11を通じて、類似症例検索要求を発行する際には、医師によって、検査画像19内の対象病変OLを含む領域が関心領域ROIとして指定される。類似症例検索要求には、対象病変OLを含む検査画像19と、指定された関心領域ROIに対応する検査画像19内の領域情報(例えば、検査画像19内の座標情報)が添付される。類似症例検索サーバ17は、類似症例検索要求を受け付けると、検査画像19の画像データと、領域情報に基づいて、関心領域ROIを特定する。そして、関心領域ROIの特徴量を算出する。特徴量を算出後、類似症例検索サーバ17は、症例DBサーバ16から、症例を1件ずつ読み出して、関心領域ROIと症例病変CLとの特徴量を比較して、類似症例を検索する。そして、類似する複数件の症例に関する情報がリスト化された類似症例リストを作成し、これを検索結果として診療科端末11に配信する。
 診療科端末11、オーダ管理端末14、検査画像DBサーバ15、症例DBサーバ16、類似症例検索サーバ17は、パーソナルコンピュータ、サーバ用コンピュータ、ワークステーションといったコンピュータをベースに、オペレーティングシステム等の制御プログラムや、クライアントプログラム又はサーバプログラム等のアプリケーションプログラムをインストールして構成される。
 図8に示すように、各DBサーバ15~17や各端末11、14を構成するコンピュータは、基本的な構成は同じであり、それぞれ、CPU(Central Processing Unit)41、メモリ42、ストレージデバイス43、通信I/F44、及び入出力部46を備えている。これらはデータバス47を介して接続されている。入出力部46は、表示部48と、キーボードやマウスなどの入力デバイス49とからなる。
 ストレージデバイス43は、例えば、HDD(Hard Disk Drive)であり、制御プログラムやアプリケーションプログラム(以下、APという)50が格納される。また、DBが構築されるサーバには、プログラムを格納するHDDとは別に、DB用のストレージデバイス43として、例えば、HDDを複数台連装したディスクアレイが設けられる。ディスクアレイは、サーバ本体に内蔵されるものでもよいし、サーバ本体とは別に設けられ、サーバ本体にケーブルやネットワークを通じて接続されるものでもよい。
 メモリ42は、CPU41が処理を実行するためのワークメモリである。CPU41は、ストレージデバイス43に格納された制御プログラムをメモリ42へロードして、プログラムに従った処理を実行することにより、コンピュータの各部を統括的に制御する。通信I/F44は、ネットワーク18との間の伝送制御を行うネットワークインタフェースである。
 診療科端末11には、AP50として、電子カルテの閲覧や編集を行う電子カルテソフトウエアや、検査画像や類似症例リストの閲覧を行うためのビューアソフトウエアなどのクライアントプログラムがインストールされている。ビューアソフトウエアは、例えば、専用ソフトウエアでもよいし、汎用的なWEBブラウザなどでもよい。
 図9に示すように、診療科端末11において、検査画像19を表示するビューアソフトウエアが起動されると、診療科端末11の表示部48Aには、GUI(Graphical User Interface)による操作機能を備えた検査画像表示画面52が表示される。診療科端末11のCPU41Aは、GUI制御部53及び検索要求発行部54として機能する。検査画像表示画面52では、検査画像19に対する関心領域ROIの指定と、類似症例検索要求の発行指示が可能になっている。GUI制御部53は、検査画像表示画面52を通じた入力デバイス49Aからの操作指示を受け付けて、受け付けた操作指示に従った画面制御を行う。また、類似症例検索要求の発行指示が入力されると、入力された発行指示がGUI制御部53から検索要求発行部54に入力される。検索要求発行部54は、類似症例検索要求に、指定された検査画像19や関心領域ROIの領域情報を添付して、類似症例検索要求を発行する。
 図10に示すように、検査画像表示画面52は、検査画像19を表示する画像表示領域52Aと、各種操作部とを有している。画像表示領域52Aには、例えば、3枚の検査画像19が並べて表示される。スクロール操作やコマ送り操作により、表示する検査画像19を切り替えることが可能である。画像表示領域52Aの上方には、検査IDを入力する入力ボックス52Bが設けられている。入力ボックス52Bに検査IDを入力すると、入力した検査IDの検査データ21が検査画像DBサーバ15からダウンロードされて、画像表示領域52Aに検査画像19が表示される。画像表示領域52Aの下方には、領域指定ボタン52C、クリアボタン52D、類似症例検索ボタン52Eが設けられている。
 領域指定ボタン52Cは、検査画像19内において関心領域ROIを指定するための操作ボタンである。領域指定ボタン52Cが、マウスのポインタ56でクリック操作されると、検査画像19内の任意の領域を指定する領域指定操作が可能になる。この状態で、ポインタ56を操作して、例えば、スプラインによって、対象病変OLを含む領域の外周を指定する。スプラインは、指定された複数の制御点を通る滑らかな曲線であり、ポインタ56で制御点を指定することによって入力される。こうした操作により対象病変OLを含む領域が、関心領域ROIとして指定される。クリアボタン52Dは、指定した関心領域ROIを取り消すためのボタンである。
 関心領域ROIは、複数個指定することが可能である。図10の例では、画像IDが「O901-01」から「O901-03」の3枚の検査画像19にそれぞれ1つずつ、No1~No3の関心領域ROIが指定されている。検査ID「901」の検査データ21において、他に関心領域ROIの指定が無ければ、1件の検査データ21に合計で3つの関心領域ROIが指定されたことになる。また、図11の例では、画像IDが「O906-01」の検査画像19内に2つの関心領域ROI(No1、No2)が指定され、画像IDが「O906-02」、「O906-03」の検査画像19内には、1つずつ関心領域ROI(No3、No4)が指定されている。No3の関心領域ROIには、2つの対象病変OL(No3、No4)が含まれている。このように、複数個の対象病変OLを含む領域を、1つの関心領域ROIとして指定してもよい。検査ID「906」の検査データ21において、他に関心領域ROIの指定が無ければ、1件の検査データ21に合計で4つの関心領域ROIが指定されたことになる。このように指定された各関心領域ROIには、異なる対象病変OLがそれぞれ1つ以上含まれる。
 図12に示すように、類似症例検索サーバ17には、AP50として類似症例検索サーバプログラムがインストールされており、プログラムを実行すると、類似症例検索サーバ17のCPU41Bは、要求受付部61、特徴量算出部62、個別類似度算出部65、総合類似度算出部66、類似症例検索部67、出力制御部69として機能する。
 要求受付部61は、診療科端末11から送信された類似症例検索要求を受け付けて、受け付けた検査画像19及び関心領域ROIの領域情報を、類似症例検索サーバ17のストレージデバイス43などに記憶する。特徴量算出部62は、受け付けた検査画像19及び領域情報に基づいて、関心領域ROIの特徴量を算出する。ここで、特徴量算出部62は、特徴量取得部として機能する。
 図13に示すように、関心領域ROIが複数有る場合には、No1の関心領域ROIの特徴量は、「RAC1」、No2の関心領域ROIの特徴量は、「RAC2」、No3の関心領域ROIは、「RAC3」というように、関心領域ROIの特徴量RACは、関心領域ROI毎に算出される。特徴量RACは、典型的な病変の画像パターンとして予め設定されている複数種類の病変パターンに対応する複数次元の数値(後述する判別器出力値)で構成される、特徴ベクトルである。典型的な病変パターンは、図6で説明したとおり、例えば、A~Hの8種類に分類される。
 図14に示すように、特徴量算出部62は、典型的な8種類の病変パターンに対応する判別器62A~62Hを有している。各判別器62A~62Hは、関心領域ROIの画像パターンに基づいて、典型的な各病変パターンに対応する数値を出力する。各判別器62A~62Hが出力する各数値は、特徴ベクトルを構成する複数次元の数値であり、ここで、各数値を判別器出力値と呼ぶ。本例では、各判別器62A~62Hに対応して判別器出力値は8種類有り、特徴ベクトルは8次元で構成される。なお、上述のとおり、典型的な病変パターンの種類は、8種類に限定されず、病変パターンの種類数に応じて、判別器の種類、特徴ベクトルの次元数も適宜決定される。
 判別器出力値は、典型的な病変パターンらしさを表すもので、関心領域ROI内に典型的な病変パターンが存在する度合いを表す値である。そのため、判別器出力値が大きいほど、関心領域ROIに典型的な病変パターンが存在する度合いが高いことを示し、判別器出力値が小さいほど、存在する度合いが低いことを示す。より詳細には、「+」の判別器出力値は、典型的な病変パターンが関心領域ROI内に存在することを示し、「-」の判別器出力値は、関心領域ROI内に存在しないことを示す。そして、「+」の出力値を示し、かつ、判別器出力値が大きいほど存在する度合いが高いことを示す。
 図14の例では、「B:空洞」の病変パターンの判別器62Bと「G:高吸収域」の判別器62Gが「+」の出力値を示し、「B:空洞」の出力値が最も大きいため、関心領域ROIは、「B:空洞」と「G:高吸収域」の病変パターンを含み、8種類の病変パターンの中では、「B:空洞」が支配的な画像パターンを持つものであることが分かる。
 なお、典型的な病変パターンに対する各判別器は、例えば、「資料名:Computer Vision and Image Understanding 88巻119頁~151頁,発行年2002年12月 Using Human Perceptual Categories for Content-Based Retrieval from a Medical Image Database  著者 Chi-Ren Shyu, Christina Pavlopoulou Avinash C.kak,and Cala E.Brodley」などに記載されているような周知の特徴量を用いて、「Ada-boost」などの機械学習アルゴリズムにより作成することができる。
 特徴量算出部62は、類似症例検索要求に添付された検査データ21において指定された複数の関心領域ROIのすべてについて、関心領域ROI毎に特徴量RACを算出する。
 図15に示すように、症例DB23内の特徴量DB23Bに格納される、各症例病変CLの特徴量CACについても、上記8種類の病変パターンに対応する特徴ベクトルで構成される。こうした特徴量CACについても、図14で説明した特徴量算出部62と同様の構成により算出され、各判別器62A~62Hの8種類の出力値から構成される8次元の特徴ベクトルで構成される。なお、各症例病変CLには、病変の種類が関連付けされて格納されている。
 図16に示すように、個別類似度算出部65は、関心領域ROIの特徴量RACと、症例病変CLの特徴量CACとの比較を行い、個別類似度ISMを算出する。具体的には、特徴量RACと特徴量CACに含まれる8次元の特徴ベクトルを比較して、個別類似度ISMを算出する。個別類似度ISMの値は、例えば、最小二乗距離または相関で算出する。前者の場合には、値が小さいほど、関心領域ROIと症例病変CLとの類似度が高く、後者の場合には、値が大きいほど、関心領域ROIと症例病変CLとの類似度が高い。
 図17に示すように、個別類似度算出部65は、1件の検査データ21に含まれる複数の関心領域ROIと、1件の症例データ24に含まれる症例病変CLとを1対1で対応させて、それぞれの特徴量RACと、特徴量CACの比較を行い、個別類似度ISMを算出する。個別類似度ISMは、関心領域ROI毎の類似度であるため、関心領域別個別類似度であり、症例病変CLと1対1で対応させて算出されるため、症例病変CL別の個別類似度でもある。個別類似度算出部65は、1つの関心領域ROIに対して、症例病変CLの個数だけ個別類似度ISMを算出する。個別類似度ISMは、すべての関心領域ROIについて算出されるので、個別類似度算出部65は、関心領域ROIの個数と、症例病変CLの個数を乗じた数だけ個別類似度ISMを算出する。
 検査IDが「O901」の検査データ21においては、No1~No3の3個の関心領域ROIが指定されており、症例IDが「C101」の症例データ24には、No1~No3の3個の症例病変CLが登録されている。「C101」の症例においては、No1~No3の3個の症例病変CLの種類は、それぞれ「B:空洞」、「F:点状影」、「E:すりガラス状陰影」で、すべて病変の種類が異なるので、3種3個の症例病変CLが登録されている。
 本例の個別類似度算出部65は、後述する総合類似度算出部66と異なり、症例病変CLの種類に着目せず、症例病変CLの個数にのみ着目して個別類似度ISMを算出する。このため、図18に示すように、「O901」の検査データ21と、「C101」の症例データ24との間では、3×3で合計9個の個別類似度ISMが算出される。
 各個別類似度ISMに付された括弧内の識別符号は、症例IDに、関心領域ROIと症例病変CLのそれぞれのシリアルナンバーを加えたものである。「C101-11」の場合には、No1の関心領域ROIと、症例IDが「C101」の症例データ24に登録された、No1の症例病変CLとの個別類似度ISMであることを示す。同様に、「C101-12」の場合には、No1の関心領域ROIと、症例IDが「C101」のNo2の症例病変CLの個別類似度ISMであることを示す。
 図19に示す、症例IDが「C102」の症例には、No1~No5の5個の症例病変CLが登録されている。「C102」の症例において、No1~No5の5個の症例病変CLの種類は、「B:空洞」、「F:点状影」、「E:すりガラス状陰影」の3種である。「B:空洞」については、同種の症例病変CL(No1~No3)が3個登録されており、No4、No5の症例病変CLについては、「F:点状影」と「E:すりガラス状陰影」がそれぞれ1個ずつ登録されているため、結局、「C102」の症例には、3種5個の症例病変CLが登録されている。
 個別類似度算出部65は、3個の関心領域ROIのそれぞれに対して、「C102」の5個の症例病変CLを、種類に拘わらず対応させるので、3×5で合計15個の個別類似度ISMを算出する。
 図20に示す症例IDが「C103」の症例は、2種4個の症例病変CLが登録されている例である。「C103」のNo1~No4の4個の症例病変CLは、種類が、「B:空洞」、「F:点状」の2種有り、「B:空洞」については、同種の症例病変CL(No1~No3)が3個で、「F:点状影」については、症例病変CL(No4)が1個である。個別類似度算出部65は、3個の関心領域ROIのそれぞれに対して、「C103」の4個の症例病変CLを、種類に拘わらず対応させて、3×4で合計12個の個別類似度ISMを算出する。
 図21の症例IDが「C104」の症例は、1種4個の症例病変CLが登録されている例である。「C104」のNo1~No4の4個の症例病変CLは、種類がすべて「B:空洞」で同種である。個別類似度算出部65は、3個の関心領域ROIのそれぞれに対して、「C103」の4個の症例病変CLを、種類に拘わらず対応させて、3×4で合計12個の個別類似度ISMを算出する。
 図22に示す症例IDが「C106」の症例は、図19に示す「C101」と同じ3種5個の症例病変CLが登録されている例であるが、内訳が異なる。すなわち、「C106」のNo1~No5の5個の症例病変CLの種類は、「B:空洞」、「F:点状影」、「E:すりガラス状陰影」の3種であり、この点については「C101」と同様である。しかし、「C106」において、「B:空洞」については、同種の症例病変CL(No1、No2)が2個登録されており、「F:点状影」についても、同種の症例病変CL(No3、No4)が2個登録されている。「E:すりガラス状陰影」については1個(No5)登録されている。個別類似度算出部65は、3個の関心領域ROIのそれぞれに対して、「C106」の5個の症例病変CLを、種類に拘わらず対応させて、3×5で合計15個の個別類似度ISMを算出する。
 図23に示すように、個別類似度算出部65は、例えば、類似症例検索サーバ17のメモリ42Bやストレージデバイス43B内に、個別類似度テーブル(以下、ISMテーブルという)71を作成して、ISMテーブル71に算出した個別類似度ISMを記録する。ISMテーブル71は、関心領域ROI毎に作成される。図23の例は、No1の関心領域ROIに対するISMテーブル71を示す。ISMテーブル71は、個別類似度ISM毎に、症例ID、病変No、病変の種類、病変画像を関連付けて記憶したテーブルである。病変画像は、症例病変CLの画像データである。すなわち、ISMテーブル71において、1件のレコードは、症例ID、病変No、病変の種類、病変画像、個別類似度ISMの5項目のデータで構成される。
 個別類似度算出部65は、最初に、ISMテーブル71に対して、各個別類似度ISMを算出した順番で記録する。各個別類似度ISMは、例えば「C101」、「C102」、「C103」というように、症例IDの番号が若い順から記録される。個別類似度ISMの値は、関心領域ROIの特徴量RACと症例病変CLの特徴量CACの相関で算出されるため、数値が大きいほど、類似度が高いことを示す。
 図24に示すように、個別類似度算出部65は、関心領域ROI毎にISMテーブル71を作成する。No1~No3の3個の関心領域ROIがある場合には、3つのISMテーブル71が作成される。この段階では、図23に示すとおり、ISMテーブル71は、症例IDの順番で各レコードが配列されている。個別類似度算出部65は、ISMテーブル71の作成が終了した場合には、ISMテーブル71を、総合類似度算出部66に引き渡す。
 図25に示すように、総合類似度算出部66は、関心領域ROI毎に作成された複数のISMテーブル71に基づいて、総合類似度TSMテーブル72(以下、TSMテーブル72という)を作成する。具体的には、図26に示すように、複数のISMテーブル71から、同一症例内の症例病変CLとの対応により算出された個別類似度ISMを1つずつ読み出し、各ISMテーブル71から読み出した複数の個別類似度ISMに基づいて、総合類似度TSMを算出する。より具体的には、総合類似度算出部66は、各ISMテーブル71から、1件ずつ読み出した個別類似度ISMを総合類似度TSMの要素として、関心領域ROIの数と症例病変CLの数に応じた数の個別類似度ISMの組み合わせを作り、組み合わせ毎に総合類似度TSMを算出する。
 個別類似度ISMは、複数の関心領域ROIと複数の症例病変CLとを1対1で対応させて算出したものであり、各関心領域ROIと各症例病変CLのどれとどれを対応させるかによって、個別類似度ISMの値は変化するので、個別類似度ISMの組み合わせパターンによって総合類似度TSMが変化する。
 個別類似度ISMの組み合わせパターンの作り方については、後述するが、本例のNo1~No3の3個の関心領域ROIと、症例IDが「C101」の3種3個の症例病変CLとを対応させる場合には、3×3=9つの個別類似度ISMが算出される。9つの個別類似度ISMは、関心領域ROI毎のISMテーブル71に3つずつ記録される。そして、総合類似度算出部66は、「C101」の9つの個別類似度ISMに関して、各ISMテーブル71が1つずつ個別類似度ISMを読み出して、個別類似度ISMの組み合わせのパターンを6通り作り、「C101-1」~「C101-6」の6通りの組み合わせに対応する6つの総合類似度TSMを算出する。各総合類似度TSMの識別符号は、症例IDに、総合類似度TSMの1~6までのシリアルナンバーを付したものである。各総合類似度TSMは、本例では3つの個別類似度ISMを足し合わせた総和である。
 各個別類似度ISMは、関心領域ROIと症例病変CLを1対1で対応させた場合の類似度であるので、総合類似度TSMが高いほど、3つの関心領域ROIと、3つの症例病変CLのそれぞれの個別類似度ISMの平均値が高いことを意味する。本例では、個別類似度ISMは、相関値で表しているため、数値が大きいほど類似度が高い。そのため、総合類似度TSMも、数値が大きいほど類似度が高いことを示す。
 本例では、6つの総合類似度TSMの中では、「C101-2」の総合類似度TSMが「2.04」で最も高い。一方、本例では「C101-3」の総合類似度TSMが「1.32」で最も低い。また、「C101-5」の総合類似度TSMには、個別類似度ISMの中に、「0.91」という最高値を示す個別類似度ISM(C101-13)を含んでいるが、最高値を含まない、「C101-2」の総合類似度TSMの方が、平均的に個別類似度ISMが高いため、数値が高くなっている。
 6つの総合類似度TSMの中で類似度が最高の総合類似度TSM(C101-2)は、No1の関心領域ROIとNo1の症例病変CLとの個別類似度ISM(C101-11)と、No2の関心領域ROIとNo3の症例病変CLとの個別類似度ISM(C101-23)と、No3の関心領域ROIとNo2の症例病変CLとの個別類似度ISM(C101-32)の総和である。そのため、「C101」の症例については、No1~No3の症例病変CLを、それぞれ、No1、No3、No2の関心領域ROIに対応させた場合が、「O901」の検査データ21に対して最も類似度が高いと評価できる。
 また、総合類似度算出部66は、総合類似度TSMを算出する個別類似度ISMの組み合わせについて、完全異種の組み合わせを作る。ここで、完全異種の組み合わせとは、同一症例内の複数の症例病変CLの中から、種類が異なる症例病変CLを種類毎に一つずつ抽出して、それらを組み合わせて構成される組み合わせを、更に関心領域ROIの一つ一つに対応させていった場合の全ての組み合わせをいう。そして、総合類似度算出部66は、完全異種の組み合わせを総合類似度TSMの算出対象と判定し、完全異種の組み合わせ以外の組み合わせについては、算出対象外と判定する。すなわち、総合類似度算出部66は、完全異種の組み合わせに対してのみ、個別類似度ISMを用いて総合類似度TSMを算出する。
 完全異種の組み合わせ以外の組み合わせには、完全同種の組み合わせと一部同種の組み合わせが含まれる。完全同種の組み合わせは、同一症例内の複数の症例病変CLの中から、種類が同じ複数の症例病変CLを一つずつ抽出して、抽出した複数の同種の症例病変CLの組み合わせを、更に関心領域ROIの一つ一つに対応させていった場合の全ての組み合わせをいう。一部同種の組み合わせは、種類が同じ2つ以上の症例病変CLと種類が異なる1つ以上の症例病変CLを少なくとも含む症例の中から、種類が同じ複数の症例病変CLと、種類が異なる1以上の症例病変CLとを一つずつ抽出して、抽出した、同種と異種が混在した複数の症例病変CLの組み合わせを、更に関心領域ROIの一つ一つに対応させていった場合の全ての組み合わせをいう。総合類似度算出部66は、各ISMテーブル71から1件ずつ読み出した個別類似度ISMに基づいて、関心領域ROIの個数と、1件の症例内の症例病変CLの個数及び種類数とに応じた数の完全異種の組み合わせを作る。
 図27~図34を用いて、個別類似度ISMの完全異種の組み合わせについて具体的に説明する。図27において、各個別類似度ISMに対して付された、「B」、「E」、「F」という符号は、各個別類似度ISMに対応する症例病変CLの種類を表す符号であり、各個別類似度ISMが、どの種類の症例病変CLを関心領域ROIに対応させて算出したものかを表すものである。
 「C101」の症例は、No1~No3の3個の症例病変CLの種類が、「B(空洞)」、「F(点状影)」、「E(すりガラス影)」ですべて異なっている。そのため、「B、E、F」の3種3個の症例病変CLを、3個の関心領域ROIに対して、「B、F、E」、「B、E、F」、「F、B、E」、「F、E、B」・・・・というように対応させて算出される個別類似度ISMの完全異種の組み合わせのパターンの数は、3個の関心領域ROIと3個の症例病変CLの順列の数と同じになる。すなわち、33=3×2×1=6通りとなる。
 図19で示した「C102」の症例のように、「B」が3個で、「E」、「F」がそれぞれ1個の3種5個の症例病変CLを有している場合には、No1~No3の関心領域ROIに対応する各ISMテーブル71には、図28に示すように、5個ずつ合計15個の個別類似度ISMが記録される。「C102」の症例のように、同種の症例病変CLが複数個存在する場合には、総合類似度算出部66は、B1~B3の同種の3個の症例病変CLを区別して、区別した同種の症例病変CL毎に完全異種の組み合わせを作る。
 「C102」の症例のように、B1~B3の同種の症例病変CLを複数個有している場合には、「C102」の症例における完全異種の組み合わせの数は、症例病変CLの種類が異なる「C101」の症例のように順列の数とはならず、図29に示すように18通りとなる。総合類似度TSMはTSM(C102-1)~(C102-18)の18個算出される。
 「C102」において、「C101」のように順列の数を求めると、「C102」の5個の症例病変CLから、関心領域ROIの数である3個を取り出したときの並べ方となるから、53=5×4×3=60通りとなる。完全異種の組み合わせの数は、この順列の数から、図30に示す「B1、B2、B3」、「B1、B3、B2」という完全同種の組み合わせや、「B1、B2、F」、「B1、B2、F」という一部同種の組み合わせを除外した数である。総合類似度算出部66は、図29に示す完全異種の組み合わせに対しては、総合類似度TSMの算出対象外と判定し、図30に示す完全異種の組み合わせ以外の組み合わせに対しては、算出対象外と判定する。
 また、図31に示す「C106」の症例のように、「C102」と同じ3種5個の症例病変CLを含む場合でも、3種5個の症例病変CLの内訳が異なると、完全異種の組み合わせの数も変化する。「C106」において、3種5個の症例病変CLの内訳は、「B」が2個(B1、B2)、「F」が2個(F1、F2)、「E」が1個である。そのため、「B1、F1、E」、「B1、F2、E」といった完全異種の組み合わせの数は24通りとなり、総合類似度算出部66は、図31に示す完全異種の24通りの組み合わせに対しては総合類似度TSMの算出対象と判定して、C106-1~C106-24の24個の総合類似度TSMを算出する。一方、図32に示すように、「B1、B2、F1」、「B1、F1、F2」という一部同種の組み合わせのように、完全異種の組み合わせ以外の組み合わせに対しては算出対象としない。
 このように、完全異種の組み合わせのみを総合類似度TSMの算出対象とすることにより、総合類似度TSMの算出対象となる組み合わせが絞り込まれるため、総合類似度算出部66の処理時間や処理負荷が軽減されて、類似症例の検索時間を短縮することができる。また、上述のとおり、本発明において、複数の関心領域ROIに着目して類似症例を検索しなければならない主要な理由は、空洞、点状影、すりガラス状陰影という3種の対象病変の存在が疾患特定の根拠となる結核などの非がん疾患の診断に利用するためである。そのため、異種の対象病変OLを含む複数の関心領域ROIに類似する症例として、異種の症例病変CLを含む症例を類似症例として検索することが必要になる。
 非がん疾患の診断に利用するという目的から考えれば、複数の関心領域ROIは、それぞれ異種の対象病変OLを含むように指定される場合が多く、そうだとすると、そもそも同種の症例病変CLを含む症例を抽出する必要性は低いので、総合類似度TSMの算出対象を完全異種の組み合わせに限定しても、実際の診断において影響は少ない。
 このような理由から、さらに、総合類似度算出部66は、図33に示す、2種4個の症例病変CLを含む「C103」の症例や、図34に示す、1種4個の症例病変CLを含む「C104」の症例については、症例病変CLの種類数が、関心領域ROIの数に満たないため、総合類似度TSMの算出対象としない。「C104」のように症例病変CLが1種の場合には、異種の症例病変CLを検索するという、検索目的から考えれば、検索対象にする必要は低い。
 また、「C103」のように症例病変CLの種類数が、関心領域ROIの数に満たない症例を対象外とする理由も、検索対象とする必要性は低い。つまり、非がん疾患の診断において、医師が関心領域ROIを3個指定するのは、3種の症例病変CLを含む症例の検索を望むからであり、関心領域ROIの数未満の種類数の症例病変CLが抽出されても、診断との関係では必要性は低いと考えられるからである。なお、本例において、2種の症例病変CLを含む「C103」の症例が算出対象外となるのは、本例が、3個の関心領域ROIを前提としているためであり、関心領域ROIが2個の場合には、当然ながら、「C103」についても総合類似度TSMの算出対象となる。
 図35に示すように、以上の内容をまとめると、総合類似度算出部66は、症例IDが「C101」、「C102」、「C105」、「C106」の症例のように、症例病変CLの種類数が、関心領域ROIの数以上ある症例についてのみ、総合類似度TSMを算出し、類似症例の検索対象とする。総合類似度TSMは、完全異種の組み合わせに対してのみ算出される。反対に、「C103」、「C104」のように、症例病変CLの登録数が関心領域ROIの数に満たない症例(本例では症例病変CLが3個未満)については、総合類似度TSMを算出せず、検索の対象外とする。
 図36に示すように、総合類似度算出部66は、検索対象となる複数の症例について算出した総合類似度TSMを、TSMテーブル72に記録する。TSMテーブル72において、1件のレコードは、症例ID、総合類似度TSM、個別類似度ISMの組み合わせパターンの3項目のデータで構成される。組み合わせパターンは、各総合類似度TSMの算出元となった個別類似度ISMの組み合わせである。総合類似度算出部66は、作成したTSMテーブル72を、類似症例検索部67に引き渡す。図12において、類似症例検索部67には、代表値判定部67Bが設けられている。
 図37に示すように、代表値判定部67Bは、TSMテーブル72において、複数の総合類似度TSMの中から、症例毎に代表値を判定する。「C101」の症例の場合には、6個の総合類似度TSMの中から、類似度が最も高い(相関値が最大)、総合類似度TSM(C101-2)が代表値と判定される。同様に、「C102」の症例の場合には、18個の総合類似度TSMの中から、類似度が最も高い(相関値が最大)、総合類似度TSM(C102-3)が代表値と判定される。
 図38に示すように、代表値判定部67Bは、代表値判定を、TSMテーブル72内の全症例について実施する。これにより、TSMテーブル72には、症例毎に1つずつ判定された代表の総合類似度TSMのみが抽出される。
 図39に示すように、類似症例検索部67は、代表の抽出が行われたTSMテーブル72において、総合類似度TSMが高い順に各レコードをソートする。TSMテーブル72において、総合類似度TSMは、1件の症例から1つだけ抽出されているので、各レコードをソートすることは、症例を総合類似度TSMが高い順にソートすることを意味する。これにより、各症例に対して順位付けが行われて、TSMテーブル72の上位に、類似度が高い類似症例が抽出される。
 類似症例検索部67には、リスト作成部67A(図12参照)が設けられている。リスト作成部67Aは、TSMテーブル72に基づいて、図40に示す類似症例リスト74を作成する。類似症例リスト74は、検索結果表示画面76内に表示される。類似症例リスト74は、複数件の類似症例に関する情報をリスト化したものである。検索結果表示画面76は、類似症例検索サーバ17が、類似症例検索要求の要求元である診療科端末11に対して検索結果を配信する画面の例である。
 リスト作成部67Aは、TSMテーブル72から、総合類似度TSMを抽出して、総合類似度TSMが高い順に、類似症例を配列した類似症例リスト74を作成する。類似症例リスト74の表示項目には、各総合類似度TSMに関して、総合類似度TSMの値、総合類似度TSMに基づく順位、症例ID、及び総合類似度TSMに関する内訳情報を含む。本例では、内訳情報は、総合類似度TSMの算出要素となる各個別類似度ISMを算出するための関心領域ROIと症例病変CLとの対応関係である。内訳情報には、症例病変CLの種類(病変種類)も含まれている。
 類似症例リスト74において、症例IDが「C106」の総合類似度TSMは、図39に示す「C106-1」の総合類似度TSMに対応する。「C106-1」の総合類似度TSMの算出要素となる個別類似度ISMの組み合わせパターンは、個別類似度ISM(C106-11)、(C106-22)、(C106-33)である。類似症例リスト74に表示される内訳情報は、これらの各個別類似度ISMを算出するための関心領域ROIと症例病変CLの対応関係である。個別類似度ISM(C106-11)は、No1の関心領域ROIと、「C106」のNo1の症例病変CLであり、個別類似度ISM(C106-22)は、No2の関心領域ROIと、「C106」のNo2の症例病変CLであり、個別類似度ISM(C106-33)は、No3の関心領域ROIと、「C106」のNo3の症例病変CLである。対応関係には、種類も表示されるため、個別類似度ISMの組み合わせパターンが、どのような種類の組み合わせかを確認することができる。
 さらに、類似症例リスト74の表示項目には、No1、No2、No3の症例病変CLの病変画像が含まれる。病変画像は、例えば、ISMテーブル71から読み出される。さらに、類似症例リスト74の上方には、No1~No3の関心領域ROIを含む検査画像19が表示される。
 類似症例リスト74は、例えば、上位6位までの症例が表示される。もちろん、画面スクロールなどの操作により6位以下の順位を表示できるようにしてもよい。また、上位10位まで表示するというように、同時に表示できる表示件数を変更できるようにしてもよい。
 出力制御部69(図12参照)は、このように作成される検索結果表示画面76に関して、例えば、XML(Extensible Markup Language)などのマークアップ言語によってWEB配信用のXMLデータを作成して、これを検索結果として、要求元の診療科端末11に配信する制御を行う。XMLデータを受信した診療科端末11は、WEBブラウザがXMLデータに基づいて検索結果表示画面76を再現して、表示部48Aに表示する。これにより、類似症例リスト74を含む検索結果表示画面76が医師の閲覧に供される。
 以下、上記構成の作用について、図41及び図42を参照しながら説明する。図41に示すように、診療科10の医師は、診療科端末11で検査画像DBサーバ15にアクセスして、検査科12に依頼した検査の検査データ21を取得する(S1100)。診療科端末11は、検査データ21を表示部48Aに表示する(S1200)。取得した検査データ21に含まれる検査画像19は、図10に示す検査画像表示画面52に表示される。検査画像表示画面52を通じて、医師により検査画像19内に関心領域ROIの指定が行われる。診療科端末11は、医師の指定操作により、複数の関心領域ROIの指定を受け付ける(S1300)。
 関心領域ROIの指定において、結核やびまん性汎細気支炎などの非がん疾患の診断の際には、病変の種類が異なる異種の対象病変OLをそれぞれに含むように複数の関心領域ROIが指定される。関心領域ROIの指定が終了すると、類似症例検索ボタン52Eが操作される。これにより、診療科端末11は、検索指示を受け付ける(S1400)。検索指示を受け付けると、検索要求発行部54が、検査画像19と領域情報とが添付された類似症例検索要求を発行して、これが類似症例検索サーバ17に送信される(S1500)。
 類似症例検索サーバ17は、類似症例検索要求を受信すると、要求受付部61がこれを受け付ける(S2100)。そして、特徴量算出部62は、検査画像19と関心領域ROIの領域情報とに基づいて、関心領域ROI毎の特徴量を算出する(S2200)。この後、類似症例検索が実行される(S2300)。
 図42に示すように、類似症例検索では、まず、個別類似度算出部65が、症例DBサーバ16から、1件の症例データ24を読み出す(S2310)。そして、個別類似度算出部65は、検査データ21内の複数の関心領域ROIと、1件の症例データ24に含まれる症例病変CLとを1対1で対応させて、個別類似度ISMを算出する(S2320)。症例病変CLが複数有る場合には、症例病変CL毎に個別類似度ISMを算出する。個別類似度算出部65は、算出した個別類似度ISMをISMテーブル71に記録して、関心領域ROI毎にISMテーブル71を作成する(S2330)。こうした処理を、1件の症例データ24について行った後、次の症例データ24に対しても同様に行う。そして、複数件の症例データ24、例えば、症例DB23内のすべての症例データ24について個別類似度ISMの算出及びISMテーブル71の作成が終了するまで(S2340でN)、同様の処理を繰り返す。
 総合類似度算出部66は、関心領域ROI毎に作成した複数のISMテーブル71に基づいて、TSMテーブル72を作成する(S2350)。TSMテーブル72の作成においては、総合類似度算出部66は、まず、症例毎の症例病変CLの種類数に基づいて検索対象とする症例を判定する(S2350)。この判定により、図34に示す症例IDが「C104」のように、症例病変CLが1種類しか無い症例や、図33に示す「C105」のように、症例病変CLの種類数が関心領域ROIの数未満の症例が検索対象外と判定される。
 次に、総合類似度算出部66は、検索対象の症例について、症例毎に個別類似度ISMの完全異種の組み合わせを作成して(S2352)、作成した完全異種の組み合わせに対してのみ総合類似度TSMを算出する(S2353)。総合類似度算出部66は、「C102」の症例については、図29に示す、18通りの完全異種の組み合わせを作り、18通りの組み合わせを算出対象と判定する。そして、判定した算出対象の組み合わせ毎に18個の総合類似度TSM(C102-1~C102-18)を算出する。一方、総合類似度算出部66は、図30に示す、完全同種の組み合わせや一部同種の組み合わせについては算出対象外と判定して、総合類似度TSMを算出しない。
 同様に、総合類似度算出部66は、「C106」の症例については、図31に示すように24通りの完全異種の組み合わせを作り、24通りの組み合わせに対応する24個の総合類似度TSM(C106-1~C106-24)を算出する。一方、総合類似度算出部66は、図32に示す、完全同種の組み合わせや一部同種の組み合わせについては算出対象外と判定して、総合類似度TSMを算出しない。そして、総合類似度算出部66は、算出した総合類似度TSMを症例毎にTSMテーブル72に記録する(S2354)。
 類似症例検索部67は、作成されたTSMテーブル72に基づいて、類似症例リスト74を作成する(S2360)。リスト作成においては、まず、代表値判定部67Bは、図38に示すように、TSMテーブル72内の症例毎の複数の総合類似度TSMの中から、代表値を判定して(S2361)、各症例の代表の総合類似度TSMのみを抽出したTSMテーブル72を作成する。この後、類似症例検索部67は、TSMテーブル72において、総合類似度TSMが高い順に症例をソートする(S2362)。これにより、TSMテーブル72の上位に、類似度が高い類似症例が抽出される。
 リスト作成部67Aは、TSMテーブル72に基づいて、所定順位までの症例を抽出して、総合類似度TSMが高い順に、類似症例を配列した類似症例リスト74を作成する(S2363)。
 図41において、出力制御部69は、リスト作成部67Aが検索結果として作成した類似症例リスト74を含む検索結果表示画面76を、配信用のXMLデータに変換して、診療科端末11に配信する(S2400)。診療科端末11は、類似症例リスト74を含むXMLデータを受信して(S1600)、XMLデータに基づいて検索結果表示画面76(図30)を再現して、表示部48Aに表示する(S1700)。
 検査データ21に複数の対象病変OLが含まれる場合において、検査データ21に類似する類似症例を検索する場合には、各対象病変OLを含む複数の関心領域ROIと複数の症例病変CLのそれぞれの特徴量に着目することに加えて、個々の特徴量を総合的に考慮した検索を行うことが望まれる。例えば、ある症例において、1つの症例病変CLが1つの関心領域ROIに対して高い類似度を示していたとしても、他の症例病変CLが他の関心領域ROIの類似度が非常に低い場合には、その症例は、少なくとも複数の関心領域ROIに着目すべき場合には、類似症例として適切ではない。
 本発明は、関心領域ROI毎に症例病変CLとの個別類似度ISMを算出し、算出した個別類似度ISMに基づいて総合類似度TSMを算出して類似症例検索を行っている。総合類似度TSMは、複数の個別類似度ISMが平均的に高い症例を、類似度が高い症例と評価する指標である。総合類似度TSMに基づいて類似症例を検索することにより、複数の対象病変OLを含む検査データ21に対して、類似度の高い適切な類似症例を検索することができる。
 従来技術は、1つの関心領域ROIの特徴量にのみ着目した類似症例検索しか行われていないため、複数の関心領域ROIがある場合の類似症例検索においては適切な類似症例を検索することができない。対して、本発明は、総合類似度TSMに基づいて類似症例を検索するから、複数の関心領域ROIがある場合の類似症例検索において、従来技術と比較して有用性の高い技術を提供することができる。
 空洞影、点状影、すりガラス影という3種の対象病変OLに着目する必要がある結核や、気管支異常影と点状影の2種の対象病変OLに着目する必要があるびまん性汎細気支炎など、非がん疾患では、複数の対象病変OLが現れていることが疾患の特定の根拠となる場合がある。このように複数の関心領域ROIの特徴量に注意を払わなければならない非がん疾患の診断において、本発明は有用である。
 また、本発明は、総合類似度TSMについて、完全異種の組み合わせのみを総合類似度TSMの算出対象とすることにより、総合類似度TSMの算出対象となる組み合わせを絞り込んでいる。そのため、総合類似度算出部66の処理時間や処理負荷が軽減されて、類似症例の検索時間を短縮することができる。図29及び図30に示した、「C102」の3種5個の症例病変CLと、3個の関心領域ROIとを対応させる例で説明すると、完全異種の組み合わせ以外の組み合わせを含めた順列の数は、53=5×4×3=60通りとなるところ、完全異種の組み合わせは18通りとなる。そのため、「C102」に関する組み合わせの数だけで単純に計算すると、総合類似度TSMの算出処理時間を約1/3にすることができる。
 また、このように、総合類似度TSMの完全異種の組み合わせに限定しても、類似症例検索を非がん疾患の診断に利用する場合には、実際の診断における影響は少ない。というのも、類似症例検索を非がん疾患の診断に利用する場合には、複数の関心領域ROIは、それぞれ異種の対象病変OLを含むように指定されるため、同種の症例病変CLを含む症例を抽出する必要性は低い。そのため、総合類似度TSMの算出対象を完全異種の組み合わせに限定しても、必要な類似症例については抽出することができるからである。
 なお、本例では、総合類似度算出部66が総合類似度TSMの算出対象外と判定する症例についても、個別類似度ISMの算出対象としているが、総合類似度算出部66が算出対象外と判定する症例については、個別類似度ISMの算出対象外としてもよい。例えば、本例においては、図31に示す2種4個の症例病変CLを有する「C103」の症例や、1種4個の症例病変CLを有する「C104」の症例のように、症例病変CLの種類数が、関心領域ROIの数未満の症例については、図35に示すように、検索対象外と判定される。このような症例については、個別類似度算出部65が、個別類似度ISMの算出の際に検索対象外と判定して、個別類似度ISMも算出しない。個別類似度算出部65は、症例病変CLの種類数と関心領域ROIの数とに基づいて、検索対象とするか否かを判定する。こうすれば、無駄な処理時間が減るため、検索時間をさらに短縮できる。
 また、本例では、総合類似度TSMを、複数の個別類似度ISMを足し合わせた総和としたが、掛け合わせた値でもよい。
 なお、本例では、代表値判定部67Bにより、症例毎に算出される複数の総合類似度TSMの中から、症例毎に代表値を判定して、代表値のみに基づいて、類似症例検索を行っている。代表値を判定することにより、TSMテーブル72に記録される総合類似度TSMの件数が減るなど、検索処理において取り扱うデータが減るため、処理時間を短縮できるという効果がある。読影レポートに記載された、症例病変CLに関する医師の所見などの診断結果は症例ごとに存在するため、代表値判定により、検索結果を症例単位で提示することで、類似症例に基づく診断を適切にかつ効率よく行うこともできる。しかし、代表値を判定せずに、類似症例検索を行ってもよい。代表値を判定しない場合には、類似症例リスト74において、同じ症例の複数の総合類似度TSMが表示されることがあり得る。これら複数の総合類似度TSMは、算出要素となる個別類似度ISMの組み合わせパターンが異なるものであるため、1つの症例について視点を変えて参照することができる。
 また、類似症例リスト74の表示項目には、類似症例の症例IDと総合類似度TSMに加えて、総合類似度TSMに関する内訳情報が表示される。内訳情報によって、総合類似度TSMの算出要素となる各個別類似度ISMを算出するための関心領域ROIと症例病変CLの対応関係を確認することができる。この対応関係を表示することにより、総合類似度TSMの算出要素となる各個別類似度ISMは、複数の関心領域ROIと複数の症例病変CLのどれとどれを対応させて算出したものかを把握することができる。また、類似症例リスト74には、検査画像19や症例病変CLの病変画像も表示されるので、画像パターンについても比較参照しやすく、画像パターン同士の類似度を直感的に判断しやすい。
 総合類似度TSMに関する内訳情報としては、図43に示す類似症例リスト75のように、総合類似度TSMの算出要素となる個別類似度ISMの値を表示してもよい。総合類似度TSMの値に加えて、その内訳となる個別類似度ISMの値を表示すれば、各関心領域ROIと各症例病変CLの対応の中で、個別類似度ISMが高いものや低いものがどれかを確認できるため、便利である。例えば、医師が、複数の関心領域ROIのうちの1つを重視したい場合には、重視する関心領域ROIの個別類似度ISMの値を見て、適切な類似症例を探し出すといったことができる。
 また、総合類似度TSMは、複数の個別類似度ISMが平均的に高い症例を、類似度が高い症例と評価する指標である。そのため、類似症例リストでは、複数の個別類似度ISMの平均値が高い症例の順位が高く、平均値が低い症例については、仮に1つの個別類似度ISMが突出して高くても、順位は低くなる。これに対して、類似度に関する医師の主観的な評価は、平均値よりも、特定の症例病変CLと特定の関心領域ROIとに関して医師が抱いた印象に大きく影響される場合があるため、医師の主観的な評価と総合類似度TSMに基づく客観的な評価(順位付け)が乖離することもある。
 このような乖離が生じた場合でも、類似症例リスト75のように、総合類似度TSMに加えて、個別類似度ISMの値が表示されれば、個別類似度ISMを確認することにより、医師が自らの主観的な評価を検証することができる。また、個別類似度ISMの値が表示されれば、医師が個別類似度ISMの値を勘案することにより、類似症例リスト75による、総合類似度TSMに基づく客観的な評価に対して、医師が主観的な評価に基づいて修正を加えながら、個々の診断に適した適切な類似症例を探し出すといったこともできる。
 また、本発明では、異種の対象病変OLをそれぞれ含む複数の関心領域ROIが指定されることを想定して、総合類似度TSMの算出対象を、複数の個別類似度ISMの組み合わせに関して完全異種の組み合わせに限定している。そのため、仮に、同種の対象病変OLをそれぞれ含む複数の関心領域ROIが指定された場合(例えば「B:空洞」が2個)でも、同種の複数の関心領域ROIに対して、同種の症例病変CLを対応させて算出した個別類似度ISMの組み合わせは総合類似度TSMの算出対象外となるため、類似症例リストにおいて類似症例として抽出されることはない。
 類似症例として抽出される可能性があるのは、同種(例えば「B:空洞」が2個)の複数の関心領域ROIに対して、少なくとも1つが同種の症例病変CL(例えば「B:空洞」と「F:点状影」や、「B:空洞」と「E:すりガラス状陰影」)を対応させて算出した個別類似度ISMの組み合わせである。このような場合において、「B」と「F」というように、異種の関心領域ROIと症例病変CLを対応させた場合の個別類似度ISMの値は無視できる程度に低くなると考えるため、総合類似度TSMに影響するのは、結局、同種の関心領域ROIと症例病変CLを対応させた場合の個別類似度ISMの値である。類似症例リスト75のように、総合類似度TSMの値に加えて、個別類似度ISMの値を表示すれば、特定の関心領域ROIに対する個別類似度ISMが高い症例病変CLを探し出しやすい。
 総合類似度算出手順に関して、上記実施形態では、症例全件分の症例病変CLについての個別類似度ISMをすべて求めた後、症例毎に、個別類似度ISMの完全異種の組み合わせを作成して、総合類似度TSMを算出している。総合類似度算出手順としては、この他に次の方法がある。まず、症例全件分ではなく、1件の症例について、その症例に含まれる症例病変CLについての個別類似度ISMを算出する。次に、個別類似度ISMを算出した1件の症例に関して、算出した個別類似度ISMの完全異種の組み合わせを作成して、総合類似度TSMを算出する。そして、1件の症例について総合類似度TSMの算出が終了した後、次の症例について、個別類似度ISMの算出、完全異種の組み合わせの作成、総合類似度算出の各手順を実行する。これを症例全件分について繰り返す。
 どちらの方法にしても、1件の症例に関する総合類似度算出手順を考えれば、個別類似度ISMを算出した後に、個別類似度ISMについての完全異種の組み合わせを作成し、完全異種の組み合わせについてのみ総合類似度TSMを算出するという手順については共通である。本発明は、どちらの方法でもよく、当然ながら、本発明の類似症例検索装置、方法及びプログラムには、どちらの方法も包含される。
[第2実施形態]
 第1実施形態においては、症例病変CLの種類数が、関心領域ROIの数未満の症例については、検索の対象外としたが、図44~図46に示す第2実施形態のように、症例病変CLの種類数が、関心領域ROIの数未満であっても、種類数が複数の症例については、検索対象に含めてもよい。例えば、図44に示すように、第2実施形態においては、3個の関心領域ROIが指定されている場合において、第1実施形態(図35参照)では対象外とされた症例IDが「C103」の症例のように、種類数が2個の症例が検索対象に含まれる。
 なお、症例IDが「C104」の症例は、種類が1つであるので、第2実施形態においても検索対象外となる。なぜならば、総合類似度TSMは、複数の関心領域ROIに対応する複数種類の個別類似度ISMの和であるため、種類数が「1」のものは、総合類似度TSMを算出できないからである。
 類似症例として抽出される症例については、症例病変CLの種類数が関心領域ROIの数以上あることが好ましい。類似症例検索を非がん疾患の診断に利用する場合には、関心領域ROIの数は、対象病変OLの種類数に対応するものだからである。しかしながら、症例DB23において、症例病変CLの種類数が多い症例の蓄積数が少ない場合がある。この場合に、関心領域ROIの数に満たない種類数の症例を、すべて検索対象外とすると、検索対象が少な過ぎて、適切な類似症例が検索できないおそれもある。症例病変CLの種類数が関心領域ROIの数に満たない症例であっても、個々の症例病変CLについて見ると、診断に際して有用な価値を持つ場合もある。第2実施形態は、症例病変CLの種類数が複数有ることを前提に、種類数が関心領域ROIの数よりも少ない症例については検索対象に含めて、類似症例として抽出することを可能にする形態である。なお、特徴量の算出方法や個別類似度ISMの算出方法は、上記実施形態と同様であるので、説明を省略し、相違点を中心に説明する。
 図45に示すように、個別類似度算出部65は、No1~No3の関心領域ROIに2種の症例病変CLを対応させて、個別類似度ISMを算出する。図35に示すハッチングの欄は、対応する症例病変CLが無いため、個別類似度ISMが算出されない箇所を示す。総合類似度算出部66は、3個の関心領域ROIに対して、2種4個の症例病変CLを対応させて算出した個別類似度ISMの完全異種の組み合わせを作り、その組み合わせを総合類似度TSMの算出対象と判定する。「C103」の症例では、完全異種の組み合わせは、「B1、F」、「B2、F」などである。これらは、種類数と関心領域ROIの数とが異なるが、同種の症例病変CLとの対応により算出した同種の個別類似度ISMを含んでいないため、完全異種の組み合わせである。一方、「B1、B2」などの完全同種の組み合わせは対象外と判定される。図45に示すとおり、「C103」において、完全異種の組み合わせは、18通りとなり、総合類似度TSMは、「C103-1」~「C103-18」の18個となる。
 類似症例検索部67は、図46に示す類似症例リスト81を作成する。ここで、類似症例リスト81において、「C103」の症例の総合類似度TSM(1.67」)は、2種類の個別類似度ISMの合計であり、3種類の個別類似度ISMの合計である他の症例の総合類似度TSMの値と同列に比較できない。そのため、類似症例検索部67は、算出要素となる個別類似度ISMの数が異なる総合類似度TSMを比較できるように、正規化を行う。正規化された類似度である正規化値に基づいて、各症例の順位付けが行われる。正規化は、例えば、総合類似度TSMを、それぞれ個別類似度ISMの数で割る処理である。
 まず、症例IDが「C103」については、個別類似度ISMの数が2個であるため、総合類似度TSM(「1.67」)を2で割った値(「0.835」)が正規化値となる。一方、個別類似度ISMが3個である他の症例は、総合類似度TSMを3で割った値が正規化値となる。例えば、症例IDが「C106」の場合には、総合類似度TSMの値が「2.34」であるので、これを3で割った値「0.78」が正規化値となる。
 本例によれば、症例DB23内において、症例病変CLの種類数が関心領域ROIの数未満であっても、症例病変CLの種類数が複数有る症例については、検索対象に含めることができるため、症例DBへの症例の登録件数が少ない場合でも、それらを有効利用して類似症例検索を行うことができる。
 また、第2実施形態の変形例として、関心領域ROIの指定数の変更に応じて、類似症例検索サーバ17が再検索を行うようにしてもよい。例えば、1回目の類似症例検索要求は、関心領域ROIの指定数を2で発行する。類似症例検索サーバ17は、その指定に基づいて類似症例検索を行って、類似症例リスト74を検索結果として配信する。医師は、その検索結果を見て、新たに関心領域ROIを追加指定して、2回目の類似症例検索要求を発行する。類似症例検索サーバ17は、追加された指定を含む要求に基づいて、類似症例検索を行って、類似症例リスト74を検索結果として配信する。こうすれば、医師は、検索結果を見ながら、必要に応じて検索条件を変更できるので、適切な類似症例を探しやすい。関心領域ROIの指定数は、増やすだけでなく、減らすことが可能であってもよい。
 また、類似症例検索サーバ17は、再検索に備えて、ISMテーブル71やTSMテーブル72など、類似症例検索の中間処理や処理結果として作成されたデータを保存しておくとよい。再検索の際にこうしたデータを利用すれば、検索時間を短縮できる。
 上記第1実施形態及び第2実施形態において、総合類似度TSMを複数の個別類似度ISMの単なる合計値とした例で説明したが、算出要素となる個別類似度ISMの値に基づいて、重み付け処理を行ってもよい。重み付け処理の例としては、例えば、個別類似度ISMが閾値以上である場合に、正の重み付け係数を掛けて、総合類似度TSMを高くしたり、反対に、個別類似度ISMが閾値未満である場合に、負の重み付け係数を掛けて、総合類似度TSMを低くする例が考えられる。関心領域ROIに酷似した症例病変CLが1つでもある場合には、その症例が診断において価値を持つ場合もある。正の重み付けを行うことで、そうした症例を類似症例として探しやすくなる。
 また、上記各実施形態では、症例病変CLの種類について、予め種類判定を行い、判定した種類情報を症例DB23に格納した例で説明したが、関心領域ROIと同様に、検索の時点で特徴量CACに基づく種類判定を行ってもよい。もちろん、こうすると検索処理に時間が掛かるため、予め種類情報を格納しておくことが好ましい。
 [第3実施形態]
 上記各実施形態では、関心領域ROIに含まれる対象病変OLの種類を判定することなく症例病変CLと対応させて個別類似度ISMを算出して、類似症例検索を行っている。しかし、図47~図51に示す第3実施形態のように、関心領域ROIに含まれる対象病変OLと症例病変CLについて、病変の種類の判定を行って、同種の病変同士のみを対応させて個別類似度ISMを算出して、類似症例検索を行ってもよい。図6に示したように、病変パターンは典型的には病変の種類によって区別される。そのため、特徴量を算出した段階で、特徴量に基づいて病変の種類を判定することができる。第3実施形態は、このような病変の種類判定を利用する形態である。
 図47に示すように、第3実施形態では、類似症例検索サーバ17には、病変種類判定部86が設けられる。図48に示すように、病変種類判定部86は、特徴量算出部62が算出した関心領域ROIの特徴量RACに基づいて、関心領域ROIに含まれる対象病変OLの種類を判定する。病変種類判定部86は、例えば、各判別器62A~62Hが出力する判別器出力値のうち、最大の判別器出力値を示す判別器に対応する病変の種類を、関心領域ROIに含まれる対象病変OLの種類と判定する。本例の場合には、「B:空洞」の判別器62Bの判別器出力値が最大であるため、対象病変OLの種類は「B:空洞」と判定される。
 図49に示すように、個別類似度算出部65は、関心領域ROIと、各症例病変CLとの個別類似度ISMを算出する際に、病変の種類が同種のもの同士についてのみ個別類似度ISMを算出し、病変の種類が異なるものについては、個別類似度ISMを算出しない。本例においては、No1の関心領域ROIの種類は「B:空洞」であるので、「C101」の症例の中から、種類が「B:空洞」である、No3の症例病変CLとの個別類似度ISMのみ算出する。1件の症例の中に、関心領域ROIと同種の症例病変CLが複数個登録されている場合には、複数個の個別類似度ISMが算出される。また、関心領域ROIと同種の症例病変CLが1個も登録されていない場合には、その症例については、個別類似度ISMは算出されない。
 これにより、個別類似度算出部65の算出処理時間を減らすことができる。また、個別類似度ISMの数が減るため、図50に示すように、総合類似度TSMの算出対象となる完全異種の組み合わせの数も減る。3種5個の「C102」の症例の場合には、第1実施形態で説明したように種類を判定しないと、完全異種の組み合わせが18通り有るのに対して、本例のように種類判定をすると、完全異種の組み合わせは3通りとなる。総合類似度算出部66は、この組み合わせに対応して、C102-1~C102-3の3つの総合類似度TSMを算出する。
 このように、第3実施形態は、病変の種類を区別せずに個別類似度ISMを算出する第1実施形態と比較して、個別類似度ISM及び総合類似度TSMの算出時間が短くなる。また、ISMテーブル71やTSMテーブル72のサイズも小さくなるため、メモリの作業領域も少なくて済む。そのため、類似症例検索サーバ17のCPU41Bに掛かる負荷が軽減されるため、検索時間の短縮化を図ることができる。
 ただし、予め病変の種類を判定して同種のものについてのみ個別類似度ISMをする態様は、病変の種類の判定精度が低い場合には、類似症例として検索すべき症例病変CLが抜け落ちる、いわゆる検索漏れの懸念もある。特に、図11で示すように、複数の対象病変OLを1つの関心領域ROIとして指定する場合には、種類判定によって、複数の対象病変OLのいずれかに偏った判定が下されることになる。そのため、第5実施形態については、病変の種類の判定精度を見極めた上で実施することが好ましい。
[第4実施形態]
 図51及び図52に示す第4実施形態は、類似症例検索サーバ17ではなく、診療科端末11において関心領域ROIの特徴量を算出する形態である。第3実施形態のように診療科端末11で関心領域ROIの特徴量を算出してもよく、この場合には、類似症例検索サーバ17には特徴量算出部62を設けずに、図12に示す、個別類似度算出部65、総合類似度算出部66、類似症例検索部67など、特徴量算出部62以外の構成が設けられる。
 図51に示すように、診療科端末11には、特徴量算出部62と同様の特徴量算出部88が設けられる。特徴量算出部88は、例えば、診療科端末11にインストールされるソフトウエアをCPU41Aが実行することによって実現される。特徴量算出部88は、検査画像19を含む検査データ21及びGUI制御部53を通じて入力される関心領域ROIの領域情報に基づいて、特徴量RACを算出する。検索要求発行部54は、関心領域ROIに対応する画像と、算出した特徴量RACとを添付して、類似症例検索要求を発行する。
 図52に示すように、類似症例検索要求は、診療科端末11から類似症例検索サーバ17に送信される。類似症例検索サーバ17は、受信した類似症例検索要求に基づいて、類似検索を行い、検索結果を診療科端末11に配信する。第7実施形態においては、類似症例検索サーバ17の要求受付部61が特徴量取得部として機能する。
 上記各実施形態では、本発明の類似症例検索装置を、診療科端末11からの要求に基づいて類似症例検索を行う類似症例検索サーバ17の形態で説明したが、もちろん、類似症例検索サーバ17を使わずに、診療科端末11が症例DBサーバ16にアクセスして類似症例を検索するというように、診療科端末11に類似症例検索機能を設けてもよい。この場合には、診療科端末11が類似症例検索装置となる。
 また、上記各実施形態では、類似症例検索サーバ17と症例DBサーバ16を、別々のサーバで構成しているが、これらを統合して1つのサーバで構成してもよい。このように、1つのサーバに複数の機能を統合したり、機能毎にサーバを分離してもよい。
 なお、コンピュータシステムのハードウェア構成は種々の変形が可能である。例えば、類似症例検索サーバ17に関して、処理能力や信頼性の向上を目的として、ハードウェアとして分離された複数台のサーバコンピュータで構成することも可能である。このように、コンピュータシステムのハードウェア構成は、処理能力、安全性、信頼性など要求される性能に応じて適宜変更することができる。さらに、ハードウェアに限らず、症例DB23やAP50などのプログラムについて、安全性や信頼性の確保を目的として、二重化したり、あるいは、複数のストレージデバイスに分散して格納することももちろん可能である。
 また、上記各実施形態では、類似症例検索サーバ17について、1つの医療施設内で利用される形態で説明したが、複数の医療施設が利用できる形態としてもよい。
 具体的には、上記各実施形態では、類似症例検索サーバ17は、診療科端末11など1つの医療施設内に設置されるクライアント端末がLANを介して通信可能に接続され、クライアント端末からの要求に基づいて、類似症例検索に関するアプリケーションサービスを提供する形態である。複数の医療施設で利用可能にするには、類似症例検索サーバ17は、例えば、インターネットや公衆通信網などのWAN(Wide Area Network)を介して、複数の医療施設に設置されるクライアント端末と通信可能に接続される。そして、類似症例検索サーバ17は、複数の医療施設のクライアント端末からの要求を受け付けて、各クライアント端末に対して類似症例検索に関するアプリケーションサービスを提供する。
 この場合の類似症例検索サーバ17の設置場所や運営主体は、例えば医療施設とは別のデータセンタでもよいし、複数の医療施設のうちの1つでもよい。また、WANを利用する場合には、情報セキュリティを考慮して、VPN(Virtual Private Network)を構築したり、HTTPS(Hypertext Transfer Protocol Secure)などのセキュリティレベルの高い通信プロトコルを使用することが好ましい。
 本発明は、上記各実施形態に限らず、本発明の要旨を逸脱しない限り種々の構成を採り得ることはもちろんである。例えば、本例では、検査画像の例として、CT、MRI、単純X線画像を例としたが、マンモグラフィーや内視鏡など他のモダリティで撮影した検査画像に対しても適用することができる。また、上述の種々の実施形態や種々の変形例を適宜組み合わせることも可能である。また、本発明は、本発明を実現するプログラムに加えて、プログラムを記憶する記憶媒体にも及ぶ。

Claims (20)

  1.  症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する類似症例検索装置において、
     前記検査画像を1枚以上含む検査データにおいて、前記検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、前記関心領域毎の特徴量を取得する特徴量取得部と、
     前記関心領域と、前記症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた前記関心領域の特徴量と前記症例病変の特徴量とを比較して、前記関心領域毎の個別類似度をそれぞれ算出する個別類似度算出部と、
     複数の前記関心領域に対してそれぞれ算出した複数の前記個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の前記症例病変と前記複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、前記個別類似度を用いて総合類似度を算出する総合類似度算出部と、
     前記完全異種の組み合わせに対応する前記総合類似度に基づいて、前記類似症例を検索する類似症例検索部と、を備えている類似症例検索装置。
  2.  前記総合類似度算出部は、前記完全異種の組み合わせに関して、前記関心領域の個数と、前記症例病変の種類数及び個数とに応じた数の組み合わせを作り、前記完全異種の組み合わせ毎に前記総合類似度を算出することを特徴とする請求項1に記載の類似症例検索装置。
  3.  同一症例内に含まれる複数の異種の前記症例病変のうち少なくとも1種に関して、同種の前記症例病変が複数個存在する場合には、前記総合類似度算出部は、同種の複数個の前記症例病変を区別して、区別した同種の前記症例病変毎に前記完全異種の組み合わせを作ることを特徴とする請求項1に記載の類似症例検索装置。
  4.  前記類似症例検索部は、前記総合類似度に基づいて、複数件の前記類似症例に関する情報をリスト化した類似症例リストを作成することを特徴とする請求項1に記載の類似症例検索装置。
  5.  前記類似症例リストにおいて、前記類似症例は、前記総合類似度順にソートされることを特徴とする請求項4に記載の類似症例検索装置。
  6.  前記類似症例リストの表示項目には、前記総合類似度の値と前記総合類似度に関する内訳情報とを含み、前記内訳情報は、前記個別類似度を算出するための前記関心領域と前記症例病変との対応関係を含むことを特徴とする請求項4に記載の類似症例検索装置。
  7.  前記類似症例リストには、前記総合類似度の値に加えて、前記総合類似度の算出要素となる複数の前記個別類似度の値が表示されることを特徴とする請求項6に記載の類似症例検索装置。
  8.  前記類似症例リストには、前記関心領域と前記症例病変のそれぞれの画像が表示されることを特徴とする請求項4に記載の類似症例検索装置。
  9.  前記症例データベースには、前記症例病変の前記病変の種類の情報が格納されていることを特徴とする請求項1に記載の類似症例検索装置。
  10.  前記総合類似度は、前記完全異種の組み合わせに含まれる複数の前記個別類似度の総和であることを特徴とする請求項1に記載の類似症例検索装置。
  11.  前記個別類似度算出部は、前記関心領域毎に、複数の前記症例病変を対応させて算出した複数の前記個別類似度を記録した個別類似度テーブルを作成することを特徴とする請求項1に記載の類似症例検索装置。
  12.  前記総合類似度算出部は、前記関心領域毎の複数の前記個別類似度テーブルのそれぞれから、1つずつ前記個別類似度を読み出して、読み出した複数の前記個別類似度を要素として、前記完全異種の組み合わせを作ることを特徴とする請求項11に記載の類似症例検索装置。
  13.  前記個別類似度算出部は、前記症例のうち、少なくとも前記症例病変の種類が複数の症例については前記個別類似度の算出対象とし、前記症例病変の種類が1種類の症例については前記個別類似度の算出対象外とすることを特徴とする請求項1に記載の類似症例検索装置。
  14.  前記総合類似度算出部は、前記総合類似度を算出する要素となる前記個別類似度の値に応じて、前記総合類似度に対して重み付け処理を行うことを特徴とする請求項1に記載の類似症例検索装置。
  15.  前記重み付け処理は、前記個別類似度が閾値以上の場合に、前記総合類似度を高くする処理であることを特徴とする請求項14に記載の類似症例検索装置。
  16.  前記類似症例検索部は、前記症例病変の種類数が、前記関心領域の数未満の前記症例については、検索の対象外とすることを特徴とする請求項1に記載の類似症例検索装置。
  17.  前記類似症例検索部は、前記症例病変の種類数が、前記関心領域の数未満の前記症例についても、前記症例病変の種類数が複数有る場合には、検索の対象とすることを特徴とする請求項1に記載の類似症例検索装置。
  18.  1つの前記関心領域と、1件の症例に含まれる複数の前記症例病変とを対応させることにより、複数の前記総合類似度が算出された場合に、複数の前記総合類似度の中から、1つの代表値を判定する代表値判定部を有しており、
     前記類似症例検索部は、前記代表値に基づいて、類似症例を検索することを特徴とする請求項1に記載の類似症例検索装置。
  19.  症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する類似症例検索方法において、
     前記検査画像を1枚以上含む検査データにおいて、前記検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、前記関心領域毎の特徴量を取得する特徴量取得ステップと、
     前記関心領域と、前記症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた前記関心領域の特徴量と前記症例病変の特徴量とを比較して、前記関心領域毎の個別類似度をそれぞれ算出する個別類似度算出ステップと、
     複数の前記関心領域に対してそれぞれ算出した複数の前記個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の前記症例病変と前記複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、前記個別類似度を用いて総合類似度を算出する総合類似度算出ステップと、
     前記完全異種の組み合わせに対応する前記総合類似度に基づいて、前記類似症例を検索する類似症例検索ステップと、
     を備えている類似症例検索方法。
  20.  症例画像を1枚以上含む症例が複数件登録された症例データベースの中から、患者の診断に用いる検査画像に類似する類似症例を検索する処理をコンピュータに実行させる類似症例検索プログラムにおいて、
     前記検査画像を1枚以上含む検査データにおいて、前記検査画像内に存在する病変である対象病変を含むように指定される関心領域であり、異なる対象病変をそれぞれ1つ以上含む複数の関心領域に関して、前記関心領域毎の特徴量を取得する特徴量取得ステップと、
     前記関心領域と、前記症例画像内に存在する複数件の症例病変とを1対1で対応させて、対応させた前記関心領域の特徴量と前記症例病変の特徴量とを比較して、前記関心領域毎の個別類似度をそれぞれ算出する個別類似度算出ステップと、
     複数の前記関心領域に対してそれぞれ算出した複数の前記個別類似度に基づいて、総合類似度を算出する総合類似度算出部であり、同一症例内の種類が異なる複数の前記症例病変と前記複数の関心領域をそれぞれ組み合わせて構成される完全異種の組み合わせに対してのみ、前記個別類似度を用いて総合類似度を算出する総合類似度算出ステップと、
     前記完全異種の組み合わせに対応する前記類似症例を検索する類似症例検索ステップと、
     を前記コンピュータに実行させることを特徴とする類似症例検索プログラム。
PCT/JP2015/056370 2014-03-27 2015-03-04 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム WO2015146520A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201580016673.9A CN106164903B (zh) 2014-03-27 2015-03-04 类似病例检索装置、以及类似病例检索方法
US15/275,851 US10839511B2 (en) 2014-03-27 2016-09-26 Similar case search device, similar case search method, and non-transitory computer readable medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014066286A JP6099593B2 (ja) 2014-03-27 2014-03-27 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム
JP2014-066286 2014-03-27

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/275,851 Continuation US10839511B2 (en) 2014-03-27 2016-09-26 Similar case search device, similar case search method, and non-transitory computer readable medium

Publications (1)

Publication Number Publication Date
WO2015146520A1 true WO2015146520A1 (ja) 2015-10-01

Family

ID=54195050

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/056370 WO2015146520A1 (ja) 2014-03-27 2015-03-04 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム

Country Status (4)

Country Link
US (1) US10839511B2 (ja)
JP (1) JP6099593B2 (ja)
CN (1) CN106164903B (ja)
WO (1) WO2015146520A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113257431A (zh) * 2021-06-18 2021-08-13 武汉泰乐奇信息科技有限公司 一种基于虚拟人虚拟病例的自然人病例生成方法和系统
CN117009829A (zh) * 2023-10-07 2023-11-07 成都华栖云科技有限公司 一种录播课程的相似度识别方法及装置

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102656542B1 (ko) * 2015-12-22 2024-04-12 삼성메디슨 주식회사 초음파 영상들을 디스플레이하는 방법 및 장치.
JP6913697B2 (ja) * 2016-06-24 2021-08-04 ベックマン コールター, インコーポレイテッド 画像アトラスシステム及び方法
AU2017204494B2 (en) 2016-09-01 2019-06-13 Casio Computer Co., Ltd. Diagnosis assisting device, image processing method in diagnosis assisting device, and non-transitory storage medium having stored therein program
JP6390746B2 (ja) * 2016-09-01 2018-09-19 カシオ計算機株式会社 診断支援装置、及び診断支援装置における画像処理方法、並びにプログラム
JP6855850B2 (ja) * 2017-03-10 2021-04-07 富士通株式会社 類似症例画像検索プログラム、類似症例画像検索装置及び類似症例画像検索方法
JP6930302B2 (ja) * 2017-09-01 2021-09-01 富士通株式会社 診断支援用コンピュータプログラム、診断支援装置及び診断支援方法
US11139080B2 (en) 2017-12-20 2021-10-05 OrthoScience, Inc. System for decision management
JPWO2020044736A1 (ja) * 2018-08-31 2021-08-10 富士フイルム株式会社 類似度決定装置、方法およびプログラム
CN112638262B (zh) * 2018-08-31 2024-05-14 富士胶片株式会社 相似度确定装置、方法及程序
JP7091475B2 (ja) 2018-11-27 2022-06-27 富士フイルム株式会社 類似度決定装置、方法およびプログラム
DE112019005902T5 (de) * 2018-11-27 2021-09-30 Fujifilm Corporation Ähnlichkeitsbestimmungsvorrichtung, ähnlichkeitsbestimmungsverfahren und ähnlichkeitsbestimmungsprogramm
CN110457425B (zh) * 2019-08-05 2024-05-03 腾讯科技(深圳)有限公司 一种案例存储方法、装置、设备及存储介质
KR102516154B1 (ko) * 2019-12-31 2023-03-30 주식회사 코어라인소프트 정량화 파라미터를 이용하는 의료 영상 분석 시스템 및 유사 증례 검색 시스템, 및 그 시스템에서 실행되는 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000285141A (ja) * 1999-01-27 2000-10-13 Ricoh Co Ltd 画像検索装置,画像分類装置およびそれらの装置としてコンピュータを機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
WO2013065090A1 (ja) * 2011-11-04 2013-05-10 パナソニック株式会社 類似症例検索装置および類似症例検索方法
WO2015015699A1 (ja) * 2013-08-01 2015-02-05 パナソニック株式会社 類似症例検索装置、症例データベース、類似症例検索方法、およびプログラム

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5544352A (en) * 1993-06-14 1996-08-06 Libertech, Inc. Method and apparatus for indexing, searching and displaying data
JP4767759B2 (ja) * 2006-06-02 2011-09-07 富士フイルム株式会社 読影レポート作成装置
JP5153281B2 (ja) * 2007-09-28 2013-02-27 キヤノン株式会社 診断支援装置及びその制御方法
CN101373479A (zh) * 2008-09-27 2009-02-25 华中科技大学 一种乳腺x线摄片计算机图像检索方法及系统
JP5094775B2 (ja) 2009-03-31 2012-12-12 富士フイルム株式会社 症例画像検索装置、方法およびプログラム
US8582802B2 (en) * 2009-10-09 2013-11-12 Edgenet, Inc. Automatic method to generate product attributes based solely on product images
JP5618787B2 (ja) * 2010-11-29 2014-11-05 キヤノン株式会社 レポート作成支援装置及びその作成支援方法、並びにプログラム
JP5852970B2 (ja) * 2011-01-31 2016-02-03 パナソニック株式会社 症例検索装置および症例検索方法
CN102958425B (zh) * 2011-06-30 2015-09-23 松下电器产业株式会社 类似病例检索装置以及类似病例检索方法
JP4979842B1 (ja) * 2011-06-30 2012-07-18 パナソニック株式会社 類似症例検索装置および類似症例検索方法
WO2013018363A1 (ja) * 2011-08-04 2013-02-07 パナソニック株式会社 類似症例検索装置および類似症例検索方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000285141A (ja) * 1999-01-27 2000-10-13 Ricoh Co Ltd 画像検索装置,画像分類装置およびそれらの装置としてコンピュータを機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
WO2013065090A1 (ja) * 2011-11-04 2013-05-10 パナソニック株式会社 類似症例検索装置および類似症例検索方法
WO2015015699A1 (ja) * 2013-08-01 2015-02-05 パナソニック株式会社 類似症例検索装置、症例データベース、類似症例検索方法、およびプログラム

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113257431A (zh) * 2021-06-18 2021-08-13 武汉泰乐奇信息科技有限公司 一种基于虚拟人虚拟病例的自然人病例生成方法和系统
CN113257431B (zh) * 2021-06-18 2021-09-28 武汉泰乐奇信息科技有限公司 一种基于虚拟人虚拟病例的自然人病例生成方法和系统
CN117009829A (zh) * 2023-10-07 2023-11-07 成都华栖云科技有限公司 一种录播课程的相似度识别方法及装置

Also Published As

Publication number Publication date
JP6099593B2 (ja) 2017-03-22
JP2015191286A (ja) 2015-11-02
CN106164903B (zh) 2019-06-18
CN106164903A (zh) 2016-11-23
US10839511B2 (en) 2020-11-17
US20170011187A1 (en) 2017-01-12

Similar Documents

Publication Publication Date Title
JP6099593B2 (ja) 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム
JP6099592B2 (ja) 類似症例検索装置及び類似症例検索プログラム
JP6109778B2 (ja) 類似症例検索装置、類似症例検索方法、及び類似症例検索プログラム
US10229497B2 (en) Integration of medical software and advanced image processing
US8837794B2 (en) Medical image display apparatus, medical image display method, and medical image display program
Aisen et al. Automated storage and retrieval of thin-section CT images to assist diagnosis: system description and preliminary assessment
US8189888B2 (en) Medical reporting system, apparatus and method
JP4905967B2 (ja) 類似症例検索装置、方法、およびプログラム
JP5383431B2 (ja) 情報処理装置、情報処理方法及びプログラム
US7933437B2 (en) Information processing apparatus and its control method and data processing system
JP4976164B2 (ja) 類似症例検索装置、方法、およびプログラム
JP5661890B2 (ja) 情報処理装置、情報処理方法及びプログラム
JP2004005364A (ja) 類似画像検索システム
KR20140024788A (ko) 고급 멀티미디어 구조화 리포트
JP2008200373A (ja) 類似症例検索装置、方法、およびプログラム、ならびに、類似症例データベース登録装置、方法、およびプログラム
JP5631914B2 (ja) データベース検索装置、方法、及び、プログラム
JP6845071B2 (ja) 自動レイアウト装置および自動レイアウト方法並びに自動レイアウトプログラム
WO2012032441A2 (en) Clinical state timeline.
JP2008217426A (ja) 症例登録システム
EP3471106A1 (en) Method and system for supporting clinical decisions
WO2020105416A1 (ja) 医療情報表示装置、方法およびプログラム、並びに医療情報を表示するためのグラフィックユーザインターフェース

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15769033

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15769033

Country of ref document: EP

Kind code of ref document: A1