Biomagnetic field measuring method and apparatus

Description Submit all comments and votes

BACKGROUND OF THE INVENTION

The present invention relates to biomagnetic field measuring method and apparatus for measuring a biomagnetic field generated by a nerve action of brain as well as a myocardial action of heart of a living body by using a plurality of fluxmeters each consisting of a highly sensitive superconducting quantum interference device (SQUID).

In addition to a magnetic field generated by a current dipole, a magnetic field due to a volume current flowing in the living body is enumerated as a biomagnetic field. Measurement of a normal component (Bz:z component in the Cartesian coordinate system or B.sub.r:radius component in the polar coordinate system) is considered to be hardly affected by the volume current. In conventional techniques, the plane of a detection coil connected to a SQUID is disposed in parallel to the body surface to measure B.sub.z or B.sub.r which is a normal component vertical to the body surface. Results of the biomagnetic field measurement are displayed in the form of a temporal change waveform of the measured field component or an isomagnetic field map (contour map) for connecting points at which magnitudes of the magnetic field component measured at desired time points are equal to each other. Various analysis methods have been proposed which analyze a magnetic field source participating in generation of the biomagnetic field from the obtained isomagnetic field curve and in typical one of them, analysis is carried out by replacing the magnetic field source with a current dipole.

An isomagnetic field map of a normal component (B.sub.z or B.sub.r) of the magnetic field generated by a current dipole is of a pattern having a source pole of magnetic field and a sink pole of magnetic field at positions which are separate from each other from the center where a magnetic field source (current dipole) is positioned. The magnitude, position and direction of the magnetic field source (current dipole) are analyzed in accordance with magnitudes of magnetic field at the two poles and a distance therebetween.

In a first prior art (H. Hosaka and D. Cohen: J. Electrocardiology, 9 (4), pp. 426-432 (1976)), a method is employed for displaying current sources distributed in the myocardium by using an isomagnetic field map of measured normal component B.sub.z with the aim of promoting visibleness of direction and intensity of currents in the myocardium and according to this method, an arrow map is contrived for expressing a current vector J (x, y) defined by equation (1) on measuring points by using an arrow. In the following description, Gothic characters are used to indicate vectors. J(x, y)=(.differential.B.sub.z(x, y)/.differential.y)e.sub.x-(.differential.B.sub.z(x, y)/.differential.x)e.sub.y (1)

In equation (1), e.sub.x designates a unit vector in x direction and e.sub.y designates a unit vector in y direction. This prior art, however, encounters a problem that when a plurality of current sources exist, it is difficult to discriminate the individual current sources from each other on the basis of the isomagnetic field map of normal component B.sub.z.

In a second prior art (K. Tukada et al: Review of the Scientific Instruments, 66(10), pp. 5085-5091 (1995)), for the sake of visualizing a plurality of distributed current sources, the normal component (Bz or B.sub.r) is not detected but tangential components B.sub.x and B.sub.y are measured by using a detection coil whose plane is disposed vertically to the body surface. Each of the measured tangential components B.sub.x and B.sub.y is displayed in the form of an isomagnetic field map. The tangential components B.sub.x and B.sub.y measured according to the second prior art are considered to be affected by the volume current but in an isomagnetic field map of two-dimensional vector magnitude B.sub.xy obtained by synthesizing B.sub.x and B.sub.y measured at time point t pursuant equation (2), a peak can always be obtained directly above a current dipole and therefore, even when a plurality of current dipoles exist, individual current dipoles can be separated for visualization. |B.sub.xy(x, y, t)= {square root over ( )}{(B.sub.x(x, y, t)).sup.2+(B.sub.y(x, y, t)).sup.2} (2)

In a third prior art (Y. Yoshida et al: Tenth International Conference on Biomagnetism, Santana Fe, N.Mex., Feb. 17 (1996)), a normal component and two tangential components of a biomagnetic field are detected by using a vector magnetic field sensor consisting of three detection coils having coil planes which are orthogonal to each other, detection results of the magnetic field components are converted in terms of the Cartesian coordinate system to determine Cartesian coordinate system components B.sub.x, B.sub.y and B.sub.z, and an isomagnetic field map of the normal component B.sub.z and an isomagnetic field map of two-dimensional vector magnitude B.sub.xy are displayed, respectively.

