Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/ShozenD/AutocorrelationShell.jl
Julia implementation of Autocorrelation Wavelets.
https://github.com/ShozenD/AutocorrelationShell.jl
Last synced: 4 days ago
JSON representation
Julia implementation of Autocorrelation Wavelets.
- Host: GitHub
- URL: https://github.com/ShozenD/AutocorrelationShell.jl
- Owner: ShozenD
- License: other
- Created: 2020-09-30T11:44:23.000Z (about 4 years ago)
- Default Branch: master
- Last Pushed: 2021-04-18T05:33:43.000Z (over 3 years ago)
- Last Synced: 2024-10-10T02:15:22.964Z (about 1 month ago)
- Language: Julia
- Homepage:
- Size: 42.7 MB
- Stars: 6
- Watchers: 2
- Forks: 1
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
[](https://ShozenD.github.io/AutocorrelationShell.jl/stable)
[](https://github.com/ShozenD/AutocorrelationShell.jl/actions)
[](https://gitlab.com/BoundaryValueProblems/autocorrelation-shell/-/commits/master)
[](https://codecov.io/gh/ShozenD/AutocorrelationShell.jl)
**Warning**: `AutocorrelationShell.jl` is deprecated and will be part of `WaveletsExt.jl` from April 23, 2021. Please follow the following steps to remove this package from your REPL and migrate to `WaveletsExt.jl`.
```
julia> ]
(@v1.x) pkg> rm AutocorrelationShell
(@v1.x) pkg> update
(@v1.x) pkg> add WaveletsExt
```
This package is a [Julia](https://github.com/JuliaLang/julia) implementation of Autocorrelation wavelet transforms (ACWT). AC wavelets are a special case of the stationary wavelet transform. This package includes the 1D autocorrelation wavelet transform, 2D autocorrelation wavelet transform, and autocorrelation wavelet packet transform.
Signal representations using autocorrelation wavelets are redundant and non-orthogonal. Some desirable properties of autocorrelation wavelet transforms are symmetry without losing vanishing moments, edge detection and characterization capabilities, and shift invariance. Autocorrelation wavelets can be used as a tool for data analysis such as time series analysis and image analysis.
## Authors
This package was first translated from Matlab by Rishi Subramanian, and was extended by Christina Chang, and currently maintained by Shozen Dan under the supervision of Professor Naoki Saito at University of California, Davis.
## Installation
The package is part of the official Julia Registry. It can be install via the Julia REPL
```
(@1.x) pkg> add AutocorrelationShell
```
or
```
julia> using Pkg; Pkg.add("AutocorrelationShell")
```
## Usage
Load `AutocorrelationShell.jl` with the `Wavelets.jl` package
```{julia}
using Wavelets, AutocorrelationShell
```
## 1D Autocorrelation Wavelet Transform
```{julia}
# Forward 1D Autocorrelation Wavelet Transform
y = acwt(x, wavelet(WT.db4))
# Inverse 1D Autocorrelation Wavelet Transform
iacwt(y)
```
### Example
Perform forward autocorrelation wavelet transform on the vector x
```{julia}
# Decompose signal
y = acwt(x, wavelet(WT.db4))
# Display decomposition
wiggle(y)
```
Result:
## 2D Autocorrelation Wavelet Transform
```{julia}
# Forward 2D Autocorrelation Wavelet Transform
y = acwt(img, wavelet(WT.db4))
```
The `acwt` function performs a forward wavelet transformation on 2D signals such as images. It returns a 4 dimensional tensor with the dimensions (num_row, num_col, levels_of_decomp_row, levels_of_decomp_col).
```{julia}
# Inverse 2D Autocorrelation Wavelet Transform
iacwt(y)
```
The `iacwt` function is the inverse function of `acwt`. It takes an array of autocorrelation wavelet coefficients and reconstructs the original signal.
### Example
```{julia}
X = load(../test/pictures/boat.jpg)
X = Float64.(Gray.(X))
Y = acwt(X, wavelet(WT.db4))
# Revert to original signal
Z = iacwt(Y)
```
## Autocorrelation Wavelet Packet Transform
```{julia}
# Autocorrelation Wavelet Packet Transform
acwpt(x, wavelet(WT.db4))
```
The `acwpt` function computes the autocorrelation wavelet packet transform for 1 dimensional signal. It returns a binary tree where the root node contains the original signal, and each child node contains a vector of 1 dimensional autocorrelation wavelet transform coefficients.
### Example
```{julia}
X₁ = randn(4); # length 4 random signal
y = acwpt(X₁, wavelet(WT.db4))
```