How to make tif output file with Coordinate Reference System (CRS)?

Hello,

I am working on tif files to produce the median matrix for them, I am using the following cod for that.

How I can make the output as .tif with keep the coordinate refrence system (CRS)?

unzip Test3.zip
% Warning: Cannot open file "" for reading.
files = dir('*.tif');
filenames = fullfile({files.folder},{files.name});
N = numel(files);
sz = zeros(N,2);
for ii = 1:N
    files(ii).data = imread(filenames{ii});
    sz(ii,:) = size(files(ii).data);
end

datablock = cat(3, files.data);
image_median = median(datablock, 3, 'omitnan');
[min(image_median(:)), max(image_median(:))]
image(image_median)

Answer 1

[Walter Roberson]

Make sure the CRS is a "map raster reference object" and use geotiffwrite from the Mapping Toolbox.

[Rafat]

Hello Walter,

Many thanks for that, but I am beginner on that, it's not easy for me!! sorry.

Cheers