Method for checking the thickness conformity of a radome

Abstract

Method for checking the thickness conformity of a radome. The invention relates to a method for checking the thickness conformity of a radome, comprising: - the definition of a reference radome, comprising a first mesh (M1) of first points, each associated with an average (μk) of thicknesses of conforming radomes, and a first control window (Fk) encompassing the first point and adjacent first points, - the measurement at each second point of a second mesh modeling the radome to be checked, of a thickness at this second point, associated with a second control window, - the checking of each second point, in which if a distance between a vector of measurements of the second points of the associated second control window and a vector of averages (μk) associated with the first points of the corresponding first control window (Fk), with a first threshold defining a confidence interval, is greater than a first threshold,The second point is non-compliant. (Figure 3)

Description

Title of the invention: Method for checking the thickness conformity of a radome

[0001] The present invention relates to the field of automobiles, and in particular to the field of conformity control of radomes protecting radars on board vehicles, these radars being used by advanced driver assistance systems.

[0002] Motor vehicles with advanced driver assistance systems, such as emergency braking devices or autonomous driving software modules, are equipped with radars, enabling the detection of obstacles or objects outside these vehicles.

[0003] In such a vehicle, several radars are, for example, positioned at the front, rear, and sides of the vehicle. Each radar is protected from the vehicle's external environment by a radome. As a reminder, a radome is a waterproof protective enclosure used to shield the radar; it is made of a lightweight, strong, non-metallic, and waterproof material such as fiberglass, but whose essential property is to minimize the attenuation of the emitted and received radar signal; its reflectivity, which is the ability to reflect radar waves, must therefore be very low. This reflectivity is, of course, related to the thickness of the radome.

[0004] It is understood that it is essential to control the thickness of such a radome at every point in the field of vision of the radar it protects, in order to guarantee the reliability of the measurements provided by the radar and consequently the proper functioning of the advanced driver assistance systems.

[0005] Today, the conformity check of the thickness of such a radome is carried out by thickness measurements taken (carefully) with a probe at twelve points within the radar's field of view, supplemented by measurements of the reflectivity or absorption by the radome of the radar waves arriving at it. In particular, precise measurements of radar wave reflectivity are used, performed in anechoic chambers (to eliminate any interference with the radar). However, the cost of these conformity certification measurements for thickness and reflectivity in anechoic chambers is high, even though these measurements are nevertheless absolutely essential for quality control.

[0006] There is therefore a need to design a method for checking the thickness conformity of a radome, which does not require an anechoic chamber and is inexpensive, particularly in terms of time and computing resources.

[0007] To this end, the invention proposes a method for checking the thickness conformity of a radome, characterized in that it comprises: - a step of defining a statistical model of a conformal radome, called a reference radome, comprising a first mesh of first points, the first mesh characterizing a radar field of view associated with the reference radome, a step in which each first point of the first mesh is associated with: - an average of a series of representative data on conformal radome thicknesses, associated with radar field-of-view points of these conformal radomes, which correspond to the first point, and - a first surface control window encompassing the first point, and the first points of the first mesh adjacent to the first point, the first surface control window being associated with a covariance matrix between the representative data series of conformal radome thicknesses whose means are each associated with one of the first points of the first surface control window, - a measurement step at each second point of a second mesh modeling a radar field of view of the radome to be controlled, of a quantity representative of a thickness of the radome to be controlled associated with this second point, to which is associated a second surface control window encompassing said second point and second points of the second mesh adjacent to said second point, - a control step for each second point of the second mesh, comprising a comparison of a distance, which is a function of the covariance matrix, between on the one hand a vector of measurements of the second points of the second surface control window associated with the second point to be controlled and on the other hand a vector of means associated with the first points of the first surface control window corresponding to the first point located in the same location in the first mesh as the second point to be controlled in the second mesh, the second surface control window associated with the second point to be controlled being geometrically identical to the corresponding first surface control window, with a first threshold defining a confidence interval, and a qualification of said second point to be controlled as compliant or non-compliant depending on the result of the comparison, - a step involving counting the number of second-point inspections deemed compliant or the number of second-point inspections deemed non-compliant, and - a step of comparing the number of second points counted with a second threshold, resulting in a qualification of the radome to be checked as compliant or non-compliant.

