## 2 Materials and Methods
### 3D Modelling
#### 2.4.1 3D Model of Houses
The BC is combined with the information on the division of the layer into zones to capture the importance of the royal streets in the city (Carrera, Herradores, San Agustin, del Agua) as seen in Figure 10.  
Figure 10: Importance of each street in the city using BC and population strata zones. The most important streets are in red tones, followed by the streets in cyan tones, and concluding with the streets in purple tones. Figure 9: Probabilistic graphical model for house generation. The application of the modified SS algorithm is shown in Figure 11. To define the initial split for a block B, the user defines a perpendicular distance \(d_{\textit{offset}}\) from block contour C(B) to an inwards offset contour C\({}^{\prime}\)(B). In our case, \(d_{\textit{offset}}\) is defined by the interior rectangular side of the plot. The contour C\({}^{\prime}\)(B) is calculated by the application of the straight skeleton to C(B) [43]. The arcs of the skeleton application specify the division of the region between C(B) and C\({}^{\prime}\)(B) into a set of strip polygons that are referred to _a-strips_. These polygons are of three types: empty space, walls, or building land as seen in Figure 11a. Strips from empty spaces and walls are fixed throughout the algorithm. The _a-strips_ from building land are first joined based on collinearity. Then, triangular regions are transferred from adjacent strips with lower BC to those with higher BC as seen in Figure 11b. These corrected strips are referred to as _\(\beta\)-strips_. Finally, the _\(\beta\)-strips_ are subdivided using the length of plot fronts as seen in Figure 11c, and the modified SS algorithm ends. The detailed pseudocode of the algorithm can be found in Appendix A. To perform the division of the interior of the blocks, the RS method in [41] has been used, which performs an adaptive spatial partitioning of the block boundary using oriented bounding boxes (OBBs). In our case, the division of the blocks reaches the walls of the block or the edge of the houses, as shown in Figure 11. This block division corresponds to the internal walls of the blocks in [PERSON]'s map. Small constructions in the map are placed on the rectangles beside the internal walls that define the block's subdivisions, as shown in Figure 11d. Upon completion of block subdivision, a collection of plots is obtained, each defined by a specific spatial location and a basic house description denoted by its _fcode_ and _ftype_. To produce a complete geometric model of each house, additional geometric features are incorporated, such as its shape evolution (e.g., elongation or with protrusions in various configurations) as described in [40], as well as whether or not an attached construction is present (a common occurrence in [PERSON]'s map). The model also includes details about textures used for facades, wooden elements, or roofs. By incorporating these features, we produce an implicit 3D model of each house, as outlined in Table 4. We intend to preserve both the geometric and textural features independently of any software tools employed to obtain explicit 3D reconstructions. Figure 11: **(a)** Block initial strip subdivision (grey is reserved for walls, white is reserved for empty space, and orange is reserved for building land); **(b)** subdivision after BC; **(c)** parcel subdivision; and **(d)** interior block subdivision, with small parcels for interior houses.