| Commit message (Collapse) | Author | Age | Files | Lines |
... | |
| |
|
|
|
|
| |
Also update the general config for the new "galactic/snr" component.
|
|
|
|
|
|
|
|
| |
The necessary but missing functionalities to simulate the Galactic SNRs
emission maps are implemented, and this new emission component is ready
for testing.
Also fix a typo in "utils/grid.py".
|
|
|
|
|
|
|
|
| |
The "map_grid_to_healpix()" maps the generated coordinate grid to a
HEALPix map. Note that only effective HEALPix pixels are returned
instead of a full HEALPix map.
TODO: Try to avoid the explicit for loop to optimize the speed.
|
|
|
|
|
|
|
| |
* make_coordinate_grid(): Make a Cartesian coordinate grid according to
the specified center coordinate, size and grid resolution.
* make_grid_ellipse(): Also generate an (rotated) ellipse and place into
the coordinate grid.
|
|
|
|
| |
NOTE: very preliminary, and missing many necessary functions!
|
|
|
|
| |
NOTE: very preliminary and may change at any time!
|
| |
|
| |
|
|
|
|
|
|
|
|
| |
Change the "__init__()" method to load the user configurations (i.e., by
invoke the "read_userconfig()" method) instead of loading extra general
configurations.
Method "getn()": also check the validity for the parameter "sep".
|
|
|
|
| |
No functional changes.
|
| |
|
| |
|
| |
|
|
|
|
|
|
| |
* Return None if specified config is None or not exist
* Raise ValueError if specified config is non-string
* Update comments
|
| |
|
|
|
|
|
|
|
|
|
| |
These functions will be used to draw the image of a SNR or a radio halo.
These functions are copied from the `scikit-image` project.
Thanks!
scikit-image: http://scikit-image.org/docs/dev/api/skimage.draw.html
|
|
|
|
|
|
|
|
|
| |
After this switch to `pkg_resources`, this package can also be installed
as an EGG (a zip archive), which is more generic.
References:
* http://setuptools.readthedocs.io/en/latest/pkg_resources.html#basic-resource-access
* http://setuptools.readthedocs.io/en/latest/setuptools.html#including-data-files
|
|
|
|
|
|
| |
Thanks D. A. Green for the Galactic SNRs catalog (294 objects)!
http://www.mrao.cam.ac.uk/surveys/snrs/
|
| |
|
|
|
|
|
|
| |
Since the logging is not avaiable at the moment, we can only use the
plain `print` for debugging, if the environment variable "DEBUG_FG21SIM"
is defined.
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* Mask the regions where the dust absorption is too high to well
determine the true Halpha absorption.
There are regions (e.g., Galactic plane) have too large dust
absorption which cause float overflow when applying absorption
correction to the Halpha map.
* Fix a bug on the wrong parameter passed to "_calc_ratio_a()" method
in "_simulate_frequency()"
* Move "_correct_dust_absorption()" invokation from "simulate()"
to "_simulate_frequency()"
Already tested and OK!
|
|
|
|
|
|
|
|
| |
* healpix2hpx(), hpx2healpix(): Remove the "header" parameter, thus the
"data" parameter can only be either the filename or a HDU;
* healpix2hpx(): Remove the check on "ORDERING", since
"read_fits_healpix()" always return the HEALPix data in RING ordering;
* Small updates to the log messages and comments.
|
| |
|
|
|
|
|
|
|
|
|
| |
* bin/fg21sim: check the "common/components" config to simulate the
enabled components;
* bin/fg21sim: add Galactic free-free component support;
* galactic/freefree.py: fix missing attributes.
NOTE: simulated free-free maps are WRONG!! Need check/debug!
|
|
|
|
|
|
|
| |
* New class "galactic.FreeFree" to simulate free-free emission
* Add new config section "galactic/freefree"
NOTE: current untested
|
|
|
|
|
|
|
| |
* Merge the upgrade/downgrade resolution w.r.t Nside to the map loading
methods;
* Update the variable names;
* Check the status of smallscales addition.
|
| |
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
| |
The missing pixels in the reprojected HEALPix map are filled with the
averages of their 8 neighboring pixels (excluding the NaN's if any).
Also add the "--inpaint" argument to the executable script.
|
|
|
|
| |
Also fix typos
|
|
|
|
|
|
|
| |
The mask array has values 0, 1, and 2, which indicate the different
statuses of the reprojected HEALPix pixels.
Also update the executable script to use the new return results.
|
| |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
The "zea2healpix()" function reprojects the maps in ZEA
(zenithal/azithumal equal area) projection to HEALPix data in Galactic
frame with RING ordering.
The other two helper functions "_image_to_healpix()" and "_convert_wcs()"
are almost copied from the "reproject" project [1]. Thanks!
TODO:
* Add some more logging
* Implement the "inpaint" argument to inpaint the HEALPix map if exists
missing pixels
[1] reproject: https://github.com/astrofrog/reproject
|
|
|
|
|
|
|
|
| |
* utils/fits.py: hack the dtype to ignore the byteorder (FITS data use
big endianness, e.g., dtype(">f4"))
* utils/healpix.py: explicit convert the dtype and log the dtype
* bin/healpix2hpx, bin/hpx2healpix: remove the --float argument
* other minor fixes/updates
|
|
|
|
| |
Also remove a unused variable in "galactic/synchrotron.py"
|
| |
|
|
|
|
|
|
| |
This function wraps on the `healpy.read_map()`, but reset the data array
to its original dtype in FITS file, as well as return the FITS header in
`astropy.io.fits.Header` object.
|
| |
|
| |
|
|
|
|
|
|
| |
NOTE:
This "write_fits_healpix()" function only implements the most common
case of the HEALPix FITS style.
|
|
|
|
|
|
|
|
|
|
|
|
| |
* Add new executable "bin/fg21sim"
* galactic/synchrotron: update to use "configs.get_path()"
* galactic/synchrotron: create output dir if not exists
* galactic/synchrotron: add logging support
* galactic/synchrotron: append FITS extension to filename
* galactic/synchrotron: pass the basic test
TODO:
* "output()" needs fixes with the FITS header
|
| |
|
| |
|
| |
|