JP2001104265A - Biological magneic field measurement equipment - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide biological magnetic field measurement equipment that performs averaging processing quickly. SOLUTION: Method for measuring biological magnetic field comprising displaying means 8-1 for the measured magnetic field using plural number of magnetic sensors by repeated measurement of the magnetic field generated from the heart of a biological body, and displaying at the same time in the same screen the measured magnetic field wave from that is obtained by indicating time interval including peak position of repeatedly measured magnetic field wave from by pre-determined magnetic sensor, performing averaging processing for each measured magnetic field wave form over the time interval including peak positions in magnetic field wave form measured by each magnetic sensor and performing averaging processing for magnetic field wave form measured by any one or plural number of magnetic sensor, and the measured magnetic field wave form that is measured by any one or plural number of magnetic sensor. Thereby, peak detection efficiency is improved.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a biomagnetic field measuring method and a biomagnetic field suitable for measuring a biomagnetic field generated by an electric current in a living body, such as a neural activity of a living brain or a cardiac myocardial activity. The present invention relates to a measuring apparatus, and more particularly to an averaging processing method for a magnetic field waveform.

[0002]

2. Description of the Related Art The distribution of a weak magnetic field generated from a living body is measured using a superconducting quantum interference device (SQUID), which is a conventional magnetic sensor, and the position of an active current inside the living body is estimated from the measurement result. A multi-channel biomagnetic imaging device for imaging the distribution is known. Such a conventional example is disclosed in, for example, Japanese Patent Laid-Open No. 4-31933.
No. 4, JP-A-5-146416, and the like.

[0003]

The above-mentioned prior art relates to the principle of operation of a biomagnetic imaging apparatus, and does not disclose any technical problems or solutions for implementation in the disclosed contents. Further, the conventional example relates to a biological activity current generated inside the brain, and does not specifically disclose other parts. An object of the present invention is to provide a biomagnetic field measurement method and a biomagnetic field measurement device that can easily measure magnetic field strengths at a plurality of measurement positions and have good operations. SUMMARY OF THE INVENTION An object of the present invention is to provide a biomagnetic field measuring apparatus that can easily display a magnetic field diagram of a living body from the magnetic field strengths at a plurality of measurement positions and that is easy to operate. It is an object of the present invention to provide a biomagnetic field measurement device including a biomagnetic field map display device that can be executed.

[0004] Another object of the present invention is to confirm whether a desired waveform appears in one or a plurality of measured magnetic field waveforms and whether the averaging process can be performed or not, and to perform a averaging process with a quick operation. It is an object of the present invention to provide a magnetic field measurement device, and in particular, to provide a biomagnetic field measurement device that improves the efficiency of peak detection. Another object of the present invention is to provide a display capable of reconfirming the result of the peak position used in the averaging process during or after the averaging process on the magnetic field waveform, thereby enabling a precise analysis of the biomagnetic field. It is to provide a device. It is another object of the present invention to provide a biomagnetic field measuring apparatus capable of measuring a magnetic field generated from a living body simultaneously with an electrocardiogram.

[0005]

SUMMARY OF THE INVENTION The present invention relates to a method for measuring a magnetic field emitted from a living body of a subject at a plurality of positions, a method for displaying a plurality of measurement positions on a display screen, and a method for displaying the displayed measurement positions. And a step of displaying information on the magnetic field strength at the selected position, wherein the plurality of measurement positions are selected in the display state of the information on the magnetic field strength. It is displayed together with the presence or absence. The present invention is characterized in that a channel item corresponding to the channel displayed on the analysis data display section is provided on the operation area display section. Another characteristic is that other measurement items are provided. The present invention specifically provides the following biomagnetic field measuring apparatus.

The present invention relates to a biomagnetic field measuring apparatus having a biomagnetic field display device for displaying a magnetic field diagram by measuring a biomagnetic field at a plurality of positions. The magnetic field diagram display device shows processing function items including measurement. A processing function display unit, an analysis data display unit for displaying a measured biomagnetic field and a processing magnetic field obtained by processing the same, and an operation area display unit for displaying operation items corresponding to the processing magnetic field for the biomagnetic field. A channel corresponding to each magnetic sensor is displayed on the analysis data display section, a selected channel item is displayed on the operation area display section, and the selected channel item is Means for selecting channel items from the displayed selected channel items, the channel items being configured as channel items associated with the channels displayed on the data display section; And, selecting the measurement item of the processing function display unit, the channel on the analysis data display unit, and the selected channel item provides biomagnetic field measuring device with a means for displaying the channel item.

[0007] Preferably, the channel and the channel item are composed of correlated rows and columns. Preferably, a means is provided for identifying a channel item consisting of a row or a column selected from the selected channel items or one channel item from another. Preferably, the pull-down menu of the measurement item includes “adjustment value file”, “fully automatic adjustment”, “offset voltage V _off adjustment”,
It has sub pull-down menus such as "manual adjustment", "adjustment value file", "measurement panel", "automatic waveform diagnosis", and "AFA offset adjustment". Preferably, measurement items for instructing start and stop are displayed on the operation area display unit, and measurement is started by a start instruction and stopped by a stop instruction. Preferably, the analysis data display unit includes an automatic waveform diagnostic means for issuing an identification signal when the display screen is out of the designated range, as an inappropriate channel for the measurement.

The present invention relates to a biomagnetic field measuring apparatus having a biomagnetic field map display device for measuring a biomagnetic field at a plurality of positions and displaying a magnetic field map as a processing result. A processing function display unit that shows processing function items including the processing biomagnetic field, and an analysis data display unit that displays a measured biomagnetic field and a processing magnetic field obtained by processing the biomagnetic field. Select and execute the “Measurement panel” submenu, and select “Single waveform display”, “Grid map display”, “Equimagnetic field diagram” and “Propagation time diagram” in the “Data analysis” menu in the analysis data display section. And a biomagnetic field measuring device for displaying the following.

In the present invention, before performing the averaging process, the digital filter or the analog filter is used to perform a filtering process to remove the fluctuation of the baseline of the magnetic field waveform due to respiration, the frequency of the commercial power source of the baseline of the magnetic field waveform (for example,
(50 Hz), and the accuracy of the averaging process can be improved. Further, in the present invention, as a method of peak detection, a point where a peak is maximum in a certain time range (maximum peak detection), a point where a peak is minimum (minimum peak detection), and a positive cross point with a threshold (from the threshold) By adopting four types of peak detection methods to detect a point of change from a small value to a large value) and a negative cross point with a threshold value (a point of change from a value larger than the threshold to a smaller value), the peak detection rate can be reduced. Improve. In the present invention, a mark is turned on at the peak or cross point of the magnetic field waveform used for the averaging process, and the appropriateness of the result of the process can be confirmed again on the display screen. Further, in the present invention, the measurement circuit body of the electrocardiogram is arranged outside the shield room, and only the electrodes and wiring portions made of a non-magnetic material for detecting the potential of the living body are arranged inside the shield, so that the magnetic noise is reduced by the electrocardiogram. It is possible to measure the biomagnetic field generated from the living body simultaneously with the electrocardiogram.

Further, a biomagnetic field measuring apparatus according to the present invention comprises a plurality of magnetic sensors for measuring a magnetic field emitted from a living body, and a display means for displaying the magnetic fields measured by the plurality of magnetic sensors. The magnetic field waveform obtained by averaging the measured magnetic field waveform and the measured magnetic field waveform in the time range including the peak position of the measured magnetic field waveform are displayed on the same screen of the display means. Having a means for storing the measured magnetic field waveform and the magnetic field waveform after the averaging process. The magnetic field waveform after the averaging process is A first mark indicating a peak position is displayed, and a second mark indicating the peak position of the measured magnetic field waveform not used in the averaging process is displayed. The time range including the peak position is displayed so as to be superimposed on the measured magnetic field waveform, and the time range including the peak position indicates a start time and an end time of the averaging process, respectively. As shown, the first and second lines intersecting the time axis of the measured magnetic field waveform, and a line substantially parallel to the time axis indicating the time range between the start point and the end point
It is also characterized in that it is superimposed on the measured magnetic field waveform and that the averaging process is performed after filtering the measured magnetic field waveform.
A biomagnetic field measuring apparatus according to the present invention is a biomagnetic field measuring apparatus comprising: a plurality of magnetic sensors for repeatedly measuring a magnetic field emitted from a living body; and display means for displaying a magnetic field measured by the plurality of magnetic sensors. A magnetic field waveform repeatedly measured by each magnetic sensor of the plurality of magnetic sensors by designating a time range including a peak position appearing in a magnetic field waveform repeatedly measured by a predetermined magnetic sensor among the plurality of magnetic sensors. The averaging process is performed for each magnetic field waveform repeatedly measured by each of the magnetic sensors over the time range including the peak position appearing in the magnetic sensor, and an arbitrary one or a plurality of magnetic sensors among the plurality of magnetic sensors perform The magnetic field waveform obtained by the averaging process on the magnetic field waveform repeatedly measured and the arbitrary A magnetic field waveform repeatedly measured by one or a plurality of magnetic sensors is displayed on the same screen of the display means, and a second position indicating a peak position appearing in the magnetic field waveform after the averaging process is performed. 1 is displayed, a second mark indicating a peak position appearing in the measured magnetic field waveform not used in the averaging process is displayed, and the time range including the peak position is displayed. It is also characterized in that it is displayed so as to be superimposed on the measured magnetic field waveform, and that the averaging process is performed after filtering is performed for each magnetic field waveform repeatedly measured by the magnetic sensors.

Further, the biomagnetic field measuring apparatus of the present invention comprises a plurality of magnetic sensors for repeatedly measuring a magnetic field emitted from a living body, and a display means for displaying the magnetic field measured by the plurality of magnetic sensors. In the device,
By specifying a time range including a peak position appearing in the magnetic field waveform measured by a predetermined magnetic sensor among the plurality of magnetic sensors, the time range including the peak position appears in the magnetic field waveform measured by each magnetic sensor of the plurality of magnetic sensors. The averaging process is performed for each magnetic field waveform measured by each of the magnetic sensors over the time range including the peak position to be measured, and the averaging process is performed by any one or a plurality of magnetic sensors among the plurality of magnetic sensors. A magnetic field waveform obtained by the averaging process on the magnetic field waveform and a magnetic field waveform measured by the arbitrary one or more magnetic sensors are displayed on the same screen of the display means.

Further, the biomagnetic field measuring apparatus of the present invention comprises a plurality of magnetic sensors for measuring a magnetic field generated from the heart of a living body,
Display means for displaying a magnetic field waveform measured by the plurality of magnetic sensors.
An analysis data display unit for displaying a magnetic field waveform obtained by processing the measured magnetic field waveform and the measured magnetic field waveform; and an operation for displaying data analysis items of the measured magnetic field waveform An operation area display section for displaying items, displaying the measured magnetic field waveform as rows or columns of a plurality of channels on the analysis data display section, and averaging a desired channel from among the plurality of channels The averaging process is performed on a magnetic field waveform measured by the plurality of magnetic sensors with reference to the designated channel as a reference channel for processing. With reference to FIG. 38, the biomagnetic field measuring apparatus of the present invention is summarized as follows. Enter a number in the text box 103a of the averaging object 103 and enter a threshold 18 for peak detection.
0 is set, and the peak detection method is set with the toggle button 103b. The averaging time width is set by inputting a number in the text box 103c.
Enter a number in d and set the time width before peak detection. Press the button 103e to perform the calculation. After peak detection, 1
The values of 03c and 103d are displayed in the waveform 106a displayed in the waveform display area 106 according to the averaging range B1, the maximum time B2, and the minimum time B3. Is displayed. The averaging process is executed by pressing the button 103f. An averaging waveform 107a is displayed in the averaging display area 107. The averaging process is performed on all the waveforms of the channels not displayed in the region 106 or the region 107 at the same time. As a result, an averaging process can be performed while confirming the measured magnetic field waveform, and a biomagnetic field measurement device capable of performing accurate analysis of the biomagnetic field with a quick operation can be realized.

