DECaPS 3D Dust Map

<h1>DECaPS 3D Dust Map</h1>

<p>A three-dimensional map of dust reddening, covering the Southern Galactic plane (<code>239° < l < 6°, |b| < 10°</code>).</p>

<p>We use DECaPS optical and near-infrared photometry, in conjunction with complementary photometric and astrometric data, to infer distances and reddenings to 709 million stars. These stars trace the reddening along different lines of sight, allowing us to build up a map of reddening in 3D.</p>

<p>The map is structured as a set of sightlines, each of which contains multiple samples of the cumulative dust reddening as a function of distance. Each sightline is identified by a nested pixel index at a HEALPix <code>N<sub>side</sub> = 8192</code>. Within each sightline, cumulative reddening is given at discrete distances, spaced evenly in distance modulus. </p>

<h2>Quality Assurance Information</h2>
<p>Quality assurance information is given for each pixel, including:</p>
<ul>
<li>Whether the fit converged in the pixel</li>
<li>Whether the pixel needed to be infilled due to lack of stars</li>
<li>The minimum reliable distance modulus in the pixel</li>
<li>The maximum reliable distance moduli in the pixel</li>
</ul>

<p>Alongside the quality assurance information, we also provide the number of stars whose PDF on distance and reddening is used to inform the line of sight fit towards each pixel. Note that this will not sum to the total number of stars used in the construction of the map (709 million) because some stars contribute to multiple pixels based on our Gaussian weighting scheme. </p>

<p>Unlike the "Bayestar" 3D dust map from Green et al. 2019 (whose reddening is provided in an arbitrary reddening unit), the DECaPS 3D dust reddening is given in units of <code>E(B-V)</code> in mags. When combined with Bayestar, the DECaPS map enables extinction corrections over the entire Galactic plane <code>|b| < 10°</code>.</p>

<p>The 3D map is described in more detail in <strong>Zucker, Saydjari, & Speagle et al. 2025.</p></strong>

<h2>Downloading and Querying the map</h2>
<p>This map is included in the Python <strong><a href="https://dustmaps.readthedocs.io/en/latest/" target="_blank">dustmaps</a></strong> package, which is available through <code>pip</code>:</p>

<pre><code>pip install dustmaps</code></pre>

<p>Note that we generate 100 samples of the reddening. In the data product <strong><code>decaps_mean_and_samples.h5</code></strong>, we calculate the mean from the 100 samples and additionally provide five random samples, all stored as <code>float16</code> to reduce memory usage (30 GB in size). For those with limited disk space, in <strong><code>decaps_mean.h5</code></strong>, we provide the option to only download the mean map (no samples), which is significantly smaller in size (7 GB).</p>