【数模】奇异值分解SVD和图形处理

Shier833_Ww

已于 2023-09-06 22:47:33 修改

阅读量1.6k

点赞数 1

分类专栏：数学建模文章标签： matlab 算法数学建模图像处理奇异值分解 SVD

于 2023-08-05 13:12:15 首次发布

本文链接：https://blog.csdn.net/shie_ww/article/details/131990287

版权

数学建模专栏收录该内容

13 篇文章 6 订阅

订阅专栏

介绍奇异值分解在图形压缩中的运用，并将简单介绍下Matlab对于图形和视频的处理

一、奇异值分解介绍

1.1 基本概念

奇异值分解(Singular Value Decomposition，以下简称SVD)是线性代数中一种重要的矩阵分解：
- U和V都是正交矩阵
- ∑是奇异值矩阵（和A的大小相同，其主对角元素是从大到小排列的，这些对角元素称为奇异值，其他位置元素为0）

1.2 求∑的步骤

在这里插入图片描述

1.2.1 U的计算

在这里插入图片描述

1.2.2 V的计算

在这里插入图片描述

1.2.3 ∑的计算

在这里插入图片描述

1.3 利用SVD对数据进行"降维"

举例：
此处的“降维”：更准确而言是使得矩阵的秩减小，矩阵大小并未改变。
保留原矩阵的特征比例

二、Matlab对奇异值分解

矩阵压缩

2.1 Matlab自带分解：[U,S,V] = svd(A)

在这里插入图片描述

% PPT_example.m文件
A = [4 0 1 6;0 0 5 1;2 1 3 2] % A : 3*4
[U,S,V] = svd(A)
% 注意：U*S*(V的转置) == A
% U:3*3 S：3*4 V ： 4*4
U*S*V' - A 
% 因为浮点数计算的缘故，所以会有非常微小的偏差
% 1.0e-14 * 0.2665 = 0.000000000000002665
U(:,1:2)*S(1:2,1:2)*V(:,1:2)'
U(:,1:3)*S(1:3,1:3)*V(:,1:3)'  % 就是A

2.2 自定义的mysvd函数

取自清风老师自定义的，不同于Matlab自带的：使用奇异值分解将矩阵A压缩到指定的特征比例

% mysvd.m文件
function [compress_A] = mysvd(A, ratio)
% 输入变量
	% A：要压缩的m*n维的矩阵
	% ratio：(至少)要保留原矩阵的特征比例（100%表示不压缩）
% 输出变量
	% compress_A：压缩后的矩阵

在这里插入图片描述

三、Matlab图形的处理

图片压缩

3.1 RGB模式

RGB色彩（三原色）：R代表Red（红色），G代表Green（绿色），B代表Blue（蓝色）。
- 计算机定义颜色时R、G、 B三种成分的取值范围是0‐255，0表示没有刺激量，255表示刺激量达最大值。
- R、G、B均为255时就合成了白光，R、G、B均为0时就形成了黑色。
- 灰色图片：R、G、B三原色完全相同；否则称为彩色图片。

3.2 图片压缩的函数

函数作用：利用SVD函数对图形进行压缩
输入变量
- photo_address：要压缩的图片存放的位置（建议输入完整的路径）
- save_address：将压缩后的图片保存的位置（建议输入完整的路径）
- ratio：要保留原矩阵的特征比例（100%表示不压缩）
- greycompress: 如果该值等于1，则会彩色的原图片转换为灰色图片后再压缩；默认值为0，表示不进行转换
输出变量：无（不需要输出，因为函数运行过程中已经将图片保存了）

% photo_compress.m文件
function []= photo_compress(photo_address, save_address,
ratio, greycompress)

3.3 SVD压缩后的效果

在这里插入图片描述

四、视频分离为图片

4.1 `video2photo.m`代码文件

Matlab对视频的要求较高，有的原视频属性详细情况不带帧数，可以采用转化器转个视频格式（如迅捷视频转换器），否则可能会报如下错误：
- 错误使用 VideoReader/init (line 619)
- 无法确定所需的编解码器
video中的对象讲解

%% 读取视频
video_file='迅捷视频转换器转换后的新闻联播.mp4';
video=VideoReader(video_file);
frame_number = video.NumberOfFrames; %视频的总帧数

%% 分离图片
for i=1:30:frame_number % 这里演示的是每30帧数保存一次（如果每一帧都全部保存的话所要花费的时间很长）
	image_name=strcat('C:\Users\hc_lzp\Desktop\数学建模视频录制\第13讲.奇异值分解\代码和例题数据\将视频转成图片\图片\image_',num2str(i),'.jpg');
	Photo=read(video,i); %读出所在帧的图片对象
	imwrite(Photo,image_name); %将图片保存到指定的位置
end

4.2 分离效果

4.3 批量处理文件夹内的所有图片

在这里插入图片描述

folder_name = 'C:\Users\hc_lzp\Desktop\数学建模视频录制\第13讲.奇异值分解\代码和例题数据\压缩文件夹内的所有图片\葫芦娃七兄弟';
dirOutput=dir(fullfile(folder_name, '*.jpg'));
files={dirOutput.name};
n = length(files);
ratio = 0.9;
for i = 1:n
	photoname = files(i);
	name = photoname{1};
	photo_address = fullfile(folder_name, name);
	disp(photo_address)
	save_address = fullfile(folder_name, strcat('compress_',name));
	photo_compress(photo_address, save_address, ratio)
end

五、SVD的评价和应用

5.1 SVD（奇异值分解）评价

优点：简化数据，去除噪声点，对数据降维（减少秩）；
缺点：数据的转换可能难以理解；
适用于数据类型：数值型。

通过SVD对数据的处理，可对原始数据进行精简，实际上是去除了噪声和冗余信息，以此达到了优化数据的目的。

5.2 SVD的另外两个重要应用

潜在语义索引：最早的SVD应用之一就是信息检索，称利用SVD的方法为潜在语义检索（LSI）或隐形语义分析（LSA）
- 有兴趣可以去阅读吴军老师的《数学之美》。
推荐系统：SVD的另一个应用就是推荐系统，较为先进的推荐系统先利用SVD从数据中构建一个主题空间，然后再在该空间下计算相似度，以此提高推荐的效果。