[0013]

FIG. 1 shows a schematic configuration of an embodiment of a biomagnetism measuring apparatus according to the present invention. In order to remove the influence of environmental magnetic noise, the biomagnetic measurement device is installed in the magnetically shielded room 1. A subject (a subject composed of a living body) 2 is measured while lying on the bed 3.
The body surface of the subject (in the case of the chest, it is generally the plane parallel to the chest wall)
Is substantially parallel to the plane of the bed 3, and this plane is parallel to the xy plane of the rectangular coordinate system (x, y, z). Although the subject's chest is curved and inclined, it is assumed to be approximately parallel for ease of explanation. A dewar 4 filled with a liquid He as a refrigerant is arranged above the chest of the subject 2, and the dewar is a superconducting quantum interference device (SQ).
UID = Superconducting Quant
um Interference Device) and a plurality of magnetic sensors including a detection coil connected to the SQUID. Liquid He is magnetic shield 1
Is supplied continuously from the automatic replenishing device 5 outside the device.

The output of the magnetic sensor outputs a voltage which has a specific relationship with the strength of the biomagnetic field (which can also be considered as a magnetic flux density) which is generated from the subject 2 and detected by the detection coil. FLL (Flux Locked l
(op) circuit 6. This FLL circuit 6
A change in the biomagnetic field (biomagnetism) input to the SQUID is canceled via a feedback coil so that the output of the QUID is kept constant (this is called a magnetic field lock). By converting the current passed through the feedback coil into a voltage, a voltage output having a specific relationship with a change in the biomagnetic field signal can be obtained. Since the detection is performed via the feedback coil as described above, a weak magnetic field can be detected with high sensitivity. The output voltage is input to an amplifier / filter / amplifier (AFA) 7 whose output is sampled, A / D converted, and taken into a computer 8.

The computer 8 comprises a personal computer, 8-1 is a display unit thereof, and 8-2 is a keyboard.
And 8-3 indicate mice. The mouse 8-3 is used to select a processing target by moving a cursor on the screen. This operation can also be performed by operating the keyboard. The input gain (I _gain ) and the output gain (O _gain ) of the AFA 7 are adjustable, and the AFA 7 is a low-pass filter (LPF) that passes a frequency signal equal to or lower than the first reference frequency. It includes a high-pass filter (HPF) that passes a frequency signal equal to or higher than the low second reference frequency and a notch filter (BEF) that cuts off the commercial power frequency. The computer 8 can perform various types of processing, and the processing results can be displayed on the display unit 8-1. Note that the computer 8 shown in FIG. 1 shows one embodiment, and the present invention is not limited to this.
For example, a device provided with a display having a touch panel, or a device using another coordinate pointing device, for example, a trackball or a joystick, instead of a mouse, may be used.
In some cases, a computer connected via a public telephone line may be used.

As the SQUID, for example, a DC SQUID is used as an example. When an external magnetic field is applied to the SQUID, a DC bias current ( _Ibias ) is supplied to the SQUID so that a voltage (V) corresponding to the external magnetic field is generated. When the external magnetic field is represented by a magnetic flux Φ, a characteristic curve of V with respect to Φ, that is, a Φ-V characteristic curve is given by a periodic function. Prior to the measurement, an operation of adjusting the offset voltage (V _OFF ) of the FLL circuit 6 to make the DC voltage of the Φ-V characteristic curve zero level is performed. Furthermore, A
The offset voltage (A _OFF ) of the AFA 7 is adjusted so that when the input of the FA 7 is zero, the output becomes zero. When a large magnetic field is applied to the SQUID from the outside, the magnetic field is trapped on the device of the SQUID and its normal operation is not performed. In that case, the SQUID can be heated to a normal state once, and then the heating can be stopped to remove the trapped magnetic field.
In this case, the SQUID heating operation is called a heat flash.

FIG. 2 shows the arrangement of the magnetic sensors. The tangential component of the biomagnetic field (biological surface,
That is, there is a coil that detects a component that is substantially parallel to the xy plane) and a coil that detects a normal component of the biomagnetic field (a component that is orthogonal to the living body plane, that is, the xy plane). As the coil for detecting the tangential component of the biomagnetic field, two coils whose coil surfaces are respectively oriented in the x direction and the y direction are used, and as the coil for detecting the normal component of the biomagnetic field, the coil surface is z.
Oriented coils are used. Plurality of magnetic sensors 20-1 to 20-8, 21-1 to 21-8, 22-1 to 22-1
22-8, 23-1 to 23-8, 24-1 to 24-8,
25-1 to 25-8, 26-1 to 26-8 and 27-1
27-9 represent the biological surface, that is, xy, as shown in FIG.
They are arranged in a matrix on a plane substantially parallel to the plane. Although the number of magnetic sensors may be arbitrary, in FIG. 2, the number of magnetic sensors is 8 × 8 = 64 because the matrix of magnetic sensors is composed of 8 rows and 8 columns. Each magnetic sensor is
As shown in FIG. 2, the longitudinal direction is the biological surface, that is, xy.
It is arranged so as to coincide with the direction (z direction) perpendicular to the plane. In this embodiment, the bed surface and the XY surface of the sensor are parallel to each other. However, it is better to approach the body, and it is possible to tilt the bed in order to increase the measurement accuracy. However, since the human body, which is the subject, is always moving, if the human body is brought into close contact with the human body, this movement will move the detection unit, and it will be difficult to perform highly accurate detection.

[0018] Figure 3 shows the structure of a sensor for detecting the respective magnetic sensors, a normal component B _z of the biomagnetic field. In the figure, a coil made of a superconducting wire (Ni-Ti wire) is arranged so that its coil surface faces the z direction. This coil comprises a combination of two coils 10 and 11 having opposite directions. The coil 10 closer to the subject 2 is a detection coil, and the coil 11 farther away is a reference coil for detecting external magnetic field noise. You. External magnetic field noise originates from a signal source that is farther than the subject, so the noise signal is detected by both the detection coil 10 and the reference coil 11. On the other hand, the magnetic field signal from the subject is weak. Therefore, the biomagnetic signal is detected by the detection coil 10, but the reference coil 11 is hardly sensitive to the biomagnetic signal. For this reason, since the detection coil 10 detects the biomagnetic field signal and the external magnetic field noise signal, and the reference coil 11 detects the external magnetic field noise signal, the difference between the signals detected by both coils is used to obtain the S / N ratio. Measurement of a high biomagnetic field becomes possible. These coils are connected to the SQUID input coil through the superconducting wire of the mounting board mounted with SQUID12, whereby the normal direction component B _z of the detected biomagnetic field signal is transmitted to the SQUID.

[0019] Figure 4 shows the configuration of a sensor for detecting the respective magnetic sensors, the tangential components B _x and B _y of the biomagnetic field. In the figure, a planar coil is used in a sensor for detecting a biomagnetic component in a tangential direction. That is, the detection coils 10 ′ and 10 ″ and the reference coils 11 ′ and 11 ″ consist of planar coils, which are respectively arranged in first and second planes that are spaced apart in the z-direction. The two faces perpendicular to each other of these coils .4 prism connected to the input coil of the mounting substrate similarly SQUID12 'and 12 "and a normal component, the sensor 13 and By components for these B _x component detection A detection sensor 14 is attached,
This sensor capable of detecting B _x component and the B _y component is formed. The tangent components Bx and By are shown in FIG.
Besides detected using the magnetic sensor shown in, the normal component B _z obtained by the magnetic sensor of FIG. 3 x, it may be obtained by partially differentiating the y. In this case it can be detected by measuring both a single magnetic sensor tangential components B _x, and B _y and the normal component B _z is.

FIG. 5 shows the positional relationship between the magnetic sensor and the chest 30 which is the measurement section of the subject 2. The indicated points represent the intersections of the rows and columns on the matrix shown in FIG. 2, that is, the measurement points of the subject 2, that is, the measurement positions. Each of these measurement positions is also called a channel. As can be seen from the figure, in this embodiment, the height direction of the subject 2 is the y direction, and the lateral direction of the subject 2 is the x direction.

FIG. 6 shows the measurement results of the biomagnetic field at the respective measurement positions shown in FIG. This measurement result shows the time-varying biomagnetic field waveform obtained by performing the above processing based on the signals detected by the magnetic sensors corresponding to each measurement position for each corresponding channel in the matrix. Things. In this embodiment, since each channel is provided at a position where the magnetic field generated by the heart muscle can be detected, the waveform in FIG. 6 shows a magnetocardiogram waveform. A waveform obtained by measuring a magnetic field generated from the heart muscle is called a magnetocardiogram waveform. As shown in FIG. 6, when displaying the measurement data measured for each channel corresponding to the position for each channel,
This is called a grid map display. FIG. 6 shows the waveform of the measurement result of the magnetocardiogram for a certain healthy person. Here (a) shows a magnetocardiogram of tangential components B _x, (b) is a magnetocardiogram waveform of the tangential component B _y, and (c) is a normal component B _z magnetocardiogram waveform respectively.

[0022] Figure 7 shows a magnetocardiogram particular 2 channels squeezed measured tangential components B _x for a healthy person. The solid line shows the magnetocardiogram waveform of one channel, and the dotted line shows the magnetocardiogram waveform of another channel. The time zone T1 during which the ventricle of the heart is depolarized, that is, the time of each waveform peak in the systolic QRS wave is shown as t _Q , t _R , and t _S , respectively. The time zone of the T wave, which is the cardiac repolarization process (diastole), is indicated as T2. When displaying the measured data, in addition to the grid map display, single-channel data may be displayed row by row or column by column,
Data of all channels or a plurality of channels may be displayed in a superimposed manner. The former is called a single waveform display, and the latter is called a superimposed waveform display.

With respect to the obtained magnetocardiogram waveform data, the results can be displayed by performing an averaging (averaging) process or creating an isomagnetic map or a time integral diagram based on the biomagnetic field. . For example, referring to FIG.
A predetermined time (t _OFF ) goes back from the time when the rising portion of the S-wave coincides with the threshold value SL, and a time t2 when a predetermined time T3 elapses from the time t1 when the time goes back.
The data up to is added a predetermined number of times. This is averaging, and a predetermined time T3 is called an averaging time, and t _OFF is called an offset time. The magnetocardiogram waveform data may be integrated over a predetermined time range. A map formed by connecting points (channels) having the same time integration value is called an isochronous integration diagram. A map formed by connecting points having the same magnetocardiogram waveform signal value is called an isomagnetic map. Since each channel is set roughly, the interval between the isomagnetic lines, that is, the magnetic field intensity difference is set in advance, and linear interpolation is performed between the channels to draw isomagnetic lines, thereby creating a diagram more suitable for diagnosis. be able to. Referring to the magnetocardiogram waveform shown in FIG. 7, the time from time t1 to the peak position of the QRS wave (time t _R ) is called a propagation time, and a map formed by connecting points having the same propagation time is equal propagation. Call it a time diagram. The setting level of the threshold value SL can be changed. At time t1, this is
Although the time when the rising part of the QRS wave matches the threshold SL is determined as a reference, the determination may be based on the time when the falling part of the QRS wave matches the threshold SL. The time point t1 further includes the peak position time point of the QRS wave (t
_(R time point) and the time point may be determined as a reference. The measured biomagnetic signal is generated by an electrophysiological phenomenon in a living body, and its source is approximated by a current dipole model. The current dipole of the magnetic field generation source is synthesized and displayed on the isomagnetic diagram, which is called a magnetic field source display.

A series of operations from registration of a subject to measurement of data of the registered subject and analysis of the measured data are performed on a display screen displayed on a display 8-1. It is done while watching. Therefore, prior to the description of the series of operations, the layout of the display screen will be described first.

FIG. 8 shows a basic layout of a display screen displayed on the display 8-1 in FIG. The upper part of the display screen is a title bar part 80 arranged in order from the top.
1, occupied by a menu bar section 802 and a toolbar section 803 in which icons are arranged. Each of the above units can be considered as a display area or area. These arrangements are commonly displayed on the display screen for other processing purposes, such as registration and reading of a subject, measurement of a magnetic field, and analysis of measured data. This increases ease of use and reduces the time for measurement and processing. The center of the display screen is occupied by a subject information section 804, an analysis data section 805 for displaying analysis data such as a diagram and a waveform, and an operation area section 806 arranged in order from left to right. The lower portion is occupied by a status bar 807, which is a message bar portion 807-1 arranged on the left side and displaying a guide message regarding the next operation.
And a date and time display section 807-2 arranged on the right side of the display. Note that the message bar section 807-1 and the date and time display 807-2 may be one display area.

On the display screen in this embodiment, a title bar section 801 at which the name of this system is always displayed at the top, a menu bar section 802 for performing basic operations of the system, and a menu bar section 802 Since the toolbar unit 803 that enables frequently used operations is provided, the user does not need to search the operation area every time the display screen changes, and always knows the currently operating system by looking at the upper part of the display screen. be able to. Moreover, since the upper part of the display screen is the first part to look at, assuming that a person reads a sentence, it is possible to provide basic items in the operation of this system at the top. ， Improved usability in a natural way. In the center of the display screen, an analysis data display section 805, which is the main body of the display screen, is provided in the center to improve the visibility, and on the right side, an operation area section unique to the display screen. By providing the 806, the display screen and the operation area in the right-hand operation are arranged in the same manner, so that the operation can be performed without a sense of incongruity. Therefore, even if this display screen is adopted as a display screen with a touch panel, the right hand operating the operation area section 806 does not disturb the analysis data display section 805. Similarly, since the subject information section 804 having only a confirmation function is provided on the left side of the analysis data display section 805, the operation can be performed while always confirming the patient. In addition, since the left position is the farthest position in the right hand operation, even if adopted for a display screen with a touch panel, the visibility of the display is not affected.

The analysis data display section 805 and the operation area section 806 are replaced by the subject list and the data list of the subject only when they are displayed (FIG. 24). When the subject list screen (FIG. 24) is displayed on the subject information section 804, information on the subject on which the cursor is placed in the subject list in the screen is always displayed. When analysis data such as a diagram or a waveform is displayed in the analysis data portion (see FIGS. 25 to 3)
4), the information of the subject who obtained the displayed analysis data is always displayed. This makes it possible to clearly know the relationship between the displayed analysis data and the subject from which the analysis data was obtained. In this manner, on the display screen of this system, the subject information section 804 is always displayed at the fixed position (left side) of the display screen, similar to the menu bar section 802. There is no need to search for the subject information area each time it changes, and it can be known by always looking at the predetermined position (left side) on the display screen.

In the title bar portion, the name of the frame, specifically, the name "Multichannel MCG S
system "is displayed (FIGS. 24 to 34). In the operation area 806, operation elements such as buttons and text boxes are arranged. The menu bar part is a part for selecting the operation menu, and the menu is “File (F)”, “Edit (E)”, “List (L)”, “Data measurement (Q)”, “Data analysis (A)” And "Help (H)"
And are arranged according to the order of operation.

FIG. 9 shows the contents of each operation menu in the menu section on the display screen. The contents of these menus are displayed as pull-down menus by clicking the corresponding menu buttons. Therefore, when the operation menu is not required, only the keywords for calling the respective menus are compactly displayed on the menu bar. Can be set widely. And
When an operation menu is required, the keywords arranged in accordance with the operation procedure can be selected and displayed from a menu bar portion and can be operated. At this time, since the keywords are arranged according to an arrangement of characters (from left to right), the operation can be instructed in a natural manner.

The pull-down menu of "File (F)" is a page layout dialog box (not shown).
To close the “Page Layout (U)” for setting the page layout, “Preview (V)” for preview before printing, “Print (P)” for printing data, and Multichannel MCG System. The item includes “End of magnetocardiography system (X)”.
The “edit (E)” pull-down menu includes an item “delete (D)”. The measurement data of the subject in the subject list on the subject list screen (FIG. 24) is displayed in the data list on the same screen, but the item "Delete (D)" is displayed in the cursor in the data list. This is for deleting the data in which is placed, and the data must be selected in advance.
Clicking this menu opens a confirmation dialog box asking if you want to delete it. If you need to delete it, click “O” in the dialog box.
When the "K" button is clicked and the deletion is to be canceled, the "Cancel" button is clicked. The display and operation of the confirmation dialog box can reduce erroneous operations of deleting measurement data of the subject by mistake.

The “List (L)” pull-down menu includes “Registration (R)”, “List (L)”, “Delete (D)”,
It includes items of “search (S)” and “cancel (X)”. When the “Register (R)” button in the pull-down menu is clicked, a subject registration dialog box shown in FIG. 10 is opened. This is used when registering data about the subject. Items that can be registered are the registration date, the subject's ID number up to a predetermined digit, name, date of birth, height, weight, sex, disease And the subject's comment indicating the classification information of the subject. When the "Register" button in the dialog box is clicked, the data entered for registration is registered, and all of the input boxes are cleared so that re-entry is possible, and the "Cancel" button is clicked. Then, all the input boxes are cleared, and if the "end" button is clicked, the subject registration dialog box is closed. When "List (L)" in the pull-down menu is clicked, the subject list screen (FIG. 24) is displayed. “Delete (D)” in the pull-down menu is for deleting the subject on which the cursor is placed in the subject list on the subject list screen (FIG. 24), and is clicked. Before the deletion, a confirmation dialog box is displayed to confirm whether the file can be deleted in the same way as when "Delete (D)" in the pull-down menu is clicked. If deletion is required, the dialog box appears. The “OK” button is clicked, and the “Cancel” button is clicked to cancel the deletion. When the subject is deleted, all data in the data list for the subject is also deleted.

When "Search (S)" in the pull-down menu is clicked, a search dialog box shown in FIG. 11 is displayed. The subject and the data are searched, and only the searched subject and data are displayed on the subject list screen. The subjects and data to be searched are all subjects and all data. Subject name, registration date, gender, data type,
Data is retrieved by giving the measurement date, diagnosis result, comment, laboratory technician, body position, etc. as keywords. A compound search combining a plurality of the above keywords can be executed.
“Release (X)” in the pull-down menu is used to release the display of only the retrieved subjects and data and display all subjects and all data. The pull-down menu of “Data measurement (Q)” is “Adjustment value file (F)”, “Fully automatic adjustment (A)”, “V _OFF adjustment (V)”, “Manual adjustment (M)”, “Adjustment value file”. (F) "," measurement panel (P) "," automatic waveform diagnosis (D) ", and" AFA offset adjustment (O) ". When "adjustment file (F)" is clicked, a sub-pull-down menu containing "open (O)", "overwrite save (S)", and "save as (A)" is displayed. “Open (O)” in the sub pull-down menu displays the system adjustment screen (FIG. 34), opens the specified adjustment value file, and sets the system adjustment values, ie, I _bias and V _OFF in the system. Used for “Save (S)” is used to open a confirmation dialog box and overwrite the current adjustment value to the currently opened adjustment value file. “Save as (A)” is used to change the name of the current adjustment value data and save it in another adjustment value file.

When "Fully automatic adjustment (A)" is selected from the pull-down menu, the system adjustment screen (FIG. 34)
Is displayed, and the bias current I _bias and the offset voltage V _OFF of the FLL circuit 6 are automatically adjusted according to the flow of FIG. When “V _OFF adjustment (V)” is selected,
The system adjustment screen (FIG. 34) is displayed, and FIG.
Offset voltage V _{OFF of the} FLL circuit 6 according to the flow of
Is automatically adjusted. When "manual adjustment (M)" is selected, a manual adjustment dialog box shown in FIG. 12 is opened. The operator selects the channel using the scroll bar and mouse, to change the bias current I _bias and offset voltage V _OFF. If there is no problem with the entered value, click the "OK" button to set the value. If you click the "Cancel" button, the changes are invalidated and the dialog box is closed.

When the "measurement panel (P)" is clicked, a data measurement screen including a grid map (FIG. 25) is displayed. When "automatic waveform diagnosis (D)" is clicked, an automatic diagnosis dialog box shown in FIG. 13 is opened. The automatic waveform diagnosis automatically checks the amplitude, the median value, and the period of the Φ-V characteristic curve in the display screen of FIG. 34, and updates the display of the Φ-V characteristic curve when displaying or displaying the same curve. At the time of updating, the channel outside the designated range is notified to the operator as a state unsuitable for measurement. If a channel in an inappropriate state for measurement is detected, an error dialog box is opened and an error message is displayed to notify the operator. However, the Φ-V characteristic curve of the channel is displayed in another color and the The status may be indicated. When the minimum amplitude value is checked, when the difference between the maximum value and the minimum value in the Φ-V characteristic curve is smaller than the minimum amplitude value, it is determined that the state is inappropriate for measurement. When the median value is specified, it is determined that the case where the absolute value of the average of the maximum value and the minimum value is larger than the specified median value is inappropriate for measurement. When the period is specified, the measurement is not performed when the period of Φ in the Φ-V characteristic curve is out of the range specified by the first text box (lower limit) and the second text box (upper limit). It is assumed that it is in an appropriate state. To enable automatic diagnosis of the minimum amplitude, median value and period, click the check box to the left of each item with the mouse so that the X mark is displayed, and enter the desired value in the corresponding text box. You just need to enter The parameters of the automatic diagnosis input in this manner are enabled by pressing the “OK” button, and the dialog box is closed. If the "Cancel" button is pressed, the entered parameters are invalidated and the dialog box is closed. "AF
"A offset adjustment (O)" is used when adjusting the offset voltage AOFF of the AFA 7, and "A
When "offset adjustment (O)" is clicked, a data measurement screen (FIG. 25) belonging to the single waveform display is displayed.
The pull-down menu of “Data analysis (A)” is “Averaging (A)”, “Single waveform display (W)”, “Overlap waveform display (M)”, “Grid map display (G)”, “Equal magnetic field”. Diagram (B), "time integration diagram (T)", "propagation time diagram (P)", "magnetic field source display (S)", "line mode (L)", and "filling mode (F)""including.

When "Averaging (A)" is clicked, an averaging screen (FIG. 27) is displayed. When "Single waveform display (W)" is clicked, a single waveform display character screen (FIG. 2) is displayed.
8) When the “overlay waveform display (M)” is clicked, the overlay waveform display screen (FIG. 29) is displayed. When the “grid map display (G)” is clicked, the grid map display screen (FIG. 30) is displayed. A check mark (x) is displayed on the left side of the menu to indicate that the selection has been made.
Clicking on “Equimagnetic field diagram (B)” displays the isomagnetic diagram screen (FIG. 31), and clicking on “Time integral diagram (T)” displays the isochronous diagram screen (FIG. 32). , When the “Propagation time diagram (P)” is clicked, the equal propagation time diagram screen (Fig. 3
3) are displayed, and a check mark (x) is displayed on the left side of the menu to indicate that the selection has been made. Also,
When "Magnetic field source display (S)" is clicked, an inverted triangle mark is displayed on the left side of the menu, and when displaying the isomagnetic map, the magnetic field source approximated by the current dipole is superimposed and displayed (not shown) . When the line mode (L) is selected on the screen for isomagnetic diagram, isochronous integration diagram, isopropagation time diagram, and magnetic field source display, the lines are displayed without painting, and the "filling mode (F ) "Is displayed with the line gaps filled. A check mark is displayed on the left side of the menu to indicate the current mode.

