|
I.INTRODUCTIONThe Solar Irradiance flux experiment a spectral filtering process as the irradiance flux propagates from the source to the sensor. When using radiometer sensors for optical measurement both the atmospheric transmittance, sensor’s filter and the responsivity of the detector produces a spectral filtering, see Fig. 1 The total response of the radiometer is the integral of the convolution from the spectral Irradiance source reaching to the sensor with the Optical Transfer Function (OTF) that is composed by the sensor filter and detector responsivity. For MetSiS and DREAM SIS radiometers the sensor filter are composed by Schott filter and an interference filter. Schott filters were added to the optical system to improve the behavior of interference filters. The OTF is dependent on the next parameters: temperature, angle of incidence of the incoming radiation and the incoming light distribution (diffusive or collimated). Any element of the optical system contributes uncertainty to the measurement, being the interference filter the most affecting in the optoelectronic chain by these parameters. Due the characterization of interference filter under different conditions of temperature and angular orientation is difficult, so modelling the behavior of the interference filters is helpful in order to assess how much uncertainty this optical element introduces into the final OTF of the radiometer instrument [1], [2]. We describe a method to find an equivalent interference filter based in the characteristic matrix approach, the theoretical model calculate filter properties as a function of the number of layers, the layers thickness and the type of material every layer is made of. Combining experimental measurement of optical transmittance with the theoretical model is possible to get an equivalent filter using inverse methods. II.EXPERIMENTAL SETUPThis part describes the experimental setup to perform the optical transmittance measurement with different orientation of the filter respect to the incident light. The setup is divided in four parts; see Fig.2:
III.MULTILAYER THIN FILM MODELINGThe Transfer Matrix method calculates the filter properties using the characteristic matrix approach, the method is presented in detail in most optical coating textbooks[3], [4],[5]. In the equations shown below n is the refractive index, n1 and ns are the refractive index from external material and substrate respectively, d is the layer thickness, φ is the angle of the incident light, ψ is the refraction angle in the layer and k0 is For s-polarization the jth maxtrix is represented by: and for p-polarization: The characteristic matrix describing the multilayer is: Table 1, show the relation between the multilayer matrix and the optical property of transmittance for the filter, the total transmittance of the filter will be the average between the transmittance for s and p polarizations. Table 1.Amplitude transmission and transmittance relations of a multilayer film.
Table 2.Parameters value of the cost function when the optimization was stopped.
IV.RETRIEVAL OF PHYSICAL PARAMETERSRetrieval of physical parameters is a problem of optimization or inverse problem where the model described in the previous section is feedback with the physical parameters to produce the best fit of the theoretical transmittance respect to the experimental one. To solve the optimization problem two algorithms has been used: (1) A heuristic optimization using a genetic algorithm (GA) performed by choosing the nth best samples of the complete offspring at every iteration, and (2) a local optimization based on the Levenberg–Marquardt algorithm is used to refine the best offspring in the nth iteration [6], [7]. The first one is based on the genetic processes of biological organism. They usually work with a population of individuals, each representing a possible solution of the problem. The best individuals, according with a fitness function, are selected to create new individual recombining and modifying the information stored in the old ones. The new individuals are then inserted in the population with some of the old ones generating a new offspring. This process continues until a stop condition is reached, and the best individuals found for the algorithm are considered the solution of the problem. Matlab software and toolbox has been used to implement this algorithm [8]. A.PROBLEM ENCODINGThe following terms are referred to any GA:
Some information about the physical parameters was provided by the manufacturer like the number of layer, the type of material used and the total physical thickness for the High & Low-Index-Layers. In this inverse problem the information from manufacturer was used, so the individual only has two chromosomes the layer thickness for material of High/Low refraction index.
The first term in (8) is the root-mean-square (RMS) error between the theoretical and the measured optical transmittance, the second and third terms are the difference between the physical thickness and the calculated thickness for the High/Low-Index-Layers respectively. α and β are constant used to weight these terms. B.PROCEDURE OF GA OPTIMIZATIONFor simplification, all calculations were performed with constant refractive index and at normal incident of light. The population number was chosen between 50-60 individuals. The population initialization can be set up randomly or it is possible to use different strategies [9], [10] to choose new individuals like suppose the individual consists on quarter wave layers (QW) or similar allowing a random variation around the reference wavelength (9). The next step are common to any simple GA: The number of generations N for the GA was establishing in 300-400. C.REFINE GA SOLUTIONThe Best Individual from GA during N generations is refined using a local search based on Levenberg–Marquardt algorithm, the cost function is the RMS error. A Maximum of 1000 iterations was established as stop condition for this method. D.STOP CONDITIONThe refine solution is introduced into the population of the GA and the process of optimization is continued until the next N generation. Then, the best individual is improved using Levenberg–Marquardt algorithm again. The process is stopped manually when the cost function value is below or near to 0.5; if the fitting solution is successful the optimization is finished, otherwise the optimization process continues. V.RESULTSFor these paper three filters has been optimized. The quality of the optimization has been successful as shown in the figures below: VI.CONCLUSIONS AND FUTURE PLANSCombining the Matrix formulation for multilayer filter, global and local search algorithm a reverse-engineering software procces has been developed and validated succesfully to find an equivalent filter from experimental measurement of interference filters using. We plan continue developing the algorithm to change the algorithm from a mono-objective cost function to a multi-objective optimization that include several cost function dependent on physical parameters. ReferencesReid E. Basher and W. Andrew Matthews,
“Problems in the use of interference filters for spectrophotometric determination of total ozone,”
Journal of Applied Meteorology,
(1977). Google Scholar
H. Guerrero,
“Development of miniaturized instrumentation for Planetary Exploration and its application to the Mars MetNet Precursor Mission,”
EGU, Viena, Austria
(2010). Google Scholar
O.Stenzel,
“The Physics of thin film optical spectra: An introduction,”
Springer-Verlag, Berlin Heidelberg New York
(2005). Google Scholar
H. Angus Macleod,
“Thin-Film Optical Filters,”
CRC Press Taylor & Francis Group2010). Google Scholar
Stéphane Larouche and Ludvik Martinu,
“OpenFilters: open-source software for the design, optimization, and synthesis of optical filters,”
APPLIED OPTICS, 47
(13),
(2008). Google Scholar
Wolfram MathWorld,
“Genetic Algorithm,”
http://marthworld.wolfram.com Google Scholar
Wolfram MathWorld,
“Levenberg-Marquardt Method,”
http://marthworld.wolfram.com Google Scholar
Matlab,
“Global Optimization Tool box,”
Google Scholar
Uwe B. Schallenberg,
“Antireflection design concepts with equivalent layer,”
APPLIED OPTICS, 45
(7),
(2006). Google Scholar
Ulrike Schulz, Uwe B. Schallenberg, and Norbert Kaiser,
“Antireflection coating design for plastic optics,”
41
(13),
(2002). Google Scholar
|