【数字图像处理matlab】RGB转HIS,HIS转RGB

【数字图像处理matlab】RGB转HIS,HIS转RGB

读入一张RGB图像,将其转换为HIS,再将HIS转换为RGB。

function F=color_processing(image)
%调用代码-----------------------------------------
%imageR=imread('flower.jpg');  
%image=im2double(imageR);color_processing(image);
%------------------------------------------------

%RGB->HIS变换,rgb2hsv()
%从三维数组中提取三幅分量图像
R1=image(:,:,1);
G1=image(:,:,2);
B1=image(:,:,3);
I=(R1+G1+B1)/3;        %亮度分量,范围[0,1]
m=min(min(R1,G1),B1);
S=1-3*m./(R1+G1+B1);    %饱和度分量,范围[0,1]
theta=acos(((R1-G1)+(R1-B1))./(2*((R1-G1).^2+((R1-B1).*(G1-B1))).^(1/2))); %弧度 
H=theta;   %色度分量,以角度表示,范围是[0,1](弧度除以2*pi后)
if B1>G1
    H=2*pi-theta;       
end     
if S==0
    H=0;
end
H=H/(2*pi);
Image1=cat(3,H,S,I);
subplot(1,3,1),imshow(image);
title('RGB原图')
subplot(1,3,2),imshow(Image1);
title('转换后的HSI图像')

%HIS->RGB变换,Image1转为imdouble型

H1=Image1(:,:,1)*2*pi;  %弧度
S1=Image1(:,:,2);
I1=Image1(:,:,3);

R1 = zeros(size(Image1, 1), size(Image1, 2)); 
G1= zeros(size(Image1, 1), size(Image1, 2)); 
B1 = zeros(size(Image1, 1), size(Image1, 2));

i= find( (0 <= H1) & (H1 < 2*pi/3)); 
B1(i) = I1(i) .* (1 - S1(i)); 
R1(i) = I1(i) .* (1 + S1(i) .* cos(H1(i)) ./ cos(pi/3 - H1(i))); 
G1(i) = 3*I1(i) - (R1(i) + B1(i)); 

i = find( (2*pi/3 <= H1) & (H1 < 4*pi/3) ); 
R1(i) = I1(i) .* (1 - S1(i)); 
G1(i) = I1(i) .* (1 + S1(i) .* cos(H1(i) - 2*pi/3) ./ cos(pi - H1(i))); 
B1(i) = 3*I1(i) - (R1(i) + G1(i)); 

i = find( (4*pi/3 <= H1) & (H1 <= 2*pi)); 
G1(i) = I1(i) .* (1 - S1(i)); 
B1(i) = I1(i) .* (1 + S1(i) .* cos(H1(i) - 4*pi/3) ./cos(5*pi/3 - H1(i))); 
R1(i) = 3*I1(i) - (G1(i) + B1(i)); 

Image2 = cat(3, R1, G1, B1); 
Image2 = max(min(Image2, 1), 0); 
subplot(1,3,3),imshow(Image2);
title('转换回的RGB图像') 


end


测试:
在这里插入图片描述

好的,针对你的问题,我可以给出以下Python算法实现: 1. RGBHIS ```python import numpy as np def rgb_to_his(rgb): # 将RGB换为[0,1]范围内的浮点数 r, g, b = rgb[0]/255.0, rgb[1]/255.0, rgb[2]/255.0 cmax, cmin = max(r, g, b), min(r, g, b) # 计算亮度值(I) intensity = (cmax + cmin) / 2.0 if cmax == cmin: return np.array([0.0, 0.0, intensity]) # 计算饱和度(S) chroma = cmax - cmin if intensity < 0.5: saturation = chroma / (cmax + cmin) else: saturation = chroma / (2.0 - cmax - cmin) # 计算色相(H) if r == cmax: hue = (g - b) / chroma elif g == cmax: hue = 2.0 + (b - r) / chroma else: hue = 4.0 + (r - g) / chroma hue *= 60.0 if hue < 0: hue += 360.0 return np.array([hue, saturation, intensity]) ``` 2. HISRGB ```python import numpy as np def his_to_rgb(his): hue, saturation, intensity = his[0], his[1], his[2] if saturation == 0: r, g, b = intensity, intensity, intensity else: if hue == 360: hue = 0 hue /= 60.0 i = int(hue) f = hue - i p = intensity * (1.0 - saturation) q = intensity * (1.0 - saturation * f) t = intensity * (1.0 - saturation * (1.0 - f)) if i == 0: r, g, b = intensity, t, p elif i == 1: r, g, b = q, intensity, p elif i == 2: r, g, b = p, intensity, t elif i == 3: r, g, b = p, q, intensity elif i == 4: r, g, b = t, p, intensity else: r, g, b = intensity, p, q # 将RGB换为[0,255]范围内的整数 return np.array([int(r*255.0), int(g*255.0), int(b*255.0)]) ``` 以上就是RGB和HIS之间的换算法实现,希望能够对你有所帮助。
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值