The "Help (H)" pull-down menu includes "Table of Contents (C)", "Search by Keyword (S)", and "Version Information (A)". Used to search for topics by keyword and to open the version dialog box. In the toolbar 803, “subject registration” (808), “subject list” (809), “print” (810), “preview” (811), “system adjustment” (812), and “data measurement” (813) and “data analysis” (814). Although these are not shown, they are linked to the functions of the menu, so that the frequently used items from the pull-down menu can be selected. That is, "subject registration"
(808) is “Registration (R)” of “List (L)”,
The “subject list” (809) is “list (L)” of “list (L)”, and “print” (810) is “print (P)” of “file (F)” and “preview” (81
1) “Preview (V)” of “File (F)”
“System adjustment” (812) includes “manual adjustment (M)” of “data measurement (Q)” and “data measurement” (81).
3) "Measurement panel (P)" of "Data measurement (Q)"
And "data analysis" (814) respectively correspond to "grid map display (G)" of "data analysis (A)". As described above, a menu linked to an icon can be selected by simply clicking the icon. Therefore, operations that are frequently used can be easily operated simply by clicking the icon adjacent to the analysis data display section, and can be operated with icons that are easy to recognize in a shorter time than the operation of the menu bar section. The icon may be selected by the user, or may be automatically displayed on the toolbar 803 according to the frequency of use (the number of times). Next, a series of operations from registration of a subject to data measurement of the registered subject and analysis of the measured data, including the adjustment operation of the system, are shown in FIGS. This will be described with reference to FIG.

