图像 行程编码与解码的matlab实现

%注意，如果已经运行了一遍该程序，那么此时RLEcode矩阵的行数就确定了，再次运行该程序时，如果选择了一个编码后产生
%的矩阵没有原先生成的RLEcode矩阵大的图像，那么此时编码生成的RLEcode矩阵大小是先前比较大的RLEcode矩阵的大小，
%对于多出来的那些行，其值保留为原先RLEcode矩阵的值，所以此时再解码会出现错误。
%所以最好是每次要运行该程序时都先关闭该程序或清空一下工作区的变量，然后再启动程序。

% open a picture
[filename, filepath] = uigetfile({'*.jpg;*.ppm; jpeg *.;*.bmp;*.png'},'Choose Input Image');
if isequal(filename,0) || isequal(filepath,0)
    disp('User pressed cancel')
    return
else
    fullfp = fullfile(filepath, filename);
end
I = imread(fullfp);   %代表要处理的图像
IGRAY=rgb2gray(I);

%将灰度图转换成二值图
% tt = graythresh(IGRAY);
% IGRAY = imbinarize(IGRAY,tt);
subplot(1,2,1); imshow(IGRAY,[]); title('编码前图像');
[m, n]=size(IGRAY);
c=IGRAY(1,1);
RLEcode(1,1:3)=[1 1 c];%建立数组RLEcode，其中元素排列形式为[行程起始行坐标、行程列坐标、灰度值]
t=2;
for k=1:m %进行行程编码，矩阵RLEcode存储编码后的图像信息
    for j=1:n
        if(not(and(k=&#