STCI 0.1.0__tar.gz

stci-0.1.0/Download_Euclid.py

@@ -0,0 +1,77 @@
+from __future__ import annotations
+
+from pathlib import Path
+
+import astropy.units as u
+import numpy as np
+from astropy.coordinates import SkyCoord
+from astropy.io import fits
+from astroquery.esa.euclid.core import EuclidClass
+
+
+BANDS = ("VIS", "NIR_Y", "NIR_J", "NIR_H")
+
+
+def _query(ra: float, dec: float) -> str:
+    return (
+        "SELECT file_name, file_path, filter_name, product_id "
+        "FROM dr1.mosaic_product "
+        "WHERE (product_type like '%Mer%Mosaic%') "
+        "AND ((instrument_name='VIS' AND filter_name='VIS') "
+        "OR (instrument_name='NISP' AND filter_name IN ('NIR_Y','NIR_J','NIR_H'))) "
+        "AND category='SCIENCE' AND fov IS NOT NULL AND "
+        f"INTERSECTS(CIRCLE('ICRS',{ra},{dec},0.0166667),fov)=1 "
+        "ORDER BY product_id ASC"
+    )
+
+
+def _is_all_zero(path: Path) -> bool:
+    with fits.open(path, memmap=False) as hdul:
+        data = next(hdu.data for hdu in hdul if hdu.data is not None)
+    return np.count_nonzero(np.nan_to_num(data, nan=0.0)) == 0
+
+
+def EUC_download(ra, dec, size, path, cred):
+    """
+    Download Euclid DR1 VIS/Y/J/H FITS cutouts.
+
+    Parameters
+    ----------
+    ra, dec : float
+        Target coordinates in degrees.
+    size : float
+        Cutout radius in arcsec.
+    path : str or Path
+        Output directory.
+    cred : str or Path
+        Euclid credentials file.
+    """
+
+    outdir = Path(path)
+    outdir.mkdir(parents=True, exist_ok=True)
+    coord = SkyCoord(float(ra), float(dec), unit="deg")
+
+    euclid = EuclidClass(environment="IDR")
+    euclid.login(credentials_file=str(cred))
+
+    rows = euclid.launch_job(_query(float(ra), float(dec)), verbose=False).get_results()
+    outputs = {}
+
+    for band in BANDS:
+        for i, row in enumerate(rows[rows["filter_name"] == band], start=1):
+            outfile = outdir / f"{band}_tile{i}.fits"
+            euclid.get_cutout(
+                file_path=f"{row['file_path']}/{row['file_name']}",
+                instrument="None",
+                id=str(row["product_id"]),
+                coordinate=coord,
+                radius=u.Quantity(float(size), u.arcsec),
+                output_file=str(outfile),
+            )
+            if outfile.stat().st_size > 0 and not _is_all_zero(outfile):
+                outputs[band] = outfile
+                break
+        if band not in outputs:
+            raise RuntimeError(f"No non-zero {band} cutout found.")
+
+    return outputs

stci-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Tian Li
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

stci-0.1.0/PKG-INFO