[0008] Thanks to the invention, measurements in an anechoic chamber are only carried out to ensure the conformity of the conforming radomes used subsequently during the definition step of the reference radome. This definition step is only carried out once, the Subsequent steps can then be repeated as many times as needed on radomes to be inspected, for example, at the end of a production line. These radomes to be inspected are the subject of the measurement step, which is carried out using a measurement system in the factory and outside of an anechoic chamber. Each radome to be inspected is then compared to the reference radome, previously constructed based on measurements taken with the same measurement system as the radome to be inspected, i.e., outside of an anechoic chamber. These measurements are taken on radomes certified as compliant in thickness and used to construct the reference radome.

[0009] It should be noted that, compared to a control method requiring a reference radome defined by deep learning, the reference radome constructed according to the invention requires much less computing resources and much less measurement (a simple sample of 30 conforming radomes is sufficient).

[0010] In the method according to the invention, the measurement vector corresponds of course to a vector of representative quantities measured for the second points of the second surface control window, during the measurement step.

[0011] The representative quantity of a thickness of the radome to be monitored is, for example, an average reflectivity over a frequency band traversed by a radar sending electromagnetic waves onto the radome to be monitored, and the representative data for the thicknesses of conformal radomes are each characteristic of such an average reflectivity. This average reflectivity is therefore representative of the reflectivity of radar waves in this frequency band, specific to the radome to be monitored. Alternatively, the representative quantity of a thickness of the radome to be monitored is an average absorbance of radar waves traversing a frequency band.

[0012] Furthermore, the distance between the vector of measurements and the vector of means is, for example, a Mahalanobis distance. The first threshold is, for example, an upper control limit corresponding to a confidence interval based on a Fisher transform applied to the Mahalanobis distance. The second point to be checked during the inspection step is considered compliant when the distance is less than or equal to this first threshold. This confidence interval is commonly accepted at 99.8% in quality control. Other transforms can, of course, be used, for example, a χ² or beta distribution.

[0013] The second threshold is, for example, an upper control limit corresponding to a confidence interval based on a Poisson distribution. This confidence interval is, for example, 99.8% with regard to the objective regulatory non-compliance requirement for thickness targeted in production; this requirement is obviously very low.

[0014] According to an optional and advantageous feature of the method for checking the thickness conformity of a radome according to the invention, the representative data conforming radome thicknesses correspond at least in part to representative measurements of conforming radome thicknesses, and the step of defining a statistical conforming radome model includes a step of removing, in these representative measurements of conforming radome thicknesses, at least one outlier measurement associated with a point in a radar field of view of a conforming radome corresponding to a first point of the first mesh, and replacing said at least one outlier measurement with an average of non-outliers measurements associated with points in radar fields of view of other conforming radomes corresponding to said first point of the first mesh.

[0015] Since the conforming radomes have been certified as conforming in thickness in an anechoic chamber and by thickness measurements carefully taken with a probe, the radome remains somewhat fragile. Therefore, an aberrant measurement at a point on a measured conforming radome can be replaced by a non-aberrant pseudo-measurement corresponding to the average of measurements taken at the same location on the other conforming radomes, without calling into question the conformity of the conforming radome. This makes it possible to correct any measurement errors that could degrade the representativeness of the reference radome and also reduces the variability of the measurements.

[0016] The removal of said at least one outlier uses, for example, a Grubbs test. Of course, other types of outlier detection tests can be used.

[0017] According to another optional and advantageous feature of the method for checking the thickness conformity of a radome according to the invention, the data representing the thicknesses of conforming radomes correspond partly to measurements or averages of measurements representing the thicknesses of N conforming radomes, and partly to virtual measurements associated with L virtual radomes determined during a step of determining these L virtual radomes, comprising a phase of generation by random draws representative of normal laws, of series of values ​​each associated with a first distinct point of the first mesh.

[0018] Indeed, to obtain a reasonable first Mahalanobis distance threshold at a 99.8% confidence level between the measurements of the radome to be tested and the control radome, the latter must be composed of a sufficient number of compliant radomes. However, since compliant radomes must be certified using an anechoic chamber and delicate thickness measurements performed with a probe, few certified compliant radomes are actually available—typically thirty—to construct the control radome. Given the surface area of ​​the control window, which is defined by measurement points adjacent to the point being tested, when a "large" area is chosen, such as 25 points (a 5x5 square), one encounters a problem with the Fisher transform statistic. associated with the Mahalanobis distance, a small number of degrees of freedom in the denominator then induces a very large distance threshold value, naturally impairing the detection power of a surface non-conformity. In this case, the reference radome must be enriched with L virtual radomes having the statistical characteristics of the N conforming radomes of the reference radome, L being large enough to establish a 99.8% confidence level with respect to the Fisher transform with a reasonable initial Mahalanobis distance threshold.

