DEVICES AND METHODS FOR DETERMINING A MAGNETIC FIELD

Abstract

A method and device for determining values of a magnetic field component of a magnetic vector field.

A method for determining values of a magnetic field component of a magnetic vector field, comprising: determining first distribution data comprising values of the magnetic field component, for a first predetermined area defined along a predetermined surface; determining second distribution data comprising second values of the component of the magnetic field for a second predetermined area defined along a second predetermined surface, wherein the first and the second predetermined surfaces are parallel;
wherein determining second distribution data comprises manipulation of the first distribution data based on making use of intrinsic physical properties of the magnetic field; and associated device.

Claims

1. A method for determining values of a magnetic field component of a magnetic vector field, comprising: determining first distribution data comprising values of said magnetic field component, for a first predetermined area defined along a first predetermined surface; determining second distribution data comprising second values of said component of said magnetic field for a second predetermined area defined along a second predetermined surface, wherein said first and said second predetermined surfaces are parallel; wherein determining second distribution data comprises manipulation of said first distribution data based on making use of intrinsic physical properties of said magnetic field.

2. The method according to claim 1, wherein said magnetic field component comprises the magnitude of the projection of said magnetic field vector on an axis or on a surface, or the magnitude of the magnetic field vector.

3. The method according to claim 1, wherein determining second distribution data of said component comprises performing a Fourier transformation of said first distribution data resulting in Fourier transformed data, followed by performing data manipulation on said Fourier transformed data resulting in manipulated Fourier transformed data, followed by an inverse Fourier transformation of said manipulated Fourier transformed data.

4. The method according to claim 3, wherein manipulation of said Fourier transformed data comprises multiplying said Fourier transformed data with a factor which is a function of spatial frequencies corresponding to a first and a second direction, said first and second direction being orthogonal and defining said first predetermined surface.

5. The method according to claim 4, wherein said factor is a function of a magnitude of a spatial frequency vector determined by said first and second direction.

6. The method according to claim 4, wherein said factor comprises an exponential function.

7. The method according to claim 4, whereby said factor is a function of a distance along a third direction, said third direction being orthogonal on said first and said second direction, between said first predetermined surface and said second predetermined surface.

8. The method according to claim 6, wherein said exponential function comprises said distance along said third direction in its exponent.

9. The method according to claim 6, wherein said exponential function comprises a magnitude of the spatial frequency vector determined by said first and second directions in its exponent.

10. The method according to claim 1, wherein determining distribution data comprising values of said component of a magnetic field, for a first predetermined area defined along a predetermined surface, comprises measuring measurement values of said component by means of a magnetic field camera.

11. The method according to claim 10, wherein determining distribution data comprising values of a component of a magnetic field, for a first predetermined area defined along said first predetermined surface, further comprises modeling said measurement values of said component based on a predetermined model and/or predetermined input parameters.

12. The method according to claim 1, wherein said distribution data of said component comprises non-zero values for said component of said magnetic field corresponding to a location at an outer border of said first predetermined area.

13. The method according to claim 1, further comprising: generating additional distribution data of said component, said additional distribution data comprising expected values for said component of said magnetic field in an extension area, said extension area adjacent to said predetermined area and along said first surface; and determining second distribution data for an extended set of distribution data, said extended set of distribution data comprising said first distribution data and said additional distribution data.

14. The method according to claim 13, wherein said values of said component comprised in said additional distribution data are set to be monotonously decreasing to zero in said extension area when moving from an outer boundary of said predetermined area, away from said predetermined area, towards an outer border of said extension area.

15. A device for determining values of a magnetic field component of a magnetic vector field, comprising: a means for determining first distribution data comprising values of said magnetic field component, for a first predetermined area defined along a predetermined surface; a means for determining second distribution data comprising second values of said component of said magnetic field for a second predetermined area defined along a second predetermined surface, wherein said first and said second predetermined surfaces are parallel; wherein said means for determining second distribution data is adapted for manipulating said first distribution data based on making use of intrinsic physical properties of said magnetic field.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0082] The disclosure will be further elucidated by means of the following description and the appended figures.

[0083] FIGS. 1 to 8 illustrate embodiments of the present invention. FIGS. 3 to 8 illustrate examples of how additional distribution data for the component of the magnetic field vector can be generated, according to embodiments of the present invention.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0084] The present disclosure will be described with respect to particular embodiments and with reference to certain drawings but the disclosure is not limited thereto but only by the claims. The drawings described are only schematic and are non-limiting. In the drawings, the size of some of the elements may be exaggerated and not drawn on scale for illustrative purposes. The dimensions and the relative dimensions do not necessarily correspond to actual reductions to practice of the disclosure.

