This code was used to create the 3D meshes from Gaia DR2 as announced in

A map of the Milky Way out to 3000 parsecs (10 thousand light years)

It is documented across several blog posts linked below.

Making a galaxy map, Part 1: Data

dr2_make_star_list_0.1.py

Making a galaxy map, Part 2: Slicing

dr2_make_slices_ob.py

Making a galaxy map, Part 3: Meshing

dr2_make_isosurfaces_all.py

Making a galaxy map, Part 4: Enclosure

dr_region_structure.py

dr2_make_iso_file.py

dr2_peaks_to_json.py

dr2_make_ranges.py

dr2_list_ionizing_ranges.py

dr2_make_mesh_sets.py

Making a galaxy map, Part 5: Overlays

Making a galaxy map, Part 6: Resources

Adding these related blog post references for completeness.