In a fourth prior art (K. Tsukada et al: Tenth International Conference on Biomagnetism, Santana Fe, N.Mex., Feb. 17 (1996)), two tangential components B.sub.x and B.sub.y of a biomagnetic field are detected and an isomagnetic field map based on |B.sub.XY|=|B.sub.x+B.sub.y| is compared with an isomagnetic field map based on a normal component B.sub.z.

Available as diagrams for indicating measurement results of electrical physiological phenomena in a living body are a magnetoencephalogram (MEG) obtained through measurement using a magnetoencephalogram and an electrocardiogram (ECG) obtained through measurement using an electrocardiograph. In measurements of the electrocardiogram, a body surface potential map for mapping an electrocardiographic figure by using a plurality of electrodes is of a well-known technique. The MEG or the body surface potential map is depicted in the form of an isopotential map for connecting isopotential points.

In a fifth prior art (T. J. Montague et al: Circulation 63, No. 5, pp. 1166-1172 (1981)), an isointegral map obtained by integrating a temporal change waveform of an output of each one of a plurality of electrodes over a desired time interval is depicted as a body surface potential map.

In the following description, "biomagnetic field" means "magnetic field generated from a living body", "cardiac magnetic field measurement" means "measurement of a magnetic field generated from the heart", and "cardiac magnetic wavaform" means "waveform indicated by a magnetocardiogram (MCG) obtained through cardiac magnetic field measurement". Further, "encephalic magnetic field measurement" means "measurement of a magnetic field generated from the brain" and "encephalic magnetic waveform" means "waveform indicated by a magnetoencephalogram (MEG) obtained through encephalic magnetic field measurement".

Each of the conventional isomagnetic field maps of the respective components has inherent features. In the presence of a single current dipole, the position, magnitude and direction of a current source can be analyzed with ease by using the isomagnetic field map of normal component B.sub.z. On the other hand, the isomagnetic field map of two-dimensional vector magnitude B.sub.xy obtained from measurement results of tangential components B.sub.x and B.sub.y features that even in the presence of a plurality of current dipoles, individual current dipoles can easily be discriminated from each other. But, for detection of a magnetic field, coils are required to be provided in x and y directions and the number of coils is doubled as compared to detection of only the normal component B.sub.z. In vector measurement for measuring all the components B.sub.x, B.sub.y and B.sub.z, the number of required coils is tripled as compared to detection of only the normal component B.sub.z. Accordingly, the magnetic field sensor consisting of a detection coil and a SQUID is increased in number and in addition, the signal processing circuit and the like are also increased in number, raising a problem that the biomagnetic field measuring system becomes an expensive one. Further, the first prior art is disadvantageous in that arrows are merely indicated on measuring points and detailed distribution states of current sources are hardly discriminated.

From the isomagnetic field map indicated in terms of a biomagnetic field component, the position, magnitude and direction, at a desired time point, of a current source in a living body can be analyzed and detailed information about changes in position, magnitude and direction of the current source can be known. Conventionally, dynamic changes in various kinds of information pieces are captured by using many figures displayed on or delivered to the apparatus so as to diagnose a disease. In the prior arts, however, many diagrams or maps indicating various kinds of information pieces are needed for diagnosis and abnormality of changes in various kinds of information pieces is known empirically. As will be seen from the above, in the prior arts, the processing of displaying, on a single map, systematic information as to what magnitude of current flows through which portion of a living body and as to which region an abnormal bio-current passes through is not executed. In the case of the body surface potential map, an isointegral technique was reported. This isointegral map was drawn by connecting between the same integral values over a desired time interval (for example, an interval during which waves of Q, R and S are generated and an interval during which S to T waves are generated). The advantage of this isointegral map is that information of the heart can be obtained from only a single electrocardiographic figure. But, in the isopotential map when the current source in the heart is assumed to be a single current dipole, a figure results disadvantageously in which an positive peak and a negative peak do not exist immediately above the current dipole but exist at a position which is separate from a point immediately above the current dipole. Further, when the position of the current dipole remains unchanged but the direction of the current dipole changes, the anode and cathode peak positions change, raising a problem that when potential is integrated, correspondence between the current source and the peak of an integral value is impaired. Like the case of the electrocardiogram, mere integration of a component of a biomagnetic field obtained through biomagnetic field measurement faces a problem that the peak position of the biomagnetic field component does not correspond to the position of the current source. Further, with only the isointegral map obtained from the electrocardiogram, because of the individual difference such as the position and size of the internal organs, making it difficult to accurately determine abnormality such as a disease by simply gathering from the isointegral map.

