《数字图像处理原理与实践(MATLAB版)》一书之代码Part2[通俗易懂]

《数字图像处理原理与实践(MATLAB版)》一书之代码Part2

大家好,又见面了,我是全栈君。

本文系《数字图像处理原理与实践(MATLAB版)》一书之代码系列的Part2(P43~80),代码运行结果请參见原书配图,建议下载代码前阅读下文:

关于《数字图像处理原理与实践(MATLAB版)》一书代码公布的说明

 

http://blog.csdn.net/baimafujinji/article/details/40987807

P44

i = imread(‘theatre.jpg’);
i = rgb2gray(i);
i = double(i);

out1 = log(1+i)/0.065;
out2 = log(1+i)/0.035;
out1(find(out1>255)) = 255;
out2(find(out2>255)) = 255;
out1 = uint8(out1);
out2 = uint8(out2);

subplot(221), imshow(out1), title(‘image, p = 0.065’);
subplot(222), imhist(out1), title(‘histgram, p = 0.065’);
subplot(223), imshow(out2), title(‘image, p = 0.035’);
subplot(224), imhist(out2), title(‘histgram, p = 0.035’);

P47

i = rgb2gray(imread(‘theatre.jpg’));
i = double(i);
y1 = 255*(i/255).^2.5;
y2 = 255*(i/255).^0.4;
y1 = uint8(y1);
y2 = uint8(y2);
subplot(221), imshow(y1), title(‘p =2.5’);
subplot(222), imhist(y1), title(‘p =2.5’);
subplot(223), imshow(y2), title(‘p =0.4’);
subplot(224), imhist(y2), title(‘p =0.4’);

P48

i = rgb2gray(imread(‘theatre.jpg’));
y1 = double(i);
y1 = 255*(y1/255).^2.5;
y2 = uint8(y1);
y3 = imadjust(y2, [ ], [ ], 0.4);
subplot(131), imshow(i), title(‘original image’);
subplot(132), imshow(y2),title(‘power = 2.5’);
subplot(133), imshow(y3),title(‘gamma = 0.4’);

P49

i = imread(‘theatre.jpg’);
i = rgb2gray(i);
i = double(i);

y1 = 1.5.^(i*0.070)-1;
y2 = 1.5.^(i*0.050)-1;
y1(find(y1>255)) = 255;
y2(find(y2>255)) = 255;
y1 = uint8(y1);
y2 = uint8(y2);

subplot(221), imshow(y1), title(‘c=0.070’);
subplot(222), imhist(y1), title(‘c=0.070’);
subplot(223), imshow(y2), title(‘c=0.050’);
subplot(224), imhist(y2), title(‘c=0.050’);

P52

i = imread(‘theatre.jpg’);
i = rgb2gray(i);
L = imadjust(i,[ ],[50/255;150/255]);
J = imadjust(L,[50/255;150/255 ],[20/255;230/255]);
subplot(221), imshow(L), title(‘low contrast’);
subplot(222), imhist(L), title(‘low contrast’);
subplot(223), imshow(J), title(‘gray stretch’);
subplot(224), imhist(J), title(‘gray stretch’);

P54

i = rgb2gray(imread(‘theatre.jpg’));
LC = imadjust(i,[ ],[50/255;150/255]);
figure(1), subplot(221), imshow(LC);
title(‘low contrast’);
figure(1),subplot(222), imhist(LC);
title(‘low contrast’);
HE1 = histeq(LC);
figure(1), subplot(223), imshow(HE1);
title(‘histogram equalization’);
figure(1),subplot(224), imhist(HE1);
title(‘histogram equalization’);

P56

img = rgb2gray(imread(‘theatre.jpg’));
img_ref = rgb2gray(imread(‘rpic.jpg’));
[hgram, x] = imhist(img_ref);
J = histeq(img, hgram);
subplot(2,3,1), imshow(img), title(‘original image’);
subplot(2,3,4), imhist(img), title(‘original image’);
subplot(2,3,2), imshow(img_ref), title(‘reference image’);
subplot(2,3,5), imhist(img_ref), title(‘reference image’);
subplot(2,3,3), imshow(J), title(‘output image’);
subplot(2,3,6), imhist(J), title(‘output image’);

P64-1

I = imread(‘apostles.jpg’);
I = double(I);
B = zeros(size(I));
H = size(I);
move_x = 100;
move_y = 150;
B(move_y + 1:H(1), move_x+1:H(2), 1:H(3))=…
I(1:H(1)-move_y, 1:H(2) – move_x, 1:H(3));
subplot(1,2,1),subimage(uint8(I))
title(‘原图像’)
subplot(1,2,2),subimage(uint8(B))
title(‘平移变换’);

P64-2

I = imread(‘apostles.jpg’);
se=translate(strel(1),[150 100]);
B = imdilate(I,se);
figure;
subplot(1,2,1),subimage(I);
title(‘原图像’);
subplot(1,2,2),subimage(B);
title(‘平移变换’);

P66

