lightstack 0.1.0__tar.gz → 0.1.2__tar.gz

lightstack-0.1.2/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Andressa Wille, Thallis Pessi
+
+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.

lightstack-0.1.2/PKG-INFO

@@ -0,0 +1,64 @@
+Metadata-Version: 2.1
+Name: lightstack
+Version: 0.1.2
+Summary: Tools for building and processing multi-filter astrophysical datacubes
+Author: Andressa Wille, Thallis Pessi
+License: MIT License
+        
+        Copyright (c) 2026 Andressa Wille, Thallis Pessi
+        
+        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.
+        
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy
+Requires-Dist: astropy
+Requires-Dist: matplotlib
+Requires-Dist: reproject
+Requires-Dist: regions
+Requires-Dist: scipy
+Requires-Dist: photutils
+
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.19393439.svg)](https://doi.org/10.5281/zenodo.19393439)
+
+<img src="https://raw.githubusercontent.com/AndressaWille/lightstack/master/images/logo.png" width="120">
+
+# Lightstack
+
+Tools for building and processing multi-filter astrophysical datacubes.
+
+## Features
+
+- Cropping FITS images  
+- Building photometric datacubes  
+- PSF matching  
+
+(Currently supports JWST and HST)
+
+## Installation
+
+```bash
+pip install lightstack
+```
+
+## Documentation
+
+https://lightstack.readthedocs.io/en/latest/
+
+See the tutorial notebook!

lightstack-0.1.2/README.md

@@ -0,0 +1,27 @@
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.19393439.svg)](https://doi.org/10.5281/zenodo.19393439)
+
+<img src="https://raw.githubusercontent.com/AndressaWille/lightstack/master/images/logo.png" width="120">
+
+# Lightstack
+
+Tools for building and processing multi-filter astrophysical datacubes.
+
+## Features
+
+- Cropping FITS images  
+- Building photometric datacubes  
+- PSF matching  
+
+(Currently supports JWST and HST)
+
+## Installation
+
+```bash
+pip install lightstack
+```
+
+## Documentation
+
+https://lightstack.readthedocs.io/en/latest/
+
+See the tutorial notebook!

{lightstack-0.1.0 → lightstack-0.1.2}/pyproject.toml

@@ -1,7 +1,8 @@
 [project]
 name = "lightstack"
-version = "0.1.0"
+version = "0.1.2"
 description = "Tools for building and processing multi-filter astrophysical datacubes"
+readme = "README.md"
 authors = [{name = "Andressa Wille"}, {name="Thallis Pessi"}]
 dependencies = [
     "numpy",

{lightstack-0.1.0 → lightstack-0.1.2}/src/lightstack/datacube.py

@@ -6,7 +6,7 @@ from astropy.wcs import WCS
 
 from reproject import reproject_interp, reproject_exact
 
-from .utils import find_ext, infer_filter
+from .utils import find_ext, infer_filter, get_pixel_scale_wcs
 
 
 def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
@@ -24,7 +24,7 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
     method : str, optional
         Reprojection method:
         - "interp" (default): faster, interpolates values
-        - "exact": slower, conserves flux   --> but reproject_exact has precision issues with resolutions below ~0.05 arcsec, so the results may not be accurate.
+        - "exact": slower   --> but reproject_exact has precision issues with resolutions below ~0.05 arcsec, so the results may not be accurate.
 
     crop : int, optional
         Number of pixels to remove from each border after reprojection
@@ -36,6 +36,7 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
         List in the form [(aligned_fits_path, filter_name), ...].
     """
 
+
     # Choose reprojection method
     if method == "exact":
         reproj_func = reproject_exact
@@ -54,6 +55,9 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
         ref_wcs = WCS(ref_header)
         shape_out = hdul_ref[ext_ref].data.shape
 
+    # Pixel scale of reference
+    scale_out = get_pixel_scale_from_wcs(ref_wcs)
+
     aligned_list = []
 
     # Loop over all FITS
@@ -65,20 +69,30 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
                 print(f"No data extension found in {fpath}. Skipping.")
                 continue
 
-            data = hdul[ext].data
+            data = hdul[ext].data.astype(float)
             header = hdul[ext].header
             wcs_in = WCS(header)
 
             unit = header.get('BUNIT', 'unknown')
             print(f"Filter {filt}: unit = {unit}")
 
+            # Pixel scale of input
+            scale_in = get_pixel_scale_from_wcs(wcs_in)
+
             # Reproject
             data_aligned, footprint = reproj_func(
                 (data, wcs_in),
                 ref_wcs,
                 shape_out=shape_out)
 
-            # Crop border (interp method)
+            # Area correction
+            area_ratio = (scale_out / scale_in)**2
+            data_aligned *= area_ratio
+
+            print(f"Pixel scale in/out: {scale_in:.4f} -> {scale_out:.4f}")
+            print(f"Applied area correction: {area_ratio:.4f}")
+
+            # Crop border
             if crop > 0:
                 data_aligned = data_aligned[crop:-crop, crop:-crop]
 
@@ -90,10 +104,13 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
 
             # Header
             header_aligned = wcs_out.to_header()
+            header_aligned['BUNIT'] = unit
 
             ref_filter = infer_filter(ref_file)
             header_aligned.add_history(
                 f"Reprojected to {ref_filter} using {method} method (reproject package)")
+            header_aligned.add_history(
+                "Flux corrected by pixel area ratio")
 
             if crop > 0:
                 header_aligned.add_history(f"Cropped {crop} pixels from each border")
@@ -105,10 +122,16 @@ def align_reproject_fits(fits_list, ref_file, method="interp", crop=1):
             fits.PrimaryHDU(data_aligned, header=header_aligned).writeto(
                 out_name, overwrite=True)
 
+            # Flux sanity check
+            sum_in = np.nansum(data)
+            sum_out = np.nansum(data_aligned)
+            print(f"Flux ratio (out/in): {sum_out/sum_in:.4f}")
+
             aligned_list.append((out_name, filt))
             print(f"Saved: {out_name}")
 
     return aligned_list
+
     
     
 def build_datacube(aligned_fits_files, reference_file, output_path):

{lightstack-0.1.0 → lightstack-0.1.2}/src/lightstack/psf.py

@@ -125,8 +125,14 @@ def resample_psf(
             raise ValueError("Provide either zoom_factor OR both psf_pixel_scale and target_pixel_scale")
         zoom_factor = psf_pixel_scale / target_pixel_scale
 
+    if zoom_factor > 1:
+    	raise ValueError("Upsampling not recommended.")
+    
     # Resample
     psf_resampled = zoom(psf, zoom_factor, order=order)
+    
+    # Area correction
+    psf_resampled /= zoom_factor**2
 
     # Ensure odd shape
     if make_odd_shape:
@@ -143,6 +149,7 @@ def resample_psf(
         header["HISTORY"] = "PSF resampled using scipy.ndimage.zoom"
         if target_pixel_scale is not None:
             header["CDELT1"] = (target_pixel_scale / 3600, "deg/pix")
+            header["CDELT2"] = (target_pixel_scale / 3600, "deg/pix")
 
     # Save
     fits.PrimaryHDU(psf_resampled, header=header).writeto(

{lightstack-0.1.0 → lightstack-0.1.2}/src/lightstack/utils.py

@@ -210,7 +210,7 @@ def MJy_sr_to_jy(aligned_list):
 
 def get_pixel_scale(fits_path):
     """
-    Compute pixel scale in arcsec/pixel using WCS.
+    Compute pixel scale in arcsec/pixel using WCS. Assumes square pixels and no significant distortion.
     """
     with fits.open(fits_path) as hdul:
         ext = find_ext(hdul)
@@ -221,3 +221,11 @@ def get_pixel_scale(fits_path):
 
     pixscale = proj_plane_pixel_scales(wcs)[0] * 3600.0
     return pixscale
+    
+    
+def get_pixel_scale_from_wcs(wcs):
+    """
+    Compute pixel scale in arcsec/pixel from a WCS object.
+    Assumes square pixels and no significant distortion.
+    """
+    return proj_plane_pixel_scales(wcs)[0] * 3600.0

lightstack-0.1.2/src/lightstack.egg-info/PKG-INFO

@@ -0,0 +1,64 @@
+Metadata-Version: 2.1
+Name: lightstack
+Version: 0.1.2
+Summary: Tools for building and processing multi-filter astrophysical datacubes
+Author: Andressa Wille, Thallis Pessi
+License: MIT License
+        
+        Copyright (c) 2026 Andressa Wille, Thallis Pessi
+        
+        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.
+        
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy
+Requires-Dist: astropy
+Requires-Dist: matplotlib
+Requires-Dist: reproject
+Requires-Dist: regions
+Requires-Dist: scipy
+Requires-Dist: photutils
+
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.19393439.svg)](https://doi.org/10.5281/zenodo.19393439)
+
+<img src="https://raw.githubusercontent.com/AndressaWille/lightstack/master/images/logo.png" width="120">
+
+# Lightstack
+
+Tools for building and processing multi-filter astrophysical datacubes.
+
+## Features
+
+- Cropping FITS images  
+- Building photometric datacubes  
+- PSF matching  
+
+(Currently supports JWST and HST)
+
+## Installation
+
+```bash
+pip install lightstack
+```
+
+## Documentation
+
+https://lightstack.readthedocs.io/en/latest/
+
+See the tutorial notebook!

{lightstack-0.1.0 → lightstack-0.1.2}/src/lightstack.egg-info/SOURCES.txt

@@ -1,3 +1,4 @@
+LICENSE
 README.md
 pyproject.toml
 src/lightstack/__init__.py

lightstack-0.1.0/PKG-INFO

@@ -1,12 +0,0 @@
-Metadata-Version: 2.1
-Name: lightstack
-Version: 0.1.0
-Summary: Tools for building and processing multi-filter astrophysical datacubes
-Author: Andressa Wille, Thallis Pessi
-Requires-Dist: numpy
-Requires-Dist: astropy
-Requires-Dist: matplotlib
-Requires-Dist: reproject
-Requires-Dist: regions
-Requires-Dist: scipy
-Requires-Dist: photutils

lightstack-0.1.0/README.md

@@ -1,18 +0,0 @@
-<img src="images/logo.png" width="120">
-
-# Lightstack
-
-
-This is a code with functions for:
-- cropping FITS images
-- building photometric datacubes
-- PSF matching
-
-(only JWST and HST for now)
-
-pip install git+https://github.com/AndressaWille/lightstack.git
-
-
-Documentation: https://lightstack.readthedocs.io/en/latest/
-
-See the tutorial notebook!

lightstack-0.1.0/src/lightstack.egg-info/PKG-INFO

@@ -1,12 +0,0 @@
-Metadata-Version: 2.1
-Name: lightstack
-Version: 0.1.0
-Summary: Tools for building and processing multi-filter astrophysical datacubes
-Author: Andressa Wille, Thallis Pessi
-Requires-Dist: numpy
-Requires-Dist: astropy
-Requires-Dist: matplotlib
-Requires-Dist: reproject
-Requires-Dist: regions
-Requires-Dist: scipy
-Requires-Dist: photutils