Metadata-Version: 2.1
Name: ps-eor
Version: 0.19
Summary: Foreground modeling/removal and Power Spectra generation
Home-page: https://gitlab.com/flomertens/ps_eor
License: GPL-3.0+
Author: Florent Mertens
Author-email: flomertens@gmail.com
Requires-Python: >=3.7,<4.0
Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
Classifier: Programming Language :: Python :: 3
Classifier: Programming Language :: Python :: 3.10
Classifier: Programming Language :: Python :: 3.7
Classifier: Programming Language :: Python :: 3.8
Classifier: Programming Language :: Python :: 3.9
Requires-Dist: GPy (>=1.9,<2.0)
Requires-Dist: astropy (>=4.0,<6.0)
Requires-Dist: click (>=7.0,<9.0)
Requires-Dist: configparser (>=4.0,<5.0)
Requires-Dist: corner (>=2.0,<3.0)
Requires-Dist: dynesty (>=1.2.2,<2.0.0)
Requires-Dist: emcee (>=3.0.2,<4.0.0)
Requires-Dist: fast-histogram (>=0.9,<0.10)
Requires-Dist: h5py (>=3.0,<4.0)
Requires-Dist: healpy (>=1.12,<2.0)
Requires-Dist: importlib-metadata (>=1.0,<2.0); python_version < "3.8"
Requires-Dist: joblib (>=1.0,<2.0)
Requires-Dist: libpipe (>=0.1.5)
Requires-Dist: matplotlib (>=3,<4)
Requires-Dist: pyfftw (>=0.12)
Requires-Dist: reproject (>=0.7,<0.8)
Requires-Dist: scikit-learn (>=0.20)
Requires-Dist: scipy (>=1.4,<2.0)
Requires-Dist: tables (>=3.2,<4.0)
Requires-Dist: torch (>=1.9.0,<2.0.0)
Requires-Dist: ultranest (>=3.0.0,<4.0.0)
Project-URL: Repository, https://gitlab.com/flomertens/ps_eor
Description-Content-Type: text/markdown

Python package for foreground removal and power spectra calculation:
===================================================================

Installation
------------

From git (development version):

    git clone https://gitlab.com/flomertens/ps_eor.git
    pip install -f ./ps_eor

