{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Illustrate wavelet \"filters\"\n",
"2019-04-02, Jeff Fessler, University of Michigan"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {},
"outputs": [],
"source": [
"using MIRT\n",
"using LinearAlgebra: I"
]
},
{
"cell_type": "code",
"execution_count": 31,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"isapprox(T' * T, I) = true\n"
]
},
{
"data": {
"image/svg+xml": [
"\n",
"\n"
]
},
"execution_count": 31,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"M,N = 16,8\n",
"T = Aodwt((M,N), level=1)[1]\n",
"T = Matrix(T)\n",
"@show isapprox(T'T, I)\n",
"jim(T)"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
"\n"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"tmp = reshape(T', M, N, :)\n",
"jim(tmp, padval=1, ncol=N)"
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"isapprox(T2' * T2, I) = true\n"
]
},
{
"data": {
"image/svg+xml": [
"\n",
"\n"
]
},
"execution_count": 33,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"T2 = Aodwt((M,N), level=2)[1]\n",
"T2 = Matrix(T2)\n",
"@show isapprox(T2'T2, I)\n",
"jim(T2)"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {},
"outputs": [
{
"data": {
"image/svg+xml": [
"\n",
"\n"
]
},
"execution_count": 34,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"tmp = reshape(T2', M, N, :)\n",
"jim(tmp, padval=1, ncol=N)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 1.1.0",
"language": "julia",
"name": "julia-1.1"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.1.0"
}
},
"nbformat": 4,
"nbformat_minor": 2
}