[0085] Furthermore, the terms first, second, third and the like in the description and in the claims, are used for distinguishing between similar elements and not necessarily for describing a sequential or chronological order. The terms are interchangeable under appropriate circumstances and the embodiments of the disclosure can operate in other sequences than described or illustrated herein.

[0086] Moreover, the terms top, bottom, over, under and the like in the description and the claims are used for descriptive purposes and not necessarily for describing relative positions. The terms so used are interchangeable under appropriate circumstances and the embodiments of the disclosure described herein can operate in other orientations than described or illustrated herein.

[0087] Furthermore, the various embodiments, although referred to as preferred are to be construed as exemplary manners in which the disclosure may be implemented rather than as limiting the scope of the disclosure.

[0088] A method is described for determining values of a magnetic field component of a magnetic vector field, comprising: [0089] determining first distribution data comprising values of the magnetic field component, for a first predetermined area defined along a first predetermined surface; [0090] determining second distribution data comprising second values of the component of the magnetic field for a second predetermined area defined along a second predetermined surface, wherein the first and the second predetermined surfaces are parallel;
wherein determining second distribution data comprises manipulation of the distribution data based on making use of intrinsic physical properties of the magnetic field.

[0091] According to preferred embodiments of the present invention, determining second distribution data of the component comprises performing a Fourier transformation of the first distribution data resulting in Fourier transformed data, followed by performing data manipulation on the Fourier transformed data resulting in manipulated Fourier transformed data, followed by an inverse Fourier transformation of the manipulated Fourier transformed data.

[0092] Below, a derivation is presented demonstrating this aspect.

[0093] It is known to persons skilled in the art of physics that a magnetic vector field {right arrow over (B)}=(B.sub.x,B.sub.y,B.sub.z) can be expressed as

{right arrow over (B)}=(B.sub.x,B.sub.y,B.sub.z)=(.sub.x.sup.,.sub.x.sup.,.sub.x.sup.)=,Equation 1

where is the magnetic potential, and is the del operator, also called nabla operator known to the skilled person.

[0094] One set of intrinsic physical properties of a harmonic potential field are the Green's Identities, which are known by persons skilled in the art of mathematics. Consider U to be a closed region in three-dimensional space, and U the boundary surface of this region. Consider and to be harmonic continuous functions with continuous partial derivatives of first and second orders in the region U. From Green's Second Identity it then follows that

[00001] $\begin{matrix} \frac{1}{4 .Math.} .Math._{U} .Math. (.Math. \frac{}{n} - .Math. \frac{}{n}) .Math. S = 0, & Equation .Math. .Math. 2 \end{matrix}$

where

[00002] $\frac{}{n}$

is the directional derivative of in the direction of the outward pointing normal {right arrow over (n)} to the surface element dS.

[0095] If P is a point inside the closed region U, then it follows from Green's Third Identity that

[00003] $\begin{matrix} (P) = \frac{1}{4 .Math.} .Math._{U} .Math. (\frac{1}{r} .Math. \frac{}{n} - .Math. \frac{}{n} .Math. \frac{1}{r}) .Math. S, & Equation .Math. .Math. 3 \end{matrix}$

where r is the distance between the surface portion dS and the point P.

[0096] It is known to persons skilled in the art of physics, that the magnetic potential of Equation 1 is a harmonious function with continuous partial derivatives of first and second orders in a region U in three-dimensional space where no sources of magnetic field are present.

[0097] Consider a Cartesian coordinate system XYZ, in which one measures the magnetic field on a XY plane at a certain position z.sub.0 along the Z axis. It is supposed that all sources of magnetic field are located at z<z.sub.0, that means, below the measurement surface. Suppose that one wants to know the magnetic field at a point P with coordinates (x,y,z.sub.0+z), that means, at a distance z above the measurement surface. One defines a region R in space enclosed by a disk with radius p lying in the measurement plane and a half sphere with the same centre and radius as the disk, extending from the disk circumference in the positive Z direction. It will be apparent to those skilled in the art that in this closed region Equation 3 is valid, whereby the function is understood to be the magnetic potential. The integral in Equation 3 can thereby be expressed as the sum of a first integral over the disk lying in the XY plane, and a second integral over the hemisphere. Next it is considered that the value of p evolves to infinity, in which case it is easily shown that the second integral, namely the one over the hemisphere, evolves to the value zero, and that the integration surface U reduces to the complete XY plane (x,y,z.sub.0). In the limit .fwdarw., Equation thus becomes:

