Publication number | US20060002602 A1 |

Publication type | Application |

Application number | US 11/207,775 |

Publication date | Jan 5, 2006 |

Filing date | Aug 22, 2005 |

Priority date | Apr 7, 1998 |

Also published as | EP0964362A1, EP1672551A2, EP1672551A3, US6853740, US7046836, US20050094860 |

Publication number | 11207775, 207775, US 2006/0002602 A1, US 2006/002602 A1, US 20060002602 A1, US 20060002602A1, US 2006002602 A1, US 2006002602A1, US-A1-20060002602, US-A1-2006002602, US2006/0002602A1, US2006/002602A1, US20060002602 A1, US20060002602A1, US2006002602 A1, US2006002602A1 |

Inventors | Hiroyuki Shinbata |

Original Assignee | Canon Kabushiki Kaisha |

Export Citation | BiBTeX, EndNote, RefMan |

Patent Citations (51), Referenced by (2), Classifications (15) | |

External Links: USPTO, USPTO Assignment, Espacenet | |

US 20060002602 A1

Abstract

An image processing method is intended to extract an irradiation area in an image easily and accurately.

A calculation unit **202 **calculates secondary difference values of density values representing respective areas in three areas of a rectangular shape arranged in parallel on an image, determined by a calculation area input unit **200 **and a calculation area determination unit **201, **between adjacent rectangular areas and the secondary difference values thus calculated are stored in a memory unit **203**. Then a judgment unit **204 **judges one end point of an irradiation area from the secondary difference values thus stored and an irradiation area determination unit **205 **determines the irradiation area, based on a plurality of end points thus judged.

Claims(14)

means for determining a plurality of areas arranged in a predetermined direction on an image and each having a predetermined shape;

means for calculating a secondary difference value of density values representing the respective areas in said plurality of areas; and

means for judging one end point of an irradiation area from said secondary difference values calculated by said calculating means.

a step of detecting an irradiation end, based on a density distribution in each area, for a plurality of areas in a desired direction in said image; and

a step of evaluating the result of said detection, based on the result of irradiation ends detected for each of said plurality of areas.

a first calculation step of calculating first values representing a gradient of radiographic image data in a direction on the radiographic image data;

a second calculation step of calculating second values representing a degree of change of the gradient of the radiographic image data based on the first values; and

a judgment step of judging an edge point of an irradiation area from the second values.

first calculation means which calculate first values representing a gradient of a radiographic image data in a direction on the radiation image data;

second calculation means which calculate second values representing a degree of change of the gradient of the radiographic image data based on the first values; and

judgment means which judge an edge point of an irradiation area from the secondary values.

a first calculation step of calculating first values representing a gradient of a radiographic image data in a direction on the radiographic image data;

a second calculation step of calculating second values representing a degree of change of the gradient of the radiographic image data based on the first values; and

a judgment step of judging an edge point of an irradiation area from the secondary values.

a first calculation step of calculating first values representing a gradient of a radiographic image data in a direction on the radiographic image data;

a second calculation step of calculating second values representing a degree of change of the gradient of the radiographic image data based on the first values;

a detecting step of calculating coordinates for edge points of an irradiation area from the second values;

a calculating step of calculating an evaluation value for evaluating positional relations among candidates; and

a judging step of judging whether radiographing is performed by an imaging device having an irradiation diaphragm function in a state of irradiation diaphragm or in a state of no irradiation diaphragm, based on the evaluation value.

a first calculation step of calculating first values representing a gradient of a radiographic image data in a direction on the radiation image data;

a second calculation step of calculating second values representing a degree of change of the gradient of the radiographic image data based on the first values;

a detecting step of calculating coordinates for edge points of an irradiation area from the second values;

a calculating step of calculating an evaluation value for evaluating positional relations among candidates; and

a judging step of judging whether radiographing is performed by an imaging device having an irradiation diaphragm function in a state of irradiation diaphragm or in a state of no irradiation diaphragm, based on the evaluation value.

first calculation means which calculate first values representing a gradient of a radiographic image data in a direction on the radiation image data;

second calculation means which calculate second values representing a degree of change of the gradient of the radiographic image data based on the first values;

detecting means which calculate coordinates for edge points of an irradiation area from the second values;

calculating means which calculate an evaluation value for evaluating positional relations among candidates; and

judging means which judge whether radiographing is performed by an imaging device having an irradiation diaphragm function in a state of irradiation diaphragm or in a state of no irradiation diaphragm, based on the evaluation value.

Description

