clear; clc; close all; A=imread('1.jpg'); I=rgb2gray(A);%将彩色图像转为灰度图像 [M,N] = size(I);%将图像转为二维矩阵 I1 = I(:);%转为一维向量 P = zeros(1,256); %获取各符号的概率; for i = 0:255 P(i+1) = length(find(I1 == i))/(M*N); end k = 0:255; dict = huffmandict(k,P); %根据灰度级k和概率数组P生成Huffman字典 enco = huffmanenco(I1,dict); %哈夫曼编码 deco = huffmandeco(enco,dict); %哈夫曼解码 Ide = col2im(deco,[M,N],[M,N],'distinct'); %把向量重新转换成图像块; subplot(1,2,1);imshow(I);title('原始图像'); subplot(1,2,2);imshow(uint8(Ide));title('解码图像'); B=length(enco);%原始图像比特长度 sumcode=length(deco);%编码后比特长度 CR=sumcode/B;%计算压缩率 disp(['原始图像 Bit: ',num2str(B),' bit']); disp(['压缩图像 Bit: ',num2str(sumcode),' bit']); disp(['压缩率: ',num2str(CR)]);

时间: 2024-02-07 16:02:52 浏览: 145
好的,我进行修改后的代码如下: ```matlab clear; clc; close all; img = imread('1.jpg'); grayImg = rgb2gray(img); %将彩色图像转为灰度图像 [M, N] = size(grayImg); %将图像转为二维矩阵 I1 = grayImg(:); %转为一维向量 probabilities = zeros(1,256); %获取各符号的概率 for i = 0:255 probabilities(i+1) = length(find(I1 == i))/(M*N); end symbols = 0:255; dict = huffmandict(symbols, probabilities); %根据灰度级symbols和概率数组probabilities生成Huffman字典 encoded = huffmanenco(I1, dict); %哈夫曼编码 decoded = huffmandeco(encoded, dict); %哈夫曼解码 decodedImg = col2im(decoded, [M, N], [M, N], 'distinct'); %把向量重新转换成图像块 subplot(1,2,1);imshow(grayImg); title('原始图像'); subplot(1,2,2);imshow(uint8(decodedImg)); title('解码图像'); originalBitLength = length(I1)*8; %原始图像比特长度 compressedBitLength = length(encoded); %编码后比特长度 compressionRatio = originalBitLength/compressedBitLength; %计算压缩率 disp(['原始图像 Bit: ',num2str(originalBitLength),' bit']); disp(['压缩图像 Bit: ',num2str(compressedBitLength),' bit']); disp(['压缩率: ',num2str(compressionRatio)]); ```
阅读全文

相关推荐

clear;clc;close all; img=imread('flower.tif'); gray=rgbimage2gray(img); %灰度化 %加入噪声 gray_noise=imnoise(gray,'salt & pepper',0.2); % 自适应中值滤波 f1 = adaptive_median_filter(gray_noise,11); if(size(img, 3) == 3) % Check if the image is a truecolor image f1 = gray2rgb(f1,img); end figure('color',[1,1,1]); subplot(221) imshow(img) title("原图") subplot(222) imshow(gray_noise) title("gray with noise") subplot(224) imshow(f1); title("自适应中值滤波") function f = adaptive_median_filter (g, Smax) % 判断邻域是否合理 if (Smax <= 1) || (Smax/2 == round(Smax/2)) || (Smax ~= round(Smax)) error ('SMAX must be an odd integer > 1.') end % f = g; f(:) = 0; % 标记是否已处理过 alreadyProcessed = false (size(g)); % 开始自适应滤波 for k = 3:2:Smax zmin = ordfilt2(g, 1, ones(k, k),'symmetric'); zmax = ordfilt2(g, k * k, ones(k, k), 'symmetric'); zmed = medfilt2(g, [k k], 'symmetric'); % 判断是否进入进程B processUsingLevelB = (zmed > zmin) & (zmax > zmed) & ~alreadyProcessed; % 若g不是脉冲,保留原值 zB = (g > zmin) & (zmax > g); outputZxy = processUsingLevelB & zB; %若是脉冲,用Zmed替换 outputZmed = processUsingLevelB & ~zB; f (outputZxy) = g(outputZxy); f (outputZmed) = zmed(outputZmed); % 已处理记录 alreadyProcessed = alreadyProcessed | processUsingLevelB; % 是否退出 if all (alreadyProcessed (:)) break; end end % 大于窗口尺寸后,Zxy替换成Zmed输出 f (~alreadyProcessed) = zmed (~alreadyProcessed); end function img_gray=rgbimage2gray(img) % 灰度变换,公式:f(x,y)=0.2989R+ 0.5870G + 0.1140B img_gray = img(:,:,1)*0.2989+ img(:,:,2)*0.5870+ img(:,:,3)*0.1140; end function img_rgb=gray2rgb(img_gray,img) % 将灰度图像转化为RGB图像 img_rgb = zeros(size(img)); img_rgb(:,:,1) = img_gray; img_rgb(:,:,2) = img_gray; img_rgb(:,:,3) = img_gray; end带有下标的赋值维度不匹配。 出错 Untitled13>gray2rgb (line 75) img_rgb(:,:,1) = img_gray; 出错 Untitled13 (line 10) f1 = gray2rgb(f1,img);