[00004] $\begin{matrix} (x, y, z_{0} + .Math. .Math. z) = \frac{1}{4 .Math.} .Math._{-}^{+} .Math._{-}^{+} .Math. .Math. (\frac{1}{r} .Math. \frac{(x^{}, y^{}, z_{0})}{z^{}} - (x^{}, y^{}, z_{0}) .Math. \frac{}{z^{}} .Math. \frac{1}{r}) .Math. x^{} .Math. y^{}, & Equation .Math. .Math. 4 \end{matrix}$

where r={square root over ((xx).sup.2+(yy).sup.2+(z.sub.0+zz).sup.2)}, and where z>0.

[0098] Equation 4 contains a term with

[00005] $\frac{(x^{}, y^{}, x_{0})}{z^{}},$

which can be eliminated in the following way. By adding together Equation 2 and Equation 3, one obtains:

[00006] $\begin{matrix} (P) = \frac{1}{4 .Math.} .Math._{U} .Math. [(+ \frac{1}{r}) .Math. \frac{}{n} - .Math. \frac{}{n} .Math. (+ \frac{1}{r})] .Math. S . & Equation .Math. .Math. 5 \end{matrix}$

[0099] One now defines the point P as the mirror image of the point P with respect to the plane z=z.sub.0, that means with P having coordinates (x,y,z.sub.0z), and one defines =1/, where ={square root over ((xx).sup.2+(yy).sup.2+(z.sub.0zz).sup.2)}. With this definition following conditions are satisfied: firstly,

[00007] $+ \frac{1}{r} = 0$

on the the XY surface where z=z.sub.0; secondly,

[00008] $+ \frac{1}{r}$

vanishes on the hemisphere in the limit .fwdarw.; and thirdly, is harmonic. In the limit .fwdarw. Equation 5 then becomes:

[00009] $\begin{matrix} (x, y, z_{0} + .Math. .Math. z) = - \frac{1}{4 .Math.} .Math._{-}^{+} .Math._{-}^{+} .Math. .Math. ((x^{}, y^{}, z_{0}) .Math. \frac{}{z^{}} [\frac{1}{r} - \frac{1}{}]) .Math. x^{} .Math. y^{} . & Equation .Math. .Math. 6 \end{matrix}$

When the derivative in this equation is calculated and z is evolving to the XY plane, one obtains:

[00010] $\begin{matrix} (x, y, z_{0} + .Math. .Math. z) = - \frac{.Math. .Math. z}{2 .Math.} .Math._{-}^{+} .Math._{-}^{+} .Math. .Math. \frac{(x^{}, y^{}, z_{0})}{{[{(x - x^{}]}^{2} + {(y - y^{})}^{2} + .Math. .Math. z^{2}]}^{3 / 2}} .Math. x^{} .Math. y^{}, & Equation .Math. .Math. 7 \end{matrix}$

where z>0. Equation 7 can be used to calculate the potential at a point (x,y,z.sub.0+z) based on measurements in a plane (x,y,z.sub.0).

[0100] In order to be useable in practice, Equation 7 can be considered in the Fourier domain. It is noted that Equation 7 has the form of a two-dimensional convolution:

(x,y,z.sub.0+z).sub..sup.+.sub..sup.+(x,y,z.sub.0)u(xx,yy,z)dxdy,Equation 8

where

[00011] $\begin{matrix} u (x, y, .Math. .Math. z) = \frac{.Math. .Math. z}{2 .Math.} .Math. \frac{1}{{(x^{2} + y^{2} + .Math. .Math. z^{2})}^{3 / 2}} . & Equation .Math. .Math. 9 \end{matrix}$

[0101] A convolution in the spatial domain corresponds to a multiplication in the Fourier domain. Therefore Equation 8 can be written as

{circumflex over ()}(k.sub.z,k.sub.y,z+z)={circumflex over ()}(k.sub.x,k.sub.y,z.sub.0){circumflex over (u)}(k.sub.x,k.sub.y,z),Equation 10

