Python package to perform resolved SED fitting using publically availabe SED fitting tools (Bagpipes and Dense Basis currently, Prospector, Beagle, etc planned), and display results.
- Interfaces with PIXEDFit, Voronoi binning, or single pixel binning to perform spatially resolved analysis.
- Portable file format - all information about a galaxy and SED fitting results are stored in an HDF5 file, allowing easy transfer, or the ability to run different parts of processing on different computers - e.g. create object on desktop, fit on a remote cluster and then analyse the results on a laptop.
- SED fitting using Bagpipes and Dense Basis with full support for all options, including custom priors, custom SFHs, dust laws, etc. Can also fit with a single pixel, or with a single SED for the whole galaxy. Multi-threading fitting is supported.
- Tools to model PSFs, either from stacking stars, or directly from WebbPSF. Can PSF homogenize internally, or accept pre-homogenized imaging.
- Interactive web-based viewer to display results of SED fitting, including interactive RGB imaging, segmentation map, binning maps, spatially resolved property maps (stellar mass density, SFR density etc), per bin/pixel SEDs, SFHs and corner plots, as well as interfacing with FITSMAP. Can also measure quantities directly from the photometry, such as Beta slopes, MUV, D4000 break etc and produce spatially resolved maps.
- Testing with Simulated galaxies from hydro sims, or generated parametrically, using the Synthesizer package. Allows comparison of recovered parameters like stellar mass, dust, SFH etc with true values. Works with FLARES, SPHINX20 and other sims which can be loaded via the Synthesizer package.
- Morphological fitting with pysersic, including multi-component models and complex priors. Pyautogalaxy and pygalfitm are partially supported.
- Radial/annular analysis of property maps, including SED fitting, radial profiles of properties, and annular SED fitting.
- Direct conversion of cutouts with modelled morphology into components for Pandeia - e.g. to plan JWST observations.
- More interactivity in the viewer - e.g. to select regions for further analysis - e.g. to measure the SED of a region, or to fit a model to it, based on drawing an aperture, or grouping pixels (implemented for EAZY SED Fitting.)
- Additional SED fitting tools, and improvement of pixel-by-pixel case. Prospector is in progress, and Beagle is planned.
- Better tools for MetaClass for population analysis - spatially resolved main sequence, outshining etc.
- Furhter reionization estimates - maps of fesc, zi_ion etc
- More radial/annular tools - annular SED fitting, radial profiles of properties etc.
- Bug fixes, documentation and examples.
Here are a couple of examples of some of the GUI features, including an interface to plot and compare the SED fitting results - SED plots, resolved parameter maps, SFH, corner plots.
This example shows an example of the interactive interface for EAZY SED fitting, which lets you place manual apertures and regions and fit their photometry.
To install, clone the repository and run pip install . in the root directory. For best results, also install our fork of Bagpipes from here.
There is a notebook (examples/intro.ipynb) which goes over the basic features of EXPANSE from the beginning of loading in the data. All the plots for the EXPANSE paper are available in scripts/outshining_paper_figure.ipynb, which show more advanced usage. The examples/documentation are a work in progress.
Please get in touch if you're interested in using the code, I'm happy to chat!
To lauch the viewer, run expanse-viewer in the terminal after installation.
The paper describing EXPANSE is currently in preparation. In the meantime, please cite this repository if you use it in your work.