clear all; close all; clc; Ia=imread('/Users/aa/Desktop/截屏2023-05-31 16.55.11.png'); I=rgb2gray(Ia); [m n]=size(I); I=double(I); w=fspecial('gaussian',[5 5]);%%高斯滤波 I=imfilter(I,w,'replicate'); imshow(uint8(I))%%sobel边缘检测 w=fspecial('sobel'); Iw=imfilter(I,w,'replicate');%求横边缘 w=w'; Ih=imfilter(I,w,'replicate');%求竖边缘 I=sqrt(Iw.^2+Ih.^2);%平方和在开方。 imshow(uint8(I))%%下面是非极大抑制 new_edge=zeros(m,n); for i=2:m-1%避开边缘像素值检测 for j=2:n-1 Mx=Iw(i,j); My=Ih(i,j); if My~=0 o=atan(Mx/My); %边缘的法线弧度 elseif My==0 && Mx>0 o=pi/2; else o=-pi/2; end %Mx处用My和img进行插值 adds=get_coords(o); %边缘像素法线一侧求得的两点坐标,插值需要 M1=My*I(i+adds(2),j+adds(1))+(Mx-My)*I(i+adds(4),j+adds(3)); %插值后得到的像素,用此像素和当前像素比较 adds=get_coords(o+pi); %边缘法线另一侧求得的两点坐标,插值需要 M2=My*I(i+adds(2),j+adds(1))+(Mx-My)*I(i+adds(4),j+adds(3)); %另一侧插值得到的像素,同样和当前像素比较 isbigger=(Mx*I(i,j)>M1)*(Mx*I(i,j)>=M2)+(Mx*I(i,j)<M1)*(Mx*I(i,j)<=M2); %如果当前点比两边点都大置1 if isbigger new_edge(i,j)=I(i,j); end end end imshow(uint8(new_edge)) %%下面是滞后阈值处理 up=120; %上阈值 low=100; %下阈值 set(0,'RecursionLimit',10000); %设置最大递归深度 for i=1:m for j=1:n if new_edge(i,j)>up &&new_edge(i,j)~=255 %判断上阈值 new_edge(i,j)=255; new_edge=connect(new_edge,i,j,low); end end end imshow(new_edge==255)