FIG. 16 shows a flow of the entire operation. When the power of the computer 8 is turned on (S-1), the operating system is started up and the program start icon is displayed on the display unit 8-1. Displayed (S-
2). From that icon, Multichannel
When the program icon of the MCG System is selected (S-3), the subject list screen shown in FIG. 24 is displayed instead (S-4). In the system according to this embodiment, a subject list screen shown in FIG. 24 is displayed as an initial screen at the start of the system. The reason for this is that the relationship between the subject and the measurement or analysis data of the subject is extremely important, and therefore, this system manages data on the subject information as keywords. That is, measurement data and analysis data cannot be managed without subject information. Therefore, in this system,
On the subject list screen, first, if the subject is registered or registered, the subject is specified. Then, in the case of a new measurement, the process shifts to measurement. To identify. It is to be noted that a display screen for waiting time at system startup may be provided prior to the subject list screen, and a display screen serving as a table of contents of the system may be provided.

The subject list screen shown in FIG. 24 will be described. The subject information section is occupied on the left side. The subject list is displayed at the top of the entire right side, and the data list is displayed at the bottom. Items displayed in the subject information section are the same as those described with reference to FIG. The items in the subject list are ID (subject ID)
No.), name, registration date (date of data registration), number of measurements (number of times data measurement was performed), date of birth, age, height, weight, comment (comment about subject)
And so on. The subject list can be scrolled with a vertical scroll bar, and the subject list items can be scrolled with a horizontal (horizontal) scroll bar. The row of the selected subject is highlighted. Items in the data list for the selected subject
ID, data type (raw (raw) data or averaging (averaging)), sampling interval (sampling interval of signal when data measurement is performed in milliseconds), sampling time (seconds), Classification (disease classification information), Date and Time (date and time at which data measurement was performed), comments (comments on data), and the like are included. The data list can be scrolled with a vertical scroll bar, and the data list items can be scrolled with a horizontal (horizontal) scroll bar. The row of the selected data is highlighted.

According to the subject list screen, one line of information of each subject is displayed on the subject list. As a result, the information of the subjects arranged vertically can be clearly separated, so that the discrimination can be improved. For example, an erroneous operation of selecting another subject by mistake can be reduced. The information of each subject can be scrolled with a horizontal (horizontal) scroll bar,
Since the information of the selected subject is arranged vertically in the vertically long subject information section for each item, the visibility is not impaired. In this case, the items in the data list of each subject are moved back and forth (left and right)
It is also possible to improve the visibility by making it movable. Further, by displaying the information of each subject on one line, many subjects can be seen at a time, so that the number of times of scrolling with the vertical scroll bar can be reduced. In addition, the data can be displayed in the lower data list by a simple operation of simply placing the cursor on the subject list from which the user wants to select data on the desired subject from the subject list and clicking the data. Moreover, since the subject list and the data list are arranged vertically, movement of the line of sight can be reduced, so that the relevance can be easily recognized. In addition, the size of the data list can be freely changed by a simple operation of moving the cursor to the upper part of the area and dragging, so that the size can be freely set according to the number of data lists. it can.

In step S-5, a desired subject line is selected from the subject list on the subject list screen. In the case of averaging processing described later, a row of raw data in the data list is always selected. Thereafter, the flow is branched into four by the menu (S-
6). According to one of the branches, the submenu “End of magnetocardiography (X)” in the menu of “File (F)” is selected, and in this case, end processing such as closing the window is performed (S-7). ), Whereby the system is shut down (S-8). After that, the power of the computer 18 is turned off.
FF is set (S-9), and all the processing ends. According to the rest of the branches, averaging processing (S-10), data analysis (S-11), and data measurement (S-12) are performed.
The averaging process can be executed by selecting a submenu “Averaging (A)” in the menu “Data analysis (A)”. The data analysis is performed by selecting “Single waveform display (W)”, “Overlap waveform display (M)”, “Grid map display (G)” in the “Data analysis (A)” menu.
It can be executed by selecting any one of the submenus of “isomagnetic diagram (B)”, “time integral diagram (T)”, “propagation time diagram (P)”, and “magnetic field source display (S)”. . Further, data measurement can be executed by selecting a submenu of “measurement panel (P)” in the menu of “data measurement (Q)”. After the completion of steps 10, 11 and 12, the flow returns to step S-4. The details of the subject selection in step S-5, the averaging process in step S-10, the data analysis in step S-11, and the data analysis in step S-12 are described below with reference to FIGS. This will be described in more detail.

FIG. 17 shows a flow of subject selection in step S-5 in FIG. In the case of subject selection, the flow is branched into four according to menu selection or subject selection. One of the branches is specified by selecting a subject (S-5).
1) This is the case where the subject selection ends. According to another branch, the "Search (S)" submenu of the "List (L)" menu is selected. by this,
A search dialog box shown in FIG. 11 is opened (S-5).
-2) The subject search conditions are input using this dialog box (S-5-3). As a result, the subject is searched (S-5-4), and based on this, the display contents of the subject list on the subject list screen shown in FIG. 24 are changed (S-5-5). According to yet another branch,
The submenu "Release (X)" of the menu "List (L)" is selected. In this case, the selected subject is returned to the list of all subjects (S-5-6), and the display contents of the subject list are changed (S-5-5). According to the remaining one of the branches, the "Registration (R)" submenu of the "List (L)" menu is selected. in this case,
The subject registration dialog box shown in FIG. 10 is opened (S-5-7), and subject information is input (S-5).
8). In these steps, the input termination is determined until all subjects have completed the input (S-5).
9) When the input is completed, the subject list is updated (S-5-5).

In this embodiment, a plurality of input data or operation instructions to be input are displayed by displaying a pull-down menu on the display screen, except for the character input such as the name and address of the subject, and the selection target is displayed. The input operation is performed by instructing the specific target with a mouse from among. As a result, most operations can be performed by mouse operation, so that a comfortable operation environment can be provided to an operator who is not accustomed to the keyboard, and the input / operation time can be reduced.
As for the plurality of selection targets of the pull-down menu, selection targets that can be input / operated in this device or its input / operation state are set and displayed in advance, so that erroneous input / erroneous operation can be reduced. Further, in this embodiment, since the cursor can be put on the input area and input can be made via the keyboard, the degree of freedom of input by the operator is secured. Further, in this embodiment, it is assumed that a character is input using a keyboard. However, it is also possible to display a dialog of the keyboard at the time of input and operate the mouse to input the dialog. Further, a handwriting input dialog may be displayed and handwriting input may be performed by operating the mouse. Further, a touch panel may be provided in the display unit, and input / operation may be performed on a screen via a fingertip or an input pen. Thus, the operability of input / operation can be remarkably improved.

FIG. 18 shows a flow of data measurement in step S-12 of FIG. First, a grid map of a magnetocardiogram waveform is displayed as an initial screen as shown in FIG. 25 (S-12-1). In the figure, channel selection, waveform monitor ON-OFF, F
Lock-unlock of the LL circuit 6, automatic adjustment of the offset voltage of the AFA 7, and operation for heat flash can be respectively performed. Further, setting of signal sampling conditions, setting of scale of waveform display, and setting of AFA parameters can be performed. It is possible. The channels consist of 64 8 × 8 channels, and all channels can be selected by clicking the “Select All Channels” button or by dragging the channel matrix diagonally from end to end. Also, by selecting a channel matrix in row units or column units or dragging one of the items (hereinafter referred to as “channel item”), a channel can be selected in row units or column units. In any case, based on the selected channel item, the magnetocardiogram waveform of the channel related to that item is displayed on the analysis data display section.
In the case of selection in units of rows or columns, the selection is displayed as shown in FIG. In this case, the selected waveform is displayed in full scale on the time axis. That is,
When all 64 channels have been selected, the analysis data display section is divided into upper, lower, left, and right (cells) as shown in FIG. As shown in FIG. 26, the analysis data section is divided into upper and lower parts to form a familiar graph form in which the time axis is left and right, thereby giving a display form giving priority to visibility.

As for the monitoring of the waveform, when the "ON" button is pressed, for example, the acquisition of the signal and the updating of the waveform are repeated every specified time from 0.5 sec to 2 sec, and the heart of the subject is monitored. The magnetic signal is monitored. Also,
When the "OFF" button is pressed, the updating of the waveform is stopped.
For FLL, click the “Lock” button or “Unlo
By clicking the "ck" button, the magnetic field lock can be performed on the 64 SQUID sensors, and the lock can be released. In that case, if one button is pressed,
The state is maintained until the other is pressed. This avoids a malfunctioning state that is not selected.
Clicking the "AFA offset adjustment" button automatically adjusts the offset voltage. Also,
If you click the "Heat Flash" button, you will see Figure 14
The heat flash operation dialog box shown in (1) is opened. Select a channel with the mouse or the arrow keys and
When the "K" button is clicked, a heat flash operation is executed for the SQUID of the selected channel. Clicking the “Cancel” button closes the dialog box and ends the process.

For the sampling time (measurement time) and interval, click the corresponding text box with an inverted triangle to open a pull-down menu of selectable numerical values, and select a desired number from the pull-down menu. it can.
The selectable numbers are, for time, for example, 1s
ec, 5 sec, 10 sec, 30 sec, 1 min, and 2 min.
sec, 0.5 msec, 1.0 msec, 2.0 ms
ec, 4.0 msec, 5.0 msec and 10.0 m
sec. Time can be from 1 sec to 2 if necessary
The selection may be made up to about 4 hours. "Time" in the "scale" box means a time scale in msec units, that is, a horizontal scale, and "signal" means a scale of the A / D converted signal, that is, a vertical scale. As for these, similarly to the selection of the sampling time and interval, a desired numerical value is selected from a pull-down menu opened by clicking the corresponding text box. The AFA parameters are input gain I _gain , output gain O _gain , low-pass filter (LPF) frequency (reference frequency), notch filter (BEF) frequency,
Includes high-pass filter frequency (reference frequency). Similarly, a desired number or character is selected from a pull-down menu opened by clicking the corresponding text box. The number or letter is I _gain
Is, for example, 1, 2, 5, 10, 20, 50,
100, 200, 500 and 1000, and O _gain is, for example, 1, 10 and 100, and LPF
About, for example, 30Hz, 50Hz, 80Hz,
100 Hz, 200 Hz, 400 Hz and 1 kHz, BEF is Off, 50 Hz and 60 Hz, and HPF is, for example, 0.05 H
z, 0.1 Hz and Thru. It is. These may be input from the keyboard instead of selection.

