PointOperation04_Matlab灰度gray图intensity windowing

最新推荐文章于 2023-03-01 17:36:55 发布

DBzs

最新推荐文章于 2023-03-01 17:36:55 发布

阅读量372

点赞数

分类专栏：图像处理 matlab编程文章标签： matlab

本文链接：https://blog.csdn.net/weixin_43290709/article/details/106162752

版权

图像处理同时被 2 个专栏收录

22 篇文章 2 订阅

订阅专栏

matlab编程

17 篇文章 0 订阅

订阅专栏

1 计算原理

一种 clamp 操作，然后将一定范围内的intensity线性拉伸到某个指定区间。
记某像素点的值为p，对各像素点做以下操作：
$\begin{cases} a& {p < a}\\ M\dfrac{p - a}{b - a} & {a \leq p \leq b}\\ M& {p > b} \end{cases}$
则将图像window在 $[a, b]$ 内，最大intensity为 $M$ .

2 代码实现

%% reading in the image
orgIMG = rgb2gray(imread("lena.jpg"));
figure(1);
subplot(1,2,1);
imshow(orgIMG);
title("Original IMG");
subplot(1,2,2);
imhist(orgIMG);
axis([0 255 0 800]);
title("The Histogram of the original img");

%% Doing intensity windowing
% Define parameters
a_min = 20;    % a_min <- 0~a_max
a_max = 220;    % a_max <- a_min~255
M = 220;        % intensity


orgIMG = double(orgIMG);

smaller = find(orgIMG <= a_min);
larger = find(orgIMG >= a_max);
middle = find(orgIMG > a_min & orgIMG < a_max);

orgIMG(smaller) = a_min;
orgIMG(larger) = M;
orgIMG(middle) = M.*(orgIMG(middle)-a_min)./(a_max-a_min);

orgIMG = uint8(orgIMG);

%% Draw the new img
figure(2);
subplot(1,2,1);
imshow(orgIMG);
title("Modified IMG");
subplot(1,2,2);
imhist(orgIMG);
axis([0 255 0 800]);
title("The Histogram of the Modified img");