Translator Disclaimer
1 October 2004 Scaled topometry in a multisensor approach
Author Affiliations +
Reliable real-time surface inspection of extended surfaces with high resolution is needed in several industrial applications. With respect to an efficient application to extended technical components such as aircraft or automotive parts, the inspection system has to perform a robust measurement with a ratio between depth resolution and lateral extension of less than 10–6. This ratio is at least 1 order beyond the solutions that are offered by existing technologies. The concept of scaled topometry consists of a systematic combination of different optical measurement techniques with overlapping ranges of resolution systematically to receive characteristic surface information with the required accuracy. In such a surface inspection system, an active algorithm combines measurements on several scales of resolution and distinguishes between local fault-indicating structures with different extensions and global geometric properties. The first part of this active algorithm finds indications of critical surface areas in the data of every measurement and separates them into different categories. The second part analyzes the detected structures in the data with respect to their resolution, and decides whether a further local measurement with a higher resolution has to be performed. The third part positions the sensors and starts the refined measurements. The fourth part finally integrates the measured local dataset into the overall data mesh. We have constructed a laboratory setup capable of measuring surfaces with extensions up to 1500×1000×500 mm3 (in x, y, and z directions, respectively). Using this measurement system we are able to separate the fault-indicating structures on the surface from the global shape, and to classify the detected structures according to their extensions and characteristic shapes simultaneously. The level of fault-detection probability is applicable by input parameter control.



Real-time surface inspection of extended surfaces is needed in a variety of industrial applications. For instance, in all branches of transportation industries (aircraft, ship, car, etc.), objects containing surfaces with extensions of several meters have to be handled and controlled. These objects, like the part of a vertical tail of an aircraft shown in Fig. 1(a), are either assembled from a number of subparts or, like the engine hood of a car [Fig. 1(b)], are produced as one part.

Fig. 1

Examples of extended technical surfaces. The part of a vertical tail of an aircraft shown in (a) with dimensions 3×2 m is assembled from a number of subparts. The engine hood of a car in (b) with dimensions 1.5×1.5 m is built from one part.


Though the objects in Fig. 1 are approximately from the same scale, they illustrate different stages of an industrial production process, which is shown in principle in Fig. 2. For both kinds of objects, quality control has to proof the overall shape to assure the intermediate and final assembly of the vehicle. That means that the objects have to fulfill certain tolerances concerning their geometric characteristics to guarantee a close and well adjusted connection of the parts. But the microstructure also has to be examined, to screen out defects that lower stability and reliability. Furthermore, both types of deviations can influence the optical appearance and, in this sense, reduce the quality of the product. If this kind of quality control is done in every stage of the production process (Fig. 2), and if the control measurements are integrated into the production line, the overall product quality will be increased.

Fig. 2

Principle of a controlled industrial production process. Every part is subjected to a quality control process, in which the ones that are out of the tolerances and cannot be modified sufficiently are sorted out and replaced by others. The tested parts are assembled to each other until they form a new part, which is subjected anew to an adapted quality control. This process is iterated to the point at which the final product is constructed.


The desired shape of an object is determined by the location of the characteristic surface structures to each other. If, for instance, the object under examination has been built up using a CAD model, the desired shape is given by this model. Deviations from the model can either influence the overall shape, the waviness (lower frequency contributions), the roughness (higher frequencies), or a mixture of them all.1 For the detection of deviations from the desired shape, a measurement system is needed that is capable of measuring the whole object surface. The resolution has to be high enough, but it is not the critical factor. In contrast, high-resolution measurements are needed to analyze the surface for microstructural defects that are strongly localized. And a combination of both should be available to detect extended defects of small amplitudes.

For the verification of the global shape, an area-related resolution (ratio between depth resolution and lateral extension, see Sec. 2.2) of 10−4 to 10−5 is sufficient in most cases. With respect to an efficient application to quality control, an inspection system that is able to analyze an object surface for shape deviations and microstructural defects at the same time has to perform a robust measurement with an area-related resolution of less than 10−6. Speaking of efficiency, the measurement system must be fast enough to be integrated into the production line without too much delay in the production cycle.

