It's so simple, I am definitely sure you can do it yourself, which by the way is good to practice your programming skills.

The pseudo code reads something like that, assuming the array to encode is named in_buf (not actual matlab syntax):

cur_val = in_buf[0]

i = 1, k = 1, o = 0

while i < length(in_buf) {

  if in_buf[i] == cur_val {

    k += 1

  }

  else {

    out_buf[o] = k

    out_buf[o+1] = curval

    o += 2

    cur_val = in_buf[i]

    k = 1

  }

  i += 1

}

Sure you can do more fancy stuff using find and cumsum, which is left as exercise :-)

Thank u but what about decoding

Sorry can u check this code

function [im] = rlemain (infile,outfile)

% Example of use:

% out=rlemain('input_image.dcm','output_image.dcm');

I = dicomread(infile);

I = mat2gray(I);

I=im2uint8(I);

figure('Name','Input image');

imshow(I);

level=graythresh(I);

bw=I;

a=bw'; a=a(:); a=a';

a=double(a);

rle(1)=a(1); m=2; rle(m)=1;

for i=1:length(a)-1

if a(i)==a(i+1)

rle(m)=rle(m)+1;

else

m=m+1; rle(m)=1; %dynamic allocation and initialization of next element of rle

end

end

figure;imshow(bw);title('RLE image');

dicomwrite(bw,'Rle image.dcm');