imicpe 1.0.0__tar.gz → 1.0.2__tar.gz

{imicpe-1.0.0 → imicpe-1.0.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: imicpe
-Version: 1.0.0
+Version: 1.0.2
 Summary: Toolbox for Maths,Signal,Image Teaching @ CPE
 Author-email: Marion Foare <marion.foare@cpe.fr>, Eric Van Reeth <eric.vanreeth@cpe.fr>, Arthur Gautheron <arthur.gautheron@cpe.fr>
 License: MIT License

{imicpe-1.0.0 → imicpe-1.0.2}/pyproject.toml

@@ -14,7 +14,7 @@ mypkg = ["*.txt", "*.mat", "*.npy"]
 #[tool.setuptools_scm]
 
 [project]
-version="1.0.0"
+version="1.0.2"
 name = "imicpe"
 authors = [
   { name="Marion Foare", email="marion.foare@cpe.fr" },

imicpe-1.0.2/src/imicpe/_version.py

@@ -0,0 +1 @@
+__version__="1.0.2"

{imicpe-1.0.0 → imicpe-1.0.2}/src/imicpe/cs/l1.py

@@ -32,16 +32,15 @@ def l1(opreg,A,At,z,x0,lam):
     
     ### init ###
     dim = x0.ndim
-    match opreg:
-        case 'id':
-            G = Id
-            Gt = Id
-        case 'gradient':
-            G = D
-            Gt = Dt
-        case 'laplacien':
-            G = L
-            Gt = Lt
+    if opreg == 'id':
+        G = Id
+        Gt = Id
+    elif opreg == 'gradient':
+        G = D
+        Gt = Dt
+    if opreg == 'laplacien':
+        G = L
+        Gt = Lt
     
     # operator norms
     lipA = opNorm(A,At,dim,x0)

{imicpe-1.0.0 → imicpe-1.0.2}/src/imicpe/cs/operators.py

@@ -9,7 +9,7 @@ import pywt
 ############################################################
 def Id(x):
     """
-    Id       Opérateur identité
+    Opérateur identité
    
     Args:
         X (numpy.ndarray)               signal 1D
@@ -28,7 +28,7 @@ def Id(x):
 # gradient
 def D(x):
     """
-    D       Calcule le gradient par différences finies à droite.
+    Calcule le gradient par différences finies à droite.
     Autrement dit, D(x) calcule le produit matriciel Dx.
    
     Args:
@@ -52,7 +52,7 @@ def D(x):
 
 def Dt(x):
     """
-    Dt       Calcule l’adjoint gradient par différences finies à droite.
+    Calcule l’adjoint gradient par différences finies à droite.
     Autrement dit, Dt(x) calcule le produit matriciel D'x.
    
     Args:
@@ -76,7 +76,7 @@ def Dt(x):
 # laplacian
 def L(x):
     """
-    L        Calcule la dérivée seconde d’un signal, ou le laplacien dans le cas d’une image.
+    Calcule la dérivée seconde d’un signal, ou le laplacien dans le cas d’une image.
     Autrement dit, L(x) calcule le produit matriciel Lx.
    
     Args:
@@ -99,7 +99,7 @@ def L(x):
 
 def Lt(x):
     """
-    Lt       Calcule l’adjoint du laplacien.
+    Calcule l’adjoint du laplacien.
     Autrement dit, Lt(x) calcule le produit matriciel L'x.
    
     Args:
@@ -126,7 +126,7 @@ def Lt(x):
 ############################################################
 def fwt(x,wavelet,level):
     """
-    P       Calcule la transformée en ondelettes directe 1D.
+    Calcule la transformée en ondelettes directe 1D.
 
     Args:
     x (numpy.ndarray)               signal 1D
@@ -143,12 +143,14 @@ def fwt(x,wavelet,level):
     return coeff_arr
 
 
-def iwt(x,wavelet,level,coeffs_slices=None,coeffs_shapes=None,):
+def iwt(x,wavelet,level):
     """
-    P       Calcule la transformée en ondelettes inverse 1D.
+    Calcule la transformée en ondelettes inverse 1D.
 
     Args:
-    x (numpy.ndarray)               Vecteur des coefficients d'ondelettes                                          
+    x (numpy.ndarray)               Vecteur des coefficients d'ondelettes                
+    wavelet (string)                Nom de l'ondelette mère (voir la librairie pywt)     
+    level (int)                     Niveau de décompisition                                     
 
     Returns: 
     (numpy.ndarray)                 Signal correspondant aux coefficients d'ondelettes donnés par x
@@ -156,6 +158,7 @@ def iwt(x,wavelet,level,coeffs_slices=None,coeffs_shapes=None,):
 
     J = level
     N = len(x)
+    _, coeffs_slices, coeffs_shapes = pywt.ravel_coeffs(pywt.wavedec(np.ones(N), wavelet=wavelet, level=J, mode='periodization'))
     
     if coeffs_shapes is None:
         # compute coeffs size at each level
@@ -178,10 +181,12 @@ def iwt(x,wavelet,level,coeffs_slices=None,coeffs_shapes=None,):
 
 def fwt2(x,wavelet,level):
     """
-    P       Calcule la transformée en ondelettes directe 2D.
+    Calcule la transformée en ondelettes directe 2D.
 
     Args:
-    x (numpy.ndarray)               image NON vectorisée 2D
+    x (numpy.ndarray)               image NON vectorisée 2D        
+    wavelet (string)                Nom de l'ondelette mère (voir la librairie pywt)     
+    level (int)                     Niveau de décompisition                    
 
     Returns: 
     (numpy.ndarray)                 Vecteur des coefficients de la décomposition en ondelettes de x
@@ -192,12 +197,14 @@ def fwt2(x,wavelet,level):
 
     return coeff_arr
 
-def iwt2(x,wavelet,level,coeffs_slices,coeffs_shapes):
+def iwt2(x,wavelet,level):
     """
-    P       Calcule la transformée en ondelettes inverse 2D.
+    Calcule la transformée en ondelettes inverse 2D.
 
     Args:
-    x (numpy.ndarray)               Vecteur des coefficients d'ondelettes                                          
+    x (numpy.ndarray)               Vecteur des coefficients d'ondelettes            
+    wavelet (string)                Nom de l'ondelette mère (voir la librairie pywt)     
+    level (int)                     Niveau de décompisition                                                          
 
     Returns: 
     (numpy.ndarray)                 Image correspondante aux coefficients d'ondelettes donnés par x
@@ -205,6 +212,8 @@ def iwt2(x,wavelet,level,coeffs_slices,coeffs_shapes):
 
     J = level
     N = len(x)
+
+    _, coeffs_slices, coeffs_shapes = pywt.ravel_coeffs(pywt.wavedec2(np.ones((int(np.sqrt(N)),int(np.sqrt(N)))), wavelet=wavelet, level=J, mode='periodization'))
     
     coeffs = pywt.unravel_coeffs(x, coeffs_slices, coeffs_shapes,
                                       output_format='wavedec2')
@@ -220,7 +229,7 @@ def iwt2(x,wavelet,level,coeffs_slices,coeffs_shapes):
 ############################################################
 def opNorm(op,opt,dim,xn):
     """
-    opNorm      Calcule la norme de l'opérateur OP, dont 
+    Calcule la norme de l'opérateur OP, dont 
     l'opérateur transposé est OPT, en dimension DIM
 
     Args:
@@ -264,7 +273,7 @@ def opNorm(op,opt,dim,xn):
 
 def matNorm(M):
     """
-    matNorm      Calcule la norme de la matrice M
+    Calcule la norme de la matrice M
 
     Args:
         M (numpy.ndarray)               matrice dont on souhaite calculer la norme

{imicpe-1.0.0 → imicpe-1.0.2}/src/imicpe/cs/tikhonov.py

@@ -29,16 +29,15 @@ def tikhonov(opreg,A,At,z,x0,lam):
     
     ### init ###
     dim = x0.ndim
-    match opreg:
-        case 'id':
-            G = Id
-            Gt = Id
-        case 'gradient':
-            G = D
-            Gt = Dt
-        case 'laplacien':
-            G = L
-            Gt = Lt
+    if opreg == 'id':
+        G = Id
+        Gt = Id
+    elif opreg == 'gradient':
+        G = D
+        Gt = Dt
+    if opreg == 'laplacien':
+        G = L
+        Gt = Lt
     
     # operator norms
     lipA = opNorm(A,At,dim,x0)

{imicpe-1.0.0 → imicpe-1.0.2}/src/imicpe.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: imicpe
-Version: 1.0.0
+Version: 1.0.2
 Summary: Toolbox for Maths,Signal,Image Teaching @ CPE
 Author-email: Marion Foare <marion.foare@cpe.fr>, Eric Van Reeth <eric.vanreeth@cpe.fr>, Arthur Gautheron <arthur.gautheron@cpe.fr>
 License: MIT License

imicpe-1.0.0/src/imicpe/_version.py

@@ -1 +0,0 @@
-__version__="1.0.0"