Where {umlaut over ()}(k.sub.x,k.sub.y,z+z) is the two-dimensional Fourier transform of the magnetic potential in the XY plane at z=z.sub.0+z. The notation {circumflex over (F)} denotes the Fourier transform of the function F.

[0102] Equation 9 can also be written as

[00012] $\begin{matrix} u (x, y, .Math. .Math. z) = - \frac{1}{2 .Math.} .Math. \frac{}{.Math. .Math. z} .Math. \frac{1}{r}, & Equation .Math. .Math. 11 \end{matrix}$

where r={square root over (x.sup.2+y.sup.2+z.sup.2)}. It can be shown that the Fourier transform of the function 1/r is given by:

[00013] $\begin{matrix} = 2 .Math. .Math. \frac{^{- .Math. k .Math. .Math. .Math. .Math. z}}{.Math. k .Math.}, | k | 0. & Equation .Math. .Math. 12 \end{matrix}$

[0103] The Fourier transform (k.sub.x,k.sub.y,z) of the function u(x,y,z) is then calculated as follows:

[00014] $\begin{matrix} \begin{matrix} \hat{u} (k_{x}, k_{y}, .Math. .Math. z) = .Math. \frac{1}{2 .Math.} .Math. \frac{}{.Math. .Math. z} .Math. \\ = .Math. \frac{}{.Math. .Math. z} .Math. \frac{^{- .Math. k .Math. .Math. .Math. .Math. z}}{.Math. k .Math.} \\ = .Math.^{- .Math. k .Math. .Math. .Math. .Math. z}, .Math. .Math. z > 0, \end{matrix} .Math. .Math. where .Math. .Math. .Math. k .Math. = \sqrt{k_{x}^{2} + k_{y}^{2}} . & Equation .Math. .Math. 13 \end{matrix}$

[0104] Equation 10 was derived for the magnetic potential . It is now shown that it is also valid for each magnetic field component B.sub.x, B.sub.y and B.sub.z separately. From Equation 1, Equation 10 and Equation 13 it follows that

[00015] $\begin{matrix} Equation .Math. .Math. 14 \\ and \\ \begin{matrix} .Math. (k_{x}, k_{y}, z_{0} + .Math. .Math. z) = .Math. - \frac{\overline{} .Math. \overline{}}{z} .Math. (k_{x}, k_{y}, z_{0} + .Math. .Math. z) \\ = .Math. - .Math. k .Math. .Math. \hat{} (k_{x}, k_{y}, z_{0} + .Math. .Math. z) \\ = .Math. - .Math. k .Math. .Math. \hat{} (k_{x}, k_{y}, z_{0}) .Math. \hat{u} (k_{x}, k_{y}, .Math. .Math. z) \\ = .Math. - .Math. (k_{x}, k_{y}, z_{0}) .Math. \hat{u} (k_{x}, k_{y}, .Math. .Math. z) \\ = .Math. .Math. (k_{x}, k_{y}, z_{0}) .Math. \hat{u} (k_{x}, k_{y}, .Math. .Math. z) \\ = .Math. .Math. (k_{x}, k_{y}, z_{0}) .Math.^{- .Math. k .Math. .Math. .Math. .Math. z} \end{matrix} & Equation .Math. .Math. 1 \end{matrix}$

[0105] In Equation 14 and Equation 15 the following properties were used, which are known to persons skilled in the art:

[00016] $\begin{matrix} = {({ik}_{x})}^{n} .Math. \hat{} .Math. .Math. = {({ik}_{y})}^{n} .Math. \hat{} .Math. .Math. = {.Math. k .Math.}^{n} .Math. \hat{} & Equation .Math. .Math. 16 \end{matrix}$

[0106] Equation 14 and Equation 15 show that the result of Equation 10 can be applied to each component of the vector field derived from , as summarized here:

custom-character (k.sub.x,k.sub.yz.sub.0+z)=(k.sub.x,k.sub.y,z.sub.0)e.sup.|k|z.Equation 17

[0107] It follows from Equation 17 that the magnetic field in the XY plane at z=z.sub.0+z can be obtained from the magnetic potential in the XY plane at z=z.sub.0 by first performing a Fourier transform, then multiplying by e.sup.|k|z, and then performing an inverse Fourier transform.

