load('peppers_RGB.mat');
imshow(RGB);
title('Peppers Image (RGB)');

t = Tiff('myfile.tif','w');  

tagstruct.ImageLength = size(RGB,1); 
tagstruct.ImageWidth = size(RGB,2);
tagstruct.Photometric = Tiff.Photometric.RGB;
tagstruct.BitsPerSample = 8;
tagstruct.SamplesPerPixel = 3;
tagstruct.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky; 
tagstruct.Software = 'MATLAB'; 

setTag(t,tagstruct)

write(t,RGB);
close(t);

load('peppers_YCbCr.mat');
imshow(Y);
title('Peppers Image (Y Component)');

t = Tiff('myfile_YCbCr.tif','w');  

tagstruct.ImageLength = size(Y,1);
tagstruct.ImageWidth = size(Y,2);
tagstruct.SampleFormat = 1; % uint
tagstruct.Photometric = Tiff.Photometric.YCbCr;
tagstruct.BitsPerSample = 8 ;
tagstruct.SamplesPerPixel = 3;
tagstruct.YCbCrSubSampling = [1,1];
tagstruct.Compression = Tiff.Compression.None;  
tagstruct.PlanarConfiguration = Tiff.PlanarConfiguration.Chunky; 
tagstruct.Software = 'MATLAB'; 

setTag(t,tagstruct)

write(t,Y,Cb,Cr)
close(t)