@@ -0,0 +1,110 @@
+Metadata-Version: 2.4
+Name: STCI
+Version: 0.1.0
+Summary: SpaceTelescopeColorImage tools for astronomy cutouts.
+Author: Tian Li
+License-Expression: MIT
+Project-URL: Homepage, https://github.com/astroskylee/STCI
+Project-URL: Repository, https://github.com/astroskylee/STCI
+Keywords: astronomy,euclid,fits,image-processing,color-image
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: astropy
+Requires-Dist: astroquery
+Requires-Dist: numpy
+Requires-Dist: pillow
+Requires-Dist: tifffile
+Dynamic: license-file
+
+# STCI
+
+STCI stands for **SpaceTelescopeColorImage**.
+
+STCI is a small astronomy image-composition tool for making display-ready
+color JPEGs from space telescope cutouts and other aligned mono images.
+
+## Installation
+
+```bash
+pip install STCI
+```
+
+The installed Python import module is `STCI`.
+
+## Example usage
+
+`mk_colorimg` creates one color JPEG from either a 3-channel RGB array or three mono images in `(R, G, B)` order. For Euclid-style color images, the default mapping is `NIR_J`, `NIR_Y`, `VIS`.
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(
+    [
+        "cutout_H.fits",    # R channel
+        "cutout_VIS.fits",  # G channel
+        "cutout_Y.fits",    # B channel
+    ],
+    output_jpg="target_mtf_vis_y_h.jpg",
+    input_mode="raw",
+)
+```
+
+For a NumPy RGB cube:
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(rgb_array, output_jpg="target_color.jpg", input_mode="normalized")
+```
+
+## Download a Euclid Color Image
+
+`Euclidimg` downloads Euclid DR1 `VIS`, `NIR_Y`, `NIR_J`, and `NIR_H` FITS cutouts, then renders one color JPEG using the `NIR_J / NIR_Y / VIS` channel order.
+
+```python
+from STCI import Euclidimg
+
+result = Euclidimg(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_color",
+    cred="Euclid/cred.txt",
+    output_jpg="EUCLJ032251.92-394609.8.jpg",
+)
+
+print(result["jpg"])
+print(result["fits"])
+```
+
+Arguments:
+
+- `ra`, `dec`: target coordinates in degrees.
+- `size`: cutout radius in arcsec. For example, `size=5.0` makes a `10" x 10"` image.
+- `path`: output directory for the FITS files and JPEG.
+- `cred`: Euclid credentials file passed to `astroquery.esa.euclid`.
+- `output_jpg`: optional JPEG filename written inside `path`.
+
+The returned dictionary contains the selected FITS paths and the final JPEG path. If the first overlapping mosaic tile for a band is empty or all zero, the downloader tries the next matching tile.
+
+## Download Euclid FITS Only
+
+```python
+from Download_Euclid import EUC_download
+
+fits_paths = EUC_download(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_fits",
+    cred="Euclid/cred.txt",
+)
+```

stci-0.1.0/README.md

@@ -0,0 +1,84 @@
+# STCI
+
+STCI stands for **SpaceTelescopeColorImage**.
+
+STCI is a small astronomy image-composition tool for making display-ready
+color JPEGs from space telescope cutouts and other aligned mono images.
+
+## Installation
+
+```bash
+pip install STCI
+```
+
+The installed Python import module is `STCI`.
+
+## Example usage
+
+`mk_colorimg` creates one color JPEG from either a 3-channel RGB array or three mono images in `(R, G, B)` order. For Euclid-style color images, the default mapping is `NIR_J`, `NIR_Y`, `VIS`.
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(
+    [
+        "cutout_H.fits",    # R channel
+        "cutout_VIS.fits",  # G channel
+        "cutout_Y.fits",    # B channel
+    ],
+    output_jpg="target_mtf_vis_y_h.jpg",
+    input_mode="raw",
+)
+```
+
+For a NumPy RGB cube:
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(rgb_array, output_jpg="target_color.jpg", input_mode="normalized")
+```
+
+## Download a Euclid Color Image
+
+`Euclidimg` downloads Euclid DR1 `VIS`, `NIR_Y`, `NIR_J`, and `NIR_H` FITS cutouts, then renders one color JPEG using the `NIR_J / NIR_Y / VIS` channel order.
+
+```python
+from STCI import Euclidimg
+
+result = Euclidimg(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_color",
+    cred="Euclid/cred.txt",
+    output_jpg="EUCLJ032251.92-394609.8.jpg",
+)
+
+print(result["jpg"])
+print(result["fits"])
+```
+
+Arguments:
+
+- `ra`, `dec`: target coordinates in degrees.
+- `size`: cutout radius in arcsec. For example, `size=5.0` makes a `10" x 10"` image.
+- `path`: output directory for the FITS files and JPEG.
+- `cred`: Euclid credentials file passed to `astroquery.esa.euclid`.
+- `output_jpg`: optional JPEG filename written inside `path`.
+
+The returned dictionary contains the selected FITS paths and the final JPEG path. If the first overlapping mosaic tile for a band is empty or all zero, the downloader tries the next matching tile.
+
+## Download Euclid FITS Only
+
+```python
+from Download_Euclid import EUC_download
+
+fits_paths = EUC_download(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_fits",
+    cred="Euclid/cred.txt",
+)
+```

stci-0.1.0/STCI.egg-info/PKG-INFO

