WO2004081595A1 - Method and system for processing a multi-channel measurement of magnetic fields - Google Patents
Method and system for processing a multi-channel measurement of magnetic fields Download PDFInfo
- Publication number
- WO2004081595A1 WO2004081595A1 PCT/FI2004/000144 FI2004000144W WO2004081595A1 WO 2004081595 A1 WO2004081595 A1 WO 2004081595A1 FI 2004000144 W FI2004000144 W FI 2004000144W WO 2004081595 A1 WO2004081595 A1 WO 2004081595A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- volume
- signal
- measurement
- source
- series
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
Definitions
- the present invention relates to a new and advanced method and system for processing a multichannel measurement signal.
- the present invention relates to measurement geometry that measures a sourceless, irrotational vector field, which measurement geometry consists of a sensor arrangement; the invention also relates to the formation of a basis for signal space and the utilisation of the basis in the pre-processing and analysis of data measured by the measurement geometry.
- the representation of measurement signals or signal vectors measured by multi-channel measuring devices in some natural signal space basis is useful from the point of view of the analysis and illustration of the measured phenomenon.
- one such basis can be obtained by modeling the source that produced the magnetic field by a current distribution whose total amplitude is the smallest possible one among the distributions that are capable of producing the measured signal vector and which current distribution has been limited within a conductor area that describes the object to be examined. This kind of method is called a minimum norm estimate and it has been described e.g.
- the source distribution that produced the vector field has also been modelled by means of a multi-pole development, but in that case, the source model has usually not been formed to be the basis for signal space.
- the formation of a basis based on a multi-pole development has been presented in the doctoral thesis "Interpretation of Neuromagnetic Measurements: Modeling and Statistical Considerations", Matti Ha.mala.inen., 1987, in which the basis is designed to describe the magnetic field outside an area containing the source distribution.
- the multi-pole development has also been used for source modeling in magnetocardiography, but in that case one has not formed a basis for signal space of the development. This is apparent e.g.
- the SSP projection method in which the sub-space caused by typical external interference sources is determined and the measured signals are projected against this orthogonally to the space.
- the projection also diminishes the actual interesting signal, in case the signal vectors produced by the object to be examined are not orthogonal to the interference sub-space, and furthermore, the SSP eliminates the interference completely only in cases in which the external interferences really belong to a predetermined interference sub-space. Also the illustration of the spatial distribution of the signals suffers from the projection.
- the SSP method has been described e.g. in patent publication FI925461 and in publication "Signal-space projection method for separating MEG or EEG into components", M.A. Uusitalo and R.J.
- reference sensors are installed in the measuring device in such places with respect to other sensors that they can be considered to measure solely external interferences without observing any signal from the object to be measured.
- the signals associated with the external interferences can be calculated and reduced from the signals measured by the signal channels, in which case just the signal relating to the object to be measured is remaining. This kind of method is justified only when the reference sensors are capable of measuring all the information relating to the external interferences in the area of the signal sensors without, however, measuring the signal associated with the object at all.
- the method is based on an assumption according to which external interferences in the area of the measuring device are uniform, whereas the signals produced by the object to be measured are weakening very fast as the distance grows.
- the method based on the reference signals is described e.g. in patent publication O9641209.
- one problem that distorts the signals is the possible movement of the object to be examined during the measurement. Lately, e.g. in magnetoencephalography there have been made excitation response measurements to testees who cannot keep their heads immovable during the measurement. An apparent solution is to reject from the signals to be averaged signals that correspond to such moments of time during which the head has been situated too far from some reference point. In this kind of method, the signal-to-noise ratio weakens due to the rejection of the responses particularly when the head moves much.
- the distortion caused by a movement can be cor- rected either by taking the into account the movement in the source modeling when using a distorted averaged signal, or by making the movement correction directly to the signals prior to averaging these.
- the latter method is better in that sense that as the outcome, a signal vector in an illustrative form is obtained.
- minimum norm estimate is used by determining the source distribution associated with every response to be averaged and by calculating ' from this a signal corresponding to some fixed location of head.
- the problem with the method is the computational slowness of the minimum norm estimate, the distorting effect of the external interferences, and the fact that the possibly great distance of the object to the measurement sensors may distort the result.
- the transformation of signal vectors from one measurement geometry to another is often necessary e.g. when comparing different measurements with one another.
- the question can be e.g. about the transformation of each measured signal vector into a signal vector of such a measuring device in which the measurement sensors are situated in determinate positions on the surface of some standard object.
- the vector field must be divided into basis function components, based on which the virtual signals can be calculated by forming the basis of the vector field for the virtual set of sensors and by using estimated components to calculate the virtual signals.
- a minimum norm estimate basis has been used, in which case the aforementioned regularisation and interference problems hamper the outcome.
- the use of the minimum norm estimate for calculating virtual signals has been de- scribed e.g. in the publication "Transformation of Multichannel Magnetocardiographic Signals to Standard Grid Form", Jussi Numminen et al, IEEE Transactions on Biomedical Engineering, Vol. 42, No. 1, 1995.
- the modeling of the sources that produced the measured signal is typically implemented by parametrising the source model and by calculating the parameters' e.g. by non-linear minimisation of the smallest square sum of the error in such a manner that the model and the measured signal correspond to each other as well as possible.
- the minimisation algorithm one must perform the calculation of the pa- rameterised model for each measurement sensor specifically, which results in a big number of arithmetic operations when the number of sensors is big.
- One further problem is the difficulty of setting a good initial guesstimate for the model that would fasten the convergence of the minimisation algorithm towards the correct solution.
- the objective of the invention is to eliminate the drawbacks referred to above, or at least to significantly alleviate them.
- One specific objective of the invention is to disclose a new type of method by which the measured signal can be divided into signals associated with the object to be examined and into signals that can be classified as external interferences by means of a computationally very easy operation without the reference sensors.
- the division of the signal space into these two parts can be performed without foreknowledge of the signal sources, so the method functions correctly, although e.g. the signal sources that are outside the research object and cause interferences would change in a non-foreseeable manner.
- the signals due to the non-idealities of the measuring device can be separated from the interesting signals and interference signals because they are mainly not included in either one of the aforementioned parts of the signal space.
- the objective of the invention is to disclose a faster and more accurate method of compensation for movement which is used to correct the distortion of the signal vector caused by a movement of an object to be measured by at the same time correcting the distortion caused by external interferences without any additional calculation.
- the objective of the invention is to disclose a calculation method of virtual signals based on the fact that in the method, a vector field can be presented as components of the spatial series develop- ment of the field, which components can be limited to relate solely to the object to be examined.
- the components are used to represent the measured signal as a signal vector of any sensor geometry.
- the objective of the invention is to disclose a method of source modeling in which the source distribution of an object to be examined can be interpreted by means of the measured basis function values, and in some cases it can be solved analytically based on these values.
- the objective of the invention is to disclose a calibration method of a multi-channel measuring device that measures an irrotational vector field, which calibration method is capable of utilising any measured vector field without knowledge of the sources that produced the field.
- this method for the parameters of the system to be calibrated, such values are sought for that the measured signal mainly falls in that sub-space of signal space which is spanned by the basis functions of the irrotational, sourceless vector field.
- the invention is based on the utilisation of two different sets of vectors in the modeling of a signal that originates from two different sources.
- the requirement is that the vector field to be measured is irrotational and sourceless.
- the vectors of the signal space in question consist of vector field samples taken in different places of the physical space or of their linear combinations.
- the present invention relates to a method which is used to form such a basis in which the information contained in the whole vector field can be expressed with a number of basis functions as small as possible, in the selection of which one has taken into account the irrotationality and sourcelessness of the vector field.
- the vector field can be represented in a form enabling a more enhanced, accurate and versatile mathematical operation of the field than before in many signal processing and analytical applications.
- the invention relates to a method the basic idea of which is to limit the degrees of freedom of the signal vectors of a measuring .device that measures an irrotational, sourceless vector field based on the information on the irrotationality and sourcelessness of the field to be measured.
- the number of possible linearly independent signal vectors can be considerably reduced based on the limitations set by Maxwell's equations.
- the signal to be measured is divided into components of two different series developments in such a manner that the series developments have different convergence characteristics in that sense that the signals produced by a first source, e.g. an object to be measured, and a second source, e.g. external interference sources, are included in the convergence areas of different series developments.
- a representation is obtained for the signal vectors in which the interesting signals and external interferences have different basis vectors, the co-ordinates corresponding to which can be unambiguously solved based on the measured signal.
- the solution is performed by minimising e.g. by means of a linear operation the difference between a model based on series develop- ments and the measured signal.
- the model describes an irrotational and sourceless vector field in the area of a measuring device sufficiently perfectly
- the difference in question is due to a measurement error which results solely from the non- idealities of the measuring device.
- the method can divide the signals into three parts: into the interesting signals of the irrotational and sourceless vector field, into the signals of the irrotational and sourceless vector field that can be classified as external interferences, and into signals not included in the series development model and resulting from the non-idealities of the measuring device, such as e.g. noise and sources in the measurement area.
- the present invention has the advantage that the method functions as an automatic interference elimination method in which the measurement sensors are not separately classified into signal channels and reference channels measuring merely interference signals, instead all the channels are in the same position when calculating the basis for signal space.
- the automatism of the interference elimination means that the method can model external interferences irrespective of the structure of the interference field in relation to the structure of the field produced by the object to be examined, or of the variation of the structure of the interference field with in respect of time, in contrast e.g. to the projection method in which the aforementioned signals are projected against some predetermined interference sub-space into an orthogonal space.
- the signal produced by an object to be examined weakens, if it is not orthogonal to the determined interference sub-space. This is eliminated in the present invention.
- the method in accordance with the invention enables one to correct the signal distortions that are due to the movement of an object to be examined faster and more accurately than by previous methods, which is a basic requirement of a real-time correction for movement.
- the correction for movement is performed by representing the movement of an object to be examined as a movement of the measurement sensors in the vicinity of a stationary object in such a manner that the variable geometry is taken into account when calculating the basis for signal space. In that case, the series development components corresponding to the object to be examined always are commensurable irrespective of the position of the object.
- the advantage of the invention is based on the fact that the calculation of the series development components requires quite few computational operations and on the fact that the external interferences are modelled when calculating the values of the series devel- opment components.
- this enables one to quickly calculate dependable values for the series development components corresponding to an interesting source area, although the distance and position of the object with respect to the measurement sensors would considerably vary and although the field structure of the external interferences would vary during the measurement .
- the series development components describe the field produced by the sources everywhere in the measurement volume, they can be used to calculate a signal which would be measured by any fictitious sensor disposed in the same space.
- Virtual signals of this kind are useful e.g. when there is a wish to represent the signal produced by an object to be examined always in the same geometry irrespective of the fact of where the object was really situated during the measurement in relation to the measurement sensors.
- the proportion of measurement errors is smaller than in original signals, as in the method one has attenuated the portion that is not included in the set model of an irrotational, sourceless vector field.
- the invention has the advantage of separating external interferences from the signals of an object to be examined.
- the method in accordance with • the invention enables one to identify individual sensors the signal measured by which clearly differs from the series development model that is based on the signals of the whole set of sensors.
- the method enables even a real-time quality control of measurement channels and an automatic rejection based on the signal analysis.
- the co-ordinates of the signals with respect to a basis in accordance with the invention can be re-calculated after the rejection of channels that operate in a non-desired manner, which enables one to obtain a more accurate result for the co-ordinates.
- the calibration of a measuring device or of a measurement geometry becomes easier than before.
- the calibration method in accordance with the invention does not need a predetermined calibration source, instead it suffices that the measuring device is known to measure a sourceless and irrotational vector field.
- the calibration is performed by setting the calibration coefficients and geometric parameters of the measuring device so that the measured signal vectors fit in the sub-space spanned by the signal vectors that are determined by the series developments.
- the calibration is performed by minimising the difference between the measured signal vectors and the sub-space calculated for the signal vectors that were produced by an irrotational vector field, in contrast to the conventional methods which minimise the difference between the signal vectors and the signal vectors produced by a source that is known beforehand as accurately as possible.
- the measured signals it is possible to use e.g. any signals produced by the sources in the environment, in which case the calibration can be performed completely automatically without any preparations of the measurement situation.
- the method is more accurate and dependable than the previous ones as the inadequacy of the information concerning the geometry of an object to be measured does not result in an error of the calibration.
- the invention can also be utilised in source modeling because the field components calculated irrespective of the external interferences and corresponding to the target area to be examined represent the continuous current distribution of the target area as a finite number of components so that each component describes some feature concerning the shape of a current distribution.
- the components can be used to check the hypotheses concerning the current distribution by comparing the modelled components to the theoretically calculated components of the distributions in accordance with the hypotheses. In some cases, this kind of comparison is easier than when using as the measurement result, the measured signal vector values, as it is done in several conventional methods. In some cases, it is possible to analytically calculate from the series development components the geometry of the current distribution without nonlinear minimisation algorithms, which makes the determination of the distribution clearly faster than in the previous methods .
- the invention also provides a new way of designing the geometry or set of sensors of a measuring device by optimising the capability thereof to dependably measure the series development components of an interesting object up to an order as high as possible by maximising at the same time the difference of the signal spaces associated with the object of the measuring device in question and the external interferences, such as e.g. the angle between the spaces.
- the external interferences such as e.g. the angle between the spaces.
- FIG. 1 shows one generalised embodiment of the present invention
- Fig. 2 is a schematic representation illustrating one embodiment of the present invention.
- Fig. 1 the space has been schematically divided into five areas by means of concentric spherical shells.
- the sensors (si, s2, s3, s4, s5, s ⁇ ) of the multichannel measuring device are situated in area 3, which is left in between the second and third innermost spheres.
- the sensors are connected to an analysis device 10.
- area 2 there is a current loop I in which represents an interesting source of magnetic field situated in the measurement object.
- a current loop situated in area I out represents a source that produces external magnetic interferences.
- the representations of the field produced by these two sources are examined by such series developments the development point of which is the common midpoint of the spheres inside area 1.
- the field produced by source I in can be represented as a Taylor's development which is convergent in area 1, which thus is such an area each point of which is closer to the development point than the closest source point.
- the field produced by source I. n can also be represented as a convergent, asymptotic series development outside area 2, i.e. in an area in which all the points are farther from the development point than any source point.
- An asymptotic series development is herein used to mean e.g. a series development that is based on spherical harmonics, the terms of which are propor- tional to the negative integer powers r "n of the distance r calculated from the development point.
- the same remarks also apply to a magnetic field produced by an external interference source I out : the Taylor ' s series development thereof is convergent everywhere in areas 1, 2 and 3, whereas the asymptomatic series development is convergent in area 5.
- the signals coming from an interesting object can be separated from interferences coming from external sources based on the fact that the external interferences are mainly shown in the sub-space spanned by the basis vectors that are based on the terms of the Taylor's series, whereas the interesting signals are included in the sub-space determined by the basis vectors associated with the asymptomatic series.
- Fig. 2 Let us think about a situation in which from the space, a point f 0 is chosen in such a manner that the smallest and biggest distance of the measurement points from this are r w ⁇ n ja r max , respectively.
- the measurement points are designated in Fig. 2 as points on a curved surface.
- the vector field produced by source 1 can be expressed in the area of the measurement points as a development of the M series with respect to point r 0 , but not as a development of the T series, and correspondingly, the vector field produced by source 2 can be expressed in the area of the measurement points as a T development with respect to point r 0 , but not as an M development.
- the signal vectors produced by sources 1 and 2 can be separated from one another, if the measured signal vector is presented in a basis whose basis functions correspond to the independent components of the T and M development with respect to the chosen development point. It must be noted that the development point does not need to be the same for different series developments and that there are several different alternatives for the developments.
- the aforementioned series development components form a linearly independent basis for signal space, then the interesting signals can be separated from the interferences independently.
- regularisation in case the susceptibility to interference of the basis is big.
- the susceptibility to interference is used to mean the product of the norms of the matrix and of its reverse matrix.
- the linear independence and the susceptibility to interference of the basis depend on the geometry of the measurement arrangement and on the number of the components to be modelled in relation to the number of measurement points.
- the invention provides a basis for the signal vectors of an irrotational vector field that is more practical than before and in which the regularisation means substantially the selection of the number of series development components to be included. This is easy to solve because of the fast attenuation of the components as the order of the development increases.
- the following describes the division of the measured signal vectors into an object to be examined and into external interferences using as the example the Taylor's series development as well as the multi- pole development of the spherical harmonics.
- the vector field being measured functions the magnetic flux density.
- the development point for the series developments, can be selected in such- a manner that the sources associated with the object to be measured are situated inside an r radial sphere, and the sources associated with the external interferences are situated outside an R radial sphere when the origin of both of the spheres is the development point in question.
- the Taylor's series development is convergent inside an R radial sphere when R is the distance from the development point of the series to the closest singular point, herein to a source.
- the signal vector ⁇ 0M produced by a source s out situated at a distance r out ° R can be expressed with respect to the development point by means of the Taylor series of the magnetic field, but the Taylor development is not convergent for a signal vector ⁇ in produced by a source s in situated at a distance r in • r.
- ⁇ in it is possible to use e.g. the multi -pole development.
- the signal measured on channel k is the surface integral of a normal component of the magnetic flux density over the measurement area.
- the magnetic scalar potential is expressed by means of so-called even and odd spherical harmonics.
- the magnetic flux density is the negative gradient of the potential in question, so herein ⁇ r ' can be expressed in the form:
- angles and ⁇ are angles of bearing of the spherical co-ordinates.
- the basis vectors corresponding to the multi-pole development can be formed by the same prin- cipe as in the case of the Taylor basis, whereby a
- the signal vectors associated with an object to be examined are of the form :
- This equation describes a signal vector measured by an ideal measuring device from an irrotational and sourceless vector field.
- a non-ideal measuring device produces into the model of an irrotational and sourceless vector field an inaudible meas- urement error , so by practical measuring devices the foregoing equation has to be written in the form:
- the series development components can be solved by minimising the difference between the measured signal and the model by the method of the smallest square sum, which in this case corresponds to the multiplying of the measured signal vector by the pseudo inverse of the matrix S :
- the solution of an ideal case is unambiguous, in case the dimension of the measured signal vector ⁇ is bigger than the number of vector fields of the basis S and in case the basis S is linearly independent.
- These requirements are typically realised by present- day multi-channel measuring devices, if the limitations set by the Maxwell's equations are taken into account.
- a sufficiently big number of series development components must be included in the basis S .
- the measured signal vector can always be divided into the part produced by an object to be measured, into the part produced by external interferences, and into the part not included in the irrotational vector field (measurement error) :
- the field of application of the invention is wide as in the method, from the measured signal vector, a natural representation is calculated containing - , ' ) substantially the same information as the signal vector, but in a form independent of the sensor geometry and typically in a lower-dimensional form than the measured signal vectors. Furthermore, the components associated with an object to be examined and with the external interferences are separate, unlike in the signal vector, which is the sum of an interesting signal and of interferences .
- interference sources external of the testee 's head are associated with the measurement situation, the structure of the magnetic field produced by which changes as a function of time.
- Interference sources of this kind can be e.g. the magnetic particles in a moving testee 's clothes and the testee ' s heart.
- the magnetic field measurements of the brain are usually made as repetition measurements in which the response produced by the brain and relating to the desired stimulus is averaged several times in a row to improve the signal-to-noise ratio.
- the averaged signal vector to be obtained as the outcome is the mean value of the signal vectors produced by a source that was situated in different places and is physiologically invariable from one stimulus to another, and thus the average signal cannot be used in the analysis as such.
- the elimination of external interferences by a projection method based on a predetermined interference sub-space is out of the question, since the field structure of the interferences has varied in a non- foreseeable manner.
- interference elimination based on separate reference sensors is not justified either, since the testee can time to time be far from both the actual signal and reference sensors, whereby a division- into an interesting signal and interferences that is based on shear uniformity of the field or on a signal amplitude is not very successful.
- the interference elimination can be performed by the series development method described above by averaging the component estimates associated with an interesting object, in the calculation of which one has taken into account the place and position of the testee with respect to the measuring device at different moments of time .
- the motion of a testee is modelled in such a manner that the motion of the head is represented as a motion of a set of sensors in the environment of an immovable head.
- basis matrixes S m (t) and S ⁇ (t) are calculated, from which a basis S (t) is formed, based on which it is possible to calculate, by the pseudo inverse, an estimate x(t) in accordance with equation 15.
- One analysis method includes the use of multi-pole components in source modeling, which is particularly useful when the source distribution asso- ciated with an object to be examined can be solved analytically from the estimated components. This is successful at least in a case in which the field is produced by one current dipole inside a spherical volume conductor.
- the outcome can also be represented as a virtual signal vector of any sensor geometry by calculating the basis S m of the sensor geometry in question, whereby the virtual signal vector is obtained in accordance with equation 16 by using an averaged estimate x m .
- the geometric parameters and the calibration coefficients are known as accurately as possible. In the method of the invention, these are found out by solely utilising the information on the fact that the measuring device measures an irrotational vector field. In that case, each measured signal vector shall be included in the sub-space S, and systematic deviations from this condition refer to an incomplete piece of information on the calibration and geometry of the measuring device.
- the calibration can be performed by setting the calibration coefficients and the geometric parameters so that the sub- space S calculated based on these corresponds to the sub-space spanned by the measured signal vectors.
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
Description
Claims
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006505618A JP5014783B2 (en) | 2003-03-14 | 2004-03-11 | Method and device for processing multi-channel measurements of magnetic fields |
AT04719479T ATE556328T1 (en) | 2003-03-14 | 2004-03-11 | METHOD AND SYSTEM FOR PROCESSING A MULTI-CHANNEL MEASUREMENT OF MAGNETIC FIELDS |
CA2518674A CA2518674C (en) | 2003-03-14 | 2004-03-11 | Method and system for processing a multi-channel measurement of magnetic fields |
EP04719479A EP1606639B1 (en) | 2003-03-14 | 2004-03-11 | Method and system for processing a multi-channel measurement of magnetic fields |
US11/220,415 US7463024B2 (en) | 2003-03-14 | 2005-09-07 | Method and device for processing a multi-channel measurement of magnetic fields |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FI20030392 | 2003-03-14 | ||
FI20030392A FI115324B (en) | 2003-03-14 | 2003-03-14 | A method and system for processing a multichannel measurement signal |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/220,415 Continuation US7463024B2 (en) | 2003-03-14 | 2005-09-07 | Method and device for processing a multi-channel measurement of magnetic fields |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2004081595A1 true WO2004081595A1 (en) | 2004-09-23 |
Family
ID=8565816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/FI2004/000144 WO2004081595A1 (en) | 2003-03-14 | 2004-03-11 | Method and system for processing a multi-channel measurement of magnetic fields |
Country Status (7)
Country | Link |
---|---|
US (1) | US7463024B2 (en) |
EP (1) | EP1606639B1 (en) |
JP (1) | JP5014783B2 (en) |
AT (1) | ATE556328T1 (en) |
CA (1) | CA2518674C (en) |
FI (1) | FI115324B (en) |
WO (1) | WO2004081595A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005078467A1 (en) * | 2004-02-13 | 2005-08-25 | Elekta Ab (Publ) | A method for interference suppression in a measuring device |
WO2006114473A1 (en) * | 2005-04-28 | 2006-11-02 | Elekta Ab (Publ.) | Method and device for interference suppression in electromagnetic multi-channel measurement |
WO2012004458A1 (en) * | 2010-07-06 | 2012-01-12 | Elekta Ab (Publ). | Method for adjusting interference signal space in bio-magnetic field measurements |
US10307105B2 (en) | 2012-01-24 | 2019-06-04 | Elekta Ab | Method and device for recognizing and removing undesired artifacts in multichannel magnetic field or electric potential measurements |
WO2021260272A1 (en) | 2020-06-25 | 2021-12-30 | Megin Oy | Magnetoencephalography apparatus and method |
US11213255B2 (en) | 2017-07-03 | 2022-01-04 | Ricoh Company, Ltd. | Information processing device, information processing method, and recording medium storing information processing program |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5361131B2 (en) * | 2007-01-03 | 2013-12-04 | エレクタ アクチボラゲット | Analysis of multi-channel measurement data using orthogonal virtual channels |
EP2456357B1 (en) * | 2009-07-20 | 2014-09-10 | Koninklijke Philips N.V. | Apparatus and method for influencing and/or detecting magnetic particles |
US20140355769A1 (en) | 2013-05-29 | 2014-12-04 | Qualcomm Incorporated | Energy preservation for decomposed representations of a sound field |
US9922656B2 (en) | 2014-01-30 | 2018-03-20 | Qualcomm Incorporated | Transitioning of ambient higher-order ambisonic coefficients |
US10770087B2 (en) | 2014-05-16 | 2020-09-08 | Qualcomm Incorporated | Selecting codebooks for coding vectors decomposed from higher-order ambisonic audio signals |
US20180306874A1 (en) * | 2016-01-07 | 2018-10-25 | University Of Florida Research Foundation, Inc. | Constructing a map of a multidimensional field using measurement data from one dimension of the field |
US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
WO2019133997A1 (en) | 2017-12-31 | 2019-07-04 | Neuroenhancement Lab, LLC | System and method for neuroenhancement to enhance emotional response |
US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
WO2020040168A1 (en) * | 2018-08-22 | 2020-02-27 | 旭化成エレクトロニクス株式会社 | Magnetic field measurement device, magnetic field measurement method, and magnetic field measurement program |
WO2020056418A1 (en) | 2018-09-14 | 2020-03-19 | Neuroenhancement Lab, LLC | System and method of improving sleep |
CN109725268B (en) * | 2018-12-10 | 2021-01-05 | 北京卫星环境工程研究所 | Efficient measurement system and measurement method for spacecraft magnetic information |
WO2020138170A1 (en) * | 2018-12-26 | 2020-07-02 | 旭化成エレクトロニクス株式会社 | Magnetic field measuring device |
US11497425B2 (en) | 2019-03-08 | 2022-11-15 | Asahi Kasei Microdevices Corporation | Magnetic field measurement apparatus |
US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
US11454679B2 (en) * | 2020-01-20 | 2022-09-27 | Asahi Kasei Microdevices Corporation | Magnetic field measuring apparatus, magnetic field measuring method and recording medium with magnetic field measuring program recorded thereon |
JP7525297B2 (en) | 2020-05-08 | 2024-07-30 | 旭化成エレクトロニクス株式会社 | Magnetic field measurement device, magnetic field measurement method, and magnetic field measurement program |
JP2022111838A (en) * | 2021-01-20 | 2022-08-01 | 旭化成エレクトロニクス株式会社 | Magnetic field measurement device, magnetic field measurement method, and magnetic field measurement program |
JP2023131647A (en) * | 2022-03-09 | 2023-09-22 | 株式会社リコー | Magnetic measurement device, magnetic measurement processing device, and control method of magnetic measurement processing device |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4492923A (en) * | 1982-06-21 | 1985-01-08 | The United States Of America As Represented By The Secretary Of The Navy | Apparatus for measuring the spatial scalar variation of a magnetic field with vector magnetic sensors on a moderately stable moving platform |
US4977896A (en) * | 1989-05-26 | 1990-12-18 | Biomagnetic Technologies, Inc. | Analysis of biological signals using data from arrays of sensors |
US5136242A (en) * | 1990-02-22 | 1992-08-04 | Siemens Aktiengesellschaft | Method for measuring the field pattern of electrical or magnetic fields emanating from a biological power source using a sensor arrangement for obtaining chronological map of the position of the biological power source |
EP0560621A2 (en) * | 1992-03-13 | 1993-09-15 | General Electric Company | System and method for noiseless measurement of a biomagnetic field |
FI925461A (en) | 1992-11-30 | 1994-05-31 | Risto Juhani Ilmoniemi | Method and apparatus for distinguishing between brain excitation responses and spontaneous function and different components of heart-measured signals |
WO1996041209A1 (en) | 1995-06-07 | 1996-12-19 | Ctf Systems Inc. | Method and system for obtaining higher order gradiometer measurements with lower order gradiometer |
US5644229A (en) * | 1994-11-07 | 1997-07-01 | U.S. Philips Corporation | Method of determining the spatial field distribution |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5020538A (en) * | 1988-08-12 | 1991-06-04 | Sam Technology, Inc. | Low noise magnetoencephalogram system and method |
JPH06189927A (en) * | 1992-10-19 | 1994-07-12 | Hitachi Ltd | Device and method for presumption of bioelectric current distribution |
JP3291858B2 (en) * | 1993-08-26 | 2002-06-17 | 株式会社島津製作所 | Magnetic detector |
JP2833506B2 (en) * | 1995-02-22 | 1998-12-09 | 株式会社島津製作所 | Magnetic detector |
FI112399B (en) * | 1998-08-28 | 2003-11-28 | Neuromag Oy | Method and apparatus for eliminating background noise from multichannel detector devices |
JP4068763B2 (en) * | 1999-06-21 | 2008-03-26 | 株式会社島津製作所 | Biological signal measuring device |
-
2003
- 2003-03-14 FI FI20030392A patent/FI115324B/en active IP Right Grant
-
2004
- 2004-03-11 EP EP04719479A patent/EP1606639B1/en not_active Expired - Lifetime
- 2004-03-11 WO PCT/FI2004/000144 patent/WO2004081595A1/en active Application Filing
- 2004-03-11 JP JP2006505618A patent/JP5014783B2/en not_active Expired - Lifetime
- 2004-03-11 CA CA2518674A patent/CA2518674C/en not_active Expired - Lifetime
- 2004-03-11 AT AT04719479T patent/ATE556328T1/en active
-
2005
- 2005-09-07 US US11/220,415 patent/US7463024B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4492923A (en) * | 1982-06-21 | 1985-01-08 | The United States Of America As Represented By The Secretary Of The Navy | Apparatus for measuring the spatial scalar variation of a magnetic field with vector magnetic sensors on a moderately stable moving platform |
US4977896A (en) * | 1989-05-26 | 1990-12-18 | Biomagnetic Technologies, Inc. | Analysis of biological signals using data from arrays of sensors |
US5136242A (en) * | 1990-02-22 | 1992-08-04 | Siemens Aktiengesellschaft | Method for measuring the field pattern of electrical or magnetic fields emanating from a biological power source using a sensor arrangement for obtaining chronological map of the position of the biological power source |
EP0560621A2 (en) * | 1992-03-13 | 1993-09-15 | General Electric Company | System and method for noiseless measurement of a biomagnetic field |
FI925461A (en) | 1992-11-30 | 1994-05-31 | Risto Juhani Ilmoniemi | Method and apparatus for distinguishing between brain excitation responses and spontaneous function and different components of heart-measured signals |
US5644229A (en) * | 1994-11-07 | 1997-07-01 | U.S. Philips Corporation | Method of determining the spatial field distribution |
WO1996041209A1 (en) | 1995-06-07 | 1996-12-19 | Ctf Systems Inc. | Method and system for obtaining higher order gradiometer measurements with lower order gradiometer |
Non-Patent Citations (9)
Title |
---|
CARL FRIEDRICH GAUSS, WILHELM WEBER, RESULTATE AUS DEN BEOBACHTUNGEN DES MAGNETISCHEN VER- EINS IM JAHRE 1838 |
GULLBERG G.T. ET AL.: "Backprojection filtering algorithms for reconstruction of vector and second order tensor fields", NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD, IEEE, 2000, 15 October 2000 (2000-10-15) - 20 October 2000 (2000-10-20), pages 15/277 - 15/281, XP010557058 * |
JUSSI NUMMINEN ET AL.: "Transformation of Multichannel Magnetocardiographic Signals to Standard Grid Form", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, vol. 42, no. 1, 1995, XP000556782, DOI: doi:10.1109/10.362916 |
K. UUTELA ET AL.: "Detecting and Correcting for Head Movements in Neuromagnetic Measurements", NEUROIMAGE, vol. 14, 2001, pages 1424 - 1431 |
M. BURGHOFF ET AL.: "Comparability of Measurement Results Obtained with Multi-SQUID- Systems of Different Sensor Configurations", IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY, vol. 7, no. 2, 1997, XP011501200, DOI: doi:10.1109/77.622136 |
M.A. UUSITALO, R.J. ILMONIEMI: "Signal-space projection method for separating MEG or EEG into components", MEDICAL & BIOLOGICAL ENGINEERING & COMPUTING, vol. 35, 1997, pages 135 - 140 |
MATTI HAMALAINEN ET AL.: "Magnetoencephalography - theory, instrumentation, and applications to noninva- sive studies of the working human brain", REVIEWS OF MODERN PHYSICS, vol. 65, no. 2, 1993 |
MATTI HAMALAINEN, INTERPRETATION OF NEUROMAGNETIC MEASUREMENTS: MODELING AND STATISTICAL CONSIDERATIONS, 1987 |
TERENCE SABAKA ET AL.: "A comprehensive model of the quiet-time, near-Earth magnetic field: phase 3", GEOPHYSICAL JOURNAL INTERNATIONAL, vol. 151, 2002, pages 32 - 68 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7649351B2 (en) | 2004-02-13 | 2010-01-19 | Elekta Ab (Publ) | Method for interference suppression in a measuring device |
WO2005078467A1 (en) * | 2004-02-13 | 2005-08-25 | Elekta Ab (Publ) | A method for interference suppression in a measuring device |
JP2007522465A (en) * | 2004-02-13 | 2007-08-09 | エレクタ アクチボラゲット(パブル) | Method for suppressing interference in measuring apparatus |
JP4669847B2 (en) * | 2004-02-13 | 2011-04-13 | エレクタ アクチボラゲット(パブル) | Method for suppressing interference in measuring apparatus |
US7933727B2 (en) | 2005-04-28 | 2011-04-26 | Elekta Ab (Publ) | Method and device for interference suppression in electromagnetic multi-channel measurement |
EP1880226A1 (en) * | 2005-04-28 | 2008-01-23 | Elekta AB (publ) | Method and device for interference suppression in electromagnetic multi-channel measurement |
WO2006114473A1 (en) * | 2005-04-28 | 2006-11-02 | Elekta Ab (Publ.) | Method and device for interference suppression in electromagnetic multi-channel measurement |
EP1880226A4 (en) * | 2005-04-28 | 2014-11-05 | Elekta Ab | Method and device for interference suppression in electromagnetic multi-channel measurement |
WO2012004458A1 (en) * | 2010-07-06 | 2012-01-12 | Elekta Ab (Publ). | Method for adjusting interference signal space in bio-magnetic field measurements |
US9642554B2 (en) | 2010-07-06 | 2017-05-09 | Elekta Ab (Publ) | Method for adjusting interference signal space in biomagnetic field measurements |
US10307105B2 (en) | 2012-01-24 | 2019-06-04 | Elekta Ab | Method and device for recognizing and removing undesired artifacts in multichannel magnetic field or electric potential measurements |
US11213255B2 (en) | 2017-07-03 | 2022-01-04 | Ricoh Company, Ltd. | Information processing device, information processing method, and recording medium storing information processing program |
WO2021260272A1 (en) | 2020-06-25 | 2021-12-30 | Megin Oy | Magnetoencephalography apparatus and method |
Also Published As
Publication number | Publication date |
---|---|
FI20030392A (en) | 2004-09-15 |
US20060031038A1 (en) | 2006-02-09 |
CA2518674C (en) | 2019-04-16 |
CA2518674A1 (en) | 2004-09-23 |
FI20030392A0 (en) | 2003-03-14 |
JP5014783B2 (en) | 2012-08-29 |
FI115324B (en) | 2005-04-15 |
ATE556328T1 (en) | 2012-05-15 |
EP1606639A1 (en) | 2005-12-21 |
JP2006520892A (en) | 2006-09-14 |
US7463024B2 (en) | 2008-12-09 |
EP1606639B1 (en) | 2012-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7463024B2 (en) | Method and device for processing a multi-channel measurement of magnetic fields | |
US8838225B2 (en) | Analysis of multi-channel measurement data using orthogonal virtual channels | |
US7933727B2 (en) | Method and device for interference suppression in electromagnetic multi-channel measurement | |
JP2006520892A5 (en) | ||
US4977896A (en) | Analysis of biological signals using data from arrays of sensors | |
JP5894156B2 (en) | Method for adjusting interference signal space in biomagnetic field measurement | |
US12038488B2 (en) | Measuring apparatus, measuring method and recording medium | |
JP2018004286A (en) | Signal processing device, signal processing method, signal processing program, and magnetic field measurement system | |
JPH04215084A (en) | Method for measuring field pattern of electric field or magnetic field | |
Nolte et al. | The effect of artifact rejection by signal-space projection on source localization accuracy in MEG measurements | |
EP0477434B2 (en) | Analysis of biological signals using data from arrays of sensors | |
US7502720B2 (en) | Method and device for using a multi-channel measurement signal in source modelling | |
CN114002634A (en) | Calibration device and method for magnetic field-voltage coefficient of multi-channel magnetocardiogram detection system | |
Değirmenci et al. | Anisotropic conductivity imaging with MREIT using equipotential projection algorithm | |
Nurminen et al. | Effects of sensor calibration, balancing and parametrization on the signal space separation method | |
US11668772B2 (en) | Magnetic field measurement device, magnetic field measurement method, and recording medium having recorded thereon magnetic field measurement program | |
Di Rienzo et al. | Theoretical lower error bound for comparative evaluation of sensor arrays in magnetostatic linear inverse problems | |
JP2022153083A (en) | Measurement device, method for measurement, and program | |
JPH05220124A (en) | Biomagnetism measuring instrument | |
JP2022153075A (en) | Measurement device, method for measurement, and program | |
JP2022153186A (en) | Measurement device, method for measurement, and program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): BW GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 11220415 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2518674 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2006505618 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2004719479 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 2004719479 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 11220415 Country of ref document: US |