Most of the existing technologies used for the measurement of extended and complex surfaces divide the surface into local areas of smaller size, in which a number of local measurements is performed.2 3 4 5 6 7 The local areas must be chosen with a sufficient overlap. For the registration of the datasets, either markers on the object surface are used or the coordinates of the local sensors in the global system have to be known very well. There are at least two important problems to this approach.

  • • The markers have to be applied onto the surface of the object. In most cases this has to be done by hand and takes time. Furthermore, the application to the surface is an interaction with the latter and may not be nondestructive or wanted.

  • • There might be an ambiguity of the registration, especially if open surfaces are chosen and no absolute position determination of the markers is performed.

As an alternative approach, we propose the measurement of the global surface at once in a single measurement, followed by the detection of critical surface areas. For the complete characterization of those critical areas, local measurements with a higher absolute resolution are performed. The integration of all measurement data leads to a description of the surface with a locally increased resolution.


Measurement Concept for Extended Object Surfaces

To overcome the lack of a measurement technique for the investigation of extended object surfaces in an acceptable time with the necessary resolution, the concept of scaled topometry8 9 10 has been introduced. To get a fast and robust measurement system, we have to focus on optical sensors. Because of several reasons, fringe projection systems have been chosen out of a variety of different optical sensor systems (Sec. 2.1). Due to the area-related resolution that characterizes a fringe projection system and limits its absolute resolution on extended surfaces (Sec. 2.2), a scale-independent algorithm for detection and classification of surface structures based on specific input parameters11 12 has been developed (Sec. 2.4). The analysis leads to a multitude of data with different resolutions that have to be merged to deliver a complete surface description according to these input parameters (Sec. 2.5).

2.1 Fringe Projection Systems

For the realization of fast and robust measurements, not all optical sensor systems are suitable. To achieve a fast registration of measurement data, the only competitive systems are those that use a 2-D matrix [e.g., a charge-coupled device (CCD) array] to grab data from the surface in a single exposure. Moreover, it is difficult to use coherent interferometric systems, because they are very sensitive to environmental influences (like vibration, background light, etc.). Among the remaining systems using incoherent illumination, we focus on fringe projection systems, because they are not only robust but also adjustable to a wide range of measurement areas.

The principle of a fringe projection system is described in Fig. 3(a). It basically consists of a projection matrix (using a liquid crystal device (LCD) or digital mirror device (DMD) array) that is capable of projecting fringes onto an object surface and an associated camera [using a CCD or complimentary metal-oxide semiconductor (CMOS) array], which takes images of the surface with the projected fringes. Projector and camera direction include the triangulation angle θ. Typically sinusoidal fringe patterns with different wavelengths g are used for illumination to enable a precise phase measurement. Those patterns form a grid with a period g from which the height variation Δz in camera direction can be calculated [Fig. 3(b)]. From the setup geometry (projector and camera position relative to each other, the imaging optics of both, etc.) and the color information of the pixels in the corresponding images, 3-D surface information can be derived. This is done by the evaluation of the phase information, which allows the identification of the fringe order with the surface position.13 14 15

Fig. 3

Principle of a fringe projection system. (a) The system basically consists of a projection matrix and an associated camera separated by an angle θ. (b) Supposing a telecentric camera position and using this triangulation angle, the relations between the fringe period g and the height variation Δz in the camera direction can be calculated. (c) During the measurement, sinusoidal fringes are projected onto the surface. The digitized projection has to take care of the sampling theorem, i.e., the period of the sine wave has to contain a minimum number of pixels to be unambiguous.


From Fig. 3(c) it can be seen that the sinusoidal fringes must have a minimum wavelength g, because the sampling theorem has to be fulfilled. This means that the discretization of the sine wave has to be performed using as many pixels (three in case of a sine wave) for a period to make the representation unambiguous.

