
X = imread('img14.tif');


R = double( X(:,:,1) );
G = double( X(:,:,2) );
B = double( X(:,:,3) );


cmap = [0:255; 0:255; 0:255]'/255;

% Separate out color planes
%pxmwrite('img14.1','g',255,R+1,cmap);
imwrite(R,'img14.1','pgm');
%pxmwrite('img14.2','g',255,G+1,cmap);
imwrite(G,'img14.2','pgm');
%pxmwrite('img14.3','g',255,B+1,cmap);
imwrite(B,'img14.3','pgm');

image(X);
axis('image');

Nclasses = 2;
MaxSamples = 200;

disp('Specify upper left and lower right points of the first training region.')
GetRegion(R,G,B,'data1',MaxSamples);

disp('Specify upper left and lower right points of the second training region.')
GetRegion(R,G,B,'data2',MaxSamples);

fid = fopen('info_file','w');
fprintf(fid,'%d \n3\n',Nclasses);
fprintf(fid,'data1 %d\n',MaxSamples);
fprintf(fid,'data2 %d\n',MaxSamples);
fclose(fid);