- [0001]1. Field of the Invention
- [0002]The present invention relates to an image processing method, an image processing apparatus, and a storage medium using secondary differences for recognition of an irradiation field.
- [0003]2. Related Background Art
- [0004]The recent progress in the digital technology makes it common to convert a radiographic image to a digital image signal, subject this digital image signal to image processing, and display it on a CRT or the like or print it out. In photography of the radiographic image it is common practice to diaphragm the irradiation area so as to radiate the radiation only into a necessary area of a human body on the humane grounds or for the reason of preventing scattering from unnecessary areas so as to prevent lowering of contrast. Normally, a processing parameter is determined from distribution of density values of the image prior to the image processing and the image processing is then carried out based on the parameter thus determined. However, if the irradiation area is not limited, unwanted information, in a sense, outside an area of interest is used for the determination of the image processing parameter, which poses a problem that appropriate image processing is not effected.
- [0005]It is thus necessary to extract the irradiation area and determine the image processing parameter from information of only the area of interest. As an extraction method of the irradiation area, for example, Japanese Patent Application (Laid-Open) No. 05-007579 describes a method for dividing the image area into small areas and extracting the irradiation area, based on values of variance in the small areas. Further, for example, Japanese Patent Publication No. 6-90412 describes a method for expressing a change of image densities between a predetermined number of pixels near an end portion of the image area by an approximate equation substantially of a linear equation and extracting the irradiation area, based on the difference between assumed image density values according to this approximate equation and actual density values indicated by a sample image signal.
- [0006]These methods are predicated on the premise that an object area in the photographic image includes the irradiation area (an area actually irradiated), and a process of judging whether the object area is an area including the irradiation area (which will also be referred to as an “area with an irradiation diaphragm”) or an area not including the irradiation area (which will also referred to as an “area without the irradiation diaphragm”) is carried out as a preliminary operation before execution of these methods.
- [0007]Methods for judging presence/absence of the irradiation diaphragm include methods for comparing an average, a median, or the like of densities in a central area of the image with an average of densities of the object area and judging that the object area is an area with the irradiation diaphragm, if the average of densities of the object area is not more than a predetermined value, as described in U.S. Pat. No. 5,091,970, for example.
- [0008]The method of above Japanese Patent Application (Laid-Open) No. 05-007579 had to obtain the density variance in each small area and thus had the problem of increased computation complexity and computation time. Even within the area of interest, density values vary suddenly, for example, at the edge of the lang field and a change rate of densities could become higher there than at the edge of the irradiation area. Particularly, the variances become high in the area where the lang is in contact with the ribs at the edge of the lang. This caused candidate points for the irradiation area to be extracted from outside the edge of the irradiation area as well and this posed another problem that the judgment became difficult.
- [0009]In the method of above Japanese Patent Publication No. 6-90412, a foot area outside the irradiation area was assumed for calculation of the linear approximate equation and there arose the problem that the method was invalid if this foot area was not able to be extracted well. Further, it was also assumed that the foot area outside the irradiation area included two areas of a gentle portion and a quickly rising portion and this posed another problem that the method was invalid where the entire foot area was gentle or where there were three or more changes of density gradient.
- [0010]The present invention has been accomplished in order to solve the problems described above and an object of the invention is, therefore, to extract the irradiation area accurately.
- [0011]An image processing method of the present invention is an image processing method comprising:
- [0012]a step of determining a plurality of areas arranged in a predetermined direction on an image and each having a predetermined shape;
- [0013]a step of calculating a secondary difference value of density values representing the respective areas in the plurality of areas; and
- [0014]a step of judging one end point of an irradiation area from the secondary difference values calculated in the calculating step.
- [0015]In the conventional image judgment method as described in above U.S. Pat. No. 5,091,970 etc., for example, where the photographic image is one in a state in which the subject overlaps with the edge of the image, the average density value at the end of the image varies depending upon the area of the subject over the edge of the image and upon the transmittance of the radiation. Therefore, the presence/absence of the irradiation diaphragm was sometimes misjudged in the area at the end of the image.
- [0016]When the radiant dose of the irradiation diaphragm was small, the density difference also became small between the central part of the image and the object area and the presence/absence of the irradiation diaphragm in the object area was misjudged in some cases.
- [0017]The present invention has been accomplished to eliminate the drawbacks described above and an object of the invention is, therefore, to make it possible to accurately judge whether an image area as a processed object is an area including the irradiation area.
- [0018]An image processing method of the present invention is an image processing method for extracting an irradiation area in an input image, the image processing method comprising:
- [0019]a step of detecting an irradiation end, based on a density distribution in each area, for a plurality of areas in a desired direction in the image; and
- [0020]a step of evaluating the result of the detection, based on the result of irradiation ends detected for each of the plurality of areas.
- [0021]Another image processing method of the present invention is an image processing method for judging whether an object area in an image includes an irradiation area, the method comprising:
- [0022]a secondary difference value acquisition step of acquiring secondary difference values from one-dimensional image data of the object area;
- [0023]an irradiation end extraction step of extracting a coordinate of an end point of the irradiation area from the secondary difference values acquired in the secondary difference value acquisition step;
- [0024]a comparison step of comparing the coordinate extracted in the irradiation end extraction step with a coordinate of an end point of the irradiation area included in the image, the coordinate being obtained preliminarily; and
- [0025]a judgment step of judging whether the object area includes the irradiation area, based on the result of the comparison in the comparison step.
- [0026]The other objects, arrangements, and effects of the present invention will become more apparent with the description of the embodiments of the present invention which follows.
- [0027]
FIG. 1 is a block diagram of an irradiation area extraction device according to Embodiment 1; - [0028]
FIG. 2 is a block diagram to show the structure inside an irradiation area extraction unit; - [0029]
FIG. 3 is a flow chart to show processing procedures in the irradiation area extraction unit; - [0030]
FIG. 4 is a structural diagram to show calculation areas; - [0031]
FIG. 5A andFIG. 5B are structural diagrams to show a radiographic image and density values of one line on the radiographic image; - [0032]
FIG. 6A andFIG. 6B are characteristic diagrams to show density values of one line onFIG. 5A and values of secondary difference on this line; - [0033]
FIG. 7 is a block diagram to show the structure of an image judgment device of Embodiment 2-1; - [0034]
FIG. 8 is a diagram for explaining an example of image data inputted into the above image judgment device; - [0035]
FIG. 9 is a flow chart for explaining a processing program carried out in the above image judgment device; - [0036]
FIG. 10 is a block diagram to show the structure of an image judgment device of Embodiment 2-2; - [0037]
FIG. 11 is a flow chart for explaining a processing program carried out in the above image judgment device; - [0038]
FIG. 12 is a block diagram to show the structure of an image judgment device of Embodiment 3-1; - [0039]
FIG. 13 is a diagram for explaining an example of image data inputted into the above image judgment device; - [0040]
FIG. 14 is a flow chart for explaining a processing program carried out in the above image judgment device; - [0041]
FIG. 15 is a block diagram to show the structure of an angle extraction device and an area extraction device in Embodiment 4-1; - [0042]
FIG. 16 is a flow chart of a processing procedure sequence in the angle extraction device and the area extraction device; - [0043]
FIG. 17 is a diagram for explaining the processing in the angle extraction device and the area extraction device; - [0044]
FIG. 18 is a flow chart of a processing procedure sequence for obtaining an approximate center in an area in Embodiment 4-4; - [0045]
FIG. 19 is a block diagram to show the structure of an area extraction device in Embodiment 4-7; - [0046]
FIG. 20 is a flow chart of a processing procedure sequence in the area extraction device of Embodiment 4-7; - [0047]
FIG. 21 is a block diagram to show the structure of an area extraction device in Embodiment 4-9; and - [0048]
FIG. 22 is a flow chart of a processing procedure sequence in the area extraction device of Embodiment 4-9. - [0049]
FIG. 1 is a block diagram to show the structure of an irradiation area extraction device according to an embodiment of the present invention. In the figure, reference numeral**101**designates an image input unit for accepting input of an image from the outside,**102**an irradiation area extraction unit for extracting an irradiation area irradiated with radiation from the image inputted at the image input unit**101**,**103**an image processing unit for performing image processing based on the irradiation area extracted at the irradiation area extraction unit**102**, and**104**a display unit for displaying an image subjected to the processing at the image processing unit**103**. - [0050]
FIG. 2 is a block diagram to show the inside structure of the irradiation area extraction unit**102**, in which reference numeral**200**designates a calculation area input unit for accepting input of a direction, a start point, and an end point for determination of calculation areas and**201**a calculation area determination unit for determining, based on the input information at the calculation area input unit**200**, areas of calculation carried out by a calculation unit**202**from the input image supplied from the image input unit**101**. Numeral**202**denotes a calculation unit for calculating primary difference values and a secondary difference value, described hereinafter, from the calculation areas determined at the calculation area determination unit**201**,**204**a judgment unit for judging irradiation area ends,**203**a memory unit for storing the values calculated at the calculation unit**202**and the irradiation area ends judged at the judgment unit**204**, and**205**an irradiation area determination unit for determining the irradiation area from the irradiation area ends stored in the memory unit**203**and judged at the judgment unit**204**. - [0051]
FIG. 3 is a flow chart to show the flow of the processing at the irradiation area extraction unit**102**.FIG. 4 is a diagram to show the calculation areas determined at the calculation area determination unit**201**.FIG. 5A is a diagram to show a radiographic image andFIG. 5B is a diagram to show density values on line X_{0 }to X_{3 }ofFIG. 5A . The abscissa indicates positions on the line X_{0 }to X_{3 }and the ordinate density values on the line.FIG. 6A is an enlarged view of part ofFIG. 5B .FIG. 6B is a plot of secondary difference values calculated at the calculation unit**202**against points on the above line between X_{0 }and X_{3}. - [0052]The operation will be described below.
- [0053]In
FIG. 1 , an image from an external device not illustrated is inputted through the image input unit**101**into the irradiation area extraction unit**102**. The irradiation area extraction unit**102**extracts the irradiation area from the input image and delivers the information of the irradiation area to the image processing unit**103**. The image processing unit**103**performs the image processing based on the information of the irradiation area extracted at the irradiation area extraction unit**102**and the image through the image processing is displayed on the display unit**104**. - [0054]The processing at the irradiation area extraction unit
**102**ofFIG. 2 will be described next according to the flow ofFIG. 3 . The irradiation area extraction unit**102**receives the input of the calculation direction, the calculation start point, and the calculation end point necessary for the determination of the calculation areas at the calculation area input unit**200**(step S**301**). It is, however, noted that the input herein is required only for an initial input screen and it does not have to be set for the next input screen and after. Next, the calculation area determination unit**201**determines the calculation areas, based on the information from the calculation area input unit**200**. The calculation areas herein are three areas A, B, and C of a rectangular shape arranged in parallel on the image, as illustrated inFIG. 4 , and are used as calculation areas for calculation of the secondary difference value at a calculation point indicated by (x,y). These three areas A, B, and C have an equal area which is determined by two parameters of “a” and “b” and the distance between the areas is represented by “d”. These parameters “a”, “b”, “d” are determined experimentally. The calculation area determination unit**201**preliminarily determines the calculation areas for all calculation points and the result, together with information of duplicate calculation areas, is stored in the memory unit**203**(step S**302**). - [0055]Then the calculation unit
**202**calculates representatives of density values in the above calculation areas A, B, and C according to a method described hereinafter and further calculates the primary difference values and the secondary difference value. The primary difference values are two values, which are a value “e” resulting from subtraction of the density representative in the area B from the density representative in the area C and a value “f” resulting from subtraction of the density representative in the area A from the density representative in the area B. The secondary difference value is a value resulting from subtraction of the value “f” from the value “e”. The primary difference values and secondary difference value calculated in this way are stored in the memory unit**203**. From the duplicate information of the calculation areas stored in the memory unit**203**, the primary difference values for duplicate calculation areas are obtained from the values stored in the memory unit**203**and the secondary difference value thereof is calculated using them (steps S**303**and S**304**). - [0056]Then the judgment unit
**204**judges the irradiation area ends from the primary difference values and the secondary difference values stored in the memory unit**203**.FIG. 6B shows the secondary difference values calculated at the calculation unit**202**. The secondary difference values are negatively large at the points X_{1 }and X_{2 }where the densities vary suddenly, and the aforementioned primary difference values f take positive values in a density increasing direction (when seen from the side of the point X_{0}) but take negative values in a density decreasing direction. From this property the judgment unit**204**judges that a point where the above secondary difference value is minimum and the primary difference values are positive is a candidate for an end point of the irradiation area. If there are plural candidates, the judgment unit determines that a first appearing candidate is an end point of the irradiation area (step S**305**). Then the irradiation area end point is stored in the memory unit**203**(step S**306**). - [0057]Further, the judgment unit
**204**judges whether the processing is completed for all the directions supplied from the calculation area input unit**200**and, if not, the processing is repeated from step S**302**(step S**307**). After the irradiation area end points are determined in all the directions, the irradiation area determination unit**205**determines the irradiation area. Here, a line passing the end point of the irradiation area and being perpendicular to the calculation direction of the secondary difference value (for example, the direction from X_{0 }to X_{3 }inFIGS. 5A and 5B ) is calculated for all the end points of the irradiation area and an area surrounded by these lines obtained is determined as an irradiation area (step S**308**). - [0058]The present embodiment presents the effects of less computational complexity and capability of readily extracting the end point of the irradiation area, because the calculation unit
**202**is arranged to calculate the representatives of density values in the calculation areas A, B, C and calculate the secondary difference value based on the representative values. - [0059]Since the device of the present embodiment has the irradiation area determination unit
**205**for determining the irradiation area, the present embodiment has such an effect that the irradiation area can be extracted even if the irradiation area is rectangular, polygonal, or circular. The present embodiment further has such an effect that when the irradiation area is polygonal, an almost circumscribed area can be extracted by selecting directions of lower dimensions than it. - [0060]Described below are methods for calculating the representative value in each of the above calculation areas A, B, C. Simple methods for such calculation include a method for calculating an average of all density values (all pixel values) in a calculation area, a method for sorting all the density values in a calculation area according to the density values and using a density value of a pixel at an intermediate position, i.e., a median, and so on. There are also methods for using the limited number of pixels in each calculation area, instead of all the pixels, by sampling or the like according to the necessity and using an average of density values of the limited number of pixels or a median of the limited number of pixels.
- [0061]Here, when the method for calculating the representative of density values in the calculation area is the method using the average of all the density values in the calculation area, the method has the effect of being resistant to noise. When the representative is determined by the method for sorting all the density values in the calculation area according to the density values and using the density value of the pixel at the intermediate position, the method has the effects of simple calculation and being resistant to noise. Further, when the representative is determined by the method using the average of the limited number of points in the calculation area or the method for sorting the density values of the limited number of points and using the density value at the intermediate position as a representative value, the methods have the effect of capability of further decreasing the computational complexity.
- [0062]Another approach for determining a representative of each calculation area A, B, C shown in
FIG. 4 will be described below. First, let us assume that there are the calculation areas having the vertical height of “a” and the horizontal height of “b”, the distance between the areas is represented by “d”, and each calculation point is indicated by (x,y), as illustrated inFIG. 4 . - [0063]In the present embodiment, where the representative value of each area A, B, or C is S(A), S(B), or S(C), respectively, the secondary difference value SS(X) is calculated according to Eq. (1) below.

*SS*(*X*)=*S*(*A*)−2×*S*(*B*)+*S*(*C*) (1) - [0064]Next described is a method for calculating the representative value according to this approach where the representative of each area is indicated by S(X) in general. First, a pixel value at each calculation point (x,y) is defined as f(x,y) and values F
**1**(*x*) are calculated according to Eq. (2) below.

*F***1**(*x*)=∫_{0}^{a}*f*(*x,y*)*dy*(2) - [0065]Here, the integration range is the vertical height “a” of the above rectangular region and for simplicity of the description thereof Eq. (2) is expressed so as to perform the integration in the range of “0” to “a”.
- [0066]Further, values F
**2**(*x*) are next calculated according to Eq. (3) below.

*F***2**(*x*)=min{*F***1**(*x+x***1**)−*h*(*x***1**)|*x***1**ε*K}*(3) - [0067]Then S(x) is calculated using the values F
**2**(*x*), according to Eq. (4) below.

