I=imread（‘G:\tupian\renwu1.jpg‘）;%讀入圖像數據，此時I為三維數組
I=rgb2gray（I）;%轉化為灰度圖像，I變為二維數組
I=double（I）;%轉化成double型
[mn]=size（I）;%獲得I的大小
c=256;%定義常量

R?=?ones（mn）;%預定義RGB三個維度的數組
G?=?ones（mn）;
B?=?ones（mn）;

for?i=1:m%根據灰度等級，計算偽彩色相應的數據
????for?j=1:n
????????if?I（ij）<=c/4
????????????R（ij）=0;
????????????G（ij）=4*I（ij）;
????????????B（ij）=c;
????????elseif?I（ij）<=c/2
????????????R（ij）=0;
????????????G（ij）=c;
????????????B（ij）=-4*I