I = imread(‘apostles.jpg’);
[height, width, dim]=size(I);
%水平镜像变换
tform = maketform(‘affine’,[-1 0 0;0 1 0; width 0 1]);
B=imtransform(I, tform, ‘nearest’);
%垂直镜像变换
tform2 = maketform(‘affine’, [1 0 0; 0 -1 0; 0 height 1]);
C=imtransform(I, tform2, ‘nearest’);
subplot(1,3,1),imshow(I);
title(‘原图像’);
subplot(1,3,2),imshow(B);
title(‘水平图像’);
subplot(1,3,3),imshow(C);
title(‘垂直图像’);

P67

A = imread(‘apostles.jpg’);
A = double(A);
figure(1), imshow(uint8(A));
H = size(A);
figure(2),B(1:H(1),1:H(2),1:H(3))=A(H(1):-1:1,1:H(2),1:H(3));%垂直镜像
imshow(uint8(B));
figure(3),C(1:H(1),1:H(2),1:H(3))=A(1:H(1),H(2):-1:1,1:H(3));%水平镜像
imshow(uint8(C));

P69

I = imread(‘apostles.jpg’);
tform = maketform(‘affine’,[0 1 0; 1 0 0; 0 0 1]);%定义转置矩阵
B = imtransform(I, tform, ‘nearest’);
subplot(1,2,1),imshow(I)
title(‘原图像’);
subplot(1,2,2),imshow(B)
title(‘转置图像’);

P74

I = imread(‘C:\apostles.jpg’);
A = imresize(I, 1.5, ‘nearest’);
B = imresize(I, 1.5, ‘bilinear’);
C = imresize(I, 1.5, ‘bicubic’);
subplot(2,2,1), imshow(I), title(‘original’);
subplot(2,2,2), imshow(A), title(‘nearest’);
subplot(2,2,3), imshow(B), title(‘bilinear’);
subplot(2,2,4), imshow(C), title(‘bicubic’);

P80

I = imread(‘apostles.jpg’);
A = imrotate(I, 30, ‘nearest’);%旋转30度,最邻近插值
figure(1),imshow(A)
B = imrotate(I, 45, ‘bilinear’,’loose’);%旋转45度,二次线性插值
figure(2),imshow(B)


(代码公布未完,请待兴许…)

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/116517.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 【AC大牛陈鸿的ACM总结贴】【ID AekdyCoin】人家当初也一样是菜鸟

    【AC大牛陈鸿的ACM总结贴】【ID AekdyCoin】人家当初也一样是菜鸟 acm总结帖_ByAekdyCoin     各路大牛都在中国大陆的5个赛区结束以后纷纷发出了退役帖,总结帖,或功德圆满,或死不瞑目,而这也许又会造就明年的各种“炸尸”风波。为了考虑在发退役贴以后明年我也成为“僵尸”的可能性,于是改名曰“总结贴”,不提比赛细节,不提比赛流水账,权当是大学本科生涯中acm生活的点滴记录……   (1)入门篇甲…

  • 如何查看Linux版本信息?

    如何查看Linux版本信息?这里所谓的Linux版本信息,包括Linux内核版本信息和Linux系统版本信息。下面分别介绍>>>(笔者在Ubuntu12.04中进行操作并截图如下)一、查看Linux内核版本信息方法1:登陆Linux,在终端输入cat/proc/version方法2:登陆Linux,在终端输入uname-a方法3:登陆Linux,在终端输入uname

  • Redis 6379端口无法telnet登陆

    Redis 6379端口无法telnet登陆一、现象:Redis所在服务器无法telnet 二、问题原因  redis-server默认开放的127.0.0.1IP地址需要修改为0.0.0.0才可以被其他机器访问三、解决办法(Windows)  1,在Redis安装目录下找到配置文件 redis.windows.config  2,修改bind127.0.0.1改为bind0.0.0.0    3,重启Redis…

  • linux常用的20个命令面试_linux常用命令详解pdf

    linux常用的20个命令面试_linux常用命令详解pdf什么是linux多用户,多任务,支持多线程和多CPU的操作系统linux的应用领域:免费,稳定,高效的,一般运行在大型服务器上用xshell连接虚拟机的步骤:1.setup设置虚拟机IP为10.10.10.10重启网卡:servicenetworerestart2.在Windows系统打开网络和共享中心,更改适配器设置,把vmnet1的ipv4的IP改成10.10.10.13.打开xshell,输入ssh10.10.10.10/根目录:一般根目录下只存放目录,有且只有一个根目

  • H3C交换机常用配置命令大全

    H3C交换机常用配置命令大全H3C交换机常用配置命令大全一.用户配置:<H3C>system-view[H3C]superpasswordH3C设置用户分级密码[H3C]undosuperpassword删除用户分级密码[H3C]localuserbigheap1234561Web网管用户设置,1(缺省)为管理级用户…

  • nginx安装及部署

    nginx安装及部署下载官方网站:https://nginx.org/en/download.htmlWindows下安装安装下载后解压(切记不能含有中文路径!!),文件结构如图(我解压的路径就有中文,记得拷贝放置于英文目录下即可!):启动两种方法:1)直接双击该目录下的”nginx.exe”,即可启动nginx服务器;2)命令行进入该文件夹,执行startngin…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号