clc; clear; %shuzituxiang32 close all; pic=imread('C:\Users\NADPH\Desktop\未爆弹4.jpg'); gray_pic=rgb2gray(pic); binary_img=1-im2bw(pic,0.73); %灰度图转换成二值图像,直接进行灰度反转,让图形区域置1 for i=202:265 binary_img(188,i)=1; end fill_hole=imfill(binary_img,'holes'); fill_hole=bwareaopen(fill_hole,10); [B,L]=bwboundaries(fill_hole,'noholes'); out_result=regionprops(fill_hole,'Extent','Centroid','boundingbox'); %Extent:各连通区域像素点与最小边界像素点比值 centroids = cat(1, out_result.Centroid); %各连通区域质心 draw_rect=cat(1,out_result.BoundingBox); %各连通区域最小边界矩形 figure; subplot(2,2,1);imshow(pic,[]);title('原图'); subplot(2,2,2);imshow(binary_img,[]);title('二值化'); subplot(2,2,3);imshow(fill_hole,[]);title('区域填充'); subplot(2,2,4);imshow(fill_hole,[]);title('圆和矩形检测'); hold on; for i=1:size(out_result) rectangle('position',draw_rect(i,:),'EdgeColor','y','LineWidth',2); %绘出各连通区域最小边界矩形 if out_result(i).Extent>0.95 text(centroids(i,1)-20, centroids(i,2)-10,'矩形','Color','b','FontSize',9); text(centroids(i,1)-32, centroids(i,2)+10,num2str(out_result(i).Extent),'Color','b','FontSize',8); elseif out_result(i).Extent>0.76&&out_result(i).Extent<0.80 text(centroids(i,1)-20, centroids(i,2)-10,'圆形','Color','b','FontSize',9); text(centroids(i,1)-32, centroids(i,2)+10,num2str(out_result(i).Extent/(pi/4)),'Color','b','FontSize',8); end end hold on; for j=1:length(B) boundary=B{j}; plot(boundary(:,2),boundary(:,1),'r','LineWidth',2); end hold off;

zip

大家在看

recommend-type

自动化-ACS800变频器知识培训(0619)[1]专题培训课件.ppt

自动化-ACS800变频器知识培训(0619)[1]专题培训课件.ppt
recommend-type

贝叶斯分类.docx

适合初学者理解的贝叶斯分类的r代码,任何编程的背后都是理论的支撑,当初花了一天半编的该代码,欢迎指正。
recommend-type

IPC-7351 使用说明

IPC-7351 软件,零件封装库制作标准软件的中文使用说明。
recommend-type

子程序参数传递学习总结.docx

关于kuka编程知识的最新总结,全局子程序与局部子程序
recommend-type

三菱FX3U-485ADP-MB通讯三种变频器程序 已实现测试的变频器:施耐德ATV312, 三菱E700,台达VFD-M三款变

三菱FX3U-485ADP-MB通讯三种变频器程序 已实现测试的变频器:施耐德ATV312, 三菱E700,台达VFD-M三款变频器,支持rtu的协议的变频器都可实现。 需要硬件:FX3UPLC,FX3U-485ADP-MB通信扩展模块,施耐德ATV312变频器或台达vfd-m变频器或三菱E700变频器,fx3u-cnv-bd 。 通过modbus rtu通讯方式 ,可以实现控制正反转,启动停止,触摸屏直接频率设定,以及对频率电流,运行状态的监控。 反馈及时,无延迟,使用方便。 内容包含plc和触摸屏程序,参数设置,接线及教程。 这里有三种变频器程序,可以通过三菱FX3U-485ADP-MB通信扩展模块实现测试。已经测试过的变频器包括施耐德ATV312、三菱E700和台达VFD-M,只要支持rtu协议的变频器都可以使用。 为了实现这个功能,您需要以下硬件设备:FX3UPLC、FX3U-485ADP-MB通信扩展模块、施耐德ATV312变频器或台达VFD-M变频器或三菱E700变频器,以及fx3u-cnv-bd。 通过modbus rtu通信方式,您可以实现控制正反转、启动停止,还可

最新推荐

recommend-type

图像去雾算法及matlab程序实例

