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 12:02:52 浏览: 134
RAR

完整版 可直接运行 MATLAB实现图像去噪 滤波 锐化 边缘检测 源程序代码 含测试图片.rar

star5星 · 资源好评率100%
好的,我进行修改后的代码如下: ```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;

最新推荐

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

tornado-6.4.1-cp38-abi3-musllinux_1_2_i686.whl

tornado-6.4.1-cp38-abi3-musllinux_1_2_i686.whl
recommend-type

tornado-6.1-cp36-cp36m-manylinux2014_aarch64.whl

tornado-6.1-cp36-cp36m-manylinux2014_aarch64.whl
recommend-type

基于java的ssm停车位短租系统程序答辩PPT.pptx

基于java的ssm停车位短租系统程序答辩PPT.pptx
recommend-type

Aspose资源包:转PDF无水印学习工具

资源摘要信息:"Aspose.Cells和Aspose.Words是两个非常强大的库,它们属于Aspose.Total产品家族的一部分,主要面向.NET和Java开发者。Aspose.Cells库允许用户轻松地操作Excel电子表格,包括创建、修改、渲染以及转换为不同的文件格式。该库支持从Excel 97-2003的.xls格式到最新***016的.xlsx格式,还可以将Excel文件转换为PDF、HTML、MHTML、TXT、CSV、ODS和多种图像格式。Aspose.Words则是一个用于处理Word文档的类库,能够创建、修改、渲染以及转换Word文档到不同的格式。它支持从较旧的.doc格式到最新.docx格式的转换,还包括将Word文档转换为PDF、HTML、XAML、TIFF等格式。 Aspose.Cells和Aspose.Words都有一个重要的特性,那就是它们提供的输出资源包中没有水印。这意味着,当开发者使用这些资源包进行文档的处理和转换时,最终生成的文档不会有任何水印,这为需要清洁输出文件的用户提供了极大的便利。这一点尤其重要,在处理敏感文档或者需要高质量输出的企业环境中,无水印的输出可以帮助保持品牌形象和文档内容的纯净性。 此外,这些资源包通常会标明仅供学习使用,切勿用作商业用途。这是为了避免违反Aspose的使用协议,因为Aspose的产品虽然是商业性的,但也提供了免费的试用版本,其中可能包含了特定的限制,如在最终输出的文档中添加水印等。因此,开发者在使用这些资源包时应确保遵守相关条款和条件,以免产生法律责任问题。 在实际开发中,开发者可以通过NuGet包管理器安装Aspose.Cells和Aspose.Words,也可以通过Maven在Java项目中进行安装。安装后,开发者可以利用这些库提供的API,根据自己的需求编写代码来实现各种文档处理功能。 对于Aspose.Cells,开发者可以使用它来完成诸如创建电子表格、计算公式、处理图表、设置样式、插入图片、合并单元格以及保护工作表等操作。它也支持读取和写入XML文件,这为处理Excel文件提供了更大的灵活性和兼容性。 而对于Aspose.Words,开发者可以利用它来执行文档格式转换、读写文档元数据、处理文档中的文本、格式化文本样式、操作节、页眉、页脚、页码、表格以及嵌入字体等操作。Aspose.Words还能够灵活地处理文档中的目录和书签,这让它在生成复杂文档结构时显得特别有用。 在使用这些库时,一个常见的场景是在企业应用中,需要将报告或者数据导出为PDF格式,以便于打印或者分发。这时,使用Aspose.Cells和Aspose.Words就可以实现从Excel或Word格式到PDF格式的转换,并且确保输出的文件中不包含水印,这提高了文档的专业性和可信度。 需要注意的是,虽然Aspose的产品提供了很多便利的功能,但它们通常是付费的。用户需要根据自己的需求购买相应的许可证。对于个人用户和开源项目,Aspose有时会提供免费的许可证。而对于商业用途,用户则需要购买商业许可证才能合法使用这些库的所有功能。"
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【R语言高性能计算秘诀】:代码优化,提升分析效率的专家级方法

![R语言](https://www.lecepe.fr/upload/fiches-formations/visuel-formation-246.jpg) # 1. R语言简介与计算性能概述 R语言作为一种统计编程语言,因其强大的数据处理能力、丰富的统计分析功能以及灵活的图形表示法而受到广泛欢迎。它的设计初衷是为统计分析提供一套完整的工具集,同时其开源的特性让全球的程序员和数据科学家贡献了大量实用的扩展包。由于R语言的向量化操作以及对数据框(data frames)的高效处理,使其在处理大规模数据集时表现出色。 计算性能方面,R语言在单线程环境中表现良好,但与其他语言相比,它的性能在多
recommend-type

在构建视频会议系统时,如何通过H.323协议实现音视频流的高效传输,并确保通信的稳定性?

要通过H.323协议实现音视频流的高效传输并确保通信稳定,首先需要深入了解H.323协议的系统结构及其组成部分。H.323协议包括音视频编码标准、信令控制协议H.225和会话控制协议H.245,以及数据传输协议RTP等。其中,H.245协议负责控制通道的建立和管理,而RTP用于音视频数据的传输。 参考资源链接:[H.323协议详解:从系统结构到通信流程](https://wenku.csdn.net/doc/2jtq7zt3i3?spm=1055.2569.3001.10343) 在构建视频会议系统时,需要合理配置网守(Gatekeeper)来提供地址解析和准入控制,保证通信安全和地址管理
recommend-type

Go语言控制台输入输出操作教程

资源摘要信息:"在Go语言(又称Golang)中,控制台的输入输出是进行基础交互的重要组成部分。Go语言提供了一组丰富的库函数,特别是`fmt`包,来处理控制台的输入输出操作。`fmt`包中的函数能够实现格式化的输入和输出,使得程序员可以轻松地在控制台显示文本信息或者读取用户的输入。" 1. fmt包的使用 Go语言标准库中的`fmt`包提供了许多打印和解析数据的函数。这些函数可以让我们在控制台上输出信息,或者从控制台读取用户的输入。 - 输出信息到控制台 - Print、Println和Printf是基本的输出函数。Print和Println函数可以输出任意类型的数据,而Printf可以进行格式化输出。 - Sprintf函数可以将格式化的字符串保存到变量中,而不是直接输出。 - Fprint系列函数可以将输出写入到`io.Writer`接口类型的变量中,例如文件。 - 从控制台读取信息 - Scan、Scanln和Scanf函数可以读取用户输入的数据。 - Sscan、Sscanln和Sscanf函数则可以从字符串中读取数据。 - Fscan系列函数与上面相对应,但它们是将输入读取到实现了`io.Reader`接口的变量中。 2. 输入输出的格式化 Go语言的格式化输入输出功能非常强大,它提供了类似于C语言的`printf`和`scanf`的格式化字符串。 - Print函数使用格式化占位符 - `%v`表示使用默认格式输出值。 - `%+v`会包含结构体的字段名。 - `%#v`会输出Go语法表示的值。 - `%T`会输出值的数据类型。 - `%t`用于布尔类型。 - `%d`用于十进制整数。 - `%b`用于二进制整数。 - `%c`用于字符(rune)。 - `%x`用于十六进制整数。 - `%f`用于浮点数。 - `%s`用于字符串。 - `%q`用于带双引号的字符串。 - `%%`用于百分号本身。 3. 示例代码分析 在文件main.go中,可能会包含如下代码段,用于演示如何在Go语言中使用fmt包进行基本的输入输出操作。 ```go package main import "fmt" func main() { var name string fmt.Print("请输入您的名字: ") fmt.Scanln(&name) // 读取一行输入并存储到name变量中 fmt.Printf("你好, %s!\n", name) // 使用格式化字符串输出信息 } ``` 以上代码首先通过`fmt.Print`函数提示用户输入名字,并等待用户从控制台输入信息。然后`fmt.Scanln`函数读取用户输入的一行信息(包括空格),并将其存储在变量`name`中。最后,`fmt.Printf`函数使用格式化字符串输出用户的名字。 4. 代码注释和文档编写 在README.txt文件中,开发者可能会提供关于如何使用main.go代码的说明,这可能包括代码的功能描述、运行方法、依赖关系以及如何处理常见的输入输出场景。这有助于其他开发者理解代码的用途和操作方式。 总之,Go语言为控制台输入输出提供了强大的标准库支持,使得开发者能够方便地处理各种输入输出需求。通过灵活运用fmt包中的各种函数,可以轻松实现程序与用户的交互功能。