*S*(*X*)=max{*F***2**(*x−x***1**)+*h*(*x***1**)|*x***1**εK} (4) - [0068]Here, h(x) is a function represented by Eq. (5) below and K is a domain of definition thereof.
$\begin{array}{cc}\begin{array}{cc}h\left(x\right)=0,& -b/4\le x\le b/4\\ =-\infty ,& \mathrm{otherwise}\end{array}\hspace{1em}& \left(5\right)\end{array}$ - [0069]Then each of the left and right irradiation ends is obtained from a value of x taking minSS(X) in the left half and the right half of the image. The upper and lower ends can also be obtained in a similar fashion.
- [0070]In other words, the above approach is nothing but an operation in which projection of pixel values is made with respect to a predetermined direction (the vertical direction) of the area, values obtained by the projection are smoothed using a one-dimensional morphology filter, and the secondary difference value is calculated with the distance d.
- [0071]With this calculation method, the irradiation area can be detected with less errors than in the case of the above methods using the average and median.
- [0072]The functional blocks
**101**to**104**,**200**to**205**inFIG. 1 andFIG. 2 may be configured on a hardware basis or may be configured in a microcomputer system comprised of CPU, memory, and so on. When they are configured in the microcomputer system, the above memory comprises the memory medium according to the present invention and this memory medium stores a program for carrying out the processing illustrated in the flow chart ofFIG. 3 . This memory medium can be selected from semiconductor memories such as ROM, RAM, and the like, optical disks, magneto-optical disks, magnetic media, and so on, or it may be substantiated in the form of either one of a CD-ROM, a floppy disk, a magnetic tape, a non-volatile memory card, and so on. - [0073]Since the device according to the present embodiment is arranged to determine the calculation areas comprised of the plural areas of the predetermined shape arranged in the predetermined direction, calculate the secondary difference value of the density values representing the respective areas in the plural areas, and judge one end point of the irradiation area from the secondary difference value thus calculated, the present embodiment presents the effects of capability of decreasing the computational complexity and in turn decreasing the computation time and capability of extracting the end points of the irradiation area with accuracy.
- [0074]Since the device of the present embodiment is arranged to determine the irradiation area from the plural end points of the irradiation area judged, it can enjoy the effect of capability of extracting the irradiation area even if the irradiation area is rectangular, polygonal, or circular. When the irradiation area is polygonal, the device presents the effect of capability of extracting the almost circumscribed area by selecting directions of lower dimensions than it.
- [0075]Since the density value representing each area in the plural areas is an average density value in each area, the device can enjoy the effect of capability of extracting the end points of the irradiation area with accuracy even from an image with noise.
- [0076]The present embodiment is carried out, for example, by an image judgment device
**2100**illustrated inFIG. 7 . This image judgment device**2100**is provided with a judgment unit**2110**, a control unit**2120**for carrying out control of the operation of the judgment unit**2110**, and a program memory**2130**to which the control unit**2120**makes access, as illustrated inFIG. 7 . - [0077]The judgment unit
**2110**is composed of a secondary difference value calculation circuit**2201**for calculating the secondary difference value from data of an object area in an input image, a left end point extraction circuit**2202**for extracting a left end point of an irradiation area included in the object area, based on the secondary difference value calculated at the secondary difference value calculation circuit**2201**, a right end point extraction circuit**2203**for extracting a right end point of the irradiation area included in the object area, based on the secondary difference value calculated at the secondary difference value calculation circuit**2201**, and a diaphragm presence/absence judgment circuit**2204**for judging whether the object area is an area with an irradiation diaphragm or an area without an irradiation diaphragm, from the left end point extracted at the left end point extraction circuit**2202**and the right end point extracted at the right end point extraction circuit**2203**. - [0078]Here, for example, an image G as illustrated in
FIG. 8 is supplied as an input image to the judgment unit**2110**. This input image G is a two-dimensional, radiographic, thoracic part image obtained by photographing the thoracic part with the irradiation diaphragm by use of an imaging device with the irradiation diaphragm function. - [0079]In above
FIG. 8 , “X_{a}” and “X_{b}” represent positions of end points (irradiation area ends) of the irradiation area preliminarily obtained with respect to the horizontal axis H, and “A” to “D” are areas at the edges of the image. - [0080]These areas A to D are object areas for which whether the irradiation diaphragm is present or absent (i.e., whether the irradiation area is present or absent) is judged. In this example, the areas A and B are areas with the irradiation diaphragm while the areas C and D are areas without the irradiation diaphragm.
- [0081]Various processing programs for controlling the operation of the judgment unit
**2110**are preliminarily stored in the program memory**2130**. - [0082]Specifically, for example, a processing program according to the flow chart as illustrated in
FIG. 9 is preliminarily stored in the program memory**2130**and this processing program is read and executed by the control unit**2120**, so as to operate the judgment unit**2110**as follows. - [0083]First, when the area C is an object area, the secondary difference value calculation circuit
**2201**calculates the secondary difference values SS(x) from data of the area C according to Eq. (21) below (step S**2301**).

*SS*(*x*)=*f*(*x−d*)−2×*f*(*x*)+*f*(*x+d*) (21) - [0084]In this Eq. (21), “f(x)” indicates image data of a one-dimensional line which crosses the area C in the horizontal direction, “x” coordinate values thereof, and “d” a constant indicating a difference distance.
- [0085]Next, using the secondary difference values SS(x) calculated at the secondary difference value calculation circuit
**2201**, the left end point extraction circuit**2202**extracts the left end point x**1**of the irradiation area included in the area C according to Eq. (22) below (step S**2302**).

*SS*(*x***1**)=min{*SS*(*x*)|0*≦x≦{overscore (x)}}*(22) - [0086]{overscore (x)}: coordinates on the horizontal axis between the left and right ends of the object area (coordinates between xa and xb)
- [0087]On the other hand, using the secondary difference values SS(x) calculated at the secondary difference value calculation circuit
**2201**, the right end point extraction circuit**2203**extracts the right end point x**2**of the irradiation area included in the area C according to Eq. (23) below (step S**2303**).

*SS*(*x***2**)=min{*SS*(*x*)|*{overscore (x)}≦x*≦Length} (23) - [0088]In this Eq. (23), “Length” indicates the length of the area C along the horizontal axis.
- [0089]Then the diaphragm presence/absence judgment circuit
**2204**obtains orders of approximation of the left end point x**1**extracted at the left end point extraction circuit**2202**and the right end point x**2**extracted at the right end point extraction circuit**2203**with respect to the left and right irradiation end points preliminarily obtained (the positions X_{a }and X_{b }of the irradiation ends along the horizontal axis H) and, if the orders of approximation are high, the diaphragm presence/absence judgment circuit**2204**judges that the left end point x**1**and the right end point x**2**are the left and right end points of the irradiation area (the area irradiated directly), i.e., that the area C is an area without the irradiation diaphragm. On the other hand, if the orders of approximation are low, it is judged that the area C is an area with the irradiation diaphragm (step S**2304**). - [0090]Specifically, using a constant “e” indicating an approximation width, if the below conditions are satisfied the area C is judged as an area without the irradiation diaphragm; otherwise the area C is judged as an area with the irradiation diaphragm.

*−e≦x***1**−*X*_{a}*≦e*and

*−e≦x***2**−*X*_{b}*≦e* - [0091]When in an object area the irradiation ends are detected similarly to the irradiation ends in the adjacent area, the object area is judged as an area irradiated with radiation as the adjacent area was.
- [0092]Steps S
**2301**to S**2304**are carried out for each of the other areas A, B, D, in a similar fashion to the area C. - [0093]As described above, the present embodiment is arranged to calculate the secondary difference values SS(x) from the data of the object area subjected to the judgment of presence/absence of the radiation diaphragm in the two-dimensional input image G and judge whether the object area is an area with or without the irradiation diaphragm, using the secondary difference values SS(x).
- [0094]Construction of the device in this structure permits boundary points between the directly irradiated area and the other areas to be extracted with accuracy even from photographic image obtained when a subject with low transmittances of radiation is photographed. Therefore, presence or absence of the irradiation diaphragm can be judged with accuracy for the object areas including the irradiation area. The presence or absence of the irradiation diaphragm can also be judged for an object area in a photographic image in which a portion with low radiation transmittances, such as the abdominal part or the like, overlaps with the edge of the image.
- [0095]It may also be contemplated that on the occasion of extracting the left end point x
**1**and the right end point x2 in the object area at the left end point extraction circuit**2202**and at the right end point extraction circuit**2203**, another condition for the secondary difference values SS(x) used in the extraction, for example such a condition that the secondary difference values SS(x) are not more than a fixed threshold, is added. - [0096]The present embodiment is arranged to apply, for example, the projection in the object area represented by Eq. (24) below to Eq. (21) in above Embodiment 2-1.

