The Beauvais Cathedral Project ∗
Total Page:16
File Type:pdf, Size:1020Kb
The Beauvais Cathedral Project ∗ Peter K. Allen, Alejandro Troccoli, Benjamin Smith, Ioannis Stamos†and Stephen Murray‡ Department of Computer Science Department of Art History and Archaeology‡ Columbia University, New York, NY 10027 Department of Computer Science, Hunter College, CUNY† Abstract Saint-Pierre in Beauvais France [9], a unique and still used church which is a prime example ofhigh gothic Preserving cultural heritage and historic sites is an architecture (fig 1). A team ofarchitectural histo- important problem. These sites are subject to ero- rians, computer scientists, and engineers has begun sion, vandalism, and as long-lived artifacts, they have to study the fragile structure of the tallest medieval gone through many phases ofconstruction, damage cathedral in France . The thirteenth-century Gothic and repair. It is important to keep an accurate cathedral at Beauvais collapsed three times in the record ofthese sites using 3-D model building tech- Middle Ages and continues to pose problems for its nology as they currently are, so preservationists can long-term survival. Our group is creating a highly ac- track changes, foresee structural problems, and allow curate three-dimensional model based on laser scans a wider audience to ”virtually” see and tour these ofthe interior and exterior ofthe building. This in- sites. Due to the complexity ofthese sites, building formation will be used to examine the weaknesses in 3-D models is time consuming and difficult, usually the building and propose remedies; visualize how the involving much manual effort. This paper discusses building looked when it was first built; and to serve new methods that can reduce the time to build a as the basis for a new collaborative way of teach- model using automatic methods. Examples ofthese ing about historic sites both in the classroom and on methods are shown in reconstructing a model ofthe the Internet. The building is included on the World Cathedral ofSaint-Pierre in Beauvais, France. Monuments Fund’s Most Endangered List. Although the cathedral survived the heavy incendi- 1 Introduction ary bombing that destroyed much ofBeauvais dur- Preserving cultural heritage and historic sites is an ing World War II, the structure is as dangerous as important problem. The advent ofnew digital 3D it is glorious, being at risk from flaws in its original scanning devices has provided new means to pre- design, compounded by differential settlement and serve these sites digitally, and to preserve the historic with stresses placed on its flying buttresses from gale record by building geometric and photorealistic 3D force winds off the English Channel. The winds cause models. A recent international conference (Virtual the buttresses to oscillate and already weakened roof and Augmented Architecture ‘01 [6]) has highlighted timbers to shift. Between the 1950s and 1980s nu- the need for simple and automatic methods to create merous critical iron ties were removed from the choir rich models ofhistoric environments. buttresses in a damaging experiment. A temporary tie-and-brace system installed in the 1990s may have A number ofother projects have addressed this and made the cathedral too rigid, increasing rather than similar problems including [17, 8, 5, 4, 3, 11]. This decreasing stresses upon it. Although the cathedral paper discusses new methods we have developed to has been intensively studied, there continues to be a recover complete geometric and photometric mod- lack ofconsensus on how to conserve the essential vi- els oflarge sites and to automate this process. In sual and structural integrity ofthis Gothic wonder. particular, we discuss new methods for data abstrac- With its five-aisled choir intersected by a towered tion and compression through segmentation, 3D to transept and its great height (keystone 152.5 feet 3D registration (both coarse and fine), and 2D to 3D above the pavement), Beauvais Cathedral, commis- texture mapping ofthe models with imagery. sioned in 1225 by Bishop Milon de Nanteuil, provides The testbed for our methods is the Cathedral of an extreme expression ofthe Gothic enterprise. ∗This work was supported in part by NSF grant IIS-01- 21239 and the Samuel Kress Foundation. point cloud is aligned relative to the other clouds at the position and orientation corresponding to the po- sition occupied by the physical surface it represents on the actual structure Although the point clouds may be registered manually, it is very time consum- ing and error-prone. Manually visualizing millions ofsmall points and matching them is quite impre- cise and difficult as the number ofscans increases. When possible, it is a common practice to use spe- cially designed targets/fiducials to help during the registration phase. In our case, however, it was al- most impossible to place targets higher than 2.5 me- ters above the ground, requiring us to develop an automatic registration method. Our registration method is a three step process. The first step is an automatic pairwise registration be- tween two overlapping scans. The pairwise registra- tion matches 3-D line segments extracted from over- lapping range scans to compute the correct transfor- mation. The second step is a global registration step that tries to align all the scans together using over- lapping pairs [15]. The third step is a multi-image Figure 1: Cathedral of Saint-Pierre simultaneous ICP algorithm [2] that does the final fine registration ofthe entire data set. 3.1 Segmenting the Scans 2DataAcquisition Our previously developed range segmentation algo- rithm [13, 14] automatically extracts planar regions The surveying process involved capturing range and and linear features at the areas of intersection of intensity data. On site work started in June 2001 neighboring planar structures. Thus, a 3–D range using a Cyrax 2400 scanner. The whole interior and scan is converted to a set ofbounded planes and a 1/3 ofthe exterior were surveyed, but due to techni- set offinite lines. The extracted 3–D intersection cal problems, the scanning process was aborted and lines are very accurate because their orientation and resumed in June 2002, this time with a new Cyrax length depends on all the points ofboth planes that 2500 scanner. Over 200 range images were acquired, intersect. 120 interior and 100 exterior scans, most ofthem sampled at 1000 by 1000 points. Intensity images The pairwise registration algorithm efficiently com- were captured by a 5 megapixel digital camera, which putes the best rigid transformation (R, T ) between a was freely placed. pair ofoverlapping scans S1 and S2. This transfor- mation has an associated grade g(R, T )thatequals The scanner uses a time-of-flight laser technology to the total number of line matches after the transfor- measure the distance to points on an object. Data mation is applied. Note that the grade is zero ifthere from the scanner comprises point clouds, each point is no overlap between the scans. comprising four coordinates, (x, y, z) and a value rep- resenting the amplitude ofthe laser light reflected In a typical scanning session tens or hundreds of back to the scanner. The amplitude is dependent on range scans need to be registered. Calculating all the reflectance ofthe material and the surfaceorien- possible pairwise registrations is impractical because tation. it leads to a combinatorial explosion. In our system, the user is providing a list ofoverlapping pairs of 1 3 Registration of Range Scans scans . All pairwise transformations are computed. Then, one ofthe scans is chosen to be the anchor scan In order to acquire data describing an entire struc- Sa. Finally, all other scans S are registered with re- ture, such as the Beauvais Cathedral, multiple scans spect to the anchor Sa. In the final step, we have must be taken from different locations with differ- the ability to reject paths ofpairwise transformation ent scanner orientations. To reconstruct the origi- that contain registrations oflower confidence. nal structure, the different scans must be registered 1Note that the list does not have to be a complete list of together correctly. When correctly registered, each all possible overlaps. In more detail, the rigid transformations (Ri,Ti)and are searched. Then the RSV distance (to the model their associated grades g(Ri,Ti) are computed be- point) for each of them is evaluated to get the closest tween each pair ofoverlapping scans. In this man- point. ner a weighted undirected graph is constructed. The Once correspondences are made, least-squares is typ- nodes ofthe graph are the individual scans, and the ically used to find the correct transformation matri- edges are the transformations between scans. Finally ces for the data points. However, in the presence the grades g(Ri,Ti) are the weights associated with ofoutliers, least-squares can be unstable. Accord- each edge. More than one path ofpairwise transfor- ingly, an M-estimator is used to weight the data mations can exist between a scan S and the anchor points. The procedure is iterative, using a conjugate- Sa. Our system uses a Dijkstra-type algorithm in or- gradient search to find the minimum. der to compute the most robust transformation path Figure 3b shows the reduction oferror vs. iteration from S to Sa.Ifp1 and p2 are two different paths for a known test data set. The data set of Beauvais from S to Sa,thenp1 is more robust than p2,if Cathedral contains over 100 scans, and it requires the cheapest edge on p1 has a larger weight than the significant computational resources and time to reg- cheapest edge of p2. This is the case because the cheapest edge on the path corresponds to the pair- ister these scans with full resolution; therefore, those wise transformation of lowest confidence (the smaller scans are sub-sampled down to 1/25 oftheir original the weight the smaller the overlap between scans).