Method for determining median line

Disclosed is a method for automatically determining the median line of a selected region of a human body so as to simplify complex imaging required hithero for a magnetic resonance imaging apparatus, and, in the method, an artifact of the imaged selected region of the human body is utilized to automatically determine the median line. For example, two artifacts 1010 and 1020 intersecting each other appear from the superior sagittal sinus as shown in FIG. 1 when the direction of phase encoding is changed. By subtracting one of these artifacts 1010 and 1020 from the other and calculating the absolute value of the result of subtraction, the point 1030 of intersection of the artifacts determines the location of the superior sagittal sinus which is the source of apperance of the artifacts, so that the median line can be automatically determined on the basis of the result of determination of location of the superior sagittal sinus.

BACKGROUND OF THE INVENTION 
This invention relates to a magnetic resonance imaging (abbreviated 
hereinafter as an MRI) apparatus utilizing the phenomenon of magnetic 
resonance, and more particularly to a method in which a human body is 
imaged by the MRI apparatus to determine a selected region of the human 
body on the image so as to determine the median line of the selected 
region of the human body. 
In an MRI apparatus, a magnetic field is applied to a human body to induce 
the phenomenon of magnetic resonance of protons in a selected region of 
the human body, and the resultant magnetic resonance signal is detected 
and measured to display a tomographic image of the selected region of the 
human body. When the mode of application of the magnetic field (referred 
to hereinafter as a sequence) is changed in this MRI apparatus, the 
corresponding resonance signal can be obtained, so that, for example, the 
blood vessel only can be selectively imaged. 
When a desired tomographic image of a selected region of a human body is to 
be obtained by this MRI apparatus, a plurality of tomographic images are 
first taken, and, after determining the desired slice direction on the 
basis of the anatomical information derived from these tomographic images, 
the most suitable sequence is selected so as to take the desired 
tomographic image. Thus, it is required to take a plurality of tomographic 
images for the purpose of positioning and also to determine the slice 
direction in order to obtain the desired tomographic image. 
In our earlier invention described in patent application Ser. No. 07/870295 
filed in the United States on Apr. 17, 1992, a method for automatically 
determining a sequence used in a usual diagnosis (a routine diagnosis) and 
a slice direction is registered beforehand, so that the desired 
tomographic image can be obtained by a single manipulation. In regard to 
the positioning, the patent application describes that characteristic 
regions of a human body are detected for the purpose of positioning and 
describes also a method for positioning, for example, the OM line of the 
head or the sagittal section of a transverse image of the neck. 
Besides the method described above, a method in which the median line of a 
human body commonly used as an index in, for example, image measurement by 
the Hough transform is reported in IEEE TRANSACTIONS ON MEDICAL IMAGING, 
VOL. 10, NO. 1, MARCH 1991. 
The above report refers to a method for detecting the median line in a 
transverse image and a coronal image of the head of a human body. 
According to the reported method, the head image is subjected to filtering 
(the Sobel operation) so as to be converted into an edge emphasized head 
image, and its binary image is subjected to the Hough transform so as to 
provide the desired median line. 
When a different sequence is used in the MRI apparatus imaging a human 
body, a correspondingly different resonance signal is generated from the 
human body. This feature is widely utilized in the research and 
development of various sequences. That is, not only detection of a proton 
density image but also implementation of an angio sequence for selective 
measurement of a blood vessel and an ultra-high speed sequence for imaging 
with a very short imaging period of time are widely attempted. 
Under such a background, MRI apparatuses have been improved from year to 
year to be able to operate with a higher performance, and the applicable 
range has also been progressively widened. However, with the improvement 
in the performance of the MRI apparatuses, they have inevitably the 
problem of requiring a complicated manipulation. 
When a tomographic image of a desired region of a human body is required by 
a doctor for the purpose of diagnosis, a plurality of tomographic images 
are first taken, and, on the basis of the anatomical information derived 
from those tomographic images, the doctor determines the desired slice 
direction so as to obtain the desired tomographic image. Among the 
information used by the doctor during the process of imaging, the median 
line [which is the intersection line between the surface of the desired 
region of the human body and a perpendicular plane (a superior sagittal 
section) including a sagittal direction of a horizontal line penetrating 
the human body, that is, a plane (a median section) dividing the human 
body into its right and left halves] is very important information for 
accurately obtaining the tomographic image of the desired region of the 
human body. 
The Hough transform described in the above paper as a method of determining 
the median line uses the sine and cosine functions for the determination. 
Therefore, the disclosed method used for the determination of the median 
line has the problem that an increase in the number of factors to be 
processed results in a corresponding increase in the period of time 
required for the processing. 
SUMMARY OF THE INVENTION 
With a view to solve the prior art problems described above, it is an 
object of the present invention to provide a method of determining the 
median line by a simple procedure. 
Imaging by an MRI apparatus is broadly classified into two cases. In one 
case, a desired region of an object to be inspected, for example, a human 
body to be diagnosed is supposed to be stationary, and, in the other case, 
the desired region of the objected to be inspected, for example, the human 
body to be diagnosed is supposed to be moving. 
In the case of imaging the desired region of the human body which is 
supposed to be stationary, the image will not be appreciably adversely 
affected by the flow of blood flowing through a small-diameter blood 
vessel, such as, a capillary vessel. However, when the desired region 
includes a large-diameter blood vessel, an artifact tends to appear on the 
image. 
This artifact is affected by the factors including the blood flow velocity, 
the slice thickness and the imaging period of time and appears in the 
phase encoding direction. In the present invention, the artifact generated 
from such a moving region of the human body is utilized so as to determine 
the median line according to the following three steps (1), (2) and (3): 
(1) Utilizing the tendency of appearance of the artifact in the phase 
encoding direction, two images are taken according to two sequences having 
respectively different directions of phase encoding, and the point of 
intersection between these artifacts is determined as the location of the 
large-diameter blood vessel. 
(2) Without the T2 weighting, the interhemispheric fissure is blackend on 
the image, and the image density of that portion of the image becomes 
lower than that of the remaining portion of the brain. (A T1 weighted 
image, a proton density image, etc. are also blackened.) Utilizing this 
tendency, the location of the interhemispheric fissure is determined by 
tracing on the basis of the image density. 
(3) Utilizing the fact that the median line exists near the center of 
gravity of the head of the human body, the median line is determined. 
The method of the present invention will now be more concretely described. 
(1) When a transverse image of the head is taken according to a sequence 
which does not restrict the blood flow, a linear artifact extending in the 
phase encoding direction from the superior sagittal sinus having an 
especially large diameter in the brain appears, as shown in FIGS. 1A-1C. 
This superior sagittal sinus lies behind the median line and is located at 
a position near the outer skin layer of the skull. Therefore, by 
determining the location of this superior sagittal sinus, the information 
regarding the position of one end of the median line can be acquired. 
(2) In the brain, the interhemispheric fissure is the region corresponding 
to the median line (plane), as shown in FIG. 8. As described already, 
without the T2 weighting, the resonance signal from the interhemispheric 
fissure is not emphasized, resulting in blackening of that portion on the 
image. Also, the superior sagittal sinus is located near the terminating 
end of the interhemispheric fissure. Utilizing the above fact, the 
interhemispheric fissure is traced from the location of the superior 
sagittal sinus selected as a trace starting point, and an approximate 
straight line of the trace line is denoted as the median line. 
(3) The median line is located on the plane dividing the human body into 
its right and left halves. Therefore, in the step (2) for determining the 
median line of the head of the human body, the trace starting point is set 
at the center of gravity of the head region. 
In the present invention, the object to be inspected is not limited to the 
human body.

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
FIG. 2 is a block diagram showing the structure of an MRI apparatus to 
which the present invention is applied. Referring to FIG. 2, the MRI 
apparatus includes a static magnetic field generation system 2010 for 
generating a homogeneous static magnetic field, a radio-frequency magnetic 
field generation system 2020 for radiating a radio-frequency magnetic 
field, and a gradient magnetic field generation system 2030 capable of 
changing the magnetic field strength in the directions of the X, Y and Z 
axes independently of one another. The magnetic field generation systems 
2010, 2020 and 2030 are controlled to induce the phenomenon of nuclear 
magnetic resonance in a human body to be diagnosed, and a signal 
measurement system 2040 receives an electromagnetic field generated from 
the human body being diagnosed. A processor 2050 implements image 
reconstruction on the basis of the received data, and the result of image 
reconstruction is displayed on a CRT 2060. A sequence controller 2070 
controls the magnetic field generation systems 2010, 2020, 2030, the 
signal measurement system 2040 and the processor 2050. The sequence 
controller 2070 includes therein a memory for storing a sequence program 
and sequence parameters. The processor 2050 carries out processing for 
determining a desired region of the human body according to the method of 
the present invention. 
The MRI apparatus carries out positioning of the desired region according 
to the chain oblique method (a method in which, when a line is drawn on an 
image, a slice position is set in a direction aligning with that line and 
perpendicular with respect to the image plane.) FIGS. 3A-3D show an 
imaging procedure for taking a transverse image of the head of the human 
body by the MRI apparatus. 
The patient is laid almost straight on a bed, but his axial section does 
not always accurately register with that of a desired tomographic image. 
When a highly accurate (distortion-free) tomographic image is required, 
imaging is desirably done according to the procedure shown in FIG. 3. A 
first slice direction 3010 is set to align with the direction of 
transverse imaging as shown in FIG. 3A, and a transverse image of the head 
of the patient is taken in that direction. Then, a second slice direction 
3020 is set to be perpendicular with respect to the direction of the 
median line 3021 of the transverse image as shown in FIG. 3B, and a 
coronal image is taken. Then, a third slice direction 3030 is set to align 
with the median line of the coronal image as shown in FIG. 3C, and imaging 
is carried out. Finally, a fourth slice direction is set to align with the 
OM line 3040 as shown in FIG. 3D, and the desired tomographic image is 
taken. 
In order to continuously implement the series of imaging described above, 
automatic positioning of the median line is indispensably required. 
How to determine the median line by approximation by a plurality of 
procedures will now be described by taking the case of automatic 
positioning of the median line of the transverse image of the head as an 
example. These procedures are as follows: 
(1) A procedure in which the location of the superior sagittal sinus lying 
behind the median line is determined by utilization of an artifact so as 
to determine the median line. 
(2) A procedure in which the interhemispheric fissure is traced from the 
location of the superior sagittal sinus so as to determine the median line 
by approximation. 
(3) A procedure in which the center of gravity of the tomographic image is 
based to trace the interhemispheric fissure so as to determine the median 
line by approximation. 
The automatic positioning of the median line on the transverse image of the 
head is as described above. A coronal image of the head is also based to 
determine the median line according to the following procedure: 
(4) A procedure in which the coronal image of the head is based to trace 
the interhemispheric fissure from the location of the center of gravity of 
the tomographic image so as to determine the median line by approximation. 
These procedures will be described in more detail. 
(1) Procedure in which an artifact is utilized so as to determine the 
location of the superior sagittal sinus: 
Because of the blood flow through the superior sagittal sinus, an artifact 
appears when the image is taken by a sequence which does not restrict the 
blood flow. Because this artifact appears in the direction of phase 
encoding, artifacts intersecting each other can be produced by changing 
the direction of application of phase encoding. In the illustrated 
embodiment of the present invention, the point of intersection between 
these intersecting artifacts is determined as the locaton of the superior 
sagittal sinus. 
FIGS. 4A and 4B show examples of the imaging sequence. This sequence is 
that commonly called the spin echo method. FIGS. 4A and 4B differ from 
each other in that the direction of application of a gradient magnetic 
field 4010 for phase encoding purpose and the direction of application of 
a gradient magnetic field 4020 for signal reading purpose are exchanged. 
FIGS. 1A and 1B show the images taken according to these two sequences 
respectively. Artifacts 1010 and 1020 appear on these two images 
respectively. Because these two images are taken on the same region, one 
of the artifacts is subtracted from the other to obtain the absolute value 
as shown in FIG. 1C. In this case, an undesirable position deviation of 
the images may result due to, for example, distortions of the static and 
gradient magnetic fields. In such a case, the subtraction is preferably 
carried out after execution of various kinds of position deviation 
correction processing already known in the art. The values of pixels of 
the image shown in FIG. 1C are integrated in both the longitudinal and 
lateral directions, and the points where the integrated values are largest 
in both the axial directions are selected as the positions of the 
artifacts, so that the point of intersection 1030 can be determined. This 
point of intersection 1030 is determined as the location of the superior 
sagittal sinus. 
(2) Procedure for determining the median line by tracing the 
interhemispheric fissure from the location of the superior sagittal sinus: 
According to this procedure, the location of the superior sagittal sinus 
determined by the procedure (1) is based to trace the interhemispheric 
fissure, and a straight line obtained by approximating the trace line by 
the method of least squares is determined as the median line. FIGS. 5, 6 
and 7 are flow charts of this procedure, and FIG. 8 shows the median line 
determined by the procedure. 
(Process 1) Separation of the head region from the background region: 
First, a threshold for separating the head of the patient from the 
background on the image is detected. Generally, the histogram of the image 
obtained by imaging by the MRI apparatus is in the form having two ridge 
portions. In the histogram, the ridge portion having a lower frequency 
distribution or density represents the background of the image, while the 
ridge portion having a higher density represents the region of the human 
body. Also, the maximum value of the ridge portion representing the 
background of the image provides the maximum value of the histogram. In 
this process, the value of the portion (the valley) intermediate between 
the two ridge portions in the histogram is used as the threshold. 
FIG. 5 is a flow chart of a sequence of steps for calculating the threshold 
from the histogram. The process includes a step 5010 of setting various 
initial values, a step 5020 of reading image data, a step 5030 of 
calculating the maximum value of the image data, and steps 5040, 5050, 
5060, 5070 and 5080 of preparing the histogram. The process further 
includes a step 5090 of caculating the maximum value of the histogram 
prepared by the above steps, steps 5100, 5110 and 5120 of detecting the 
valley between two ridge portions of the histogram, and a step 5130 of 
calculating the threshold on the basis of the detected valley. When the 
valley does not exist in the histogram, the threshold is set in a step 
5140. 
Then, the edge line of the head relative to the background is determined. 
FIG. 6 is a flow chart of a sequence of steps for determining the edge 
line of the head on the basis of the threshold. The process includes a 
step 6010 of setting various initial values, steps 6020, 6030, 6040, 6050, 
6060 and 6070 in which points having values larger than the threshold are 
searched in the direction of from the left end toward the right end of the 
image so as to determine the edge points, and steps 6080, 6090, 6100, 
6110, 6120 and 6130 in which points having values larger than the 
threshold are searched in the direction of from the right end toward the 
left end of the image so as to determine the edge points. The zone inside 
these edge points is determined as the head region, while that outside the 
edge points is determined as the background region. 
(Process 2) Trace of the interhemispheric fissure from the location of the 
superior sagittal sinus: 
FIG. 7 is a flow chart of a sequence of steps for tracing the 
interhemispheric fissure from the location of the superior sagittal sinus. 
The process includes a step 7010 of reading the image and a step 7020 of 
determining the threshold together with the edge points according to the 
method described already. In a step 7030, the range for tracing the 
interhemispheric fissure is determined. In this range determination, the 
position where the edge point is first detected on the Y line extending 
downward from the top of the image is selected as a final trace point 8010 
as shown in FIG. 8. The trace starts from the location of the superior 
sagittal sinus 8020 shown in FIG. 8. The tracing operation utilizes the 
fact that the interhemispheric fissure is blackened on the T1 weighted 
image. Thus, in steps 7040, 7050, 7060, 7070, 7080, 7090, 7100 and 7110, 
pixels having the minimum image density within the extent of the width of 
a plurality of pixels in the X direction are detected on a line connecting 
trace points. This manner of tracing is repeated until the 
interhemispheric fissure is traced up to the final trace point 8010 
thereby deciding a trace line 8030 as shown in FIG. 8. 
(Process 3) Determination of the median line on the basis of the trace 
line: 
On the basis of the decided trace line, an approximate straight line 8040 
as shown in FIG. 8 is calculated by approximation according to the method 
of least squares, and this line 8040 is denoted as the median line. The 
formula for calculating this straight line by approximation according to 
the method of least squares will now be described. 
Calculation of the straight line by approximation according to the method 
of least squares: 
EQU Approximate straight line: X=aY+b 
EQU a=(.intg.A(y)Ydy.intg.dy=.intg.A(y)dy.intg.Ydy)/(.intg.(Y.multidot.Y)dy.int 
g.dy-.intg.Ydy.intg.Ydy) 
EQU b=(.intg.A(y)dy.intg.(Y.multidot.Y)dy-.intg.(A(y)Y)dy.intg.Ydy)/(.intg.(Y.m 
ultidot.Y)dy.intg.dy-.intg.Ydy.intg.Ydy) 
In the above equations, A(y): the X coordinate value corresponding to the Y 
coordinate value on the line tracing the interhemispheric fissure. 
(3) Procedure for determining the median line by tracing the 
interhemispheric fissure on the basis of the location of the center of 
gravity of the tomographic image: 
In this procedure, a point located near the center of gravity of the 
tomographic image is selected as the trace starting point. 
First, the head region is separated from the background region on the image 
according to the procedure (2) described above. Then, the center of 
gravity of the head region is calculated according to the following 
formula: 
Center of gravity: 
EQU X=.intg..intg.x*f(x,y)dxdy/.intg..intg.f(x,y)dxdy 
EQU Y=.intg..intg.y*f(x,y)dxdy/.intg..intg.f(x,y)dxdy 
In the above equations: 
x: the x coordinate value on the image, 
y: the y coordinate value on the image, and 
f(x, y): image data. The head region=1, and the background region=0. 
The interhemispheric fissure exists near the center of gravity of the head 
region. Therefore, the point where the pixels having the minimum density 
are located in a rectangular range within the extent of the width of a 
plurality of pixels located at a position near the center of gravity of 
the head region is detected, and the trace is started from that point. The 
manner of tracing and the manner of calculating the straight line by 
approximation according to the method of least squares are similar to 
those described already in the procedure (2). In this case, however, the 
trace proceeds downward from the position of the center of gravity of the 
head region. 
(4) Procedure for determining the median line on a coronal image of the 
head: 
FIG. 9 shows a coronal image of the head. In the coronal image of the head, 
the interhemispheric fissure does not traverse the tomographic image and 
exists in the upper part of the tomographic image as shown. 
In order to determine the median line, the head region is first separated 
from the background region in a manner similar to the procedure described 
already, and the edge line and the point of the center of gravity 9010 of 
the head are determined. Then, the top point 9020 of the edge line is 
determined. The determination of the top point of the edge line is such 
that the Y coordinate value first detected during a downward trace from 
the top of the image is selected as the edge top point 9020. A point 9060 
for determining a trace starting point 9030 has a Y coordinate value 
between the center of gravity 9010 and the edge top point 9020, and the X 
coordinate is the coordinate of the center of gravity. The pixels having 
the minimum density are detected from a rectangular range near the point 
9060 so determined so as to determine the trace starting point 9030. The 
interhemispheric fissure is traced from the trace starting point 9030 on 
the basis of the detected image density accordng to a procedure similar to 
that described already. A trace termination line 9040 is selected to be 
located beneath the edge top point 9020 but above the trace starting point 
9030. This is because the blood vessel runs in the vicinity of the edge 
top point 9020, and the blood flowing through this blood vessel may have 
such a high concentration which will obstruct the otherwise successful 
tracing depending on the method of imaging. The trace termination line 
9040 is selected to be located at such a position which avoids an 
undesirable trouble. A trace line 9050 is finally obtained as a result of 
tracing from the trace starting point 9030 to the trace termination line 
9040, and a straight line obtained by approximating the trace line 9050 
according to the method of least squares is selected as the median line. 
It will be understood from the foregoing description of the present 
invention that the median line and the median plane can be automatically 
detected, so that, in the case of imaging to take a desired tomographic 
image of a human body by means of positioning on the basis of the median 
line (plane), the positioning can be continuously attained. 
Further, because of the fact that the median line (plane) that is the index 
for implementing the desired accurate positioning can be automatically 
detected, the accuracy of determination of the imaging position can be 
improved.