Abstract:
A method for quantitative measurement of surface accuracy of an area is provided. This comprises directing a monochromatic flat light wave towards a predefined surface area, recording an image of the reflected light with a camera and lens system focused on said surface area, and deducing surface accuracy parameters from the recorded image. The method is characterized in that said surface accuracy parameters are determined by obtaining a Fourier transform of the recorded image. Then, fitting predetermined Fourier components to a Fourier spectrum of said Fourier transform, wherein said Fourier components are determined along the major and across the minor elongation axes of the Fourier transform as at least a large Gaussian component, and a peak of the spectrum. Followed by, determining surface accuracy parameters of said surface area from said Fourier components.

Description:
TECHNICAL FIELD 
       [0001]    The embodiments herein relates to a method for in-situ optical quantitative measurement of surface accuracy of an area during polishing process. 
       BACKGROUND 
       [0002]    Surface modifications, especially polishing, are a required process in many industrial production areas. These include:
       Automotive, e.g. polishing of moulds for plastic parts such as boot lids and light covers and lenses.   Engine parts such as cam axes.   Medical implants, e.g. the surface of artificial hip joints.   Optics such as injection moulds for safety glasses and contact lenses.       
 
         [0007]    Such polished surfaces can be several square meters in size and require a surface accuracy of 1 micrometer or less. Today, the sample is moved from the polishing process to a laboratory where the surface is measured at small areas at a time in a time-consuming process. Furthermore, most of the polishing is done manually, where the manual polisher establishes the quality by looking at the tool/mould, whereas the automatic procedure using robot or laser has, in general, to be checked by taking the tool/mould out of the process chamber. 
         [0008]    In order to check the surface accuracy, several methods are known in the art. 
         [0009]    U.S. Pat. No. 6,084,671 discloses a method and apparatus for surface analysis using Gaussian beam profiles. A Gaussian beam is shaped with a lens system and illuminates a surface. The reflected light is distorted by surface irregularities. The reflected beam is then observed through a series of spatial filters which basically remove the original Gaussian beam shape. The measured beam with difference in residual beam shape is the recorded with a camera. This solution requires that the surface is illuminated by a very specific beam pattern in the form a Gaussian shape. Further, the illumination of the surface is uneven over the surface in that the edges of the studied area is less intensely illuminated than the center and any statistical measure of the surface structure is therefore biased towards the center. 
         [0010]    US 2005/0046870 disclose a method for characterizing surfaces by illuminating it at an angle to show shadowgraph of larger structure on the surface. This solution describes and discusses geometrical reflections only. It does not relate to illuminating polished surfaces with monochromatic light. Further, the illumination is not directed perpendicularly to the surface. 
         [0011]    US 2012/107971 discloses a polishing pad assembly for a chemical mechanical polishing apparatus. The polishing pad has a polishing surface and a surface opposite the polishing surface for attachment to a platen, and a solid light-transmissive window formed in the polishing pad. The light-transmissive window is more transmissive to light than the polishing pad. The light-transmissive window has a light-diffusing bottom surface. This solution requires a chemical transparent top layer over the surface to produce interference with the surface reflected wave. This is a single spot measurement instrument and technique which measures interference between two overlapping surfaces. 
         [0012]    U.S. Pat. No. 4,873,434 discloses a scanning optical microscope which comprises a source of optical radiation and means for focusing radiation from the source into an interrogating spot on a surface under examination. The spot is deflected about a point on the surface, and the surface topography is measured by measuring the amplitude and/or phase of the radiation reflected from the surface at the spot position. This technique measures one spot at a time. Each spot is handled separately and does not give any direct quantitative analysis of a larger area. 
         [0013]    U.S. Pat. No. 5,917,191 discloses a method for measuring surface topography characterized by making multiple scans of the surface with a laser scanning unit and utilizing the multiple scans to create representations of the surface&#39;s topography. The surface topography data can also be used to calculate the compressive or tensile stress caused by a thin film applied to the surface of a semiconductor wafer. The apparatus of the solution scans a laser beam across a surface in an x direction, and detects displacements of a reflected portion of the laser beam in a z direction. A pair of photo detectors is used to translate z direction displacements of the reflected beam into analog signals which are digitized and input into a microcomputer for analysis. The multiple scans of the surface are preferably accomplished by placing the work piece on a pedestal which can be rotated to various angular positions. This technique is scanning point by point to give point topographic data. It neither describes nor teaches quantitative analysis over a large area simultaneously. 
       SUMMARY 
       [0014]    It is therefore an object of the embodiments herein to provide a method for quantitative measurements of surface quality during polishing processes, wherein the method would eliminate the drawbacks of the prior art. 
         [0015]    The embodiments herein relates to a method and apparatus for determining quantitative accuracy of a surface, wherein a beam of parallel monochromatic light is sent perpendicularly towards the surface of a sample. Quantitative parameters defining the surface topography are then deduced by simultaneously recording the intensity of the reflected light at each pixel in a camera image focused with a camera lens on the surface area under investigation. The area intensity image is then Fourier-transformed from the surface aperture dimensions to scattering angle dimensions. Image surface parameters are then calculated from the intensity Fourier transform and coupled to quantitative surface parameters. 
         [0016]    The embodiments herein describes a method and apparatus for in-situ optical quantitative measurement of surface accuracy of an area during polishing process, wherein a monochromatic flat light wave is directed towards a predefined surface area and the image of the reflected light is recorded with a camera and lens system focused on the surface area. The flat light wave is directed perpendicularly or at an angle less than 6 degrees towards the surface. Further, the camera and lens system may be positioned perpendicularly to or at an angle less than 6 degrees towards the surface. The difference in angle between the reflected light entering the camera and the optical axis of the camera and lens system may be less than 6 degrees. 
         [0017]    The surface accuracy parameters are then deduced from the recorded image. The surface accuracy parameters may be determined with a Fourier transform of the recorded image, predetermined Fourier components are fitted to the Fourier spectrum of the surface area, and surface accuracy parameters are determined from the Fourier components. Such Fourier components may then be determined along the major and across the minor elongation axes of the Fourier transform as a large Gaussian component; the peak of the spectrum; a wide Gaussian component as foundation to the peak; and component peak outside the central maximum. 
         [0018]    Further components may be fitted to the major and minor axes as the quota of the major and minor axes of large Gaussian components; the direction of the major axis; and the direction of the component peak outside the maximum centre. Surface accuracy parameters may then be deduced such as the large Gaussian component along the major axis divided with the central peak maximum, thus defining the quality of the polishing; the wide Gaussian component near the peak divided by the large Gaussian component defining the gloss of the surface; the large Gaussian component along the major elongation axis divided with the large Gaussian component along the minor axis defining the symmetry of the surface; the direction of the major axis showing the direction of structural lines of the surface, and the size of a component peak outside the centre maximum peak indicating structural lines. 
         [0019]    The method and apparatus are defined in number of wavelengths rather than geometrical size, and are not limited to a specific optical wavelength even though this is used as example. 
         [0020]    The method and apparatus have been verified and tested for the cases where the focus pixel area is less than 100 wavelengths in diameter, and the surface observing imaging area is larger than 4000 wavelengths in diameter. These are not theoretical limitations but rather practical due to present limits on number of pixels of camera. 
         [0021]    Embodiments herein are especially useful for determining surface parameters over a large area. It is robust and mounted to make measurements inside a process chamber. The working distance is long, 70-800 mm, which is much longer than for microscope objective systems used presently. 
         [0022]    The industrial use ranges from decision on process route and method in automatic polishing for automatic, medical and optical industrial tools, moulds and products; to objective decision making in manual polishing processes. The apparatus can be mounted on a robot arm, outside a laser process chamber or on the floor in manual polishing workshops. The method comprises taking a single image of a large part, in the case discussed below 4×4 mm, of the surface area and is therefore quick and much more insensitive to vibrations compared to single point measurements and interferometry instruments. The total surface area can then also be measured by stepping over the total area in steps of, in the case discussed below, 2-4 mm. 
         [0023]    The surface accuracy may be defined as how the surface is topographically deviating from being flat in the same plane as the light wave. 
         [0024]    Some specific advantages of some embodiments herein are:
       Long working distance: can be placed 70-800 mm from the target area;   Large surface simultaneously: fast measurement;   Quantitative data over large area as feedback to the polishing process.       
 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0028]    Embodiments will now be described by way of example with reference to the accompanying drawings, of which: 
           [0029]      FIG. 1  illustrates an apparatus according to embodiments herein; 
           [0030]      FIG. 2  illustrates a geometry of the system according to embodiments herein; 
           [0031]      FIG. 3  illustrates the intensity of the Fourier transform of the image. 
       
    
    
     DETAILED DESCRIPTION 
       [0032]      FIG. 1  depicts a monochromatic light source (laser or laser diode)  1  at wavelength λ emits light into a beam expander  2 . It may here be noted that a beam expander expands the laser beam so that the central area where the measurements are taken can be considered as having the same phase and amplitude on a flat plane, i.e. there is a monochromatic flat light wave hitting the surface. The flat wave front light from the beam expander is reflected from a beam splitter  3  to an object surface  4 . The light is reflected from the surface through the beam splitter into a focusing lens  5  and recorded with a camera  6 . The camera used in the example is a CCD camera; having a camera lens of 75 mm. The diameter of the focusing spot at the surface in the example is 2 micrometer. 
         [0033]    The geometry of the system is shown in  FIG. 2 . The object surface  4  is limited to inside of the focusing area spot. The surface has a reference plane  7  and the incoming flat wave plane is represented with a line  8 . The normal  9  to the surface reference plane  7  has a tilt angle φ to the normal  10  of the incoming flat wave  8 . The incoming light from direction  11  at x=x 1 , hits the object surface  4  at point  12  x 1 , z 1 +Δz 1 . The phase of the light  13  emitted from point  12  in a direction θ to the normal  10  of the incoming light wave plane  8  is calculated. The phase center is a position (x=0, z=0) in the coordinate system defined by the incoming flat wave  8  and the normal  10  at a position where it crosses the surface reference plane  7 . 
         [0034]      FIG. 3  shows the intensity of the Fourier transform of the image in a cut through the major axis of the elliptical Gaussian noise distribution of the transform. The notations a 1 -a 4  are components discussed below. 
         [0035]    The phase of the incoming light at z=0 is zero at all points x along the object surface. The discussion below is concentrated to the light scatter inside one focusing spot, i.e. same as one recorded pixel on the CCD. We also assume that the emitting source at position (x 1 , z 1 +Δz 1 ) is much smaller than one wavelength and thus works as a spherical scatter point of the incoming wave. 
         [0036]    The light coming in through x=x 1  and emitted at the angle θ is delayed by: 
         [0000]      Δ L=z   1   +Δz   1   +Δl  
 
         [0037]    Where Δz 1  is the offset of the surface from a plane and 
         [0000]    
       
         
           
             
               z 
               1 
             
             = 
             
               
                 x 
                 1 
               
               × 
               
                 tan 
                  
                 
                   ( 
                   ϕ 
                   ) 
                 
               
             
           
         
       
       
         
           
             
               Δ 
                
               
                   
               
                
               l 
             
             = 
             
               
                 
                   z 
                   1 
                 
                 + 
                 
                   Δ 
                    
                   
                       
                   
                    
                   z 
                 
               
               
                 cos 
                  
                 
                   ( 
                   Θ 
                   ) 
                 
               
             
           
         
       
     
         [0038]    In addition, there is a geometrical path difference Ag, between the wave emitted from position (x=0, z=0) and the wave emitted from the surface and passing z=0 at x=x 1 +Δx 1 . This path can be calculated as 
         [0000]      Δ g =( x   1   +Δx   1 )×sin(Θ)
 
         [0000]      Where 
         [0000]      Δ x   1 =( z   1   +Δz   1 )×tan(Θ)
 
         [0039]    The phase difference between the wave coming in through and emitted at the phase reference position (x=0, z=0) and the corresponding wave at (x=x 1 , z=z 1 +Δz 1 ) at the line going through position (x=0, z=0) perpendicular to the emitting angle θ is therefore 
         [0000]    
       
         
           
             
               δ 
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               2 
                
               π 
                
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       L 
                        
                       
                         ( 
                         x 
                         ) 
                       
                     
                   
                   - 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       g 
                        
                       
                         ( 
                         x 
                         ) 
                       
                     
                   
                 
                 λ 
               
             
           
         
       
     
         [0040]    This function can be further extended by introducing a curvature of the focal surface f(x) to become 
         [0000]    
       
         
           
             
               δ 
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               2 
                
               π 
                
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       L 
                        
                       
                         ( 
                         x 
                         ) 
                       
                     
                   
                   - 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       g 
                        
                       
                         ( 
                         x 
                         ) 
                       
                     
                   
                   + 
                   
                     f 
                      
                     
                       ( 
                       x 
                       ) 
                     
                   
                 
                 λ 
               
             
           
         
       
     
         [0041]    The complex voltage of the signal reflected at position (x=x 1 , z=z 1 +Δz 1 ) is 
         [0000]        V   1 ( x   1 )= A   1   e   −iδ     1     (x     1     )    
         [0000]    where A 1  is the amplitude of the signal. 
         [0042]    This is readily extended to two dimensions (x, y). The focal area is tilted at angles (φ x , φ y ) and the offset from a plane is the two dimensional function Δz(x, y) The complex voltage of the reflected light from the focusing area in the direction (θ x , θ y ) is then the coherent sum of all complex voltages as: 
         [0000]        V (focalarea)=Σ x Σ y   V ( x,y )
 
         [0043]    The intensity emitted in direction (θ x , θ y ) from the focal area S(X, Y) is 
         [0000]        I   X,Y (α,β)=           V ( X,Y )* V ( X,Y )*           T  
 
         [0000]    T is the averaging time for the intensity, X, Y is the position of the focal area on focal plane over the sample, α is the zenith angle and β is the azimuth angle of the emitted light. 
         [0044]    In the same manner we now define the tilt angles φ x  and φ y  as zenith X and azimuth ξ angles to the normal of the focal plane. We also set the focal plane to be parallel to the incoming wave front. 
         [0045]    According to some embodiments herein, the camera is angled to look onto the sample in the same direction as the light from the beam expander, as shown in  FIG. 1 . The zenith angle is then α=0, or (θ x =0, θ y =0). In this case the phase offset for the position (x, y) in the plane where the x-axis in the direction of the azimuth angle φ y  can be written as: 
         [0000]    
       
         
           
             
               δ 
                
               
                 ( 
                 
                   x 
                   , 
                   y 
                 
                 ) 
               
             
             = 
             
               2 
                
               π 
                
               
                 
                   
                     x 
                     * 
                     
                       tan 
                        
                       
                         ( 
                         
                           ϕ 
                           x 
                         
                         ) 
                       
                     
                   
                   + 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       z 
                        
                       
                         ( 
                         
                           x 
                           , 
                           y 
                         
                         ) 
                       
                     
                   
                   + 
                   
                     f 
                      
                     
                       ( 
                       
                         x 
                         , 
                         y 
                       
                       ) 
                     
                   
                 
                 λ 
               
             
           
         
       
     
         [0046]    The complex voltage emitted from the focal area is then 
         [0000]    
       
         
           
             
               V 
                
               
                 ( 
                 focalarea 
                 ) 
               
             
             = 
             
               
                 Σ 
                 x 
               
                
               
                 Σ 
                 y 
               
                
               
                 A 
                  
                 
                   ( 
                   
                     x 
                     , 
                     y 
                   
                   ) 
                 
               
                
               
                  
                 
                   
                     - 
                     i 
                   
                    
                   
                       
                   
                    
                   2 
                    
                   π 
                    
                   
                     
                       
                         x 
                         * 
                         tan 
                          
                         
                             
                         
                          
                         
                           ϕ 
                           x 
                         
                       
                       + 
                       
                         Δ 
                          
                         
                             
                         
                          
                         
                           z 
                            
                           
                             ( 
                             
                               x 
                               , 
                               y 
                             
                             ) 
                           
                         
                       
                       + 
                       
                         f 
                          
                         
                           ( 
                           
                             x 
                             , 
                             y 
                           
                           ) 
                         
                       
                     
                     λ 
                   
                 
               
             
           
         
       
     
         [0047]    Assuming that the material is the same at each pixel x, y, and the surface variations Δz and f are small over the focal area, then A is constant over the focal area and this becomes: 
         [0000]    
       
         
           
             
               V 
                
               
                 ( 
                 focalarea 
                 ) 
               
             
             = 
             
               A 
               * 
               
                 Σ 
                 x 
               
                
               
                 Σ 
                 y 
               
                
               
                  
                 
                   
                     - 
                     i 
                   
                    
                   
                       
                   
                    
                   2 
                    
                   π 
                    
                   
                     
                       
                         x 
                         * 
                         tan 
                          
                         
                             
                         
                          
                         
                           ϕ 
                           x 
                         
                       
                       + 
                       
                         Δ 
                          
                         
                             
                         
                          
                         
                           z 
                            
                           
                             ( 
                             
                               x 
                               , 
                               y 
                             
                             ) 
                           
                         
                       
                       + 
                       
                         f 
                          
                         
                           ( 
                           
                             x 
                             , 
                             y 
                           
                           ) 
                         
                       
                     
                     λ 
                   
                 
               
             
           
         
       
     
         [0048]    The intensity emitted from the focal area at position X, Y can then be written as 
         [0000]    
       
         
           
             
               I 
               
                 X 
                 , 
                 Y 
               
             
             = 
             
               
                 〈 
                 
                   
                     A 
                     2 
                   
                    
                   
                     Σ 
                     x 
                   
                    
                   
                     Σ 
                     y 
                   
                    
                   
                      
                     
                       
                         - 
                         i 
                       
                        
                       
                           
                       
                        
                       2 
                        
                       π 
                        
                       
                         
                           
                             x 
                             * 
                             tan 
                              
                             
                                 
                             
                              
                             
                               ϕ 
                               x 
                             
                           
                           + 
                           
                             Δ 
                              
                             
                                 
                             
                              
                             
                               z 
                                
                               
                                 ( 
                                 
                                   x 
                                   , 
                                   y 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             f 
                              
                             
                               ( 
                               
                                 x 
                                 , 
                                 y 
                               
                               ) 
                             
                           
                         
                         λ 
                       
                     
                   
                   * 
                   
                     Σ 
                     x 
                   
                    
                   
                     Σ 
                     y 
                   
                    
                   
                      
                     
                       i 
                        
                       
                           
                       
                        
                       2 
                        
                       π 
                        
                       
                         
                           
                             x 
                             * 
                             tan 
                              
                             
                                 
                             
                              
                             
                               ϕ 
                               x 
                             
                           
                           + 
                           
                             Δ 
                              
                             
                                 
                             
                              
                             
                               z 
                                
                               
                                 ( 
                                 
                                   x 
                                   , 
                                   y 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             f 
                              
                             
                               ( 
                               
                                 x 
                                 , 
                                 y 
                               
                               ) 
                             
                           
                         
                         λ 
                       
                     
                   
                 
                 〉 
               
               T 
             
           
         
       
     
         [0049]    Here A 2  is the intensity reflected by a perfect mirror and the second term is the attenuation of this by out of phase interference depending on the tilt angle φ x , the surface variation Δz(x, y) and the form factor f(x, y). Note that the interference also can be constructive resulting in Speckle radiation from the focal surface area. 
         [0050]    Assuming that the material is the same over the whole sample and that the amplitude A thus can be considered the same everywhere, the direct returned intensity for each focal area is only dependent on the local surface geometry:
       1. tilt angle φ x ,   2. surface variation Δz(x, y), and   3. form factor f(x, y).       
 
         [0054]    The focal area is the same as the Airy disc of the camera objective. The camera lens system is constructed so that an Airy disc is imaged onto one CCD pixel. In this case therefore the observed intensity image of the sample is a direct measure of the relative local surface geometry of the focal areas of the sample. Nearby pixels with the same intensity form image structure, caused by e.g. a similar tilt angle. 
         [0055]    The embodiments herein are not limited to zenith angle of zero degrees. A small angle between the incoming and reflected beams is possible if this angle is small compared to the diffraction angle of size scale under consideration. In practice, this means that the difference in angle between the illuminating monochromatic light and the camera lens system optical axis should be less than 6 degrees. 
       Example 
       [0056]    A tool surface is being polished by a machine or a human. The polishing process requires measurements of the surface area representing the surface roughness in order to be able to determine the way to proceed with the polishing process in order to reach a predefined smoothness. In this case, the target in blocks of 4×4 mm with 2000×2000 pixels at focal resolution of 2000 nm is studied. The light is emitted from a laser diode at the wavelength of 635 nm. 
         [0057]    An image of the 4×4 mm surface area on the target surface is recorded with a CCD camera with 2000×2000 gray scale pixels. The image in position (i, j) at the surface is I i,j (X, Y) is Fourier transformed and the absolute value of the Fourier transform is stored as dB in a matrix: 
         [0000]        S   i,j ( u,v )=20*log 10 (∥ F ( I ( X,Y )∥)
 
         [0058]    Here, u and v are the spatial frequency in the X and Y directions, respectively. 
         [0059]    S i,j (u, v) is a two-dimensional spatial spectrum of the surface area geometry. This is now studied and the following surface structure parameters are extracted as:
       1. The central peak (S 1 ) at (u=0, v=0) is the total reflected intensity of the area.   2. There is a large bulge (S 3 ), which corresponds to the randomly distributed intensity variations.   3. The peak S 1  has a pedestal (S 2 ), which corresponds to the nearby pixels having similar intensity.   4. Isolated peak (S 4 ) correspond to a large-scale linear structure.   5. S 3  is elliptical (asymmetry), indicating that the intensity is more ordered in one direction along the structure angle φ.       
 
         [0065]    We can now connect these intensity structure parameters to the surface geometry parameters above.
       A perfect mirror with tilt angle φ x =0, Δz=0, and f=0 everywhere would have a single peak S 1  only.   Linear lines, i.e. caused by the previous honing step, would show as isolated peaks S 4 .   A one-directional polishing would show more random features in one direction compared to the elliptical distributed noise S 3  with major axis position angle φ gauss .   A rough surface would have a large variation in tilt angles φ x  and show a large Gaussian distribution S 3 .   A surface where the tilt angles φ x  are similar for most focal areas would show smaller S 3  and larger S 2 . This would look like a glossy surface, but not quite a mirror.   Also, the gloss can be more in one direction and more mirror-like in the other at an angle φ gloss .   Finally, when the focal area geometry becomes flat, i.e. both Δz and f are small (&lt;λ/10), the surface would turn into a mirror.       
 
         [0073]    The surface parameters are then further reduced as cuts are made through the Fourier transform S(u,v) along the major axis position angle φ gauss  and perpendicular to this. Gaussian profiles are fitted to the features described above and denoted as: 
         [0074]    major axis cut: a 1   p , a 2   p , a 3   p , and a 4   p  ( FIG. 3 ); 
         [0075]    minor axis cut: a 1   m , a 2   m , a 3   m.    
         [0076]    The maximum intensity in the center of the transform S is denoted as amax. Note that these values are in dB. The structural parameters are then calculated as:
       1. amax−a 3   p =Gaussian noise part, i.e. the amount of wide band noise compared to the total returned intensity.   2. (amax−a 3   p )/(amax−a 3   m )=the asymmetry of the Gaussian noise.   3. φ gauss =the angle of the asymmetry structure.   4. a 4   p −a 3   p =the amount of structural lines compared to Gaussian noise.   5. a 2   p −a 3   p =gloss.       
 
         [0082]    Then, these parameters can either be mapped over the sample surface, or used individually for 4×4 mm spot tests of the sample. 
         [0083]    The parameters can then be fed back to the polishing process to determine subsequent polishing steps. 
         [0084]    To perform the method actions in the apparatus for determining quantitative accuracy of a surface, the apparatus may in some embodiments comprise or be connectable to a processing circuitry, which may also be referred to as processing module or processor. The processing circuitry may be configured to perform the image processing and analysis steps of the method actions as described in the embodiments herein. 
         [0085]    The embodiments for determining quantitative accuracy of a surface may be implemented through one or more processors, such as, e.g. the processing circuitry comprised in or connectable to the apparatus depicted in  FIG. 1 , together with computer program code for performing the functions and actions of the embodiments herein. The program code mentioned above may also be provided as a computer program product, for instance in the form of a data carrier carrying computer program code or code means for performing the embodiments herein when being loaded into the processing circuitry. The computer program code may e.g. be provided as pure program code in the processing circuitry or on a server and downloaded to the processing circuitry. The carrier may be one of an electronic signal, optical signal, radio signal, or computer readable storage medium, such as, e.g. electronic memories like a RAM, a ROM, a Flash memory, a magnetic tape, a CD-ROM, a DVD, a Blueray disc, etc. 
         [0086]    Thus, the apparatus may further comprise or be connectable a memory, which may be referred to or comprise one or more memory modules or units. The memory may be arranged to be used to store executable instructions and data to perform the methods described herein when being executed in the processing circuitry. Those skilled in the art will also appreciate that the processing circuitry and the memory described above may refer to a combination of analog and digital circuits, and/or one or more processors configured with software and/or firmware, e.g. stored in the memory, that when executed by the one or more processors such as the processing circuitry perform the method as described above. One or more of these processors, as well as the other digital hardware, may be included in a single application-specific integrated circuit (ASIC), or several processors and various digital hardware may be distributed among several separate components, whether individually packaged or assembled into a system-on-a-chip (SoC). 
         [0087]    Thus, a computer program, comprising instructions which, when executed on at least one processor, e.g. the processing circuitry described above, cause the at least one processor to carry out the method for determining quantitative accuracy of a surface as described above is presented. Also, a carrier containing the computer program, wherein the carrier is one of an electronic signal, optical signal, radio signal, or computer readable storage medium, is presented. 
         [0088]    The embodiments herein are not limited to the above described preferred embodiments. Various alternatives, modifications and equivalents may be used. Therefore, the above embodiments should not be construed as limiting.