EP2191444A1 - Verfahren zur tomographischen rekonstruktion einer messgrösse durch gezielte bestimmung der messpunkte - Google Patents
Verfahren zur tomographischen rekonstruktion einer messgrösse durch gezielte bestimmung der messpunkteInfo
- Publication number
- EP2191444A1 EP2191444A1 EP08801207A EP08801207A EP2191444A1 EP 2191444 A1 EP2191444 A1 EP 2191444A1 EP 08801207 A EP08801207 A EP 08801207A EP 08801207 A EP08801207 A EP 08801207A EP 2191444 A1 EP2191444 A1 EP 2191444A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- distribution
- measured
- variable
- auxiliary
- values
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000009826 distribution Methods 0.000 claims abstract description 158
- 230000008859 change Effects 0.000 claims abstract description 24
- 239000000446 fuel Substances 0.000 claims description 26
- 238000012937 correction Methods 0.000 claims description 16
- 230000007704 transition Effects 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 claims description 4
- 239000000523 sample Substances 0.000 description 53
- 238000005457 optimization Methods 0.000 description 52
- 239000011159 matrix material Substances 0.000 description 49
- 239000013598 vector Substances 0.000 description 18
- 239000010410 layer Substances 0.000 description 13
- 230000005484 gravity Effects 0.000 description 9
- 230000036961 partial effect Effects 0.000 description 9
- 230000002829 reductive effect Effects 0.000 description 9
- 238000013459 approach Methods 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 239000000047 product Substances 0.000 description 6
- 238000010606 normalization Methods 0.000 description 5
- 230000005624 perturbation theories Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 3
- 230000003197 catalytic effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000010439 graphite Substances 0.000 description 3
- 229910002804 graphite Inorganic materials 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 229940050561 matrix product Drugs 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000009827 uniform distribution Methods 0.000 description 3
- 108010014172 Factor V Proteins 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000010936 titanium Substances 0.000 description 2
- 229910052719 titanium Inorganic materials 0.000 description 2
- 229920000557 Nafion® Polymers 0.000 description 1
- 239000004809 Teflon Substances 0.000 description 1
- 229920006362 Teflon® Polymers 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000002346 layers by function Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- CFQCIHVMOFOCGH-UHFFFAOYSA-N platinum ruthenium Chemical compound [Ru].[Pt] CFQCIHVMOFOCGH-UHFFFAOYSA-N 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/08—Measuring current density
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
Definitions
- the invention relates to a method for the tomographic reconstruction of a measured variable from an auxiliary variable accessible to the measurement, a method for determining optimum measuring locations for the auxiliary variable and apparatus for carrying out the method.
- auxiliary variable is measured, which is physically related to the distribution of the measured variable and which is more accessible to a measurement. From the values of the auxiliary quantity, the distribution of the measured quantity in the interior of the sample is then concluded by tomographic reconstruction. For example, in computer tomography, the density distribution of matter in the
- Body, the measure, and the absorption of guided by the body X-ray radiation is used as an auxiliary size.
- Tomographic reconstructions require measurements of auxiliary size in many locations. This requires either a large number of detectors to simultaneously measure values of the auxiliary size in many locations, or a lot of measurement time to successively measure the values at many locations with one or more movable detectors. Task and solution
- a method for the tomographic reconstruction of the distribution of a measured variable within a sample from values of an auxiliary variable was found. 5
- the values of the auxiliary size can be measured outside the sample. However, they can also be expressly measured within the sample, for example, by introducing a corresponding probe into the sample.
- the values of the auxiliary variable are weighted on the basis of a relevance measure.
- a measure of relevance means any function which assigns a numerical value to a location, which is a measure of the significance of a value of the auxiliary variable measured at this location with regard to the distribution of the measured variable, and which at the same time does not alone constitute a measure of is the distance of the location from a reference point of the sample.
- this reference point may be the center of the sample.
- the specific design of the relevance measure depends on what type of information is requested in the specific application.
- a location within the meaning of this invention is in a frame of reference in which the sample rests. Thus, it is irrelevant, for example, for a measurement of the auxiliary size at two different locations, whether the detector used for this measurement is moving (or physically present a second time at the location of the second measurement) or whether the sample is being moved, in particular translated or rotated.
- a relevance measure does not have to be explained in the entire room. For example, it may only be explained on a set of locations that satisfy a constraint imposed by the use case.
- the relevance measure can be explained, for example, only on the quantity of those places which are technically or financially accessible for a measurement of the auxiliary quantity, so that it is expedient to select relevant locations only from this quantity.
- a measure of relevance does not necessarily have to be describable by an analytical formula, nor should it be continuous.
- the values of the relevance measure can be determined exclusively by numerical simulation on the basis of a model of the sample. Such numerical simulations can also be used by the person skilled in the art in order to determine a relevance measure for a specific application. If, for example, the sample is to be examined for the presence of certain defects, it can introduce these defects into an undisturbed model of the sample, study its effect on the values of the auxiliary size and configure the relevance measurement to give particular value to the locations at which measurements are taken are particularly important for reliable detection of defects.
- measured values of the auxiliary variable from locations at which they can only be measured inaccurately or only with difficulty can be weighted with the value zero and thus completely ignored during the tomographic reconstruction.
- the cost of data acquisition can be reduced. If it is known in advance by the relevance measure at which locations measurements with regard to the distribution of the measured variable are particularly worthwhile, measurements can advantageously be limited to these locations. If the values of the auxiliary variable are measured with several stationary detectors, fewer detectors are required for the same quality of the reconstruction. If the values of the auxiliary variable are measured with one or more movable detectors which approach the measuring locations in each case, a shorter measuring time is required for the same quality of the reconstruction. This allows dynamic processes in the sample to be studied on a shorter time scale.
- the optimal distribution of measurement locations can be used, which was determined by an optimization method, or this optimal distribution can at least enter into such a relevance measure.
- a distribution can be described as density, for example.
- the invention also relates to a method for determining the optimal distribution of measurement sites at which for the purpose of tomographic reconstruction a measurand within the sample is to measure an auxiliary variable accessible to the measurement.
- the distribution of measurement locations is optimized to the extreme value of a function, which includes a convolution of this distribution of locations with a function f (S) of the operator S, which assigns a distribution of the auxiliary variable to each possible distribution of the measured variable.
- the invention also relates to a method for the tomographic reconstruction of the distribution of a measured variable within a sample from values of an auxiliary variable, in which the thus optimized distribution of measuring locations is used as relevance measure.
- f (S) can also contain an explicit location dependence f (S, x). may include S + S with f (S), in particular the multiplication of the adjoint operator to S. In general, it is advantageous if f (S) contains the multiplication of an operator with the adjoining operator.
- the operator S represents the tomographic problem to be solved. In addition to the physical relationship between measured variable and auxiliary variable, it also contains the dimensionality of the spaces in which the distribution of the auxiliary variable is measured or the distribution of the measured variable is sought.
- the function f (S) can also be derived from an operator derived from the S operator
- Hang operator For example, if S is a linear operator, it can be written as a matrix that can be expressed as UWV ⁇ by means of singular value decomposition .
- the matrix may be discrete or continuously indexed.
- U is an orthogonal matrix
- W a diagonal matrix with the singular values
- V a unitary matrix.
- a function f (U) of the operator U is therefore also a function f (S): Since U can be calculated from S, f (U (S)) can also be computed into a form f (S).
- the distribution of measurement sites may, for example, be embodied in a function which, in its domain of definition, assigns to each location a local density of measurement locations per unit volume, area or length along a curve.
- the notion of distribution of measurement locations is not limited to such functions, and particularly not to those that vary continuously with location. It comprises any assignment which indicates to a location in its domain of definition whether or even with which priority the auxiliary variable is to be measured at exactly this location in order to achieve an optimal tomographic reconstruction of the measured variable.
- the convolution of the distribution of measurement locations with the function f (S) is understood to mean any integral and sum over more than one location x, which contains an association of the distribution at location x with this function f (S). In this case, the location dependence of f (S), which is always implicitly given over S, must be taken into account. This notion of convolution also explicitly extends to functions f (S), which are matrices with discrete or continuous indices.
- the defined target may also have some influence on the choice of the operator who enters the function f (S). If the operator S is used in the above linear example, then the convolution contains a measure of how strongly changes occur. within the sample to measured values of the auxiliary variable. If the operator U is used, the convolution contains a measure of how much information a measured value of the auxiliary variable contains over the entire distribution of the measured variable.
- the functional is generally given the distribution of measurement sites as an argument and provides a scalar value W.
- this value can be a measure of a desired advantageous property or a combination of such properties, such as for the achievable in the tomographic reconstruction accuracy ,
- the invention likewise relates to a method for determining the optimal distribution of measurement locations at which an auxiliary variable accessible to the measurement is to be measured for tomographic reconstruction of a measured variable within the sample, the distribution of measuring locations being optimized to the extreme value of a function contains at least one singular value of an operator which is a product of a power of the distribution of
- the power does not have to be an integer, but can also be rational or real.
- the product may include the square root of the distribution of measurement locations.
- Distribution of locations with a function f (S) of the operator S which assigns a distribution of the auxiliary size to each possible distribution of the measured variable, is derived in the case of a linear operator in the special description part.
- Sites of measurement correlate directly with how relevant a measurement at each location x where this distribution is explained is for the quality of tomographic reconstruction defined by the functional.
- the functional is usually determined by the optimization goal specified in the specific task. Right.
- the distribution is not only to be interpreted as a prescription for the density of measuring points to be realized in space, but also as a measure of relevance for places x with respect to that function, in the optimization of which it was obtained.
- the optimization method is carried out iteratively, wherein the distribution of measurement locations in the transition from the iteration n-1 to the iteration n by a pointwise product of the iteration n-1 of this distribution with a orthogonal to this distribution, in particular orthonormal , Correction function is changed.
- the objective pursued by the optimization method is defined according to the invention by a functional which contains the distribution of measurement locations. Since, in the transition from the iteration n-1 to the iteration n, the iteration n-1 of the distribution is assumed to be known, the insertion of a convolution of the iteration n-1 with the correction function into the functional determines the determination of the iteration n on the determination attributed to the correction function.
- the functional used for the optimization can vary from one iteration to the next. This allows a flexible trade-off between conflicting optimization goals, such as absolute and relative accuracy.
- a homogeneous distribution in space or even a determined on the basis of a priori information about the distribution of the measured variable distribution of measurement sites into consideration for example, a homogeneous distribution in space or even a determined on the basis of a priori information about the distribution of the measured variable distribution of measurement sites into consideration.
- the relationship between the distribution of the measured variable within the sample and the values of the auxiliary variable is advantageously linearized by an operating point in the space of the possible distributions of the measured variable. Then the distribution of the measurand and the values
- the auxiliary size is then linked by a linear operator S, even if the physical relationship between these quantities is non-linear. For such operators, the singular value decomposition of S provides a statement about how strongly a change in the distribution of the measured variable affects the values of the auxiliary variable.
- the auxiliary size is usually only limited
- the functional advantageously contains a measure for the smallest eigenvalue of the operator S + S, and thus a measure for the smallest singular value of the operator S, which distributes the measured variable represents a distribution of the auxiliary size, or for the change of this value.
- ⁇ 0 is the largest eigenvalue of S + S, to which eigenvectors in the space of the auxiliary size distributions are taken into account.
- the relative error thus depends on the ratio between the largest and the smallest eigenvalue.
- the functional advantageously contains a measure of the ratio between the largest and the smallest singular value of the operator, which maps a distribution of the measured variable to a distribution of the auxiliary size, or for the change of this relationship.
- the distribution of the auxiliary variable is made up of the distribution of measuring locations sought in the optimization method and the measured values which the auxiliary variable assumes at these measuring locations.
- the least possible absolute error and the lowest possible relative error can be conflicting optimization goals for the determination of the distribution of the measurement sites.
- the linear approximation also provides information on the extent to which tomographic reconstruction is possible at all. If the operator S + S (where S + is the operator adjoined to S) has a null space, then there are distributions of the measurand that do not derive any distribution of the values of the values
- Auxiliary size can be determined. The same applies if, in order to improve the absolute or the relative accuracy of the reconstruction, eigenvectors for certain eigenvalues of SS + are not taken into account.
- the spatial resolution with which the distribution of the measured variable can be reconstructed depends on the number N of eigenvectors considered by SS + :
- d w is the dimension of the area in which the distribution of the measurand is sought (line: 1, area: 2, volume: 3).
- is accordingly the length, the area or the volume of this area.
- an a priori expected distribution of the measured variable is selected as the operating point. If, for example, the sample is a technical device, such as a fuel cell, the operating point may be the distribution of the measured variable during trouble-free operation.
- the extreme value of the function is determined under the secondary condition that a convolution of the distribution of the measuring locations with a location-dependent weighting function assumes a predetermined value.
- the extreme value of the function is determined under the secondary condition that a convolution of the distribution of the measurement locations with itself assumes a predetermined value.
- a convolution of the distribution of the measurement sites with themselves is understood in particular to mean an integer or real valued power of this distribution.
- the predefined value can be used to control whether a clustering of Measuring locations or a uniform distribution of sites is preferred.
- a concentration of measuring locations can be advantageous if the measuring locations are difficult to access or it is expensive to move mobile measuring devices over longer distances.
- Auxiliary size and sample are expressly suitable both for the method for tomographic reconstruction and for the optimization method. Also, a device that can perform one of the methods, so that at the same time perform the other.
- a relevance measure is selected whose value is a predetermined maximum of at most three quarters, preferably half and more preferably at most one quarter of all points within any finite distance from the center of the sample on which it is explained Threshold exceeds.
- Threshold exceeds.
- the value of the relevance measure is set in relation to its mean value on a set of points, is to be construed as multiplying all values of the relevance measure by any other than 0 Constants and / or addition of any constant to all values of the relevance measure needs to be fulfilled.
- the effect according to the invention does not depend on absolute values of the relevance measure, but on the relation of those values which the relevance measure assigns to different possible measurement locations for the auxiliary variable. This relation is neither obtained by multiplying all values by a constant other than 0 nor by
- the relevance measure exceeds or falls short of a value, this refers to the amount of the relevance measure. Since a relevance measure is not negative in common usage, this is not indicated in the following for reasons of readability.
- the relevance measure can be a measure of how strongly the value of the auxiliary variable measured at the evaluated location reacts to changes in the measured variable.
- the relevance measure can also be a measure of the size of the resolvable structures in the reconstructed distribution of the measured variable, which are changed when the value of the auxiliary variable measured at the evaluated location changes.
- the total reconstruction m of the distribution of the measured variable can be written as a linear combination of individual contributions, each of which is based solely on a measured value H k of the auxiliary variable.
- the measured values H k of the auxiliary variable to which these contributions are based must be taken into account when determining the distribution m of the measured variable.
- the measured values H k which are particularly sensitive to changes in the distribution m of the measured variable, are not the same, which contribute individual contributions with particularly fine structures to the reconstructed distribution m.
- the relevance measure can advantageously also be a measure of how much information about the entire distribution of the measured variable is contained in the value of the auxiliary variable measured at the evaluated location. Such a measure can be particularly meaningful, for example, when using only a few measured values Hk of the auxiliary variable to obtain a quick overview of the measured variable m in the entire sample volume investigated.
- the matrix S in the above-mentioned relevance measure R k can be written as a singular value decomposition UWV ⁇ , where U is an orthogonal matrix, W is a diagonal matrix with the singular values and V is a unitary matrix. Then the function is
- no values of the auxiliary variable are taken into account which were measured at locations for which the relevance measure is a predetermined one
- Threshold falls below. This selection can take place before the actual data acquisition and thus lead to a saving of required detectors and / or measurement time. Surprisingly, however, it has been found that even a selection after the data acquisition, that is to say a discarding of already recorded measured values of locations rated as less relevant by the relevance measure, is the case
- Quality of reconstruction can improve. If the relevant information about the distribution of the measured variable already exists, it can be nullified by adding irrelevant information. This is due to the fact that the solution of an overdetermined system of equations is fundamentally a decision between several contradictory solutions.
- the overdetermination is caused, for example, by unavoidable measurement inaccuracies.
- the equation system would formally have more equations than unknown ones.
- the surplus equations would be linearly dependent, so that the system of equations in the final effect would not be overdetermined.
- the predetermined threshold below which a measuring location is disregarded, usually has at least the magnitude of the average value of the relevance measure for all measuring locations. Their exact value can be determined, for example, by parameter optimization.
- the aim of such an optimization may be, for example, a maximum accuracy of the tomographic reconstruction or an adjustable ratio between accuracy and effort for data acquisition.
- the predetermined threshold may be 1.5 times, preferably 2 times, the average value of the relevance measure for a given set of measurement locations. This requirement may, for example, be that, due to experimental or financial constraints, measurements are usefully possible only at a certain amount of measurement sites.
- the locations at which the values of the auxiliary variable are measured lie on a cuboid surface, in particular on a cube surface, around the sample. Then, for example, the measuring locations on each individual surface of the cuboid can be approached by a separate, two-dimensional detector.
- the locations at which the values of the auxiliary variable are measured can also lie on an elliptical surface, in particular on a spherical surface, around the sample. Such a shape has maximum symmetry.
- a single detector movable along a spiral path, for example, on the surface can be used to approach all auxiliary size measurement sites.
- auxiliary size there is no restriction on the places where the auxiliary size is measured.
- an electrical current density is selected as the measured variable. This creates a long-range magnetic field that permeates most materials and is therefore measurable outside the sample, where it is accessible without destructive interference with the sample.
- a magnetic field is selected as an auxiliary variable.
- the electric current density j is also subject to restrictions due to Maxwell's equations.
- the magnetic field H in the outer space of the sample is uniquely determined by the current density j within the sample. Since the current density j is related to the electrical field E via the electrical conductivity ⁇ , for which restrictions also apply according to Maxwell's equations, not every current density distribution j in the sample is permissible.
- the current density distribution j sought at N points in the sample can thus be expressed by a set N c of independent variables with N C ⁇ N, with the aid of the constraint that it must be permissible according to Maxwell's equations.
- N R Due to the relative errors of the measuring instruments used to measure the auxiliary magnitude magnetic field, there is an N R for which W kk ⁇ O for k> N R can be assumed.
- N c is determined by the mesh size at which the sample volume is used to set up the linear
- a fuel cell is selected as a sample.
- An operating fuel cell is not or hardly accessible direct measurements of the current density in its interior, since it must be changed for this, which in turn leads to difficult to predict changes in the current density.
- the current density can be determined destructively and without influence as a measured variable.
- a fuel cell which provides only a very low voltage of about one volt, arise in technically relevant applications high currents of the order of 50 A and more. These currents produce magnetic fields that can be measured very well outside the fuel cell.
- each of these devices comprises in each case at least one measuring device for the auxiliary variable and an evaluation unit for
- the measuring device can be arranged outside the sample. But it can also be arranged within the sample; for example, it may be a probe that has been introduced into the sample. It may be able to measure the auxiliary size at different locations in the room. However, it may also be fixed, and means may be provided, such as means for rotating or translating the sample, which at different times act on the values of the measurand from different areas on the sample to the value of the auxiliary quantity at the physical location of the meter to let. For simplicity, therefore, a location will be considered below in a frame of reference in which the sample rests.
- the first embodiment of the device is characterized by an evaluation unit which is able to apply the relevance measure to location coordinates.
- This embodiment offers full flexibility for carrying out the method: it can be decided in advance at which locations at all measurements be performed. However, it is also possible to record a large collection of values H k for the auxiliary variable, with the relevance dimension being evaluated only after the fact in which the values were recorded.
- the evaluation unit may include a list of location coordinates at which the values of the auxiliary variable are determined when carrying out the method.
- these location coordinates can be selected in advance with the aid of the relevance measurement and, for example, introduced into the evaluation unit by means of a data carrier. Therefore, in addition to the location coordinates, the list may also contain the associated values of the relevance measure.
- the values of the relevance measure are at least 1.5 times, preferably at least 2 times, as large as at any other location at which the relevance measure is explained.
- This selection of location coordinates records the small fraction of the possible measurement points, which contains the most information about the desired distribution of the measured variable within the sample.
- the values of the relevance measure are at least 1.5 times, preferably at least 2 times, as large as its average value over all locations at which the relevance measure is explained. This selection of location coordinates also captures the small fraction of the possible measuring points, which contains the most information about the desired distribution of the measured variable within the sample.
- FIG. 3 Values of the relevance measure R k * on a cuboid around the model
- FIG. 4 logarithms of the singular values of the matrix S, which includes the 5536 locations considered in FIG. 3a (curve A); Logarithms of the singular values of the matrix S, which only comprises the 868 locations considered in FIG. 3b (curve B).
- FIG. 5 Logarithm of the relative precision with which the magnetic field measurement has to be performed (ordinate), above the desired spatial resolution of the current distribution j to be reconstructed in centimeters (abscissa).
- FIG. 6 Dependence of the inverse condition of the matrix S (ordinate) on the proportion of possible measurement locations for the auxiliary variable that comprises this matrix S (abscissa).
- FIG. 7 Development of the largest (partial image a) and the smallest (partial image b) eigenvalues of [S + S] in carrying out the optimization method according to the invention using the example of a fuel cell with the current distribution in the interior as the measured variable and the magnetic field in the external space as the auxiliary variable.
- FIG. 8 Representation of the 7% of the measuring sites originally distributed uniformly on the outer surface around the fuel cell, which with the positioning and weighting shown here (area of the dots) the absolute accuracy of the tomographic reconstruction by a factor of 2.77 and the relative accuracy around the Factor 2.8 compared to the original one
- FIG. 1 outlines a model of an experimental direct methanol fuel cell, on which the method according to the invention was tested simulatively.
- the heart of the cell is the Membrane Electrode Unit (MEA). It consists of a porous layer, a catalytic layer, an electrolytic layer (Nafion), a second catalytic layer and a second porous layer.
- MEA Membrane Electrode Unit
- the catalytic layers contain platinum or platinum-ruthenium and are each embedded in porous graphite.
- the MEA consists of 5 layers, it only has a thickness of 0.6 mm.
- the MEA is framed by two nonmagnetic steel plates, a graphite layer, a thin titanium layer (integrated into the steel plates) and a Teflon frame (omitted in Figure 1).
- the negative pole and at the bottom right in FIG. 1 the positive terminal of the fuel cell is connected.
- the operating point is selected with 60 A current at 0.4 A voltage drop so that the electrical conductivity is homogeneous over the entire MEA.
- the cell has a terminal voltage of 0.7 V when de-energized.
- FIG. 2 shows, by way of example, the standard of some components of column vectors from the matrix U for the fuel cell from FIG. 1 in the form of two-dimensional images, wherein the degree of blackening of the screening in the image indicates the absolute value of the standard.
- the parameter i indicates the column index of the vector.
- Magnetic field values is a linear combination of column vectors from the matrix U. Therefore, the exemplary representation of some of these column vectors in image form clarifies that different column vectors contribute information portions with strongly different spatial resolution to the final reconstruction. Therefore, it makes sense to choose a measure of relevance, which is a measure of the size of the
- the grid spacing of the rated locations is about 0.6 cm, so a total of 5536 locations were rated.
- the most relevant sites for the reconstruction according to the relevance measure R k are located on the front and on the rear side of the fuel cell. This is due to the fact that the true diameters of cylindrical flows can best be determined at these levels.
- the condition of the matrix S improves by the factor 1, 4.
- the logarithms of the singular values of the matrix S which includes the locations considered in FIG. 3a, (upper curve A)
- the logarithms of the singular values of the matrix S which only includes the 868 locations considered in FIG. 3b (lower curve B)
- large singular values decrease more rapidly than small ones, the reduction of the amount of data to be recorded via the magnetic field by a factor of 6.3 does not have to be paid for with a loss of resolution in the tomographically reconstructed current distribution.
- the logarithm of the relative precision with which the magnetic field measurement has to be carried out (ordinate) is above the desired spatial resolution the current distribution j to be reconstructed is plotted in centimeters (abscissa).
- the lower curve A with circles as symbols is based on a matrix S which comprises all 5536 locations from FIG. 3a.
- the upper curve B with asterisks as symbols is based on a matrix S which comprises only the 868 locations from FIG. 3b.
- the drastic reduction of the recorded measuring points leads to a given spatial resolution being achieved with more unreliable measured values for the auxiliary magnetic field or, for a given measuring accuracy for the magnetic field, achieving a better spatial resolution in the reconstructed current distribution.
- Cost of data acquisition and quality of tomographic reconstruction can be selected, comparable to the lossy compression of image or audio files.
- the freely selectable parameter Nc determines the fineness of the discretization and thus the accuracy with which each of the N c measured values m is spatially localized.
- S kj are the elements of a matrix representation S of the linear operator S.
- S k j is an abbreviation for S x , j.
- W is the curve, area or volume within which measurements are taken.
- a coordinate (i or j) corresponds to an index of the desired component [S + S] y of the matrix product [S + S].
- the other coordinate is set according to the location x over which is integrated in the convolution.
- p (x) can be understood not only as a density function but also as a measure of relevance. This assigns a weight to each location x, with which measured values obtained in this location should enter into the tomographic reconstruction, in order to maximize the quality of this reconstruction with regard to the functional used in the optimization of p (x).
- the physically meaningful standardization becomes
- the A k are the amplitudes of the ⁇ -functions, which represent their weights in the density function p (x).
- the integral notation makes it possible, for example, to interpolate measured values between the locations X k at which they were actually recorded and to smooth out the density function p (x) obtained therefrom.
- [S + S] of p (x) can be calculated. Due to the specific task, a mathematical condition for the eigenvalue spectrum of [S + S] is given, which should be optimally fulfilled. As a rule, such a condition exists in the form of a function which depends either on the whole matrix [S + S] or on one or more of its elements. In this functional, the above representation, in which [S + S] is expressed as an integral over p (x), is used. The functional of [S + S] thus becomes a functional of p (x), that is, a functional of the sought distribution of measuring locations at which the values of the auxiliary variable are expediently to be determined.
- each location x in the convolution is assigned a measure of how strongly changes in the measured variable m within the sample affect the auxiliary value H (x) measured at this location, ie how relevant is the location x for the tomographic reconstruction.
- each location x in the convolution can be assigned a measure of how much information a value H (x) 0 of the auxiliary variable measured at that location over the total distribution m ( x) contains the measured variable.
- T n is the unitary transformation that diagonalizes [S + S] nI .
- the matrix product in the integral can be simplified so that
- ⁇ n (x) is the correction function.
- the density function p (x) is thus changed at the transition from the iteration n-1 to the iteration n by a pointwise product of itself with the correction function.
- the correction function ⁇ "(x) is required to be orthonormal to ⁇ n- i (x), ie
- any such functional which contains an eigenvalue of the matrix [S + S] also contains at least one element of this matrix.
- the unitary transformations of T [S + S] are represented as linear combinations of the elements of [S + S] by the unitary transformation T n . Since both these elements and the elements of T n can in turn be expressed as a convolution of the distribution p (x) of measurement locations with a concatenation of the operators S and S + , as described above, each functional contains an eigenvalue of the matrix [ S + S], thus automatically including such a convolution.
- the eigenvalues of [S + S] are the singular values of ⁇ J p (x) • S.
- the invention therefore also relates to a method for determining the optimal distribution of measuring locations at which an auxiliary variable accessible to the measurement is to be measured within the sample for tomographic reconstruction of a measured variable, the distribution of measuring locations being optimized to the extreme value of a function which contains at least one singular value of an operator which contains a product of a power of the distribution of measurement locations with the operator S, which assigns a distribution of the auxiliary variable to each possible distribution of the measured variable.
- the potency in this case is 1 A.
- the optimization goal is a maximum absolute accuracy
- the goal is to make the smallest, ie the Nc'th, eigenvalue of [S + S] as large as possible. Therefore, the change of this eigenvalue should be maximized.
- An example of a functional describing the step from the iteration n-1 to the iteration n with respect to this optimization goal is
- the first sum is the change of the Nc'th eigenvalue at the transition from the iteration n-1 to the iteration n in the 1st order perturbation theory.
- ⁇ 0 and ⁇ i are Lagrange parameters that should be chosen so that the constraints are met.
- This functional can be generalized to the effect that when optimizing the smallest eigenvalue of [S + S], the effects of the changes made to [S + S] (via the change of the density function p (x), which causes [S + S]) to the other eigenvalues of [S + S]. be taken into account. This can be avoided in particular that occurs during the iterations degeneracy. In the case of degeneracy, the change of eigenvalues from one iteration to the next can no longer be approximated in order of perturbation theory as described above.
- the first two summation terms are a measure of the change in the ratio between the largest and the smallest eigenvalues of [S + S] and thus also a measure of the ratio between the largest and the smallest singular value of S.
- the change was the same as in the optimization goal the maximum absolute accuracy in 1st order perturbation theory. hert.
- ⁇ n (x) has to be chosen in order to steer the center of gravity ⁇ s in the desired direction.
- the optimization goal is to be formulated as a functional, which depends on this focus ⁇ s.
- the focus .lambda..sub.s can be, as described above, n (x) described by the expansion coefficients ⁇ v ausdrü- CKEN that ⁇ , the composition of the correction function of the basic functions of d v.
- these evolution coefficients are ultimately the only free parameters. Secondary conditions can be considered in such a function by terms with Lagrange parameters.
- Each such functional which contains the center of gravity ⁇ s, also depends on at least one element of the matrix [S + S], since the eigenvalues depend on the matrix elements.
- These matrix elements can in turn, as described above, be expressed as a convolution of the distribution p "(x) of measurement locations with a function f (S) of the operator S (see the above representation for the elements [S + S], j ). Therefore, every functional that contains an eigenvalue of [S + S], and certainly every functional that contains a centroid of a set of such eigenvalues, automatically contains such a convolution.
- the goal is to make the center of gravity ⁇ s > u of the smallest eigenvalues of [S + S] and thus also the change ⁇ s , u of this center of gravity as large as possible.
- the smallest eigenvalues of [S + S] are considered only those eigenvalues of [S + S] which are at most a factor p larger than the smallest eigenvalue of [S + S].
- the smallest possible values of p should therefore be sought.
- the value of p is typically on the order of 1.1.
- the first sum is proportional to the change of ⁇ s , u in the transition from the iteration n-1 to the iteration n in the first order perturbation theory.
- the value of this function should be maximized as a function of ⁇ v .
- ⁇ o and ⁇ are Lagrange parameters that are to be chosen so that the constraints are met.
- the optimum obtained is thus a function of these Lagrangian parameters.
- the end result is the optimum obtained for those values of the Lagrange parameters that satisfy the constraints.
- the iterations must be terminated as soon as the optimum of F n over a number of iterations does not yield a positive value for ⁇ s jU , because then no improvement can be expected.
- the value of ⁇ inserted in ⁇ s > u determines the strength of the correction and thus controls the speed of the optimization process. Since, due to the orthonormality condition for the correction function ⁇ n (x), this function must change its sign, for each such correction function ⁇ n (x) there is a maximum ⁇ max > 0, at which point the approach described above for the transition from Iteration n-1 no longer works on the iteration n.
- this goal can be achieved, for example, via the functional
- ⁇ s u and ⁇ s 0 are respectively the centers of gravity of the smallest and the largest considered eigenvalues of [S + S]. Analogous to the amount of the smallest in
- the eigenvalues considered are only eigenvalues considered to be the largest eigenvalues, which are at most a factor q below the largest eigenvalue. Typical values for q are in the order of 0.9. The smaller q and the larger the bandwidth of those eigenvalues, the largest eigenvalues are considered, the more stable is the focus of this amount of largest eigenvalues and the less accurate it is at the same time determined.
- FIG. 7a shows the development of the largest and FIG. 7b the development of the 5 smallest eigenvalues of [S + S] for an operator S, which is typical for the problem of the tomographic reconstruction of the current distribution in a fuel cell (measured variable) from the outside of the fuel cell Magnetic field values (auxiliary size).
- the optimization procedure was performed iteratively, and the optimization goal was to increase the heavy point ⁇ s > u of the smallest eigenvalues of [S + S] as large as possible. Plotted are the eigenvalues over the number of iterations.
- the eigenvalues were determined after each iteration n by performing a singular value decomposition of the matrix [S + S] n resulting from the respective distribution p n (x) of measurement locations.
- the smallest eigenvalue of [S + S] is increased during the iterations by a factor 5 of about 8, which seems to be a typical order of magnitude for the said application on the fuel cell.
- the number of measuring locations can be drastically reduced for the exemplary embodiment of the optimization method according to the invention for the fuel cell shown in FIG.
- p (x) of measurement locations on an outer surface, which has a distance of 1 cm from the fuel cell in each dimension 5.
- the distribution was determined using the optimization method according to the invention for a fixed number NH of 6,700 measurement locations, wherein a uniform distribution of the measurement locations over the entire outer surface was selected as the starting value for the first iteration. After optimization, all measuring locations were rejected where p (x) given threshold undercut.
- the threshold value chosen in this case was the average value of p (x).
- the remaining measurement sites corresponded to a reduced density P R (X) of measurement sites, each of which received a ⁇ function in this density. Retaining the number and positions of these ⁇ functions, the reduced density P R (X) has now been reduced to the physically meaningful
- Layer 1 is steel-titanium
- layer 2 is the membrane-electrode assembly (MEA)
- layer 3 is graphite
- layer 4 is steel.
- the remaining measuring locations are marked by dots, the area contents of the dots being greater, the greater the weight assigned to the corresponding measuring location by the second optimization.
- the number of measuring locations could be reduced by 93% and at the same time the absolute accuracy increased by a factor of 2.77, with the remaining measuring locations being given very different weights by the second optimization have received. At the same time, the relative accuracy could be increased by a factor of 2.8.
- the practical implementation of the measurement can be further simplified by taking into account only the measuring sites on the one side surface where the vast majority of measurement sites remaining after the first optimization are concentrated. Then enough a single scanning detector that processes this side surface to tomographically determine the three-dimensional flow distribution within the fuel cell.
- a scanning detector which processes the measuring locations one after the other, can be used in this exemplary embodiment of a fuel cell because time-independent stationary operating states are set there.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102007044257A DE102007044257A1 (de) | 2007-09-17 | 2007-09-17 | Verfahren zur tomographischen Rekonstruktion einer Messgröße aus einer der Messung zugänglichen Hilfsgröße |
PCT/DE2008/001384 WO2009036726A1 (de) | 2007-09-17 | 2008-08-20 | Verfahren zur tomographischen rekonstruktion einer messgrösse durch gezielte bestimmung der messpunkte |
Publications (2)
Publication Number | Publication Date |
---|---|
EP2191444A1 true EP2191444A1 (de) | 2010-06-02 |
EP2191444B1 EP2191444B1 (de) | 2017-07-05 |
Family
ID=40112839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP08801207.5A Not-in-force EP2191444B1 (de) | 2007-09-17 | 2008-08-20 | Verfahren zur tomographischen rekonstruktion einer messgrösse durch gezielte bestimmung der messpunkte |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP2191444B1 (de) |
DE (2) | DE102007044257A1 (de) |
WO (1) | WO2009036726A1 (de) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102012016837A1 (de) | 2012-08-27 | 2014-02-27 | Forschungszentrum Jülich | Vorrichtung und Verfahren zur qualitativen Bestimmung des Betriebszustandes eines Prüfobjekts |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19812285A1 (de) * | 1998-03-20 | 1999-09-23 | Philips Patentverwaltung | Bildgebendes Verfahren für medizinische Untersuchungen |
DE102006035741A1 (de) * | 2005-07-28 | 2007-05-03 | Karl-Heinz Dr. Hauer | Verfahren zur Bestimmung der Stromdichteverteilung in einem Brennstoffzellenstack |
EP1780672A1 (de) * | 2005-10-25 | 2007-05-02 | Bracco Imaging, S.P.A. | Verfahren zur Registrierung von Bild-Algorithmen zum Durchführen der Bildregistrierung, Programm zur Registrierung, das den Algorithmus verwendet und Verfahren zur Behandlung von biomedizinischen Bildern zur Verminderung von durch Objektbewegung verursachten Bildgebungsartefakten |
-
2007
- 2007-09-17 DE DE102007044257A patent/DE102007044257A1/de not_active Withdrawn
-
2008
- 2008-08-20 DE DE112008003156T patent/DE112008003156A5/de not_active Withdrawn
- 2008-08-20 EP EP08801207.5A patent/EP2191444B1/de not_active Not-in-force
- 2008-08-20 WO PCT/DE2008/001384 patent/WO2009036726A1/de active Application Filing
Non-Patent Citations (1)
Title |
---|
See references of WO2009036726A1 * |
Also Published As
Publication number | Publication date |
---|---|
WO2009036726A1 (de) | 2009-03-26 |
DE102007044257A1 (de) | 2009-03-19 |
EP2191444B1 (de) | 2017-07-05 |
DE112008003156A5 (de) | 2010-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2165158B1 (de) | Messverfahren, sensoranordnung und messsystem | |
DE102013220933B3 (de) | Shim-Verfahren mit Festlegung der Zielfeldverteilung durch Optimierung in einem Parameterraum reduzierter Dimension | |
DE102009014723B4 (de) | Kontrastabhängige Regularisierungsstärke bei der iterativen Rekonstruktion von CT-Bildern | |
EP3195258A1 (de) | Verfahren und system zur bestimmung der lokalen qualität von aus volumendaten extrahierten oberflächendaten | |
DE102005028475B4 (de) | Verfahren und Vorrichtung zur Bestimmung von Koeffizienten eines Diffusionstensors mittels magnetischer Resonanz | |
DE102009050662A1 (de) | Verfahren und Vorrichtung zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung für eine Magnetresonanzanlage | |
DE102010001577A1 (de) | Verfahren zur Reduktion von Verzerrungen in der Diffusionsbildgebung und Magnetresonanzanlage | |
DE102009028104A1 (de) | Bildrekonstruktionsverfahren für energiereiches Doppelenergie-CT-System | |
DE10317384A1 (de) | Computertomographie | |
EP2454601A1 (de) | Verfahren und anordnung zum rekonstruieren der quelle eines elektromagnetischen feldes | |
DE69724681T2 (de) | Vorrichtung und Verfahren zur Analyse von Kristallstruktur | |
EP3649614B1 (de) | Verfahren zur bestimmung von unsicherheiten in messdaten aus einer messung eines objekts | |
DE602004012290T2 (de) | Verfahren zum automatischen Shimmen für die Kernspinresonanzspektroskopie | |
DE102009010501A1 (de) | Verfahren zur Rauschreduktion von CT-Bilddaten und Bildbearbeitungssystem | |
DE602005004506T2 (de) | Entfernen störender Geräteeinflüsse von einer Beugungsstruktur durch Entfaltung, wobei eine beugungswinkelabhängige Gerätefunktion verwendet wird | |
DE102012223789B4 (de) | Ermittlung einer Phasendifferenz-Karte | |
DE102019127795A1 (de) | Schaltung und ein Verfahren zum Bestimmen einer Lage eines Magneten und Joystick | |
DE102017205113A1 (de) | Ermitteln der Pose einer Röntgeneinheit relativ zu einem Objekt anhand eines digitalen Modells des Objekts | |
EP2191444B1 (de) | Verfahren zur tomographischen rekonstruktion einer messgrösse durch gezielte bestimmung der messpunkte | |
WO2014053256A1 (de) | Isotrope rekonstruktion von bilddaten | |
DE102017210700A1 (de) | Verfahren zum automatisierten Quantifizieren eines Analyten sowie NMR-Messgerät zur Durchführung des Verfahrens | |
DE112019006092T5 (de) | Lose gekoppeltes inspektions- und metrologiesystem zur überwachung eines produktionsprozesses mit hohem volumen | |
DE60037235T2 (de) | Verfahren zur ermittlung eines wahren spektrums aus einem gemessenen spektrum | |
DE4309958C1 (de) | Verfahren und Vorrichtung zur ortsauflösenden Magnetresonanzuntersuchung eines Meßobjekts | |
DE102014226052A1 (de) | Ermittlung einer phasenbeschreibenden Karte mit Hilfe eines min-cut/max-flow-Verfahrens |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20100409 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MT NL NO PL PT RO SE SI SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL BA MK RS |
|
RIN1 | Information on inventor provided before grant (corrected) |
Inventor name: REISSEL, MARTIN Inventor name: SCHMIDT, ULI Inventor name: LUSTFELD, HANS Inventor name: HIRSCHFELD, JULIAN Inventor name: STEFFEN, BERNHARD |
|
DAX | Request for extension of the european patent (deleted) | ||
17Q | First examination report despatched |
Effective date: 20140520 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
INTG | Intention to grant announced |
Effective date: 20170223 |
|
RIN1 | Information on inventor provided before grant (corrected) |
Inventor name: HIRSCHFELD, JULIAN Inventor name: LUSTFELD, HANS Inventor name: SCHMIDT, ULI Inventor name: STEFFEN, BERNHARD Inventor name: REISSEL, MARTIN |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MT NL NO PL PT RO SE SI SK TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D Free format text: NOT ENGLISH |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: REF Ref document number: 907071 Country of ref document: AT Kind code of ref document: T Effective date: 20170715 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D Free format text: LANGUAGE OF EP DOCUMENT: GERMAN |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 502008015433 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 10 |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20170705 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: HR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: NO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20171005 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20171005 Ref country code: IS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20171105 Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20171006 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 502008015433 Country of ref document: DE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: MM4A |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20170831 |
|
26N | No opposition filed |
Effective date: 20180406 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170820 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170820 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 11 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170831 Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO Effective date: 20080820 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CY Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20170705 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20170705 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20200824 Year of fee payment: 13 Ref country code: FR Payment date: 20200819 Year of fee payment: 13 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: IT Payment date: 20200831 Year of fee payment: 13 Ref country code: CH Payment date: 20200825 Year of fee payment: 13 Ref country code: AT Payment date: 20200818 Year of fee payment: 13 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R081 Ref document number: 502008015433 Country of ref document: DE Owner name: FACHHOCHSCHULE AACHEN, DE Free format text: FORMER OWNER: FACHHOCHSCHULE AACHEN, 52428 JUELICH, DE Ref country code: DE Ref legal event code: R081 Ref document number: 502008015433 Country of ref document: DE Owner name: FACHHOCHSCHULE AACHEN, DE Free format text: FORMER OWNERS: FACHHOCHSCHULE AACHEN, 52428 JUELICH, DE; FORSCHUNGSZENTRUM JUELICH GMBH, 52428 JUELICH, DE |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MM01 Ref document number: 907071 Country of ref document: AT Kind code of ref document: T Effective date: 20210820 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20210820 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210831 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210831 Ref country code: AT Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210820 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IT Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210820 Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210820 Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20210831 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20220720 Year of fee payment: 15 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 502008015433 Country of ref document: DE |