*f*(*x*)=∫_{b}^{c}*f*(*x,y*)*dy*(24) - [0097]In this Eq. (24), “f(x,y)” represents the image data of the object area in the input image G, and “x” and “y” are coordinates on the horizontal axis and on the vertical axis, respectively. Further, “b” and “c” indicate the object area.
- [0098]When the projection in the object area represented by above Eq. (24) is used in above Eq. (21) in Embodiment 2-1, it can achieve the effects similar to those in the case where steps S
**2301**to S**2304**described above are carried out in the averaged state of the data of the object area. Namely, misjudgment of presence or absence of the irradiation diaphragm can be prevented from being caused by alteration of coordinates of the object area due to influence of scattered rays, noise, and so on. Therefore, presence/absence of the irradiation diaphragm in the object area can be judged with better accuracy. - [0099]The present embodiment is applied, for example, to an image judgment device
**2400**as illustrated inFIG. 10 . - [0100]This image judgment device
**2400**has the structure similar to that of the image judgment device**2100**in Embodiment 2-1 described above, but is different in the internal structure of the judgment unit**2110**therefrom. - [0101]In the image judgment device
**2400**ofFIG. 10 , elements operating similarly to those of the image judgment device**2100**ofFIG. 7 are denoted by the same reference symbols and the detailed description thereof is omitted herein. Only the different structure from Embodiment 2-1 will be described in detail herein. - [0102]The judgment unit
**2140**is composed of, as illustrated inFIG. 10 , a coordinate indication circuit**2404**, a secondary difference value calculation circuit**2401**for calculating the secondary difference values from data of an object area (either one of the areas A to D in the input image G ofFIG. 8 in this example) in the input image according to coordinates indicated by the coordinate indication circuit**2404**, a left end point extraction circuit**2402**for extracting a left end point of the irradiation area included in the object area, based on the secondary difference values calculated at the secondary difference value calculation circuit**2401**, and a right end point extraction circuit**2403**for extracting a right end point of the irradiation area included in the object area, based on the secondary difference values calculated at the secondary difference value calculation circuit**2401**. - [0103]The judgment unit
**2140**further has a memory circuit**2405**for storing coordinates of a plurality of left and right end points extracted for a plurality of lines by the left end point extraction circuit**2402**and the right end point extraction circuit**2403**, an average coordinate calculation circuit**2406**for calculating an average coordinate at each end of the coordinates of the left and right end points stored in the memory circuit**2405**, a first judgment circuit**2407**for judging whether the object area is an area with the irradiation diaphragm or an area without the irradiation area, using the averages of the coordinates of the left and right end points calculated at the average coordinate calculation circuit**2406**, a coordinate variance calculation circuit**2408**for calculating variances of the coordinates of the left and right end points according to the judgment result of the first judgment circuit**2407**, and a second judgment circuit**2409**for again judging whether the object area is an area with the irradiation diaphragm or an area without the irradiation diaphragm, based on the variances calculated at the coordinate variance calculation circuit**2408**. - [0104]The coordinate indication circuit
**2404**is configured to give an indication of a coordinate for calculation of the secondary difference value at the secondary difference value calculation circuit**2401**thereto after the left end point extraction circuit**2402**and the right end point extraction circuit**2403**extract the left and right end points. - [0105]A processing program, for example, according to the flow chart as illustrated in
FIG. 11 , for controlling the operation of the judgment unit**2140**is preliminarily stored in the program memory**2130**and this processing program is read and executed by the control unit**2120**, whereupon the judgment unit**2140**operates as follows. - [0106]First, supposing the area C is an object area, the secondary difference value calculation circuit
**2401**calculates the secondary difference values SSi(x) from the data of the area C according to Eq. (25) below for each coordinate i indicated by the coordinate indication circuit**2404**(step S**2501**).

*SSi*(*x*)=*fi*(*x−d*)−2*×fi*(*x*)+*fi*(*x+d*) (25) - [0107]In this Eq. (25), “fi(x)” represents image data of a one-dimensional line crossing the area C in the horizontal direction, “x” coordinates thereof, and “i” a coordinate of the horizontal line (a coordinate of the one-dimensional data) indicated by the coordinate indication circuit
**2404**. Further, “d” is a constant indicating a difference distance. - [0108]Using the secondary difference values SSi(x) calculated at the secondary difference value calculation circuit
**2401**, the left end point extraction circuit**2402**then extracts the left end point xLi according to Eq. (26) below (step S**2502**).

*SSi*(*xLi*)=min{*SSi*(*x*)|0*≦x≦{overscore (x)}}*(26) - [0109]{overscore (x)}: coordinates on the horizontal axis between the left and right ends of the object area (coordinates between xa and xb)
- [0110]Then the memory circuit
**2405**stores the left end point xLi extracted at the left end point extraction circuit**2402**(step S**2503**). - [0111]Using the secondary difference values SSi(x) calculated at the secondary difference value calculation circuit
**2401**, the right end point extraction circuit**2403**also extracts the right end point xRi according to Eq. (27) below (step S**2504**).

*SSi*(*xRi*)=min{*SSi*(*x*)|*{overscore (x)}≦x*≦Length} (27) - [0112]In this Eq. (27), “Length” indicates the length of the area C along the horizontal axis.
- [0113]Then the memory circuit
**2405**stores the right end point xRi extracted at the right end point extraction circuit**2403**(step S**2505**). - [0114]After completion of the storage of the left end point xLi and the right end point xRi in the memory circuit
**2405**, the coordinate indication circuit**2404**gives another instruction of a coordinate i of a new horizontal line to the secondary difference value calculation circuit**2401**(step S**2506**). - [0115]According to this indication, the processing from step S
**2501**is carried out again. This loop processing is executed before the coordinate indication circuit**2404**gives an indication of the end to the secondary difference value calculation circuit**2401**. - [0116]After steps S
**2501**to S**2506**are carried out repeatedly and the processing is terminated according to the end indication from the coordinate indication circuit**2404**, the memory circuit**2405**is in a storage state of the left end points xLi and the right end points xRi corresponding to the coordinates i indicated during the processing by the coordinate indication circuit**2404**to the secondary difference value calculation circuit**2401**. - [0117]Then the average coordinate calculation circuit
**2406**calculates the averages of the left end points xLi and the right end points xRi stored in the memory circuit**2405**(step S**2507**). - [0118]This obtains the average of the left end points xLi (left end point average) XL and the average of the right end points xRi (right end point average) XR.
- [0119]Next, similar to the irradiation diaphragm judgment circuit
**2204**ofFIG. 7 , the first judgment circuit**2407**compares the left end point average XL and the right end point average XR calculated at the average coordinate calculation circuit**2406**with the left and right irradiation end points preliminarily obtained (the positions X_{a }and X_{b }of the irradiation ends on the horizontal axis H), and if the following relations are satisfied it judges that the left and right end points indicated by the left end point average XL and the right end point average XR are the left and right end points of the irradiation area i.e., that the area C is an area with the irradiation area (or an area without the irradiation diaphragm).

−*e≦XL−X*_{a}*≦e*and

−*e≦XR−X*_{b}*≦e*

On the other hand, if either one is not met, the judgment circuit judges that the area C is an area without the irradiation area (or an area with the irradiation diaphragm) (step S**2508**). - [0120]The next steps S
**2509**and S**2510**are carried out only if the first judgment circuit**2407**judges that the area is an area without the irradiation diaphragm. - [0121]Namely, the coordinate variance calculation circuit
**2408**calculates the variances of the coordinates of the left and right end points, using the left end points xLi and the right end points xRi stored in the memory circuit**2405**and the left end point average XL and the right end point average XR calculated at the average coordinate calculation circuit**2406**(step S**2509**). - [0122]This operation provides the variance VL of the coordinates of the left end points and the variance VR of the coordinates of the right end points.
- [0123]The second judgment circuit
**2409**compares the variances VL and VR of the coordinates of the left and right end points obtained at the coordinate variance calculation circuit**2408**with a predetermined threshold TH. When the following relations are satisfied, the second judgment circuit**2409**judges that the area C is an area with the irradiation area (or an area without the irradiation diaphragm).

VL<TH and

VR<TH

On the other hand, if either one is not met, the second. judgment circuit**2409**judges that the area C is an area without the irradiation area (or an area with the irradiation diaphragm) (step S**2510**). - [0124]Steps S
**2501**to S**2510**are also carried out for each of the other areas A, B, D in the similar fashion to the area C. - [0125]As described above, the present embodiment is arranged to carry out the detection of the irradiation ends on the plural lines in the object area while the coordinate indication circuit
**2404**gives the indications of the coordinates i of one-dimensional data lines crossing the object area in the horizontal direction to the secondary difference value calculation circuit**2401**. - [0126]Since this structure is arranged to carry out the detection of the irradiation ends on the plural lines in the object area, the present embodiment is more unlikely to be affected by the scattered rays and the noise and can perform the judgment with better accuracy than in the case of the detection of the irradiation ends on a certain line.
- [0127]Since the variances of the irradiation ends are also used for the judgment of presence/absence of the irradiation diaphragm by provision of the second judgment circuit
**2409**, for example, when the image end portions are exposed to radiation because of the influence of the scattered rays, so as to result in making the irradiation ends extracted in the image end portions coincident with the irradiation ends preliminarily obtained, the degree of spread of the scattered rays can be quantified and thus misjudgment of presence/absence of the irradiation diaphragm can be prevented for sure. - [0128]The present embodiment is arranged to calculate the secondary difference values SS(x) according to above Eq. (21) and to extract the left end point x
**1**of the irradiation area included in the object area, using the secondary difference values SS(x), for example, similar to Embodiment 2-1 described above. At this time the present embodiment also uses the sign of the primary difference value S(x**1**) given by Eq. (28) below.

*S*(*x***1**)=*f*(*x***1**)−*f*(*x***1**−*d*) (28) - [0129]Specifically, for example, if the sign of the primary difference value S(x
**1**) is “negative” and if above Eq. (22) is satisfied, the left end point x**1**is regarded as a left end point of the irradiation area. - [0130]For extracting the right end point x
**2**of the irradiation area included in the object area, the present embodiment also uses the sign of the primary difference value S(x**2**) given by Eq. (29) below.

*S*(*x***2**)=*f*(*x***2***+d*)−*f*(*x***2**) (29)