[0108] Although in the derivation of Equation 17 it is assumed that z>0, which means that the field is calculated at distances farther away from the magnetic field source, it can also be applied with z<0, that means, to calculate the magnetic field at distances closer to the magnetic field source. A condition is that the position where the field is calculated is still in a region free from magnetic field sources. In practice this means that it is possible to measure the magnetic field at a certain distance from a magnet, and that the field can be calculated at positions closer to the magnet, even up to the magnet surface. Since existing magnetic field sensors and magnetic field camera devices often have a minimum measurement distance from the magnet, it is an advantage of the present invention to be able to calculate the field at closer distances to the magnet.

[0109] Although Equation 17 was derived for the case of a Cartesian coordinate system, the same principle can be applied in other coordinate systems, such as the cylindrical coordinate system. First, one considers the case in which any component of the magnetic field is recorded on a flat disk or ring surface in a cylindrical coordinate system (R, , Z). Since the predetermined area is in a flat surface, it can be projected onto a plane in the Cartesian coordinate system. Thereby it is noted that the data grid in the cylindrical coordinate system may not transform to a regular grid in the Cartesian coordinate system. However, the transformed grid can be made regular again by interpolation methods which are well known to persons skilled in the art. To the obtained data, Equation 17 can be applied, in order to obtain the magnetic field at another distance in the Z-direction. The coordinates of the obtained data are then transformed back to the cylindrical coordinate system. Thereby, according to preferred embodiments, another interpolation step is performed in order to obtain the data values at the original grid points in the cylindrical coordinate system.

[0110] In order to obtain suitable boundary conditions on the outer periphery of the transformed disk surface in Cartesian coordinates, it is preferred to apply extrapolation methods according to methods described in the present disclosure.

[0111] The person skilled in the art will recognize that the above method equally applies to a ring surface, where only the area between an inner radius and an outer radius is considered. Thereby, the inner disk surface, on which no data points are present, may be extrapolated according to methods described in the present invention.

[0112] The principle of Equation 17 can be used also for data on a curved cylinder surface, in order to calculate the magnetic fields at radial distances different from this surface. This is already apparent from Equation 3, which says that the magnetic field in a closed region can be deduced from the field on the boundary of that region. For the case where the predetermined area is a curved cylinder surface with radius R.sub.0, consider the region bound by two cylinders, one with radius R.sub.0 and one with radius R.sub.1 with R.sub.0<R.sub.1, and both extending from z.sub.0 to +z.sub.0 in the Z-direction. In analogy to the given derivation for the Cartesian coordinate system, it is apparent that the integrand of Equation vanishes on the outer cylinder and on the planar ring surfaces in the limit where R.sub.1,z.sub.0.fwdarw., which means that the field at any point in the region R>R.sub.0 can be derived from the field at R.sub.0 using an expression similar to Equation 17, expressed in cylindrical coordinates.

[0113] Furthermore, it will be apparent to those skilled in the art that this method is equally applicable to determine the field at a radius R<R.sub.0, assuming that no magnetic field sources are present in the region between R and R.sub.0.

[0114] When the distribution of the measured magnetic field component on the boundaries of the predetermined area is not approaching zero, because for example the magnetic object is larger than the measurement area of the magnetic field camera, it can be needed to extend the predetermined measurement area in order to make sure that the measured distribution is sufficiently approaching zero at its boundaries, so that the methods of the present invention can be applied. For such case, it may be needed to measure a larger area by stepping the measurement area in the first and/or second directions, and at each new location measure the distribution in a new predetermined area which is adjacent to the previous predetermined area, after which all determined distributions are stitched together as to result in one large predetermined area, on the borders of which the magnetic field has sufficiently evolved towards zero as to apply the extrapolation methods described in the present disclosure, with the aim of determining second distribution data comprising second values of the component of the magnetic field for a second predetermined area defined along a second predetermined surface, wherein the first and the second predetermined surfaces are parallel. Methods and devices for accomplishing this are for example described in European patent application EP12188521.4 filed on 15 Oct. 2012 by the applicant of the present application, which are hereby incorporated by reference. Such a device can be described as a device for determining a magnetic field distribution of a magnet along a main surface of the magnet, the device comprising: [0115] a. an arrangement of at least two independent magnetic field camera modules being arranged in a fixed relative position with respect to each other, each magnetic field camera module being adapted for measuring a magnetic field distribution to which it is exposed by means of a respective detection surface; [0116] b. a means for providing a predetermined relative movement between the main surface and the arrangement to thereby scan the magnetic field distribution of the magnet along the main surface.

