US20030228040A1 - Anatomical characteristic position detecting system, and object structure measurement system - Google Patents

Anatomical characteristic position detecting system, and object structure measurement system Download PDF

Info

Publication number
US20030228040A1
US20030228040A1 US10/455,417 US45541703A US2003228040A1 US 20030228040 A1 US20030228040 A1 US 20030228040A1 US 45541703 A US45541703 A US 45541703A US 2003228040 A1 US2003228040 A1 US 2003228040A1
Authority
US
United States
Prior art keywords
artificial image
shape information
characteristic position
artificial
anatomical
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
US10/455,417
Other versions
US7315635B2 (en
Inventor
Akira Oosawa
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Fujifilm Corp
Original Assignee
Fuji Photo Film Co Ltd
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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Assigned to FUJI PHOTO FILM CO., LTD. reassignment FUJI PHOTO FILM CO., LTD. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: OOSAWA, AKIRA
Publication of US20030228040A1 publication Critical patent/US20030228040A1/en
Assigned to FUJIFILM CORPORATION reassignment FUJIFILM CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: FUJIFILM HOLDINGS CORPORATION (FORMERLY FUJI PHOTO FILM CO., LTD.)
Application granted granted Critical
Publication of US7315635B2 publication Critical patent/US7315635B2/en
Expired - Fee Related legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • 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/60Analysis of geometric attributes

Definitions

  • This invention relates to an anatomical characteristic position detecting system and an object structure measurement system for aiding in diagnosing anatomical structure of an object on the basis of a medical image.
  • the profile of the cardiothorax is automatically detected by template matching on the plane of the polar coordinate system by the use of a template which is substantially similar to the average cardiothorax profile (a reference profile), and the cardiothoracic ratio is calculated on the basis of the automatically detected profile of the cardiothorax.
  • the primary object of the present invention is to provide an anatomical characteristic position detecting system and an object structure measurement system which make it feasible to detect anatomical characteristic positions and to measure structural physical values of an object on the basis of information on the detected characteristic positions even if the image has a part which is unclear to know the anatomical structure of the object.
  • an anatomical characteristic position detecting system comprising
  • a shape information obtaining means which obtains shape information for use in generating an artificial image corresponding to an input medical image of an object
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means.
  • an object structure measurement system comprising
  • a shape information obtaining means which obtains shape information for use in generating an artificial image corresponding to an input medical image of an object
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means
  • a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means.
  • the “artificial image corresponding to an input medical image of an object” means an artificial image representing an object substantially conforming in structure to the object of the input medical image.
  • the “structural physical value of an object” is a value of the anatomic structure of the object which can be geometrically measured or calculated, and may be, for instance, the angle, the length, the area, the roundness, the flatness or a value calculated by a combination of two or more these factors.
  • the “structural physical value of an object” may be a combination of the maximum widths of the thorax and the heart or the cardiothoracic ratio (the ratio of the width of the heart to the width of the thorax)
  • the “shape information obtaining means” may be a means which obtains shape information for use in generating an artificial image of the object by the use of a technology of AAM (Active Appearance Model).
  • AAM Active Appearance Model
  • the technology of AAM is a method of generating an artificial image using shape information and texture information, and, for example, generates an artificial image by warping a density pattern represented by the texture information into a shape represented by the shape information.
  • the shape information and the texture information can be represented by a linear sum of eigen vectors of the shape or the texture of an object obtained by carrying out statistical processing on a plurality of images each representing the structure of the object as samples. By controlling a parameter obtained by integrating the eigen vectors, substantially any shape or substantially any texture can be obtained.
  • the shape information of the object of the input medical image can be obtained by generating an artificial image representing an object substantially conforming in structure to the object of the input medical image. This can be accomplished by searching for a value of the parameter which minimizes the difference between corresponding pixels of the input medical image and the artificial image by changing the shape and the texture of the artificial image by controlling the parameter. (This technology will be described in detail later.)
  • the shape information clearly represents the shape of the object independently of clearness of the anatomic structure of the object represented by the texture information. Accordingly, when an artificial image substantially conforming to an image of the object including an unclear part is generated, anatomical characteristic positions of the object can be clearly detected on the basis of the shape information obtained at that time.
  • the shape information obtaining means may obtain a plurality of pieces of shape information for use in generating an artificial image corresponding to the input medical image of the object by a plurality of anatomical structural regions of the object.
  • anatomical structural region means, for instance, when the object is a chest, the thorax, the heart and the like, but may include a background together with such an anatomical structural region. Further, the “anatomical structural region” as used here means a part of the anatomical structural region such as the thorax, the heart or the like.
  • an anatomical characteristic position detecting system comprising
  • an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures an are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data,
  • an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means, and
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means.
  • an object structure measurement system comprising
  • an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data,
  • an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means
  • a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means.
  • an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, and an anatomical characteristic position of the object is detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means, unlike in those in accordance with the first and second aspects of the present invention, where shape information for use in generating an artificial image substantially conforming in structure to the object of an input medical image is obtained.
  • a computer program for causing a computer to perform the function of the anatomical characteristic position detecting system or the object structure measurement system in accordance with the present invention or each of the means of the systems may be recorded in a computer readable medium so that the computer can perform the function when loaded with the recording medium.
  • the computer readable medium is not limited to any specific type of storage devices and includes any kind of device, including but not limited to CDs, floppy disks, RAMs, ROMs, hard disks, magnetic tapes and internet downloads, in which computer instructions can be stored and/or transmitted. Transmission of the computer code through a network or through wireless transmission means is also within the scope of this invention.
  • computer code/instructions include, but are not limited to, source, object and executable code and can be in any language including higher level languages, assembly language and machine language.
  • the medical image need not be limited to a radiation image but may any image so long as it represents anatomical object structure.
  • the medical image may be an image obtained by a CT system or an MRI system.
  • image should be broadly interpreted to include image data as well as a visible image actually displayed.
  • an expression “to obtain an image” includes to merely obtain image data on the basis of which an image is displayed.
  • anatomical characteristic positions of an object can be detected even if the image has a part which is unclear to detect the anatomical object structure, since shape information for use in generating an artificial image corresponding to an input medical image of an object is obtained, and the anatomical characteristic positions of the object are detected on the basis of the shape information obtained.
  • the shape information when the shape information is obtained by a plurality of anatomical structural regions of the object, the shape information can be finely fitted to the input medical image and the accuracy of detecting the anatomical characteristic position of the object is improved.
  • anatomical characteristic positions of an object can be detected and a structural physical value of the object can be measured even if the image has a part which is unclear to detect the anatomical object structure, since shape information for use in generating an artificial image corresponding to an input medical image of an object is obtained, the anatomical characteristic positions of the object are detected on the basis of the shape information obtained and a structural physical value of the object is measured on the basis of the anatomical characteristic position of the object detected.
  • the shape information when the shape information is obtained by a plurality of anatomical structural regions of the object, the shape information can be finely fitted to the input medical image and the accuracy of measuring the structural physical value of the object is improved.
  • anatomical characteristic positions of an object can be detected even if the image has a part which is unclear to detect the anatomical object structure as in the anatomical characteristic position detecting system in accordance with the first aspect of the present invention, since an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, and the anatomical characteristic positions of the object are detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means.
  • anatomical characteristic positions of an object can be detected and a structural physical value of the object can be measured even if the image has a part which is unclear to detect the anatomical object structure as in the object structure measurement system in accordance with the second aspect of the present invention, since an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, the anatomical characteristic positions of the object are detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means and a structural physical value of the object is measured on the basis of the anatomical characteristic positions of the object detected.
  • FIG. 1 is a schematic view showing cardiothoracic ratio measuring system in accordance with a first embodiment of the present invention
  • FIG. 2 is a flow chart for illustrating the steps of generating an artificial image
  • FIGS. 3A to 3 C show landmarks attached to chest images on the outer edge of the lungs
  • FIG. 4 shows an image representing an average shape of chests
  • FIG. 5 shows warping of the landmarks attached on the outer edge of the lungs to the average object shape
  • FIGS. 6A to 6 C are views respectively showing different eigen shapes represented by landmarks
  • FIGS. 7A to 7 C are views respectively showing different eigen textures
  • FIG. 8 is a view showing the partial artificial images
  • FIG. 9 is a view showing a chest radiation image including a part whose anatomical structure is unclear
  • FIG. 10 is a view showing the landmarks and the profile of the cardiothorax obtained by joining the landmarks with lines, and
  • FIG. 11 is a schematic view showing cardiothoracic ratio measuring system in accordance with a second embodiment of the present invention.
  • a cardiothoracic ratio measuring system in accordance with a first embodiment of the present invention comprises a shape information obtaining means 10 which obtains shape information K′ by generating an artificial image P′ representing a chest substantially conforming in structure to an input chest radiation image P by the use of shape information K′ and texture information T′ obtained by changing a parameter obtained by integrating the eigen vectors of shape and texture, a cardiothorax profile detecting means 20 (an anatomical characteristic position detecting means) which detects a profile of the cardiothorax of the object on the basis of the shape information K′, and a cardiothoracic ratio calculating means (a measurement means) 30 which calculates the cardiothoracic ratio of the object on the basis of information d on the position of the profile of the cardiothorax detected by the cardiothorax profile detecting means 20 .
  • the shape information 10 obtains shape information K′ according to a technology of AAM to be described later by generating an artificial image P′ corresponding to the input chest radiation image P (representing a chest substantially conforming in structure to the input chest radiation image P) by the use of shape information K′ and texture information T′ obtained by changing a parameter obtained by integrating the eigen vectors of shape and texture.
  • AAM Active Appearance Model
  • images stored in a database in which a plurality of different images each representing an anatomical object structure are stored and anatomical characteristic shapes (landmarks) extracted from the images stored in the database are used as teacher data and applied to the AAM.
  • a statistical model of object structure (a model of change of shape, a model of change of texture and a model of correlation between the shape and the texture) is made on the basis of the teacher data, and by the use of the statistical model, substantially any artificial image (in shape and texture) corresponding to the input image can be obtained.
  • the images stored in the database be actually photographed in view of making a statistical model of the object structure, they may be artificial images made at random with the characteristic basic object structure held in the case where a sufficient number of images are not stored in the database. It is necessary to store also an image of an object having an abnormal structure so that an artificial image corresponding to an image of a diseased object such as those suffered from cardiomegaly can be generated.
  • step #1 an average shape and an average texture are generated.
  • FIGS. 3A to 3 C show landmarks attached to chest images S 1 , S 2 and S 3 on the outer edge of the lungs.
  • a plurality of black dots indicated at M in FIGS. 3A to 3 C are landmarks.
  • the points employed as the landmarks are anatomically equivalent points in the images shown in FIGS. 3A to 3 C.
  • Average object shape is made on the basis of the shapes of the objects represented by the landmarks.
  • FIG. 4 shows an image Sa representing an average shape of chests.
  • the parallel displacement components be normalized in advance. For example, in the case of a chest radiation image, the upper ends of the lungs are aligned with the horizontal center line
  • FIG. 5 shows warping of the landmarks attached on the outer edge of the lungs to the average object shape (indicated Ka in FIG. 5).
  • the amount of shift of each land mark is approximated by the use of the following two-dimensional quintic polynomial approximation.
  • x and y are coordinates of each landmark
  • x′ and y′ are the amounts of shift of each landmark to the average object shape to which the landmark is warped
  • n is the order
  • a ij and b ij are coefficients.
  • the coefficients in polynomial approximation are obtained by a least squares method.
  • the amount of shift of each pixel to the average object shape is calculated for each image and the pixel is warped.
  • the values of pixels which are warped to positions which are represented by not integers but figures including a decimal are obtained by first-order approximation on the basis of four pixels around the pixel. That is, the value of the pixel is distributed to the four pixels around the pixel in proportion to the distances from the pixel to the four pixels after the warping. Then average texture are obtained from the images warped to the average object shape.
  • step #2 a model of change of shape is generated.
  • An eigen shape is calculated by carrying out a principal component analysis on the shape of the object by the use of each shape of the images stored in the database and the average shape.
  • eigen vectors can be obtained by a method disclosed in [Matthew Turk, Alex Pentl and, Eigenfaces for recognition, Journal of Cognitive Neuroscience, vol. 3, num 1, 1991]. Any shape can be approximated as follows by the use of the linear sum of eigen shapes.
  • x is a shape factor (x 1 , y 1 , . . . , x i , y i , . . . , x n , y n ),
  • x ave is an average shape vector
  • P S is a natural vector of shape
  • b s is a set of shape factors.
  • FIGS. 6A to 6 C are views respectively showing different eigen shapes Ke 1 , Ke 2 and Ke 3 represented by landmarks M.
  • the solid line Ka shows the average shape.
  • step #3 a model of change of texture is generated.
  • An eigen texture is calculated by carrying out a principal component analysis on the texture of the object by the use of each texture of the images stored in the database and the average texture.
  • g is a factor of gray level (a pixel value of texture) (g 11 , g 21 , . . . , g ij , . . . , g nn ), g ave is a vector of the average gray level, P g is a natural vector of gray level, and b g is a set of gray level factors.
  • FIGS. 7A to 7 C are views respectively showing different eigen textures Te 1 , Te 2 and Te 3 .
  • step #4 a model of correlation between the shape and the texture is generated.
  • the shape/texture of each of the images stored in the database can be expressed by coefficients b s and b g which are values obtained by projection onto each eigen shape/eigen texture.
  • W s is a diagonal matrix and a weight on each shape coefficient.
  • step #5 an artificial image is generated.
  • An artificial image is generated for an input radiation image in the following manner.
  • An appearance parameter calculated on the basis of the average shape and images in the average shape is set as an initial value.
  • An artificial image obtained on the basis of the value of the appearance parameter obtained in step 2 is selected as an artificial image corresponding to the input radiation image.
  • the database is not essential to the object structure measurement system.
  • the artificial image P′ is divided into three anatomical structural regions, the right lung region, the left lung region and the heart region as shown in FIG. 8 and three partial artificial images P 1 ′, P 2 ′ and P 3 ′ corresponding to the respective regions are generated, whereby three pieces of shape information K 1 ′, K 2 ′ and K 3 ′ for the respective anatomical structural regions are obtained. That is, each of the images stored in the database is divided into the three anatomical structural regions, and the above models are generated by the regions.
  • ⁇ D w 1
  • P 1 , P 2 and P 3 are partial images corresponding to the respective anatomical regions of the input chest radiation image P
  • ⁇ (Pi ⁇ Pi′) is the total of the difference between the values of corresponding pixels of the partial image and the partial artificial image
  • w 1 , w 2 and w 3 are weight from 0 to 1.
  • the profile of the cardiothorax from which the widths of the thorax and the heart necessary to calculate the cardiothoracic ratio can be obtained by obtaining pieces of shape information K 1 ′, K 2 ′ and K 3 ′ for the respective anatomical structural regions. (Though the profile of the upper and lower portions of the heart is not accurately obtained, it involves no problem since the width is only measured.) Generally, by dividing an artificial image into a plurality of regions, the degree of freedom in adjusting the artificial image is increased and an artificial image more accurately corresponding to the input medical image can be generated.
  • the cardiothorax profile detecting means 20 detects the profile of the cardiothorax in this manner on the basis of the pieces of shape information K 1 ′, K 2 ′ and K 3 ′ for the respective anatomical structural regions of the object obtained upon generation of the partial artificial images P 1 ′, P 2 ′ and P 3 ′.
  • the shapes represented by the respective pieces of shape information K 1 ′, K 2 ′ and K 3 ′ are expressed by, for instance, landmarks.
  • the profile of the cardiothorax can be derived from the landmarks by joining the landmarks with lines, for instance, by smoothing or first-order approximation.
  • FIG. 10 shows the landmarks and the profile of the cardiothorax obtained by joining the landmarks with lines by the above method.
  • the cardiothorax ratio calculating means 30 detects the maximum width L 2 of the thorax and the maximum width L 1 of the heart and calculates the cardiothorax ratio (L 1 /L 2 ).
  • the profile of the object can be obtained even for an input image including an unclear portion by generating an artificial image corresponding to the input image and obtaining shape information, and accordingly, the structural physical values can be accurately measured.
  • the cardiothoracic ratio measuring system of this embodiment comprises an artificial image storage means 40 in which a plurality of pieces of artificial image data representing different object structures are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data, an artificial image selection means 50 which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image P out of the images represented by the plurality of pieces of image data stored in the artificial image storage means 40 , a cardiothorax profile detecting means 20 (an anatomical characteristic position detecting means) which detects a profile of the cardiothorax of the object on the basis of the shape information K′ attached to the artificial image data selected by the artificial image selection means 50 , and a cardiothoracic ratio calculating means (a measurement means) 30 which calculates the cardiothoracic ratio of the object on the basis of information d
  • an artificial image P′ representing a chest substantially conforming in structure to the chest represented by the input medical image P is selected out of the artificial images which have been generated and stored in the artificial image storage means 40 , and the cardiothoracic ratio of the object is calculated on the basis of the shape information K′ attached to the artificial image P′ selected by the artificial image selection means 50 , unlike in the cardiothoracic ratio measuring system of the first embodiment, where shape information is obtained by generating an artificial image for the input image P.
  • the artificial image P′ may be divided into anatomical structural regions, e.g., the right lung region, the left lung region and the heart region and partial artificial images corresponding to the respective regions may be selected.

Abstract

Shape information for use in generating an artificial image corresponding to an input medical image of an object is obtained, and an anatomical characteristic position of the object is detected on the basis of the shape information obtained.

Description

    BACKGROUND OF THE INVENTION
  • 1. Field of the Invention [0001]
  • This invention relates to an anatomical characteristic position detecting system and an object structure measurement system for aiding in diagnosing anatomical structure of an object on the basis of a medical image. [0002]
  • 2. Description of the Related Art [0003]
  • It has been carried out to measure a structural part such as a bone or an organ by the use of a medical image. For example, when diagnosing cardiomegaly, the widths of the thorax and the heart are measured by the use of a chest radiation image and the ratio of the latter to the former (cardiothoracic ratio) is calculated. [0004]
  • Conventionally, the widths of the thorax and the heart have been manually measured by applying a scale to the chest radiation image. Recently, an image is often digitized and handled as digital image data, and there is a demand for automatically calculating the cardiothoracic ratio on the basis of the digital image data. [0005]
  • In order to meet this demand, there have been proposed various technologies of automatically calculating the cardiothoracic ratio. To calculate the cardiothoracic ratio, it is necessary to identify the profiles of the thorax and the heart (cardiothorax). As the method of detecting the profile of the cardiothorax, there have been known a method disclosed in U.S. Pat. Nos. 5,790,690 and 6,011,862 where the profile of the thorax (a region including the left and right lungs), which is a closed region, is detected and a method proposed by Mr. Katsuragawa where the profile of the heart is detected by the use of a deformation ellipse model function. Further, there has been proposed in U.S. Laid-Open No. 20020057828 a method where the image of the chest is expressed in a polar coordinate system, the profile of the cardiothorax is automatically detected by template matching on the plane of the polar coordinate system by the use of a template which is substantially similar to the average cardiothorax profile (a reference profile), and the cardiothoracic ratio is calculated on the basis of the automatically detected profile of the cardiothorax. [0006]
  • However, it has been found that by the method disclosed in U.S. Laid-Open No. 20020057828, the error of the detected cardiothorax profile from the actual profile is too large to obtain an accuracy necessary to diagnosis. Accordingly, this applicant has proposed in Japanese Unexamined Patent Publication No.2003-6661 a more accurate method where the cardiothorax profile is divided into a plurality of regions and the regions are recombined after edges on borders of regions are detected. [0007]
  • However, the method disclosed in Japanese Unexamined Patent Publication No. 2003-6661 is disadvantageous in that edges on borders cannot be accurately found when fat is on the heart due to heart disease or edges are not clear due to noise. [0008]
  • SUMMARY OF THE INVENTION
  • In view of the foregoing observations and description, the primary object of the present invention is to provide an anatomical characteristic position detecting system and an object structure measurement system which make it feasible to detect anatomical characteristic positions and to measure structural physical values of an object on the basis of information on the detected characteristic positions even if the image has a part which is unclear to know the anatomical structure of the object. [0009]
  • In accordance with a first aspect of the present invention, there is provided an anatomical characteristic position detecting system comprising [0010]
  • a shape information obtaining means which obtains shape information for use in generating an artificial image corresponding to an input medical image of an object, and [0011]
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means. [0012]
  • In accordance with a second aspect of the present invention, there is provided an object structure measurement system comprising [0013]
  • a shape information obtaining means which obtains shape information for use in generating an artificial image corresponding to an input medical image of an object, [0014]
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means, and [0015]
  • a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means. [0016]
  • The “artificial image corresponding to an input medical image of an object” means an artificial image representing an object substantially conforming in structure to the object of the input medical image. [0017]
  • The “structural physical value of an object” is a value of the anatomic structure of the object which can be geometrically measured or calculated, and may be, for instance, the angle, the length, the area, the roundness, the flatness or a value calculated by a combination of two or more these factors. When the object is a chest, the “structural physical value of an object” may be a combination of the maximum widths of the thorax and the heart or the cardiothoracic ratio (the ratio of the width of the heart to the width of the thorax) [0018]
  • The “shape information obtaining means” may be a means which obtains shape information for use in generating an artificial image of the object by the use of a technology of AAM (Active Appearance Model). The technology of AAM is a method of generating an artificial image using shape information and texture information, and, for example, generates an artificial image by warping a density pattern represented by the texture information into a shape represented by the shape information. The shape information and the texture information can be represented by a linear sum of eigen vectors of the shape or the texture of an object obtained by carrying out statistical processing on a plurality of images each representing the structure of the object as samples. By controlling a parameter obtained by integrating the eigen vectors, substantially any shape or substantially any texture can be obtained. The shape information of the object of the input medical image can be obtained by generating an artificial image representing an object substantially conforming in structure to the object of the input medical image. This can be accomplished by searching for a value of the parameter which minimizes the difference between corresponding pixels of the input medical image and the artificial image by changing the shape and the texture of the artificial image by controlling the parameter. (This technology will be described in detail later.) [0019]
  • The shape information clearly represents the shape of the object independently of clearness of the anatomic structure of the object represented by the texture information. Accordingly, when an artificial image substantially conforming to an image of the object including an unclear part is generated, anatomical characteristic positions of the object can be clearly detected on the basis of the shape information obtained at that time. [0020]
  • In the anatomical characteristic position detecting system and the object structure measurement system in accordance with the first and second aspects of the present invention, the shape information obtaining means may obtain a plurality of pieces of shape information for use in generating an artificial image corresponding to the input medical image of the object by a plurality of anatomical structural regions of the object. [0021]
  • The “anatomical structural region” as used here means, for instance, when the object is a chest, the thorax, the heart and the like, but may include a background together with such an anatomical structural region. Further, the “anatomical structural region” as used here means a part of the anatomical structural region such as the thorax, the heart or the like. [0022]
  • In accordance with a third aspect of the present invention, there is provided an anatomical characteristic position detecting system comprising [0023]
  • an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures an are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data, [0024]
  • an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means, and [0025]
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means. [0026]
  • In accordance with a fourth aspect of the present invention, there is provided an object structure measurement system comprising [0027]
  • an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data, [0028]
  • an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means, [0029]
  • an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means, and [0030]
  • a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means. [0031]
  • That is, in the anatomical characteristic position detecting system and the object structure measurement system in accordance with the third and fourth aspects of the present invention, an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, and an anatomical characteristic position of the object is detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means, unlike in those in accordance with the first and second aspects of the present invention, where shape information for use in generating an artificial image substantially conforming in structure to the object of an input medical image is obtained. [0032]
  • A computer program for causing a computer to perform the function of the anatomical characteristic position detecting system or the object structure measurement system in accordance with the present invention or each of the means of the systems may be recorded in a computer readable medium so that the computer can perform the function when loaded with the recording medium. A skilled artisan would know that the computer readable medium is not limited to any specific type of storage devices and includes any kind of device, including but not limited to CDs, floppy disks, RAMs, ROMs, hard disks, magnetic tapes and internet downloads, in which computer instructions can be stored and/or transmitted. Transmission of the computer code through a network or through wireless transmission means is also within the scope of this invention. Additionally, computer code/instructions include, but are not limited to, source, object and executable code and can be in any language including higher level languages, assembly language and machine language. [0033]
  • The medical image need not be limited to a radiation image but may any image so long as it represents anatomical object structure. For example, the medical image may be an image obtained by a CT system or an MRI system. [0034]
  • The “image” as used here should be broadly interpreted to include image data as well as a visible image actually displayed. For example, an expression “to obtain an image” includes to merely obtain image data on the basis of which an image is displayed. [0035]
  • In the anatomical characteristic position detecting system in accordance with the first aspect of the present invention, anatomical characteristic positions of an object can be detected even if the image has a part which is unclear to detect the anatomical object structure, since shape information for use in generating an artificial image corresponding to an input medical image of an object is obtained, and the anatomical characteristic positions of the object are detected on the basis of the shape information obtained. [0036]
  • In the anatomical characteristic position detecting system in accordance with the first aspect of the present invention, when the shape information is obtained by a plurality of anatomical structural regions of the object, the shape information can be finely fitted to the input medical image and the accuracy of detecting the anatomical characteristic position of the object is improved. [0037]
  • In the object structure measurement system in accordance with the second aspect of the present invention, anatomical characteristic positions of an object can be detected and a structural physical value of the object can be measured even if the image has a part which is unclear to detect the anatomical object structure, since shape information for use in generating an artificial image corresponding to an input medical image of an object is obtained, the anatomical characteristic positions of the object are detected on the basis of the shape information obtained and a structural physical value of the object is measured on the basis of the anatomical characteristic position of the object detected. [0038]
  • In the object structure measurement system in accordance with the second aspect of the present invention, when the shape information is obtained by a plurality of anatomical structural regions of the object, the shape information can be finely fitted to the input medical image and the accuracy of measuring the structural physical value of the object is improved. [0039]
  • In the anatomical characteristic position detecting system in accordance with the third aspect of the present invention, anatomical characteristic positions of an object can be detected even if the image has a part which is unclear to detect the anatomical object structure as in the anatomical characteristic position detecting system in accordance with the first aspect of the present invention, since an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, and the anatomical characteristic positions of the object are detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means. [0040]
  • In the object structure measurement system in accordance with the fourth aspect of the present invention, anatomical characteristic positions of an object can be detected and a structural physical value of the object can be measured even if the image has a part which is unclear to detect the anatomical object structure as in the object structure measurement system in accordance with the second aspect of the present invention, since an artificial image representing an object substantially conforming in structure to the object of an input medical image is selected out of the images represented by the plurality of pieces of image data which have been generated and stored in the artificial image storage means, the anatomical characteristic positions of the object are detected on the basis of the shape information attached to the artificial image selected by the artificial image selection means and a structural physical value of the object is measured on the basis of the anatomical characteristic positions of the object detected.[0041]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a schematic view showing cardiothoracic ratio measuring system in accordance with a first embodiment of the present invention, [0042]
  • FIG. 2 is a flow chart for illustrating the steps of generating an artificial image, [0043]
  • FIGS. 3A to [0044] 3C show landmarks attached to chest images on the outer edge of the lungs,
  • FIG. 4 shows an image representing an average shape of chests, [0045]
  • FIG. 5 shows warping of the landmarks attached on the outer edge of the lungs to the average object shape, [0046]
  • FIGS. 6A to [0047] 6C are views respectively showing different eigen shapes represented by landmarks,
  • FIGS. 7A to [0048] 7C are views respectively showing different eigen textures,
  • FIG. 8 is a view showing the partial artificial images, [0049]
  • FIG. 9 is a view showing a chest radiation image including a part whose anatomical structure is unclear, [0050]
  • FIG. 10 is a view showing the landmarks and the profile of the cardiothorax obtained by joining the landmarks with lines, and [0051]
  • FIG. 11 is a schematic view showing cardiothoracic ratio measuring system in accordance with a second embodiment of the present invention.[0052]
  • DESCRIPTION OF THE PREFERRED EMBODIMENTS
  • In FIG. 1, a cardiothoracic ratio measuring system in accordance with a first embodiment of the present invention comprises a shape [0053] information obtaining means 10 which obtains shape information K′ by generating an artificial image P′ representing a chest substantially conforming in structure to an input chest radiation image P by the use of shape information K′ and texture information T′ obtained by changing a parameter obtained by integrating the eigen vectors of shape and texture, a cardiothorax profile detecting means 20 (an anatomical characteristic position detecting means) which detects a profile of the cardiothorax of the object on the basis of the shape information K′, and a cardiothoracic ratio calculating means (a measurement means) 30 which calculates the cardiothoracic ratio of the object on the basis of information d on the position of the profile of the cardiothorax detected by the cardiothorax profile detecting means 20.
  • Operation of the cardiothoracic ratio measuring system of this embodiment will be described, hereinbelow. [0054]
  • When a chest radiation image P including a part U whose anatomical structure is unclear as shown in FIG. 9 is input, the [0055] shape information 10 obtains shape information K′ according to a technology of AAM to be described later by generating an artificial image P′ corresponding to the input chest radiation image P (representing a chest substantially conforming in structure to the input chest radiation image P) by the use of shape information K′ and texture information T′ obtained by changing a parameter obtained by integrating the eigen vectors of shape and texture.
  • The technology of AAM according to which an artificial image is generated will be described, hereinbelow. [0056]
  • The technology of AAM has been proposed by Cootes et al of the university of Manchester and is disclosed in “Active Appearance Model (AAM)”[T. F. Cootes, and C. J. Taylor, Active Appearance Models, Proc, European Conference on Computer Vision, Vol.2, pp.484-498, Springer, 1998]. [0057]
  • In this embodiment, images stored in a database in which a plurality of different images each representing an anatomical object structure are stored and anatomical characteristic shapes (landmarks) extracted from the images stored in the database are used as teacher data and applied to the AAM. A statistical model of object structure (a model of change of shape, a model of change of texture and a model of correlation between the shape and the texture) is made on the basis of the teacher data, and by the use of the statistical model, substantially any artificial image (in shape and texture) corresponding to the input image can be obtained. [0058]
  • Though it is preferred that the images stored in the database be actually photographed in view of making a statistical model of the object structure, they may be artificial images made at random with the characteristic basic object structure held in the case where a sufficient number of images are not stored in the database. It is necessary to store also an image of an object having an abnormal structure so that an artificial image corresponding to an image of a diseased object such as those suffered from cardiomegaly can be generated. [0059]
  • The steps of generating an artificial image corresponding to an input image will be described on the case where a chest radiation image is input, by way of example, hereinbelow. [0060]
  • In [0061] step #1, an average shape and an average texture are generated.
  • Anatomical characteristic shapes of the objects are first extracted from the images stored in the database as land marks. FIGS. 3A to [0062] 3C show landmarks attached to chest images S1, S2 and S3 on the outer edge of the lungs. A plurality of black dots indicated at M in FIGS. 3A to 3C are landmarks. The points employed as the landmarks are anatomically equivalent points in the images shown in FIGS. 3A to 3C. Average object shape is made on the basis of the shapes of the objects represented by the landmarks. FIG. 4 shows an image Sa representing an average shape of chests. When the shape of the object is extracted in the form of landmarks, it is preferred that the parallel displacement components be normalized in advance. For example, in the case of a chest radiation image, the upper ends of the lungs are aligned with the horizontal center line
  • Then in order to warp the images in the database to the average object shape thus made, the amount of shift of each landmark on each image to that of the average object shape is calculated. FIG. 5 shows warping of the landmarks attached on the outer edge of the lungs to the average object shape (indicated Ka in FIG. 5). The amount of shift of each land mark is approximated by the use of the following two-dimensional quintic polynomial approximation. [0063]
  • x′=x+Δx  (1)
  • y′=y+Δy  (2) Δ x = i = 0 n j = 0 n - i a ij · x i · y i ( 3 ) Δ y = i = 0 n j = 0 n - i b ij · x i · y i ( 4 )
    Figure US20030228040A1-20031211-M00001
  • wherein x and y are coordinates of each landmark, x′ and y′ are the amounts of shift of each landmark to the average object shape to which the landmark is warped, n is the order, and a[0064] ij and bij are coefficients. The coefficients in polynomial approximation are obtained by a least squares method.
  • By the use of the above formulae, the amount of shift of each pixel to the average object shape is calculated for each image and the pixel is warped. The values of pixels which are warped to positions which are represented by not integers but figures including a decimal are obtained by first-order approximation on the basis of four pixels around the pixel. That is, the value of the pixel is distributed to the four pixels around the pixel in proportion to the distances from the pixel to the four pixels after the warping. Then average texture are obtained from the images warped to the average object shape. [0065]
  • At this time point, the average shape of the images stored in the database, images warped to the average shape and the average texture are obtained. [0066]
  • In step #2, a model of change of shape is generated. [0067]
  • An eigen shape is calculated by carrying out a principal component analysis on the shape of the object by the use of each shape of the images stored in the database and the average shape. For example, eigen vectors can be obtained by a method disclosed in [Matthew Turk, Alex Pentl and, Eigenfaces for recognition, Journal of Cognitive Neuroscience, vol. 3, [0068] num 1, 1991]. Any shape can be approximated as follows by the use of the linear sum of eigen shapes.
  • x=x ave +P s ·b s  (5)
  • wherein x is a shape factor (x[0069] 1, y1, . . . , xi, yi, . . . , xn, yn), xave is an average shape vector, PS is a natural vector of shape, and bs is a set of shape factors. FIGS. 6A to 6C are views respectively showing different eigen shapes Ke1, Ke2 and Ke3 represented by landmarks M. In FIGS. 6A to 6C, the solid line Ka shows the average shape.
  • In step #3, a model of change of texture is generated. [0070]
  • An eigen texture is calculated by carrying out a principal component analysis on the texture of the object by the use of each texture of the images stored in the database and the average texture. [0071]
  • g=g ave +P g ·b g  (6)
  • wherein g is a factor of gray level (a pixel value of texture) (g[0072] 11, g21, . . . , gij, . . . , gnn), gave is a vector of the average gray level, Pg is a natural vector of gray level, and bg is a set of gray level factors. FIGS. 7A to 7C are views respectively showing different eigen textures Te1, Te2 and Te3.
  • In step #4, a model of correlation between the shape and the texture is generated. The shape/texture of each of the images stored in the database can be expressed by coefficients b[0073] s and bg which are values obtained by projection onto each eigen shape/eigen texture. An appearance parameter for controlling the shape and the gray level is calculated by carrying out a principal component analysis on the coefficients as follows. b = ( W s · b s b g ) = ( W s · P s T · ( x - x ave ) P g T · ( g - g ave ) ) ( 7 )
    Figure US20030228040A1-20031211-M00002
  • wherein W[0074] s is a diagonal matrix and a weight on each shape coefficient.
  • b=Q·c  (8)
  • wherein Q is an eigenvector and c is an appearance parameter for controlling the shape and the gray level (c is 0 in average.) By varying the appearance parameter, substantially any artificial image of the object can be obtained from the images stored in the database. [0075]
  • In step #5, an artificial image is generated. An artificial image is generated for an input radiation image in the following manner. [0076]
  • 1. An appearance parameter calculated on the basis of the average shape and images in the average shape is set as an initial value. [0077]
  • 2. Artificial images are generated while changing the parameter from an appearance parameter for an uppermost eigen vector and an appearance parameter which minimizes the difference between the input radiation image and the generated artificial image is searched for. [0078]
  • 3. An artificial image obtained on the basis of the value of the appearance parameter obtained in step 2 is selected as an artificial image corresponding to the input radiation image. [0079]
  • In accordance with the above method, when elements such as the average shape/texture, the eigen shape/eigen texture and the eigenvector Q of the object, have been obtained in advance, substantially any artificial image of the object can be generated by simply controlling the appearance parameter c. [0080]
  • In this embodiment, it is assumed that the above elements have been obtained for the chest radiation image (a medical image of the object) and an artificial image P′ corresponding to the input chest radiation P can be generated and the shape information can be obtained by simply controlling the appearance parameter c. Accordingly, the database is not essential to the object structure measurement system. [0081]
  • In this embodiment, the artificial image P′ is divided into three anatomical structural regions, the right lung region, the left lung region and the heart region as shown in FIG. 8 and three partial artificial images P[0082] 1′, P2′ and P3′ corresponding to the respective regions are generated, whereby three pieces of shape information K1′, K2′ and K3′ for the respective anatomical structural regions are obtained. That is, each of the images stored in the database is divided into the three anatomical structural regions, and the above models are generated by the regions. Then the values of appearance parameters c1, c2 and c3 defined by the following formulae (9) to (11) which minimize the difference A D defined by the following formula (12) are obtained and the partial artificial images P1′, P2′ and P3′ are generated.
  • b1=Q1·c1  (9)
  • b2=Q2·c2  (10)
  • b3=Q3·c3  (11)
  • ΔD=w1|Δ(P1−P1′|+w2|Δ(P2−P2′|+w3|Δ(P3−P3′|  (12)
  • wherein P[0083] 1, P2 and P3 are partial images corresponding to the respective anatomical regions of the input chest radiation image P, Δ (Pi−Pi′) is the total of the difference between the values of corresponding pixels of the partial image and the partial artificial image, and w1, w2 and w3 are weight from 0 to 1.
  • The profile of the cardiothorax from which the widths of the thorax and the heart necessary to calculate the cardiothoracic ratio can be obtained by obtaining pieces of shape information K[0084] 1′, K2′ and K3′ for the respective anatomical structural regions. (Though the profile of the upper and lower portions of the heart is not accurately obtained, it involves no problem since the width is only measured.) Generally, by dividing an artificial image into a plurality of regions, the degree of freedom in adjusting the artificial image is increased and an artificial image more accurately corresponding to the input medical image can be generated.
  • The cardiothorax [0085] profile detecting means 20 detects the profile of the cardiothorax in this manner on the basis of the pieces of shape information K1′, K2′ and K3′ for the respective anatomical structural regions of the object obtained upon generation of the partial artificial images P1′, P2′ and P3′. The shapes represented by the respective pieces of shape information K1′, K2′ and K3′ are expressed by, for instance, landmarks. The profile of the cardiothorax can be derived from the landmarks by joining the landmarks with lines, for instance, by smoothing or first-order approximation. FIG. 10 shows the landmarks and the profile of the cardiothorax obtained by joining the landmarks with lines by the above method.
  • When detected position information d is input into the cardiothorax ratio calculating means [0086] 30, the cardiothorax ratio calculating means 30 detects the maximum width L2 of the thorax and the maximum width L1 of the heart and calculates the cardiothorax ratio (L1/L2).
  • In accordance with the conventional technology where a profile of the object is extracted from a medical image, when the input image includes an unclear portion, the profile of the object cannot be detected to cause errors and even if the profile of the object can be detected, the calculated structural physical values contains a large error. To the contrast, in the cardiothoracic ratio measuring system of this embodiment, the profile of the object can be obtained even for an input image including an unclear portion by generating an artificial image corresponding to the input image and obtaining shape information, and accordingly, the structural physical values can be accurately measured. [0087]
  • A cardiothoracic ratio measuring system in accordance with a second embodiment of the present invention will be described with reference to FIG. 11, hereinbelow. As shown in FIG. 11, the cardiothoracic ratio measuring system of this embodiment comprises an artificial image storage means [0088] 40 in which a plurality of pieces of artificial image data representing different object structures are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data, an artificial image selection means 50 which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image P out of the images represented by the plurality of pieces of image data stored in the artificial image storage means 40, a cardiothorax profile detecting means 20 (an anatomical characteristic position detecting means) which detects a profile of the cardiothorax of the object on the basis of the shape information K′ attached to the artificial image data selected by the artificial image selection means 50, and a cardiothoracic ratio calculating means (a measurement means) 30 which calculates the cardiothoracic ratio of the object on the basis of information d on the position of the profile of the cardiothorax detected by the cardiothorax profile detecting means 20.
  • In the cardiothoracic ratio measuring system of this embodiment, an artificial image P′ representing a chest substantially conforming in structure to the chest represented by the input medical image P is selected out of the artificial images which have been generated and stored in the artificial image storage means [0089] 40, and the cardiothoracic ratio of the object is calculated on the basis of the shape information K′ attached to the artificial image P′ selected by the artificial image selection means 50, unlike in the cardiothoracic ratio measuring system of the first embodiment, where shape information is obtained by generating an artificial image for the input image P.
  • Also in this embodiment, the artificial image P′ may be divided into anatomical structural regions, e.g., the right lung region, the left lung region and the heart region and partial artificial images corresponding to the respective regions may be selected. [0090]

Claims (7)

What is claimed is:
1. An anatomical characteristic position detecting system comprising
a shape information obtaining means which obtains shape information for use at least in generating an artificial image corresponding to an input medical image of an object, and
an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means.
2. An anatomical characteristic position detecting system as defined in claim 1 in which the shape information obtaining means obtains a plurality of pieces of shape information for use in generating an artificial image corresponding to the input medical image of the object by a plurality of anatomical structural regions of the object.
3. An object structure measurement system comprising
a shape information obtaining means which obtains shape information for use at least in generating an artificial image corresponding to an input medical image of an object,
an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information obtained by the shape information obtaining means, and
a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means.
4. An object structure measurement system as defined in claim 3 in which the shape information obtaining means obtains a plurality of pieces of shape information for use in generating an artificial image corresponding to the input medical image of the object by a plurality of anatomical structural regions of the object.
5. An object structure measurement system as defined in claim 3 in which the structural physical value is a cardiothoracic ratio.
6. An anatomical characteristic position detecting system comprising
an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures an are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data,
an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means, and
an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means.
7. An object structure measurement system comprising
an artificial image storage means in which a plurality of pieces of artificial image data representing different object structures are stored with each piece of artificial image data attached with shape information on the shape of the object represented by the artificial image data,
an artificial image selection means which selects an artificial image representing an object substantially conforming in structure to the object of an input medical image out of the images represented by the plurality of pieces of image data stored in the artificial image storage means,
an anatomical characteristic position detecting means which detects an anatomical characteristic position of the object on the basis of the shape information attached to the artificial image selected by the artificial image selection means, and
a measurement means which measures a structural physical value of the object on the basis of the anatomical characteristic position of the object detected by the anatomical characteristic position detecting means.
US10/455,417 2002-06-06 2003-06-06 Anatomical characteristic position detecting system, and object structure measurement system Expired - Fee Related US7315635B2 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2002165104A JP4138371B2 (en) 2002-06-06 2002-06-06 Anatomical feature position detecting device, recording medium, subject structure measuring device, and recording medium
JP(PAT)165104/2002 2002-06-06

Publications (2)

Publication Number Publication Date
US20030228040A1 true US20030228040A1 (en) 2003-12-11
US7315635B2 US7315635B2 (en) 2008-01-01

Family

ID=29706678

Family Applications (1)

Application Number Title Priority Date Filing Date
US10/455,417 Expired - Fee Related US7315635B2 (en) 2002-06-06 2003-06-06 Anatomical characteristic position detecting system, and object structure measurement system

Country Status (2)

Country Link
US (1) US7315635B2 (en)
JP (1) JP4138371B2 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060171586A1 (en) * 2004-11-08 2006-08-03 Bogdan Georgescu Method of database-guided segmentation of anatomical structures having complex appearances
US20070263769A1 (en) * 2006-05-10 2007-11-15 Stefan Roell Device, method and computer program product for generation of an exposure plan
US20080101667A1 (en) * 2006-10-25 2008-05-01 Rcadia Medical Imaging Ltd. Method and system for the presentation of blood vessel structures and identified pathologies
US20080101674A1 (en) * 2006-10-25 2008-05-01 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures and pathologies
US7873194B2 (en) * 2006-10-25 2011-01-18 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures and pathologies in support of a triple rule-out procedure
US7940977B2 (en) 2006-10-25 2011-05-10 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures to identify calcium or soft plaque pathologies
US7940970B2 (en) 2006-10-25 2011-05-10 Rcadia Medical Imaging, Ltd Method and system for automatic quality control used in computerized analysis of CT angiography
JP2014018376A (en) * 2012-07-17 2014-02-03 Canon Inc Image processing device and method
CN106991697A (en) * 2017-04-17 2017-07-28 心医国际数字医疗系统(大连)有限公司 The ambition ratio measuring method of chest digitized video
US20180108156A1 (en) * 2016-10-17 2018-04-19 Canon Kabushiki Kaisha Radiographing apparatus, radiographing system, radiographing method, and storage medium
US10922894B2 (en) * 2016-06-06 2021-02-16 Biodigital, Inc. Methodology and system for mapping a virtual human body

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8276091B2 (en) * 2003-09-16 2012-09-25 Ram Consulting Haptic response system and method of use
JP4795658B2 (en) * 2004-03-23 2011-10-19 富士フイルム株式会社 Difference image acquisition method, difference image acquisition device, and program thereof
JP4690204B2 (en) * 2006-01-16 2011-06-01 富士フイルム株式会社 Image reproduction apparatus and program thereof
JP4891144B2 (en) * 2007-05-08 2012-03-07 キヤノン株式会社 Image search apparatus and image search method
CN101502420B (en) * 2008-02-05 2013-05-22 Ge医疗系统环球技术有限公司 Image processing method and equipment
US8077956B2 (en) * 2008-10-28 2011-12-13 Carestream Health, Inc. Orientation detection for chest radiographic images
JP2010220908A (en) 2009-03-25 2010-10-07 Sony Corp Image processor, image processing method, and program
JP5405285B2 (en) * 2009-12-15 2014-02-05 株式会社東芝 X-ray equipment
US10565723B2 (en) 2016-02-02 2020-02-18 Samsung Electronics Co., Ltd. Systems and methods for drift correction
US10579162B2 (en) 2016-03-24 2020-03-03 Samsung Electronics Co., Ltd. Systems and methods to correct a vehicle induced change of direction

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5790690A (en) * 1995-04-25 1998-08-04 Arch Development Corporation Computer-aided method for automated image feature analysis and diagnosis of medical images
US5862249A (en) * 1995-12-01 1999-01-19 Eastman Kodak Company Automated method and system for determination of positional orientation of digital radiographic images
US5943435A (en) * 1997-10-07 1999-08-24 Eastman Kodak Company Body part recognition in radiographic images
US20020057828A1 (en) * 2000-11-06 2002-05-16 Fuji Photo Film Co., Ltd. Apparatus for automatically setting measurement reference element and measuring geometric feature of image
US20020085743A1 (en) * 2000-04-04 2002-07-04 Konica Corporation Image processing selecting method, image selecting method and image processing apparatus
US20030210813A1 (en) * 2002-05-13 2003-11-13 Fuji Photo Film Co., Ltd. Method and apparatus for forming images and image furnishing service system
US6678703B2 (en) * 2000-06-22 2004-01-13 Radvault, Inc. Medical image management system and method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002109550A (en) * 2000-09-29 2002-04-12 Fuji Photo Film Co Ltd Cardiothoracic contour detection method and cardiothoracic ratio calculation method
JP3679990B2 (en) * 2000-10-31 2005-08-03 株式会社東芝 Medical image processing apparatus and method
JP2003006661A (en) 2001-06-22 2003-01-10 Fuji Photo Film Co Ltd Thoracic contour detector

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5790690A (en) * 1995-04-25 1998-08-04 Arch Development Corporation Computer-aided method for automated image feature analysis and diagnosis of medical images
US6011862A (en) * 1995-04-25 2000-01-04 Arch Development Corporation Computer-aided method for automated image feature analysis and diagnosis of digitized medical images
US5862249A (en) * 1995-12-01 1999-01-19 Eastman Kodak Company Automated method and system for determination of positional orientation of digital radiographic images
US5943435A (en) * 1997-10-07 1999-08-24 Eastman Kodak Company Body part recognition in radiographic images
US20020085743A1 (en) * 2000-04-04 2002-07-04 Konica Corporation Image processing selecting method, image selecting method and image processing apparatus
US6678703B2 (en) * 2000-06-22 2004-01-13 Radvault, Inc. Medical image management system and method
US20020057828A1 (en) * 2000-11-06 2002-05-16 Fuji Photo Film Co., Ltd. Apparatus for automatically setting measurement reference element and measuring geometric feature of image
US20030210813A1 (en) * 2002-05-13 2003-11-13 Fuji Photo Film Co., Ltd. Method and apparatus for forming images and image furnishing service system

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7876934B2 (en) * 2004-11-08 2011-01-25 Siemens Medical Solutions Usa, Inc. Method of database-guided segmentation of anatomical structures having complex appearances
US20060171586A1 (en) * 2004-11-08 2006-08-03 Bogdan Georgescu Method of database-guided segmentation of anatomical structures having complex appearances
US20070263769A1 (en) * 2006-05-10 2007-11-15 Stefan Roell Device, method and computer program product for generation of an exposure plan
US8238636B2 (en) 2006-05-10 2012-08-07 Siemens Aktiengesellschaft Device, method and computer-readable storage medium for generating an exposure plan to implement radiation therapy
US20110222661A1 (en) * 2006-05-10 2011-09-15 Stefan Roell Device, method and computer program product for generation of an expousre plan
US7995826B2 (en) * 2006-05-10 2011-08-09 Siemens Aktiengesellschaft Device, method and computer program product for generation of an exposure plan
US7940970B2 (en) 2006-10-25 2011-05-10 Rcadia Medical Imaging, Ltd Method and system for automatic quality control used in computerized analysis of CT angiography
US7940977B2 (en) 2006-10-25 2011-05-10 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures to identify calcium or soft plaque pathologies
US7873194B2 (en) * 2006-10-25 2011-01-18 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures and pathologies in support of a triple rule-out procedure
US7983459B2 (en) 2006-10-25 2011-07-19 Rcadia Medical Imaging Ltd. Creating a blood vessel tree from imaging data
US7860283B2 (en) * 2006-10-25 2010-12-28 Rcadia Medical Imaging Ltd. Method and system for the presentation of blood vessel structures and identified pathologies
US20080101674A1 (en) * 2006-10-25 2008-05-01 Rcadia Medical Imaging Ltd. Method and system for automatic analysis of blood vessel structures and pathologies
US8103074B2 (en) 2006-10-25 2012-01-24 Rcadia Medical Imaging Ltd. Identifying aorta exit points from imaging data
US20080101667A1 (en) * 2006-10-25 2008-05-01 Rcadia Medical Imaging Ltd. Method and system for the presentation of blood vessel structures and identified pathologies
JP2014018376A (en) * 2012-07-17 2014-02-03 Canon Inc Image processing device and method
US10922894B2 (en) * 2016-06-06 2021-02-16 Biodigital, Inc. Methodology and system for mapping a virtual human body
US20180108156A1 (en) * 2016-10-17 2018-04-19 Canon Kabushiki Kaisha Radiographing apparatus, radiographing system, radiographing method, and storage medium
US10861197B2 (en) * 2016-10-17 2020-12-08 Canon Kabushiki Kaisha Radiographing apparatus, radiographing system, radiographing method, and storage medium
CN106991697A (en) * 2017-04-17 2017-07-28 心医国际数字医疗系统(大连)有限公司 The ambition ratio measuring method of chest digitized video

Also Published As

Publication number Publication date
JP2004008419A (en) 2004-01-15
US7315635B2 (en) 2008-01-01
JP4138371B2 (en) 2008-08-27

Similar Documents

Publication Publication Date Title
US7315635B2 (en) Anatomical characteristic position detecting system, and object structure measurement system
US9524582B2 (en) Method and system for constructing personalized avatars using a parameterized deformable mesh
US8345927B2 (en) Registration processing apparatus, registration method, and storage medium
US8023708B2 (en) Method of segmenting anatomic entities in digital medical images
Bosch et al. Automatic segmentation of echocardiographic sequences by active appearance motion models
US8682054B2 (en) Method and system for propagation of myocardial infarction from delayed enhanced cardiac imaging to cine magnetic resonance imaging using hybrid image registration
US7616789B2 (en) Method and apparatus for forming images and image furnishing service system
US7382907B2 (en) Segmenting occluded anatomical structures in medical images
US20010048757A1 (en) Method and apparatus for matching positions of images
US8165359B2 (en) Method of constructing gray value or geometric models of anatomic entity in medical image
US7876945B2 (en) Method for processing slice images
US7221786B2 (en) Method for automatic construction of 2D statistical shape model for the lung regions
CN100577103C (en) Image processing device and method which use two images
US20220207771A1 (en) Heart Position Estimation
Ecabert et al. Modeling shape variability for full heart segmentation in cardiac computed-tomography images
Hansegard et al. Constrained active appearance models for segmentation of triplane echocardiograms
US8433159B1 (en) Compressed target movement model using interpolation
Orderud et al. Real-time tracking of the left ventricle in 3D echocardiography using a state estimation approach
US7697739B2 (en) Method, apparatus and program for image processing, and abnormal shadow detection
US20060149511A1 (en) Method of segmenting a three-dimensional data set allowing user corrections
US8867809B2 (en) Image processing method
CN111127488A (en) Method for automatically constructing patient anatomical structure model based on statistical shape model
Foroughi et al. Elastic registration of 3d ultrasound images
JP2002109550A (en) Cardiothoracic contour detection method and cardiothoracic ratio calculation method
Ecabert et al. Automatic whole heart segmentation in CT images: method and validation

Legal Events

Date Code Title Description
AS Assignment

Owner name: FUJI PHOTO FILM CO., LTD., JAPAN

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:OOSAWA, AKIRA;REEL/FRAME:014143/0969

Effective date: 20030519

AS Assignment

Owner name: FUJIFILM CORPORATION, JAPAN

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:FUJIFILM HOLDINGS CORPORATION (FORMERLY FUJI PHOTO FILM CO., LTD.);REEL/FRAME:018904/0001

Effective date: 20070130

Owner name: FUJIFILM CORPORATION,JAPAN

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:FUJIFILM HOLDINGS CORPORATION (FORMERLY FUJI PHOTO FILM CO., LTD.);REEL/FRAME:018904/0001

Effective date: 20070130

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

FPAY Fee payment

Year of fee payment: 4

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
STCH Information on status: patent discontinuation

Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362

FP Lapsed due to failure to pay maintenance fee

Effective date: 20160101