As the channels, in addition to the 64 channels, for example, auxiliary channels of 16 channels are prepared, and for example, an electrocardiographic waveform may be obtained in the auxiliary channels. The waveform of the tenth channel as a reference channel is displayed at the bottom of FIG. 25, which is the electrocardiographic waveform obtained on the tenth channel among the auxiliary channels. The magnetocardiogram waveform generally contains magnetic noise, while the electrocardiogram waveform does not contain such noise.
Therefore, the displayed magnetocardiogram waveform is compared with the electrocardiogram waveform of the reference channel to obtain information as to whether the magnetocardiogram waveform contains magnetic noise. Of course, the electrocardiographic waveform may be obtained from any of the regular 64 channels instead of the auxiliary channel. Also, other than the electrocardiographic waveform, an electroencephalogram, a blood flow waveform, a blood pressure waveform, or the like may be used. Further, the electrocardiographic waveform of the pregnant woman and the magnetocardiographic waveform of the fetus may be compared. As the reference waveform, not only a reference waveform of one channel but also a reference waveform of a plurality of channels may be displayed. Further, the reference channel may be used to input not only a signal from a living body but also various control signals for maintenance or the like.

Returning to the flow of FIG. 18, step S-1
In 2-2, the monitor channel is selected in the manner described above (S-12-2), and when the lock button of the FLL is pressed, the magnetic fields of all SQUIDs are locked (step S-12-3). In this state, the setting of the sampling time and the signal, which are the measurement parameters, and the setting of the AFA parameters are performed (step S-12-4). With respect to the setting, this setting can be omitted from the next time by using the set condition. As a result, it is not necessary to set the conditions every time, so that the setting time can be reduced. The setting conditions may be recorded with a name or the like so that the setting conditions can be called. "Start" in the "Measure" box
When the button is pressed, the measurement is started, "Measurement in progress" is displayed as shown in FIG. 15, and a progress bar indicating the progress of the measurement is displayed (step S-12).
5). In the progress bar of this embodiment, the bar extends in the bar graph format from left to right according to the progress of the processing. , A pie chart, etc. Further, the progress bar is displayed at a predetermined position on the measurement screen while leaving the measurement screen around. As a result, since the contents of the display screen do not change significantly, erroneous operations can be reduced.

When the measurement is started, the displayed signal waveform is fixed as it is, and the progress bar is displayed on the fixed display screen. Updating of the progress bar is repeated, for example, every second until the set time ends (S-12-6). "Cancel" in the "Measure" box
When the button is pressed, the measurement stops. When the measurement is completed, the screen shown in FIG. 26 is displayed, and the waveform is confirmed (S-12-7). Thereafter, the necessity of saving the data is determined (S-12-8), and if saving is necessary, the menu "File (F)"-"Save (S)" is selected, and the signal is saved. The data is added to the data list of the subject (S-12-9). Thereafter, it is determined whether the measurement is necessary again, including the case where the storage is not necessary (S-12-10). If necessary, the above steps are repeated. If not, all the steps of the data measurement are completed. In this case, the menu is selected so that the screen is displayed as shown in FIG. FIG. 26 shows an example in which the channel in the second column is selected as the channel.

In FIG. 26, at the bottom of the analysis data portion, there is a scroll bar 262 through which a scroll box 261 moves. The scroll box 261 is movable between the left and right ends of the scroll bar 262, and the width w of the scroll box represents a time scale. The time width between the left and right ends of the scroll bar 262 indicates the measurement time. Therefore, the displayed waveform is an enlarged waveform of a part of the waveform generated during the measurement time corresponding to the time scale w of the scroll box 261. It is. Accordingly, the operator indicates how long (the width of the scroll box 261) the waveform currently displayed in the analysis data section is within the measurement time (the width of the scroll bar 262), and the waveform indicates the measurement time. Since the first half or the second half can be grasped by the human eye, visibility can be improved. Further, since the position of the scroll box 261 and the waveform of the analysis data section are linked, by moving the cursor while dragging the cursor to the scroll box 261, the display area of the analysis data section is moved for a predetermined time. May be viewed. By doing so, the scroll bar 262 can be treated as a table of contents, for example, the waveform of a predetermined time can be easily confirmed, and the measurement content can be confirmed in detail.

FIG. 19 shows the flow of the averaging process in step S-10 of FIG. In the averaging process, in order to remove noise of each data measured in each channel, the data is added for each channel and the average value is calculated. The following processing is performed to set a reference time for the averaging processing of each channel, and the averaging processing of each channel is executed. First, raw data of the specified subject is read (S-10-1). This is performed by selecting a row in the list data of FIG. 24 in which "data type" is raw data. Thus, the magnetocardiogram waveform of the channel in the first column shown in FIG. 27 is displayed as an initial display (S-10-2). Next, a display channel is selected (S-10-3). FIG. 27 shows an example in which the channel in the second column is selected. Thereafter, the setting channel is specified (S-10-4). This is performed by clicking the text box of “channel” in the “averaging condition” box in the operation area portion or by clicking the mouse on the area displaying the waveform of the desired channel. In this case, clicking the triangle button in the text box increases the channel number, and clicking the inverted triangle button decreases the channel number. FIG. 27 shows an example in which the specified channel is the channel in the second column and second row. By this channel designation, a channel serving as a reference time for the averaging process is specified. In specifying this channel, it is preferable to select a channel having the most typical or easy-to-understand waveform. If there is no channel having a good waveform in the selected row or column, it is possible to start over from the step (S-10-3) again.

When one channel in the analysis data display section is specified, a threshold cursor 271 is displayed at the position of the designated channel as shown in FIG. Thus, the specified channel can be visually confirmed. In FIG. 27, three slider cursors 273 to 275 movable in the left and right directions are displayed at the upper part of the analysis data display unit 805. These are automatically displayed at the same time when the display screen of FIG. 27 is displayed. Is done. In step S-10-5, a threshold value, an offset time, and an averaging time, which are averaging conditions, are set by clicking corresponding text boxes in the “averaging condition” box. Click the triangle button to increase the number; click the inverted triangle button to decrease the number. The threshold is set by selecting a number representing the threshold in the "Threshold" text box. As a result, the threshold value cursor 271 automatically moves to the position corresponding to the numeral. The moving position in that case can be clearly confirmed by the cursor line. The change (selection) of the number indicating the threshold value in the “threshold value” text box and the movement of the threshold value cursor 271 are interlocked with each other. Therefore, the threshold value can also be set by moving the threshold value cursor 271. It is possible. The slider cursor 273 indicates the position (reference time) at which the rising portion of the waveform matches the threshold value set by the slider cursor 271, and the position indicated by the cursor can be clearly confirmed by the cursor line. The slider cursor 273 moves following the reference time position so as to coincide with the reference time position that changes according to the change of the threshold value. When a number representing the offset time is selected in the “Offset” text box and a number representing the averaging time is selected in the “Time” text box, the slider cursors 274 and 275 are positioned at the corresponding positions.
Moves. The movement position in that case can be clearly confirmed by the cursor line of those slider cursors. The movement of the slider cursors 274 and 275 is linked to the selection of numbers in the "Offset" and "Time" text boxes. For this reason, the offset time and the averaging time can be set by moving the slider cursor.

In step S-10-6, whether the averaging should be performed automatically or manually as the averaging mode is set (S-10-7). Then, click the “Set” button (S-10-
8), an addition process is performed. Clicking the "Cancel" button cancels all averaging conditions. As the addition processing, a time t (reference time) exceeding the threshold value is searched for in each channel (S-10).
-9), and then a waveform is displayed around time t (that is, t-50 so that time t is positioned at the center of the screen)
The waveform from msec to t + 50msec is displayed)
(S-10-10), and it is determined whether or not cancel is selected in the manual mode (S-10-).
11) If the cancel is not selected, the addition of the waveform is executed (S-10-12). Step S-10
The steps from -9 to S-10-12 are repeated by the number of additions for each channel, and the added data is divided by the number of additions (S-10-13), and the averaging process is completed. I do. Thus, FIG.
In the display screen of the averaging process shown in FIG. 5, the various conditions of the averaging process can be input / operated from both the operation area portion and the analysis data display portion 805. Set
It is possible to provide the operator with various setting methods, such as setting detailed values in the operation area section, and to shorten the time for setting conditions. In particular, in this embodiment, the analysis data display unit 8
Since a specific channel can be designated as a reference for the averaging process by specifying a specific channel from the plurality of channels 05 through the cursor, erroneous operations are reduced and operability is improved. In addition, a threshold cursor and a slider cursor are displayed near the designated channel to enable visual recognition. Further, since various conditions for the averaging process can be input / operated by using a threshold cursor or a slider cursor arranged near the channel, the movement of the operator's line of sight is reduced, and the visualization is performed in accordance with the waveform display. Erroneous operation can be reduced and operability can be improved. These setting conditions, for example, the latest setting conditions may be stored and displayed as the next setting conditions, or each setting condition may be stored with a name and stored.

Another embodiment for performing the averaging process on the display screen of the averaging process shown in FIG. 27 is shown in FIG.
5 to FIG. 42 will be described in detail. FIG. 35 shows a display screen of the averaging processing program. The display screens shown in FIGS. 35 to 42 show the analysis data section 805 and the operation area section 806 shown in FIG. 8 in detail. The screen is mainly divided into five object groups. When the object group is represented in accordance with the procedure of the analysis processing, the operation area portion includes a data read object 101, a graph display control object 102,
An averaging object 103, an all-channel overview display object 104, and a file save and end object 105 are arranged. An original waveform (a measured magnetic field waveform is referred to as an original waveform) display area 106 is arranged at an upper part of the display screen, and an averaging display area 107 is arranged at a lower part of the display screen. Hereinafter, embodiments of the operation method will be described according to the operation procedure.

An embodiment of reading a file will be described with reference to FIG. Button LOAD101 to read the file
When a is pressed, a file selection dialog 108 is displayed, the icon is moved to the position of the file name to be selected, the file name is selected, and the file selection can be set.
In this embodiment, the name is test1. The file with the dat name is selected, and is clearly displayed in the highlight 108a. To select all files, select All button 10 in selection button 108b
When 8b-1 is selected, all files are highlighted and all files can be selected. When a plurality of files are selected, a plurality of windows shown in FIG. 1 are displayed.
After selecting one or more files as described above, the user presses the OK button 108b-2 to move to the next mode. If you open the file selection dialog 108 by mistake,
The user can select the Cancel button 108b-3 to return to the display screen shown in FIG. The number of data to be read for the selected file is displayed in the data point input window 101b.
Can be set by inputting a numerical value into. Set the sampling frequency of the data of the selected file to Frequen
A selection is made with the cy toggle 101c. In this embodiment, 1k
Hz sampling frequency is selected. Freq
If the method of manual selection is cumbersome as in the case of the uency toggle 101c, for example, information on the sampling frequency at the time of measurement may be recorded in the selected file together with the measurement data, and the sampling frequency may be determined automatically.