@@ -0,0 +1,110 @@
+Metadata-Version: 2.4
+Name: STCI
+Version: 0.1.0
+Summary: SpaceTelescopeColorImage tools for astronomy cutouts.
+Author: Tian Li
+License-Expression: MIT
+Project-URL: Homepage, https://github.com/astroskylee/STCI
+Project-URL: Repository, https://github.com/astroskylee/STCI
+Keywords: astronomy,euclid,fits,image-processing,color-image
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: astropy
+Requires-Dist: astroquery
+Requires-Dist: numpy
+Requires-Dist: pillow
+Requires-Dist: tifffile
+Dynamic: license-file
+
+# STCI
+
+STCI stands for **SpaceTelescopeColorImage**.
+
+STCI is a small astronomy image-composition tool for making display-ready
+color JPEGs from space telescope cutouts and other aligned mono images.
+
+## Installation
+
+```bash
+pip install STCI
+```
+
+The installed Python import module is `STCI`.
+
+## Example usage
+
+`mk_colorimg` creates one color JPEG from either a 3-channel RGB array or three mono images in `(R, G, B)` order. For Euclid-style color images, the default mapping is `NIR_J`, `NIR_Y`, `VIS`.
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(
+    [
+        "cutout_H.fits",    # R channel
+        "cutout_VIS.fits",  # G channel
+        "cutout_Y.fits",    # B channel
+    ],
+    output_jpg="target_mtf_vis_y_h.jpg",
+    input_mode="raw",
+)
+```
+
+For a NumPy RGB cube:
+
+```python
+from STCI import mk_colorimg
+
+mk_colorimg(rgb_array, output_jpg="target_color.jpg", input_mode="normalized")
+```
+
+## Download a Euclid Color Image
+
+`Euclidimg` downloads Euclid DR1 `VIS`, `NIR_Y`, `NIR_J`, and `NIR_H` FITS cutouts, then renders one color JPEG using the `NIR_J / NIR_Y / VIS` channel order.
+
+```python
+from STCI import Euclidimg
+
+result = Euclidimg(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_color",
+    cred="Euclid/cred.txt",
+    output_jpg="EUCLJ032251.92-394609.8.jpg",
+)
+
+print(result["jpg"])
+print(result["fits"])
+```
+
+Arguments:
+
+- `ra`, `dec`: target coordinates in degrees.
+- `size`: cutout radius in arcsec. For example, `size=5.0` makes a `10" x 10"` image.
+- `path`: output directory for the FITS files and JPEG.
+- `cred`: Euclid credentials file passed to `astroquery.esa.euclid`.
+- `output_jpg`: optional JPEG filename written inside `path`.
+
+The returned dictionary contains the selected FITS paths and the final JPEG path. If the first overlapping mosaic tile for a band is empty or all zero, the downloader tries the next matching tile.
+
+## Download Euclid FITS Only
+
+```python
+from Download_Euclid import EUC_download
+
+fits_paths = EUC_download(
+    ra=50.7163333,
+    dec=-39.7693889,
+    size=5.0,
+    path="euclid_fits",
+    cred="Euclid/cred.txt",
+)
+```

stci-0.1.0/STCI.egg-info/SOURCES.txt

@@ -0,0 +1,11 @@
+Download_Euclid.py
+LICENSE
+README.md
+STCI.py
+Tian_color.py
+pyproject.toml
+STCI.egg-info/PKG-INFO
+STCI.egg-info/SOURCES.txt
+STCI.egg-info/dependency_links.txt
+STCI.egg-info/requires.txt
+STCI.egg-info/top_level.txt

stci-0.1.0/STCI.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

stci-0.1.0/STCI.egg-info/requires.txt

@@ -0,0 +1,5 @@
+astropy
+astroquery
+numpy
+pillow
+tifffile

stci-0.1.0/STCI.egg-info/top_level.txt

@@ -0,0 +1,3 @@
+Download_Euclid
+STCI
+Tian_color

stci-0.1.0/STCI.py

@@ -0,0 +1,21 @@
+"""Public entry point for STCI."""
+
+from pathlib import Path
+
+from Tian_color import *  # noqa: F401,F403
+
+
+def Euclidimg(ra, dec, size, path, cred, output_jpg="Euclid_color.jpg"):
+    """Download Euclid FITS cutouts and render one color JPEG."""
+
+    from Download_Euclid import EUC_download
+
+    outdir = Path(path)
+    fits_paths = EUC_download(ra, dec, size, outdir, cred)
+    jpg_path = outdir / output_jpg
+    mk_colorimg(
+        [fits_paths["NIR_J"], fits_paths["NIR_Y"], fits_paths["VIS"]],
+        output_jpg=jpg_path,
+        input_mode="raw",
+    )
+    return {"fits": fits_paths, "jpg": jpg_path}