[X,map] = imread('trees.tif');
I = ind2gray(X,map);
imshow(I)
I2 = roifill;
% select the region to fill
figure, imshow(I2)