Surface with multiple textures

This example demonstrates one possible method for displaying a 3D surface with multiple textures.

Thanks to Emmanuel Reynaud and Luis Gutierrez for providing the gorgeous coral model for this demo. You can find the data on FigShare: https://zenodo.org/records/13380203

More information on the methods used to generate this model can be found in L. Gutierrez-Heredia, C. Keogh, E. G. Reynaud, Assessing the Capabilities of Additive Manufacturing Technologies for Coral Studies, Education, and Monitoring. Front. Mar. Sci. 5 (2018), doi:10.3389/fmars.2018.00278.

A bit about 3D models

A standard way to define a 3D model (mesh, or Surface in napari) is by listing vertices (3D point coordinates) and faces (triplets of vertex indices - each face is a triangle in 3D space). Meshes are often stored in “Wavefront” (.obj) files, which may have companion material (.mtl) files that describe some shading properties (base color, shinyness, etc.) for different parts of the model.

In some cases, the color of a vertex is given by a single point value that is then colormapped on the fly (vertex_values). In other cases, each vertex or face may be assigned a specific color (vertex_colors). These methods are demonstrated in Surface with texture and vertex_colors.

In the case of “photorealistic” models, the color of each vertex is instead determined by mapping a vertex to a point in an image called a texture using 2D texture coordinates in the range [0, 1]. The color of each individual pixel is smoothly interpolated (sampled) on the fly from the texture (the GPU makes this interpolation very fast).

Napari does not (yet) support models with multiple textures or materials. If the textures don’t overlap, you can display them on separate meshes as shown in this demo. If the textures do overlap, you may instead be able to combine the textures as images. This relies on textures having the same texture coordinates, and may require resizing the textures to match each other.

Tags: visualization-nD

import os

import matplotlib.pyplot as plt
import pooch
from vispy.io import imread, read_mesh

import napari

Download the model

download = pooch.DOIDownloader(progressbar=True)
doi = '10.5281/zenodo.13380203'
tmp_dir = pooch.os_cache('napari-surface-texture-example')
os.makedirs(tmp_dir, exist_ok=True)
data_files = {
    'mesh': 'PocilloporaDamicornisSkin.obj',
    # "materials": "PocilloporaVerrugosaSkinCrop.mtl",  # not yet supported
    'Texture_0': 'PocilloporaDamicornisSkin_Texture_0.jpg',
    'GeneratedMat2': 'PocilloporaDamicornisSkin_GeneratedMat2.png',
}
print(f'downloading data into {tmp_dir}')
for file_name in data_files.values():
    if not (tmp_dir / file_name).exists():
        print(f'downloading {file_name}')
        download(
            f'doi:{doi}/{file_name}',
            output_file=tmp_dir / file_name,
            pooch=None,
        )
    else:
        print(f'using cached {tmp_dir / file_name}')