The graph display control will be described with reference to FIG. The waveform 106a displayed in the original waveform display area 106 in FIG. 37 is the test 1. dat file in Tr
ig. The waveform of the channel selected by the ch button 102a (80 channels in this embodiment) is displayed. In this embodiment, in addition to the 64 magnetic sensors shown in FIG.
The second lead electrocardiogram waveform (FIG. 43) which has an external input of the channel and is measured simultaneously with the magnetocardiogram waveform on the 80th channel
FIG. 5 shows a display screen on which averaging processing has been performed using the waveform of an electrocardiogram. Naturally, Trig. The channel selected by the ch button 102a is not limited to the electrocardiographic waveform of 80 channels, and any one of the magnetocardiogram waveforms shown in FIG. 6 can be selected. Dsp. The channel button 102b can select a channel (80 channels in this embodiment) displayed in the averaging display area 107. In the present embodiment, the channels selected in the original waveform display area 106 and the averaging display area 107 are the same, but need not necessarily be the same. The graph display area (time area width) of the original waveform display area 106 and the averaging display area 107 can be changed by inputting numerical values into the text box width 102c and the text box offset 102d.

The text box width 102c is
Original waveform display area 106 and averaging display area 107
The time width (unit: ms) displayed on the horizontal axes 106b and 107b is set in the text box offset102.
d can set the time to skip (the unit is ms). For example, if the data width of 4000 ms is to be read by skipping 1000 ms from the data point at time 0 of the data, the value of the text box offset 102 d is set to 1000 m.
s, and the numerical value in the text box width 102c is 4000 ms. In the present embodiment, the numerical value of the text box offset102d is input as 1 ms,
Set the value of the text box width102c to 4000
An example in which ms is input is shown. Original waveform display area 1 of the present embodiment
Although 06 and the vertical axis of the averaging display area 107 are displayed with data values converted into digital data by the A / D converter, they may be displayed using the values of the intensity of the magnetic field and the potential. Further, in this embodiment, the display width of the vertical axis is automatically scaled to a range in which the maximum value and the minimum value of the data within the time width selected in the display section can be displayed.

A method for setting the peak detection will be described with reference to FIG. The peak detection is performed by the averaging object 103. First, a number is input into the text box Threshold103a to set a peak detection threshold (180 in this embodiment). Next, a detection method is set by a toggle button Method 103b. As a detection method, a method of detecting a maximum peak having a value exceeding a threshold value (peak +, described in FIG. 38) and a method of detecting a minimum peak having a value smaller than the threshold value (peak)
−), A method of detecting a cross point with the threshold when the value changes from a value smaller than the threshold to a larger value (cross +), a method of detecting a cross point with the threshold when the value changes from a value larger than the threshold to a smaller value (cross +) ( cross-) are prepared. In the present embodiment, the method (peak +) for detecting the maximum peak is selected. Further selection of the averaging time width is available in the text box FrameLe.
This can be done by writing a number in the ngth 103c. By inputting a number into the text box DataLength 103d, the time width before the detection peak can be set. To perform peak detection based on the set contents (103a, 103b, 103c, 103d), the button CALC. PEAK1
Press 03e to perform calculations. As soon as the peak detection is completed, the set FrameLength103c and Dat
The value of aLength 103d is displayed on the original waveform 106a displayed on the original waveform display area 106 on the averaging display range display section B1, the averaging display range maximum time B2,
The setting range is displayed on the original waveform 106a by the averaging display range minimum time B3.

In addition, white circles C1, C2, C3 and C4 are all displayed at the peak points extracted in the waveform on the display screen. The number 1 is displayed on the circle C1 so that the number of the detected peak can be recognized. In the present embodiment, the value of FrameLength 103c is 1
At 400, the value of DataLength 103d is set at 700. The original waveform 106a displayed on the original waveform display area 106 and the averaging display range display section B1,
If the user wants to perform the averaging process while confirming the averaging display range maximum time B2 and the averaging display range minimum time B3, the averaging process is performed by pressing the button NEXT 103f, and the averaging result is displayed in the averaging display area 107. Is displayed as an averaging waveform 107a. In FIG. 38, a circle C
1 shows an example in which the averaging process is performed by pressing the button NEXT103f after observing the waveform 1 and the averaging waveform 107a is represented by a circle C because the averaging is performed once.
The same waveform as the waveform centered at 1 is displayed. In the averaging process, the displayed original waveform display area 106 is displayed.
Alternatively, averaging processing is also performed on all waveforms of channels not displayed in the averaging display area 107 at the same time. In FIG. 38, 1 is displayed in the original waveform display area 106.
Although only one waveform is displayed in one waveform and the averaging display area 107, a plurality of waveforms may be displayed in the original waveform display area 106 and the averaging display area 107, respectively.

FIG. 39 shows a case where the second peak is not selected as the averaging data by pressing the button SKIP 103g. Further, after pressing the button SKIP 103g in FIG. 39, a text box “wid” for setting the waveform display time width which was 4000 in FIG. 38 is set.
The case where the number of th102c is changed to 8000 is shown. Therefore, the averaging waveform 107a displayed in the averaging display area 107 is the first waveform selected in FIG. 38 as it is. If the user wants to interrupt the peak detection process, by pressing the button CLEAR103h, the averaging waveform 107a displayed in the averaging display area 107 disappears, and nothing is displayed in the averaging display area 107. For the peak selected for the averaging process, the circle C1 is changed from white to a black circle, and for the peak not selected for the averaging process, the mark is changed to a double circle. Of course, the shape and color of the circles are not limited to white, black, and circles. By leaving the averaging selection state of the peak point as a history as described above, it is possible to confirm again after the averaging process is completed. It is also possible to select and remove the peak averaging waveform.

FIG. 40 shows an embodiment in which the averaging process shown in FIG. 38 or 39 is continued, and the same process is performed on the sixth detected peak. Since the text box Threshold103a is set to 180 and the toggle button Method103b is set to the method of detecting the maximum peak (peak +), the peak after the sixth one is not detected because the peak value is below 180. When a peak is detected in the next time width, the value of the text box offset102d is set to, for example, 800.
If 0 is input, the same averaging process can be performed by displaying the next 8000 ms to 16000 ms. As described above, among the circles C1 to C6, black circles C1, C3, C4, C5, and C6 indicate peaks selected as the averaging process, and are outlined with double circles. C2 indicates one not selected for the averaging process.

FIG. 41 shows a button View Average 104 in the all channel summary display object 104.
A waveform display 109 of the averaging processing result of all channels when b is pressed is shown. With this operation, the waveform 109a on which averaging processing has been performed for all channels can be displayed in accordance with the position of the sensor. Similarly, the button View Sou
By pressing the rc 104a, the waveforms of all channels before processing can be displayed.

FIG. 42 illustrates a waveform file saving process after the averaging process. To save the file after averaging, click the button SAVE1
Press 05a to display the save dialog 110.
After the display, if a new file name is to be created, the user can enter a file name in the text box 110a and press the OK button 110b to save the file under an arbitrary file name. If you want to overwrite the existing file, highlight the file name you want to use and press the button OK11
Press 0b to save the averaging result as the same file. If the save dialog 110 is displayed by mistake or the file is not saved to a file, the user can return to the display screen shown in FIG. 41 without saving the file by pressing the button Cancel 110c. When all processes are completed, click the button QUIT10
By pressing 5b, the averaging process application can be terminated.

The above is the procedure and the screen display method when the averaging process is performed. When performing the above-mentioned averaging process, if the following preprocessing is performed, appropriate averaging process can be performed. If the waveform is fluctuating due to low-frequency magnetic field noise such as respiratory noise, the signal of one channel or the signals of all channels should be passed through a high-pass filter as pre-processing for averaging. . When only one channel has passed through the high-pass filter, it is also possible to perform peak detection from the signal that has passed through the high-pass filter, and apply the above-described averaging process to the signals of all channels from the peak detection result. The type of filter used for the high-pass filter is, for example, a high-pass filter having a cutoff frequency in the range of 0.1 to 3 Hz, an IIR (infinite impulse response) filter, or an F-type filter.
It is preferable to use a digital filter such as an IR (finite impulse response) filter. In the case where commercial power supply noise (for example, 50 Hz or 60 Hz) is a problem as noise and low frequency respiratory noise is a problem, a comb-shaped notch filter (removing an integer multiple of 50 Hz or 60 Hz) is removed. Filter). Comb notch filters are 0Hz, 50H
A notch filter having a bandwidth of about 1 Hz centered on a frequency such as z, 100 Hz, 150 Hz... If the above filter processing is used as preprocessing of the averaging processing, the averaging processing with less noise can be performed.

Next, when the above averaging process is performed on an electrocardiogram, a configuration for simultaneously measuring the magnetocardiogram and the electrocardiogram will be described with reference to FIG. The electrocardiograph main body 121 is arranged outside the shield room to avoid magnetic noise. The wiring 122-1 of the electrocardiograph connected to the electrocardiograph main body 121 is inserted into the shield room, and a limb lead wire 122-2 is connected to the tip of the wiring 122-1 of the electrocardiograph. A carbon electrode 122-3 that does not easily generate magnetic noise is disposed at the end of the limb induction wire 122-2. In the present embodiment, the description has been made using the limb induction wire 122-2 and the carbon electrode 122-3, but the present invention is not limited to these. For example, a carbon electrode may be used for the chest induction wire (12 lead), The electrode portion may be formed using a magnetic material. Further, the electrocardiograph main body 121 may be arranged in a rack in which the FLL circuit 6 and the amplifier / filter / amplifier 7 are stored. The waveform of the electrocardiograph is stored as digital data on the computer 8 at the same time as the waveform of the electrocardiogram, and an arbitrary number of electrocardiogram or electrocardiogram waveforms can be displayed in real time simultaneously with the acquisition of the electrocardiogram and electrocardiogram data. . Real-time waveform display is 1
The screen is not limited to the screen of one computer 8, but can be selectively displayed on a plurality of screens (computer display, television monitor, etc.) simultaneously inside and outside the shield room.

FIG. 44 shows two embodiments showing the relationship between the arrangement of the electrocardiograph main body 121 and the biomagnetic field measuring device. FIG.
Reference numeral 4 denotes a positional relationship between the FLL circuit 6 and the amplifier / filter / amplifier 7 arranged outside the seal room 1 shown in the embodiment of the biomagnetic field measuring apparatus shown in FIG. In the embodiment shown in FIG. 44A, the FLL circuit 6 and the amplifier / filter / amplifier 7 are built in the box 124 as one unit, and the FLL circuit 6 and the amplifier / filter / amplifier 7 are turned on when the power switch 125 is turned on. Are turned on all at once. An electrocardiograph main body 121 is placed above the box 124, connected to the wiring 122-1 and introduced into the shield room. The wiring 122-1 and the wiring connected between the FLL circuit 6 and the cryostat 4 may be wired inside the shield room through the same hole provided in the wall of the shield room, or may be connected to the electrocardiograph. If magnetic noise due to electromagnetic waves or the like induced in the cable becomes a problem, the wiring 122 is passed through different holes provided in the wall of the shield room.
-1 and the wiring connected between the FLL circuit 6 and the cryostat 4 may be separately wired inside the shield room. The electrocardiograph main body 121 is provided with a roll paper 123 for printing an electrocardiogram waveform. What is printed on the roll paper 123 is not limited to the electrocardiogram waveform, but may print a magnetocardiogram waveform, or a roll paper unit dedicated to printing the electrocardiogram waveform may be installed in the box 124. FIG.
In the embodiment of FIG. 4 (b), the FLL circuit 6, the amplifier / filter / amplifier 7, and the electrocardiograph main body 121 are built in the box 124 as one unit.
11 shows a flow of data analysis in step S-11. Data analysis is intended to display various types of waveforms and diagrams to obtain the information necessary for diagnosis. Select the menu shown in Fig. 9 to selectively display various types of waveforms and diagrams. Can be displayed. That is, if "single waveform display (W)" of "data analysis (A)" is selected, the single waveform screen shown in FIG. 28 is displayed (S-11-2) and "data analysis (A)".
If "Superimposed waveform display (M)" is selected, the superimposed waveform screen shown in FIG. 29 is displayed (S-11-3), "Data analysis (A)"
If "grid map display (G)" is selected, the grid map waveform screen shown in FIG. 30 is displayed (S-11-4), and "isomagnetic diagram (B)" of "data analysis (A)" is selected. If the isomagnetic diagram screen shown in FIG. 31 is (S-11-5) and “Propagation time diagram (P)” of “Data analysis (A)” is selected, FIG.
If the isopropagation time diagram screen shown in FIG. 2 is displayed (S-11-6), and the "time integration diagram (T)" of "data analysis (A)" is selected, the isotime integration diagram screen shown in FIG. (S-11-7)
Each of them is displayed, and the system is terminated by selecting “End of magnetocardiography (X)” in “File (F)”.

