灰度图像 伪彩色处理 matlab

function P=psedocolor（I）
G=im2double（I）;
I=im2uint8（I）;
I=fix（I/4）;
figureimshow（I[]）
idxr1=find（I<=23）;
idxr2=find（（I>=24&I<=38））;
idxr3=find（（I>=39&I<=55））;
idxr4=find（（I>=56&I<=63））;
red=zeros（size（G））;
red=G;
red（idxr2）=（1/16）*I（idxr2）-23/16;
red（idxr3）=1;
red（idxr4）=-（1/16）*I（idxr4）+71/16;
red=mat2gray（red）;
figureimshow（red[]）
idxg1=find（（I>=12&I<=27））;
idxg2=find（（I>=28&I<=44））;
idxg3=find（（I>=45&I<=59））;
idxg4=find（（I<=12）|（I>=60&I<=63））;
green=zeros（size（G）