单波段伪彩色密度分割matlab编程.docxVIP

clc;clear all;close all; %单波段伪彩色密度分割 image=imread(全色图像 1.jpg);%读入图片 img1=rgb2gray(image); figure; imshow(img1); % [a,b]=size(img1); imgrgb=imread(真彩色图像 1.jpg); figure; imshow(imgrgb); % figure % imhist(img1); ind1=find(img135); ind2=find(img1=35img141); ind3=find(img1=41img152); ind4=find(img1=52img163); ind5=find(img1=63img192); ind6=find(img1=92img196); ind7=find(img1=96img1115); ind8=find(img1=115img1140); ind9=find(img1140img1170); ind10=find(img1=170); resR=img1; resG=img1; resB=img1; resR(ind1)=22; resG(ind1)=47; resB(ind1)=18; resR(ind2)=17; resG(ind2)=24; resB(ind2)=53; resR(ind3)=16; resG(ind3)=25; resB(ind3)=64; resR(ind4)=19; resG(ind4)=28; resB(ind4)=83; resR(ind5)=27; resG(ind5)=45; resB(ind5)=125; resR(ind6)=81; resG(ind6)=123; resB(ind6)=60; resR(ind7)=101; resG(ind7)=146; resB(ind7)=79; resR(ind8)=113; resG(ind8)=153; resB(ind8)=100; resR(ind9)=115; resG(ind9)=156; resB(ind9)=142; resR(ind10)=213; resG(ind10)=222; resB(ind10)=159; image2(:,:,1)=resR; image2(:,:,2)=resG; image2(:,:,3)=resB; figure; imshow(image2);