[0019] In this step of determining these L virtual radomes, assuming, for example, that the normal distributions are centered and reduced, the generation phase is followed by a correction phase for the series of values ​​generated by random sampling. This correction ensures that the mean of each series of corrected values ​​associated with each first point of a first surface control window of the first mesh is zero, and that the covariances between the series of corrected values, each associated with a distinct first point of the first surface control window, form the identity matrix. This correction phase eliminates residual means and covariances, making it easier to then work on the associated series of values ​​to determine the mean and covariance values ​​of the N conforming radomes of the reference radome.

[0020] The correction phase is then followed by a phase of generating virtual measurements, each associated with one of the L virtual radomes, in which the virtual measurements associated with the same point of the L virtual radomes corresponding to a first point of the first control window are obtained as a function of the series of corrected values ​​associated with this first point, an average of measurements representative of the thicknesses of the N conforming radomes corresponding to this first point, and the covariances between the series of data representative of the thicknesses of the N conforming radomes each associated with a first distinct point of the first surface control window.Alternatively, it is possible to work on all the series of values ​​associated with all the first points of the first mesh simultaneously, rather than working first surface control window by first surface control window, if sufficient computing resources are available during the construction stage of the test radome.

[0021] Other features and advantages of the invention will become apparent from the following description on the one hand, and from several illustrative and non-limiting examples of embodiments given by reference to the accompanying schematic drawings on the other hand, in which:

[0022] [Fig. 1] represents a first step of a conformity control method according to the invention, for the thickness of a radome, in an embodiment of the invention, this first step being the definition of a test radome,

[0023] [Fig.2] represents measurement data on meshes corresponding to radar fields of view through conforming radomes comprising at least part of the reference radome, used during the first stage of the conformity check process of the [Fig.1],

[0024] [Fig.3] represents data associated with a radar field of view of the radome witness as modeled at the end of this first step,

[0025] [Fig.4] represents other steps in the conformity control process of [Fig.1], And

[0026] [Fig.5] represents measurement data on a mesh corresponding to a field radar vision through an unknown radome to be checked, during a second step of the conformity check process of the [Fig.l].

[0027] A method for checking the thickness conformity of a radome according to the invention, shown in [Fig. 1], comprises a first step 10 of defining a statistical model of a conforming radome, called a reference radome. This first step 10 is performed only once, and the reference radome resulting from this first step 10 can then be used to check several radomes.

[0028] This first step 10 comprises a first substep 102 of reflectivity measurements of conformal radomes, i.e., radomes each capable of housing a radar without interfering with its detection of obstacles located outside the radar. In other words, the thickness at every point of the field of view of a conformal radome is within a tolerance range suitable for radar detection of obstacles through the conformal radome.

[0029] The conformity of these radomes has for example been established by precision measurements in an anechoic chamber and previously by conformity measurements with the "probe" at a few points (usually 12) of the radar field of view.

[0030] In this first substep 102, the Rhode & Schwarz® QAR (Quality Automotive Radome Tester) system is used, for example, on N = 30 conformal radomes, each having an overall ellipsoidal radar field of view. N is, of course, a natural number strictly greater than 1. This system provides, for each i-th of the N conformal radomes, i being between 1 and N, measurements of the reflectivity or absorption of a radar wave through the i-th conformal radome, at a plurality of points, corresponding to K points of a mesh M of the overall elliptical radar field of view, shown in [Fig. 2]. The mesh M is geometrically the same for the N conformal radomes. Of course the geometric shape of the mesh depends on the radar, its field of view and the shape of the type of radome whose conformity we want to control, other forms of surface mesh of the radar's field of view are therefore also usable in the invention.

[0031] The system performs a frequency scan of the radar waves used for these measurements and provides an average over these frequencies of the reflectivity at each point k of the mesh M for each i-th conformal radome. Thus, at the output of this first step 102, for each i-th conformal radome, we have K reflectivity measurements rk>i with k between 1 and K (maximum number of points in the mesh M) and i between 1 and N, rkji being representative of a reflectivity over a frequency band. K is approximately 27400 in this embodiment of the invention.

