| Commit message (Collapse) | Author | Age | Files | Lines |
... | |
| |
|
|
|
|
|
|
|
|
|
|
| |
NOTE: However, there are two significant problems:
* Cannot figure out why ``nb.optional(nb.types.UniTuple(nb.int64, 2))``
does NOT work. Therefore, make ``shape`` as mandatory parameter
instead of optional.
* Cannot figure out multi-dispatch that allows both int and float types
for ``r``, ``c``, ``r_radius`` and ``c_radius``. Thus only support
the int type for the moment.
|
| |
|
|
|
|
|
| |
Remove the original "make_coordinate_grid()" function, and rename the
"make_coordinate_grid_fast()" to override the former one.
|
|
|
|
|
| |
Change the parameters of "make_coordinate_grid_fast()", make it share
the same parameter list as "make_coordinate_grid()".
|
|
|
|
|
| |
There are some limitations with the numba JIT, which are commented in
the code.
|
|
|
|
| |
Also fix a wrong notation in "random.py"
|
|
|
|
|
|
|
|
|
|
| |
To optimize "map_grid_to_healpix()" in `grid.py` which uses
"healpy.ang2pix()" using `numba.jit`, implement the latter's
functionality explicitly with `numba.jit` support.
The implementation simply mimic the corresponding functions
"ang2pix_ring_z_phi()" and "pix2ang_ring_z_phi()" in HEALPix's
"src/C/subs/chealpix.c". Thanks!
|
|
|
|
|
| |
Also place "_calc_hpx_row_idx()" before "_calc_hpx_indicies()", which is
required by `numba`.
|
|
|
|
|
|
|
|
| |
Optimize "make_coordinate_grid()" using "numba.jit", which improves the
speed a bit, since the other two functions (i.e., "make_grid_ellipse()",
and "map_grid_to_healpix()") are much slow and time-consuming.
TODO: optimize "make_grid_ellipse()" and "map_grid_to_healpix()"
|
|
|
|
|
|
|
|
|
| |
* Optimize the "Fnu_to_Tb()" and "Sb_to_Tb()" functions by explicitly
calculating the conversions, avoiding the slow `astropy.unit`
conversions (which is rather slow). The new fast functions are named
as "Fnu_to_Tb_fast()" and "Sb_to_Tb_fast()".
* Optimize the new functions using `numba.jit` further.
* Add `numba` to the dependencies.
|
|
|
|
|
| |
Currently implement "spherical_uniform()" to generate random points
uniformly distributed on the spherical surface.
|
|
|
|
|
|
|
|
| |
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.
|
| |
|
| |
|
| |
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
* 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.
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
| |
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.
|
| |
|
|
|
|
|
|
|
| |
By saving the existing formatter and setting for the new handlers, the
configured format styles are kept. Otherwise, the newly created
handlers does not have a formatter, thus the default format style gets
used.
|
|
|
|
| |
Also improve the comments.
|
|
|
|
|
| |
"logging.basicConfig()" does NOT accept parameter "filemode" if the
corresponding paramter "filename" NOT specified.
|
| |
|
|
|
|
| |
Also fix a bug about parameter "stream".
|
|
|
|
| |
Also fix a bug about "append_comment".
|
|
|
|
|
|
| |
This "setup_logging()" utility can setup the logging module with the
given configuration dict, and can also update the logging configs by
specify the additional parameters.
|
| |
|
|
* healpix2hpx: reorganize HEALPix data in HPX projection
* hpx2healpix: recover the HEALPix data from HPX projection
Thanks: M. Calabretta: WCSLIB: HPXcvt
TODO: Add test codes
|