When a radio button (circular button in the drawing) in the operation area is clicked on each screen, a screen of a waveform or a diagram designated by the click is displayed instead. For this reason, in FIG. 20, the branch portion is not "branch by menu" but "branch by menu or radio button". Therefore, according to this embodiment, various analysis data can be obtained only by clicking the radio button in the operation area without selecting the menu of FIG. 9, so that the operation time can be reduced and Erroneous operation can be reduced and operability can be improved.

In FIGS. 28 to 30, the “magnetic flux density” in the “scale” box is the value of the full scale on the plus side and the minus side with reference to the zero level (the unit is picotesla (pT)). The value is selected from the pull-down menu that opens when you click the triangle button in the text box. 28 to 33, "display component"
By clicking the radio box in the box, the waveform of the normal component or the waveform of the tangent component can be selected and displayed on the screen.

In FIG. 28, the waveform of each channel selected in the channel is displayed with the left end of the analysis data portion adjusted to the offset time. According to this analysis data, it is possible to compare the shapes and sizes of the waveforms of the channels arranged and displayed in the analysis data section in the vertical direction. Similarly, in FIG.
8, the waveforms arranged vertically are displayed in a superimposed manner, and the shapes and sizes of the waveforms can be compared. Also, in FIG. 30, all the channels are displayed on the basis of the offset time as in FIGS. 28 and 29. Therefore, the operator can select and analyze the number of channels as needed.

In FIG. 31, a vertically elongated magnetic field strength index box 310 is arranged at the right end of the analysis data section. The magnetic field strength index boxes have different colors from each other.
It is divided into compartments. This is to improve the visual angle (color) recognizability by distinguishing the intensity range of the magnetic field indicated by each island pattern on the isomagnetic diagram screen shown in FIG. 31 by the type of color. That is, the center position 311 in the longitudinal direction of the magnetic field strength index box 310 is a position where the magnetic field strength is zero, and the sections above the center position are referred to as the first to sixth sections in order from the center position. Then, for example, the first section is in the magnetic field intensity range of 0 to 2 pT, the second section is in the magnetic field intensity range of 2 to 4 pT, and the third section is 4 to 6 pT.
The fourth section corresponds to a magnetic field strength range of 6-8 pT, the fifth section corresponds to a magnetic field strength range of 8-10 pT, and the sixth section corresponds to a magnetic field strength range of 10-12 pT. I have. The same applies to the section below the center position. However, the section above the center position indicates the magnetic field strength in the plus direction, and the section below the center position indicates the magnetic field strength in the minus direction. The isomagnetic diagram shown in FIG. 31 is displayed in different colors according to the magnetic field strength in accordance with the definition of the correspondence between the magnetic field strength range and the color in the magnetic field strength index box 310. In addition, as a color, the plus side of the magnetic field strength may be a warm color system, the minus side may be a cool color system, and the center may be yellow. Thereby, the strength of the magnetic field can be recognized in color, so that the visibility can be improved. Moreover, according to this embodiment, the magnetic field strength index box 310 is located near the analysis data section.
Is provided, the color of the comparison object, that is, the color given to the map and the predetermined color of the magnetic field strength index box 310 can be confirmed while comparing without greatly moving the line of sight, so that the level of the strength of the magnetic field can be confirmed. And the relationship between color and color can be clearly determined. In this embodiment, the magnetic field strength index box 310 is provided at the right end of the analysis data section, but may be provided near the analysis data section, for example, at the upper, lower, or left side.

In FIG. 31, “number of maps” in the “reconstruction parameter” box indicates the number of displayed isomagnetic maps.
The “maximum value” indicates the magnetic field strength corresponding to both ends of the magnetic field strength index box 310, and the “interval” indicates the magnetic field strength index box 3.
10 means the magnetic field range corresponding to the length of each section.
Its value can be selected by clicking the triangle or inverted triangle button in the corresponding text box. The ECG waveform of the reference channel and two map time selection cursors 311 and 312 are displayed at the bottom of the analysis data display section. The two map time selection cursors 3
Divided lines having the same interval are displayed between 11 and 312, and the number of these lines matches the number of maps selected by the map number selection. Also, two map time selection cursors 3
11 and 312 can move their positions independently in the left-right direction. When the movement changes the distance between the cursors, the distance between the dividing lines also changes, but the distance between the dividing lines is always the same. Of course, it is also possible to provide a cursor for each division line and set each one separately. In FIG. 31, the number of displayed isomagnetic maps is 16, but these diagrams are diagrams at the time when the dividing line on the electrocardiogram waveform is located.
For each map, the time is also displayed so that it is possible to know when the map is at that time. Thus, in the same manner as described with reference to FIG. 26, the operator can determine how much the map currently displayed on the analysis data display section (the width of the two cursors 31) within the analysis time (the width of the electrocardiographic waveform).
1 and 312), and the range indicated by the map in the analysis time can be grasped at a glance, so that the visibility can be improved. Further, the range indicated by the map can be set by simply moving the two cursors with a mouse operation, so that the operation is easy. Furthermore, if the intervals between the dividing lines are freely set, various analysis environments can be provided to the operator, for example, by denying the questionable part and sparsely setting the other parts. When a check mark is displayed by clicking the check box of “current direction” in the operation area, an arrow is displayed on the isomagnetic diagram. The isomagnetic map on which the arrow is displayed is called an arrow map (not shown). For the arrow, the position indicates the position of the channel (position of the magnetic sensor), the length indicates the strength of the magnetic field, and the direction indicates the direction of the current when the direction of the magnetic field is converted into the direction of the current.

FIG. 32 shows an equal propagation time diagram. In this diagram, the starting position of the propagation time (t1 in FIG. 7) is shown.
The time can be changed by changing the position of the cursor 321 moving on the reference waveform.
Can be changed by dragging the cursor with a mouse.

FIG. 33 shows two time integral diagrams and one difference diagram. The "difference display" of the isochronous integration diagram can be easily displayed by clicking a check box and displaying a check mark. When "differential display" is checked, four cursors 331 to 334 appear on the reference waveform, and two isochronous integration diagrams are displayed on the upper left and right sides, and a difference diagram is displayed on the lower left side, as shown in the figure. The two isochronous diagrams show the magnetocardiogram waveform and the cursors 331 and 332 on the reference waveform.
And 100 ms to 140 ms and 180 ms set using the cursors 333 and 334, respectively.
ec to 240 msec based on the value integrated over the time range, and each time range is indicated by the cursor 3
31 and 332 and cursors 333 and 334 can be changed by dragging with the mouse, respectively. The difference diagram is 2
It represents the difference between the two time integral diagrams. If there is no check mark, only two cursors (for example, cursors 331 and 332) appear for the cursor, and only one isochronous integration diagram is displayed for the time integration diagram. Of course,
The integration time can be changed by changing the position of the cursor. As described above, according to this embodiment, the user clicks the check box for “difference display” to prompt the next operation.
Since the set of cursors is displayed, the operation time can be reduced without giving any doubt about the operation, and the time range can be easily set by moving the two sets of cursors with the mouse, so that the operability can be improved.

FIG. 21 shows a flow of system adjustment. The flow is branched into five by menu selection (S-1).
5). In this case, the currently displayed data is stored as it is (S-13), and the Φ-V characteristic curve shown in FIG. 34 is displayed (S-14). When “Fully automatic adjustment (A)” of “Data measurement (Q)” is selected, the adjustment values of I _bias and V _OFF are automatically calculated for the specified channel (S-16). The calculation will be described later),
The calculated adjustment value is set in the FLL circuit (S-1).
7), the flow returns to step S-14. When “V _OFF adjustment (V)” of “Data measurement (Q)” is selected,
V _OFF is calculated for the specified channel (S-
18) (Calculation of V _OFF will be described later), and then the flow proceeds to step S-17 described above. When "manual adjustment (M)" of "data measurement (Q)" is selected, FIG.
Is opened (S-19). When the operator inputs the adjustment values of I _bias and V _OFF for each channel, the input adjustment values are accepted (S-20), and the adjustment values are set to FL by pressing the “OK” button in the dialog box.
It is set in the L circuit (S-21). The dialog box is thereby closed (S-22), and the flow proceeds to step S-17. When the "adjustment value file (F)" of "data measurement (Q)" is selected, the flow is further branched into two by a sub pull-down menu (S-
22). That is, "Open (O)" is selected from the sub pull-down menu to inquire the file name to the operator, and the operator inputs the file name including the contents of FIG. 11 (S-23). The contents of the adjustment value file are FLL
The circuit is set (S-24). Also, "Data measurement (Q)"-"Adjustment value file (F)"-"Overwrite save"
Or, select “Save as (A)” and select step S
-23 is performed (S-25), and the adjustment value can be written to the adjustment value file (S-2).
6).

FIG. 22 shows a flow of the automatic calculation in step S-15 in FIG.

S-151: In the case of fully automatic adjustment, the following processing is performed for all SQUID channels one by one. The channel being processed is set to ch.

[0076] S-15-2: The calculated bias current and the offset voltage is assumed to be held in the memory of named I _BIAS and V _{OF F.} I _BIAS and V _OFF can hold values as many as the number of channels, and the initial values are all 0.
And The temporarily stored amplitude of the Φ-V characteristic curve is ΔV
And

[0077] S-15-3: For each channel ch, changes in steps of [Delta] I _bias to I _bias specified by the "scan parameters" box in Figure 34 the bias current I _b from 0 (scanning) is, [Phi Let I _{bias at} which the amplitude ΔV of the −V characteristic curve becomes the largest be the optimum bias current I _BIAS (ch) of the channel ch.

S-15-4 to 8: Bias current 0 to I
Until _b, the amplitude ΔV of the Φ-V characteristic curve is obtained by the following processing. The external magnetic field Φ applied to the SQUID of the channel ch is changed (scanned) from 0 to Φext specified in the “scan parameter” box in FIG. 34 in steps of ΔΦ, and the A / D converted signal is stored. determining the maximum value M _max and the minimum value V _min of the signal Te.

S-15-9-10: Here, maximum value M
If the difference between _max and the minimum value V _min is greater than [Delta] V I _BIAS (c
h) is the value of I _bias and the value of VOFF (ch) is the maximum value M
the average of the _max and the minimum value V _min, respectively as the ΔV in the difference between the maximum value V _max and the minimum value V _min. If Trip ΔV holds the previous value is greater than the difference between the maximum value M _max and the minimum value V _min.

S-15-11: I _BIAS (c) when the above processing is repeated until the bias current I _b becomes I _bias
h) and V _OFF (ch) are the optimum bias current and offset voltage of the SQUID channel ch.

