US3869602A - Apparatus for automatic computation of cardiothoracic ratio - Google Patents
Apparatus for automatic computation of cardiothoracic ratio Download PDFInfo
- Publication number
- US3869602A US3869602A US284558A US28455872A US3869602A US 3869602 A US3869602 A US 3869602A US 284558 A US284558 A US 284558A US 28455872 A US28455872 A US 28455872A US 3869602 A US3869602 A US 3869602A
- Authority
- US
- United States
- Prior art keywords
- cardiothoracic
- data
- heart
- data processing
- automatically measuring
- 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.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/59—Transmissivity
- G01N21/5907—Densitometers
- G01N21/5911—Densitometers of the scanning type
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
-
- 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/503—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 heart
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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/30048—Heart; Cardiac
Definitions
- G06f /42 and heart are first inferred by processing information [58] Field of Search 235/1513; 444/1; on brightness from an X-ray image and h n n rr r 128/205 R, 2.06 R, 2.06 B, 2.6 G, 2.1 R, in the inference of the profiles is corrected by contex- 419 R, DIG, 3; 178/DIG, 5, DIG, 22; tual data. processing, so that the cardiothoracic ratio 250/320, 321, 322, 323, 401 may be computed on the basis of the corrected data.
- the present invention relates to an apparatus for automatically computing the cardiothoracic ratio especially on the basis of an X-ray photograph.
- a heart disease is accompanied by an abnormal enlargement of the heart
- the universally adopted method of measurement thereof is by the cardiothoracic ratio advocated by C. S. Danzer, which is defined as a ratio in percentage of the transverse diameter of the heart to the maximal transverse internal diameter of the thorax.
- C. S. Danzer which is defined as a ratio in percentage of the transverse diameter of the heart to the maximal transverse internal diameter of the thorax.
- the CTR is usually determined with a rule or similar means applied to an X-ray or equivalent image of the chest. Even though various devices have been devised in an effort to facilitate the measurement, no substantial progress has been made in the direction of greatly saving the labor for the measurement.
- the principal operations or operating principle of the apparatus according to the present invention consist of the first step of processing brightness information included in an X-ray image by an improved method to infer the profiles of the thorax and heart, the second step of correcting an error, if any, in the inference by contextual data processing and the third step of computing the cardiothoracic ratio on the basis of the corrected data.
- Another object of the present invention is to provide an automatic cardiothoracic ratio measuring apparatus with a high performance at a reasonable cost.
- FIG. 1 is a diagram for explaining what is the cardiothoracic ratio
- FIG. 2 is a diagram showing an electrical signal generated by photo-electric conversion of an X-ray chest image along a straight line at right angles to the median line;
- Flg. 3 is a diagram showing a typical electrical signal obtained from an actual X-ray chest image
- FIG. 4 is a diagram for explaning the manner in which the peripheral points are determined in a typical electrical signal
- FIG. 5 is a diagram similar to FIG. 4 for explaining the manner in which the peripheral points are determined in a distorted electrical signal
- FIG. 6 is a diagram for explaining the manner in which the transverse diameter of the heart and the maximal transverse internal diameter of the thorax are determined; 1
- FIGS. 7A and 7B are diagrams for explaining the manner in which data obtained are corrected
- FIG. 8 is a block diagram showing an embodiment of the present invention.
- FIGS. 9 and 10 are diagrams showing in detail a part of the embodiment of FIG. 8;
- FIG. 11 shows a scanning locus
- FIG. 13 is a plan showing a film driving system
- FIG. 14 is a plan showing a frame detecting system
- FIG. 15 an enlarged plan showing a part of the frame detecting system of FIG. 14.
- the reference numeral shows a reference called the median line, numeral 11 the heart and numerals 12 and 12' lungs.
- the numerals 13 and 14 show lines drawn perpendicular to the median line on the left and right sides thereof respectively at such points where the heart is projected most to the left and right I f the lengths of t l 1e tw g lines are expressed as T3 and 14 respectively, 13 14 represents the maximal transverse diameter-of the heart.
- the numeral 15 shows the maximal transverse internal diameter of the thorax the length of which is expressed as l 5 Therefore, the
- FIG. 2 an electrical signal 23 obtained by photoelectric conversion of an X-ray image ofthe chest along a given straight line at right angles to the median line 10 is shown in a simplified form, the numerals 21 and 22 respectively showing the white and black levels.
- FIG. 2 is concerned with an ordinary image with a dark lung field, its brightness may be reversed.
- the conventional straightward method of determin-' ing the peripheral points of the heart and thorax consisted in differentiating the electrical signal 23 to obtain points 24-1, 24-2, 24-3 and 24-4 with great gradients or finding the points including 25-1, 25-2, 25-3 and 25-4 on the signal curve where it crosses a predetermined lgyel if line.
- the data thus obtained was used to obtain l3, l4 and 1 5 and then to compute the car'diothoracic ratio.
- the electrical signal 23 of FIG. 2 usually takes a more complex form as shown in FIG. 3.
- the curve of FIG. 3 shows an example ofa typical electrical signal on which numerals 31, 32, 33 and 34 indicate the peripheral points corresponding to those on the straight line 20 as determined by the observation of a person with a wealth of medical knowledge and experience.
- the points 31, 32, 33 and 34 are considered to be true peripheral points. If the first one of the conventional methods explained with reference to FIG. 2 in which an electrical signal is differentiated to find points with the greatest gradients is applied to curve of FIG. 3 it is possible to detect the points on the curve of FIG.
- peripheral points detected by the conventional methods do not agree with the correct peripheral points in many cases, sometimes even with an error for beyond the practically allowable limit. This is the very reason why, in spite of the simple engineering processes involved, the conventional methods have not been successfully commercialized.
- Points 40 and 40' show the ends of the curve 30.
- Numeral 46 shows the maximum point on the gentle part of the curve almost at the center thereof. This maximum point 46 of the curve 30 which is corresponding to the straight line 20 is found in every ordinary vX-ray image of the chest.
- FIGS. 5(A), 5(8) and 5(C) show how the peripheral points 41, 42, 43 and 44 determined according to the invention are almost free from influences of any distortions of an electrical signal. According to the conventional methods, by contrast, even if satisfactory results. may be obtained from the signal of the waveform shown in FIG. 5(A), distortions of the signal as shown in FIGS.
- a set of peripheral points thus determined will as-- merals 60-1, 60-2, 60-3, 60-n show equidistant straight lines at right angles to the median line which are equivalent to the straight line 20 moved up and down.
- the distance between the straight lines which is determined taking into consideration the required accuracy and economy may be much longer than in the television raster.
- the peripheral points 61-1, 62-1, 63-1 and 64-1; 61-2, 62-2, 63-2 and 64-2; 61-3, 62-3, 61-n, 62-n, 63-n and 64-n determined by the abovementioned method are the series of data providing the profiles of the thorax and the heart.
- the line 62-n on the straight line 60-n is known to be a peripheral point of the diaphragm 69 but not the peripheral point between the lungs and the heart, the line 62-n will be disregarded, as mentioned in detail later with reference to th e second step of operation, in gle term in ing th e length 62-1.
- 63-1, 63-2, 63-3, 65 show the length of lines perpendicular to the median line and passing through the points 6i 1 ,63-2, 63-n respectively
- the maximum length 63-J corresponding to 14 in FIG. 1 is determined without regard to 637.
- the distange between the points 61-1 and 64-1 is expressed as 65-]
- the distance between points 61-2 a r i d 64-2 as 5 and so on, the maxim n distance 65-K is determined which corresponds to 15 in FIG. 1.
- t h e ca rd i othorzgc ratiol computed from the length 62-1, 63-J and 65 -K as 62-I +63 J/65- R X 100 percent.
- the cardiothoracic ratio is obtained only after the correction of the profile at the second step of the operation which will be explained below.
- peripheral points 41, 42, 43 and 44 determined in FIGS. 4 and 5 according to the first step of operation agree well with the true peripheral points and that the results obtained are very stable. Even in this case, however, an error sometimes occurs which exceeds an allowable limit. This is mainly attributable to the fact that as mentioned above the shades of the ribs and blood vessels are indistinguishable from that of the lung field or the heart. The result is an electrical signal representing false information indicating an erroneous peripheral point. As far as the observation.
- the straight lines -1, 70-2, 70-6 in FIG. 7A correspond to the straight line 20.
- the peripheral points 72-1, 72-2, 72-6 were determined by following the first step of operation, the point 724 being a false peripheral point determined as mentioned above.
- the co-ordinates representing the points 72-1, 72-2, 72-6 are D D D respectively.
- FIG. 7B shows the case in which D does not actually exist within the range of P i A estimated on the basis of D and P
- D is again employed for a further procedure disregarding P
- Such a case occurs in the neighborhood of the diaphragms 69 and 69' in FIG. 6.
- the data series tend to depart from the median line sharply in the neighborhood of the diaphragm, and it is known that points 72-4, 72-5, 72-6, correspond to the diaphragm, so that these data may be eliminated in the process of determining the maximal transverse diameter of the heart.
- FIGS. 7A and 7B are primarily concerned with the determination of the profile of the heart, but a similar correcting process is applied to the determination of the profile of the thorax.
- the method of correction explained with reference to FIGS. 7A and 7B is onlyan example of many possible methods, and does not preclude the use of a combination of more complex techniques of prediction. The important thing is to verify the co-relation with points in the neighborhood.
- the present invention may be also effectively applied to the automatic measurements of not only the size of various internal organs but the area occupied thereby as its image is projected in a certain direction, as well as the volume of the internal organ through the combination of the above data and the ratio between them on the basis of radiographs.
- the reference numeral 800 shows a light source including a lamp 801, a light diffusion plate 802 and a box 803 serving simultaneously as a reflector.
- the light source 800 which is so designed that almost uniform brightness is obtained over the whole area ofthe light diffusion plate 802 is identical with the one used by a doctor for the observation of an X.-ray film.
- the X-ray film-810 is disposed near the light diffusion plate 802.
- Numeral 820 shows a photo-electric conversion system including a lens 821, vidicon tube 822, deflection circuit 823, video signal amplifier 824, sample and hold circuit 825, A-D converter 826 and control circuit 827.
- the lens 821, vidicon tube 822, deflection circuit 823 and video signal amplifier 824 or a part thereof are equivalent to and may be replaced by a standard type of television camera with 525 scanning lines.
- the lens 821 should preferably be provided with an automatic iris control regulated by, say, a CdS sensor in the neighborhood of the lens. Theiris of thelens is automatically adjusted to maintain the good operating conditions of the vidicon tube and the video signal amplifier even if the average shade of the film image varies with the photographing and developing conditions or an object.
- the output of the video signal amplifier 824 is applied through the sample and hold circuit 825 to the A-D converter 826 where parallel digital outputs are produced.
- the number of levels into which the shade is quantized in the digital data processing of an image depends upon the object of the data processing, but 7 bits that is 128 levels or 8 bits that is 256 levels meet the requirements of the present invention. In other words, one datum is containedin one byte.
- FIG. 9. the output of the video signal amplifier 824 is applied to the lead wire 9-1.
- the reference numeral 92 shows an operational amplifier, numerals 93 and 94 MOS FETs, numeral 95 a drive circuit for supplying a sampling signal to the gateof the MOS FET. 93.
- a timing signal for the sampling operation is introduced from the control circuit 827 and applied to the lead wire 96.
- a high input imedpance of the gate of the MOS FET 94 and capacitor 97 constitute a hold circuit, a signal held by which is taken out of the source of the MOS FET and applied through the lead wire 98 to the A-D converter 826. No explanation will be made of the A-D converter 826 which is well known.
- the control circuit 827 acts to control not only the sample and-hold circuit 825 and A-D converter 826 but the deflection circuit 823 and is shown in detail in FIG. 10.
- the numeral 101 shows a clock pulse generator, numeral 102 a countdown circuit, numeral 103 a lead wire for supplying horizontal and vertical synchronizing pulses to the deflection circuit 823 of the vidicon tube 822, and numeral 104 a dynamic shift register for a circulating storage with a circulation time equal to a horizontal scanning period, which is provided with output terminals for designating the positions of the parallel straight lines 60-1, 60-2, 60-3, 60-n in FIG. 6.
- Numeral 105 shows a shift register for determining the number of the 1 parallel straight lines 60-1, 60-2, 60-3, 60-n in FIG.
- numeral 106 a timing circuit
- numeral 107 a lead wire for introducing a data read start instruction signal of image pick-up within which the vidicon tube operates
- numeral 112 a locus of the electron beam scanning.
- the locus'in FIG. 11 is shown rotated by even though actually its aspect ratio is 3 to 4 and agrees with a standard television system as already explained.
- The-interlaced scanning ratio of 2 to l in the standard television system is not employed in this embodiment as will be clear from FIG. 10, partly because a resolution of about 128 lines is sufficient for the purpose of the present invention, and partly because it isdesired that the control circuit'827 should not unnecessarily be complicated.
- the clock pulse generator 101 produces an output at the frequency of 1 MHz, while the horizontal and vertical synchronizing I signals produced from the count-down circuit 102 have the period of 64p sec or the frequency of [5625 KHZ and 16.384 in sec or 61.035 Hz respectively. Strictly speakingQtherefore, they are different from the standard system but it is needless to say that this small difference in frequency does not affect the normal operation of the standard type of camera or monitor. This merit, however, is not an important element of the present invention.
- the above-mentioned method of obtaining the virtual scanning lines-by the standardtype of scanning and the controlling of timing of the sampling operation is only one of the possible methods.
- Another method of achieving the object of the present invention is to deflect the electron beam along the above-mentioned vertical scanning lines. In other words, what is important is to convert the brightness of an image into an electrical signal along such loci.
- the digital signal thus obtained is processed on the aforementioned principle by the data processing means 830 which comprises an interface 831, random access memory 832, processor unit 833 and read only memory 834.
- the data processing means 830 should preferably be provided with an interrupt function, whereby the processing time can be shortened.
- This data processing means may be replaced by a small digital computer which is generally called a minicomputer, in which case the program may be written not in the read only memory 834 but in the random access memory 832 with a capacity of approximately 4,000 bytes or words.
- the numeral 840 in FIG. 8 shows an output means for producing the computation results of the cardiothoracic ratio.
- Such an output means is typically represented by an input-output typewriter such as Teletype ASR33 of Teletype corporation or a CRT display terminal.
- the profile of the heart and therefore the silhouette thereof as shown in FIG. 12 which is an interim result of the data processing operation is easily printed or otherwise displayed as'desired by a doctor involved.
- the availability of this option is very important as it helps the doctor, if he so desires, diagnose a heart disease especially when the cardiothoracic ratio is abnormally high.
- the reference numeral 850 in FIG. 8 shows a film handling system including a film handling mechanism 851 and a control circuit 852.
- the types of film handled are roughly divided into a sheet and a roll.
- a film handling system for a sheet is similar in construction to an OCR system called a document reader, while a roll film handling system is like an OCR system generally called a journal tape reader.
- sheet filmsin an input hopper are fed one by one to automatically measure the cardiothoracic ratio following the above-mentioned operating steps, while this system is so arranged that on completion of a processing the processed films are housed in a stacker to make ready for the next film processing.
- a roll of film is fed from one spool to the other -by means of a pinch roller, and it is desirable that the film is maintained stationary at least during the period of time in which a frame of film is being processed or when the film is being scanned to derive information.
- the shaft 131 of the pinch roller is connected to the shaft of the motor 133 through the electromagnetic clutch 132 or the shaft of the electromagnetic clutch/brake 134 through the electromagnetic clutch 132.
- the time when the film should be stopped occurs when the center of a frame of film to be processed reaches the optical axis of the lens 821, as detected by photo'electric sensors.
- the vidicon 822, lens 821, film 810, light diffusion plate 802 and lamp 801 in FIG. 14 are identical to those shown in FIG. 8, the numerals 141 and 142 showing the photo-electric sensors.
- the numeral 810 shows a roll film which is fed in the direction of arrow 143.
- FIG. 15 The relationship between a series of frames of the film 810 and the photo-electric sensors 141 and 142 as seen from the lens is shown in FIG. 15.
- the distance 151 between the two photo-electric sensors 141 and 142 is equal to the length 152 of a frame of the recorded roll film 810.
- either the control circuit 852 is energized by a signal from the data processing means 830, or the operator in response to an alarm from the displayer indicating the completion of the processing 10 depresses a button for feed instructions, thereby to restart the feed of the" roll film 810 in the direction of arrow 143, but it is not until the next frame 8102 reaches the mid-point between the sensors 141 and 142 that the outputs of both the sensors 141 and 142 change to the white level.
- the output of the sensors is shaped and, after an AND operation thereof, the control circuit 852 is energized, whereupon the film 810 is stopped by the electromagnetic clutch/brake 134.
- the film handling system 850 may be replaced by a manual system in which the film is fed by hand sheet by sheet or frame by frame, as the occasion may be.
- a system for automatically measuring the cardiothoracic ratios comprising photo-electric conversion means (820) for effectively scanning a chest X-ray image (810) along a plurality of parallel lines perpen dicular to the median line of said chest X-ray image and converting the brightness of said X-ray image into a corresponding electrical signal, data processing means (830) for processing data obtained from said photo electric conversion means, and means (840) for recording or displaying the value of the cardiothoracic.
- a system for automatically measuring the cardiov thoracic ratios according to claim 1, in which said film handling means is manually operated.
- a system for automatically measuring the cardiothoracic ratios according to claim 1, in which said photo-electric conversion means comprises a standard type of television image pickup device or an equivalent film camera.
- a system for automatically measuring the cardiothoracic ratios according to claim 1, in which said data processing means is replaced by a general-purpose digital computer.
- a system for automatically measuring the cardio-' thoracic ratios accordingto claim 1, in which data are processed in such a manner as to display or record the profiles of the heart and/or the thorax or the silhouettes thereof as well as the cardiothoracic ratios.
- a system for automatically measuring the cardiothoracic ratios according to claim 1, further comprising film handling means for automatically feeding the film carrying said chest X-ray image.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- General Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- Geometry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Cardiology (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Biochemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Input (AREA)
Abstract
An apparatus for automatically measuring the cardiothoracic ratios speedily and with a sufficiently high accuracy, characterized in that the profiles of the thorax and heart are first inferred by processing information on brightness from an Xray image and then an error in the inference of the profiles is corrected by contextual data processing, so that the cardiothoracic ratio may be computed on the basis of the corrected data.
Description
United States Patent I [111 3,869,602 Sezaki et al. Q 1 V y 4 I Mar. 4, 1975 I 4] APPARATUS FOR AUTOMATIC 'Fo .P COMPUTATION or CARDIOTHORACIC REIGN ATENTS OR APPLICATIONS 611,258 12/1960 Canada l78/DlG. 5 RATIO 723,256 12/1965 Canada 178/DlG. 5 [75] Inventors; Nobuhiko s Kohji Ukena, both 869,706 6/1961 Great Britain 178/D1G. 5
Y of Yokohama, Japan [73] Assignee: Matsushita Electric Industrial Co., Primary Examiner.loseph F. Ruggiero 7 Ltd., Osaka, Japan Armrney, Agent, or FirmStevens, Davis, Miller & 221 Filed: Aug. 29, 1972 Mosh [21] App]. No.: 284,558
57 ABSTRACT [30] Foreign Application Priority Data Sept. 3, 1971 Japan 46-68426 An apparatus for automatically measuring the cardiothoracic ratios speedily and with a sufficiently high ac- [52] US. Cl 235/l5l.3, 128/205 R, l78/DlG. 5 curacy, characterized in that the profiles of-the thorax [51] Int. Cl. G06f /42 and heart are first inferred by processing information [58] Field of Search 235/1513; 444/1; on brightness from an X-ray image and h n n rr r 128/205 R, 2.06 R, 2.06 B, 2.6 G, 2.1 R, in the inference of the profiles is corrected by contex- 419 R, DIG, 3; 178/DIG, 5, DIG, 22; tual data. processing, so that the cardiothoracic ratio 250/320, 321, 322, 323, 401 may be computed on the basis of the corrected data.
[56] References Cited 6 Claims, 16 Drawing Figures UNITED STATES PATENTS 3,580,997 5/1971 Webb l78/DIG. 5 UX H 5E p y 824 ,;825 826 ag 8 l f AMPLIFIER ClRCUlT CONVERTER l I =80l I DEFLECTlON CONTROL 802 CIRCUIT CIRCUIT m i i FILM 1 HANDLING 850 MEANS l l CONTROL L i i, CIRCUIT 852 RANDOM L ACCESS MEMORY 83o INTER- PROCESSOR FACE UNIT 840 i OUTPUT 1 MEMORY I MEANS 833 834 PATENTED 41975 3.869.602
sum 2 9f 9 FIG. 3
FIG. 4
PATENTED 41975 snwanrs) PATENTED 4 7 SHEET 8 BF 9 OOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO O00OOOOO00000000000000OOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00000000000OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOO0O00OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000000000000OOOOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOO PATENTEDH 41975 3.869.602
sum 9 {3f 9 FIG. l3
1 APPARATUS FOR AUTOMATIC COMPUTATION OF CARDIOTHORACIC RATIO The present invention relates to an apparatus for automatically computing the cardiothoracic ratio especially on the basis of an X-ray photograph.
The recent years have seen a great increase in the number of persons in the advanced nations who have died of cancer or malignant tumor, heart diseases and other adult diseases, raising the need for improved preventive medicine. A reinforcement of the present system of mass examination may be an answer. One object of the mass examination is to find and pick up out of examinees those diseased persons who have no subjective symptoms. The achievement of this purpose requires the processing of a large volume of medical data collected. This job which is known as a screening, even though very monotonous, must be endorsed by a deep knowledge as well as a wealth of experience in the field of medicine, demanding a great time and labor of a doctor involved. Nevertheless, the shortage of the number of doctors poses a serious problem which cannot be solved in the near future.
It is common that a heart disease is accompanied by an abnormal enlargement of the heart, and the universally adopted method of measurement thereof is by the cardiothoracic ratio advocated by C. S. Danzer, which is defined as a ratio in percentage of the transverse diameter of the heart to the maximal transverse internal diameter of the thorax. (See, for example, the Japanese Medical Journal No. 2433, PP. 130 and No. 2,422, PP. 134 to PP. 135). In other words, since the cardiac enlargement is a helpful overall indicator of a heart disease, the CTR is a valuable clinical parameter for the screening process.
The CTR is usually determined with a rule or similar means applied to an X-ray or equivalent image of the chest. Even though various devices have been devised in an effort to facilitate the measurement, no substantial progress has been made in the direction of greatly saving the labor for the measurement.
It is easily imagined that the computation of the CTR will be made easily by a hypothetical method in which the brightness of an image is converted into an electrical signal by means of an appropriate photoelectric converter and portions of the signal with a large gradient are located by differentiating them or points are determined where the signal crosses a predetermined leve, in order to obtain the profiles of the thorax and heart. The past attempts made in this direction, however, did not necessarily agree with the actual results of measurements, thereby preventing the widespread use of this method. H. P. Meyers et al reported in l964 that the cardiothoracic ratio automatically computed with a flying spot scanner and a digital computer almost agreed with the actual measurements (H. P. Meyers et al. Automated Computer Analysis of Radiographic Images, Radiology Vol. 83, PP, 1,029 to 1,034, Dec. 1964). But their algorithm is based on a principle different from the one mentioned above-and is entirely empiric. Anyway, the amount of information contained in a sheet of X-ray film is so great that its complete processing requires either very expensive hardware or I greatly prolonged time.
2 speedly and with a sufficiently high accuracy without human labor. I V
The principal operations or operating principle of the apparatus according to the present invention consist of the first step of processing brightness information included in an X-ray image by an improved method to infer the profiles of the thorax and heart, the second step of correcting an error, if any, in the inference by contextual data processing and the third step of computing the cardiothoracic ratio on the basis of the corrected data.
Another object of the present invention is to provide an automatic cardiothoracic ratio measuring apparatus with a high performance at a reasonable cost.
Still another object of the invention is to automate the screening operation eitherby recording or displaying the results of measurements of the cardiothoracic ratios of many examinees which have been obtained as a result of automatically and successively processing the X-ray images of their chests, or by identifying for recording or display purposes the serial numbers of the examinees chest images with a cardiothoracic ratio amost equal to or higher than a reference or producing an alarm. signal when their cardiothoracic ratios are almost equal to or greater than a predetermined reference, thereby to save the labor of a doctor and -elimi nate not only the differences in the results of measurements, which otherwise might occur for individual doctors, but their exhaustion for improved objectivity, reliability and repeatability of the measurements.
The above and other objects, features and advantages will be made apparent by the detailed description taken in conjunction with the accompanying drawings, in which:
FIG. 1, is a diagram for explaining what is the cardiothoracic ratio;
FIG. 2 is a diagram showing an electrical signal generated by photo-electric conversion of an X-ray chest image along a straight line at right angles to the median line;
Flg. 3 is a diagram showing a typical electrical signal obtained from an actual X-ray chest image;
FIG. 4 is a diagram for explaning the manner in which the peripheral points are determined in a typical electrical signal;
FIG. 5 is a diagram similar to FIG. 4 for explaining the manner in which the peripheral points are determined in a distorted electrical signal;
FIG. 6 is a diagram for explaining the manner in which the transverse diameter of the heart and the maximal transverse internal diameter of the thorax are determined; 1
FIGS. 7A and 7B are diagrams for explaining the manner in which data obtained are corrected;
FIG. 8 is a block diagram showing an embodiment of the present invention;
FIGS. 9 and 10 are diagrams showing in detail a part of the embodiment of FIG. 8;
FIG. 11 shows a scanning locus;
heart;
FIG. 13 is a plan showing a film driving system;
FIG. 14 is a plan showing a frame detecting system; and
FIG. 15 an enlarged plan showing a part of the frame detecting system of FIG. 14.
First, preliminary explanation will be made of an X-ray image of the chest and the cardiothoracic ratio.
Referring to FIG. 1 drawn for explaining what is the cardiothoracic ratio, the reference numeral shows a reference called the median line, numeral 11 the heart and numerals 12 and 12' lungs. The numerals 13 and 14 show lines drawn perpendicular to the median line on the left and right sides thereof respectively at such points where the heart is projected most to the left and right I f the lengths of t l 1e tw g lines are expressed as T3 and 14 respectively, 13 14 represents the maximal transverse diameter-of the heart. The numeral 15 shows the maximal transverse internal diameter of the thorax the length of which is expressed as l 5 Therefore, the
' cardiothoracic ratio is given as [(E+1 I)/ij 100 percent.
In FIG. 2, an electrical signal 23 obtained by photoelectric conversion of an X-ray image ofthe chest along a given straight line at right angles to the median line 10 is shown in a simplified form, the numerals 21 and 22 respectively showing the white and black levels. Although FIG. 2 is concerned with an ordinary image with a dark lung field, its brightness may be reversed.
It is assumed in all the descriptions below that the lung field is shown dark and the white and black levels are positioned at upper and lower sides respectively. It should be noted, however, that these assumptions do not adversely affect the possibility of general application of the invention. I
The conventional straightward method of determin-' ing the peripheral points of the heart and thorax consisted in differentiating the electrical signal 23 to obtain points 24-1, 24-2, 24-3 and 24-4 with great gradients or finding the points including 25-1, 25-2, 25-3 and 25-4 on the signal curve where it crosses a predetermined lgyel if line. The data thus obtained was used to obtain l3, l4 and 1 5 and then to compute the car'diothoracic ratio.
Actually, however, partly because an X-ray image of the chest contains complicated intermediate tones with irregular variations in shade andpartly because the images of ribs and blood vessels overlap each other in the lung field, the electrical signal 23 of FIG. 2 usually takes a more complex form as shown in FIG. 3. The curve of FIG. 3 shows an example ofa typical electrical signal on which numerals 31, 32, 33 and 34 indicate the peripheral points corresponding to those on the straight line 20 as determined by the observation of a person with a wealth of medical knowledge and experience. The points 31, 32, 33 and 34 are considered to be true peripheral points. If the first one of the conventional methods explained with reference to FIG. 2 in which an electrical signal is differentiated to find points with the greatest gradients is applied to curve of FIG. 3 it is possible to detect the points on the curve of FIG. 3 corresponding to the points 24-1, 24-2, 24-3 and 24-4 on the curve of FIG. 2 but additional points which are almost indistinguishable from the true peripheral points are unnecessarily detected. On the other hand, the application to the curve 30 of the second conventional method in which points where the signal crosses a preset reference level are located has the disadvantage of the difficulty in presetting the reference level and often leads to the result that the detected peripheral points are highly variable when the shade gradient is gentle,
resulting in the failure to detect the true peripheral points depending upon the shape of the curve 30 and the reference level-setting process.
The peripheral points detected by the conventional methods do not agree with the correct peripheral points in many cases, sometimes even with an error for beyond the practically allowable limit. This is the very reason why, in spite of the simple engineering processes involved, the conventional methods have not been successfully commercialized.
The operating principle of the present invention will be now explained. Referring first to the first step of the operation, the manner is described below in which the peripheral points are determined with reference to the curve 30 of FIG. 4 showing a typical electrical signal.
When the curve 30 is followed starting at the end point 40, the first point 41 is reached where the gradient is negative. This point 41 corresponds to the point 31 of FIG. 3. Following the curve 30 further toward the point 46, a minimum point 45 is reached before point 46. Intermediate between point 45 and point 46, there always exists a point with the maximum gradient which corresponds to point 32 of FIG. 3 and may be expressed as point 42. In like manner, points 43 and 44 corresponding to points 33 and 34 respectively on the curve of FIG. 3 are determined by following the curve 30 from the other end point 40'.
The points 41, 42, 43 "and 44 'thus determined agree with the true peripheral points very well. Further, the process of determining the points is stable'and simple, and a measurement error is sufficiently within the limit of commercial applications excepting the cases men tioned later. A qualitative analysis of this advantage is shown in FIGS. 5(A), 5(8) and 5(C), FIG. 5(A) being quite identical with FIG. 4. FIGS. 5(B) and 5(C) show how the peripheral points 41, 42, 43 and 44 determined according to the invention are almost free from influences of any distortions of an electrical signal. According to the conventional methods, by contrast, even if satisfactory results. may be obtained from the signal of the waveform shown in FIG. 5(A), distortions of the signal as shown in FIGS. 5(8) and 5(C) will result in an increased error or make it impossible to detect the peripheral points. Photo-electric converters now commonly used include the flying spot scanner and the vidicon camera, although a mechanical scanning means is also available. Regardless of which means is employed, however, it is well known that the distortions as shown in FIG. 5 are a frequent phenomenon due to various causes. It is also expected that an X-ray image itself carries such distortions. Further, although not referred to in the explanation of FIG. 5, it is evident that the peripheral points determined by the method of the invention are affected very little by variations in the amplitude and DC level of an electrical signal resulting from the variations in the contrast of the X-ray image or the gain of an amplifier.
A set of peripheral points thus determined will as-- merals 60-1, 60-2, 60-3, 60-n show equidistant straight lines at right angles to the median line which are equivalent to the straight line 20 moved up and down. The distance between the straight lines which is determined taking into consideration the required accuracy and economy may be much longer than in the television raster. The peripheral points 61-1, 62-1, 63-1 and 64-1; 61-2, 62-2, 63-2 and 64-2; 61-3, 62-3, 61-n, 62-n, 63-n and 64-n determined by the abovementioned method are the series of data providing the profiles of the thorax and the heart.
Before going to the second step, explanation will be made now of the third step of operation with reference to FIG. 6 for convenience of illustration.
In order to compute the cardiothoracic ratio, it is necessary first of all to determine 13, 13 and I? as explained above with reference to Flg. l.-Assuming that 623, 62 2, 623, m show the length of the lines perpendicular to the median line and passing the po int 62-1, 62-2, 62-3, 62-n respec t iv e ly and that 62-1 is the longest of them, the length 62-1 corresponds to the length 13 shown in FIG. 1. In this particular illustration, however, since the point 62-n on the straight line 60-n is known to be a peripheral point of the diaphragm 69 but not the peripheral point between the lungs and the heart, the line 62-n will be disregarded, as mentioned in detail later with reference to th e second step of operation, in gle term in ing th e length 62-1. Similarly, assuming that 63-1, 63-2, 63-3, 65 show the length of lines perpendicular to the median line and passing through the points 6i 1 ,63-2, 63-n respectively, the maximum length 63-J corresponding to 14 in FIG. 1 is determined without regard to 637. Furthermore, if the distange between the points 61-1 and 64-1 is expressed as 65-], the distance between points 61-2 a r i d 64-2 as 5 and so on, the maxim n distance 65-K is determined which corresponds to 15 in FIG. 1.
Accordingly, t h e ca rd i othorzgc ratiol computed from the length 62-1, 63-J and 65 -K as 62-I +63 J/65- R X 100 percent. Actually, however, the cardiothoracic ratio is obtained only after the correction of the profile at the second step of the operation which will be explained below.
It was already explained that the peripheral points 41, 42, 43 and 44 determined in FIGS. 4 and 5 according to the first step of operation agree well with the true peripheral points and that the results obtained are very stable. Even in this case, however, an error sometimes occurs which exceeds an allowable limit. This is mainly attributable to the fact that as mentioned above the shades of the ribs and blood vessels are indistinguishable from that of the lung field or the heart. The result is an electrical signal representing false information indicating an erroneous peripheral point. As far as the observation. is one-dimensional, this error is inevitable, and it is easily expected that even an observer with a wealth of medical knowledge and experience will make a similar mistake if he is to observe an X-ray image of the chest through the slits perpendicular to the median line. Nevertheless he is usually able to infer the profiles of the thorax and the heart accurately, this is not only due to his knowledge and experience but to his ability to correct the superficial results of the one-dimensional observation by making reference to all the factors involved, thereby effectively inferring the original two dimensional picture. Along this line, a contextual data processing step was followed in correcting an erroneous determination of the profile according to the present invention. The operation of the present invention is based on the principle of the prediction method as shown in FIGS. 7A and 7B to correct an error and identify the diaphragm. v
The straight lines -1, 70-2, 70-6 in FIG. 7A correspond to the straight line 20. Let us assume that the peripheral points 72-1, 72-2, 72-6 were determined by following the first step of operation, the point 724 being a false peripheral point determined as mentioned above. Also, assume that the co-ordinates representing the points 72-1, 72-2, 72-6 are D D D respectively. It In a given series of data, it is generally known how to predict the probable value of D on the basis of D and D The simplest method of such a prediction is the linear prediction method in which the predicted value of D(i) is given as P i A, if D 1 and D are known where P 2 X D,,-.,, D and the symbol A shows an allowable limit which will depend on the standard shape of the heart and the spaces between the straight lines 70-1, 7.0-2, and 70-6, for example.
In the example of FIG. 7A, if D departs from the range predicted on the basis of D and D D is provisionally disregarded and P is substituted therefor. and then it is examined whether D actually exists within the range of P(5) A'estimated on the basis of D and P If so, the preceding provisional substitution is deemed reasonable and a further procedure is followed, after completing the correction of D The point 72-4' on the-straight line 70-4 shows the result thus corrected, the thick portion of the straight line indicating the range covering i A. The point 72-5' is provided for the purpose of the above-mentioned examination and meaningless once the presence of point 72-5 within the range of i A is verified.
FIG. 7B shows the case in which D does not actually exist within the range of P i A estimated on the basis of D and P In this case, the provisional substitution is not deemed reasonable and D is again employed for a further procedure disregarding P Such a case occurs in the neighborhood of the diaphragms 69 and 69' in FIG. 6. The data series tend to depart from the median line sharply in the neighborhood of the diaphragm, and it is known that points 72-4, 72-5, 72-6, correspond to the diaphragm, so that these data may be eliminated in the process of determining the maximal transverse diameter of the heart. FIGS. 7A and 7B are primarily concerned with the determination of the profile of the heart, but a similar correcting process is applied to the determination of the profile of the thorax. Incidentally, the method of correction explained with reference to FIGS. 7A and 7B is onlyan example of many possible methods, and does not preclude the use of a combination of more complex techniques of prediction. The important thing is to verify the co-relation with points in the neighborhood.
Although the above explanation of the operating principle of the present invention was with reference to' the automatic measurement of the cardiothoracic ratio, the present invention may be also effectively applied to the automatic measurements of not only the size of various internal organs but the area occupied thereby as its image is projected in a certain direction, as well as the volume of the internal organ through the combination of the above data and the ratio between them on the basis of radiographs.
It will be understood from above that it is possible to embody the present invention in various forms,'one of which will be explained below in detail with reference to FIG. 8.
In this figure, the reference numeral 800 shows a light source including a lamp 801, a light diffusion plate 802 and a box 803 serving simultaneously as a reflector. The light source 800 which is so designed that almost uniform brightness is obtained over the whole area ofthe light diffusion plate 802 is identical with the one used by a doctor for the observation of an X.-ray film. The X-ray film-810 is disposed near the light diffusion plate 802.Numeral 820 shows a photo-electric conversion system including a lens 821, vidicon tube 822, deflection circuit 823, video signal amplifier 824, sample and hold circuit 825, A-D converter 826 and control circuit 827. The lens 821, vidicon tube 822, deflection circuit 823 and video signal amplifier 824 or a part thereof are equivalent to and may be replaced by a standard type of television camera with 525 scanning lines. The lens 821, however, should preferably be provided with an automatic iris control regulated by, say, a CdS sensor in the neighborhood of the lens. Theiris of thelens is automatically adjusted to maintain the good operating conditions of the vidicon tube and the video signal amplifier even if the average shade of the film image varies with the photographing and developing conditions or an object.
The output of the video signal amplifier 824 is applied through the sample and hold circuit 825 to the A-D converter 826 where parallel digital outputs are produced. The number of levels into which the shade is quantized in the digital data processing of an image depends upon the object of the data processing, but 7 bits that is 128 levels or 8 bits that is 256 levels meet the requirements of the present invention. In other words, one datum is containedin one byte. The detail in this connection is shown in FIG. 9. In this figure, the output of the video signal amplifier 824 is applied to the lead wire 9-1. The reference numeral 92 shows an operational amplifier, numerals 93 and 94 MOS FETs, numeral 95 a drive circuit for supplying a sampling signal to the gateof the MOS FET. 93. A timing signal for the sampling operation is introduced from the control circuit 827 and applied to the lead wire 96. A high input imedpance of the gate of the MOS FET 94 and capacitor 97 constitute a hold circuit, a signal held by which is taken out of the source of the MOS FET and applied through the lead wire 98 to the A-D converter 826. No explanation will be made of the A-D converter 826 which is well known. The control circuit 827 acts to control not only the sample and-hold circuit 825 and A-D converter 826 but the deflection circuit 823 and is shown in detail in FIG. 10. In this figure, the numeral 101 shows a clock pulse generator, numeral 102 a countdown circuit, numeral 103 a lead wire for supplying horizontal and vertical synchronizing pulses to the deflection circuit 823 of the vidicon tube 822, and numeral 104 a dynamic shift register for a circulating storage with a circulation time equal to a horizontal scanning period, which is provided with output terminals for designating the positions of the parallel straight lines 60-1, 60-2, 60-3, 60-n in FIG. 6. Numeral 105 shows a shift register for determining the number of the 1 parallel straight lines 60-1, 60-2, 60-3, 60-n in FIG.
6, numeral 106 a timing circuit, numeral 107 a lead wire for introducing a data read start instruction signal of image pick-up within which the vidicon tube operates, and numeral 112 a locus of the electron beam scanning. To make it easier to understand the relationship with FIG. 6, the locus'in FIG. 11 is shown rotated by even though actually its aspect ratio is 3 to 4 and agrees with a standard television system as already explained. The-interlaced scanning ratio of 2 to l in the standard television system is not employed in this embodiment as will be clear from FIG. 10, partly because a resolution of about 128 lines is sufficient for the purpose of the present invention, and partly because it isdesired that the control circuit'827 should not unnecessarily be complicated. More in detail, the clock pulse generator 101 produces an output at the frequency of 1 MHz, while the horizontal and vertical synchronizing I signals produced from the count-down circuit 102 have the period of 64p sec or the frequency of [5625 KHZ and 16.384 in sec or 61.035 Hz respectively. Strictly speakingQtherefore, they are different from the standard system but it is needless to say that this small difference in frequency does not affect the normal operation of the standard type of camera or monitor. This merit, however, is not an important element of the present invention. The parallel straight lines 60-1, 60-2, 60-3, 60-n, as shown perpendicular to the horizontal scanning lines 112 in FIG. 11, are a locus of points where the signals are sampled and converted from analog into digital under the control of the control circuit 827, being so controlled as to recover aline of data for each field. These parallel straight lines correspond to the parallel straight lines 60-1, 60-2, 60-3, 60-n in FIG. 6. Hereinafter these parallel lines are called virtual scanning lines.
The above-mentioned method of obtaining the virtual scanning lines-by the standardtype of scanning and the controlling of timing of the sampling operation is only one of the possible methods. Another method of achieving the object of the present invention is to deflect the electron beam along the above-mentioned vertical scanning lines. In other words, what is important is to convert the brightness of an image into an electrical signal along such loci. Back to FIG. 8, the digital signal thus obtained is processed on the aforementioned principle by the data processing means 830 which comprises an interface 831, random access memory 832, processor unit 833 and read only memory 834. The capacity of the random access memory 832 5 of the program for carrying out the data processing on the afore-mentioned principle, and in this embodiment it must be about 2,000 bytes.
The data processing means 830 should preferably be provided with an interrupt function, whereby the processing time can be shortened. This data processing means may be replaced by a small digital computer which is generally called a minicomputer, in which case the program may be written not in the read only memory 834 but in the random access memory 832 with a capacity of approximately 4,000 bytes or words. The numeral 840 in FIG. 8 shows an output means for producing the computation results of the cardiothoracic ratio. Such an output means is typically represented by an input-output typewriter such as Teletype ASR33 of Teletype corporation or a CRT display terminal. In any case, the profile of the heart and therefore the silhouette thereof as shown in FIG. 12 which is an interim result of the data processing operation is easily printed or otherwise displayed as'desired by a doctor involved. The availability of this option is very important as it helps the doctor, if he so desires, diagnose a heart disease especially when the cardiothoracic ratio is abnormally high.
The reference numeral 850 in FIG. 8 shows a film handling system including a film handling mechanism 851 and a control circuit 852. The types of film handled are roughly divided into a sheet and a roll. A film handling system for a sheet is similar in construction to an OCR system called a document reader, while a roll film handling system is like an OCR system generally called a journal tape reader. In this case of a sheet film handling system, sheet filmsin an input hopper are fed one by one to automatically measure the cardiothoracic ratio following the above-mentioned operating steps, while this system is so arranged that on completion of a processing the processed films are housed in a stacker to make ready for the next film processing. In the case of a roll film handling system, on the other hand, a roll of film is fed from one spool to the other -by means of a pinch roller, and it is desirable that the film is maintained stationary at least during the period of time in which a frame of film is being processed or when the film is being scanned to derive information. Accordingly, as shown in FIG. 13, the shaft 131 of the pinch roller is connected to the shaft of the motor 133 through the electromagnetic clutch 132 or the shaft of the electromagnetic clutch/brake 134 through the electromagnetic clutch 132. The time when the film should be stopped occurs when the center of a frame of film to be processed reaches the optical axis of the lens 821, as detected by photo'electric sensors. The vidicon 822, lens 821, film 810, light diffusion plate 802 and lamp 801 in FIG. 14 are identical to those shown in FIG. 8, the numerals 141 and 142 showing the photo-electric sensors. The numeral 810 shows a roll film which is fed in the direction of arrow 143.
The relationship between a series of frames of the film 810 and the photo- electric sensors 141 and 142 as seen from the lens is shown in FIG. 15. In this figure, the distance 151 between the two photo- electric sensors 141 and 142 is equal to the length 152 of a frame of the recorded roll film 810. On completion of the processing one frame 8101, either the control circuit 852 is energized by a signal from the data processing means 830, or the operator in response to an alarm from the displayer indicating the completion of the processing 10 depresses a button for feed instructions, thereby to restart the feed of the" roll film 810 in the direction of arrow 143, but it is not until the next frame 8102 reaches the mid-point between the sensors 141 and 142 that the outputs of both the sensors 141 and 142 change to the white level. The output of the sensors is shaped and, after an AND operation thereof, the control circuit 852 is energized, whereupon the film 810 is stopped by the electromagnetic clutch/brake 134.
The film handling system 850 may be replaced by a manual system in which the film is fed by hand sheet by sheet or frame by frame, as the occasion may be.
Also, it is needless to say that instead of the digital system of data processing employed in the embodiment, an analog or appropriate hybrid system may be used to carry out the operating principles explained in FIGS. 1 to 7.
Furthermore, in place of a transparent image aimed at in the embodiment of FIG. 8, it is possible to con' struct the apparatus according to the present invention in such a manner as to handle an opaque image by slight modification of the light source 800 and the photo-electric conversion system 820.
What is claimed is:
l. A system for automatically measuring the cardiothoracic ratios comprising photo-electric conversion means (820) for effectively scanning a chest X-ray image (810) along a plurality of parallel lines perpen dicular to the median line of said chest X-ray image and converting the brightness of said X-ray image into a corresponding electrical signal, data processing means (830) for processing data obtained from said photo electric conversion means, and means (840) for recording or displaying the value of the cardiothoracic.
2. A system for automatically measuring the cardiov thoracic ratios according to claim 1, in which said film handling means is manually operated.
3. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which said photo-electric conversion means comprises a standard type of television image pickup device or an equivalent film camera.
4. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which said data processing means is replaced by a general-purpose digital computer. i
5. A system for automatically measuring the cardio-' thoracic ratios accordingto claim 1, in which data are processed in such a manner as to display or record the profiles of the heart and/or the thorax or the silhouettes thereof as well as the cardiothoracic ratios.
6. A system for automatically measuring the cardiothoracic ratios according to claim 1, further comprising film handling means for automatically feeding the film carrying said chest X-ray image.
Claims (6)
1. A system for automatically measuring the cardiothoracic ratios comprising photo-electric conversion means (820) for effectively scanning a chest X-ray image (810) along a plurality of parallel lines perpendicular to the median line of said chest X-ray image and converting the brightness of said X-ray image into a corresponding electrical signal, data processing means (830) for processing data obtained from said photo-electric conversion means, and means (840) for recording or displaying the value of the cardiothoracic ratio obtained from said data processing means, said data processing means (830) inferring intersections between the boundary lines of the thorax and the heart and said plurality of parallel lines on the basis of data obtained along each of said parallel lines, said data processing means (830) correcting erroneous data which may be contained in four series of data corresponding to the boundary lines of the thorax and the heart, whereby the cardiothoracic ratio is determined on the basis of said four series of data.
2. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which said film handling means is manually operated.
3. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which said photo-electric conversion means comprises a standard type of television image pickup device or an equivalent film camera.
4. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which said data processing means is replaced by a general-purpose digital computer.
5. A system for automatically measuring the cardiothoracic ratios according to claim 1, in which data are processed in such a maNner as to display or record the profiles of the heart and/or the thorax or the silhouettes thereof as well as the cardiothoracic ratios.
6. A system for automatically measuring the cardiothoracic ratios according to claim 1, further comprising film handling means for automatically feeding the film carrying said chest X-ray image.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP46068426A JPS5116091B2 (en) | 1971-09-03 | 1971-09-03 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| US3869602A true US3869602A (en) | 1975-03-04 |
Family
ID=13373337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US284558A Expired - Lifetime US3869602A (en) | 1971-09-03 | 1972-08-29 | Apparatus for automatic computation of cardiothoracic ratio |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US3869602A (en) |
| JP (1) | JPS5116091B2 (en) |
| CA (1) | CA979074A (en) |
| DE (1) | DE2243449B2 (en) |
| FR (1) | FR2152055A5 (en) |
| GB (1) | GB1400471A (en) |
| NL (1) | NL7211930A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3988602A (en) * | 1975-02-03 | 1976-10-26 | Goodyear Aerospace Corporation | Method and apparatus for enhancing data |
| US4101961A (en) * | 1977-02-16 | 1978-07-18 | Nasa | Contour detector and data acquisition system for the left ventricular outline |
| US4240440A (en) * | 1977-11-17 | 1980-12-23 | Siemens Gammasonics, Inc. | Method and apparatus for nuclear kymography providing a motion versus time display of the outer transverse dimensions of an organ |
| EP0066824A3 (en) * | 1981-06-08 | 1983-02-09 | General Electric Company | X-ray image digital substraction system |
| US4521808A (en) * | 1979-03-22 | 1985-06-04 | University Of Texas System | Electrostatic imaging apparatus |
| US5022063A (en) * | 1989-01-25 | 1991-06-04 | Hitachi Medical Corporation | Multiple-mode scanning and beam current control x-ray TV apparatus |
| US5065435A (en) * | 1988-11-16 | 1991-11-12 | Kabushiki Kaisha Toshiba | Method and apparatus for analyzing ventricular function |
| ES2087824A1 (en) * | 1994-09-16 | 1996-07-16 | Duque Rafael Aparicio | Improved device for the rapid determination of the cardiothoracic index by means of a graduated strip of elastic material |
| FR2790657A1 (en) * | 1999-03-12 | 2000-09-15 | Iodp | Three-dimensional contactless measurement of a tissue of living organism from digital ultrasound image by using graphic representation of profile of gray level of ultrasound image |
| US20090196481A1 (en) * | 2008-02-05 | 2009-08-06 | Huanzhong Li | Image processing method and apparatus |
| CN104732520A (en) * | 2015-01-31 | 2015-06-24 | 西安华海盈泰医疗信息技术有限公司 | Cardio-thoracic ratio measuring algorithm and system for chest digital image |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL7611419A (en) * | 1976-10-15 | 1978-04-18 | Optische Ind De Oude Delft Nv | DEVICE FOR READING AND PROCESSING INFORMATION CONTAINED IN IMAGE FRAMES, SUCH AS FORMED BY IT CONSEQUENTLY FROM MULTIPLE DIRECTIONS WITH A BASICALLY FLAT BEAM OF SHORT-WAVE RADIATION, RADIATION OF AN OBJECT. |
| JPS5530835A (en) * | 1978-08-25 | 1980-03-04 | Takashi Katoda | Method of forming heterogeneous junction in semiconductor |
| JPS57168579A (en) * | 1981-04-08 | 1982-10-16 | Fuji Photo Film Co Ltd | Method and device for reading of picture information |
| JPS5876712A (en) * | 1981-10-31 | 1983-05-09 | Shimadzu Corp | X ray checking device |
| JP6082337B2 (en) * | 2013-10-30 | 2017-02-15 | 横河電機株式会社 | Cardiothoracic ratio calculation device |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3580997A (en) * | 1968-07-05 | 1971-05-25 | Balteau Electric Corp | Video system for automatic production line inspection by x-ray |
-
1971
- 1971-09-03 JP JP46068426A patent/JPS5116091B2/ja not_active Expired
-
1972
- 1972-08-29 US US284558A patent/US3869602A/en not_active Expired - Lifetime
- 1972-09-01 CA CA150,809A patent/CA979074A/en not_active Expired
- 1972-09-01 GB GB4071272A patent/GB1400471A/en not_active Expired
- 1972-09-01 NL NL7211930A patent/NL7211930A/xx unknown
- 1972-09-01 FR FR7231125A patent/FR2152055A5/fr not_active Expired
- 1972-09-04 DE DE2243449A patent/DE2243449B2/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3580997A (en) * | 1968-07-05 | 1971-05-25 | Balteau Electric Corp | Video system for automatic production line inspection by x-ray |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3988602A (en) * | 1975-02-03 | 1976-10-26 | Goodyear Aerospace Corporation | Method and apparatus for enhancing data |
| US4101961A (en) * | 1977-02-16 | 1978-07-18 | Nasa | Contour detector and data acquisition system for the left ventricular outline |
| US4240440A (en) * | 1977-11-17 | 1980-12-23 | Siemens Gammasonics, Inc. | Method and apparatus for nuclear kymography providing a motion versus time display of the outer transverse dimensions of an organ |
| US4521808A (en) * | 1979-03-22 | 1985-06-04 | University Of Texas System | Electrostatic imaging apparatus |
| EP0066824A3 (en) * | 1981-06-08 | 1983-02-09 | General Electric Company | X-ray image digital substraction system |
| US5065435A (en) * | 1988-11-16 | 1991-11-12 | Kabushiki Kaisha Toshiba | Method and apparatus for analyzing ventricular function |
| US5022063A (en) * | 1989-01-25 | 1991-06-04 | Hitachi Medical Corporation | Multiple-mode scanning and beam current control x-ray TV apparatus |
| ES2087824A1 (en) * | 1994-09-16 | 1996-07-16 | Duque Rafael Aparicio | Improved device for the rapid determination of the cardiothoracic index by means of a graduated strip of elastic material |
| FR2790657A1 (en) * | 1999-03-12 | 2000-09-15 | Iodp | Three-dimensional contactless measurement of a tissue of living organism from digital ultrasound image by using graphic representation of profile of gray level of ultrasound image |
| US20090196481A1 (en) * | 2008-02-05 | 2009-08-06 | Huanzhong Li | Image processing method and apparatus |
| US8433112B2 (en) * | 2008-02-05 | 2013-04-30 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for processing chest X-ray images |
| CN104732520A (en) * | 2015-01-31 | 2015-06-24 | 西安华海盈泰医疗信息技术有限公司 | Cardio-thoracic ratio measuring algorithm and system for chest digital image |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS5116091B2 (en) | 1976-05-21 |
| DE2243449A1 (en) | 1973-03-22 |
| NL7211930A (en) | 1973-03-06 |
| DE2243449B2 (en) | 1975-08-14 |
| CA979074A (en) | 1975-12-02 |
| JPS4834452A (en) | 1973-05-18 |
| GB1400471A (en) | 1975-07-16 |
| FR2152055A5 (en) | 1973-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6501829B2 (en) | Grid holding frame, which provides grid information to X-ray image processing apparatus | |
| EP0314506A1 (en) | Method and apparatus for improving the alignment of radiographic images | |
| JP2001351091A (en) | Image processing apparatus, photographing apparatus, image processing system, image processing method, and storage medium | |
| CN113397577A (en) | Dynamic quality management device, dynamic quality management program, and dynamic quality management method | |
| JPH10127622A (en) | X-ray computed tomography apparatus | |
| EP1080690B1 (en) | X-ray camera | |
| US6801647B1 (en) | Method and apparatus for evaluating digital radiation image | |
| US6829389B1 (en) | Method and apparatus for evaluating image | |
| GB1400471A (en) | Apparatus for automatic computation of cardiothoracic ratio | |
| GB2096440A (en) | Radiography apparatus including a photographic film memory | |
| WO1986001912A1 (en) | Crt multi-frame photographing apparatus | |
| JPH07100060B2 (en) | Image measurement method to obtain lung contour lines | |
| JP2019005013A (en) | Analysis device and analysis system | |
| US5710841A (en) | Method for improving the contrast in photographic film materials | |
| JPH07194583A (en) | X-ray diagnostic device | |
| JP2002112985A (en) | Image diagnosis support device | |
| JPS61231476A (en) | Apparatus for converting data | |
| JPH07148143A (en) | Bone measurement method and device | |
| JPH0777749A (en) | Radiation image information reader | |
| RU2007120C1 (en) | Method roentgenographic investigating of biological objects | |
| RU177536U1 (en) | DIGITAL X-RAY FLUOROGRAPHIC CAMERA | |
| JP2644996B2 (en) | Nuclear medicine data processing equipment | |
| JPH03215256A (en) | Bone evaluation apparatus | |
| JP2771452B2 (en) | Bone measurement method and device therefor | |
| JP2732982B2 (en) | Image reading method and apparatus |