For example, if the sign of the primary difference value S(x**2**) is “negative” and if above Eq. (23) is satisfied, the right end point x**2**is regarded as a right end point of the irradiation area. - [0131]Using the left end point x
**1**and the right end point x**2**obtained in this way, the judgment of presence/absence of the irradiation diaphragm in the object area is carried out in the similar fashion to that in Embodiment 2-1 described above. - [0132]When the conditions of the signs of the primary difference values are also added to the detection of the irradiation ends of the object area, presence/absence of the irradiation diaphragm in the object area can be judged with consideration to the inclination of image data outside the irradiation area due to the scattered rays. This permits the presence/absence of the irradiation diaphragm in the object area to be judged without misjudgment and with better accuracy even if there is a portion or the like where densities change quickly in the object area.
- [0133]The present embodiment is arranged first to calculate the secondary difference values from the data of the object area, for example similar to Embodiment 2-1, but at this time, the present embodiment is arranged to calculate the secondary difference values from data obtained after the data of the object area is filtered.
- [0134]Specifically, for example, the image data of a one-dimensional line of an object area is represented by “f(x)”, the data is subjected to a filtering process according to Eqs. (30) and (31) below, and the secondary difference values are calculated from the values F
**2**obtained as a result.$\begin{array}{cc}\mathrm{F1}\left(x\right)=\mathrm{min}\left\{f\left(x+\mathrm{x1}\right)-h\left(\mathrm{x1}\right)\u2758-d\le \mathrm{x1}\le d\right\}& \left(30\right)\\ \mathrm{F2}\left(x\right)=\mathrm{max}\left\{\mathrm{F1}\left(x-\mathrm{x1}\right)+h\left(\mathrm{x1}\right)\u2758-d\le \mathrm{x1}\le d\right\}\text{}h\left(x\right)\{\begin{array}{cc}=0\text{:}& -d\le x\le d\\ =-\infty \text{:}& \mathrm{otherwise}\end{array}& \left(31\right)\end{array}$ - [0135]When Embodiment 2-1 is modified so as to calculate the second difference values after the data of the object area is smoothed by the filtering process as described above, the presence/absence of the irradiation diaphragm in the object area can be judged with better accuracy without influence of the noise, particularly, without influence of the noise on the line.
- [0136]The present embodiment is carried out, for example, by an image judgment device
**3100**as illustrated inFIG. 12 . - [0137]This image judgment device
**3100**is composed of a judgment unit**3110**, a control unit**3120**for controlling the operation of the judgment unit**3110**, and a program memory**3130**to which the control unit**3120**makes access, as illustrated inFIG. 12 . - [0138]The judgment unit
**3110**is composed of a coordinate indication circuit**3201**, a characteristic value extraction circuit**3202**for extracting a characteristic value from data of an object area in an input image according to a coordinate indicated by the coordinate indication circuit**3201**, an end point extraction circuit**3203**for extracting a coordinate of an end point of the irradiation area included in the object area, based on the characteristic value extracted at the characteristic value extraction circuit**3202**, a memory circuit**3204**for storing coordinates of end points extracted at the end point extraction circuit**3203**, and a judgment circuit**3205**for judging whether the object area is an area with the irradiation diaphragm or an area without the irradiation diaphragm, from the coordinates of the end points stored in the memory circuit**3204**. - [0139]Here, for example, an image G illustrated in
FIG. 13 is supplied as an input image to the judgment unit**3110**. This input image G is a two-dimensional, radiographic, thoracic part image obtained by photographing the thoracic part with the irradiation diaphragm by use of an imaging device with the irradiation diaphragm function. - [0140]In above
FIG. 13 , “X_{a}” and “X_{b}” represent positions of irradiation area ends on the horizontal axis X. “A” to “D” are areas at the edges of the image part. In this example, the areas A and B out of these “A” to “D” are areas with the irradiation diaphragm while the other areas C and D are areas without the irradiation diaphragm. - [0141]Various processing programs for controlling the operation of the judgment unit
**3110**are preliminarily stored in the program memory**3130**. - [0142]Specifically, for example, a processing program according to the flow chart as illustrated in
FIG. 14 is preliminarily stored in the program memory**3130**and this processing program is read and executed by the control unit**3120**to operate the judgment unit**3110**as follows. - [0143]Described below is a process of judging whether the irradiation diaphragm is present or absent at the lower edge.
- [0144]First, the coordinate indication circuit
**3201**gives indications of plural coordinates on the X-axis between X_{a }and X_{b }in the input image G of aboveFIG. 13 to the characteristic value extraction circuit**3202**(step S**3301**). - [0145]Specifically, the coordinate indication circuit
**3201**successively indicates, for example, coordinates i (=1 to 10) on the X-axis, which are ten points distributed at equal intervals between X_{a }and X_{b}. - [0146]According to each coordinate i indicated by the coordinate indication circuit
**3201**, the characteristic value extraction circuit**3202**then calculates a characteristic value of one-dimensional image data corresponding to the coordinate i, for example, the secondary difference values SSi(y) according to Eq. (41) below (step S**3302**).

*SSi*(*y*)=*fi*(*y−d*)−2×*fi*(*y*)+*fi*(*y+d*) (41) - [0147]In this Eq. (41), “fi(y)” indicates image data of a one-dimensional line crossing the object area in the vertical direction, “y” coordinates thereof, and “d” a constant indicating a difference distance.
- [0148]Using the secondary difference values SSi(y) obtained at the characteristic value extraction circuit
**3202**, the end point extraction circuit**3203**then extracts a coordinate yi of an end point (an irradiation end point) of the irradiation area included in the object area according to Eq. (42) below (step S**3303**).

*SS*(*yi*)=min{*SSi*(*y*)|0≦*y*≦Length} (42) - [0149]In this Eq. (42), “Length” represents the length of the input image G along the vertical axis (the direction of the Y-axis).
- [0150]Then the memory circuit
**3204**stores the coordinate yi of the irradiation end point obtained at the end point extraction circuit**3203**(step S**3304**). - [0151]The above processing of steps S
**3301**to S**3304**is carried out repeatedly before the coordinate i indicated by the coordinate indication circuit**3201**reaches “10”. This results in storing the coordinates y**1**, y**2**, y**3**, . . . , y**10**of the irradiation end points corresponding to the coordinates i=1, 2, 3, . . . , 10 indicated by the coordinate indication circuit**3201**, in the memory circuit**3204**. - [0152]After completion of the processing of steps S
**3301**to S**3304**up to the coordinate i=10, the judgment circuit**3205**then calculates an average of the coordinates y**1**, y**2**, y**3**, . . . , y**10**stored in the memory circuit**3204**(step S**3305**) and calculates a variance Bv of the irradiation end points (step S**3306**). - [0153]Then the judgment circuit
**3205**judges that the object area is an area without the irradiation diaphragm, if the variance Bv is not less than a predetermined threshold TH; otherwise, it is judged that the object area is an area with the irradiation diaphragm (step S**3307**). - [0154]The threshold TH is a constant which is determined experimentally.
- [0155]As described above, the present embodiment is arranged to carry out the extraction of irradiation ends of plural column lines in the object area for which presence/absence of the irradiation diaphragm is judged and judge the presence/absence of the irradiation diaphragm of the object area according to the variance of the irradiation ends. If there is an irradiation area in the object area the coordinates of the irradiation end points of the plural column lines will be aligned approximately on one horizontal axis and the variance thereof will be thus small. If there is no irradiation area in the object area the coordinates of the irradiation end points of the plural column lines will be distributed and the variance will be thus large.
- [0156]Therefore, the presence/absence of the irradiation diaphragm in the object area can be determined with accuracy because of the structure using the variance.
- [0157]Since the device is constructed to use the secondary difference value in order to extract the irradiation end, boundary points can be extracted with accuracy between an area irradiated directly and the other areas even in a photographic image obtained by photographing a subject with low transmittances of radiation. Therefore, the device of the present embodiment can judge the presence/absence of the irradiation diaphragm in the object area including the irradiation area with accuracy. In addition, the presence/absence of the irradiation diaphragm in the object area can be judged with accuracy even in a photographic image in which a portion with low radiation transmittances such as the abdominal part or the like overlaps with an end portion of the image.
- [0158]Although the present embodiment is arranged to use the secondary difference value in order to make a judgment of the irradiation diaphragm, the apparatus of the present invention does not always have to be limited to this; for example, where change of density is quick at an irradiation end, the apparatus may also be arranged to use the primary difference values or higher-order difference values. In this case, the primary difference values or the higher-order difference values are obtained from the object area and a first appearing point of a value not less than a predetermined threshold is regarded as a candidate for an irradiation end.
- [0159]Although the present embodiment is arranged to use the variance in order to make a judgment of the irradiation diaphragm, the apparatus of the present invention does not always have to be limited to this; for example, the apparatus may also be arranged to use another index that can indicate discrete degrees of frequency, such as standard deviation.
- [0160]The present embodiment is arranged to apply the projection of a constant width in an object area according to Eq. (43) below to Eq. (41) in Embodiment 3-1 described above.

*fi*(*y*)=∫_{xi−1}^{xi}*f*(*x,y*)*dx*(43) - [0161]In this Eq. (43), “fi(y)” represents image data of a one-dimensional line in the object area in above Eq. (41), and “Xi” and “Xi−1” represent coordinates indicated by the coordinate indication circuit
**3201**. - [0162]When Embodiment 3-1 is modified so as to apply the projection in the object area represented by above Eq. (43) to above Eq. (41) as described above, the modification can enjoy the same effects as in the case where aforementioned steps S
**3301**to S**3307**are carried out in an averaged state of the data of the object area. Namely, this structure can prevent misjudgment about the presence/absence of the irradiation diaphragm from being caused by change of coordinates of the object area due to influence of the scattered rays, the noise, and so on. The presence/absence of the irradiation diaphragm in the object area can be judged with better accuracy accordingly. - [0163]The present embodiment is arranged to calculate the secondary difference values SSi(y) according to above Eq. (41) and extract the irradiation end of the object area using the secondary difference values SSi(y), for example, as Embodiment 3-1 described above was. At this time, the present embodiment also uses the sign of the primary difference value Si(y) expressed by Eq. (44) below.

*Si*(*y*)=*fi*(*y*)−*fi*(*y−d*) (44) - [0164]Specifically, for example, if the sign of the primary difference value Si(y) is “negative” and if above Eq. (42) is satisfied, the point “y” is regarded as an irradiation end.
- [0165]Using the irradiation ends obtained in this way, whether the irradiation diaphragm is present or absent in the object area is judged in the similar fashion to aforementioned Embodiment 3-1.
- [0166]When Embodiment 3-1 is modified to use the sign of the primary difference value as well on the occasion of detecting the irradiation end in the object area, the presence/absence of the irradiation diaphragm in the object area can be judged, also taking the inclination of image data outside the irradiation area due to the scattered rays into consideration. Therefore, whether the irradiation diaphragm is present or absent in the object area can be judged with better accuracy without misjudgment, even if there is a quickly changing portion of density or the like within the object area.
- [0167]The present embodiment is arranged first to calculate the secondary difference values from the data of the object area, similar to Embodiment 3-1; but at this time, the data of the object area is subjected to a filtering process and the secondary difference values are calculated from the data after the filtering process.
- [0168]Specifically, for example, where the image data of a one-dimensional line of the object area is “f(x)”, it is subjected to the filtering process according to Eq. (45) and Eq. (46) below and the secondary difference values are calculated from values F
**2**obtained as a result.$\begin{array}{cc}\mathrm{F1i}\left(y\right)=\mathrm{min}\left\{\mathrm{fi}\left(y+\mathrm{y1}\right)-h\left(\mathrm{y1}\right)\u2758-d\le \mathrm{y1}\le d\right\}& \left(45\right)\\ \mathrm{F2i}\left(y\right)=\mathrm{max}\left\{\mathrm{F1i}\left(y-\mathrm{y1}\right)+h\left(\mathrm{y1}\right)\u2758-d\le \mathrm{y1}\le d\right\}\text{}h\left(x\right)\{\begin{array}{cc}=0\text{:}& -d\le x\le d\\ =-\infty \text{:}& \mathrm{otherwise}\end{array}& \left(46\right)\end{array}$ - [0169]When Embodiment 3-1 is modified so as to calculate the secondary difference values after the data of the object area is smoothed by the filtering process as described above, whether the irradiation diaphragm is present or absent in the object area can be judged with better accuracy without being affected by the noise, particularly, without being affected by the noise on the line.
- [0170]
FIG. 15 is a block diagram to show the structure of an angle extraction device and an area extraction device of two-dimensional image data according to Embodiment 4-1. InFIG. 15 , reference numeral**4101**designates a characteristic quantity calculation means for calculating a characteristic quantity from data of a one-dimensional line at a coordinate indicated by a coordinate indication means**4103**,**4102**an end point extraction means for extracting an area end portion from the characteristic quantity calculated at the characteristic quantity calculation means**4101**,**4103**a coordinate indication means for indicating a coordinate of image data to be calculated by the characteristic quantity calculation means**4101**,**4104**an end point memory means for storing a coordinate of an end point extracted at the end point extraction means**4102**, and**4105**a rotation angle indication means for indicating a rotation angle of a rotation axis onto which coordinates of end points stored in the end point memory means**4104**are projected. - [0171]Reference numeral
**4106**denotes an accumulated quantity calculation means for calculating an accumulated quantity of the coordinates of the end points projected onto the aforementioned rotation axis and stored in the end point memory means**4104**,**4107**an accumulated quantity memory means for storing the accumulated quantity calculated at the accumulated quantity calculation means**4106**,**4108**a rotation angle judgment means for judging a rotation angle from the accumulated quantity stored in the accumulated quantity memory means**4107**, and**4109**an end point judgment means for extracting an area end from the angle judged at the rotation angle judgment means**4108**and the accumulated quantity stored in the accumulated quantity memory means**4107**. - [0172]
FIG. 16 is a flow chart of a processing sequence in the angle extraction device and the area extraction device in Embodiment 4-1.FIG. 17 is a schematic diagram for explaining the processing in the angle extraction device and area extraction device of the present embodiment, in which “a” represents two-dimensional data (data within a light receiving surface in the case of copiers etc., or data in a sensor surface in the case of X-ray apparatus etc.), “b” a rectangular area (an irradiation area in the case of the X-ray apparatus or a sheet surface in the case of the copiers etc.), “c” a start point of the rotation axis, “d” an area end point, “e” a start axis of angle, “f” a rotation axis, and “g” an angle difference θn between the start axis e of angle and the rotation axis f. - [0173]The operation of the present embodiment will be described below according to the flow chart of
FIG. 16 . - [0174]The characteristic quantity calculation means
**4101**ofFIG. 15 calculates the secondary difference values according to Eq. (51) below (steps S**4101**and S**4102**). In the equation fi(x) represents image data of the i-th row line indicated by the coordinate indication means**4103**, SSi(x) the secondary difference values, and “c” a constant.

*SSi*(*x*)=*fi*(*x−c*)−2*×fi*(*x*)+*fi*(*x+c*) (51) - [0175]Then the end point extraction means
**4102**extracts a candidate for an end point of the area according to Eq. (52) below (step S**4103**). - [0176]For the left area end of the i-th row, the following equation is used.

*xi*=min{*SSi*(*x*)|0*≦x≦{overscore (x)}}*(52)

In this equation {overscore (x)} indicates a coordinate on an axis along an indicated direction within the area (the details will be described hereinafter and reference should be made to Eq. (76)). Then the end point memory means**4104**(FIG. 15 ) stores the coordinate of the area end obtained according to Eq. (52) (step S**4104**). Then steps S**4101**to S**4104**are repeated before the processing is complete for all rows indicated by the coordinate indication means**4103**(FIG. 15 ). - [0177]Next, the rotation angle indication means
**4105**(FIG. 15 ) indicates the rotation angle θn (FIG. 17 ) (step S**4105**), and the accumulated quantity calculation means**4106**(FIG. 15 ) performs processing indicated by Eq. (53) to Eq. (57) below to calculate the accumulated quantity L(θn, X) of coordinates x on the aforementioned rotation axis f (FIG. 17 ) (step S**4106**). The projection of the area end points onto the aforementioned rotation axis is obtained from projected values of candidate points onto the rotation axis rotated about the start point c(xc,yc) of the aforementioned rotation axis.

*X*_{iθn}=−(*Xi−xc*)*cos(*θn*)+(*Yi−yc*)*sin(θ*n*) (53)

Here, Xi, Yi represent coordinates of a candidate point for the area end of the i-th row, and X_{iθn }represents projected positions of the candidate points for the area ends on the rotation axis x, which are points resulting from rotation by θn about the start point c (FIG. 17 ) of the aforementioned rotation axis.$\begin{array}{cc}L\left(\theta \text{\hspace{1em}}n,x\right)=\sum _{i=1}^{m}\stackrel{.}{Z}\left(X,{X}_{i\text{\hspace{1em}}\theta \text{\hspace{1em}}n}\right)& \left(54\right)\end{array}$

Here, L(θn,X) indicates the number of projected candidate points on the axis resulting from the rotation of the aforementioned rotation axis by θn.$\begin{array}{cc}\text{\hspace{1em}}\begin{array}{cc}\mathrm{Here},Z\left(X,\mathrm{Xi}\right)=1,& x-\mathrm{d1}\le \mathrm{Xi}\le X+\mathrm{d1}\\ =0,& \mathrm{otherwise},\end{array}\text{\hspace{1em}}& \left(55\right)\end{array}$

where d1 is a constant. - [0178]Then the accumulated quantity memory means
**4107**(FIG. 15 ) stores the accumulated quantity L(θn,x) obtained by Eq. (54) (step S**4107**). Further, the processing of steps S**4105**to S**4107**is repeated for all angles θn indicated by the rotation angle indication means**4105**(FIG. 15 ). - [0179]Then the rotation angle judgment means carries out the processing of Eq. (56) to Eq. (58) below to extract the rotation angle of the area (step S
**4108**).

*n*(θ*n*)=*L*(*θn, X*max−1)+*L*(θ*n, X*max)+*L*(θ*n, X*max+1) (56)

This represents the number of overlaps of area ends at the maximum overlap point in the projection of area ends onto the rotation axis. Here, Xmax indicates a coordinate of X that satisfies the following.

*L*(*θn, X*max)=max{*L*(*θn,X*)|*Xεk*1} (57)

Here, k1 represents coordinates on the rotation axis. - [0180]After all, the rotation angle θ is an angle obtained according to the following equation.

θ=max{*n*(θ*n*)|θ*nεK}*(58)

In this equation K is an arbitrary domain of definition. - [0181]Then the end point judgment means
**4109**(FIG. 15 ) extracts Xmax satisfying Eq. (59) as a candidate for a left end point of the rectangular area from θ determined by Eq. (58) and the accumulated quantities L(θ,X) stored in the accumulated quantity memory means**4107**(Step S**4109**).

*L*(θ,*X*max)=max{*L*(θ,*X*)|*Xεk*1} (59) - [0182]Likewise, candidates are also extracted for the right, upper, and lower end points. It is, however, noted that they do not always have to be extracted, because the angle of the area is already determined.
- [0183]As described above, the present embodiment has the following effect because of the use of the secondary difference value; when the area is extracted, boundary points between the area directly irradiated with X-rays and the other areas can be extracted as end points of the area with accuracy even from an object with low X-ray transmittances.
- [0184]Further, since the angle and end points of the area are determined from the accumulation of end points, the angle and end points of the area can be extracted with accuracy and within short computation time.
- [0185]In the case of the copiers, FAX, OCR, etc., a necessary area can be extracted based on the area and angle obtained. Further, the present embodiment has the effect of capability of automatically carrying out selection of a sheet with accuracy, because the size of the area can be identified. In addition, even if a sheet is placed obliquely, correction can be made and thus the result is equivalent to that where the sheet is placed at a normal position.
- [0186]In the present embodiment, where the image data of a one-dimensional line is f(x) and the secondary difference values thereof are SS(x) defined by Eq. (1), the sign of the primary difference value S(XL) of Eq. (60) below is also added to the operation of extracting the end point XL at the end point extraction means
**4102**(FIG. 15 ).

*S*(*XL*)=*f*(*XL*)−*f*(*XL−d*) (60) - [0187]For example, for extracting the left end point, XL that makes S(XL) negative and that satisfies Eq. (61) below is regarded as a left end point.

*SS*(*XL*)=min{*SS*(*x*)|0*≦x≦{overscore (x)}}*(61) - [0188]As described above, the present embodiment can also take the inclination of the image data outside the area into consideration, because the sign of the primary difference is added to the extraction of the area end. This eliminates erroneous extraction of a quickly changing portion of density in the area and thus the area end point can be extracted with better accuracy.
- [0189]In the present embodiment, where the one-dimensional data of the i-th row is fi(x) and is subjected to the filtering process according to Eqs. (64) and (65) below and values resulting from the filtering process are defined as F
**1***i*(*x*) and F**2***i*(*x*), the characteristic quantity calculation means**4101**(FIG. 15 ) uses F**2***i*(*x*) for calculating the secondary difference values defined by Eq. (51).

*F***1***i*(*x*)=min{*fi*(*x+x***1**)−*h*(*x***1**)|−*d≦x***1***≦d}*(64)

Further, F**2***i*(*x*) is defined as follows.

*F***2***i*(*x*)=max{*F***1***i*(*x−x***1**)+*h*(*x***1**)|−*d≦x***1***≦d}*(65)

Here, h(x) is a function defined below.$\begin{array}{cc}\begin{array}{cc}h\left(x\right)=\text{\hspace{1em}}0,& -d\le \mathrm{x1}\le d\\ =\infty ,& \mathrm{otherwise}\end{array}\hspace{1em}& \left(66\right)\end{array}$ - [0190]Hence, the secondary difference values are calculated as follows.

*SSi*(*x*)=*F***2***i*(*x−c*)−2×*F***2***i*(*x*)+*F***2***i*(*x+c*) (67) - [0191]As described above, the present embodiment is arranged to smooth the one-dimensional image data for calculation of the secondary difference values by the filtering process, thereby accomplishing the effect of being not affected by the noise, particularly by the noise on the line.
- [0192]The present embodiment is an example in which the start point c (
FIG. 17 ) of the aforementioned rotation axis is placed at the barycenter of data not-less than a fixed density. The flow of processing will be described according to a flow chart of a procedure sequence for obtaining an approximate center in the area, shown inFIG. 18 . - [0193]First, an accumulation histogram of whole image data is produced in step S
**4201**and upper densities above TH % are extracted (steps S**4202**, S**4203**, S**4204**). Then the barycenter of the upper densities above TH % is regarded as an approximate barycenter of the rectangular area.$\begin{array}{cc}\stackrel{\_}{x}=\frac{\int \int \pounds \left(x,y\right)xdxdy}{\int \int \pounds \left(x,y\right)dxdy}& \left(68\right)\\ \stackrel{\_}{y}=\frac{\int \int \pounds \left(x,y\right)ydxdy}{\int \int \pounds \left(x,y\right)dxdy}& \left(69\right)\end{array}$ - [0194]Here, f(x,y) represents an image of the upper densities above TH % and {overscore (y)}, {overscore (x)} stand for the barycenter.
- [0195]As described above, since the present embodiment is arranged to calculate the barycenter of data above the fixed density value, the present embodiment can extract the approximate barycenter of the area of interest with accuracy without being affected by the peripheral data, even if the object area is in the peripheral part of the two-dimensional image data.
- [0196]Further, since the start point c (
FIG. 17 ) of the aforementioned rotation axis f is placed in the object area, degrees of change in accumulated quantities of area ends become large against rotation of the rotation axis and the present embodiment can thus enjoy the effect of capability of extracting the rotation angle and area with better accuracy. - [0197]The present embodiment is arranged as follows. When the end point extraction means
**4102**(FIG. 15 ) extracts an end point, a point is not extracted as a candidate for an end point if f maxi<Th at the point, f maxi being defined below.

*f*max*i*=max{*fi*(*x*)|*xεk}*(70)

Here, “k” indicates a domain of definition of row data and Th a fixed threshold value. - [0198]As described above, the present embodiment can eliminate end points that can be points outside the object area with high possibilities and thus can extract the area and the angle of the area with better accuracy.
- [0199]The present embodiment uses a first order differential or a higher order differential as a characteristic quantity at the characteristic quantity calculation means
**4101**(FIG. 15 ). When Di(x) represents differential values of the i-th row, the end point extraction means**4102**(FIG. 15 ) extracts a point satisfying the condition of xi≧th as a candidate for an end point on the occasion of extracting a candidate for the left end point of the i-th row. - [0200]As described above, the present embodiment has the effect of shorter computation time, because it adopts the simple differentiation.
- [0201]
FIG. 19 is a block diagram to show the structure of an area extraction device according to the present embodiment. InFIG. 19 , reference numeral**4301**designates a secondary difference value calculation means for calculating the secondary difference values of one-dimensional image data in a designated direction (for example, which is determined according to the value determined by the angle extraction device described above),**4302**a left end point extraction means for extracting a left end point of the area, based on the secondary difference values calculated at the secondary difference value calculation means**4301**, and**4303**a right end point extraction means for extracting a right end point of the area, based on the secondary difference values calculated at the secondary difference value calculation means**4301**.FIG. 20 is a flow chart of a processing procedure sequence in the area extraction device according to the present embodiment. - [0202]The flow of the processing in the present embodiment will be described according to
FIG. 20 . The secondary difference value calculation means**4301**(FIG. 19 ) calculates the secondary difference values SS(x) according to a calculation equation defined by Eq. (71) below (step S**4301**). Here, f(x) represents the one-dimensional data of a line crossing the area in the designated direction and x represents coordinates thereof. Further, “d” denotes a constant indicating a difference distance.

*SS*(*d*)=*f*(*x−d*)−2×*f*(*x*)+*f*(*x+d*) (71) - [0203]The left end point extraction means
**4302**(FIG. 19 ) extracts the left end point x1 according to Eq. (72) below (step S**4302**). Here, {overscore (x)} represents a coordinate on the horizontal axis in the area.

*SS*(*x***1**)=min{*SS*(*x*)|0*≦x≦{overscore (x)}}*(72) - [0204]Then the right end point extraction means
**4203**(FIG. 19 ) extracts the right end point x**2**according to Eq. (73) below (step S**4303**).

*SS*(*x***2**)=min{*SS*(*x*)|{overscore (x)}*≦x*≦Length} (73)

Here, “Length” indicates the length of the image data along the horizontal axis. - [0205]As described above, since the present embodiment uses the secondary difference values, it has the effect of capability of extracting the area with accuracy even if the density change is gentle at the area end.
- [0206]The present embodiment uses the projection of the image area as defined by Eq. (74) below, as the one-dimensional data f(x) in aforementioned Eq. (71).

*f*(*x*)=∫_{b}^{c}*f*(*x,y*)*dy*(74) - [0207]Here, f(x,y) represents the image data in the designated direction and x, y-coordinates on the horizontal and vertical axes, respectively (which are inclined in a certain direction). Further, “b”, “c” represent an arbitrary domain across the object area.
- [0208]As described above, since the present embodiment uses the projection in the image area, it has the same effect as in the case of averaging the data of the image area. Therefore, the present embodiment can prevent the extraction from being affected by the noise inside and outside the object area and thus can extract the area with better accuracy.
- [0209]
FIG. 21 is a block diagram to show the structure of an area extraction device according to the present embodiment. InFIG. 21 , reference numeral**4401**designates a secondary difference value calculation means for calculating the secondary difference values of one-dimensional image data in a designated direction (for example, which is determined according to the value determined by the angle extraction device described above) indicated by a coordinate indication means**4404**,**4402**a left end point extraction means for extracting a left end point of the area, based on the secondary difference values calculated at the secondary difference value calculation means**4401**, a right end point extraction means for extracting a right end point of the area, based on-the secondary difference values calculated at the secondary difference value calculation means**4401**, and**4404**the coordinate indication means for indicating a coordinate of one-dimensional data for calculation of the secondary difference values at the secondary difference value calculation means**4401**after extraction of the right area end at the right end point extraction means**4403**. - [0210]Numeral
**4405**denotes a memory means for storing coordinates of right end points and left end points extracted at the left end point extraction means**4402**and at the right end point extraction means**4403**, and**4406**an average coordinate calculation means for calculating an average coordinate for each end point from the left and right ends stored in the memory means**4405**. - [0211]
FIG. 22 is a flow chart of a processing procedure sequence of the area extraction device according to the present embodiment. The operation will be described according toFIG. 22 . - [0212]The secondary difference value calculation means
**4401**(FIG. 21 ) calculates the secondary difference values SSi(x) according to a calculation equation defined by Eq. (75) below (step S**4401**). Here, fi(x) represents one-dimensional data of a line crossing the area in the designated direction, x coordinates thereof, and “i” a coordinate of the row in the designated direction indicated by the coordinate indication means**4404**. Here, “d” represents a constant indicating a difference distance.

*SSi*(*x*)=*fi*(*x−d*)−2*×fi*(*x*)+*fi*(*x+d*) (75) - [0213]Then the left end point extraction means
**4402**(FIG. 21 ) extracts the left area end XLi according to Eq. (76) below (step S**4402**). - [0214]Here, {overscore (x)} represents a coordinate on the axis along the designated direction in the area.

*SSi*(*XLi*)=min{*SSi*(*x*)|0*≦x≦{overscore (x)}}*(76) - [0215]Then the memory means
**4405**(FIG. 21 ) stores the left area end XLi extracted at the left end point extraction means**4402**(step S**4403**). - [0216]Next, the right end point extraction means (
FIG. 21 ) extracts the right area end XRi according to Eq. (77) below (step S**4404**).

*SSi*(*XRi*)=min{*SSi*(*x*)|*{overscore (x)}≦x*≦Length} (77) - [0217]Here, “Length” represents the length of the image data along the horizontal axis in the predetermined direction.
- [0218]Then the memory means
**4405**(FIG. 21 ) stores the right area end XRi extracted at the right end point extraction means**4403**(step S**4404**). - [0219]Next, the coordinate indication means
**4404**(FIG. 21 ) indicates a coordinate of a row in a new designated direction and step S**4401**to step S**4405**are repeated. An indication of the end of the loop is issued from the coordinate indication means**4404**(FIG. 21 ) (step S**4406**). - [0220]Next, after completion of the above extraction of the left and right end points, the average coordinate calculation means
**4406**(FIG. 21 ) calculates an average of the coordinates of the left and right end points stored in the memory means**4405**(step S**4407**). - [0221]As described above, since the present embodiment uses the secondary difference values, the boundary points of the area can be extracted with accuracy even if the density values vary gently at the boundary of the area.
- [0222]Further, since the area end is determined from the average of the area ends on plural lines, the present embodiment has the effects of being more resistant to the noise and capable of extracting the area with higher accuracy than in the case using a single point.
- [0223]In the present embodiment, where the image data of a one-dimensional line in the designated direction is expressed by f(x) and the secondary difference values thereof by SS(x) defined by Eq. (78) below, the sign of the primary difference value S(XL) defined by Eq. (79) below is added to the extraction of the left end point XL at the left end point extraction means
**4302**(FIG. 19 ) or**4402**(FIG. 21 ).

*SS*(*x*)=*f*(*x−d*)−2×*f*(*x*)+*f*(*x+d*) (78)

*S*(*XL*)=*f*(*XL*)−*f*(*XL−d*) (79) - [0224]For example, for extracting a left end point, XL that makes s(XL) negative and that satisfies Eq. (80) below is regarded as a left end point.

*SS*(*XL*)=min{*SS*(*x*)|0*≦x≦{overscore (x)}}*(80)

*S*(*XR*)=*f*(*XR+d*)−*f*(*XR*) (81) - [0225]Likewise, where the right end point is extracted at the right end point extraction means
**4303**(FIG. 19 ) or**4403**(FIG. 21 ), XR that makes S(XR) defined by above Eq. (81) negative and that satisfies Eq. (82) below is regarded as a right area end.

*SS*(*XR*)=min{*SS*(*x*)|*{overscore (x)}≦x*≦Length} (82) - [0226]As described above, since the present embodiment is arranged to add the sign of the primary difference to the extraction of the area end, it can also take the inclination of the image data outside the area into consideration. Therefore, erroneous extraction can be prevented at a quickly changing portion of density in the area, so that the area can be extracted with better accuracy.
- [0227]In the present embodiment, where the one-dimensional data in the designated direction is defined by f(x) and values resulting from the filtering process according to Eqs. (83), (84) below are defined by F
**1**(*x*), F**2**(*x*), the values F**2**(*x*) are used for calculation of the secondary difference values at the secondary difference value calculation means**4301**(FIG. 19 ) or**4401**(FIG. 21 ).

*F***1**(*x*)=min{*f*(*x+x***1**)−*h*(*x***1**)|−*d≦x***1***≦d*) (83)

Further, F**2**(*x*) is defined as follows.

*F***2**(*x*)=max{*F***1**(*x−x***1**)+*h*(*x***1**)|−*d≦x***1***≦d*} (84) - [0228]Here, h(x) is a function defined as follows.
$\begin{array}{cc}\begin{array}{cc}h\left(x\right)=\text{\hspace{1em}}0,& -d\le \mathrm{x1}\le d\\ =\infty ,& \mathrm{otherwise}\end{array}\hspace{1em}& \left(85\right)\end{array}$ - [0229]As described above, since the present embodiment is arranged to smooth the one-dimensional image data for the computation of the secondary difference values by the filtering process, it has the effect of being not affected by the noise, particularly, by the noise on the line.
- [0230]While the present invention has been described with respect to what is presently considered to be the preferred embodiments, it is to be understood that the invention is not limited to the disclosed embodiments. To the contrary, the invention is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.

Patent Citations

Cited Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|

US4536797 * | Oct 12, 1983 | Aug 20, 1985 | Rca Corporation | Television receiver with auxiliary on-screen display |

US4710822 * | Aug 21, 1986 | Dec 1, 1987 | Konishiroku Photo Industry Co., Ltd. | Image processing method |

US4859850 * | Jan 12, 1988 | Aug 22, 1989 | Fuji Photo Film Co., Ltd. | Irradiation field recognizing method, and method of adjusting image processing conditions using the same |

US4970393 * | Nov 25, 1987 | Nov 13, 1990 | Fuji Photo Film Co., Ltd. | Irradiation field recognizing method, and method of adjusting image processing conditions using the same |

US4992663 * | Mar 17, 1989 | Feb 12, 1991 | Fuji Photo Film Co., Ltd. | Method of judging the correctness or incorrectness of a prospective contour point of an irradiation field |

US5046147 * | Oct 17, 1989 | Sep 3, 1991 | Fuji Photo Film Co., Ltd. | Radiation image processing method |

US5081580 * | Apr 20, 1989 | Jan 14, 1992 | Fuji Photo Film Co., Ltd. | Method for recognizing an irradiation field, and method for judging the correctness or incorrectness of prospective contour points of an irradiation field |

US5091970 * | Jul 25, 1990 | Feb 25, 1992 | Fuji Photo Film Co., Ltd. | Method for judging the presence or absence of limited irradiation field |

US5115476 * | Jan 24, 1991 | May 19, 1992 | Fuji Photo Film Co., Ltd. | Edge finding method and apparatus |

US5123054 * | Oct 19, 1990 | Jun 16, 1992 | Fuji Photo Film Co., Ltd. | Abnormal pattern detecting apparatus |

US5151947 * | Nov 26, 1990 | Sep 29, 1992 | Konica Corporation | Method and apparatus for processing radiation image |

US5495536 * | May 8, 1991 | Feb 27, 1996 | Sandia Corporation | Image processing system and method for recognizing and removing shadows from the image of a monitored scene |

US5533142 * | Oct 14, 1993 | Jul 2, 1996 | Fuji Photo Film Co., Ltd. | Method for adjusting read-out and processing conditions for Magen images |

US5557687 * | Feb 15, 1995 | Sep 17, 1996 | Fuji Photo Film Co., Ltd. | Abnormal pattern detecting apparatus |

US5659402 * | Jan 3, 1995 | Aug 19, 1997 | Mita Industrial Co., Ltd. | Image processing method and apparatus |

US5668889 * | Jun 7, 1995 | Sep 16, 1997 | Fuji Photo Film Co., Ltd. | Apparatus for determining an image position, and method for adjusting read-out conditions and/or image processing conditions for a radiation image |

US5732149 * | May 9, 1997 | Mar 24, 1998 | Konica Corporation | Apparatus for extracting an irradiation field region from a radiation image |

US5838758 * | Mar 13, 1995 | Nov 17, 1998 | Vivid Technologies | Device and method for inspection of baggage and other objects |

US5850465 * | Aug 28, 1991 | Dec 15, 1998 | Fuji Photo Film Co., Ltd. | Abnormnal pattern detecting or judging apparatus, circular pattern judging apparatus, and image finding apparatus |

US5850836 * | Sep 22, 1995 | Dec 22, 1998 | Hologic, Inc. | Morphometric x-ray absorptiometry (MXA) |

US5875225 * | Jun 3, 1997 | Feb 23, 1999 | Siemens Aktiengesellschaft | Spiral scan computed tomography apparatus, and method for image reconstruction therein |

US5892840 * | Feb 29, 1996 | Apr 6, 1999 | Eastman Kodak Company | Method and apparatus for irradiation field detection in digital radiographic images |

US5960102 * | Sep 16, 1996 | Sep 28, 1999 | U.S. Philips Corporation | X-ray image processing method and device for performing that method in which a portion corresponding to an x-ray absorption filter is selectively processed |

US6002959 * | Dec 9, 1997 | Dec 14, 1999 | Hologic, Inc. | Morphometric x-ray absorptiometry (MXA) |

US6057847 * | May 9, 1997 | May 2, 2000 | Jenkins; Barry | System and method of image generation and encoding using primitive reprojection |

US6061465 * | Aug 29, 1997 | May 9, 2000 | Fuji Photo Film Co., Ltd. | Radiation image processing method and apparatus |

US6111582 * | Jul 23, 1997 | Aug 29, 2000 | Jenkins; Barry L. | System and method of image generation and encoding using primitive reprojection |

US6198837 * | Aug 19, 1997 | Mar 6, 2001 | Konica Corporation | Method of storing medical image information recorded in sheet shaped recording medium |

US6243485 * | Mar 12, 1998 | Jun 5, 2001 | Fuji Photo Film Co., Ltd. | Method and apparatus for recognizing irradiation fields on radiation images |

US6322505 * | Jun 8, 1999 | Nov 27, 2001 | Acuson Corporation | Medical diagnostic ultrasound system and method for post processing |

US6352511 * | Nov 2, 2000 | Mar 5, 2002 | Acuson Corporation | Medical diagnostic ultrasound system and method for post processing |

US6389310 * | Mar 2, 2000 | May 14, 2002 | Georges Demonceau | Method and apparatus for analyzing heart function using 4D ECG synchronized heart cavity tomoscintigraphy |

US6415049 * | Apr 14, 1999 | Jul 2, 2002 | Konica Corporation | Apparatus for detecting and processing a radiation image |

US6445771 * | Apr 27, 2000 | Sep 3, 2002 | Koninklijke Philips Electronics N.V. | X-ray device incorporating detection of shutter edges |

US6600831 * | Mar 21, 2000 | Jul 29, 2003 | Fuji Photo Film Co., Ltd. | Connection processing method for radiation images |

US6611615 * | Jun 25, 1999 | Aug 26, 2003 | University Of Iowa Research Foundation | Method and apparatus for generating consistent image registration |

US6671625 * | Feb 17, 2000 | Dec 30, 2003 | Vialogy Corp. | Method and system for signal detection in arrayed instrumentation based on quantum resonance interferometry |

US6740883 * | Aug 13, 1999 | May 25, 2004 | Robert Z. Stodilka | Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas |

US6768782 * | Dec 16, 2002 | Jul 27, 2004 | University Of Notre Dame Du Lac | Iterative method for region-of-interest reconstruction |

US6832001 * | Nov 15, 2000 | Dec 14, 2004 | Canon Kabushiki Kaisha | Image processing apparatus and method |

US6850596 * | Jun 12, 2001 | Feb 1, 2005 | Ge Medical Systems | Collimation device, radiology apparatus, test kit and method of testing a radiology apparatus |

US6901158 * | Sep 17, 1998 | May 31, 2005 | Canon Kabushiki Kaisha | Image discrimination apparatus and image discrimination method |

US20020146093 * | Apr 10, 2001 | Oct 10, 2002 | Williams John J. | X-ray system alignment method and apparatus |

US20030018457 * | Mar 11, 2002 | Jan 23, 2003 | Lett Gregory Scott | Biological modeling utilizing image data |

US20030095698 * | Nov 7, 2002 | May 22, 2003 | Konica Corporation | Feature extracting method, subject recognizing method and image processing apparatus |

US20040125103 * | Feb 26, 2001 | Jul 1, 2004 | Kaufman Arie E. | Apparatus and method for volume processing and rendering |

US20040179646 * | Mar 14, 2003 | Sep 16, 2004 | Jianying Li | Imaging systems and methods |

US20040197727 * | Jul 14, 2003 | Oct 7, 2004 | Orametrix, Inc. | Method and system for comprehensive evaluation of orthodontic treatment using unified workstation |

US20040223591 * | Jun 14, 2004 | Nov 11, 2004 | Dietmar Sundermann | Method of testing a radiology apparatus and program |

US20050004617 * | Apr 28, 2004 | Jan 6, 2005 | Dawant Benoit M. | Apparatus and methods of optimal placement of deep brain stimulator |

US20050069187 * | Sep 27, 2004 | Mar 31, 2005 | Konica Minolta Medical & Graphic, Inc. | Image processing method, image processing apparatus and image processing program |

Referenced by

Citing Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|

US8081811 * | Apr 11, 2008 | Dec 20, 2011 | Fujifilm Corporation | Method, apparatus, and program for judging image recognition results, and computer readable medium having the program stored therein |

US20080260226 * | Apr 11, 2008 | Oct 23, 2008 | Fujifilm Corporation | Method, apparatus, and program for judging image recognition results, and computer readable medium having the program stored therein |

Classifications

U.S. Classification | 382/132 |

International Classification | G06K9/32, G06F19/00, G06K9/00 |

Cooperative Classification | G06T7/0083, G06F19/321, G06K9/3233, G06T7/0008, G06T2207/10116, G06T5/50, G06T2207/30004 |

European Classification | G06K9/32R, G06T5/50, G06T7/00S2, G06T7/00B1P |

Rotate