SUMMARY OF THE INVENTION

An object of the present invention is to provide biomagnetic field measuring method and apparatus which can grasp the whole state of a living body portion by using maps which are greatly reduced in number as compared to the maps required in the prior arts.

Another object of the present invention is to provide biomagnetic field measuring method and apparatus which can permit analysis of a magnetic field source by measuring a vertical component B.sub.z of a biomagnetic field without increasing the number of detection coils.

According to the present invention, (1) a biomagnetic field measuring method comprises: a first step of measuring a temporal change of a component of a biomagnetic field generated from a living body by using a plurality of fluxmeters disposed externally of the living body and each including a superconducting quantum interference device (SQUID), the magnetic field component being in a first direction which is vertical to the surface of the living body; a second step of determining a temporal change of a value proportional to a root of square sum of change rates of the first-direction magnetic field component in second and third directions which cross the first direction; a third step of integrating the temporal change of the value obtained in the second step over a predetermined interval to determine an integral value, and a fourth step of displaying the integral value obtained in the third step.

According to the present invention, (2) a biomagnetic field measuring method comprises: a first step of measuring temporal changes of components of a biomagnetic field generated from a living body by using a plurality of fluxmeters disposed externally of the living body and each including a superconducting quantum interference device (SQUID), the magnetic field components being in first and second directions which are parallel to the surface of the living body; a second step of determining a temporal change of a value proportional to a root of square sum of the first-direction and second-direction magnetic field components; a third step of integrating the temporal change of the value obtained in the second step over a predetermined interval to determine an integral value; and a fourth step of displaying the integral value obtained in the third step.

Specifically, in the biomagnetic field measuring methods (1) and (2) as above, the above integral values are used through interpolation and extrapolation to display an isointegral map for connecting points at which the integral values in the above fourth step are equal to each other, the above third step of integrating the temporal change of the value obtained in the second step over a predetermined interval to determine the integral value is carried out over a plurality of predetermined intervals to determine a plurality of integral values, and computation for determining any of the ratio, the sum or the difference between the plurality of integral values is carried out. In the Cartesian coordinate system (x, y, z), the direction vertical to the body surface is defined as z axis, the first direction is defined as z direction, the second direction is defined as x direction and the third direction is defined as y direction. In the polar coordinate system (r, .theta., .phi.), the direction vertical to the body surface is defined as r axis, the first direction is defined as r direction, the second direction is defined as .theta. direction and the third direction is defined as .phi. direction.

According to the present invention, (1) a biomagnetic field measuring apparatus for measuring biomagnetic field distribution comprises: a plurality of fluxmeters disposed externally of a living body and each including a superconducting quantum interference device (SQUID) for detecting, as a signal, a biomagnetic field generated from the living body; operation processing unit for performing the operation processing of the signal; and display unit for displaying a result of the operation processing. In the biomagnetic field measuring apparatus, the fluxmeters detect a temporal change of a component of a biomagnetic field, the magnetic field component being in a first direction which is vertical to the surface of the living body, the operation processing unit performs computation for determining a temporal change of a value proportional to a root of square sum of change rates of the first-direction magnetic component in second and third directions which cross the first direction and computation for integrating the temporal change of the value over a predetermined interval to determine an integral value, and the display unit displays the integral value.

According to the present invention. (2) in the above biomagnetic field measuring apparatus, the fluxmeters detect temporal changes of components of a biomagnetic field, the magnetic field components being in first and second directions which are parallel to the surface of the living body, the operation processing unit performs computation for determining a temporal change of a value proportional to a root of square sum of the first-direction and second-direction magnetic components and