2.2 Area-Related Resolution

In optical systems, the resolution of the setup is defined by the capability to laterally distinguish two objects from each other. For those systems, the minimal distance δx between two such objects on the surface to be resolved is

Eq. (1)

δx0.61λn sin α,
where λ is the wavelength of the light and n⋅ sin(α) is the numerical aperture of the optical system with the refractive index n [Fig. 4(a)].

Fig. 4

Dependency of the lateral and height resolution on the measurement setup.


A change Δz in the height on the object surface as well can be described by a displacement Δx in the lateral direction that depends on the triangulation angle θ [Fig. 4(b)]. As a result, using Eq. (1), the height resolution δz can be described as

Eq. (2)

δz=δx tan θ0.61λn sin α tan θ.
Inserting typical values of a fringe projection system, we receive for the vertical resolution
δzFP0.61800nm1 sin 0.7 tan 0.5<1.5μm.
This is the lower bound of resolution that is theoretically achievable with a typical fringe projection setup. It could even be lowered if further influences concerning the height measurement were taken into account.

But in a fringe projection system, the limiting factor concerning the vertical and lateral resolution is mainly given by the phase error Δφ. Talking about the phase error Δφ, we have in mind the shift of the projected pixels on the surface. Such a shift of 2π would result in the original pattern shifted by a wavelength g [Fig. 3(c)]. Using the relation for the relative errors Δφ/2π=Δg/g and taking into account the geometric properties in Fig. 3, the height error δz is then given by

Eq. (3)

δz=δg sin θ=Δφg2π sin θ.
Another restriction is introduced by the sampling theorem. To display a sinusoidal fringe pattern nonambiguously, it is necessary to have a minimum number of pixels per period g. The resulting relation for the period is given by

Eq. (4)

gNxobjicam cos θ,
where N is the number of pixels per period, x obj is the width of the measurement area in the direction of the fringe modulation, and i cam is the number of camera pixels in the same direction. Inserting Eq. (4) into Eq. (3) leads to an estimation of the minimal height error

Eq. (5)

δzNΔ.φ2πicam tan θxobj.
We use the estimation for δz as a boundary condition for the achievable resolution. And as can be seen from Eq. (5), there is proportionality between δz and x obj . This relation is used to define the area-related resolution

Eq. (6)

δzA=δzANΔ.φ2πicam tan θ(Axobj),
with A as the measurement area. For typical values representing a fringe projection system, we get

Eq. (7)

δzA(FP)30.052π1024 tan 0.5>105.
This value is derived under the assumption of an ideal setup. A typical value, which is achievable in an optimized setup, is δz≈3⋅10−4.

2.3 Resolution Enhancement

From Eq. (6), it is obvious that for a calibrated system, the area-related resolution is almost constant. In Fig. 5 the consequence for two measurements on areas with different lateral extensions is demonstrated. It can be concluded that with an adequate focusing optic, a measurement system is capable of acquiring both datasets.

Fig. 5

Suppose that the global dataset (a) of an arbitrary surface has been measured by a calibrated system with an area related resolution δzA≈5⋅10−4. With an adequate focusing optic, this system is capable of acquiring the local dataset in (b), which is a part of the global surface marked by the white rectangle in (a). It is obvious that the absolute resolution Δz increases, while the area-related resolution δzA remains the same.


To improve the global resolution Δz obj for a given test object according to the input requirements (detailed information about the object shape, dimensions of smallest features to detect, etc.), successive measurements with decreasing sizes of surface areas are performed in areas of interest indicated by a detection algorithm (Sec. 2.4). The precondition for this approach is the assumption that the global measurement already contains information about the smallest features to detect. This defines minimal conditions for the relation between the resolution (lateral and vertical) of the global measurement and the feature size.

  • • The lateral extensions of the structures to be detected have to be at least in the range of the lateral resolution.

  • • The vertical extension must be bigger than the depth resolution Δz.

