Skip to main content
This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (
Help us improve Python packaging - Donate today!

Uniform interface for multiple dust reddening maps.

Project Description

A uniform interface for a number of 2D and 3D maps of interstellar dust reddening/extinction.

Supported Dust Maps

  1. Burstein & Heiles (1982; BH‘82)
  2. Chen et al. (2014)
  3. Green, Schlafly, Finbeiner et al. (2015,2018; Bayestar)
  4. Marshall et al. (2006)
  5. Planck Collaboration (2013)
  6. Sale et al. (2014; IPHAS)
  7. Schlegel, Finkbeiner & Davis (1998; SFD‘98)


Download the repository and then run:

python install –large-data-dir=/path/where/you/want/large/data/files/stored

Alternatively, you can use pip:

pip install dustmaps

Getting the Data

To fetch the data for the SFD dust map, run:

python fetch –map-name=sfd

You can download the other dust maps by changing “sfd” to “planck”, “bayestar”, “iphas”, “marshall”, “chen2014” or “bh”.

Alternatively, if you have used pip to install dustmaps, then you can configure the data directory and download the data by opening up a python interpreter and running:

>>> from dustmaps.config import config
>>> config['data_dir'] = '/path/where/you/want/large/data/files/stored'
>>> import dustmaps.sfd
>>> dustmaps.sfd.fetch()
>>> import dustmaps.planck
>>> dustmaps.planck.fetch()
>>> import dustmaps.bayestar
>>> dustmaps.bayestar.fetch()
>>> import dustmaps.iphas
>>> dustmaps.iphas.fetch()
>>> import dustmaps.marshall
>>> dustmaps.marshall.fetch()
>>> import dustmaps.chen2014
>>> dustmaps.chen2014.fetch()

Querying the Maps

Maps are queried using astropy.coordinate.SkyCoord objects. This means that any coordinate system can be used as input. For example, we can query SFD‘98 as follows:

>>> from dustmaps.sfd import SFDQuery
>>> from astropy.coordinates import SkyCoord
>>> sfd = SFDQuery()
>>> c = SkyCoord(
>>> print sfd(c)

Above, we’ve used the ICRS coordinate system (the inputs are RA and Dec). We can use other coordinate systems, such as Galactic coordinates, and we can provide coordinate arrays. The following example uses both features:

>>> c = SkyCoord(
        [75.00000000, 130.00000000],
        [-89.00000000, 10.00000000],
>>> print sfd(c)
[ 0.0146584   0.97695869]


Read the full documentation at

Release History

This version
History Node


History Node


History Node


History Node


History Node


History Node


History Node


History Node


History Node


History Node


History Node


Download Files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

File Name & Hash SHA256 Hash Help Version File Type Upload Date
(419.3 kB) Copy SHA256 Hash SHA256
Source Jan 14, 2018

Supported By

WebFaction WebFaction Technical Writing Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Rackspace Rackspace Cloud Servers DreamHost DreamHost Log Hosting