WO2012032940A1 - Dementia diagnosis support device and dementia diagnosis support method - Google Patents
Dementia diagnosis support device and dementia diagnosis support method Download PDFInfo
- Publication number
- WO2012032940A1 WO2012032940A1 PCT/JP2011/069116 JP2011069116W WO2012032940A1 WO 2012032940 A1 WO2012032940 A1 WO 2012032940A1 JP 2011069116 W JP2011069116 W JP 2011069116W WO 2012032940 A1 WO2012032940 A1 WO 2012032940A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- region
- brain
- image
- graph
- dementia
- Prior art date
Links
- 206010012289 Dementia Diseases 0.000 title claims abstract description 104
- 238000003745 diagnosis Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims description 28
- 210000004556 brain Anatomy 0.000 claims abstract description 142
- 238000004364 calculation method Methods 0.000 claims description 25
- 230000005484 gravity Effects 0.000 claims description 11
- 230000000052 comparative effect Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 abstract description 8
- 210000001769 parahippocampal gyrus Anatomy 0.000 abstract description 5
- 238000010191 image analysis Methods 0.000 abstract description 3
- 210000003140 lateral ventricle Anatomy 0.000 abstract description 2
- 230000002861 ventricular Effects 0.000 description 62
- 238000012545 processing Methods 0.000 description 22
- 206010061216 Infarction Diseases 0.000 description 16
- 230000007574 infarction Effects 0.000 description 16
- 238000004891 communication Methods 0.000 description 13
- 238000007792 addition Methods 0.000 description 9
- 208000024827 Alzheimer disease Diseases 0.000 description 8
- 230000002490 cerebral effect Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 210000003625 skull Anatomy 0.000 description 6
- 230000008034 disappearance Effects 0.000 description 5
- 206010003694 Atrophy Diseases 0.000 description 4
- 230000037444 atrophy Effects 0.000 description 4
- 208000013677 cerebrovascular dementia Diseases 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 210000004885 white matter Anatomy 0.000 description 3
- 208000024806 Brain atrophy Diseases 0.000 description 2
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 210000004884 grey matter Anatomy 0.000 description 2
- 210000001320 hippocampus Anatomy 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 210000003478 temporal lobe Anatomy 0.000 description 2
- 102000013455 Amyloid beta-Peptides Human genes 0.000 description 1
- 108010090849 Amyloid beta-Peptides Proteins 0.000 description 1
- 208000030518 CARASIL syndrome Diseases 0.000 description 1
- 206010058558 Hypoperfusion Diseases 0.000 description 1
- 208000004552 Lacunar Stroke Diseases 0.000 description 1
- 201000004810 Vascular dementia Diseases 0.000 description 1
- 206010063661 Vascular encephalopathy Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000004958 brain cell Anatomy 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000002008 hemorrhagic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007917 intracranial administration Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4088—Diagnosing of monitoring cognitive diseases, e.g. Alzheimer, prion diseases or dementia
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/501—Apparatus 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 the head, e.g. neuroimaging or craniography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
Definitions
- the present invention relates to a dementia diagnosis support apparatus and a dementia diagnosis support method that support dementia diagnosis using medical images such as CT images, MR images, and US images.
- Non-Patent Document 1 Non-Patent Document 2
- VD cerebrovascular dementia
- Alzheimer's dementia which occurs when a substance called beta amyloid develops in brain cells.
- VD is divided into multiple infarct dementia (large vessel occlusion), dementia due to small vessel lesions (small vessel occlusion), hypoperfusion cerebrovascular dementia, cerebral hemorrhagic cerebrovascular dementia, etc. It is done. Furthermore, multiple infarct dementia is further divided into multiple lacunar infarctions and Binswanger disease (progressive subcortical vascular encephalopathy).
- dementia can be divided into various types, but it is important to know which type of dementia belongs to first in order to perform appropriate treatment.
- dementia can be divided into various types, but it is important to know which type of dementia belongs to first in order to perform appropriate treatment.
- the parahippocampal gyrus is the earliest atrophy (Non-patent Document 3), but the parahippocampal gyrus has a small volume, and visual evaluation should be performed on CT images and MRI images. It was considered difficult.
- Patent Document 1 describes an apparatus for determining brain atrophy.
- the ratio of gray matter volume and white matter volume with respect to the entire brain is calculated using the brain MRI image, and the ratio between the intracranial volume and the white matter volume is also calculated. Determine atrophy.
- Non-Patent Document 3 describes an early stage Alzheimer-type dementia diagnosis support system called VSRAD (registered trademark).
- VSRAD registered trademark
- MRI images are used to perform numerical evaluation by comparing the degree of atrophy of the parahippocampal volume with that of a normal brain.
- Non-Patent Document 4 describes a method of diagnosing dementia by measuring the lateral ventricular lower angle using CT images, quantitatively evaluating the hippocampus and its surrounding atrophy.
- Patent Document 1 is a method that can be applied only to MRI images. Further, in order to use the VSRAD shown in Non-Patent Document 3 described above, it was premised on having an image (MRI image) in which the hippocampus was reflected first. Non-patent document 4 was also premised on having an image in which the lower angle of the lateral ventricle was shown.
- the present invention has been made in view of the above problems, and utilizes brain images in which evaluation target sites such as parahippocampal gyrus and lower ventricular angle are not clearly shown, and CT images widely used in clinical practice. It is an object of the present invention to provide a dementia diagnosis support apparatus and a dementia diagnosis support method that can perform image analysis of the brain and can effectively support the diagnosis of the type of dementia.
- the present invention relates to a focused density area extracting unit that extracts a pixel corresponding to a specific density range from an input brain image as a focused density area, and a distribution of the focused density area as a figure.
- This is a dementia diagnosis support device.
- a brain region extracting means for extracting a brain region from the inputted brain image is provided.
- the present invention also provides a focused density region extraction step for extracting pixels corresponding to a specific density range from the input brain image as a focused density region, and the distribution of the focused density region as a figure, a graph, or a numerical value.
- a dementia diagnosis support method comprising: an operation step for performing an operation for representing; and a display step for displaying a figure, graph, or numerical value representing the distribution of the concentration region of interest obtained by the operation means. It is.
- the present invention it is possible to analyze a brain image using a brain image in which an evaluation target site such as a parahippocampal gyrus or a lateral ventricular lower angle is not clearly shown or a CT image widely used in clinical practice, and further, dementia It is possible to provide a dementia diagnosis support apparatus and a dementia diagnosis support method capable of effectively supporting the types of diagnosis.
- the figure which shows the whole structure of the dementia diagnosis support apparatus 100 The flowchart explaining the procedure of the concentration distribution graphicalization process which the dementia diagnosis assistance apparatus 100 which concerns on this invention performs Examples of brain images (input images) by type of dementia Figure clearly showing the low density area of the input image (a) Image obtained by extracting the low density area, (b) Image obtained by removing the ventricle area from the low density area
- An example of brain region division Example showing the size of the low-concentration area from which the ventricular area has been removed as a figure for each divided area Display example in left-right overlay mode (multiple infarct type) Display example in left-right overlay mode (Alzheimer type) Display example in image difference mode Setting example of division center position Setting example of dividing line
- Other examples of brain region segmentation and other graphic display formats Examples of brain images (input images) by type of dementia Flow chart explaining the procedure of density distribution graph and digitization processing Examples of images to be processed (CT images and binarized images) Graph G1 showing CT value distribution against standardized radius by
- the image processing system 1 includes a display device 107, a dementia diagnosis support device 100 having an input device 109, an image database 111 connected to the dementia diagnosis support device 100 via a network 110, And a medical image photographing device 112.
- the dementia diagnosis support apparatus 100 is a computer that performs processing such as image generation and image analysis.
- a medical image processing apparatus installed in a hospital or the like is included.
- the dementia diagnosis support apparatus 100 includes a CPU (Central Processing Unit) 101, a main memory 102, a storage device 103, a communication interface (communication I / F) 104, a display memory 105, a mouse 108, etc.
- An interface (I / F) 106 with a device is provided, and each unit is connected via a bus 113.
- the CPU 101 calls a program stored in the main memory 102 or the storage device 103 to the work memory area on the RAM of the main memory 102 and executes it, and drives and controls each unit connected via the bus 113 to diagnose dementia.
- Various processes performed by the support device 100 are realized.
- the CPU 101 executes a density distribution graphing process (see FIG. 2), a density distribution graph / digitization process (see FIG. 15), which will be described later, and the like on the captured brain region image (hereinafter referred to as a brain image).
- the distribution of the concentration region of interest in the brain image is displayed as a figure, graph, or numerical value.
- the density distribution graphicizing process and the density distribution graph / digitizing process will be described in each embodiment described later.
- the main memory 102 is composed of ROM (Read Only Memory), RAM (Random Access Memory), and the like.
- ROM Read Only Memory
- RAM Random Access Memory
- the ROM permanently holds a computer boot program, a BIOS program, data, and the like.
- the RAM temporarily holds programs, data, and the like loaded from the ROM, the storage device 103, and the like, and includes a work area that the CPU 101 uses for performing various processes.
- the storage device 103 is a storage device that reads / writes data to / from an HDD (hard disk drive) or other recording medium, and stores programs executed by the CPU 101, data necessary for program execution, an OS (operating system), and the like. .
- As for the program a control program corresponding to the OS and an application program are stored. Each of these program codes is read by the CPU 101 as necessary, transferred to the RAM of the main memory 102, and executed as various means.
- the communication I / F 104 includes a communication control device, a communication port, and the like, and mediates communication between the dementia diagnosis support device 100 and the network 110.
- the communication I / F 104 performs communication control with the image database 111, another computer, or a medical image photographing apparatus 112 such as an X-ray CT apparatus or an MRI apparatus via the network 110.
- the I / F 106 is a port for connecting a peripheral device, and transmits / receives data to / from the peripheral device.
- a pointing device such as a mouse 108 or a stylus pen may be connected via the I / F 106.
- the display memory 105 is a buffer that temporarily stores display data input from the CPU 101.
- the accumulated display data is output to the display device 107 at a predetermined timing.
- the display device 107 includes a display device such as a liquid crystal panel and a CRT monitor, and a logic circuit for executing display processing in cooperation with the display device, and is connected to the CPU 101 via the display memory 105.
- the display device 107 displays the display data stored in the display memory 105 under the control of the CPU 101.
- the input device 109 is an input device such as a keyboard, for example, and outputs various instructions and information input by the operator to the CPU 101.
- the operator interactively operates the dementia diagnosis support apparatus 100 using external devices such as the display device 107, the input device 109, and the mouse 108.
- the network 110 includes various communication networks such as a LAN (Local Area Network), a WAN (Wide Area Network), an intranet, the Internet, etc., and communication between the image database 111, a server, other information devices, etc. and the dementia diagnosis support apparatus 100. Mediate the connection.
- LAN Local Area Network
- WAN Wide Area Network
- intranet the Internet
- the dementia diagnosis support apparatus 100 Mediate the connection.
- the image database 111 stores and stores image data captured by the medical image capturing device 112.
- the image database 111 is connected to the dementia diagnosis support apparatus 100 via the network 110, but the image database 111 is stored in, for example, the storage device 103 in the dementia diagnosis support apparatus 100. May be provided. [First embodiment] Next, a first embodiment of the operation of the dementia diagnosis support apparatus 100 according to the present invention will be described with reference to FIGS.
- the CPU 101 of the dementia diagnosis support apparatus 100 reads a program and data related to the concentration distribution graphic processing shown in FIG. 2 from the main memory 102, and executes processing based on this program and data.
- the brain image data to be calculated is acquired from the image database 111 or the like via the network 110 and the communication I / F 104 and stored in the storage device 103 of the dementia diagnosis support apparatus 100. It shall be.
- the brain image is a tomographic image (two-dimensional) or three-dimensional volume image of the head imaged by an X-ray CT apparatus, MRI apparatus, ultrasonic apparatus or the like.
- a two-dimensional axial CT image is used as a brain image, but a coronal image can also be used.
- the CPU 101 of the dementia diagnosis support apparatus 100 captures a brain image to be processed.
- FIG. 3 (a) is a head CT image 20a of a patient with multiple infarct type dementia, (b) is a head CT image 20b of a patient with Alzheimer type dementia, (c) is a head CT image 20c of a patient with mixed dementia. Indicates.
- the CPU 101 performs alignment using the skull region or the cerebrospinal fluid region (step S1). This is to match the shape of the brain image based on the brain image of each individual case.
- the skull region and cerebrospinal fluid region can be extracted from the CT value of the input image. Further, the image alignment processing is performed on the basis of a reference image given in advance or a past image of the same patient. By performing this preprocessing, it becomes possible to statistically and accurately analyze brain images of various cases.
- the CPU 101 extracts a focused density area from the preprocessed image.
- a low concentration region from which the ventricle region is removed is extracted as the concentration region of interest (step S2).
- the first CPU101 is a brain region of the preprocessed image thresholded to extract a lower concentration region than the predetermined threshold value t 2 (low concentration region). Then, regions indicated by thick lines in the images 21a, 21b, and 21c in FIG. 4, that is, the ventricle region 23 and other low-concentration regions 22 as shown in the image 24b ′ in FIG. 5 (a) are extracted. The CPU 101 further removes the ventricular region 23 from the extracted region and removes the ventricular region as shown in the image 24b of FIG.5 (b) (hereinafter, the ventricular removal low concentration region 22). Extract).
- the ventricular region has a particularly low CT value and can be easily recognized from anatomical information that it is located near the center of the brain.
- the CPU 101 executes the processing of steps S3 to S5 in FIG. 2 to display the distribution of the complex shape of the ventricular removal low concentration region 22 in a simple figure.
- the CPU 101 divides the brain region or the inside of the skull into a plurality of regions (step S3).
- the CPU 101 divides the brain region or the inside of the skull symmetrically by a center line 3a that divides at least the left brain and the right brain.
- the left-right symmetry here is based on the arrangement of the left and right brains and does not necessarily match the left-right symmetry on the image. Therefore, the center line 3a is, for example, a line connecting brain depressions by cerebral columns.
- the CPU 101 moves the center line 3a within the image plane by a predetermined angle around the center of the brain or the center point of the center line 3a (the line connecting the brain depressions by the cerebral column).
- FIG. 6 shows an example in which the brain region is divided into 16 by the center line 3a and the dividing lines 3b to 3h.
- the areas divided by the center line 3a and the dividing lines 3b to 3h are referred to as a divided area A, a divided area B,.
- the CPU 101 obtains the total area and average coordinates of the ventricular removal low concentration region 22 for each divided region (step S4). If the image to be processed is two-dimensional, the total area is obtained. If the image to be processed is a three-dimensional image, the total volume is obtained instead of the total area.
- the CPU 101 generates figures 5a to 5p having a size corresponding to the size of the total area (or total volume) calculated in step S4 and displays them in the divided areas A to B as shown in FIG.
- the display position is displayed at a position centered on the average coordinate of the ventricular removal low concentration region 22 in each of the divided regions A to B (step S5).
- step S5 the sizes of the figures 5a to 5p are the same as the total area calculated in step S4 or an area proportional to the total area (in the case of two dimensions). In the case of three dimensions, the volume is the same as the total volume or a volume proportional to the total volume.
- step S5 A display example of step S5 is shown in FIG.
- graphic display images 25a, 25b, and 25c are generated in which the graphic 5a to 5p indicating the size of the ventricular removal low concentration region 22 is displayed in each divided region of the brain image.
- the size of each figure 5a to 5p represents the size of the ventricular removal low concentration region 22 in the corresponding divided region
- the position of each figure 5a to 5p is the position of the ventricular removal low concentration region 22 in the corresponding divided region. Represents average coordinates.
- ventricular removal low concentration region 22 By generating and displaying such a graphic display image, it can be displayed in an easy-to-understand manner based on the position and size of the graphic 5a to 5p where and how much the ventricular removal low concentration region 22 is distributed in the brain region.
- the distribution of the ventricular removal low concentration region 22 is displayed as a graphic for each of the plurality of divided regions divided symmetrically, so the left-right symmetry of the ventricular removal low concentration region 22 is clearly displayed. it can.
- FIG. 7 (a) is a graphic display image 25a in the multiple infarct type
- FIG. 7 (b) is a graphic display image 25b in the Alzheimer type
- FIG. 7 (c) is a graphic display image 25c in the mixed type.
- the Alzheimer type of (b) is more symmetrical than the multiple infarct type of (a), and the focused concentration region (ventricular removal low concentration region 22) Can be seen at a glance.
- Such a difference in the distribution of the low-concentration region by type is useful information for determining the type of dementia.
- the dementia diagnosis support apparatus 100 includes a left-right superimposed display mode (see FIGS. 8 and 9), a left-right difference display mode, and an inter-image difference display mode (FIG. 10) as display modes different from the display example of FIG. Etc.). Further, as a user interface for switching between these display modes, the dementia diagnosis support apparatus 100 includes a “measurement” button 42, a “horizontal overlap” button 43, a “left / right difference” button 44, and an “image difference” button 45. It is desirable to provide software buttons including the like. Further, a “next patient number” button 46 for switching images and an “end” button 47 for ending the processing may be provided.
- the mode is switched to the horizontal overlay mode shown in FIGS.
- the CPU 101 selects the right region of the line (center line 3a in FIG. 8) connecting the brain depressions in the cerebral column among the figures 5a to 5p of the divided regions A to P calculated in the measurement mode.
- the graphic in the left area is duplicated and displayed at a position symmetrical to the center line 3a.
- the superimposed graphic hereinafter referred to as the superimposed graphic
- the original graphic in different display formats (color, pattern, graphic, etc.).
- a multiple infarct type graphic display image 25a and a left-right superimposed image 26a of the graphic display image 25a are displayed side by side.
- the figure 5p in the right area P is displayed as a superimposed figure 6p in the symmetry area A with the center line 3a as the symmetry axis
- the figure 5o in the right area O is in the symmetry area B.
- the figure 5n in the right area N is displayed as the superimposed figure 6n in the symmetrical area C
- the figure 5m in the right area M is displayed as the superimposed figure 6m in the symmetrical area D.
- the graphic 5l in the right region L is displayed as a superimposed graphic 6l in the symmetric region E
- the graphic 5i in the right region I is displayed as a superimposed graphic 6i in the symmetric region H.
- the superimposed figures 6i to 6p displayed at the left and right symmetrical positions are displayed in a different display format from the original figures 5a to 5p.
- the mode is switched to the left / right difference mode.
- the CPU 101 sets the regions 5a to 5p of each divided region calculated in the measurement mode to regions that are bilaterally symmetrical with respect to the line connecting the brain depressions by the cerebral tandem (center line 3a in FIG. 8).
- a difference in size of each figure (hereinafter, left-right difference) is taken, and a figure indicating the left-right difference in size (hereinafter, referred to as a left-right difference graphic) is displayed (not shown).
- the shape of the left-right difference graphic is arbitrary, but the size is the same as the size of the left-right difference or a size proportional to the size of the left-right difference.
- the display position of the left / right difference graphic is the average position of the display positions of the respective figures at the symmetrical positions (average coordinates in each divided region of the ventricular removal low concentration region 22). Or you may display a left-right difference figure in the position where the big figure is displayed among each figure in a left-right symmetric position.
- the mode is switched to the inter-image difference mode.
- the CPU 101 determines the distribution of the ventricular removal low concentration region 22 of the input image and the distribution of the ventricular removal low concentration region 22 of the comparison target image in the same manner as Steps S1 to S5 described above.
- Each of these images is measured, and the difference between the ventricular removal low concentration regions 22 (hereinafter referred to as inter-image differences) is taken for each corresponding divided region, and a figure corresponding to the size of the inter-image difference (hereinafter referred to as the inter-image difference).
- This is called an inter-image difference graphic (corresponding to 7b to 7o in FIG.
- the shape of the inter-image difference graphic 7b to 7o is arbitrary, but the size thereof is the same as the size of the inter-image difference or a size proportional to the size of the inter-image difference. In order to prevent confusion, it is desirable to display the inter-image difference graphics 7b to 7o and the original graphics 5a to 5p and 5A to 5P in different display formats (colors, patterns, graphics, etc.).
- the display positions of the inter-image difference graphics 7b to 7o are the average positions of the respective figures in the corresponding areas of the past image 25p and the current image 25P (for example, the average position of the display position of the graphic 5a and the display position of the graphic 5A). To do.
- FIG. 10 is an example of a display screen in the inter-image difference mode.
- a coronal image is a measurement target.
- a graphic display image 25p for the past image of the same patient and a graphic display image 25P for the current image are displayed side by side.
- An inter-image difference image 28 obtained by taking the inter-image difference of the distribution of the ventricular removal low-concentration region 22 is displayed for each corresponding divided region of these graphic display images 25p and 25P.
- the corresponding inter-image difference graphic is not displayed. The same applies to the figure 5p and the figure 5P.
- the image for which the difference is taken is not limited to the past image and the current image, and may be a reference image that is stored in advance and a comparison image that is input.
- the reference image is preferably stored in the storage device 103 of the dementia diagnosis support apparatus 100 for each type of dementia.
- a division number input field 48 may be provided as a user interface for designating the division number of the divided area.
- the CPU 101 divides the brain region into the division number input in the division number input field 48 in the above-described step S3 (FIG. 2).
- the CPU 101 of the dementia diagnosis support apparatus 100 extracts the brain region from the input brain image, and the attention belonging to a specific concentration range from the extracted brain region A concentration region (here, a low concentration region from which the ventricle region has been removed; a ventricular removal low concentration region 22) is extracted.
- the CPU 101 displays the distribution of the ventricular removal low-concentration region 22 as simple graphics 5a to 5p in the original image.
- the graphics 5a to 5p are respectively displayed in the divided areas A to P obtained by dividing the brain area based on the symmetry of the left and right brains.
- the sizes of figures 5a to 5p indicate the size of the ventricular removal low concentration area 22 in each divided area A to P, and the display positions of figures 5a to 5p are the ventricular removal in each divided area A to P.
- the average position of the low concentration region 22 is used.
- the distribution of the ventricular removal low concentration region 22 is displayed in a graphic form for each of the divided regions A to P divided symmetrically, so that the left / right symmetry of the ventricular removal low concentration region 22 is increased. Can be clearly displayed. Therefore, if a doctor or the like uses the dementia diagnosis support apparatus 100 of the present invention for diagnosis of dementia, the distribution of the ventricular removal low-concentration region 22 can be confirmed in an easy-to-understand manner, and the discrimination of the type of dementia is particularly effective. Can help. In addition, since the ventricular removal low concentration region 22 can be easily extracted from the CT image, it becomes possible to perform dementia diagnosis using CT images widely used in clinical practice.
- the dementia diagnosis support apparatus 100 displays a graphic displayed on either the left or right brain symmetrically with respect to the center line 3a dividing the left and right brain.
- the ventricular removal low concentration region 22 is distributed symmetrically compared to the case of multiple infarct dementia, so pay attention to the left-right symmetry of the ventricular removal low concentration region 22 It is also effective for the diagnosis of dementia type.
- the difference in size (left / right difference) of each figure in the left / right symmetrical area with respect to the center line 3a dividing the left / right brain is taken.
- a right-and-left differential graphic having a size proportional to the size is generated and displayed.
- a graphic display image is displayed for each of the reference brain image (for example, the past image of the patient) and the comparative brain image (for example, the current image of the patient), and the reference brain image and the comparative brain image are also displayed.
- the number of divisions of the brain region is 16, but the number is not limited to 16.
- the center of the division may be not only the center of gravity of the brain region but also any point on the line connecting the brain depressions by the cerebral column.
- it may be divided into portions effective for diagnosis from an anatomical viewpoint. In this case as well, it is desirable to divide the brain region so that at least the symmetry of the left and right brains can be seen.
- a handle H0 for setting the center point of the division
- a center line for setting the division
- the handles H0, H1, H2a, and H2b can be moved by an operation using a pointing device such as the mouse 108.
- the center line can be manually set based on the judgment of a doctor or the like. If the initial position of the handle H1 is a straight line connecting the brain depressions by the cerebral tandem lines detected by the CPU 101, the setting of the center line becomes easy and flexible.
- Handle H0 moves on the center line.
- the position of the handle H0 moved by the operation of the operator is set as the position of the center point.
- Handles H2a and H2b are handles that are moved relative to one another, and a straight line connecting the positions of the handles H2a and H2b and the center point (the position of the handle H0) is set as a dividing line.
- the positions of the handles H2a and H2b are interlocked so as to be always symmetrical with respect to the center line.
- H2a and H2b for setting a division line are provided in a number corresponding to the number of divisions.
- the dividing lines 3b to 3h of the region may be not only a straight line but also a curved line.
- a handle h1, h2, h3, h4 may be provided on each dividing line as a user interface 292 for setting the shape of the dividing line.
- the handles h1, h2, h3, and h4 can be moved by an operation using a pointing device such as the mouse 108.
- the shape of the line (partition line) connecting the handle H2a and the center line can be transformed into an arbitrary curve.
- a plurality of handles similar to h1, h2, h3, and h4 may be provided on each dividing line.
- it is desirable that the corresponding handles are interlocked so that the left and right dividing lines are symmetrical with respect to the center line.
- the dividing line can be set to an appropriate position and shape according to the cross section (axial, coronal) of the image, or the dividing line can be set to an appropriate position and shape based on anatomical information.
- the graphic representing the distribution of the ventricular removal low concentration region 22 is not limited to a circle.
- the graphics 5A to 5F may be rectangular or other shapes. Also, the pattern and color of the figure are arbitrary.
- FIG. 14 (a) is a Binswanger type, (b) is a multiple infarct type, (c) is an Alzheimer type, and (d) is a typical brain image of mixed type dementia. A region surrounded by a central line in each image in FIG. 14 is a ventricle region.
- low-concentration regions other than ventricular region 23 are densely distributed throughout
- low-concentration regions other than ventricular region 23 are distributed.
- C In the Alzheimer type, low-concentration regions other than the ventricular region 23 are distributed around the edge of the brain region.
- low-concentration regions other than the ventricular region 23 are distributed. It is distributed near the ventricle and the edge of the brain area. That is, the type of dementia depends not only on the left-right symmetry of the ventricular removal low concentration region, but also on how the ventricular removal low concentration region spreads.
- the characteristic concentration distribution of such an image is graphed or digitized so that a doctor or the like can easily determine which type it belongs to. Process.
- the CPU 101 of the dementia diagnosis support apparatus 100 reads out the program and data related to the concentration distribution graph / digitization processing shown in FIG. 15 from the main memory 102, and based on this program and data, A target density area is extracted, a calculation for graphing and digitizing the distribution of the extracted target density area is performed, and a calculation result is displayed.
- the brain image data to be calculated is acquired from the image database 111 or the like via the network 110 and the communication I / F 104 and stored in the storage device 103 of the dementia diagnosis support apparatus 100. It shall be.
- the brain image is a tomographic image (two-dimensional) or three-dimensional volume image of the head imaged by an X-ray CT apparatus, MRI apparatus, ultrasonic apparatus or the like.
- a two-dimensional axial CT image is used as a brain image, but the present invention can also be applied to a coronal image.
- the CPU 101 of the dementia diagnosis support apparatus 100 captures a brain image to be processed. Further, the CPU 101 generates an addition memory in the RAM work memory, and initializes (zero clear) (step S21).
- the addition memory is a one-dimensional array memory.
- the CPU 101 extracts a brain region based on the CT value of the captured brain image (input image) (step S22). Further, the CPU 101 obtains the center of gravity of the extracted brain region (step S23). Next, the CPU 101 recognizes the ventricle region of the input image (step S24).
- the ventricular region has a particularly low CT value and can be easily recognized from anatomical information that it is located near the center of the brain. The ventricular region is removed from the brain region by the processing as described above.
- the concentration distribution is analyzed with the region excluding the ventricular region from the brain region as the processing target.
- a CT image of a region obtained by removing the ventricle region from the brain region as shown in FIG. 16 (a) or only a low concentration region of interest as shown in FIG. 16 (b) is extracted. Any of binarized binarized images may be used.
- the CT value of each pixel is a calculation target of a profile described later, and in the binarized image in FIG. 16 (b), a binary value of “1” or “0” is described later.
- Profile calculation target This binarized image represents a low density pixel value “1” and the other pixels “0” in a region obtained by removing the ventricle from the brain region.
- the CPU 101 sets a moving radius with the center of gravity (or the above-mentioned center point) of the brain region obtained in step S23 as the rotation center (step S25).
- the CPU 101 obtains a concentration distribution profile along the moving radius at a certain rotation angle ⁇ (step S26).
- CT value is a calculation target as in the CT image of FIG. 16 (a)
- the CPU 101 acquires the pixel value (CT value) at each point on the radius vector as a one-dimensional array profile.
- CT value pixel value
- a binary value is to be calculated as in the binarized image of FIG. 16B
- the pixel value “1” or “0” at each point on the radius vector is acquired as a one-dimensional array profile.
- the CPU 101 determines the radius length at each rotation angle. Normalization is performed, the profile obtained in step S26 is interpolated, and the normalized profile is calculated (step S27). That is, when the radial length of each rotational angle is normalized with reference to the radial radius at a certain rotational angle, the profile obtained at a short angle of the total length of the radial radius (distance from the origin to the brain region end) Since the number of data to be acquired (number of pixels) is small compared to the profile obtained at the long angle of the total length, this is interpolated. For interpolation, the pixel value of the interpolated pixel may be estimated from the pixel value (CT value or binary) of the adjacent pixel.
- the CPU 101 adds the standardized profile to the addition memory (step S28). At this stage, the pixel value of each point on the radius vector at a certain rotation angle ⁇ is normalized and stored in the addition memory.
- the average CT value of the entire brain region may be subtracted from the normalized profile before addition to the addition memory. In this way, even if a CT value having a large value is set as a calculation target, the added value becomes small and overflow during calculation can be prevented. In addition, the profile change is easy to understand.
- step S29 the CPU 101 rotates the moving radius by a predetermined angle d ⁇ (step S29), and executes the processing of step S26 to step S29 in all angle directions (360 degrees).
- step S30 the CPU 101 divides the data (normalized profile) acquired in the addition memory by a predetermined value to obtain smoothed data. For example, if the data in the addition memory is divided by the number of additions, the data is smoothed over 360 degrees.
- the CPU 101 displays the obtained smoothed data as a graph (step S31).
- the graph represents the cumulative density value (CT value or binary value) with respect to the normalized radius.
- FIG. 17 is a graph G1 of the smoothed data (profile) with each CT image of FIG. 14 as the calculation target
- FIG. 18 is the smoothed data with the binary image as shown in FIG. (Profile) graph G2.
- the horizontal axis represents the normalized radius Rn
- the vertical axis represents the normalized average CT value (profile representing the relationship between the distance from the origin and the CT value at each normalized radius).
- the thin line represents the Binswanger type
- the broken line represents the multiple infarct type
- the alternate long and short dash line represents the Alzheimer type
- the thick line represents the mixed data.
- the Binswanger type has a low CT value as a whole
- the multiple infarct type has a low CT value concentrated in a region close to the center of the brain (region with a small radius)
- the Alzheimer type has a brain edge. It can be seen that low CT values are distributed in a concentrated manner (regions with a large radius), and the CT value is almost constant in any region in the mixed type.
- the horizontal axis is the normalized radius Rn
- the vertical axis is the count number (smoothed data; the value obtained by adding and smoothing the number of pixels with a normalized pixel value of 1 at each radial radius over 360 degrees. ).
- the thin line represents the Binswanger type
- the broken line represents the multiple infarct type
- the alternate long and short dash line represents the Alzheimer type
- the thick line represents the mixed data.
- the Binswanger type has a large count number regardless of radius (i.e., high density and low concentration as a whole), and the multiple infarct type has a count number near the center of the brain (part with a small radius)
- the number of counts (low density pixels) in the brain edge (large radius) is large, and in the mixed type, the count number is almost constant (low density overall). The ratio is small).
- the CPU 101 of the dementia diagnosis support apparatus 100 adds the CT value over 360 degrees for each point on the radial center with the center of the brain (e.g., the center of gravity) as the origin, or 360 only low-density pixels.
- the graphs G1 and G2 representing the relationship between the distance from the center of the brain and the density value (CT value or low-density pixel) are generated by accumulating (counting) over time. Since the density distribution of the image of the brain region differs depending on the type of dementia, it can be displayed in a form that can be easily compared with graphs G1 and G2, and the type of dementia can be easily determined.
- image features may be represented by numerical values.
- FIG. 19 shows each type of binarized image.
- Binswanger type the number of low-density pixels is large and distributed at high density. The average distance over which the low density pixels are distributed is medium, and the distribution spread (dispersion) is large.
- the multiple infarct type the number of low-density pixels is small and distributed at a low density. The average distance over which the low density pixels are distributed is small, and the distribution spread (dispersion) is also small.
- the mixed type the number of low-density pixels is small and distributed at a low density. The average distance over which the low density pixels are distributed is medium, and the distribution spread (distribution) is large.
- the CPU 101 calculates the feature of the concentration distribution as a feature amount (step S32 in FIG. 15).
- the CPU 101 displays the graph G1 in FIG. 17 and the graph G2 in FIG. 18 on the display device 107, and further calculates a reference feature amount to display as a numerical index 8.
- the correlation between the graph of the profile obtained from the input image and the reference graph calculated in advance from the reference image of each dementia type may be calculated and displayed as the type determination result 9.
- average distances R av , R CT , feature quantity M, disappearance rate, ventricular rate, etc. are calculated and displayed on the display screen.
- Rn is a normalized value of the distance from the center of the radius to the end of the brain region (or inside the skull).
- the average distance R CT is a value representing the characteristic of how the density value spreads from the center of the above-mentioned moving radius for the CT image of FIG. 16 (a), and is calculated by the following equation (1).
- CT max is the maximum CT value
- CT i is the CT value at pixel i. Also, calculation is performed only for pixels corresponding to the low density region from which the ventricle is removed (t1 ⁇ CT i ⁇ t2).
- Equation (1) is calculated for each radial angle ( ⁇ ). That is, the target pixel i is a pixel on the radius vector forming the pixel ⁇ .
- a graph of the calculation result of Expression (1) is shown in G3 of FIG.
- the characteristics of the density distribution can be expressed.
- a graph (reference graph G10) of the average distance R CT ( ⁇ ) in images of various cases may be displayed together with the graph G3.
- the concentration value deviation M is obtained by the following equation (2).
- CT max is the maximum CT value
- CT i is the CT value at pixel i. Also, calculation is performed only for pixels corresponding to the low density region from which the ventricle is removed (t1 ⁇ CT i ⁇ t2).
- the disappearance rate is a numerical value indicating the ratio of the low concentration region (concentration concentration region) in the brain region, and is calculated by the following equation (3).
- the ventricular rate is a numerical value indicating the proportion of the ventricle in the brain region, and is calculated by the following equation (4).
- the CPU 101 uses the graph G1 or G2 indicating the relationship between the normalized radius calculated for the input image and the concentration value, and the reference curve pattern for each type of dementia ( A graph indicating the relationship between the normalized radius and the density value) is compared with each other to determine which pattern the input image is close to (step S33; FIG. 15).
- the correlation between the graph and the reference pattern may be quantified from the Pearson product moment correlation coefficient r in Equation (6).
- the CPU 101 of the dementia diagnosis support apparatus extracts the brain region from the input brain image, removes the ventricle region from the extracted brain region, and the concentration of interest
- the pixel distribution (CT value distribution or low density region distribution) is analyzed and displayed in a graph or numerical form.
- CPU101 sets a radial radius centered on the center of gravity of the brain region, acquires a profile of the distance from the center and the concentration value (CT value or binary value) of each point on the radial radius over 360 degrees, A profile is obtained by sequentially adding density values (CT value or binary value) at each angle at the same distance from the center. At this time, the difference in the radial length of the profile at each angle is normalized.
- the points that lack data when normalized are interpolated by CT values or binary data of adjacent points on the radius vector.
- the CPU 101 calculates smoothed data obtained by dividing the acquired normalized profile by a predetermined value, generates graphs G1 and G2 representing the cumulative density values with respect to the normalized radius from the smoothed data, and displays them on the display screen. .
- the CPU 101 calculates the average distance R av and R CT that can be calculated from the density value of the image, the deviation M of the density value, the disappearance rate, the ventricular rate, etc.
- the feature amount is calculated and displayed. This makes it possible to recognize the type and progression of dementia quantitatively, and is effective for diagnosis.
- the CPU 101 compares the graph G1 or G2 indicating the relationship between the normalized radius calculated for the input image and the concentration value with the reference graph for each type of dementia, and the input image is Find which pattern is close. As a result, it is possible to easily confirm by numerical values which type the density distribution of the input image is close to and can effectively support the diagnosis of dementia.
- the low density pixels are extracted as “1” and the others as “0” as the binarized image.
- the high density pixels are extracted and the binarized image is extracted. You may get In this case, the ventricular region having a low CT value can be excluded without performing the ventricular region recognition process in step S24.
- a CT image is described as an example, but the present invention can also be applied to an MRI image or an ultrasound image.
- the “low density area” in the present embodiment is replaced with a “high density area”.
- graphic display image of the first embodiment described above and the graph, numerical index, and type determination result of the second embodiment may be displayed in combination.
- FIG. 1 image processing system, 100 dementia diagnosis support device, 101 CPU, 102 main memory, 103 storage device, 104 communication I / F, 105 display memory, 106 I / F, 107 display device, 108 mouse, 109 input device, 20a , 20b, 20c Input image, 22 Ventricular removal low concentration region, 23 Ventricular region, 24b Ventricular removal image, A, B, ..., P division region, 3a center line, 3b-3h division line, 25a, 25b, 25c figure display image, 5a-5p figure, 26a, 26b, 26c left-right superimposed image, 6a-6p left-right superimposed figure, 25p reference image, 25P comparison image, 28 image difference image, 7b-7o image difference Figure, H0 center point setting handle, H1 center line setting handle, H2a, H2b, H3a, H3b dividing line setting handle, h1 to h4 dividing line deformation handle, G1, G2 distance from the center of brain region and concentration value , 8 numeric
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Neurology (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Psychology (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Hospice & Palliative Care (AREA)
- Quality & Reliability (AREA)
- Psychiatry (AREA)
- Radiology & Medical Imaging (AREA)
- Neurosurgery (AREA)
- Physiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
Description
大別すると、脳血管障害に関連して出現した認知症を総称した脳血管性認知症(Vascular dementia:VD)と、ベータアミロイドと呼ばれる物質が脳細胞に発生することで起こるアルツハイマー認知症がある。 Recent research has revealed that there are various types of dementia (Non-Patent Document 1, Non-Patent Document 2).
Broadly classified, there are cerebrovascular dementia (VD), which is a general term for dementia that appears in connection with cerebrovascular disorders, and Alzheimer's dementia, which occurs when a substance called beta amyloid develops in brain cells. .
そして、好ましくは入力された脳画像から脳領域を抽出する脳領域抽出手段を備える。 In order to achieve the above-described object, the present invention relates to a focused density area extracting unit that extracts a pixel corresponding to a specific density range from an input brain image as a focused density area, and a distribution of the focused density area as a figure. , A graph or a numerical value, and a display means for displaying a figure, graph or numerical value representing the distribution of the concentration area of interest obtained by the calculating means. This is a dementia diagnosis support device.
Preferably, a brain region extracting means for extracting a brain region from the inputted brain image is provided.
まず、図1を参照して、本発明の認知症診断支援装置100を適用する画像処理システム1の構成について説明する。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
First, the configuration of an image processing system 1 to which the dementia diagnosis support
濃度分布図形化処理、及び濃度分布グラフ・数値化処理については後述する各実施の形態にて説明する。 In addition, the
The density distribution graphicizing process and the density distribution graph / digitizing process will be described in each embodiment described later.
[第1の実施の形態]
次に、図2~図13を参照して、本発明に係る認知症診断支援装置100の動作の第1実施の形態について説明する。 The
[First embodiment]
Next, a first embodiment of the operation of the dementia
この前処理を行うことで、様々な症例の脳画像を統計的に正確に解析することが可能となる。 The skull region and cerebrospinal fluid region can be extracted from the CT value of the input image. Further, the image alignment processing is performed on the basis of a reference image given in advance or a past image of the same patient.
By performing this preprocessing, it becomes possible to statistically and accurately analyze brain images of various cases.
ステップS3において、CPU101は、少なくとも左脳と右脳とを分割する中心線3aによって脳領域または頭蓋骨内部を左右対称に分割する。ここでいう左右対称とは、左右脳の配置に基づくものであり、必ずしも画像上の左右対称とは一致しなくてもよい。そのため、中心線3aは、例えば大脳縦列による脳の窪みを結ぶ線とする。更に、分割領域を増やす場合には、CPU101は、脳の重心または上記中心線3a(大脳縦列による脳の窪みを結ぶ線)の中点を中心として、中心線3aを所定角度ずつ画像平面内を回転させ、所望の分割数に分割する。図6は、脳領域を中心線3a及び分割線3b~3hにて16分割した例である。図6において、中心線3a、分割線3b~3hにより分割された各領域を図の上部から半時計回りに、分割領域A、分割領域B、・・・分割領域Pと呼ぶこととする。 First, the
In step S3, the
図7に示すように、脳画像の各分割領域内に、脳室除去低濃度領域22の大きさを示す図形5a~5pがそれぞれ表示された図形表示画像25a,25b,25cが生成される。各図形5a~5pの大きさは該当する分割領域にある脳室除去低濃度領域22の大きさを表し、各図形5a~5pの位置は該当する分割領域にある脳室除去低濃度領域22の平均座標を表している。 A display example of step S5 is shown in FIG.
As shown in FIG. 7,
なお、上述の例では、脳領域の分割数を16としたが、16分割に限定されない。 This makes it possible to easily compare the distribution of the ventricular removal low-
In the above example, the number of divisions of the brain region is 16, but the number is not limited to 16.
この場合、例えば図12に示すように、分割線の形状を設定するためのユーザインターフェース292として、各分割線上にハンドルh1,h2、h3、h4を備えてもよい。
ハンドルh1,h2、h3、h4は、マウス108等のポインティングデバイスによる操作にて移動可能である。 Further, the
In this case, for example, as shown in FIG. 12, a handle h1, h2, h3, h4 may be provided on each dividing line as a
The handles h1, h2, h3, and h4 can be moved by an operation using a pointing device such as the
[第2の実施の形態]
次に、図14~図22を参照して、本発明に係る認知症診断支援装置100の動作の第2の実施の形態について説明する。 Further, the graphic representing the distribution of the ventricular removal
[Second Embodiment]
Next, a second embodiment of the operation of the dementia
以上のような処理により、脳領域から脳室領域が除去される。 Next, the
The ventricular region is removed from the brain region by the processing as described above.
処理対象とする画像は、例えば、図16(a)のように脳領域から脳室領域を除いた領域のCT画像か、または図16(b)のように着目する低濃度領域のみを抽出し、二値化した二値化画像のいずれでもよい。 In the subsequent procedure, the concentration distribution is analyzed with the region excluding the ventricular region from the brain region as the processing target.
For the image to be processed, for example, a CT image of a region obtained by removing the ventricle region from the brain region as shown in FIG. 16 (a) or only a low concentration region of interest as shown in FIG. 16 (b) is extracted. Any of binarized binarized images may be used.
また、細線はBinswanger型、破線は多発梗塞型、一点鎖線はアルツハイマー型、太線は混合型の各平滑化データを表す。 In the graph G1 of FIG. 17, the horizontal axis represents the normalized radius Rn, and the vertical axis represents the normalized average CT value (profile representing the relationship between the distance from the origin and the CT value at each normalized radius).
The thin line represents the Binswanger type, the broken line represents the multiple infarct type, the alternate long and short dash line represents the Alzheimer type, and the thick line represents the mixed data.
また、グラフ以外にも、数値にて画像の特徴を表すようにしてもよい。 As described above, the
In addition to graphs, image features may be represented by numerical values.
図19を参照すると、(a)Binswanger型では、低濃度画素の画素数が多く、高密度に分布している。そして、低濃度画素が分布する平均距離は中程度であり、分布の広がり方(分散)は大きい。(b)多発梗塞型では、低濃度画素の画素数が少なく、低密度に分布している。そして、低濃度画素が分布する平均距離は小さく、分布の広がり方(分散)も小さい。(c)アルツハイマー型では、低濃度画素の画素数が少なく、低密度に分布している。そして、低濃度画素が分布する平均距離は大きく、分布の広がり方(分散)は小さい。(d)混合型では、低濃度画素の画素数が少なく、低密度に分布している。そして、低濃度画素が分布する平均距離は中程度で、分布の広がり方(分散)は大きい。 FIG. 19 shows each type of binarized image.
Referring to FIG. 19, in the (a) Binswanger type, the number of low-density pixels is large and distributed at high density. The average distance over which the low density pixels are distributed is medium, and the distribution spread (dispersion) is large. (b) In the multiple infarct type, the number of low-density pixels is small and distributed at a low density. The average distance over which the low density pixels are distributed is small, and the distribution spread (dispersion) is also small. (c) In the Alzheimer type, the number of low-density pixels is small and distributed at a low density. The average distance over which the low density pixels are distributed is large, and the spread (distribution) of the distribution is small. (d) In the mixed type, the number of low-density pixels is small and distributed at a low density. The average distance over which the low density pixels are distributed is medium, and the distribution spread (distribution) is large.
Claims (15)
- 入力された脳画像から特定の濃度範囲に該当する画素を着目濃度領域として抽出する着目濃度領域抽出手段と、
前記着目濃度領域の分布を、図形、グラフ、または数値にて表すための演算を行なう演算手段と、
前記演算手段によって求められた前記着目濃度領域の分布を表す図形、グラフ、または数値を表示する表示手段と、
を備えることを特徴とする認知症診断支援装置。 A focused density region extracting means for extracting a pixel corresponding to a specific density range from the input brain image as a focused density region;
A calculation means for performing a calculation for expressing the distribution of the concentration region of interest by a figure, a graph, or a numerical value;
Display means for displaying a figure, graph, or numerical value representing the distribution of the concentration region of interest obtained by the computing means;
A dementia diagnosis support apparatus comprising: - 前記入力された脳画像から脳領域を抽出する脳領域抽出手段を更に備えることを特徴とする請求項1に記載の認知症診断支援装置。 2. The dementia diagnosis support apparatus according to claim 1, further comprising brain region extraction means for extracting a brain region from the inputted brain image.
- 前記着目濃度領域の分布を前記図形にて表す場合、
前記演算手段は、
前記脳画像の脳領域を、左右脳の対称性に基づき複数の分割領域に分割する分割手段と、
前記着目濃度領域抽出手段により抽出された着目濃度領域の大きさを前記分割手段により分割された分割領域毎に計測し、計測された着目濃度領域の大きさを表す大きさの図形を分割領域毎に生成する図形生成手段と、を備え、
前記表示手段は、
前記図形生成手段によって生成された図形を該当する分割領域に夫々配置した図形表示画像を生成し、表示する図形表示手段を備えることを特徴とする請求項2に記載の認知症診断支援装置。 When the distribution of the concentration region of interest is represented by the graphic,
The computing means is
A dividing unit that divides the brain region of the brain image into a plurality of divided regions based on left-right brain symmetry;
The size of the focused density area extracted by the focused density area extracting unit is measured for each divided area divided by the dividing unit, and a figure having a size representing the measured size of the focused density area is measured for each divided area. And a graphic generation means for generating
The display means includes
3. The dementia diagnosis support apparatus according to claim 2, further comprising graphic display means for generating and displaying a graphic display image in which the graphic generated by the graphic generating means is arranged in a corresponding divided area. - 前記表示手段は、
前記図形生成手段によって生成された図形を該当する分割領域に夫々表示するとともに、前記図形のうち左右脳のいずれか一方に表示される図形を左右脳を分割する中心線に対して対称に表示する左右重ね合わせ表示手段を更に備えることを特徴とする請求項3に記載の認知症診断支援装置。 The display means includes
The graphic generated by the graphic generation means is displayed in the corresponding divided area, and the graphic displayed on either the left or right brain is displayed symmetrically with respect to the center line dividing the left and right brain. 4. The dementia diagnosis support apparatus according to claim 3, further comprising a left-right superimposed display means. - 前記表示手段は、
前記図形生成手段によって生成された図形を該当する分割領域に夫々表示するとともに、左右対称となる分割領域に表示される各図形の大きさの差分を表す左右差分図形を生成し、表示する左右差分表示手段を更に備えることを特徴とする請求項3に記載の認知症診断支援装置。 The display means includes
The left and right differences are generated by displaying the graphics generated by the graphic generation means in the corresponding divided areas, and generating and displaying the left and right differential figures that represent the difference in the size of each figure displayed in the symmetrical divided areas. 4. The dementia diagnosis support apparatus according to claim 3, further comprising display means. - 前記表示手段は、
基準脳画像についての前記図形表示画像と、比較脳画像についての前記図形表示画像と、を夫々表示するとともに、
前記基準脳画像及び前記比較脳画像についての各図形表示画像の対応する各分割領域について、前記着目濃度領域の大きさの差分を算出し、その大きさを表す画像間差分図形を生成し、該当する分割領域に表示する画像間差分表示手段を更に備えることを特徴とする請求項3に記載の認知症診断支援装置。 The display means includes
Displaying the graphic display image for the reference brain image and the graphic display image for the comparative brain image, respectively;
For each corresponding divided region of each graphic display image for the reference brain image and the comparative brain image, calculate the difference in the size of the concentration region of interest, and generate an inter-image difference graphic representing the size, 4. The dementia diagnosis support apparatus according to claim 3, further comprising an inter-image difference display means for displaying in the divided area. - 前記着目濃度領域の分布を前記グラフにて表す場合、
前記演算手段は、
前記脳領域の重心を原点とする動径を設定する動径設定手段と、
前記動径上の各点について、前記原点からの距離と濃度値との関係を表すプロファイルを全角度にわたって取得するプロファイル取得手段と、
前記プロファイル取得手段によって取得したプロファイルを基準動径長に基づいて規格化し、規格化したプロファイルについてのグラフを生成するグラフ生成手段と、を備え、 前記表示手段は、
前記グラフ生成手段により生成されたグラフを表示するグラフ表示手段と、
を備えることを特徴とする請求項2に記載の認知症診断支援装置。 When representing the distribution of the concentration region of interest in the graph,
The computing means is
A radial setting means for setting a radial with the center of gravity of the brain region as the origin;
Profile acquisition means for acquiring a profile representing the relationship between the distance from the origin and the concentration value over all angles for each point on the radius vector;
Normalizing the profile acquired by the profile acquisition unit based on a reference radial length, and generating a graph for the normalized profile, and the display unit includes:
Graph display means for displaying the graph generated by the graph generation means;
The dementia diagnosis support apparatus according to claim 2, further comprising: - 前記プロファイル取得手段は、前記脳領域から脳室領域を除去した領域の画素値を対象として前記プロファイルを取得することを特徴とする請求項7に記載の認知症診断支援装置。 8. The dementia diagnosis support apparatus according to claim 7, wherein the profile acquisition unit acquires the profile for a pixel value of a region obtained by removing a ventricle region from the brain region.
- 前記プロファイル取得手段は、前記着目濃度領域に該当するか否かを示す二値を対象として前記プロファイルを取得することを特徴とする請求項7に記載の認知症診断支援装置。 8. The dementia diagnosis support apparatus according to claim 7, wherein the profile acquisition unit acquires the profile for a binary value indicating whether or not the target concentration region is met.
- 認知症のタイプ別に取得され、規格化されたプロファイルをグラフ化した基準グラフを保持する基準グラフ保持手段と、
入力された脳画像について前記グラフ生成手段により生成されたグラフと、前記基準グラフとに基づき、認知症のタイプを判定するための相関値を算出する相関値算出手段と、 前記相関値算出手段により算出された相関値を、認知症のタイプ別に表示するタイプ表示手段と、
を更に備えることを特徴とする請求項7に記載の認知症診断支援装置。 A reference graph holding means for holding a reference graph obtained by graphing a standardized profile acquired for each type of dementia;
Correlation value calculation means for calculating a correlation value for determining the type of dementia based on the graph generated by the graph generation means for the input brain image and the reference graph, and by the correlation value calculation means Type display means for displaying the calculated correlation value for each type of dementia;
8. The dementia diagnosis support apparatus according to claim 7, further comprising: - 前記着目濃度領域の分布を前記数値にて表す場合、
前記演算手段は、前記脳領域の濃度分布の特徴を示す特徴量を演算する特徴量演算手段を備え、
前記表示手段は、前記特徴量演算手段により算出された特徴量を表示する数値指標表示手段を備えることを特徴とする請求項2に記載の認知症診断支援装置。 When expressing the distribution of the concentration region of interest by the numerical value,
The calculation means includes a feature quantity calculation means for calculating a feature quantity indicating a feature of the concentration distribution of the brain region,
3. The dementia diagnosis support apparatus according to claim 2, wherein the display unit includes a numerical index display unit that displays the feature amount calculated by the feature amount calculation unit. - 前記特徴量は、
前記脳領域の重心からの前記着目濃度領域の平均距離、濃度値の偏り、脳領域における着目濃度領域の割合、及び脳領域における脳室の割合のうちいずれか一つを含むことを特徴とする請求項11に記載の認知症診断支援装置。 The feature amount is
It includes any one of an average distance of the concentration region of interest from the center of gravity of the brain region, a deviation in concentration value, a proportion of the concentration region of interest in the brain region, and a proportion of the ventricle in the brain region. 12. The dementia diagnosis support apparatus according to claim 11. - 入力された脳画像から特定の濃度範囲に該当する画素を着目濃度領域として抽出する着目濃度領域抽出ステップと、
前記着目濃度領域の分布を、図形、グラフ、または数値にて表すための演算を行なう演算ステップと、
前記演算ステップによって求められた前記着目濃度領域の分布を表す図形、グラフ、または数値を表示する表示ステップと、
を備えることを特徴とする認知症診断支援方法。 A focused density region extraction step of extracting pixels corresponding to a specific density range from the input brain image as a focused density region;
A calculation step of performing a calculation for representing the distribution of the concentration region of interest by a figure, a graph, or a numerical value;
A display step for displaying a figure, a graph, or a numerical value representing the distribution of the concentration region of interest obtained by the calculation step;
A method for supporting the diagnosis of dementia, comprising: - 前記入力された脳画像から脳領域を抽出する脳領域抽出ステップを更に備え、
前記着目濃度領域の分布を前記図形にて表す場合、
前記演算ステップは、
前記脳画像の脳領域を、左右脳の対称性に基づき複数の分割領域に分割する分割ステップと、
前記着目濃度領域抽出ステップにより抽出された着目濃度領域の大きさを前記分割ステップにより分割された分割領域毎に計測し、計測された着目濃度領域の大きさを表す大きさの図形を分割領域毎に生成する図形生成ステップと、を備え、
前記表示ステップは、
前記図形生成ステップによって生成された図形を該当する分割領域に夫々配置した図形表示画像を生成し、表示する図形表示ステップを備えることを特徴とする請求項13に記載の認知症診断支援方法。 A brain region extracting step of extracting a brain region from the input brain image;
When the distribution of the concentration region of interest is represented by the graphic,
The calculation step includes:
A division step of dividing the brain region of the brain image into a plurality of divided regions based on left-right brain symmetry;
The size of the target density region extracted by the target density region extraction step is measured for each divided region divided by the division step, and a figure having a size representing the size of the measured target density region is measured for each divided region. And a graphic generation step for generating
The display step includes
14. The dementia diagnosis support method according to claim 13, further comprising a graphic display step of generating and displaying a graphic display image in which the graphic generated by the graphic generating step is arranged in a corresponding divided area. - 前記入力された脳画像から脳領域を抽出する脳領域抽出ステップを更に備え、
前記着目濃度領域の分布を前記グラフにて表す場合、
前記演算ステップは、
前記脳領域の重心を原点とする動径を設定する動径設定ステップと、
前記動径上の各点について、前記原点からの距離と濃度値との関係を表すプロファイルを全角度にわたって取得するプロファイル取得ステップと、
前記プロファイル取得ステップによって取得したプロファイルを基準動径長に基づいて規格化し、規格化したプロファイルについてのグラフを生成するグラフ生成ステップと、を備え、 前記表示ステップは、
前記グラフ生成ステップにより生成されたグラフを表示するグラフ表示ステップと、
を備えることを特徴とする請求項13に記載の認知症診断支援方法。 A brain region extracting step of extracting a brain region from the input brain image;
When representing the distribution of the concentration region of interest in the graph,
The calculation step includes:
A radial setting step for setting a radial with the center of gravity of the brain region as the origin;
For each point on the radius vector, a profile acquisition step of acquiring a profile representing the relationship between the distance from the origin and the density value over all angles;
A graph generation step of normalizing the profile acquired by the profile acquisition step based on a reference radial length and generating a graph of the normalized profile, and the display step includes
A graph display step for displaying the graph generated by the graph generation step;
The dementia diagnosis support method according to claim 13, comprising:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201180043149.2A CN103096787B (en) | 2010-09-07 | 2011-08-25 | Dementia diagnosis support device |
JP2012532930A JP5878125B2 (en) | 2010-09-07 | 2011-08-25 | Dementia diagnosis support apparatus and dementia diagnosis support method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010-199530 | 2010-09-07 | ||
JP2010199530 | 2010-09-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012032940A1 true WO2012032940A1 (en) | 2012-03-15 |
Family
ID=45810543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/069116 WO2012032940A1 (en) | 2010-09-07 | 2011-08-25 | Dementia diagnosis support device and dementia diagnosis support method |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP5878125B2 (en) |
CN (1) | CN103096787B (en) |
WO (1) | WO2012032940A1 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015043880A (en) * | 2013-08-28 | 2015-03-12 | コニカミノルタ株式会社 | Thoracic diagnosis support system |
WO2016047683A1 (en) * | 2014-09-25 | 2016-03-31 | 大日本印刷株式会社 | Medical image display processing method, medical image display processing device, and program |
JP2019136519A (en) * | 2019-04-08 | 2019-08-22 | 大日本印刷株式会社 | Medical image display processing method, device and program |
CN110390624A (en) * | 2019-07-25 | 2019-10-29 | 科大讯飞股份有限公司 | A kind of pattern evaluation method, device, equipment and storage medium |
EP3677180A4 (en) * | 2017-08-29 | 2020-10-14 | FUJIFILM Corporation | Medical information display device, method, and program |
JP2021037257A (en) * | 2019-09-05 | 2021-03-11 | キヤノンメディカルシステムズ株式会社 | Image processing apparatus, magnetic resonance imaging apparatus, and image processing method |
WO2022071158A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, method for operating diagnosis assistance device, program for operating diagnosis assistance device, dementia diagnosis assistance method, and learned model for deriving dementia findings |
WO2022071160A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, operation method of diagnosis assistance device, operation program of diagnosis assistance device, and dementia diagnosis assistance method |
WO2022071159A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, operation method for diagnosis assistance device, operation program for diagnosis assistance device, and dementia diagnosis assistance device |
JP2022547909A (en) * | 2019-09-05 | 2022-11-16 | コリア ユニバーシティ リサーチ アンド ビジネス ファウンデーション | Method and apparatus for predicting cerebral cortex contraction rate by region based on CT image |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102067335B1 (en) * | 2018-04-26 | 2020-01-16 | 인제대학교 산학협력단 | DlEMENTIA CLASSIFICATION METHOD BY DISTANCE ANALYSIS FROM THE CENTRAL CORONAL PLANE OF THE BRAIN HIPPOCAMPUS |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002325761A (en) * | 2000-06-30 | 2002-11-12 | Hitachi Medical Corp | Image diagnosis supporting device |
JP2002336221A (en) * | 2001-05-15 | 2002-11-26 | Fuji Photo Film Co Ltd | Abnormal shadow candidate detector |
JP2006136506A (en) * | 2004-11-12 | 2006-06-01 | Hitachi Medical Corp | Image processor |
WO2007114238A1 (en) * | 2006-03-30 | 2007-10-11 | National University Corporation Shizuoka University | Apparatus for determining brain atrophy, method of determining brain atrophy and program for determining brain atrophy |
WO2009005013A1 (en) * | 2007-06-29 | 2009-01-08 | Toshinori Kato | White matter enhancing device, white matter enhancing method, and program |
JP2010110567A (en) * | 2008-11-10 | 2010-05-20 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7245754B2 (en) * | 2000-06-30 | 2007-07-17 | Hitachi Medical Corporation | image diagnosis supporting device |
-
2011
- 2011-08-25 JP JP2012532930A patent/JP5878125B2/en not_active Expired - Fee Related
- 2011-08-25 WO PCT/JP2011/069116 patent/WO2012032940A1/en active Application Filing
- 2011-08-25 CN CN201180043149.2A patent/CN103096787B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002325761A (en) * | 2000-06-30 | 2002-11-12 | Hitachi Medical Corp | Image diagnosis supporting device |
JP2002336221A (en) * | 2001-05-15 | 2002-11-26 | Fuji Photo Film Co Ltd | Abnormal shadow candidate detector |
JP2006136506A (en) * | 2004-11-12 | 2006-06-01 | Hitachi Medical Corp | Image processor |
WO2007114238A1 (en) * | 2006-03-30 | 2007-10-11 | National University Corporation Shizuoka University | Apparatus for determining brain atrophy, method of determining brain atrophy and program for determining brain atrophy |
WO2009005013A1 (en) * | 2007-06-29 | 2009-01-08 | Toshinori Kato | White matter enhancing device, white matter enhancing method, and program |
JP2010110567A (en) * | 2008-11-10 | 2010-05-20 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
Non-Patent Citations (2)
Title |
---|
HIROSHI MATSUDA: "Voxel-based Specific Regional Analysis System as an Adjunct to Diagnosis of Early Alzheimer's Disease", JAPANESE JOURNAL OF RADIOLOGICAL TECHNOLOGY, vol. 62, no. 8, 2005, pages 1066 - 1072 * |
MASUMI HATTORI: "Tobu CT Gazo o Mochiita Sokutoyo Naisokubu Ishuku no Jido Keisokuho no Kaihatsu", JAPANESE JOURNAL OF RADIOLOGICAL TECHNOLOGY, vol. 60, no. 7, 2004, pages 993 - 998 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015043880A (en) * | 2013-08-28 | 2015-03-12 | コニカミノルタ株式会社 | Thoracic diagnosis support system |
WO2016047683A1 (en) * | 2014-09-25 | 2016-03-31 | 大日本印刷株式会社 | Medical image display processing method, medical image display processing device, and program |
JP2016064004A (en) * | 2014-09-25 | 2016-04-28 | 大日本印刷株式会社 | Medical image display processing method, medical image display processing device and program |
US10285657B2 (en) | 2014-09-25 | 2019-05-14 | Dai Nippon Printing Co., Ltd. | Medical image display processing method, medical image display processing device, and program |
EP3677180A4 (en) * | 2017-08-29 | 2020-10-14 | FUJIFILM Corporation | Medical information display device, method, and program |
US11295442B2 (en) | 2017-08-29 | 2022-04-05 | Fujifilm Corporation | Medical information display apparatus displaying cavity region in brain image, medical information display method, and medical information display program |
JP2019136519A (en) * | 2019-04-08 | 2019-08-22 | 大日本印刷株式会社 | Medical image display processing method, device and program |
CN110390624A (en) * | 2019-07-25 | 2019-10-29 | 科大讯飞股份有限公司 | A kind of pattern evaluation method, device, equipment and storage medium |
CN110390624B (en) * | 2019-07-25 | 2023-05-30 | 科大讯飞股份有限公司 | Graph evaluation method, device, equipment and storage medium |
JP2021037257A (en) * | 2019-09-05 | 2021-03-11 | キヤノンメディカルシステムズ株式会社 | Image processing apparatus, magnetic resonance imaging apparatus, and image processing method |
JP2022547909A (en) * | 2019-09-05 | 2022-11-16 | コリア ユニバーシティ リサーチ アンド ビジネス ファウンデーション | Method and apparatus for predicting cerebral cortex contraction rate by region based on CT image |
JP7317306B2 (en) | 2019-09-05 | 2023-07-31 | コリア ユニバーシティ リサーチ アンド ビジネス ファウンデーション | Method and apparatus for predicting cerebral cortex contraction rate by region based on CT image |
JP7374786B2 (en) | 2019-09-05 | 2023-11-07 | キヤノンメディカルシステムズ株式会社 | Image processing device, magnetic resonance imaging device, and image processing method |
WO2022071158A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, method for operating diagnosis assistance device, program for operating diagnosis assistance device, dementia diagnosis assistance method, and learned model for deriving dementia findings |
WO2022071160A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, operation method of diagnosis assistance device, operation program of diagnosis assistance device, and dementia diagnosis assistance method |
WO2022071159A1 (en) * | 2020-10-01 | 2022-04-07 | 富士フイルム株式会社 | Diagnosis assistance device, operation method for diagnosis assistance device, operation program for diagnosis assistance device, and dementia diagnosis assistance device |
Also Published As
Publication number | Publication date |
---|---|
JP5878125B2 (en) | 2016-03-08 |
JPWO2012032940A1 (en) | 2014-01-20 |
CN103096787B (en) | 2014-12-31 |
CN103096787A (en) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5878125B2 (en) | Dementia diagnosis support apparatus and dementia diagnosis support method | |
JP5970766B2 (en) | Medical image processing apparatus, medical image processing method, and program | |
JP4937922B2 (en) | Image analysis apparatus and method | |
JP5601378B2 (en) | Medical image processing device | |
JP6036009B2 (en) | Medical image processing apparatus and program | |
US20050283070A1 (en) | Systems and methods for qualifying symmetry to evaluate medical images | |
JP2005237441A (en) | Method and apparatus for supporting diagnosis of encephalopathy | |
WO2007114238A1 (en) | Apparatus for determining brain atrophy, method of determining brain atrophy and program for determining brain atrophy | |
JP5098393B2 (en) | Region of interest determination device | |
US20130303900A1 (en) | Method and apparatus for processing of stroke ct scans | |
Qiu et al. | User-guided segmentation of preterm neonate ventricular system from 3-D ultrasound images using convex optimization | |
JP2018068400A (en) | Dynamic image processing device | |
JP6703323B2 (en) | ROI setting technology for biological image inspection | |
JP6442311B2 (en) | Technology for extracting tumor contours in nuclear medicine images | |
JP5107538B2 (en) | Diagnostic imaging support system and method | |
JP2009226043A (en) | Medical image processor and method for detecting abnormal shadow | |
US9218657B2 (en) | Method of obtaining and analyzing data from an upright MRI from the spinal region of a subject | |
JP6442309B2 (en) | Nuclear medicine image analysis technology | |
Liao et al. | Tracing the deformed midline on brain CT | |
WO2022059596A1 (en) | Skin surface analysis device and skin surface analysis method | |
JP6442310B2 (en) | Technology for extracting tumor regions from nuclear medicine images | |
KR102349360B1 (en) | Diagnosis Method and System of Idiopathic Normal Pressure Hydrocephalus Using Imaging Diagnostic Equipment | |
JP7154012B2 (en) | Image processing program, image processing apparatus, and image processing method | |
JP6286606B1 (en) | Image processing program, image processing apparatus, and image processing method | |
WO2016125349A1 (en) | Nuclear medical image analysis technique |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201180043149.2 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11823419 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012532930 Country of ref document: JP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 11823419 Country of ref document: EP Kind code of ref document: A1 |