Features just fulfilling these minimal conditions are not fully described by the global measurement. The successive local measurements are used to locally increase the resolution and achieve feature characterization.

The aim of a quality control is to ensure tolerances of parts and products (Sec. 1). To ensure tolerances, the measurement accuracy (standard deviation) has to be about ten times better. This again implies a measurement resolution that is ten times better then the standard deviation. Taking into account this rule of thumb,10 i.e.,

Eq. (8)

tolerance  10. accuracy  100. resolution,
we can conclude that with this measurement technique we can track features with tolerances in the range of the measurement resolution. In this sense, the area-related resolution of the global system can be increased locally by 2 orders of magnitude. If we assume a high-resolution optical measurement system with an area-related resolution of 10−4 [see Eq. (7)] and use Eq. (8), we can achieve a depth resolution of less than 10 μm on an area of 1 m 2.

2.4 Scale-Independent Detection Algorithm

The areas, in which surface features on a certain scale may be located, are indicated by an analyzing algorithm. It distinguishes between different feature sizes down to lateral dimensions in the magnitude of the resolution. It basically consists of two transformations, a wavelet transformation (WT) and a fractal decomposition (FD), each followed by a further analytic investigation in the areas of interest. As a result of this analyzing algorithm, we get a list of surface areas on different scales, which probably contain faults or critical features (Fig. 6).

Fig. 6

Principle of the detection algorithm. From the input dataset, according to the size, orientation, and characteristics of the features of interest, wavelet transformations on the corresponding scales are performed. These transformations are tracked for the areas, which are pointed out by a fractal analysis comparing points to their surrounding. The resulting scale map indicates local areas with features of interest. If necessary, refined measurements in the indicated areas are performed.


The analysis result on the first scale is taken as the starting point for further local measurements with a higher absolute resolution.

In the first part of the investigation, the discrete surface data z(xi,yj) is subjected to a discrete wavelet transformation

Eq. (9)

WTz,ψ(ςi, ηj,s)=i,jz(xi,yj)ψ*(xiςis,yjηjs),
with a wavelet kernel ψ and a scaling factor s. For a single value of s, we receive a matrix i×j with the same size as the matrix of the input data. It contains the surface data z convolved with the wavelet kernel ψ, i.e., the localized surface frequencies in the scale of s. In the areas indicated by the fractal decomposition, which is explained later, every point (xi,yj) is analyzed as a continuous function over the scale s (Fig. 7).

Fig. 7

(a) Discrete surface data in a matrix i×j are subjected to (b) a discrete wavelet transformation WTz,ψ(xi,yj,s). If indicated by the fractal decomposition analysis, (c) the corresponding points (xi,yj) are further analyzed as a continuous function over the scale s.


The advantages of the wavelet transformation over a simple Fourier transformation concerning the separation of the surface data into the different scales consist not only in the localization of the frequency information, which is lost in the FT. In addition, the wavelet kernel ψ can be chosen out of a variety of different kernels, and in this sense can be adapted to the surface characteristics.

On the scales of interest, which are defined by the input parameters to the analysis process, a fractal decomposition, called the fractal pyramid, is performed on the wavelet-transformed data, respectively (Fig. 8). This process consists of the determination of the fractal dimension DF. In the field of image processing, where we deal with arrays of pixels and color value as height information, we can approximate the fractal dimension by creating parallelepipeds with volumes V=a⋅r3. The number N(r) of such volumes that we need on every scale r to completely cover the surface can be used to determine the fractal dimension by box counting:

Eq. (10)

