一杯清酒邀明月
天下本无事,庸人扰之而烦耳。
posts - 3121,comments - 209,views - 578万

色相饱和度明度和三原色的相互转换

复制代码
 1 clear all;
 2 close all;
 3 clc;
 4 
 5 img=imread('lena_color.jpg');
 6 img=mat2gray(img);  %任意区间映射到[0,1];
 7 [m n dim]=size(img);
 8 imshow(img);
 9 %%图像的RGB
10 R=img(:,:,1);
11 G=img(:,:,2);
12 B=img(:,:,3);
13 
14 %%图像的RGB2HSV
15 H=zeros(m,n);   %色相角
16 S=zeros(m,n);   %饱和度
17 V=zeros(m,n);   %明度
18 for i=1:m
19    for j=1:n
20        r=R(i,j);
21        g=G(i,j);
22        b=B(i,j);
23        MAX=max([r,g,b]);
24        MIN=min([r,g,b]);
25 
26        if MAX==MIN
27             H(i,j)=0;
28        elseif MAX==r && g>=b
29             H(i,j)=60*(g-b)/(MAX-MIN);
30        elseif MAX==r && g<b
31             H(i,j)=60*(g-b)/(MAX-MIN)+360;
32        elseif MAX==g
33             H(i,j)=60*(b-r)/(MAX-MIN)+120;
34        elseif MAX==b
35             H(i,j)=60*(r-g)/(MAX-MIN)+240;
36        end
37 
38        if MAX==0
39             S(i,j)=0;
40        else
41             S(i,j)=1-MIN/MAX;
42        end
43 
44        V(i,j)=MAX;
45    end
46 end
47 
48 %%图像HSV2RGB
49 for i=1:m
50     for j=1:n
51         h=floor(H(i,j)/60);
52         f=H(i,j)/60-h;
53         v=V(i,j);
54         s=S(i,j);
55         p=v*(1-s);
56         q=v*(1-f*s);
57         t=v*(1-(1-f)*s);
58 
59         if h==0
60             R(i,j)=v;G(i,j)=t;B(i,j)=p;
61         elseif h==1
62             R(i,j)=q;G(i,j)=v;B(i,j)=p;            
63         elseif h==2
64             R(i,j)=p;G(i,j)=v;B(i,j)=t;            
65         elseif h==3
66             R(i,j)=p;G(i,j)=q;B(i,j)=v;            
67         elseif h==4
68             R(i,j)=t;G(i,j)=p;B(i,j)=v;            
69         elseif h==5
70             R(i,j)=v;G(i,j)=p;B(i,j)=q;            
71         end
72     end
73 end
74 
75 %%如果正反变换都没错的话,那么图像是不变的
76 img(:,:,1)=R;
77 img(:,:,2)=G;
78 img(:,:,3)=B;
79 figure;
80 imshow(img)
复制代码

 

posted on   一杯清酒邀明月  阅读(647)  评论(0编辑  收藏  举报
编辑推荐:
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
· 没有源码,如何修改代码逻辑?
· 一个奇形怪状的面试题:Bean中的CHM要不要加volatile?
· [.NET]调用本地 Deepseek 模型
· 一个费力不讨好的项目,让我损失了近一半的绩效!
阅读排行:
· 全网最简单!3分钟用满血DeepSeek R1开发一款AI智能客服,零代码轻松接入微信、公众号、小程
· .NET 10 首个预览版发布,跨平台开发与性能全面提升
· 《HelloGitHub》第 107 期
· 全程使用 AI 从 0 到 1 写了个小工具
· 从文本到图像:SSE 如何助力 AI 内容实时呈现?(Typescript篇)
< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5

点击右上角即可分享
微信分享提示