该算法首先将彩色图像转换为灰度图像,然后对灰度图像进行直方图均衡。下面是该算法的MATLAB程序实例: ```matlab clc; clear all; close all PS = imread('c:\\ce.png'); figure; subplot(2, 2, 1); imshow(PS); ...
recommend-type

利用MATLAB计算分形维数

I = rgb2gray(I); [G, gabout] = gaborfilter(I, 2, 4, 16, pi/16); J = fft2(gabout); A = double(J); [m, n] = size(A); B = A; C = zeros(m, n); for i = 1:m-1 for j = 1:n-1 B(i, j) = A(i+1, j+1); C(i, j)...
recommend-type

基于遗传算法的动态优化物流配送中心选址问题研究(Matlab源码+详细注释),遗传算法与免疫算法在物流配送中心选址问题的应用详解(源码+详细注释,Matlab编写,含动态优化与迭代,结果图展示),遗传

基于遗传算法的动态优化物流配送中心选址问题研究(Matlab源码+详细注释),遗传算法与免疫算法在物流配送中心选址问题的应用详解(源码+详细注释,Matlab编写,含动态优化与迭代,结果图展示),遗传算法 求解物流配送中心选址问题 源码+详细注释(Matlab编写) 有两种解决选址问题代码,说明如下: 代码一:免疫算法物流配送中心选址 模型应用场景: 1.配送中心能够配送的总量≥各揽收站需求之和 2.一个配送中心可为多个揽收站配送物,但一个快递揽收站仅由一个配送中心供应 需求点,需求点容量,配送中心数目可以根据实际随意更改(结果图如图1,2,3,4所示) 代码二:遗传算法配送中心选址 可以修改需求点坐标,需求点的需求量,备选中心坐标,配送中心个数 注:2≤备选中心≤20,需求点中心可以无限个 [new]优化与迭代过程是动态更新的喔[火]有需要的可以直接拿哈 (结果图如图5,6,7,8所示) 代码一经出不予 保证运行 可回答简单问题[托腮] ,核心关键词:遗传算法;物流配送中心选址问题;免疫算法;源码;Matlab编写;模型应用场景;需求点;配送中心;备选中心坐标;优化与迭代过程。
recommend-type

SpringBoot博客项目.zip(毕设&课设&实训&大作业&竞赛&项目)

项目工程资源经过严格测试运行并且功能上ok,可实现复现复刻,拿到资料包后可实现复现出一样的项目,本人系统开发经验充足(全栈全领域),有任何使用问题欢迎随时与我联系,我会抽时间努力为您解惑,提供帮助 【资源内容】:包含源码+工程文件+说明(如有)等。答辩评审平均分达到96分,放心下载使用!可实现复现;设计报告也可借鉴此项目;该资源内项目代码都经过测试运行;功能ok 【项目价值】:可用在相关项目设计中,皆可应用在项目、毕业设计、课程设计、期末/期中/大作业、工程实训、大创等学科竞赛比赛、初期项目立项、学习/练手等方面,可借鉴此优质项目实现复刻,设计报告也可借鉴此项目,也可基于此项目来扩展开发出更多功能 【提供帮助】:有任何使用上的问题欢迎随时与我联系,抽时间努力解答解惑,提供帮助 【附带帮助】:若还需要相关开发工具、学习资料等,我会提供帮助,提供资料,鼓励学习进步 下载后请首先打开说明文件(如有);整理时不同项目所包含资源内容不同;项目工程可实现复现复刻,如果基础还行,也可在此程序基础上进行修改,以实现其它功能。供开源学习/技术交流/学习参考,勿用于商业用途。质量优质,放心下载使用
recommend-type

基于改进蚁群算法与动态窗口法的多机器人路径规划与避障算法研究:去除冗余点、实现全局与局部实时动态规划,基于改进蚁群算法与动态窗口法的多机器人路径规划与避障算法研究:去除冗余点,实现全局与局部实时动态规

基于改进蚁群算法与动态窗口法的多机器人路径规划与避障算法研究:去除冗余点、实现全局与局部实时动态规划,基于改进蚁群算法与动态窗口法的多机器人路径规划与避障算法研究:去除冗余点,实现全局与局部实时动态规划及路径平滑,改进蚁群算法多机器人路径规划算法 改进蚁群算法+去除多余冗点(路径平滑)+全局和局部多移动机器人实时动态规划。 代码注释完全易懂,效果请看下图。 预先声明:该程序为版权所述,仅供学习参考使用,蚁群算法 动态窗口法 dwa 多机器人 路径规划算法 基于改进蚁群算法实现静态已知障碍物避障,基于动态窗口算法实现静态未知障碍物避障,动态未知障碍物避障 ,核心关键词: 改进蚁群算法; 多机器人路径规划算法; 去除多余冗点; 路径平滑; 全局和局部动态规划; 静态已知障碍物避障; 动态窗口法(DWA); 静态未知障碍物避障; 动态未知障碍物避障。,基于改进蚁群算法的机器人多路径规划与平滑技术
recommend-type

PHP集成Autoprefixer让CSS自动添加供应商前缀

标题和描述中提到的知识点主要包括:Autoprefixer、CSS预处理器、Node.js 应用程序、PHP 集成以及开源。 首先,让我们来详细解析 Autoprefixer。 Autoprefixer 是一个流行的 CSS 预处理器工具,它能够自动将 CSS3 属性添加浏览器特定的前缀。开发者在编写样式表时,不再需要手动添加如 -webkit-, -moz-, -ms- 等前缀,因为 Autoprefixer 能够根据各种浏览器的使用情况以及官方的浏览器版本兼容性数据来添加相应的前缀。这样可以大大减少开发和维护的工作量,并保证样式在不同浏览器中的一致性。 Autoprefixer 的核心功能是读取 CSS 并分析 CSS 规则,找到需要添加前缀的属性。它依赖于浏览器的兼容性数据,这一数据通常来源于 Can I Use 网站。开发者可以通过配置文件来指定哪些浏览器版本需要支持,Autoprefixer 就会自动添加这些浏览器的前缀。 接下来,我们看看 PHP 与 Node.js 应用程序的集成。 Node.js 是一个基于 Chrome V8 引擎的 JavaScript 运行时环境,它使得 JavaScript 可以在服务器端运行。Node.js 的主要特点是高性能、异步事件驱动的架构,这使得它非常适合处理高并发的网络应用,比如实时通讯应用和 Web 应用。 而 PHP 是一种广泛用于服务器端编程的脚本语言,它的优势在于简单易学,且与 HTML 集成度高,非常适合快速开发动态网站和网页应用。 在一些项目中,开发者可能会根据需求,希望把 Node.js 和 PHP 集成在一起使用。比如,可能使用 Node.js 处理某些实时或者异步任务,同时又依赖 PHP 来处理后端的业务逻辑。要实现这种集成,通常需要借助一些工具或者中间件来桥接两者之间的通信。 在这个标题中提到的 "autoprefixer-php",可能是一个 PHP 库或工具,它的作用是把 Autoprefixer 功能集成到 PHP 环境中,从而使得在使用 PHP 开发的 Node.js 应用程序时,能够利用 Autoprefixer 自动处理 CSS 前缀的功能。 关于开源,它指的是一个项目或软件的源代码是开放的,允许任何个人或组织查看、修改和分发原始代码。开源项目的好处在于社区可以一起参与项目的改进和维护,这样可以加速创新和解决问题的速度,也有助于提高软件的可靠性和安全性。开源项目通常遵循特定的开源许可证,比如 MIT 许可证、GNU 通用公共许可证等。 最后,我们看到提到的文件名称 "autoprefixer-php-master"。这个文件名表明,该压缩包可能包含一个 PHP 项目或库的主分支的源代码。"master" 通常是源代码管理系统(如 Git)中默认的主要分支名称,它代表项目的稳定版本或开发的主线。 综上所述,我们可以得知,这个 "autoprefixer-php" 工具允许开发者在 PHP 环境中使用 Node.js 的 Autoprefixer 功能,自动为 CSS 规则添加浏览器特定的前缀,从而使得开发者可以更专注于内容的编写而不必担心浏览器兼容性问题。
recommend-type

揭秘数字音频编码的奥秘:非均匀量化A律13折线的全面解析

# 摘要 数字音频编码技术是现代音频处理和传输的基础,本文首先介绍数字音频编码的基础知识,然后深入探讨非均匀量化技术,特别是A律压缩技术的原理与实现。通过A律13折线模型的理论分析和实际应用,本文阐述了其在保证音频信号质量的同时,如何有效地降低数据传输和存储需求。此外,本文还对A律13折线的优化策略和未来发展趋势进行了展望,包括误差控制、算法健壮性的提升,以及与新兴音频技术融合的可能性。 # 关键字 数字音频编码;非均匀量化;A律压缩;13折线模型;编码与解码;音频信号质量优化 参考资源链接:[模拟信号数字化:A律13折线非均匀量化解析](https://wenku.csdn.net/do
recommend-type

arduino PAJ7620U2

### Arduino PAJ7620U2 手势传感器 教程 #### 示例代码与连接方法 对于Arduino开发PAJ7620U2手势识别传感器而言,在Arduino IDE中的项目—加载库—库管理里找到Paj7620并下载安装,完成后能在示例里找到“Gesture PAJ7620”,其中含有两个示例脚本分别用于9种和15种手势检测[^1]。 关于连线部分,仅需连接四根线至Arduino UNO开发板上的对应位置即可实现基本功能。具体来说,这四条线路分别为电源正极(VCC),接地(GND),串行时钟(SCL)以及串行数据(SDA)[^1]。 以下是基于上述描述的一个简单实例程序展示如
recommend-type

网站啄木鸟:深入分析SQL注入工具的效率与限制

网站啄木鸟是一个指的是一类可以自动扫描网站漏洞的软件工具。在这个文件提供的描述中,提到了网站啄木鸟在发现注入漏洞方面的功能,特别是在SQL注入方面。SQL注入是一种常见的攻击技术,攻击者通过在Web表单输入或直接在URL中输入恶意的SQL语句,来欺骗服务器执行非法的SQL命令。其主要目的是绕过认证,获取未授权的数据库访问权限,或者操纵数据库中的数据。 在这个文件中,所描述的网站啄木鸟工具在进行SQL注入攻击时,构造的攻击载荷是十分基础的,例如 "and 1=1--" 和 "and 1>1--" 等。这说明它的攻击能力可能相对有限。"and 1=1--" 是一个典型的SQL注入载荷示例,通过在查询语句的末尾添加这个表达式,如果服务器没有对SQL注入攻击进行适当的防护,这个表达式将导致查询返回真值,从而使得原本条件为假的查询条件变为真,攻击者便可以绕过安全检查。类似地,"and 1>1--" 则会检查其后的语句是否为假,如果查询条件为假,则后面的SQL代码执行时会被忽略,从而达到注入的目的。 描述中还提到网站啄木鸟在发现漏洞后,利用查询MS-sql和Oracle的user table来获取用户表名的能力不强。这表明该工具可能无法有效地探测数据库的结构信息或敏感数据,从而对数据库进行进一步的攻击。 关于实际测试结果的描述中,列出了8个不同的URL,它们是针对几个不同的Web应用漏洞扫描工具(Sqlmap、网站啄木鸟、SqliX)进行测试的结果。这些结果表明,针对提供的URL,Sqlmap和SqliX能够发现注入漏洞,而网站啄木鸟在多数情况下无法识别漏洞,这可能意味着它在漏洞检测的准确性和深度上不如其他工具。例如,Sqlmap在针对 "http://www.2cto.com/news.php?id=92" 和 "http://www.2cto.com/article.asp?ID=102&title=Fast food marketing for children is on the rise" 的URL上均能发现SQL注入漏洞,而网站啄木鸟则没有成功。这可能意味着网站啄木鸟的检测逻辑较为简单,对复杂或隐蔽的注入漏洞识别能力不足。 从这个描述中,我们也可以了解到,在Web安全测试中,工具的多样性选择是十分重要的。不同的安全工具可能对不同的漏洞和环境有不同的探测能力,因此在实际的漏洞扫描过程中,安全测试人员需要选择合适的工具组合,以尽可能地全面地检测出应用中存在的漏洞。 在标签中指明了这是关于“sql注入”的知识,这表明了文件主题的核心所在。SQL注入是一种常见的网络攻击方式,安全测试人员、开发人员和网络管理员都需要对此有所了解,以便进行有效的防御和检测。 最后,提到了压缩包子文件的文件名称列表,其中包含了三个文件:setup.exe、MD5.exe、说明_Readme.html。这里提供的信息有限,但可以推断setup.exe可能是一个安装程序,MD5.exe可能是一个计算文件MD5散列值的工具,而说明_Readme.html通常包含的是软件的使用说明或者版本信息等。这些文件名暗示了在进行网站安全测试时,可能涉及到安装相关的软件工具,以及进行文件的校验和阅读相应的使用说明。然而,这些内容与文件主要描述的web安全漏洞检测主题不是直接相关的。
recommend-type

【GPStoolbox使用技巧大全】:20个实用技巧助你精通GPS数据处理

# 摘要 GPStoolbox是一个广泛应用于GPS数据处理的软件工具箱,它提供了从数据导入、预处理、基本分析到高级应用和自动化脚本编写的全套功能。本文介绍了GPStoolbox的基本概况、安装流程以及核心功能,探讨了如何