DF(BC)=3 log [N(r)ar3] log (r)= log [aN(r)] log (r).
The fractal dimension of a surface area A is calculated by subdividing the area in square grids of sizes r2 and counting the number of parallelepipeds with height a⋅r to completely cover a surface [Fig. 8(b)]. This is done for a number of different surface area sizes, from the global surface down to local areas with only some pixels of size. The analysis of the “fractal dimension” of every pixel in the surrounding of surface areas of decreasing sizes leads to a classification of the pixel compared to its neighbors [Fig. 8(c)]. Critical surface areas are indicated by significant changes in the behavior throughout the different scales c of the fractal decomposition. This behavior is compared to the behavior of the pixel through the scales s of the wavelet transformation [Fig. 7(c)] and gives indications for successive local measurements with a higher resolution.

Fig. 8

On (a) the wavelet transformed data on the scale s, (b) a fractal decomposition is performed, taking into account local areas of size scale c±Δc. (c) For single points, the correspondence to their surroundings (spatially closest points) is analyzed.


2.5 Integration of Data with Different Resolutions

The datasets returned from the local measurements have to be integrated into the global dataset. For this purpose, we use a software that is based on an optimization of the standard iterative closest point (ICP) algorithm.16 17 This algorithm analyzes the closest point relations and minimizes the root mean square error by translating and rotating the point clouds. It returns translation and rotation matrices for the different point clouds, and is able to merge the point clouds into one single dataset.

The example in Fig. 9 shows the integration of two datasets. The algorithms work fast and robust under two constraints.

Fig. 9

(a) Measurement data of the back side of a former German coin. (b) Measurement data with a higher resolution of a local area of the same coin. (c) Integration of the two datasets with different resolutions. The position of the local dataset relative to the global one has to be known as good, as there is an overlap of the two letters “k.”

  • • The closest local minimum for the integration of the datasets is equivalent to the global one. This means that the two datasets should not be separated too far from each other to guarantee the merging into the right areas.

  • • Low and high resolutions do not differ much from each other, i.e., there are enough data points in the area of the low-resolution dataset into which the high-resolution local dataset should be integrated, to direct the merging to a high-precision registration.

The first condition is a matter of setup and configuration of the measurement systems (see Sec. 3.1). The second one has to be regarded using the results of the detection analysis. The measurement areas of the refined measurements should not be too small. In cases where areas of just a few pixels in the global measurement are detected, it could be useful to perform an intermediate measurement as a link between global and local datasets.


Experimental Results

3.1 Laboratory Setup

In our laboratory we have built up a test system that is based on an aluminum construction with dimensions of 2×1.5×2 m 3. The bottom consists of a solid aluminum plate on which the test objects can be placed. On the upper end of the frame, a xyz-positioning system is integrated that combines three axes mounted on top of each other. The axes are used for the movement of the local sensor systems, which can be placed at the end of the z axis (Fig. 10).

Fig. 10

Laboratory setup for the measurements using the concept of scaled topometry: the aluminum construction with dimensions 2×1.5×2 m 3 contains a ground plate, on which the objects to be measured can be placed. Different kinds of local sensor systems can be mounted onto the z axis of the xyz-positioning system. As global fringe projection system, (a) a color DMD-projector and (b) a corresponding color camera are fixed on top of the setup. (c) As a local adaptable fringe projection system, a compact system with integrated LCD projector and CCD camera is used.


The local sensor can be chosen according to application-specific requirements from a variety of different sensor types. It could be a confocal microscope,18 a laser triangulation sensor (point or line),19 a microscopic fringe projection sensor,20 21 22 or, as we use it, a further fringe projection sensor, which is adjustable for a range of measurement areas23 [Fig. 10(c)].

The global measurement system consists of a macroscopic fringe projection sensor with a DMD projector [Fig. 10(a)] and a CCD camera [Fig. 10(b)]. This sensor is placed on the very top of the frame at a fixed position. After calibration, it is capable of measuring the entire surface with a definite resolution. The measurement result, consisting of a 3-D dataset and a corresponding mask excluding invalid values, is transferred together with the information of interesting feature sizes and orientation to the analyzing algorithm. The resulting scale maps are used for the repositioning of the local sensor system that measures the relevant local datasets one after the other. Because the coordinates of the global system are well known and those of the local system can be obtained by the positioning system, the different datasets can be integrated by the standard matching software PolyWorks™.24

