dustmaps3d 2.1.10

An all-sky 3D dust map based on Gaia and LAMOST.

Read this in: English | 中文

dustmaps3d

🌌 An all-sky 3D dust extinction map based on Gaia and LAMOST

📄 Wang et al. (2025), An all-sky 3D dust map based on Gaia and LAMOST
📌 DOI: 10.12149/101620

📦 A Python package for easy access to the 3D dust map
📌 DOI: 10.12149/101619

---

📦 Installation

Install via pip:

pip install dustmaps3d

Note: Installing the package does not include the data file.
The ~350 MB model data will be automatically downloaded from GitHub on first use.
⚠️ If you experience network issues when downloading from GitHub,
you can manually download the data from NADC:
🔗 https://nadc.china-vo.org/res/r101619/

---

🚀 Usage

from dustmaps3d import dustmaps3d

l = [120.0]
b = [30.0]
d = [1.5]

EBV, dust, sigma, max_d = dustmaps3d(l, b, d)

print(f"EBV: {EBV[0]:.4f} [mag]")
print(f"Dust: {dust[0]:.4f} [mag/kpc]")
print(f"Sigma: {sigma[0]:.4f} [mag]")
print(f"Max distance: {max_d.iloc[0]:.4f} kpc")

Batch example with FITS:

import numpy as np
from astropy.table import Table
from dustmaps3d import dustmaps3d

data = Table.read('input.fits')   
l = data['l'].astype(float)
b = data['b'].astype(float)
d = data['distance'].astype(float)

EBV, dust, sigma, max_d = dustmaps3d(l, b, d)

data['EBV_3d'] = EBV
data['dust'] = dust
data['sigma'] = sigma
data['max_distance'] = max_d
data.write('output.fits', overwrite=True)

---

🧠 Function Description

dustmaps3d(l, b, d, n_process=None)

Estimates 3D dust extinction and related quantities for given galactic coordinates and distances.

Input	Type	Description	Unit
l	np.ndarray	Galactic longitude	degrees
b	np.ndarray	Galactic latitude	degrees
d	np.ndarray	Distance	kpc

Returns:

Output	Type	Description	Unit
EBV	np.ndarray	E(B–V) extinction	mag
dust	np.ndarray	Dust density (d(EBV)/dx)	mag/kpc
sigma	np.ndarray	Estimated uncertainty in E(B–V)	mag
max_d	np.ndarray	Maximum reliable distance	kpc

If d contains NaN, it will be automatically replaced by the maximum reliable distance along that line of sight (max_d).

If the input d exceeds max_d, it indicates the point lies beyond the model's reliable range. The returned values in this case are extrapolated and not guaranteed to be accurate.

---

⚡ Performance

• Fully vectorized and optimized with NumPy
• On a modern personal computer, evaluating 100 million stars takes only ~100 seconds

---

📂 Data Version

This version uses data_v2.1.parquet, released under v2.1

---

📜 Citation

If you use this model or the Python package, please cite both:

• Wang, T. (2025), An all-sky 3D dust map based on Gaia and LAMOST. DOI: 10.12149/101620
• dustmaps3d: A Python package for easy access to the 3D dust map. DOI: 10.12149/101619

---

🛠️ License

This project is open-source and distributed under the MIT License.

---

📫 Contact

If you have any questions, suggestions, or encounter issues using this package,
please feel free to contact the authors via GitHub issues or email.

• Prof. Yuan Haibo: [e-mail]
• Wang Tao: [e-mail]

🔗 GitHub Repository