downloading data into /home/runner/.cache/napari-surface-texture-example
downloading PocilloporaDamicornisSkin.obj

  0%|                                              | 0.00/93.8M [00:00<?, ?B/s]
  0%|                                     | 4.10k/93.8M [00:00<43:05, 36.3kB/s]
  0%|                                     | 15.4k/93.8M [00:00<19:56, 78.4kB/s]
  0%|                                      | 48.1k/93.8M [00:00<10:32, 148kB/s]
  0%|                                       | 116k/93.8M [00:00<05:39, 276kB/s]
  0%|                                       | 254k/93.8M [00:00<03:05, 503kB/s]
  1%|▏                                      | 518k/93.8M [00:00<01:42, 909kB/s]
  1%|▍                                    | 1.07M/93.8M [00:01<00:52, 1.76MB/s]
  2%|▊                                    | 2.17M/93.8M [00:01<00:27, 3.39MB/s]
  5%|█▋                                   | 4.36M/93.8M [00:01<00:13, 6.60MB/s]
  9%|███▏                                 | 8.14M/93.8M [00:01<00:07, 11.8MB/s]
 12%|████▍                                | 11.2M/93.8M [00:01<00:05, 14.0MB/s]
 15%|█████▌                               | 14.1M/93.8M [00:01<00:04, 17.0MB/s]
 17%|██████▎                              | 15.9M/93.8M [00:01<00:05, 15.3MB/s]
 20%|███████▍                             | 18.8M/93.8M [00:02<00:04, 16.2MB/s]
 24%|████████▋                            | 22.0M/93.8M [00:02<00:04, 17.5MB/s]
 27%|█████████▊                           | 24.9M/93.8M [00:02<00:03, 19.9MB/s]
 29%|██████████▋                          | 27.0M/93.8M [00:02<00:03, 17.9MB/s]
 31%|███████████▋                         | 29.5M/93.8M [00:02<00:03, 17.6MB/s]
 34%|████████████▍                        | 31.6M/93.8M [00:02<00:03, 18.3MB/s]
 36%|█████████████▎                       | 33.7M/93.8M [00:02<00:03, 16.8MB/s]
 39%|██████████████▍                      | 36.7M/93.8M [00:03<00:03, 17.3MB/s]
 43%|███████████████▊                     | 40.2M/93.8M [00:03<00:02, 18.5MB/s]
 46%|█████████████████                    | 43.2M/93.8M [00:03<00:02, 18.4MB/s]
 49%|██████████████████▎                  | 46.3M/93.8M [00:03<00:02, 18.6MB/s]
 52%|███████████████████▍                 | 49.2M/93.8M [00:03<00:02, 18.4MB/s]
 56%|████████████████████▌                | 52.0M/93.8M [00:03<00:02, 18.2MB/s]
 59%|█████████████████████▊               | 55.3M/93.8M [00:04<00:02, 18.6MB/s]
 62%|██████████████████████▉              | 58.0M/93.8M [00:04<00:01, 20.5MB/s]
 64%|███████████████████████▋             | 60.2M/93.8M [00:04<00:01, 18.4MB/s]
 67%|████████████████████████▋            | 62.5M/93.8M [00:04<00:01, 17.5MB/s]
 69%|█████████████████████████▌           | 64.8M/93.8M [00:04<00:01, 18.6MB/s]
 72%|██████████████████████████▌          | 67.4M/93.8M [00:04<00:01, 20.5MB/s]
 74%|███████████████████████████▍         | 69.5M/93.8M [00:04<00:01, 18.3MB/s]
 77%|████████████████████████████▍        | 72.2M/93.8M [00:04<00:01, 18.0MB/s]
 80%|█████████████████████████████▊       | 75.5M/93.8M [00:05<00:00, 18.7MB/s]
 84%|███████████████████████████████▏     | 78.9M/93.8M [00:05<00:00, 22.4MB/s]
 87%|████████████████████████████████     | 81.3M/93.8M [00:05<00:00, 20.4MB/s]
 89%|████████████████████████████████▉    | 83.6M/93.8M [00:05<00:00, 20.8MB/s]
 92%|█████████████████████████████████▉   | 86.1M/93.8M [00:05<00:00, 20.0MB/s]
 95%|███████████████████████████████████  | 88.7M/93.8M [00:05<00:00, 21.5MB/s]
 97%|███████████████████████████████████▉ | 91.0M/93.8M [00:05<00:00, 19.4MB/s]
100%|████████████████████████████████████▉| 93.7M/93.8M [00:06<00:00, 18.7MB/s]
  0%|                                              | 0.00/93.8M [00:00<?, ?B/s]