3.2 Measurement of a Car Engine Hood

The measurement technique is demonstrated using a cooperative test object from the automotive industry—a car engine hood partly deformed by a traffic accident [photo in Fig. 1(b) and data in Fig. 11(a)]. Due to the car crash, the test object contains a lot of surface defects on different scales. We demonstrate the analysis for the scaled topometry at an area with small gradients in the global shape, in which several defects of different kinds are located. In this area (256×256 mm 2), a local measurement is performed [Fig. 11(b)]. The dataset is transformed by the analysis process that has been described in Sec. 2. The results of this analysis process and the consecutive refined measurements are shown in Fig. 12.

Fig. 11

Measurement data of the test object. (a) The engine hood has been measured by the global fringe projection system. (b) The area framed by the square (256×256 mm 2) in (a) has been measured with a higher resolution by a local sensor system.


Fig. 12

Two branches of the data analysis from the local area measurement (256×256 mm 2) of the selected part of the engine hood (left and right side). Two scales of the wavelet transformation (a) and (d), the result of the corresponding fractal decompositions (b) and (e), and the succeeding further refined measurement (c) and (f) are shown, respectively. Only a few features on a coarse scale (left) and a lot of surface features on a finer scale (right) are detected by the fractal pyramid.


Figure 12 shows two different scales of the surface topometry and the fault indicating structures inside them: one on the left and the other on the right side. One can see two scales of the wavelet transformation [Figs. 12(a) and 12(d)] and the result of the corresponding fractal decompositions [Figs. 12(b) and 12(e)]. On the bottom, the corresponding succeeding measurements [Figs. 12(c) and 12(f)] are depicted, respectively.

On the left side, the analysis for the scale of 1/8 of the extension of the dataset is shown. Only three structures are contained in the upper part of the dataset, out of which the black area in the fractal decomposition [Fig. 12(b)] is subjected to a refined measurement. On the right side, the results of the scale 1/64 are depicted. For further analysis, a feature was chosen out of a large amount of others [black and gray squares in the fractal decomposition of Fig. 12(e)].

The high-resolution measurements show that the groove [Fig. 12(c)] is a cut that originates from a beat or a bump with a massive object. The cut is surrounded by an area where the color is partly removed because of corrosion. The structure analyzed in the right part [Fig. 12(f)] is a convex feature, in contrast to the other one, and has been characterized as a deposition consisting of rubber.

The input parameters to the detection algorithm (detailed information about the object shape, dimensions of smallest features to detect, etc.) determine the scales on which surface feature information has to be analyzed. According to these parameters, the local measurements are performed, and the resulting datasets are integrated into a global one. As an example, we show the integration of one of the local measurements (Fig. 13).

Fig. 13

Registration of the local measurement data [Fig. 12(c)] into the refined measurement [Fig. 11(b)] and the global dataset [Fig. 11(a)]. Especially in the magnification in (b), the difference in the resolutions of the datasets can be seen.


The final dataset is a 3-D patchwork of surface information with different resolutions. Higher resolution local measurements are only performed in areas where the analysis process indicates deviations of the global shape that are out of the tolerances allowed for the tested product.



We emphasize the importance of quality control for industrial applications. Integration into the production line in every stage of the production process is desirable with regard to an efficient product output. To achieve these objectives for extended object surfaces, we demonstrate a multisensor inspection system using the concept of scaled topometry. The need for a new measurement approach originates from the fact that a single measurement is not capable of measuring the whole extended surface with a sufficient resolution in an acceptable time.