[0117] The associated method is a method for determining a magnetic field distribution of a magnet along a main surface of the magnet, the device comprising: [0118] providing the magnet; [0119] providing an arrangement of at least two independent magnetic field camera modules being arranged in a fixed relative position with respect to each other, each magnetic field camera module being adapted for measuring a magnetic field distribution to which it is exposed by means of a respective detection surface; [0120] providing a predetermined relative movement between the main surface and the arrangement to thereby scan the magnetic field distribution of the magnet along the main surface.

[0121] In order to determine the second distribution data comprising second values of the magnetic field component for this stitched distribution, the methods described in the present invention can be applied to the resulting (stitched) predetermined area.

[0122] Another method however is to apply the so-called overlap-add and overlap-save methods, which are known to persons skilled in the art of signal processing, whereby the Fourier transform, data manipulation based on intrinsic physical properties of a magnetic vector field, and inverse transform are performed on sub areas of the respective predetermined areas or combination of predetermined area and extension area.

[0123] The overlap-add and overlap-save methods can also be applied in the other cases described in the present disclosure in the following way. Instead of first determining the additional distribution data in the extended area, and after that applying data manipulation based on intrinsic physical properties of a magnetic vector field (for instance comprising applying a Fourier transform, manipulating the Fourier transformed data, and performing an inverse Fourier Transform) in order to obtain the second distribution data of the component, the Fourier Transform, data manipulations, and Inverse Fourier Transform are applied separately to the predetermined area and to (sub-regions of) the extended area in a block-wise fashion. Thereby each of the blocks can be zero padded to a certain extent in order to create overlap regions between the back-transformed blocks. The resulting back-transformed blocks are then combined in the final larger matrix, whereby the overlap regions are added together. Additionally, the predetermined area can itself also be treated block-wise, for example in the case of a large magnet where the predetermined area is relatively large and is measured in a block-wise way. This method corresponds to the overlap-add method. A person skilled in the art will readily recognize the analogous possibility of applying the overlap-save method to the same data.

[0124] Alternatively, the predetermined and extended areas can be treated block-wise or section-wise, where the manipulations are performed on each block/section separately, and the results of them added together or saved according to overlap-add and overlap-save methods, respectively.

[0125] As a further improvement on the previous method, the overlap-add and overlap-save methods can also be used on an infinite extension area, where the extrapolation is expressed as an analytical function which the Fourier transform and other operations can be analytically determined on the interval stretching from infinity to the border of the predetermined area.

[0126] Distribution data of the component of the magnetic vector field in the extension area can be determined in many ways, some of which are explained below.

[0127] FIGS. 3 to 8 illustrate examples of how additional distribution data for the component of the magnetic vector field can be generated, according to embodiments of the present disclosure. The depicted patterns illustrate the evolution of the values of the component, the z-component of the magnetic field, along to a cross-section S as depicted in FIG. 1. The cross-section can comprise an outer portion of the predetermined area only, but can extend up until the centre of the predetermined area.

[0128] According to a first embodiment, illustrated in FIG. 3, the component of the magnetic field vector is set to value zero (0) in the complete extension area. When the measured values of the component of the magnetic vector field on the outer boundary of the predetermined area differ from zero, the method may though introduce a discontinuity in the values on the outer boundary of the predetermined area, which may still introduce errors in a subsequent calculation of the second distribution data of the component of the magnetic field vector.

[0129] According to an alternative embodiment, illustrated in FIG. 8, of the present disclosure, the additional distribution data comprising the component of the magnetic field vector is simulated or modeled for the predetermined area and the simulation results or model are/is extrapolated into the extension area.

[0130] According to preferred embodiments, illustrated in FIG. 2, the magnetic field component at locations lying in the extension area may be set to the same value as the nearest point in the predetermined area. For a rectangular predetermined area, this would mean that the boundary values of the predetermined area and thus measured area are set as a fixed, constant, value throughout the extension area along a direction which is orthogonal on the outer boundary of the predetermined area (portions 20). In the corner areas which are then remaining (21), the value of the corresponding corner of the predetermined area is set as a constant value. Subsequently, a window function is applied on the predetermined area and the extension area, which has a value 1 in the predetermined or measured area and which evolves continuously to a (near to) zero value on the outer boundary of the predetermined area. For instance, the window function can be a Tukey-window (illustrated in FIG. 5) or a Planck-Taper window. Alternatively a so-called bump-function or a test-function can be used, which have the property that they evolve from a value 1 to a value zero within a limited area, whereby they can be infinitely differentiated, and whereby they thus do not introduce discontinuities in the function or any of its first or higher order derivatives, which makes them suitable for being used as a window function in this context.