100%|██████████████████████████████████████| 93.8M/93.8M [00:00<00:00, 540GB/s]
downloading PocilloporaDamicornisSkin_Texture_0.jpg

  0%|                                              | 0.00/17.3M [00:00<?, ?B/s]
  0%|                                     | 11.3k/17.3M [00:00<03:56, 73.0kB/s]
  0%|                                      | 39.9k/17.3M [00:00<02:07, 135kB/s]
  1%|▏                                     | 95.2k/17.3M [00:00<01:14, 230kB/s]
  1%|▍                                      | 205k/17.3M [00:00<00:42, 406kB/s]
  2%|▉                                      | 428k/17.3M [00:00<00:22, 756kB/s]
  5%|█▉                                    | 886k/17.3M [00:00<00:11, 1.46MB/s]
 10%|███▊                                 | 1.80M/17.3M [00:01<00:05, 2.82MB/s]
 21%|███████▋                             | 3.62M/17.3M [00:01<00:02, 5.48MB/s]
 38%|██████████████▏                      | 6.63M/17.3M [00:01<00:01, 9.55MB/s]
 57%|████████████████████▉                | 9.78M/17.3M [00:01<00:00, 12.6MB/s]
 74%|███████████████████████████▎         | 12.8M/17.3M [00:01<00:00, 14.4MB/s]
 92%|██████████████████████████████████   | 15.9M/17.3M [00:01<00:00, 15.8MB/s]
  0%|                                              | 0.00/17.3M [00:00<?, ?B/s]
100%|█████████████████████████████████████| 17.3M/17.3M [00:00<00:00, 98.8GB/s]
downloading PocilloporaDamicornisSkin_GeneratedMat2.png

  0%|                                               | 0.00/120k [00:00<?, ?B/s]
  6%|██▎                                   | 7.17k/120k [00:00<00:02, 49.2kB/s]
 32%|████████████▋                          | 38.9k/120k [00:00<00:00, 139kB/s]
 79%|██████████████████████████████▉        | 95.2k/120k [00:00<00:00, 236kB/s]
  0%|                                               | 0.00/120k [00:00<?, ?B/s]
100%|████████████████████████████████████████| 120k/120k [00:00<00:00, 542MB/s]

Load the model

Next, read the model data from the .obj file. Currently napari/vispy do not support reading material properties (.mtl files) nor separate texture and vertex indices (i.e. repeated vertices). Normal vectors read from the file are also ignored and re-calculated from the faces.

vertices, faces, _normals, texcoords = read_mesh(tmp_dir / data_files['mesh'])

Notice reading .OBJ: material properties are ignored.

Load the textures

This model comes with two textures: Texture_0 is generated from photogrammetry of the actual object, and GeneratedMat2 is a generated material to fill in parts of the model lacking photographic texture.

photo_texture = imread(tmp_dir / data_files['Texture_0'])
generated_texture = imread(tmp_dir / data_files['GeneratedMat2'])

This is what the texture images look like in 2D:

fig, axs = plt.subplots(1, 2)
axs[0].set_title(f'Texture_0 {photo_texture.shape}')
axs[0].imshow(photo_texture)
axs[0].set_xticks((0, photo_texture.shape[1]), labels=(0.0, 1.0))
axs[0].set_yticks((0, photo_texture.shape[0]), labels=(0.0, 1.0))
axs[1].set_title(f'GeneratedMat2 {generated_texture.shape}')
axs[1].imshow(generated_texture)
axs[1].set_xticks((0, generated_texture.shape[1]), labels=(0.0, 1.0))
axs[1].set_yticks((0, generated_texture.shape[0]), labels=(0.0, 1.0))
fig.show()

Create the napari layers

Next create two separate layers with the same mesh - once with each texture. In this example the texture coordinates happen to be the same for each texture, but this is not a strict requirement.

photo_texture_layer = napari.layers.Surface(
    (vertices, faces),
    texture=photo_texture,
    texcoords=texcoords,
    name='Texture_0',
)
generated_texture_layer = napari.layers.Surface(
    (vertices, faces),
    texture=generated_texture,
    texcoords=texcoords,
    name='GeneratedMat2',
)

Add the layers to a viewer

Finally, create the viewer and add the Surface layers. sphinx_gallery_thumbnail_number = 2

viewer = napari.Viewer(ndisplay=3)

viewer.add_layer(photo_texture_layer)
viewer.add_layer(generated_texture_layer)

viewer.camera.angles = (90.0, 0.0, -75.0)
viewer.camera.zoom = 75

if __name__ == '__main__':
    napari.run()