This inspection system performs a global measurement with the highest possible resolution. The data analysis returns fault-indicating structures on different scales of resolutions down to the required minimum feature extensions. If these smallest features can not be resolved by the global measurement, and if there are at least indications about them in the global dataset, consecutive local area measurements with higher absolute resolutions are performed and analyzed in the same way. The resulting datasets are merged by matching software and integrated into one overall dataset that describes the object surface completely according to the input requirements.

We illustrate the mode of operation using a car engine hood. This cooperative object is a typical surface that needs a better testing system for industrial quality control. To demonstrate the multiscale ability of our approach, we chose an object with many faults (car crash). The results are not representative concerning an industrial application, but they are convincing, especially in the background stage of this early laboratory setup.


This work has been supported by the State of Bremen as project AMST-P7 within the AMST initiative. We would like to thank Mr. R. Windecker from the Institut fu¨r Technische Optik, Universita¨t Stuttgart for the measurements with the microscopic fringe projection sensor [Figs. 12(c) and 12(f)].



D. J. Whitehouse, Handbook of Surface Metrology, Institute of Physics Publishing Ltd., Philadelphia, PA (1994).


Automatic body inspection system (ABIS), Steinbichler Optotechnik GmbH, Germany. For more details, see


H. Steinbichler, E. Heinrich No¨sekabel, and R. Ro¨sch, “Optical inspection in the production line,” in Fringe 2001—4th Intl. Workshop Automatic Process. Fringe Patt. W. Osten and W. Ju¨ptner, Eds., pp. 587–592, The Data Science Library, Elsevier, Paris (2001).


Self-calibrating 3-D measuring system, IOF Jena, Germany. For more details, see


G. Notni, “360-deg shape measurement with fringe projection—calibration and application,” in Fringe 2001-4th Intl. Workshop Automatic Process. Fringe Patt. W. Osten and W. Ju¨ptner, Eds., pp. 311–323, The Data Science Library, Elsevier, Paris (2001).


G. Wiora , “High resolution measurement of phase-shift amplitude and numeric object phase calculation,” Proc. SPIE , 4117 289 –299 (2000). Google Scholar


G. Wiora, “Optical 3D-metrology: Precise shape measurement with an extended fringe projection method,” PhD Thesis, Interdisziplina¨res Zentrum fu¨r Wissenschaftliches Rechnen (IWR), Fakulta¨t fu¨r Pysik und Astronomie, Universita¨t Heidelberg (2001).


P. Andra¨, E. Ivanov, and W. Osten, “Scaled topometry—an active measurement approach for wide scale 3D surface inspection,” in Fringe 1997 3rd Intl. Workshop Automatic Process. Fringe Patt., W. Osten and W. Ju¨ptner, Eds., pp. 179–189, Akademie-Verlag, Berlin (1997).


D. Kayser , T. Bothe , and W. Osten , “Fault detection in gray-value images of surfaces on different scales,” Proc. SPIE , 2744 110 –117 (1999). Google Scholar


W. Osten , P. Andra¨ , and D. Kayser , “Highly-resolved measurement of extended technical surfaces with scalable topometry,” Tech. Mess. , 66 (11), 413 –428 (1999). Google Scholar


D. Kayser, W. Osten, and W. P. O. Ju¨ptner, “A scale independent algorithm for the detection of fault indicating structures in range images,” in Interferometry in Speckle Light: Theory and Applications, Proc. Intl. Conf. Lausanne, P. Jacquot and J. M. Fournier, Eds., pp. 389–396, Springer Verlag, Berlin (2000).


D. Kayser, T. Bothe, and W. Osten, “An integrated measurement system for the inspection of extended surfaces in industrial quality control,” in Proc. Intl. Symp. Photon. Meas. Aachen, pp. 339–344, VDI-Verlag, Du¨sseldorf (2002).


V. Srinivasan , H. C. Liu , and M. Halioua , “Automated phase-measuring profilometry of 3-d diffuse objects,” Appl. Opt. , 23 3105 –3108 (1984). Google Scholar