[0032] It should be noted that for visibility reasons some parameter letters are indexed in the equations but are not necessarily indexed in the figures and the text.

[0033] A second substep 104, following the first substep 102, is the removal of outliers from the first substep 102, when the measurements from this first substep 102 include outliers. This step is optional but beneficial to the end customer, as it reduces the acceptable dispersion of the measurements used in the construction of the reference radome, made up of the conforming radomes, and thus improves the quality of the conformity check of the unknown radomes at the end of production.

[0034] In this embodiment of the invention, this second substep 104 uses a Grubbs test to detect outliers in the measurements from the first substep 102. For this purpose, the series of N measurements rk>i, with i ranging from 1 to N, at a k-th point of the mesh M, is approximated by a probability distribution Rk, with mean

[0035] and standard deviation _ hy^ ( .2 mk- —Vmk)

[0036] Next, for each i-th conforming radome and each k-th point of the mesh M, the following values ​​are defined:

[0037] 7 and / A1V-21Z,2 is statistically assimilable for the k- !..-L Æ'' ^amU / VZ^) The i-th point of the mesh M and over the set of N radomes follows a Student's t-distribution with N-2 degrees of freedom. By setting a risk α, for example at 10 parts per million (0.001%), we define a confidence interval 1-α. The corresponding upper control limit LSC(k) is 5.12 in this case (for N=30).

[0038] In this second substep, if > f^C(k) ai°rs the measurement provided for the k-th point of the i-th radome is an outlier, it is removed from the set of measurements from the first substep 102. Each outlier is replaced, in this set of measurements from the first substep 102, by the average r'^i of the non-outliers corresponding to the same k-th point carried out on the other conforming radomes of the set N of conforming radomes. Thus, in the case where only one outlier is detected for this k-th point, we calculates an average of the Nl non-aberrant measurements at the k-th point. At the end of this second substep 104, we thus obtain a set of measurements, some of which are corrected, associated with the N conformal radomes and the K points of the mesh M.

[0039] Of course, alternatively, in this detection of outliers, one can use fewer than N conforming radomes, for example a subset of the N conforming radomes, to perform the calculations of the values ​​^kj and the upper control limit is then adapted accordingly.

[0040] In a third substep 106 following the second substep 104, a statistical model of a conformal radome is defined, determined from the set of measurements possibly corrected, from the second substep 104. This statistical model includes a first mesh M1, represented [Fig.3], geometrically identical to the mesh M, therefore also including K points.

[0041] To each k-th point of the first mesh Ml is associated the average qk of the possibly corrected measurements provided for this k-th point at the end of the second sub-step 104:

[0042] *

[0043] With r'^ = if < LSC(k) or a mean of the non-outliers provided for the kth point otherwise.

[0044] Furthermore, a first surface control window Fk is associated with each k-th point of the first mesh Ml, consisting of the k-th point and Pl points adjacent to the k-th point, that is, near or in the vicinity of this point k. In the example in [Fig. 3], P is equal to 9, and the first surface control window is a portion of the first mesh Ml comprising the k-th point and each of its eight nearest neighbors surrounding this k-th point. Of course, P can be chosen to be equal to 25 or 49, for example, by enlarging the area of ​​the first surface control window Fk centered on the k-th point. Note that the larger the control area P, the better the "local" thickness conformity diagnosis performed by the surface control window, which relies on matrix calculations but obviously at the cost of increased computation time.

[0045] Associated to this first surface control window Fk is an inverse covariance matrix E*, which is the inverse matrix of the covariance matrix of dimensions P*P, where each element in row u and column v is equal to Cov(R'u, R'v), where R'u and R'v are the probability models to which are associated the series of N possibly corrected measures / ' \ and f \ at the u-th and v-th points ^i=làN \v^î=làN of the first surface window Fk. We therefore have:

[0046] Cov(R-.R-,)= , ArV w "u A' "v!

[0047] As will be seen in a second step 12 of the conformity control process 1, the inverse covariance matrix EÀ1 can be considered a metric associated with a Hotelling T2 distribution, to which is associated an upper control limit LSC(T2), and an acceptance threshold of 1-a' commonly equal to 99.8% confidence in Quality Control, based on a transform of the Fisher distribution with (P, NP) degrees of freedom. This Fisher distribution is defined by:

[0048] F = ((NP) / ((N-1)*P))* T2

