rqtreemesh 1.0.0

An implementation of restricted quadtree triangulation (RQT), for converting raster heightmaps to 3D meshes.

Triangulated terrain near Horseshoe Bend, Arizona.

rqtreemesh

rqtreemesh is a basic implementation of restricted quadtree triangulation (RQT), for converting raster heightmaps to 3D meshes. For an overview of RQT-based mesh generation, see

Pajarola, Renato. "Overview of Quadtree-based Terrain Triangulation and Visualization." (2002).

Implemented in C++ with a pybind11 wrapper for Python.

Quick Start

git clone --recursive this repo and run pip install . on it.

This package requires Numpy and Pillow.

Usage

Functionality is provided by the Heightmap class.

from rqtreemesh import Heightmap

Two constructors are available:

heightmap = Heightmap(array: np.ndarray, pixel_dim: float, top_left_x: float, top_left_y: float)

Parameters

array: np.ndarray

Input raster heightmap as a 2-dimensional numpy array where each dimension is a power of two. E.g. 1024 x 2048. Maximum size is 32768 x 32768.

 pixel_dim: float

The spatial dimension of one pixel in the input heightmap.

top_left_x:  float, top_left_y:  float

The x and y locations of the top left corner of the input heightmap. Note that this location is not the center of the top left pixel, but rather its top left corner.

---

Alternatively:

heightmap = Heightmap.from_geotiff(path: str)

Parameters

path: str

Path to input raster heightmap in geoTIFF format with dimesnsions a power of two. E.g. 1024 x 2048. Maximum size is 32768 x 32768. The geolocation and spatial resolution of the heightmap will be obtained from the appropriate tags in the geoTIFF.

---

To generate a mesh, use the generate_mesh method.

heightmap = Heightmap...
(vertices, triangles) = heightmap.generate_mesh(max_error: float, show_progress: bool = False) -> Tuple[np.ndarray, np.ndarray]

Parameters

max_error: float

The maximum absolute vertical distance between the input heightmap and the resulting mesh.

show_progress: bool = False

Show a progress bar while generating the mesh. Default value: False.

Returns

(vertices, triangles): Tuple[np.ndarray, np.ndarray]

A tuple containing the output mesh. vertices is an array where each row contains the x, y and z coordinates of a vertex. triangles is a connectivity matrix, i.e. an array where each row contains a triangle as three (zero-indexed) indices into the first element.

---

To provide easy compatibility with most 3D software, a simple exporter to the Wavefront Object (.obj) format is also provided:

from rqtreemesh import write_obj
write_obj(path: str, verts: np.ndarray, triangles: np.ndarray, flip_zy: bool = False)

Parameters

path: str

The output file path.

verts: np.ndarray, triangles: np.ndarray

A coordinate array and connectivity matrix, such as the output of generate_mesh.

flip_zy: bool = True

Interchange the z and y axes of the exported mesh. Default value: False.

---

Triangulation of a 2 meter per pixel resolution heightmap of East Pawnee Butte, Colorado. A maximum error of 0 meters results in a mesh consisting of 32768 triangles (left). A maximum error of 0.5 meters produces 8156 triangles (middle), and a maximum error of 1 meter produces 4074 triangles (right).