aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/utils
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim/utils')
-rw-r--r--fg21sim/utils/healpix.py34
-rw-r--r--fg21sim/utils/random.py4
2 files changed, 19 insertions, 19 deletions
diff --git a/fg21sim/utils/healpix.py b/fg21sim/utils/healpix.py
index fe436d4..f72a8a6 100644
--- a/fg21sim/utils/healpix.py
+++ b/fg21sim/utils/healpix.py
@@ -71,9 +71,9 @@ def healpix2hpx(data, append_history=None, append_comment=None):
logger.info("Loaded HEALPix data: dtype={0}, Npixel={1}, Nside={2}".format(
dtype, npix, nside))
hp_data = np.append(hp_data, np.nan).astype(dtype)
- logger.info("Calculating the HPX indices ...")
- hpx_idx = _calc_hpx_indices(nside)
- # Fix indices of "-1" to set empty pixels with above appended NaN
+ logger.info("Calculating the HPX indexes ...")
+ hpx_idx = _calc_hpx_indexes(nside)
+ # Fix indexes of "-1" to set empty pixels with above appended NaN
hpx_idx[hpx_idx == -1] = len(hp_data) - 1
hpx_data = hp_data[hpx_idx]
hpx_header = _make_hpx_header(hp_header,
@@ -126,11 +126,11 @@ def hpx2healpix(data, append_history=None, append_comment=None):
logger.info("Determined HEALPix Nside=%d" % nside)
#
npix = hp.nside2npix(nside)
- logger.info("Calculating the HPX indices ...")
- hpx_idx = _calc_hpx_indices(nside).flatten()
+ logger.info("Calculating the HPX indexes ...")
+ hpx_idx = _calc_hpx_indexes(nside).flatten()
hpx_idx_uniq, idxx = np.unique(hpx_idx, return_index=True)
if np.sum(hpx_idx_uniq >= 0) != npix:
- raise ValueError("Number of pixels does not match indices")
+ raise ValueError("Number of pixels does not match indexes")
hpx_data = hpx_data.flatten()
hp_data = hpx_data[idxx[hpx_idx_uniq >= 0]]
hp_header = _make_healpix_header(hpx_header, nside=nside,
@@ -141,7 +141,7 @@ def hpx2healpix(data, append_history=None, append_comment=None):
@nb.jit(nb.int64[:](nb.int64, nb.int64, nb.int64), nopython=True)
def _calc_hpx_row_idx(nside, facet, jmap):
- """Calculate the HEALPix indices for one row of a facet.
+ """Calculate the HEALPix indexes for one row of a facet.
NOTE
----
@@ -210,8 +210,8 @@ def _calc_hpx_row_idx(nside, facet, jmap):
@nb.jit(nb.int64[:, :](nb.int64), nopython=True)
-def _calc_hpx_indices(nside):
- """Calculate HEALPix element indices for the HPX projection scheme.
+def _calc_hpx_indexes(nside):
+ """Calculate HEALPix element indexes for the HPX projection scheme.
Parameters
----------
@@ -220,15 +220,15 @@ def _calc_hpx_indices(nside):
Returns
-------
- indices : 2D `~numpy.ndarray`
+ indexes : 2D `~numpy.ndarray`
2D integer array of same size as the input/output HPX FITS image,
- with elements tracking the indices of the HPX pixels in the
+ with elements tracking the indexes of the HPX pixels in the
HEALPix 1D array, while elements with value "-1" indicating
null/empty HPX pixels.
NOTE
----
- * The indices are 0-based;
+ * The indexes are 0-based;
* Currently only HEALPix RING ordering supported;
* The null/empty elements in the HPX projection are filled with "-1".
"""
@@ -250,7 +250,7 @@ def _calc_hpx_indices(nside):
FACETS_LAYOUT[4, :] = [-1, -1, -1, 2, 6]
#
shape = (nfacet*nside, nfacet*nside)
- indices = -np.ones(shape, dtype=np.int64)
+ indexes = -np.ones(shape, dtype=np.int64)
#
# Loop vertically facet-by-facet
for jfacet in range(nfacet):
@@ -266,9 +266,9 @@ def _calc_hpx_indices(nside):
else:
idx = _calc_hpx_row_idx(nside, facet, j)
col = ifacet * nside
- indices[row, col:(col+nside)] = idx
+ indexes[row, col:(col+nside)] = idx
#
- return indices
+ return indexes
def _make_hpx_header(header, append_history=None, append_comment=None):
@@ -510,7 +510,7 @@ def ang2pix_ring(nside, theta, phi):
Returns
-------
ipix : 1D or 1D `~numpy.ndarray`
- The indices of the pixels corresponding to the input coordinates.
+ The indexes of the pixels corresponding to the input coordinates.
The shape is the same as the input array.
NOTE
@@ -539,7 +539,7 @@ def pix2ang_ring(nside, ipix):
Parameters
----------
ipix : 1D or 2D `~numpy.ndarray`
- The indices of the HEALPix pixels in the RING ordering
+ The indexes of the HEALPix pixels in the RING ordering
Returns
-------
diff --git a/fg21sim/utils/random.py b/fg21sim/utils/random.py
index f652f6d..d825101 100644
--- a/fg21sim/utils/random.py
+++ b/fg21sim/utils/random.py
@@ -20,11 +20,11 @@ def spherical_uniform(n=1):
Returns
-------
theta : float, or 1D `~numpy.ndarray`
- The polar angles, φ ∈ [0, π]. (unit: rad)
+ The polar angles, θ ∈ [0, π]. (unit: rad)
If ``n > 1``, then returns a 1D array containing all the generated
coordinates. (unit: rad)
phi : float, or 1D `~numpy.ndarray`
- The azimuthal angles, θ ∈ [0, 2π).
+ The azimuthal angles, φ ∈ [0, 2π).
NOTE
----