S-15-12: The above processing is executed for all SQUID channels, and the fully automatic adjustment ends.

FIG. 23 shows V in step 18 of FIG.
4 shows an OFF adjustment flow.

S-18-1: In the case of VOFF adjustment, the following processing is performed for all SQUID channels one by one.

S-18-2 to 6: SQ of channel ch
The external magnetic field Φ applied to the UID is changed (scanned) from 0 to Φ _ext specified in the “scan parameter” box in FIG. 34 in steps of ΔΦ, and the maximum value V _{max of the} A / D converted signal And the minimum value _Vmin .

S-18-7: SQUID channel ch
The optimal offset voltage VOFF (ch) is the maximum value V
It is calculated as the difference between the _max and the minimum value V _min.

S-18-8: The above processing is performed for all SQUAs.
V _{OF F} adjusted by running the ID channel ends. In the above description, “scanning of the external magnetic field Φ” can be realized by giving a signal obtained by D / A conversion of a value set inside the computer to the SQUID sensor.

[0087]

According to the present invention, there is provided a method for measuring a magnetic field of a living body which is easy to measure the magnetic field strength at a plurality of measurement positions and has good operability. According to the present invention, since the channel items corresponding to the channels displayed on the analysis data display unit are provided, the measurement and display operation can be easily performed by the correspondence between the channels and the channel items. Further, according to the present invention, it is possible to improve the efficiency of peak detection and the accuracy of peak detection, and it is possible to reconfirm whether the result of peak detection processing and the result of averaging processing are appropriate. Further, according to the present invention, an electrocardiogram and a biomagnetic field can be simultaneously measured.

[Brief description of the drawings]

FIG. 1 is a schematic configuration diagram of an embodiment of a biomagnetic field measurement apparatus according to the present invention.

FIG. 2 is a perspective view showing an arrangement configuration of a magnetic sensor used in the biomagnetic field measuring apparatus of FIG. 1;

FIG. 3 is used in the biomagnetic field measurement device of FIG.
FIG. 3 is a perspective view of a single magnetic sensor that detects a normal component of a magnetic field.

FIG. 4 is used in the biomagnetic field measurement device of FIG.
FIG. 3 is a perspective view of a single magnetic sensor that detects a normal component of a magnetic field.

FIG. 5 is a diagram showing a positional relationship between a magnetic sensor and a chest of a subject in the biomagnetic field measuring apparatus in FIG. 1;

6 is a time waveform diagram of each component of a healthy person's biomagnetic field (cardiac magnetic field) measured by each magnetic sensor in the biomagnetic field measurement apparatus of FIG. 1;

FIG. 7 is a diagram showing a time waveform of a tangential component of a magnetocardiogram of two specific channels measured for a healthy person in the biomagnetic field measuring apparatus of FIG. 1;

FIG. 8 is a diagram showing a basic layout of a display screen displayed on a display unit in the biomagnetic field measuring apparatus of FIG. 1;

9 is a diagram showing an operation menu in a menu bar section of a display screen displayed on a display unit in the biomagnetic field measuring apparatus of FIG.

10 is a diagram illustrating a subject opened when “list (L)”-“registration (R)” is selected as an operation menu on a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of the registration dialog box.

11 is a search dialog box that is opened when “list (L)”-“search (S)” is selected as an operation menu on a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. FIG.

12 is a manual that is opened when “data measurement (Q)”-“manual adjustment (M)” is selected as an operation menu on a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of the adjustment dialog box.

13 is an automatic menu that is opened when “data measurement (Q)”-“measurement panel (P)” is selected as an operation menu in a display screen displayed on a display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of a diagnosis dialog box.

FIG. 14 shows a heat flash operation dialog box opened when a “heat flash” button is pressed when the display screen of FIG. 25 or FIG. 26 is displayed on the display unit in the biomagnetic field measurement apparatus of FIG. FIG.

FIG. 15 is a diagram showing a measurement progress bar displayed on a display unit during measurement in the biomagnetic field measurement device of FIG. 1;

FIG. 16 is a diagram showing a flow of an overall operation performed in the biomagnetic field measuring apparatus of FIG. 1;

FIG. 17 is a diagram showing a flow of subject selection in a subject selection step in the operation flow of FIG. 16;

FIG. 18 is a diagram showing a flow of data measurement in a data measurement step in the operation flow of FIG. 16;

FIG. 19 is a diagram showing a flow of an averaging process in an averaging process step in the operation flow of FIG. 16;

FIG. 20 is a diagram showing a flow of data analysis in a data analysis step in the operation flow of FIG. 16;

FIG. 21 is a diagram showing a flow of system adjustment performed in the biomagnetic field measurement apparatus of FIG. 1;

FIG. 22 is a view showing a flow of automatic calculation in an automatic calculation step in the system adjustment flow of FIG. 21.

FIG. 23 is a diagram showing a VOFF adjustment flow in a VOFF adjustment step in the system adjustment flow of FIG. 21.

FIG. 24 is a view showing a subject list screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 25 is a view showing a data measurement screen displayed on a display unit of the biomagnetic field measurement apparatus of FIG. 1;

FIG. 26 is a diagram showing a waveform confirmation screen displayed on the display unit of the biomagnetic field measurement apparatus of FIG. 1 when measurement is completed.

FIG. 27 is a diagram showing an averaging screen displayed on a display unit of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 28 is a diagram showing a single waveform display screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 29 is a diagram showing a superimposed waveform display screen displayed on the display unit of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 30 is a view showing a grid map display screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 31 is a view showing an isomagnetic map display screen displayed on a display unit of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 32 is a diagram showing a propagation time diagram display screen displayed on the display unit of the biomagnetic field measurement apparatus in FIG. 1;

FIG. 33 is a diagram showing a time integration diagram display screen displayed on the display unit of the biomagnetic field measurement apparatus of FIG. 1;

FIG. 34 is a view showing a system adjustment screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 35 is a diagram showing an initial screen of an averaging screen displayed on the display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 36 is a view showing a file selection screen of an averaging screen displayed on the display unit of the biomagnetic field measurement apparatus of FIG. 1;

FIG. 37 is a view showing a raw waveform display screen of an averaging screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 38 is a view showing an averaging process setting screen of the averaging screen displayed on the display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 39 is a diagram showing a screen during the averaging process of the averaging screen displayed on the display unit of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 40 is a diagram showing a screen during the averaging process of the averaging screen displayed on the display unit of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 41 is a diagram showing an entire waveform display screen of an averaging screen displayed on the display unit of the biomagnetic field measurement apparatus of FIG. 1;

FIG. 42 is a view showing a file saving processing screen of an averaging screen displayed on the display unit of the biomagnetic field measuring apparatus of FIG. 1;

FIG. 43 is a view showing the arrangement and connection of an electrocardiograph of the biomagnetic field measuring apparatus in FIG. 1;

FIG. 44 is a diagram showing an arrangement of an electrocardiograph of the biomagnetic field measuring apparatus in FIG. 1;

[Explanation of symbols]

1 ... magnetic shield room, 2 ... subject, 3 ... bed, 4
... Dewa, 5 ... Automatic replenishing device, 6 ... FLL circuit, 7 ... Amplifier / filter / amplifier, 8 ... Computer, 8-1 ... Display unit, 8-2 ... Keyboard, 8-3 ... Mouse, 2
0-1 to 20-8, 21-1 to 21-8, 22-1 to 2
2-8, 23-1 to 23-8, 24-1 to 24-8, 2
5-1 to 25-8, 26-1 to 26-8 and 27-1 to
27-8: Magnetic sensors, 10, 10 'and 10 "and 11, 11' and 11" ... Coils, 12, 12 'and 1
2 "SQUID, 13 and 14 Sensor, 30 Chest, 261 Scroll Box, 262 Scroll Bar, 271 Threshold Cursor, 273 to 275 Slider Cursor, 311, 312 and 331 to 334
... Cursor, 801, Title bar section, 802, Menu bar section, 803, Toolbar section, 804, Examinee information section, 805, Analysis data section, 806, Operation area section, 80
8 status bar section, 807-1 message bar section, 807-2 date / time display section, 808 to 814 icons, 101 ... data reading object, 101a ...
Button LOAD, 101b: Data point input window, 101c: Frequency toggle, 102: Graph display control object, 102a: Trig. c
h button, 102b... Dsp. ch button, 102c ...
Text box width, 102d: text box offset, 103: averaging object, 103a: text box Threshold
d, 103b ... Toggle button Method, 103c ...
Text box FrameLength, 103d…
Text box DataLength, 103e; Button Calc. Peak, 103f ... button NEXT,
103g ... button SKIP, 103h ... button CLEA
R, 104: all channel summary display object, 1
04a Button ViewSource, 104b Button ViewAverage, 105 File save and end object, 105a Button SAVE, 10
5b button QUIT, 106 ... original waveform display area, 10
6a: original waveform, 106b: horizontal axis, 107: averaging display area, 107a: averaging waveform, 107b
... horizontal axis, 108 ... file selection dialog, 108 a ...
Highlight, 108b ... select button, 108b-1 ... S
select All button, 108b-2 OK button, 108b-3 Cancel button, 109 waveform display, 109a waveform, 110 dialog box for saving, 1
10a text box, 110b button OK, 1
10c ... button Cancel, 120 ... gantry, 1
21: Electrocardiograph main body, 122-1: Electrocardiograph wiring, 122
-2: Limb guide wire, 122-3: Carbon electrode, 12
3 roll paper, 124 box, 125 power switch, B1 averaging range display section, B2 averaging display range maximum time, B3 averaging display range minimum time, C1, C2, C3, C4, C5 C6 ... circle.

Continued on the front page (72) Inventor Keiji Tsukada 1-280 Higashi Koikebo, Kokubunji-shi, Tokyo Inside the Central Research Laboratory, Hitachi, Ltd. (72) Inventor Hiroyuki Suzuki 882 Ma, Hitachinaka-shi, Ibaraki Pref.Hitachi, Ltd.

Claims (5)

[Claims] 1. A biomagnetic field measuring method for repeatedly measuring a magnetic field generated from a heart of a living body by a plurality of magnetic sensors and displaying a magnetic field measured by the plurality of magnetic sensors on a display means. Specifying a time range including a peak position appearing in a magnetic field waveform repeatedly measured by a predetermined magnetic sensor among the plurality of magnetic sensors; and a step of specifying the time range including the peak position appearing in a magnetic field waveform repeatedly measured by each of the plurality of magnetic sensors Performing an averaging process for each magnetic field waveform repeatedly measured by each of the magnetic sensors over the time range including the position; and repeatedly performing measurement by any one or a plurality of magnetic sensors among the plurality of magnetic sensors. A magnetic field waveform obtained by the averaging process on the applied magnetic field waveform and the arbitrary one or more magnetic fields. Biomagnetic field measuring method characterized by comprising the step of displaying a repeating measured magnetic field waveform by a sensor on the same screen of the display means. 2. The biomagnetic field measuring method according to claim 1, wherein a first mark indicating a peak position appearing in the magnetic field waveform after the averaging process is displayed. Magnetic field measurement method. 3. The biomagnetic field measuring method according to claim 1, wherein a second mark indicating a peak position appearing in the measured magnetic field waveform not used in the averaging process is displayed. A biomagnetic field measurement method characterized by the above-mentioned. 4. The biomagnetic field measuring method according to claim 1, wherein the time range including the peak position is displayed so as to overlap the measured magnetic field waveform. . 5. The biomagnetic field measurement method according to claim 1, wherein the averaging process is performed after filtering for each magnetic field waveform repeatedly measured by each of the magnetic sensors. Method.