数字信号处理Day2－小波基与规范正交化

我们有这么一张灰度图64*64

我们能够定义出4096个基，各自是某一位是0其它是1，在这样的情况下，假设我们传输图片，那么就相当于传输原始数据

如果传到一半，网络坏了。

于是，我们得到数字信号处理Day2－小波基与规范正交化

我们能够计算原图像和这图像的差距

error = I - I_approx;

distance = sqrt(sum(sum(error.*error)))

distance =

713.5559

如果我们使用Haar基

数字信号处理Day2－小波基与规范正交化

function h = haar(N)

    % Assume N is power of 2

    h = zeros(N,N);

    h(1,1:N) = ones(1,N)/sqrt(N);

    for k = 1:N-1

        p = fix(log(k)/log(2));

        q = k-(2^p);

        k1 = 2^p;

        t1 = N/k1;

        k2 = 2^(p+1);

        t2 = N/k2;

        for i=1:t2

            h(k+1,i+q*t1) = (2^(p/2))/sqrt(N);

            h(k+1,i+q*t1+t2) = -(2^(p/2))/sqrt(N);

        end

    end

然后，我们看一下仅仅接受一半的图片能够恢复成什么样子，这个东西能够用到图片压缩（PCA降维）

clear all

close all

clc

% Load image

I = double(imread('camera.jpg'));

% Arrange image in column vector

I = I(:);

% Generate Haar basis vector (rows of H)

H = haar(4096);

% Project image on the new basis

I_haar = H*I;

% Remove the second half of the coefficient

I_haar(2049:4096) = 0;

% Recover the image by inverting change of basis

I_haar = H'*I_haar;

% Rearrange pixels of the image

I_haar = reshape(I_haar, 64, 64);

% Display image

figure

imshow(I_haar,[]);

数字信号处理Day2－小波基与规范正交化

效果还能够

error = I - I_haar;

distance = sqrt(sum(sum(error.*error)))

distance =

350.6765

Haar基是正交基

怎样将普通基变成正交基呢

线性代数中的拉格慕－施密特正交化即可

function E=gs_orthonormalization(V)

% V is a matrix where each column contains the vectors spanning the space of which we want to compute the orthonormal base

% E is a matrix where each column contains an ortho-normal vector of the base of the space

[numberLines,numberColumns]=size(V);

% U is a matrix containing the orthogonal vectors (non normalized)

U=zeros(numberLines,numberColumns);

for indexColumn=1:numberColumns

    U(:,indexColumn)=V(:,indexColumn);

    for index=1:(indexColumn-1)

        R(index,indexColumn) =U(:,index)'*V(:,indexColumn);

        U(:,indexColumn)=U(:,indexColumn)-R(index,indexColumn)*U(:,index);

    end

    R(indexColumn,indexColumn)=norm(U(:,indexColumn));

    U(:,indexColumn)=U(:,indexColumn)./R(indexColumn,indexColumn);

end

E=U;

return

秒客网

数字信号处理Day2－小波基与规范正交化

相关文章