W. Osten and W. Ju¨ptner , “Measurement of displacement vector fields of extended objects,” Opt. Lasers Eng. , 24 261 –285 (1996). Google Scholar


F. Chen , G. M. Brown , and M. Song , “Overview of three-dimensional shape measurement using optical methods,” Opt. Eng. , 39 10 –22 (2000). Google Scholar


P. J. Besl and N. D. McKay , “A method for registration of 3-D shapes,” IEEE Trans. Pattern Anal. Mach. Intell. , 14 (2), 239 –256 (1992). Google Scholar


D. W. Eggert , A. W. Fitzgibbon , and R. B. Fisher , “Simultaneous registration of multiple range views for use in reverse engineering of CAD models,” Comput. Vis. Image Underst. , 69 (3), 253 –272 (1998). Google Scholar


H. J. Tiziani , R. Windecker , M. Wegner , K. Leonhardt , D. Steudle , and M. Fleischer , “Measurement and description of microstructures with consideration of material-specific characteristics,” Tech. Mess. , 66 (11), 429 –436 (1999). Google Scholar


Keyence Corporation, Osaka, Japan. For more details see


R. Windecker , M. Fleischer , and H. J. Tiziani , “Three-dimensional topometry with stereo microscopes,” Opt. Eng. , 36 (12), 3372 –3377 (1997). Google Scholar


K. Ko¨rner , R. Windecker , M. Fleischer , and H. J. Tiziani , “One-grating projection for absolute three-dimensional profiling,” Opt. Eng. , 40 (8), 1653 –1660 (2001). Google Scholar


R. Windecker , M. Fleischer , K. Ko¨rner , and H. J. Tiziani , “Testing micro devices with fringe projection and white-light interferometry,” Opt. Lasers Eng. , 36 (2), 141 –154 (2001). Google Scholar


K. Ko¨rner , U. Droste , R. Windecker , M. Fleischer , H. J. Tiziani , T. Bothe , and W. Osten , “Projection of structured light in object planes of varying depths for absolute 3-D profiling in a triangulation setup,” Proc. SPIE , 4398 23 –34 (2001). Google Scholar


InnovMetric Inc., Canada. For more information see

008410j.m1.jpg Daniel Kayser studied physics at the University of Bremen. As a member of the Institute of Solid State Physics, he received his MS degree in 1998 working on the growth of ZnSe/ZnCdSe epitaxial layers. Since 1998 he has been with the Bremer Institut fu¨r angewandte Strahltechnik as a member of the optical metrology group in the field of optical 3-D measurement and image processing.

008410j.m2.jpg Thorsten Bothe studied physics at Carl von Ossietzky University, Oldenburg. As a member of the applied optics group he received his MS degree in 1995 working on speckle pattern interferometry. Afterward, he was project manager for deformation monitoring in historical monuments by 3-D ESPI. Since 1998 he has been with Bremer Institut fu¨r angewandte Strahltechnik as member of the optical metrology group in the area of shape and deformation measurement, and has been working on his thesis.

Wolfgang Osten received the BSc from the University of Jena in 1979. From 1979 to 1984 he was a coworker of the Institute of Mechanics in Berlin. In 1983 he received the PhD degree from the University of Halle in the field of holographic interferometry. From 1984 to 1991 he was employed at the Central Institute of Cybernetics and Information Processes, Berlin, making investigations in digital image processing and computer vision. Between 1991 and 2002 he was employed at the Bremen Institute of Applied Beam Technology as the head of the Department of Optical Metrology. Since September 2002 he has been a professor at the University of Stuttgart and director of the Institut fu¨r Technische Optik. He is concerned with new concepts in industrial inspection by combining modern principles of optical metrology and image processing.

©(2004) Society of Photo-Optical Instrumentation Engineers (SPIE)
Daniel Kayser, Thorsten Bothe, and Wolfgang Osten "Scaled topometry in a multisensor approach," Optical Engineering 43(10), (1 October 2004).
Published: 1 October 2004

Back to Top