from [pypi](https://pypi.org/project/ps-eor/) (latest stable version):

    pip install ps_eor

The following packages will also be installed:

- click
- numpy
- scipy
- astropy
- matplotlib
- pytable (for saving loading in h5 format)
- pyfftw (optional, can save time)
- sklearn (for PCA)
- GPy (for GPR)
- healpy,reproject (optional)


Command-line usage:
-------------------

The ps_eor package comes with a command line tool: pstool, which can:
- Generate gridded visibility cube from fits images cubes (gen_vus_cube)
- Produce power-spectra from gridded visibility cube (make_ps)
- Combine different visibility cube (combine)
- Remove foregrounds using GPR (run_gpr)
- ... and several more

Usage of pstool is as follow:

    Usage: pstool [OPTIONS] COMMAND [ARGS]...

    Options:
      --version  Show the version and exit.
      --help     Show this message and exit.

    Commands:
      combine               Combine all datacubes listed in FILE_LIST...
      combine_sph           Combine all sph datacubes listed in FILE_LIST...
      diff_cube             Compute the difference between CUBE1 and CUBE2
      even_odd_to_sum_diff  Create SUM / DIFF datacubes from EVEN / ODD...
      gen_vis_cube          Create a datacube from image and psf fits files.
      make_ps               Produce power-spectra of datacubes FILE_I: Input...
      run_flagger           Run flagger on datacubes and save flag.
      run_gpr               Run GPR & generate power spectra FILE_I: Input...
      simu_noise_img        Compute noise fits cube from a simulated UV
                            coverage...

      simu_noise_ps         Compute noise PS from a simulated UV coverage...
      simu_uv               Compute gridded UV coverage.
      vis_to_sph            Load visibilities datacubes VIS_CUBE, transform to...


To get help for any of the pstool command, execute:

    pstool CMD --help


API usage:
----------

This package provide a simple interface for Foreground removal and PS estimation tasks which can complete the command line tool.

Below we briefly document the main tasks that can be perform with this package:


1. Loading your data: 
---------------------

    data_cube_i = datacube.CartDataCube.load_from_fits_image(files, umin, umax, theta_fov)

This will do the following steps:
- Read Fits image file
- Trim image to theta_fov
- Convert image from Jy/PSF to K, using imager_scale_factor or WSCNORMF
 attribute to get PSF "solid angle" (otherwise use Gaussian approx of the PSF)
- FFT image per frequencies to get visibilities
- Keep only non-zero visibilities between umin and umax.
and return a CartDataCube object storing the visibilities in an ungridded way.

It is possible (recommended) to save/load to an h5 format with the save()/load() method.

One can also regrid the data and make an image with the regrid() and image() method.


2. Run FG removal algorithm:
----------------------------

The main FG removal code is GPR, but PCA, GMCA (the python version) and Polynamial fitting are also implemented.

To run GPR, one do the following:

    data_cube_i = datacube.CartDataCube.load_from_fits_image(files_i, umin, umax, theta_fov)
    data_cube_v = datacube.CartDataCube.load_from_fits_image(files_v, umin, umax, theta_fov)

    eor_bin_list = pspec.EorBinList(data_cube_i.freqs)
    # Create an EoR bin 122-134 MHz with a 120-136MHz range for the FG fitting
    eor_bin_list.add_freq(1, 122, 134, 120, 136)
    eor = eor_bin_list.get(1)

    gpr_config = fitutil.GprConfig.load(gpr_config_filename)
    gpr_fit = fgfit.GprForegroundFit(gpr_config)
    gpr_res = gpr_fit.run(eor.get_slice_fg(data_cube_i), eor.get_slice_fg(data_cube_v))

This return a GprForegroundResult object which have the following attributes:
- fit: The FG model in a form of a CartDataCube object
- sub: The residual 
- pre_fit: The pre-fit FG model
- post_fit: The post-fit FG model

And the following method:
- get_fg_model(): return the GPR fg model
- get_eor_model(): return the GPR eor model

On can then save/load those CartDataCube as needed for later processing.

The CartDataCube cubes of the GPR model and residual contains the error covariance 
from the GPR model that need to be taken into account when generating the 
power spectra.

Look at gpr_config.parset and gpr_config_v.parset for examples of GPR configuration.


3. Generate Power Spectra:
--------------------------

The PS code into account automatically the error covariance of the GPR model.

It is possible to generate spatial only PS, Cylindrically averaged PS (2D) or spherically averaged PS (3D).

    # Create a PS configuration
    el = 2 * np.pi * (np.arange(data_cube.ru.min(), data_cube.ru.max(), du))
    ps_conf = pspec.PowerSpectraConfig(el)

    pb = datacube.LofarHBAPrimaryBeam()
    
    # Create a PS generation object
    ps_gen = pspec.PowerSpectraCart(eor, ps_conf, pb)

    # Create a Spatial PS, plot it and save it to a file
    ps = ps_gen.get_ps(data_cube)
    ps.plot(title='Spatial power spectra')
    plt.savefig('ps.pdf')
    ps.save_to_txt('ps.txt')

    # Create a Cylindrically averaged PS
    ps2d = ps_gen.get_ps2d(data_cube)
    ps2d.plot(title='Cylindrically averaged power spectra')
    plt.savefig('ps2d.pdf')
    ps2d.save_to_txt('ps2d.txt')

    # Create a Spherically averaged PS
    kbins = np.logspace(np.log10(ps_gen.kmin), np.log10(0.5), 10)
    ps3d = ps_gen.get_ps3d(kbins, data_cube)
    ps3d.plot(title='Spherically averaged power spectra')
    plt.savefig('ps3d.pdf')
    ps3d.save_to_txt('ps3d.txt')