[0049] In order to have a 1-a' confidence interval of 99.8% for a reasonable upper control limit LSC(T2), i.e. of the order of 100 for the K = 27400 points of the first mesh Ml, it is necessary to have at least 21 conforming radomes if we choose P=9, but at least 72 conforming radomes if we choose P=25.

[0050] The number P of points surrounding each first point of the first mesh Ml in its first surface control window Fk is chosen with regard to the computation time and cost of thickness conformity checks on the N conforming (physical) radomes; it should be noted that the number P=25 is the best compromise between cost and computation time for the sample of N=30 conforming radomes composing the control radome. It should also be noted that the larger P is, the greater the detection power for a thickness non-conformity at a point to be checked on a radome, this point to be checked corresponding to the first point.The invention is not limiting with respect to the quantity P of points forming the cloud of points neighboring the first point defining its first surface control window Fk, but a high value of P obviously requires increasing the number N of conforming radomes, which naturally has a cost (anechoic measurements + thickness conformity measurements with the probe). That said, this cost remains much lower than that which would be imposed by a training sample for a deep learning control method.

[0051] In the case where P=9, the statistical model defined in this third substep 106 constitutes the control radome used in the second step 12 of the conformity control process 1, since the N=30 control radomes used during these three substeps 102 to 106 are sufficient, and the next step of the conformity control process 1 is the second step 12.

[0052] In the case where P=25, the statistical model of the control radome must be enriched with L = 42 virtual radomes, which are determined (or constructed) during a fourth substep 108 of the first step 10 of the conformity control process 1. L is of course a natural number strictly greater than 1.

[0053] The fourth substep 108 of determining L virtual radomes is therefore an optional step which depends on the number of reference radomes used during the three Substeps 102 to 106, and especially the quantity P of points defining the size of the surface control window. It includes a first phase 180 of generation, using Monte Carlo simulations of standard normal distributions, of L values ​​for a first surface control window Fk of the first mesh Ml, and for each p-th point of the P points of this first surface control window Fk, by random sampling. These values ​​form a matrix (Xki>p) with i ranging from 1 to L and p from 1 to P.

[0054] For the first surface control window Fk, we then have a covariance matrix Sk defined in the same way as that of the statistical model defined in the third substep 106, this covariance matrix Sk being substantially equal to the identity matrix of dimensions P*P, and for each p-th point of the first surface control window Fk, a mean mkp of the measurements close to zero but not zero. The means mkp of the first surface control window Fk form a vector mk of p rows, each corresponding to a point of the first surface control window Fk.

[0055] A second phase 182 of the fourth substep 108 is then a correction phase in which a corrected matrix (X”ki>p) of values ​​associated with each first control surface window Fk is obtained, for which the mean of the L values ​​associated with each p-th point of the first control surface window Fk is strictly zero, and for which the covariance matrix, defined in the same way as that of the statistical model defined in the third substep 106, is strictly equal to the identity matrix. This corrected matrix (X”ki>p) of values ​​is obtained from the matrix (Xki>p) of values ​​obtained in the preceding first phase 180.

