a d@sdZddlmZmZmZgdZddlmZddlm Z ddl Z ddl m Z mZdd lmZmZmZdd lmZmZmZmZdd dZGdddeZGdddeZGdddeZGdddeZGdddeZGdddeZGdddeZ dS)z*1D and 2D Wavelet packet transform module.)divisionprint_functionabsolute_import)BaseNodeNode WaveletPacketNode2DWaveletPacket2DNodeNDWaveletPacketND)product) OrderedDictN)Wavelet _check_dtype)dwtidwt dwt_max_level)dwt2idwt2dwtnidwtnadcsLg}t|dD]2}fdd|Dfdd|dddD}q|S)Nrcsg|] }|qSr.0path)xr^/var/www/html/stable-diffusion-webui/venv/lib/python3.9/site-packages/pywt/_wavelet_packets.py z&get_graycode_order..csg|] }|qSrrr)yrrr r!)range)levelrr"graycode_orderirrr"rget_graycode_orders r)c@s eZdZdZdZdZddZddZd;dd Zdd"d#Zd$d%Zd?d&d'Zd(d)Zd*d+Zd,d-Zed.d/Zed0d1Zd@d2d3ZdAd5d6ZdBd7d8Zd9d:Z dS)Cra BaseNode for wavelet packet 1D and 2D tree nodes. The BaseNode is a base class for `Node` and `Node2D`. It should not be used directly unless creating a new transformation type. It is included here to document the common interface of 1D and 2D node and wavelet packet transform classes. Parameters ---------- parent : Parent node. If parent is None then the node is considered detached (ie root). data : 1D or 2D array Data associated with the node. 1D or 2D numeric array, depending on the transform type. node_name : A name identifying the coefficients type. See `Node.node_name` and `Node2D.node_name` for information on the accepted subnodes names. NcCs||_|durH|j|_|j|_|jd|_|j|_|j||_|j|_nd|_d|_d|_d|_d|_||_|jdur~d|_ nt |j |_ | dS)Nrr)parentwaveletmoder%maxlevel _maxlevelraxesdata _data_shapenpasarrayshape_init_subnodes)selfr+r1 node_namerrr__init__:s$    zBaseNode.__init__cCs|jD]}||dqdSN)PARTS _set_noder7partrrrr6Us zBaseNode._init_subnodesTcCs tdSr:NotImplementedErrorr7r>r1 overwriterrr_create_subnodeYszBaseNode._create_subnodecKsJ|||s&||dur&||S||||fi|}||||Sr:)_validate_node_name _get_noder<)r7node_clsr>r1rBkwargsnoderrr_create_subnode_base\s    zBaseNode._create_subnode_basecCs t||Sr:)getattrr=rrrrEeszBaseNode._get_nodecCst|||dSr:)setattrr7r>rHrrrr<hszBaseNode._set_nodecCs||ddSr:r<r=rrr _delete_nodekszBaseNode._delete_nodecCs0||jvr,tdddd|jD|fdS)N'Subnode name must be in [%s], not '%s'., css|]}d|VqdSz'%s'Nrrprrr qr!z/BaseNode._validate_node_name..)r; ValueErrorjoinr=rrrrDns zBaseNode._validate_node_namecs6jtj}tfddtd|dDS)zThe path to the current node in tuple form. The length of the tuple is equal to the number of decomposition levels. cs(g|] }|dj|jqSr)PART_LEN)rnrr7rrr {sz'BaseNode.path_tuple..r)rlenrXtupler$)r7ZnlevrrZr path_tupless  zBaseNode.path_tupler+cCs|dvs J|jdur|jS|jdur@|jtt|jj|jS|dkr`|jdur|j|SnB|dkr|j D]2}t ||d}|durn||}|durn|SqndSz Try to find the value of maximum decomposition level if it is not specified explicitly. Parameters ---------- evaluate_from : {'parent', 'subnodes'} )r+subnodesNr+r_) r/r1r%rminr5r,r+_evaluate_maxlevelr;rJr7 evaluate_fromr8rHr%rrrra~s$        zBaseNode._evaluate_maxlevelcCs<|jdur|jS|jdd|_|jdur6|jdd|_|jS)Nr+)rcr_)r/rar7rrrr.s   zBaseNode.maxlevelcCs|j|j dSr:)rrXrdrrrr8szBaseNode.node_namecCs |j|jkr|StddS)a2 Decompose node data creating DWT coefficients subnodes. Performs Discrete Wavelet Transform on the `~BaseNode.data` and returns transform coefficients. Note ---- Descends to subnodes and recursively calls `~BaseNode.reconstruct` on them. z$Maximum decomposition level reached.N)r%r. _decomposerUrdrrr decomposes zBaseNode.decomposecCs tdSr:r?rdrrrreszBaseNode._decomposeFcCs|js |jS||S)aX Reconstruct node from subnodes. Parameters ---------- update : bool, optional If True, then reconstructed data replaces the current node data (default: False). Returns: - original node data if subnodes do not exist - IDWT of subnodes otherwise. )has_any_subnoder1 _reconstruct)r7updaterrr reconstructszBaseNode.reconstructcCs tdSr:r?rdrrrrhszBaseNode._reconstructcCs<||||}|dur8|r8|js8|||}|S)a Returns subnode or None (see `decomposition` flag description). Parameters ---------- part : Subnode name decompose : bool, optional If the param is True and corresponding subnode does not exist, the subnode will be created using coefficients from the DWT decomposition of the current node. (default: True) N)rDrEis_emptyrf)r7r>rfsubnoderrr get_subnodes    zBaseNode.get_subnodecCsdt|}t|tr@zd|}Wnty>t|Yn0t|tr|jdurpt||j|jkrpt d|r| |d|jd||jdS|Snt|dS)a Find node represented by the given path. Similar to `~BaseNode.get_subnode` method with `decompose=True`, but can access nodes on any level in the decomposition tree. Parameters ---------- path : str String composed of node names. See `Node.node_name` and `Node2D.node_name` for node naming convention. Notes ----- If node does not exist yet, it will be created by decomposition of its parent node. zMInvalid path parameter type - expected string or tuple of strings but got %s.r*NzPath length is out of range.rT) type isinstancer\rV TypeErrorstrr.r[rX IndexErrorrm)r7rerrmsgrrr __getitem__s&     zBaseNode.__getitem__cCst|tr|jdur:t|jt||j|jkr:td|r||d|jd}|dur||d|jd||d|jd}||||jd<qt|t rt |j |_ n t ||_ t |}|j j|kr|j ||_ ntdt|dS)a Set node or node's data in the decomposition tree. Nodes are identified by string `path`. Parameters ---------- path : str String composed of node names. data : array or BaseNode subclass. NzPath length out of range.rFz9Invalid path parameter type - expected string but got %s.)rorqr.r[rrXrrrmrCrr3r4r1rdtypeZastyperprn)r7rr1rlrurrr __setitem__s*    zBaseNode.__setitem__cCs.||}|j}d|_|r*|jr*||jdS)z Remove node from the tree. Parameters ---------- path : str String composed of node names. N)r+r8rN)r7rrHr+rrr __delitem__6s  zBaseNode.__delitem__cCs |jduSr:r1rdrrrrkIszBaseNode.is_emptycCs$|jD]}||durdSqdS)NTF)r;rEr=rrrrgMs zBaseNode.has_any_subnodecs$gfdd}|j|dS)z Returns leaf nodes. Parameters ---------- decompose : bool, optional (default: True) cs<|j|jkr |js |dSs8|js8|dSdSNFT)r%r.rkappendrgrHrfresultrrcollect_s   z(BaseNode.get_leaf_nodes..collectrf)walk)r7rfr~rr|rget_leaf_nodesTs zBaseNode.get_leaf_nodesrcCsb|dur i}||g|Ri|r^|j|jkr^|jD](}|||}|dur4|||||q4dS)as Traverses the decomposition tree and calls ``func(node, *args, **kwargs)`` on every node. If `func` returns True, descending to subnodes will continue. Parameters ---------- func : callable Callable accepting `BaseNode` as the first param and optional positional and keyword arguments args : func params kwargs : func keyword params decompose : bool, optional If True (default), the method will also try to decompose the tree up to the `maximum level `. N)r%r.r;rmrr7funcargsrGrfr>rlrrrrjs"  z BaseNode.walkcCsb|dur i}|j|jkrH|jD](}|||}|dur|||||q||g|Ri|dS)a Walk tree and call func on every node starting from the bottom-most nodes. Parameters ---------- func : callable Callable accepting :class:`BaseNode` as the first param and optional positional and keyword arguments args : func params kwargs : func keyword params decompose : bool, optional (default: False) N)r%r.r;rm walk_depthrrrrrs   zBaseNode.walk_depthcCs|jdt|jS)Nz: )rrqr1rdrrr__str__szBaseNode.__str__)NT)NT)r+)F)T)F)rNT)rNT)!__name__ __module__ __qualname____doc__rXr;r9r6rCrIrEr<rNrDpropertyr]rar.r8rfrerjrhrmrtrvrwrkrgrrrrrrrrrsD       &%     rc@s>eZdZdZdZdZeefZdZd ddZd d Z d d Z dS)rz WaveletPacket tree node. Subnodes are called `a` and `d`, just like approximation and detail coefficients in the Discrete Wavelet Transform. rrrNTcCs|jt|||dSN)rFr>r1rB)rIrrArrrrCs zNode._create_subnodecCs|jrLd\}}||jdur,||j|||jdur||j|n8t|j|j|j|j d\}}||j|||j|||j||jfS)zq See also -------- dwt : for 1D Discrete Wavelet Transform output coefficients. NNNaxis) rkrEArCDrr1r,r-r0)r7data_adata_drrrres zNode._decomposecCsd\}}||j||j}}|dur2|}|durB|}|dur\|dur\tdnTt|||j|j|jd}|j dur|j |j kr|t dd|j D}|r||_ |SdS)Nrz>Node is a leaf node and cannot be reconstructed from subnodes.rcSsg|] }t|qSrslicerszrrrr r!z%Node._reconstruct..) rErrrjrUrr,r-r0r2r5r\r1)r7rirrZnode_aZnode_drecrrrrhs    zNode._reconstruct)NT) rrrrrrr;rXrCrerhrrrrrs rc@sReZdZdZdZdZdZdZeeeefZdZ dd d Z d d Z d dZ ddZ dS)rz WaveletPacket tree node. Subnodes are called 'a' (LL), 'h' (HL), 'v' (LH) and 'd' (HH), like approximation and detail coefficients in the 2D Discrete Wavelet Transform rhvrrNTcCs|jt|||dSr)rIrrArrrrCs zNode2D._create_subnodecCs|jrd\}}}}n"t|j|j|j|jd\}\}}}||j|||j|||j |||j || |j| |j | |j| |j fS)q See also -------- dwt2 : for 2D Discrete Wavelet Transform output coefficients. NNNNr0) rkrr1r,r-r0rCLLLHHLHHrE)r7data_lldata_lhdata_hldata_hhrrrres zNode2D._decomposec Csd\}}}}||j||j||j||jf\}}}} |durP|}|dur`|}|durp|}| dur| }|dur|dur|dur|durtd|jnf||||ff} t| |j |j |j d} |j dur| j |j kr| tdd|j D} |r| |_| SdS)NrSTree is missing data - all subnodes of `%s` node are None. Cannot reconstruct node.rcSsg|] }t|qSrrrrrrr r!z'Node2D._reconstruct..)rErrrrrjrUrrr,r-r0r2r5r\r1) r7rirrrrZnode_llZnode_lhZnode_hlZnode_hhcoeffsrrrrrhs@    zNode2D._reconstructcsL|jd|jd|jd|jdidfdd|Ddfdd|DfS) NhhZhlZlhZllr*csg|]}|dqS)rrrRZexpanded_pathsrrr +r!z)Node2D.expand_2d_path..csg|]}|dqSrWrrRrrrr ,r!)rrrrrV)r7rrrrexpand_2d_path$szNode2D.expand_2d_path)NT)rrrrrrrrr;rXrCrerhrrrrrrs   rcsleZdZdZfddZddZddZdd Zd d Zd d Z dddZ dddZ ddZ ddZ ZS)r aV WaveletPacket tree node. Unlike Node and Node2D self.PARTS is a dictionary. For 1D: self.PARTS has keys 'a' and 'd' For 2D: self.PARTS has keys 'aa', 'ad', 'da', 'dd' For 3D: self.PARTS has keys 'aaa', 'aad', 'ada', 'daa', ..., 'ddd' Parameters ---------- parent : Parent node. If parent is None then the node is considered detached (ie root). data : 1D or 2D array Data associated with the node. 1D or 2D numeric array, depending on the transform type. node_name : string A name identifying the coefficients type. See `Node.node_name` and `Node2D.node_name` for information on the accepted subnodes names. ndim : int The number of data dimensions. ndim_transform : int The number of dimensions that are to be transformed. csXtt|j|||d||_t|_td|jD]}d|jd|<q2||_||_ dS)N)r+r1r8)adr*) superr r9rXr r;r rVndimndim_transform)r7r+r1r8rrkey __class__rrr9JszNodeND.__init__cCsdSr:rrdrrrr6TszNodeND._init_subnodescCs |j|Sr:)r;r=rrrrEYszNodeND._get_nodecCs ||jvrtd||j|<dS)Nz invalid part)r;rUrLrrrr<\s zNodeND._set_nodecCs||ddSr:rMr=rrrrNaszNodeND._delete_nodecCs8||jvr4tddddt|jD|fdS)NrOrPcss|]}d|VqdSrQrrRrrrrThr!z-NodeND._validate_node_name..)r;rUrVlistkeysr=rrrrDds   zNodeND._validate_node_nameNTcCs|jt||||j|jdS)N)rFr>r1rBrr)rIr rrrArrrrCjs zNodeND._create_subnoder+cCs|dvs J|jdur|jS|jdur@|jtt|jj|jS|dkr`|jdur|j|Sn>|dkr|j D]*\}}|durr||}|durr|SqrdSr^) r/r1r%rr`r5r,r+rar;itemsrbrrrraos"      zNodeND._evaluate_maxlevelcsjjrddjD}ntjjjjd}|D]\}} ||q<fddjDS)rcSsi|] }|dqSr:rrrrrr r!z%NodeND._decompose..rc3s|]}|VqdSr:)rErrdrrrTr!z$NodeND._decompose..) rkr;rrr1r,r-r0rrC)r7Zcoefsrr1rrdrres zNodeND._decomposecCsdd|jD}d}|jD]*}||}|dur"|d7}|||<q"|dkrftd|jn$t||j|j|j d}|r||_ |SdS)NcSsi|] }|dqSr:rrrrrrr!z'NodeND._reconstruct..rrrr) r;rrErjrUrrr,r-r0r1)r7rirZnnodesrrHrrrrrhs" zNodeND._reconstruct)NT)r+)rrrrr9r6rEr<rNrDrCrarerh __classcell__rrrrr /s   r csBeZdZdZdfdd ZddZdfd d Zdd dZZS)ra Data structure representing Wavelet Packet decomposition of signal. Parameters ---------- data : 1D ndarray Original data (signal) wavelet : Wavelet object or name string Wavelet used in DWT decomposition and reconstruction mode : str, optional Signal extension mode for the `dwt` and `idwt` decomposition and reconstruction functions. maxlevel : int, optional Maximum level of decomposition. If None, it will be calculated based on the `wavelet` and `data` length using `pywt.dwt_max_level`. axis : int, optional The axis to transform. symmetricNr#cstt|d|dt|ts&t|}||_||_||_|durt |}|jdkrb|j|j |_d|jkrz|j ksnt d|j |_ |durt|j |j|j}nd|_ ||_dS)Nr*rz!Axis greater than data dimensions)rrr9rorr,r-r0r3r4rrUr5 data_sizerr/)r7r1r,r-r.rrrrr9s"   zWaveletPacket.__init__cCst|j|j|j|jffSr:)rr1r,r-r.rdrrr __reduce__szWaveletPacket.__reduce__TcsT|jrNtt||}|jdur@|j|jkr@|dd|jD}|rJ||_|S|jS)a  Reconstruct data value using coefficients from subnodes. Parameters ---------- update : bool, optional If True (default), then data values will be replaced by reconstruction values, also in subnodes. NcSsg|] }t|qSrrrrrrr r!z-WaveletPacket.reconstruct..)rgrrrjrr5r1r7rir1rrrrjs zWaveletPacket.reconstructnaturalcs|dvrtd||jkr.td|jgfdd}|j||d|dkrZS|dkrtd d Dt}fd d |DStd |dS)a Returns all nodes on the specified level. Parameters ---------- level : int Specifies decomposition `level` from which the nodes will be collected. order : {'natural', 'freq'}, optional - "natural" - left to right in tree (default) - "freq" - band ordered decompose : bool, optional If set then the method will try to decompose the data up to the specified `level` (default: True). Notes ----- If nodes at the given level are missing (i.e. the tree is partially decomposed) and `decompose` is set to False, only existing nodes will be returned. Frequency order (``order="freq"``) is also known as as sequency order and "natural" order is sometimes referred to as Paley order. A detailed discussion of these orderings is also given in [1]_, [2]_. References ---------- ..[1] M.V. Wickerhauser. Adapted Wavelet Analysis from Theory to Software. Wellesley. Massachusetts: A K Peters. 1994. ..[2] D.B. Percival and A.T. Walden. Wavelet Methods for Time Series Analysis. Cambridge University Press. 2000. DOI:10.1017/CBO9780511841040 rfreqInvalid order: {}KThe level cannot be greater than the maximum decomposition level value (%d)cs|jkr|dSdSryr%rzr{r%r}rrr~s  z(WaveletPacket.get_level..collectrrrcss|]}|j|fVqdSr:)rrrHrrrrT'r!z*WaveletPacket.get_level..csg|]}|vr|qSrrr)r}rrr )r!z+WaveletPacket.get_level..zInvalid order name - %s.N)rUformatr.rdictr))r7r%orderrfr~r&rrr get_levels " zWaveletPacket.get_level)rNr#)T)rT rrrrr9rrjrrrrrrrsrcsBeZdZdZdfdd ZddZdfd d Zdd dZZS)r a Data structure representing 2D Wavelet Packet decomposition of signal. Parameters ---------- data : 2D ndarray Data associated with the node. wavelet : Wavelet object or name string Wavelet used in DWT decomposition and reconstruction mode : str, optional Signal extension mode for the `dwt` and `idwt` decomposition and reconstruction functions. maxlevel : int Maximum level of decomposition. If None, it will be calculated based on the `wavelet` and `data` length using `pywt.dwt_max_level`. axes : 2-tuple of ints, optional The axes that will be transformed. smoothNr#cstt|ddt|ts&t|}||_||_t||_t t |jdkrXt ddurt jdkr|t dj|_fdd|jD}|durtt||j}nd|_||_dS)Nr*zExpected two unique axes.z8WaveletPacket2D requires data with 2 or more dimensions.csg|]}j|qSrr5raxrxrrr Sr!z,WaveletPacket2D.__init__..)rr r9rorr,r-r\r0r[r3uniquerUr4rr5rrr`r/)r7r1r,r-r.r0transform_sizerrxrr9Bs(    zWaveletPacket2D.__init__cCst|j|j|j|jffSr:)r r1r,r-r.rdrrrrZszWaveletPacket2D.__reduce__TcsT|jrNtt||}|jdur@|j|jkr@|dd|jD}|rJ||_|S|jS)" Reconstruct data using coefficients from subnodes. Parameters ---------- update : bool, optional If True (default) then the coefficients of the current node and its subnodes will be replaced with values from reconstruction. NcSsg|] }t|qSrrrrrrr kr!z/WaveletPacket2D.reconstruct..)rgrr rjrr5r1rrrrrj^s zWaveletPacket2D.reconstructrc s|dvrtd|jkr.tdjgfdd}j||d|dkrifdd DD]\\}}}||i|<qltd d d }fd d |DgD]fdd |DqS)a Returns all nodes from specified level. Parameters ---------- level : int Decomposition `level` from which the nodes will be collected. order : {'natural', 'freq'}, optional If `natural` (default) a flat list is returned. If `freq`, a 2d structure with rows and cols sorted by corresponding dimension frequency of 2d coefficient array (adapted from 1d case). decompose : bool, optional If set then the method will try to decompose the data up to the specified `level` (default: True). Notes ----- Frequency order (``order="freq"``) is also known as as sequency order and "natural" order is sometimes referred to as Paley order. A detailed discussion of these orderings is also given in [1]_, [2]_. References ---------- ..[1] M.V. Wickerhauser. Adapted Wavelet Analysis from Theory to Software. Wellesley. Massachusetts: A K Peters. 1994. ..[2] D.B. Percival and A.T. Walden. Wavelet Methods for Time Series Analysis. Cambridge University Press. 2000. DOI:10.1017/CBO9780511841040 rrrcs|jkr|dSdSryrr{rrrr~s  z*WaveletPacket2D.get_level..collectrrcsg|]}|j|fqSr)rrrrdrrr sz-WaveletPacket2D.get_level..lrr(csg|]}|vr|qSrrr)nodesrrr r!csg|]}|vr|qSrrr)rowrrr r!)rUrr.r setdefaultr)rz) r7r%rrfr~Zrow_pathZcol_pathrHr&r)r%rr}rr7rrqs.   zWaveletPacket2D.get_level)rNr)T)rTrrrrrr .sr cs:eZdZdZd fdd Zd fdd Zd d d ZZS)r a Data structure representing ND Wavelet Packet decomposition of signal. Parameters ---------- data : ND ndarray Data associated with the node. wavelet : Wavelet object or name string Wavelet used in DWT decomposition and reconstruction mode : str, optional Signal extension mode for the `dwt` and `idwt` decomposition and reconstruction functions. maxlevel : int, optional Maximum level of decomposition. If None, it will be calculated based on the `wavelet` and `data` length using `pywt.dwt_max_level`. axes : tuple of int, optional The axes to transform. The default value of `None` corresponds to all axes. rNc sJdur|durtd|dur,tj}nt|r<|f}t|}tt|t|krbtdt|}durtjdkrtdj}nt|}t t | dd||t |t st |}||_||_||_||_dur:jt|krtdj|_fdd|jD}|dur@tt||j}nd|_||_dS) Nz'If data is None, axes must be specifiedzExpected a set of unique axes.rzdata must be at least 1Dr*z9The number of axes exceeds the number of data dimensions.csg|]}j|qSrrrrxrrr r!z,WaveletPacketND.__init__..)rUr$rr3Zisscalarr\r[rr4rr r9rorr,r-r0rr5rrr`r/) r7r1r,r-r.r0rrrrrxrr9sD       zWaveletPacketND.__init__TcsT|jrNtt||}|jdur@|j|jkr@|dd|jD}|rJ||_|S|jS)rNcSsg|] }t|qSrrrrrrr r!z/WaveletPacketND.reconstruct..)rgrr rjrr5r1rrrrrjs zWaveletPacketND.reconstructcs<|jkrtd|jgfdd}|j||dS)ah Returns all nodes from specified level. Parameters ---------- level : int Decomposition `level` from which the nodes will be collected. decompose : bool, optional If set then the method will try to decompose the data up to the specified `level` (default: True). rcs|jkr|dSdSryrr{rrrr~s  z*WaveletPacketND.get_level..collectr)r.rUr)r7r%rfr~rrrrs zWaveletPacketND.get_level)rNN)T)T)rrrrr9rjrrrrrrr s ,r )rr)!r __future__rrr__all__ itertoolsr collectionsr numpyr3Z_extensions._pywtrrZ_dwtrrrZ _multidimrrrrr)objectrrrr rr r rrrrs*   ;Q~