foscat 3.6.1__py3-none-any.whl → 3.7.1__py3-none-any.whl

foscat/scat_cov2D.py

@@ -1,5 +1,5 @@
 import foscat.scat_cov as scat
-
+import numpy as np
 
 class scat_cov2D:
     def __init__(self, s0, s2, s3, s4, s1=None, s3p=None, backend=None):
@@ -11,8 +11,108 @@ class scat_cov2D:
     def fill(self, im, nullval=0):
         return self.fill_2d(im, nullval=nullval)
 
-
 class funct(scat.funct):
     def __init__(self, *args, **kwargs):
         # Impose que use_2D=True pour la classe scat
         super().__init__(use_2D=True, *args, **kwargs)
+     
+    def spectrum(self,image):
+        """
+        Compute the 1D power spectrum of a 2D image by averaging the 2D power spectrum
+        over concentric frequency rings (radial averaging), using np.bincount for efficiency.
+
+        Parameters:
+        - image : ndarray (2D), input image
+
+        Returns:
+        - freq : radial frequencies
+        - spectrum_1d : corresponding 1D power spectrum
+        """
+        import numpy as np
+        
+        # Compute the 2D Fourier Transform and shift the zero frequency to the center
+        F = np.fft.fftshift(np.fft.fft2(image))
+        power_spectrum = np.abs(F) ** 2
+
+        # Create coordinate grids and compute the radial distance from the center
+        y, x = np.indices(power_spectrum.shape)
+        center = np.array(power_spectrum.shape) // 2
+        r = np.sqrt((x - center[1])**2 + (y - center[0])**2).astype(int)
+
+        # Use np.bincount for fast summation and counting
+        sum_power = np.bincount(r.ravel(), weights=power_spectrum.ravel())
+        counts = np.bincount(r.ravel())
+
+        # Compute the mean power for each radial bin
+        spectrum_1d = sum_power / counts
+
+        return spectrum_1d
+
+    def plot_results(self,in_image,out_image,vmin=None,vmax=None,cmap='coolwarm',spec_range=None):
+        import matplotlib.pyplot as plt
+        
+        if len(out_image.shape)>2:
+            nimage=out_image.shape[0]
+            ndraw=np.min([3,nimage])
+            plt.figure(figsize=(16,12))
+            plt.subplot(2,ndraw+1,1)
+            plt.title('Original field') 
+            plt.imshow(in_image,cmap=cmap,vmin=vmin,vmax=vmax,origin='lower')
+            plt.xticks([])
+            plt.yticks([])
+            for k in range(ndraw):
+                plt.subplot(2,ndraw+1,2+k)
+                plt.title('Modeled field #%d'%(k)) 
+                plt.imshow(out_image[k],cmap=cmap,vmin=vmin,vmax=vmax,origin='lower')
+                plt.xticks([])
+                plt.yticks([])
+            plt.subplot(2,2,3)
+            plt.title('Histogram') 
+            for k in range(nimage):
+                if k==0:
+                    plt.hist(out_image[k].flatten(),bins=100,label='modeled',color='b',histtype='step',log=True,alpha=0.5)
+                else:
+                    plt.hist(out_image[k].flatten(),bins=100,color='b',histtype='step',log=True,alpha=0.5)
+            plt.hist(in_image.flatten(),bins=100,label='original',color='r',histtype='step',log=True)
+            plt.legend(frameon=0)
+            plt.subplot(2,2,4)
+            plt.title('Powerspectra') 
+            for k in range(nimage):
+                if k==0:
+                    plt.plot(self.spectrum(out_image[k]),color='b',label='modeled',alpha=0.5)
+                else:
+                    plt.plot(self.spectrum(out_image[k]),color='b',alpha=0.5)
+            plt.plot(self.spectrum(in_image),color='r',label='original')
+            plt.xscale('log')
+            plt.yscale('log')
+            plt.legend(frameon=0)
+            if spec_range is not None:
+                plt.ylim(spec_range[0],spec_range[1])
+        else:
+            plt.figure(figsize=(16,3))
+            plt.subplot(1,4,1)
+            plt.title('Original field') 
+            plt.imshow(in_image,cmap=cmap,vmin=vmin,vmax=vmax,origin='lower')
+            plt.xticks([])
+            plt.yticks([])
+            plt.subplot(1,4,2)
+            plt.title('Modeled field') 
+            plt.imshow(out_image,cmap=cmap,vmin=vmin,vmax=vmax,origin='lower')
+            plt.xticks([])
+            plt.yticks([])
+            plt.subplot(1,4,3)
+            plt.title('Histogram') 
+            plt.hist(in_image.flatten(),bins=100,label='original',color='r',histtype='step',log=True)
+            plt.hist(out_image.flatten(),bins=100,label='modeled',color='b',histtype='step',log=True)
+            plt.legend(frameon=0)
+            plt.subplot(1,4,4)
+            plt.title('Powerspectra') 
+            plt.plot(self.spectrum(in_image),color='b',label='original')
+            plt.plot(self.spectrum(out_image),color='r',label='modeled')
+            plt.xscale('log')
+            plt.yscale('log')
+            if spec_range is not None:
+                plt.ylim(spec_range[0],spec_range[1])
+            plt.legend(frameon=0)
+        
+        