[0131] According to a further embodiment, illustrated in FIG. 4, the values of the component of the magnetic vector field at the outer boundary of the predetermined area are set to be evolving exponentially towards a zero value within the extension area. Preferably the exponential factor is predetermined such that its sufficiently small value or zero value is achieved on the outer boundary of the extension area. In practice, it is enough that the magnetic field distribution on the outer boundary of the extension area has a value which is below the measurement noise. This method guarantees the continuity of the values on the outer border of the predetermined area, but does not guarantee the continuity of the first derivative thereof.

[0132] In a further preferred embodiment, the component of the magnetic vector field in the (boundary area of the) predetermined area is represented by a polynomial representation. The order of the extrapolated polynomial function can be reduced, for instance to order two (quadratic polynomial), or to order one (linear polynomial), in order not to obtain instable extrapolation values. This method guarantees that the continuity of the measurement values and of the derivatives is ascertained to the same extent as the order of the used polynomial function. In practice, it is further preferred to further apply a window-function on these extrapolated values, such that it can be guaranteed that the distribution of the component of the magnetic field vector reaches a value which is small enough or zero at the outer boundary of the extension area.

[0133] In a further preferred embodiment, the component of the magnetic vector field in the (boundary area of the) predetermined area is represented by a rational function representation. The order of the extrapolated rational function can be chosen, for instance to order 1 (1/x), order 2 (1/x.sup.2) or order 3 (1/x.sup.3). The latter case is particularly useful since it corresponds to the decline rate of the magnetic field of a magnetic dipole, which is an approximation of a magnet at large distances from the magnet field source, e.g. magnet. This method guarantees that the magnetic field approaches zero when moving further into the extension area. In practice, it can further be preferred to apply a window-function on these extrapolated values, such that it can be guaranteed that the distribution of the component of the magnetic field vector reaches a value which is small enough or zero at the outer boundary of the extension area.

[0134] In a further preferred embodiment, the component of the magnetic vector field in the (boundary area of the) predetermined area is represented by a spline-representation, illustrated in FIG. 6. The order of the extrapolated spline-function can be reduced, for instance to order two (quadratic spline), or to order one (linear spline), in order not to obtain instable extrapolation values. This method guarantees that the continuity of the measurement values and of the derivatives is ascertained to the same extent as the order of the used spline-function. In practice, it is further preferred to further apply a window-function on these extrapolated values, such that it can be guaranteed that the distribution of the component of the magnetic field vector reaches a value which is small enough or zero at the outer boundary of the extension area. This is illustrated in FIG. 7.

[0135] According to a still further embodiment of the present disclosure, the values of the component of the magnetic field at the boundary of the predetermined area can be extrapolated on the basis of a simulated magnetic field distribution for the magnetic object of which the magnetic field distribution was measured. Simulation algorithms can be used to calculate the magnetic field distribution based for instance on one or more properties of the magnetic object such as for instance shape, material, magnetization vector, position and so forth. Also here, use can be made of any of the methods disclosed for instance in EP2508906.

[0136] According to another embodiment of the present disclosure, the determination of the second distribution data of the component of the magnetic field distribution of a stitched predetermined area is performed in the following way, by using the overlap-add or overlap-save method. Thereby the Fourier transform, manipulation operations and inverse Fourier transform are performed on each individual predetermined area, which may or may not be zero padded, without applying the extrapolation methods described in the current disclosure. When combining the determined second data distributions of the component into the larger distribution, the overlap-save or overlap-add method can be applied.

[0137] According to another embodiment of the present disclosure, the overlap-add or overlap-save method is used on a finite or an infinite extension area.

DEVICES AND METHODS FOR DETERMINING A MAGNETIC FIELD

Assignee

Inventors

Cpc classification

Classification Explorer

G01R33/10

PHYSICS

Classification Explorer

G01R33/12

PHYSICS

Classification Explorer

G01R33/1207

PHYSICS

Classification Explorer

G01R33/072

PHYSICS

Classification Explorer

G01R33/1253

PHYSICS

Classification Explorer

G01R33/1215

PHYSICS

International classification

Classification Explorer

G01R33/10

PHYSICS

Classification Explorer

G01R33/12

PHYSICS

Abstract

Claims

Description