+ Ai?1Rt^RIt^RIHt^RIHtRRltRRltRRlt RRlt Rt Rt R t RR ltR tRR ltRR ltR#)aFunction of unitary fourier transform (uft) and utilities This module implements the unitary fourier transform, also known as the ortho-normal transform. It is especially useful for convolution [1], as it respects the Parseval equality. The value of the null frequency is equal to .. math:: \frac{1}{\sqrt{n}} \sum_i x_i so the Fourier transform has the same energy as the original image (see ``image_quad_norm`` function). The transform is applied from the last axis for performance (assuming a C-order array input). References ---------- .. [1] B. R. Hunt "A matrix theory proof of the discrete convolution theorem", IEEE Trans. on Audio and Electroacoustics, vol. au-19, no. 4, pp. 285-288, dec. 1971 N)_supported_float_typecnVf VPp\P!V\V)^4RR7pV#)aN-dimensional unitary Fourier transform. Parameters ---------- inarray : ndarray The array to transform. dim : int, optional The last axis along which to compute the transform. All axes by default. Returns ------- outarray : ndarray (same shape than inarray) The unitary N-D Fourier transform of ``inarray``. Examples -------- >>> input = np.ones((3, 3, 3)) >>> output = ufftn(input) >>> np.allclose(np.sum(input) / np.sqrt(input.size), output[0, 0, 0]) True >>> output.shape (3, 3, 3) orthoaxesnorm)ndimfftfftnrangeinarraydimoutarrays&& U/var/www/html/photoedit/myenv/lib/python3.14/site-packages/skimage/restoration/uft.pyufftnrs12 {llxxeSD!n7CH OcnVf VPp\P!V\V)^4RR7pV#)aN-dimensional unitary inverse Fourier transform. Parameters ---------- inarray : ndarray The array to transform. dim : int, optional The last axis along which to compute the transform. All axes by default. Returns ------- outarray : ndarray The unitary inverse nD Fourier transform of ``inarray``. Has the same shape as ``inarray``. Examples -------- >>> input = np.ones((3, 3, 3)) >>> output = uifftn(input) >>> np.allclose(np.sum(input) / np.sqrt(input.size), output[0, 0, 0]) True >>> output.shape (3, 3, 3) rr)rr ifftnr r s&& ruifftnr;s14 {llyyucT1~GDH OrcnVf VPp\P!V\V)^4RR7pV#)a$N-dimensional real unitary Fourier transform. This transform considers the Hermitian property of the transform on real-valued input. Parameters ---------- inarray : ndarray, shape (M[, ...], P) The array to transform. dim : int, optional The last axis along which to compute the transform. All axes by default. Returns ------- outarray : ndarray, shape (M[, ...], P / 2 + 1) The unitary N-D real Fourier transform of ``inarray``. Notes ----- The ``urfft`` functions assume an input array of real values. Consequently, the output has a Hermitian property and redundant values are not computed or returned. Examples -------- >>> input = np.ones((5, 5, 5)) >>> output = urfftn(input) >>> np.allclose(np.sum(input) / np.sqrt(input.size), output[0, 0, 0]) True >>> output.shape (5, 5, 3) rr)rr rfftnr r s&& rurfftnr[s2D {llyyucT1~GDH OrcnVf VPp\P!W\V)^4RR7pV#)aN-dimensional inverse real unitary Fourier transform. This transform considers the Hermitian property of the transform from complex to real input. Parameters ---------- inarray : ndarray The array to transform. dim : int, optional The last axis along which to compute the transform. All axes by default. shape : tuple of int, optional The shape of the output. The shape of ``rfft`` is ambiguous in case of odd-valued input shape. In this case, this parameter should be provided. See ``np.fft.irfftn``. Returns ------- outarray : ndarray The unitary N-D inverse real Fourier transform of ``inarray``. Notes ----- The ``uirfft`` function assumes that the output array is real-valued. Consequently, the input is assumed to have a Hermitian property and redundant values are implicit. Examples -------- >>> input = np.ones((5, 5, 5)) >>> output = uirfftn(urfftn(input), shape=input.shape) >>> np.allclose(input, output) True >>> output.shape (5, 5, 5) rr)rr irfftnr )r rshapers&&& ruirfftnrs2L {llzz'ucT1~GLH Orc\V^4#)a2-dimensional unitary Fourier transform. Compute the Fourier transform on the last 2 axes. Parameters ---------- inarray : ndarray The array to transform. Returns ------- outarray : ndarray (same shape as inarray) The unitary 2-D Fourier transform of ``inarray``. See Also -------- uifft2, ufftn, urfftn Examples -------- >>> input = np.ones((10, 128, 128)) >>> output = ufft2(input) >>> np.allclose(np.sum(input[1, ...]) / np.sqrt(input[1, ...].size), ... output[1, 0, 0]) True >>> output.shape (10, 128, 128) )rr s&rufft2rs: ! rc\V^4#)a32-dimensional inverse unitary Fourier transform. Compute the inverse Fourier transform on the last 2 axes. Parameters ---------- inarray : ndarray The array to transform. Returns ------- outarray : ndarray (same shape as inarray) The unitary 2-D inverse Fourier transform of ``inarray``. See Also -------- uifft2, uifftn, uirfftn Examples -------- >>> input = np.ones((10, 128, 128)) >>> output = uifft2(input) >>> np.allclose(np.sum(input[1, ...]) / np.sqrt(input[1, ...].size), ... output[0, 0, 0]) True >>> output.shape (10, 128, 128) )rrs&ruifft2r!s: '1 rc\V^4#)a2-dimensional real unitary Fourier transform Compute the real Fourier transform on the last 2 axes. This transform considers the Hermitian property of the transform from complex to real-valued input. Parameters ---------- inarray : ndarray, shape (M[, ...], P) The array to transform. Returns ------- outarray : ndarray, shape (M[, ...], 2 * (P - 1)) The unitary 2-D real Fourier transform of ``inarray``. See Also -------- ufft2, ufftn, urfftn Examples -------- >>> input = np.ones((10, 128, 128)) >>> output = urfft2(input) >>> np.allclose(np.sum(input[1,...]) / np.sqrt(input[1,...].size), ... output[1, 0, 0]) True >>> output.shape (10, 128, 65) )rrs&rurfft2r#s> '1 rc\V^VR7#)ai2-dimensional inverse real unitary Fourier transform. Compute the real inverse Fourier transform on the last 2 axes. This transform considers the Hermitian property of the transform from complex to real-valued input. Parameters ---------- inarray : ndarray, shape (M[, ...], P) The array to transform. shape : tuple of int, optional The shape of the output. The shape of ``rfft`` is ambiguous in case of odd-valued input shape. In this case, this parameter should be provided. See ``np.fft.irfftn``. Returns ------- outarray : ndarray, shape (M[, ...], 2 * (P - 1)) The unitary 2-D inverse real Fourier transform of ``inarray``. See Also -------- urfft2, uifftn, uirfftn Examples -------- >>> input = np.ones((10, 128, 128)) >>> output = uirfftn(urfftn(input), shape=input.shape) >>> np.allclose(input, output) True >>> output.shape (10, 128, 128) )r)r)r rs&&ruirfft2r%sD 7AU ++rcVPR,VPR,8wd^\P!\P!\P!V4^,RR7RR7,\P!\P!VR,4^,RR7, #\P!\P!\P!V4^,RR7RR7#)aReturn the quadratic norm of images in Fourier space. This function detects whether the input image satisfies the Hermitian property. Parameters ---------- inarray : ndarray Input image. The image data should reside in the final two axes. Returns ------- norm : float The quadratic norm of ``inarray``. Examples -------- >>> input = np.ones((5, 5)) >>> image_quad_norm(ufft2(input)) == np.sum(np.abs(input)**2) True >>> image_quad_norm(ufft2(input)) == image_quad_norm(urfft2(input)) True )axis).)rnpsumabsrs&rimage_quad_normr.6s4}}RGMM"--266"&&A!5B?bIIBFF FF76? #q (rM   vvbffRVVG_1;"EErc V'g VPp\VP4p\P!WR7pY\ VP Uu.uFp\^V4NK up4&\VP 4F\wrxWpPV, 8gK\P!V\\P!V^, 44)VR7pK^ V'd\PM\Pp V !V\V)^4R7p \P !V\P"4p V P%V RR7#uupi)aCompute the transfer function of an impulse response (IR). This function makes the necessary correct zero-padding, zero convention, correct fft2, etc... to compute the transfer function of IR. To use with unitary Fourier transform for the signal (ufftn or equivalent). Parameters ---------- imp_resp : ndarray The impulse responses. shape : tuple of int A tuple of integer corresponding to the target shape of the transfer function. dim : int, optional The last axis along which to compute the transform. All axes by default. is_real : bool, optional If True (default), imp_resp is supposed real and the Hermitian property is used with rfftn Fourier transform. Returns ------- y : complex ndarray The transfer function of shape ``shape``. See Also -------- ufftn, uifftn, urfftn, uirfftn Examples -------- >>> np.all(np.array([[4, 0], [0, 0]]) == ir2tf(np.ones((2, 2)), (2, 2))) True >>> ir2tf(np.ones((2, 2)), (512, 512)).shape == (512, 257) True >>> ir2tf(np.ones((2, 2)), (512, 512), is_real=False).shape == (512, 512) True Notes ----- The input array can be composed of multiple-dimensional IR with an arbitrary number of IR. The individual IR must be accessed through the first axes. The last ``dim`` axes contain the space definition. )dtype)shiftr')rF)copy)rrr0r+zerostuplerslice enumeraterollintfloorr rr r promote_types complex64astype) imp_resprris_realirpadded_dtypeirpaddedsr' axis_sizefuncout cplx_dtypes &&&& rir2tfrFXs^ mm*8>>:Nxx4H=EU8AE!QK8 9:%X^^4 ==3& &wwxBHHY]4K0L/LSWXH5 399SXXD xucT1~ /C!!.",,?J ::ju: --9sE c\P!^.V,4p\V4Fp\\ ^^4.V,\ R4.,\ ^^4.W, ^, ,,4p\P !.RO4P \V4Uu.uF qfV8XdRM^NK up4W5&K RV,V\ ^^43V,&\W1VR7V3#uupi)aReturn the transfer function of the Laplacian. Laplacian is the second order difference, on row and column. Parameters ---------- ndim : int The dimension of the Laplacian. shape : tuple The support on which to compute the transfer function. is_real : bool, optional If True (default), imp_resp is assumed to be real-valued and the Hermitian property is used with rfftn Fourier transform to return the transfer function. Returns ------- tf : array_like, complex The transfer function. impr : array_like, real The Laplacian. Examples -------- >>> tf, ir = laplacian(2, (32, 32)) >>> np.all(ir == np.array([[0, -1, 0], [-1, 4, -1], [0, -1, 0]])) True >>> np.all(tf == ir2tf(ir, (32, 32))) True Ng@)r>)grHr()r+r3r r4r5arrayreshaperF)rrr>imprridxis&&& r laplacianrNs> 88QC$J DT{ 1a[MC 5;- /5A;-4:PQ>2R R HH./77,1$K 8Kq8R "K 8  #&*D%1+$  g . 44 9s/C< )N)NN)NT)T)__doc__numpyr+ scipy.fftr _shared.utilsrrrrrrr!r#r%r.rFrNrrrTsW*1>@%P)X@@D",JFD@.F(5r