甘肃龙网

  1. 首页 > 生活 > >

【三】图像的点运算( 二 )

生活百科直方图


5 灰度阈值变换
简单来说就是设置一个阈值，把灰度图变成二值图，小于阈值的就置0，大于阈值的置1 。
实现
I = imread('rice.png');thresh = graythresh(I);bw1 = im2bw(I,thresh);bw2 = im2bw(I,130/255);subplot(1,3,1);imshow(I);title('原图像');subplot(1,3,2);imshow(bw1);title('自动选择阈值');subplot(1,3,3);imshow(bw2);title('阈值130');

【三】图像的点运算

文章插图
6 分段线性变换
就是分段去增强或抑制灰度区域 。
实现
首先自定义函数.m
function out = imgrayscaling(varargin)% IMGRAYSCALING执行灰度拉伸功能%语法：%out = imgrayscaling(I, [x1,x2], [y1,y2]);%out = imgrayscaling(X, map, [x1,x2], [y1,y2]);%out = imgrayscaling(RGB, [x1,x2], [y1,y2]);%这个函数提供灰度拉伸功能，输入图像应当是灰度图像，但如果提供的不是灰度%图像的话，函数会自动将图像转化为灰度形式 。x1，x2，y1，y2应当使用双精度%类型存储，图像矩阵可以使用任何MATLAB支持的类型存储 。[A, map, x1 , x2, y1, y2] = parse_inputs(varargin{:});% 计算输入图像A中数据类型对应的取值范围range = getrangefromclass(A);range = range(2);% 如果输入图像不是灰度图，则需要执行转换if ndims(A)==3,% A矩阵为3维，RGB图像A = rgb2gray(A);elseif ~isempty(map),% MAP变量为非空，索引图像A = ind2gray(A,map);end % 对灰度图像则不需要转换% 读取原始图像的大小并初始化输出图像[M,N] = size(A);I = im2double(A);% 将输入图像转换为双精度类型out = zeros(M,N);% 主体部分，双级嵌套循环和选择结构for i=1:Mfor j=1:Nif I(i,j)x2out(i,j) = (I(i,j)-x2)*(range-y2)/(range-x2) + y2;elseout(i,j) = (I(i,j)-x1)*(y2-y1)/(x2-x1) + y1;endendend% 将输出图像的格式转化为与输入图像相同if isa(A, 'uint8') % uint8out = im2uint8(out);elseif isa(A, 'uint16')out = im2uint16(out);% 其它情况，输出双精度类型的图像end% 输出:if nargout==0 % 如果没有提供参数接受返回值imshow(out);return;end%-----------------------------------------------------------------------------function [A, map, x1, x2, y1, y2] = parse_inputs(varargin);% 这就是用来分析输入参数个数和有效性的函数parse_inputs% A输入图像，RGB图 (3D), 灰度图 (2D), 或者索引图 (X)% map索引图调色板 (:,3)% [x1,x2] 参数组 1，曲线中两个转折点的横坐标% [y1,y2] 参数组 2，曲线中两个转折点的纵坐标% 首先建立一个空的map变量，以免后面调用isempty(map)时出错map = [];%IPTCHECKNARGIN(LOW,HIGH,NUM_INPUTS,FUNC_NAME) 检查输入参数的个数是否%符合要求，即NUM_INPUTS中包含的输入变量个数是否在LOW和HIGH所指定的范围%内 。如果不在范围内，则此函数给出一个格式化的错误信息 。iptchecknargin(3,4,nargin,mfilename);%IPTCHECKINPUT(A,CLASSES,ATTRIBUTES,FUNC_NAME,VAR_NAME, ARG_POS) 检查给定%矩阵A中的元素是否属于给定的类型列表 。如果存在元素不属于给定的类型，则给出%一个格式化的错误信息 。iptcheckinput(varargin{1},...{'uint8','uint16','int16','double'}, ...{'real', 'nonsparse'},mfilename,'I, X or RGB',1);switch nargincase 3 %可能是imgrayscaling(I, [x1,x2], [y1,y2]) 或 imgrayscaling(RGB, [x1,x2], [y1,y2])A = varargin{1};x1 = varargin{2}(1);x2 = varargin{2}(2);y1 = varargin{3}(1);y2 = varargin{3}(2);case 4A = varargin{1};%imgrayscaling(X, map, [x1,x2], [y1,y2])map = varargin{2};x1 = varargin{2}(1);x2 = varargin{2}(2);y1 = varargin{3}(1);y2 = varargin{3}(2);end% 检测输入参数的有效性% 检查RGB数组if (ndims(A)==3) && (size(A,3)~=3)msg = sprintf('%s: 真彩色图像应当使用一个M-N-3维度的数组', ...upper(mfilename));eid = sprintf('Images:%s:trueColorRgbImageMustBeMbyNby3',mfilename);error(eid,'%s',msg);endif ~isempty(map) % 检查调色板if (size(map,2) ~= 3) || ndims(map)>2msg1 = sprintf('%s: 输入的调色板应当是一个矩阵', ...upper(mfilename));msg2 = '并拥有三列';eid = sprintf('Images:%s:inColormapMustBe2Dwith3Cols',mfilename);error(eid,'%s %s',msg1,msg2);elseif (min(map(:))


上一篇：基础概念：书中的域、final关键字

下一篇：ETCD原理详细解析