[0056] For this purpose, we remove from each value Xki p of the i-th row located in the p-th column of the matrix (Xki.p), the mean mkp, whose value depends of course on the p-th column, which results in a matrix (X'ki>p) with zero mean.

[0057] Then we obtain the corrected matrix (X' 'ki>p) by multiplying the matrix (X'ki,p) with zero mean by the inverse square root of the covariance matrix Sk, mathematically this is a centering and multidimensional reduction of the Monte Carlo sampling data:

[0058] (X”ki.p) = (X’ki.p)Sk1 / 2

[0059] It should be noted that the covariance matrix Sk, which is a residual covariance matrix, is symmetric positive definite, and that its inverse square root can therefore be calculated quite easily thanks to the P eigenvectors and P eigenvalues ​​of this residual covariance matrix.

[0060] The corrected matrix (X”ki>p) obtained at the end of the second phase 182 therefore corresponds to a sample of values ​​from a rigorously standardized multi-normal probability distribution. This multidimensional sample of standardized values ​​is propagated over the first mesh Ml over the first surface control window Fk+5 not overlapping with the first surface control window Fk, then over the first surface control window Fk+i0 and so on until the entire first mesh Ml is covered.

[0061] A third phase 184 following the second phase 182 is then a virtual measurement generation phase in which a virtual measurement matrix (Yki p) of L virtual measurement lines is calculated for a p-th point of the first surface control window Fk, of which: - the average over a line is ^p, that is to say identical to the average of the N=30 possibly corrected measurements obtained at the end of the second sub-step 104 of the first step 10 of the conformity control procedure 1, and - the covariance matrix defined in the same way as that of the statistical model defined in the third sub-step 106, is identical to this one, that is to say it is

[0062] For this, we perform the product of the corrected matrix (X”ki>p) by the square root of the covariance matrix and we add to each value of this product, located on a p-th column, the mean Pp, this one depending of course on the value p;

[0063] (Yki,) = (X”Ri.p) Zl'2+ / 1 \ “... I 1

[0064] The column vector on the right of this equation has L values ​​equal to one.

[0065] This virtual measurement matrix (Yki>p) is equivalent to measurements performed on L virtual radomes. Therefore, these are virtual measurements, although referred to as measurements hereafter for simplicity, in the remainder of this application. A measurement matrix is ​​formed in a similar manner for all the first surface control windows Fk+5m of the first mesh ML

[0066] It should be noted that it is possible to generate a larger sample of measurements from the first phase of this fourth substep 108, for example by working on a matrix of 50 draws per p-th point of a first surface control window Fk, then retaining only L = 42 measurements per p-th point of the first surface control window Fk in the measurement matrix (Yki>p) at the end of the third phase 184.

[0067] A fifth substep 110 of the first step 10 of the conformity control process 1 is then the aggregation of the measurements of the virtual radomes obtained during the fourth substep 108, with the measurements from the second sub- step 104, which are possibly corrected measurements associated with N=30 conforming radomes.

[0068] During this fifth substep 110, a reference radome is defined, which is an enriched statistical model of a conformal radome determined from all the measurements possibly corrected, from the second substep 104, and if necessary from all the measurements of the virtual L radomes obtained during the fourth substep 108. This reference radome includes the first mesh Ml.

[0069] For simplicity, the same designation dk4 is used for the measurements of these two sets, which are representative data of conformal radome thicknesses:

[0070] dkji = r'ki for i varying from 1 to N=30 and k from 1 to K=27400

[0071] and dkji = Ykijk for i varying from N+l to N+L, when k corresponds to the center of one of the first surface control windows Fk considered during the fourth substep 108 of determination of virtual radomes, otherwise dkji = Yqi>p where q is the q-th point of the first mesh Ml corresponding to a first surface control window Fq considered during the fourth substep 108 of determination of virtual radomes, and to which belongs the k-th point, p being an integer between 1 and P corresponding to the k-th point in this first surface control window Fq.

[0072] Each k-th point of the mesh Ml is associated with the mean qk of the sets of measurements obtained for this k-th point, identical to the average qk of the possibly corrected measurements r'^ from the second sub-step 104:

[0073] f*k “ N+L

[0074] Furthermore, to each k-th point of the first mesh Ml, we associate the first surface control window Fk, consisting of the k-th point and the Pl points located at neighborhood of the k-th point. In this example of an implementation of the invention, P is 25.

[0075] A covariance matrix is ​​associated with this first control window Fk reverse , which is the inverse matrix of the covariance matrix of dimensions (P,P) where each element in row u and column v is equal to Cov(Zu, Zv), where Zu and Zv are the probability distributions to which we approximate the sets of N + L data representing the thicknesses of radomes conforming to the u-th and v-th points of the first surface control window Fk. Therefore:

[0076] Cov(Zu, Zv) = _l_yN+L(du}(d. N+LL^i=l ( atlj l*u ) ( avj ivf )

[0077] The inverse covariance matrix can be considered a metric associated with a Hotelling T2 distribution, to which is associated an upper control limit LSC(T2), and therefore an acceptance threshold of 1-a' equal to 99.8% confidence, based on a The Fisher distribution is a transform of (P, N+LP) degrees of freedom. This Fisher distribution is defined by:

[0078] F = ((N+LP) / ((N+L-1)*P))* T2

[0079] The upper control limit LSC(T2) of this 1-a' 99.8% confidence interval is on the order of 100, which is therefore a "reasonable" threshold value.

[0080] The first step 10 of the thickness conformity control process 1 is thus finalized.

[0081] The next step (reference A) takes place when one wants to check the conformity of a radome for which no anechoic chamber is used, and therefore for which there is less possibility of detecting poor transmission of radar waves.

[0082] This next step, represented [Fig.4], is a step 12 of reflectivity measurements on the radome to be controlled, using the same system as in the first step 10, i.e. for example the Rhode&Schwarz®QAR system.

[0083] In this measurement step 12, the system provides a reflectivity measurement xk for each point k of a second mesh M2, shown in [Fig. 5], comprising K points. The second mesh M2 is geometrically identical to the first mesh ML. xk is a representative value of the reflectivity over the frequency band traversed by the Rhode & Schwarz® QAR system.

[0084] The next step 14 is a control step for each point k of the second mesh M2. A second surface control window Gk (represented in [Fig. 5] with only 9 points for simplicity) is associated with this point k, encompassing point k and points of the second mesh M2 adjacent to point k. This second surface control window Gk is geometrically identical to the first surface control window Fk. In this embodiment of the invention, we consider the case where the second surface control window Gk has P = 25 points and the first surface control window Fk also has P = 25 points. This embodiment can easily be transposed to the case where P is equal to 9 points.

[0085] This step 14 of checking point k includes a first substep 142 of calculating a distance Dk between a vector of measurements xk+p belonging to the second surface control window Gk and the data attached to the first surface control window Fk, p varying from 0 to 24. In this example of an embodiment of the invention, a Mahalanobis distance is used:

[0086] / xk~^k

[0087] A second substep 144 of step 14 of the k-point control is then the comparison of the calculated distance Dk with an upper control limit associated with the Mahalanobis distance, LSC(T2), here set at 99.78, given a size N=30 for the reference radome, a size L=42 for the virtual radomes, and an area of ​​P=25 points for the surface control window. Hotelling's T2 law refers here to the Mahalanobis distance. Other approximations are of course usable.

[0088] In other words, if Dk is strictly greater than LSC(T2) (branch N), the point k and the second surface control window Gk associated with it is qualified as non-compliant, otherwise (branch Y) the point k is qualified as compliant.

[0089] The step following step 14 of checking each point k is then a step 16 of counting the number Nb of points qualified as non-compliant in the previous step to judge the conformity of the "overall" thickness of the surface of the radome to be checked, a surface which corresponds to the field of vision of the radar.

[0090] The next step, 18, is a comparison of the number Nb of non-compliant qualified points with an upper control limit LSC(Nb) based on the Poisson distribution. For example, an objective non-compliance rate of 0.001% is chosen with regard to the European regulatory requirement Global Safety Regulation GSR2. The confidence interval is then 99.8%. Therefore, LSC(Nb) = 2.5. A decimal limit value is deliberately chosen, even though the count Nb is strictly an integer, to leave no doubt during the decision-making process.

[0091] In other words, if the number Nb of qualified non-compliant points is strictly greater than LSC(Nb) (branch N), the radome to be checked is qualified non-compliant, otherwise (branch Y) the radome to be checked is qualified compliant.

[0092] Of course, the invention is not limited to the examples just described, and many modifications can be made to these examples without departing from the scope of the invention. In particular, measurement systems other than the Rhode&Schwarz®QAR system can be used, the number N of conforming radomes can be different from 30, the number P can also take values ​​other than 9 or 25, and it is possible to use metrics other than the Grubbs test, the Mahalanobis distance, or the Poisson distribution.

Claims

1. Demands A method (1) for checking the thickness conformity of a radome, characterized in that it comprises: - a definition step (10) of a statistical model of a conformal radome, called a reference radome, comprising a first mesh (Ml) of first points, the first mesh (Ml) characterizing a radar field of view associated with the reference radome, a step in which we associate with each first point of the first mesh (Ml): - an average (qk) of a series of data representative of the thicknesses of conformal radomes, associated with points of the radar fields of view of these conformal radomes, which correspond to the first point, and - a first surface control window (Fk) encompassing the first point, and the first points of the first mesh adjacent to the first point, the first surface control window (Fk) being associated with a covariance matrix ( , ^½) between the series of representative data of conformal radome thicknesses, the averages of which are each associated with one of the first points of the first surface control window (Fk), - a measurement step (12) for each second point of a second mesh (M2) modeling a radar field of view of the radome to be controlled, of a representative quantity (xk) of a thickness of the radome to be controlled associated with this second point, to which is associated a second surface control window (Gk) encompassing said second point and second points of the second mesh (M2) adjacent to said second point, - a control step (14) of each second point of the second mesh (M2), comprising a comparison (144) of a distance (Dk), which is a function of the covariance matrix ( , between on the one hand, a vector of measurements of the second points of the second surface control window (Gk) associated with the second point to be controlled, and on the other hand, a vector of means (qk) associated with the first points of the first surface control window (Fk) corresponding to the first point located in the same position in the first mesh (M1) as the second point to be controlled in the second mesh (M2), the second surface control window (Gk) associated with the second point to be controlled being geometrically identical to the corresponding first surface control window (Fk), with a first threshold (LSC(T2)) defining a confidence interval, and a qualification of said second point to be checked as compliant or non-compliant depending on the result of the comparison (144), - a counting step (16) of the number of second points checked qualified as compliant or the number (Nb) of second points checked qualified as non-compliant, and - a comparison step (18) of the number (Nb) of second points counted with a second threshold (LSC(Nb)), resulting in a qualification of the radome to be checked as compliant or non-compliant.

2. Method of checking (1) the thickness conformity of a radome according to claim 1, wherein the representative quantity (xk) of a thickness of the radome to be checked is an average reflectivity over a frequency band traversed by a radar sending electromagnetic waves to the radome to be checked, and the representative data (r'k>i, dk>i) of conforming radome thicknesses are each characteristic of such average reflectivity.

3. Method for checking (1) the thickness conformity of a radome according to claim 1 or 2, wherein the distance (Dk) between on the one hand the measurement vector and on the other hand the average vector (qk), is a Mahalanobis distance.

4. Method of checking (1) the thickness conformity of a radome according to any one of claims 1 to 3, wherein the first threshold (LSC(T2)) is an upper control limit corresponding to a confidence interval based on a transform of Fisher's law applied to the Mahalanobis distance, the second point to be checked during the control step being qualified as conforming when the distance (Dk) is less than or equal to this first threshold (LSC(T2)).

5. Method of checking (1) the thickness conformity of a radome according to any one of claims 1 to 4, wherein the second threshold (LSC(Nb)) is an upper control limit corresponding to a confidence interval based on a Poisson distribution.

6. A method (1) for checking the thickness conformity of a radome according to any one of claims 1 to 5, wherein the representative data (r'k>i, dk>i) of conforming radome thicknesses correspond at least in part to representative measurements (rk>i) of conforming radome thicknesses, and wherein the step of defining a statistical model of a conforming radome includes a step of removing (104), from these representative measurements (rk>i) of conforming radome thicknesses, at least one outlier measurement associated with a point in a radar field of view of a conforming radome corresponding to a first point of the first mesh (Ml), and replacing said at least one outlier measurement with an average (r'k>i) of non-outliers measurements associated with points in the radar fields of view of other conforming radomes corresponding to said first point of the first mesh (Ml).

7. Method for checking (1) the thickness conformity of a radome according to claim 6, wherein the removal of said at least one outlier measurement uses a Grubbs test.

8. Method of checking (1) the thickness conformity of a radome according to any one of claims 1 to 7, wherein the representative data (r'k>i, dk>i) of thicknesses of conforming radomes correspond partly to measurements or averages of measurements (r'k>i) representative of thicknesses of N conforming radomes, and partly to virtual measurements associated with L virtual radomes determined during a determination step (108) of these L virtual radomes, comprising a generation phase (180) by random draws representative of normal laws, of series of values ​​each associated with a first distinct point of the first mesh (Ml).

9. A method for checking (1) the thickness conformity of a radome according to claim 8, wherein the normal laws being centered and reduced, the generation phase (180) is followed by a correction phase (182) of the series of values ​​generated by random draws such that an average of each series of corrected values ​​(X”ki>p) associated with each first point of a first surface control window (Fk) of the first mesh (Ml) is zero, and that the covariances between the series of corrected values ​​each associated with a first distinct point of the first surface control window (Fk) form the identity matrix.

10. Method (1) for checking the thickness conformity of a radome according to claim 9, wherein the correction phase (182) is followed by a generation phase (184) of virtual measurements (Yki>p) each associated with one of the L virtual radomes, wherein the virtual measurements (Yki>p) associated with the same point of the L virtual radomes corresponding to a first point of the first surface control window (Fk), are obtained as a function of the series of corrected values ​​(X”ki>p) associated with this first point, of an average of measurements (r'k>i) representative of the thicknesses of the N conforming radomes corresponding to this first point, and of the covariances between the series of data representative of the thicknesses of the N conforming radomes each associated with a first distinct point of the first surface control window (Fk).