foscat/scat_cov_map.py

@@ -24,10 +24,23 @@ class funct(scat.funct):
         super().__init__(return_data=True, *args, **kwargs)
 
     def eval(
-            self, image1, image2=None, mask=None, norm=None, Auto=True, calc_var=False,out_nside=None
+        self,
+        image1,
+        image2=None,
+        mask=None,
+        norm=None,
+        Auto=True,
+        calc_var=False,
+        out_nside=None,
     ):
         r = super().eval(
-            image1, image2=image2, mask=mask, norm=norm, Auto=Auto, calc_var=calc_var,out_nside=out_nside
+            image1,
+            image2=image2,
+            mask=mask,
+            norm=norm,
+            Auto=Auto,
+            calc_var=calc_var,
+            out_nside=out_nside,
         )
         return scat_cov_map(
             r.S2, r.S0, r.S3, r.S4, S1=r.S1, S3P=r.S3P, backend=r.backend

{foscat-3.6.1.dist-info → foscat-3.7.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: foscat
-Version: 3.6.1
+Version: 3.7.1
 Summary: Generate synthetic Healpix or 2D data using Cross Scattering Transform
 Author-email: Jean-Marc DELOUIS <jean.marc.delouis@ifremer.fr>
 Maintainer-email: Theo Foulquier <theo.foulquier@ifremer.fr>
@@ -18,7 +18,7 @@ Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Requires-Python: >=3.9
 Description-Content-Type: text/markdown
-License-File: LICENCE
+License-File: LICENSE
 Requires-Dist: imageio
 Requires-Dist: imagecodecs
 Requires-Dist: matplotlib

foscat-3.7.1.dist-info/RECORD

@@ -0,0 +1,26 @@
+foscat/CNN.py,sha256=j0F2a4Xf3LijhyD_WVZ6Eg_IjGuXw3ddH6Iudj1xVaw,4874
+foscat/CircSpline.py,sha256=CXi49FxF8ZoeZ17Ua8c1AZXe2B5ICEC9aCXb97atB3s,4028
+foscat/FoCUS.py,sha256=bYQctBoTiKyJd9xtIuPgcroek4XxIEnUmdirkSxdcmQ,102659
+foscat/GCNN.py,sha256=5RV-FKuvqbD-k99TwiM4CttM2LMZE21WD0IK0j5Mkko,7599
+foscat/Softmax.py,sha256=aBLQauoG0q2SJYPotV6U-cxAhsJcspWHNRWdnA_nAiQ,2854
+foscat/Spline1D.py,sha256=rKzzenduaZZ-yBDJd35it6Gyrj1spqb7hoIaUgISPzY,2983
+foscat/Synthesis.py,sha256=3oL-WIwOzjO6jYn-7J3TfAQSiPaQTqGtrmxcGBh1Gvs,13787
+foscat/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+foscat/alm.py,sha256=PaAaq82CH3JzH520eQpqTCivxE4VeJmmc6I_9HXIhcA,32480
+foscat/backend.py,sha256=b7c6nL4vaolypzcmhzsZm7vzsJhjrkd--l5DjtESXEA,44851
+foscat/backend_tens.py,sha256=9Dp136m9frkclkwifJQLLbIpl3ETI3_txdPUZcKfuMw,1618
+foscat/loss_backend_tens.py,sha256=dCOVN6faDtIpN3VO78HTmYP2i5fnFAf-Ddy5qVBlGrM,1783
+foscat/loss_backend_torch.py,sha256=k3z18Dj3SaLKK6ZIKcm7GO4U_YKYVP6LtHG1aIbxkYk,1627
+foscat/scat.py,sha256=qGYiBIysPt65MdmF07WWA4piVlTfA9-lFDTaicnqC2w,72822
+foscat/scat1D.py,sha256=W5Uu6wdQ4ZsFKXpof0f1OBl-1wjJmW7ruvddRWxe7uM,53726
+foscat/scat2D.py,sha256=boKj0ASqMMSy7uQLK6hPniG87m3hZGJBYBiq5v8F9IQ,532
+foscat/scat_cov.py,sha256=mHGOFa6ZnmGksOSa_zQQ-92FlRs2D1Z1zmkA-22X0Tg,232476
+foscat/scat_cov1D.py,sha256=XOxsZZ5TYq8f34i2tUgIfzyaqaTDlICB3HzD2l_puro,531
+foscat/scat_cov2D.py,sha256=HVYBNtpZgRlz6OOm0LwG_w8Q4Kld4Qc1n6XYeGoxIug,4802
+foscat/scat_cov_map.py,sha256=Swt39-nYEaQkBzyX4EOAQBvUuYQpERzJ-uVxSWS2b-Y,2911
+foscat/scat_cov_map2D.py,sha256=FqF45FBcoiQbvuVsrLWUIPRUc95GsKsrnH6fKzB3GlE,2841
+foscat-3.7.1.dist-info/LICENSE,sha256=i0ukIr8ZUpkSY2sZaE9XZK-6vuSU5iG6IgX_3pjatP8,1505
+foscat-3.7.1.dist-info/METADATA,sha256=Jgds24vLXVaJ2W3FaBHysyu19qY2hrywOpn4ggy46rA,7216
+foscat-3.7.1.dist-info/WHEEL,sha256=EaM1zKIUYa7rQnxGiOCGhzJABRwy4WO57rWMR3_tj4I,91
+foscat-3.7.1.dist-info/top_level.txt,sha256=AGySXBBAlJgb8Tj8af6m_F-aiNg2zNTcybCUPVOKjAg,7
+foscat-3.7.1.dist-info/RECORD,,

{foscat-3.6.1.dist-info → foscat-3.7.1.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.42.0)
+Generator: setuptools (75.9.1)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

foscat/alm_tools.py

@@ -1,11 +0,0 @@
-import numpy as np
-
-#====================================================================================================================
-# This class is an automatic traduction of the fortran healpix software
-#====================================================================================================================
-
-
-class alm_tools():
-    def __init__(self):
-        pass
-

foscat-3.6.1.dist-info/RECORD

@@ -1,27 +0,0 @@
-foscat/CNN.py,sha256=j0F2a4Xf3LijhyD_WVZ6Eg_IjGuXw3ddH6Iudj1xVaw,4874
-foscat/CircSpline.py,sha256=DjP1gy88cnXu2O21ww_lNnsHAHXc3OAWk_8ey84yicg,4053
-foscat/FoCUS.py,sha256=YoZqPiDLQanUXGgN4RyN6S4_F2rHn273P6QQjq0hlBU,101774
-foscat/GCNN.py,sha256=5RV-FKuvqbD-k99TwiM4CttM2LMZE21WD0IK0j5Mkko,7599
-foscat/Softmax.py,sha256=aBLQauoG0q2SJYPotV6U-cxAhsJcspWHNRWdnA_nAiQ,2854
-foscat/Spline1D.py,sha256=a5Jb8I9tb8y20iM8W-z6iZsIqDFByRp6eZdChpmuI5k,3010
-foscat/Synthesis.py,sha256=3_Lq5-gUM-WmO2h15kajMES8XjRo2BGseoxvTLW_xEc,13626
-foscat/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-foscat/alm.py,sha256=IIvYOz0kAtufhlBzS_QOx891-Y30ARmofyqomOFmJDY,29335
-foscat/alm_tools.py,sha256=BMObyIMF3_v04JdvMWBMMaF2F8s2emvDKmirKDnHWDA,387
-foscat/backend.py,sha256=bN-b0CWcJXVsLCyqMhJACom0JlSStgrMWyD0uB6HqoU,40218
-foscat/backend_tens.py,sha256=9Dp136m9frkclkwifJQLLbIpl3ETI3_txdPUZcKfuMw,1618
-foscat/loss_backend_tens.py,sha256=dCOVN6faDtIpN3VO78HTmYP2i5fnFAf-Ddy5qVBlGrM,1783
-foscat/loss_backend_torch.py,sha256=k3z18Dj3SaLKK6ZIKcm7GO4U_YKYVP6LtHG1aIbxkYk,1627
-foscat/scat.py,sha256=qGYiBIysPt65MdmF07WWA4piVlTfA9-lFDTaicnqC2w,72822
-foscat/scat1D.py,sha256=W5Uu6wdQ4ZsFKXpof0f1OBl-1wjJmW7ruvddRWxe7uM,53726
-foscat/scat2D.py,sha256=boKj0ASqMMSy7uQLK6hPniG87m3hZGJBYBiq5v8F9IQ,532
-foscat/scat_cov.py,sha256=gQvGC2s93oDurBDW0A6mGyzzMETjYAuBsBsDOI-ZqUY,150642
-foscat/scat_cov1D.py,sha256=XOxsZZ5TYq8f34i2tUgIfzyaqaTDlICB3HzD2l_puro,531
-foscat/scat_cov2D.py,sha256=3gn6xjKvfKsyHJoPfYIu8q9LLVAbU3tsiS2l1LAJ0XM,531
-foscat/scat_cov_map.py,sha256=9Yymbr1MxUNY5nJd9kIEEHt1x2IoOjc0EW4kkJVtmQ4,2783
-foscat/scat_cov_map2D.py,sha256=FqF45FBcoiQbvuVsrLWUIPRUc95GsKsrnH6fKzB3GlE,2841
-foscat-3.6.1.dist-info/LICENCE,sha256=i0ukIr8ZUpkSY2sZaE9XZK-6vuSU5iG6IgX_3pjatP8,1505
-foscat-3.6.1.dist-info/METADATA,sha256=6t347EDYxBbM5qhlC5T6pwbJDBwvAnJ5VJnNppblC6A,7216
-foscat-3.6.1.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-foscat-3.6.1.dist-info/top_level.txt,sha256=AGySXBBAlJgb8Tj8af6m_F-aiNg2zNTcybCUPVOKjAg,7
-foscat-3.6.1.dist-info/RECORD,,