{"id":1022,"date":"2023-10-24T09:12:28","date_gmt":"2023-10-24T09:12:28","guid":{"rendered":"https:\/\/uneedtalk.com\/config\/tredword\/?page_id=1022"},"modified":"2023-12-01T13:12:41","modified_gmt":"2023-12-01T13:12:41","slug":"inversion-method-for-a-generic-model-of-diffusely-reflecting-opaque-surfaces","status":"publish","type":"page","link":"https:\/\/radi-cal.org\/method\/inversion-method-for-a-generic-model-of-diffusely-reflecting-opaque-surfaces\/","title":{"rendered":"Inversion method for a generic model of diffusely reflecting, opaque surfaces"},"content":{"rendered":"<div class=\"row\"  id=\"row-604656108\">\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<h1>Inversion method for a generic model of diffusely reflecting, opaque surfaces<\/h1>\n<h2>4.13.1. Background and theory<\/h2>\n<p>The following method is proposed to derive generic surface reflectance models for painted or coated opaque surfaces. The approach can be considered an intermediate model downgrade and follows the principle of fallback approaches (see section 2.2). While simplified assumptions are made for the complex-valued refractive index, the diffuse reflection based on subsurface scattering still contains spectrally resolved information. Beyond that, even the assumption of a constant refractive index is still able to model scattering events in a spectrally resolved way.<br \/>Highly resolved spectral measurement data is often unavailable for practically relevant, coated or painted materials, e.g. lammelas of shading elements. Nonetheless, the manufacturer will usually still provide integrated reflectance values for specific ranges of the global radiation spectrum. The reflectance values usually represent either the total diffuse or the total hemispherical reflectance, i.e. they exclude or include specular reflection. The values are usually determined by using incidence light with a low, near-normal direction. The commonly applied standard ASTM E 903 (ASTM E 903, 1996) advises using incidence angles in the 6 to 12 degrees range. In order to determine the solar reflectance value, the reflectance values are weighted based on a global radiation standard (see section 4.4).<br \/>As discussed in section 4.3.3, the diffuse reflection of painted or coated material is commonly the result of subsurface scattering processes on pigments submerged in a binder. It is, therefore, a reasonable approximation to apply a generic model for these materials, assuming a constant or smooth refractive index for the binder with little or no absorption. The subsurface scattering on the pigments, providing the characteristic absorption and colour of the surface, is then only modelled by using the subsurface reflectance function \ud835\udf0c\ud835\udc60\ud835\udc60, see section 4.11.<br \/>The actual shape of the spectral reflectance function \ud835\udf0c\ud835\udc60\ud835\udc60 is determined in an optimisation process. Again, simulated annealing (section 5.2) is utilized for this purpose. In the optimisation process, the reflectance measurement is performed virtually and in accordance with the relevant standard (ASTM E 903, 1996). The spectral distribution of the incidence light is modelled according to the standard the measurement data refers to. The intensity and spectral distribution of the reflected rays are tracked, and the function \ud835\udf0c\ud835\udc60\ud835\udc60 is altered until the quantities measured in the model match the provided empirical data. The optimisation is performed by moving the knots of the implemented cubic-spline subsurface reflectance function (see section 4.11). The use of the spline functions to approximate the reflectance spectra takes into account that solid-state absorption (unlike gas absorption) usually exhibits a continuous and smooth spectral dependence.<br \/>The empirical reflectance data can be provided for hemispherical reflectance or for diffuse reflectance. In the latter case, a filter is used to exclude the specular reflection (gloss) from the measurement. The same features are implemented in the method to ensure that the model measurement mimics the real-world measurement.<br \/>As the empirical data commonly only covers four different reflectance values (UV, visible, NIR and total) and the cubic spline has comparably many degrees of freedom, an additional optimisation objective is generated by setting a target for the colour of the surface. For this purpose, the spectrum of all diffusely reflected rays is converted into an RGB triple. This is performed by first converting the wavelengths into the XYZ colour space of CIE 1931 and then further into the RGB colour space Light, sun and optics &#8211; applied principles, models and methods RadiCal, D. R\u00fcdisser 106 (Fairman et al., 1997; Service, 2016; Wikipedia, 2022). Since the brightness of the colour is already determined by the visible reflectance value, only the nomalised colour value is defined as the optimisation target. A third, less significant optimisation target measures the variability of the reflectance functions. It is simply defined as the average vertical distance between successive spline knots.<br \/>The three optimisation targets and an additional precision term are formulated as loss functions to assess the quality of the fitting.<br \/>They are defined as:<\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_703847547\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"672\" height=\"353\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a18-5.jpg\" class=\"attachment-large size-large\" alt=\"\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_703847547 {\n  width: 69%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<p>The weighted loss functions determine the entire cost function for the optimisation process:<\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_2036973272\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"592\" height=\"37\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a19-4.jpg\" class=\"attachment-large size-large\" alt=\"\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_2036973272 {\n  width: 70%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<p>\ud835\udc45_\ud835\udc52\ud835\udc5f\ud835\udc5f_\ud835\udc5f\ud835\udc5a\ud835\udc60\ud835\udc5e represents the deviation for the reflectance in the given spectral ranges, \ud835\udc50\ud835\udc5c\ud835\udc59_\ud835\udc52\ud835\udc5f\ud835\udc5f is the  deviation of the nomalised colour value and \ud835\udc63\ud835\udc4e\ud835\udc5f\ud835\udc56_\ud835\udc5d\ud835\udc52\ud835\udc5b\ud835\udc4e\ud835\udc59\ud835\udc61\ud835\udc66 provides a penalty for strongly fluctuating  functions. Finally, \ud835\udc5d\ud835\udc5f\ud835\udc52\ud835\udc50_\ud835\udc5d\ud835\udc52\ud835\udc5b\ud835\udc4e\ud835\udc59\ud835\udc61\ud835\udc66 adds a penalty related to the precision of the current solution. In  order to increase the performance of the optimisation, the first optimisation cycles are performed  with fewer samples. The potential statistical error of the results is considered by adding three times  the standard error to the cost function. This very effectively prevents early solutions from being  identified as optimal solutions. The coefficients \ud835\udc50\ud835\udc56 are used to balance the four components of the  cost function. In the optimisations performed, the applied cost function weighting was: \ud835\udc501 = \ud835\udc502 = \ud835\udc504 \u2219  = 1 and \ud835\udc503 = 1  1000 \u2044 .<br \/>\nThe method has proven to provide stable results matching the targeted reflectance quantities with  good accuracy. However, given the limited restrictions for the optimisation process and the relatively  high degree of freedom, the simulation results can certainly not be considered unique. In general,  multiple different solutions will be able to meet the optimisation targets. Even if the generated  subsurface reflection includes a certain degree of randomness and deviates from the actual spectral  reflectance, it will still be a good approximation for the unknown function as it exhibits the same  integrated reflectance values within the specific, provided spectral ranges and is based on a continuous reflectance function. In commonly applied alternative methods, a constant diffuse reflection coefficient \ud835\udf0c\ud835\udc60\ud835\udc5c\ud835\udc59, which is integrated over the entire global radiation range, is generally used.<\/p>\n<h2>4.13.2. Implementation<\/h2>\n<p>The implementation is contained in the module rc_TsubsurfaceOpt. The method OPT_ASTM_E903_Reflectance can be called to perform an optimisation of subsurface reflectance function represented by a cubic spline with N Knots.<\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_853010505\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"1001\" height=\"407\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a20-4.jpg\" class=\"attachment-large size-large\" alt=\"\" srcset=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a20-4.jpg 1001w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a20-4-800x325.jpg 800w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a20-4-768x312.jpg 768w\" sizes=\"(max-width: 1001px) 100vw, 1001px\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_853010505 {\n  width: 93%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<p>As can be seen in the implementation, currently, an optimisation method performing virtual measurements according to the ASTM E 903 standard is implemented. In order to align the model with the empirical data, the optimisation method offers the same functionality as the \u201creal-world measurement\u201d method. The incidence direction, spectrum, size of the gloss filter and type of reflectance value can be defined freely. Beyond that, the colour target, the number of free knots for the spline optimisation and the polarisation state of the incidence beam have to be provided.<\/p>\n<h2>4.13.3. Testing and validation<\/h2>\n<p>The feasibility of the inversion method is demonstrated based on three different surface coatings for Venetian blinds. The materials reflect the materials of the actual testing specimens used in the full system validation (chapter 8).<\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_767799797\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"888\" height=\"387\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a21-2.jpg\" class=\"attachment-large size-large\" alt=\"\" srcset=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a21-2.jpg 888w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a21-2-800x349.jpg 800w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a21-2-768x335.jpg 768w\" sizes=\"(max-width: 888px) 100vw, 888px\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_767799797 {\n  width: 75%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<p>Figure 67 shows the material data as provided by the manufacturer. The table provides the integrated reflectance values for the total, UV, visible and NIR ranges for three different coatings\/colours. The manufacturer also specifies the testing method and spectrum used. In the product data sheet, the Light, sun and optics &#8211; applied principles, models and methods RadiCal, D. R\u00fcdisser 108 material of the coating is referred to as PE (polyester). The specifications are inconsistent, as PE is usually used for polyethylene and not polyester. As most polymers exhibit refractive indices in the range of 1.4 to 1.6 in the relevant spectral range, a constant value of 1.5 is commonly assumed to model polymers generically. Alternatively, additional information on the gloss at specified angles could be used to approximate the material&#8217;s refractive index. A precise measurement of the gloss value can be used to narrow down the refractive index range for the applied spectral range, as the gloss value is a measure for the specular reflection (see sections 4.3.2 and 4.3.3).<br \/>\nThe generic value of 1.5 is assumed for the refractive index for the inversions performed here. The absorption of the material representing the binder of the coating is considered to be negligible; therefore, the complex-valued refractive index is defined as \ud835\udc5b\u0303 = 1.5 \u2212 0\ud835\udc56.<br \/>\nIn order to model the diffuse subsurface reflection caused by the pigments, the subsurface reflectance functions \ud835\udf0c\ud835\udc60\ud835\udc60(\uf06c) are determined in the optimisation process described above. The results are depicted in Figure 68. As shown in Table 9, the integrated reflectance values for the different spectral ranges of optimised subsurface functions match the targeted values with reasonable accuracy. The deviations are in the range of a few tenths of a per cent. The reason for the deviation is that the optimisation based on smooth spline functions additionally has to consider the colour target. Without this objective, the optimisation would be able to provide almost perfectly matching reflectance values. Without additional color measurements available, the nomalised colour target was set to \ud835\udc5f\ud835\udc52\ud835\udc51 = \ud835\udc54\ud835\udc5f\ud835\udc52\ud835\udc52\ud835\udc5b = \ud835\udc4f\ud835\udc59\ud835\udc62\ud835\udc52 = 1.0 for all three materials, representing a desaturated colour (white or grey). <\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_2146997357\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"899\" height=\"500\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a22-3.jpg\" class=\"attachment-large size-large\" alt=\"\" srcset=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a22-3.jpg 899w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a22-3-719x400.jpg 719w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a22-3-768x427.jpg 768w\" sizes=\"(max-width: 899px) 100vw, 899px\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_2146997357 {\n  width: 90%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner\"  >\n<p>The subsurface reflectance function for the white coating shows strong fluctuations in the UV range. This is required to fulfil the steep drop from 80.9% reflectance in the visible to only 12.2% in the UV range (stated in the manufacturer&#8217;s datasheet). In order to assess the quality and relevance of the optimised functions, it has to be considered that the visible and lower parts of the NIR spectral range are most relevant. The UV range contributes only 3%, and the range above 1460 nm accounts for only 10% of global radiation (see Table 2).<\/p>\n<\/div><\/div>\n<div class=\"col small-12 large-12\"  ><div class=\"col-inner text-center\"  >\n\t<div class=\"img has-hover x md-x lg-x y md-y lg-y\" id=\"image_1868737727\">\n\t\t\t\t\t\t\t\t<div class=\"img-inner dark\" >\n\t\t\t<img loading=\"lazy\" decoding=\"async\" width=\"824\" height=\"422\" src=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a23-3.jpg\" class=\"attachment-large size-large\" alt=\"\" srcset=\"https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a23-3.jpg 824w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a23-3-781x400.jpg 781w, https:\/\/radi-cal.org\/method\/wp-content\/uploads\/2023\/10\/a23-3-768x393.jpg 768w\" sizes=\"(max-width: 824px) 100vw, 824px\" \/>\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n\t\t\t\t\t\t\t\t\n<style scope=\"scope\">\n\n#image_1868737727 {\n  width: 95%;\n}\n<\/style>\n\t<\/div>\n\t\n<\/div><\/div>\n<\/div>\n<style>\n#menu-main li:nth-child(4) .sub-menu{display:block !important;}\n<\/style>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"page-right-sidebar.php","meta":{"footnotes":""},"_links":{"self":[{"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/pages\/1022"}],"collection":[{"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/comments?post=1022"}],"version-history":[{"count":16,"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/pages\/1022\/revisions"}],"predecessor-version":[{"id":1967,"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/pages\/1022\/revisions\/1967"}],"wp:attachment":[{"href":"https:\/\/radi-cal.org\/method\/wp-json\/wp\/v2\/media?parent=1022"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}