How to replace the data in tif file with my own data by matlab?

Hello,

I have .tif file, how I can replace or edit its data with my own matrix of data?

I cc'd mymatrix and the .tif file.

If you have an idea, please elaborat a littile bit on explanation how to do it, since I am beginner with matlab.

Cheers

Answer 1

[Walter Roberson]

The file you have appears to be a GeoTIFF file, which you can read with readgeoraster or the older geotiffread

You can write out the modified data using geotiffwrite

geotiffinfo might also be of interest.

[Rafat]

Can you help me coding that?

[Walter Roberson]

D=unzip('Test5.zip');
%%
[D1, R] = readgeoraster(D{1,1});

DS = load(D{1,3}).image_median;
subplot(211);geoshow(double(D1), R); title('doy2019190052709_aid0001.tif')
subplot(212); geoshow(double(DS), R); title('mymatrix.mat')

ind = DS>0.9;
D1(ind) = DS(ind);
%%
figure
geoshow(double(D1), R); title('Updated TIFF file')

geotiffwrite("UPDATED.tif", D1, R)

[Rafat]

@Walter Roberson

I test it, it's perfect!!

If you please, is this 3 or 2 in the following comand line? I think it's 2?

DS = load(D{1,3}).image_median;

[Walter Roberson]

The file appears to contain a __MACOSX subdirectory.

pwd

ans = '/users/mss.system.Bh1JYo'

ls -la

total 76
drwxr-xr-x 2 mluser worker  4096 Nov  7 01:53 .
drwx------ 1 mluser worker  4096 Nov  7 01:52 ..
-rw-r--r-- 1 mluser worker 68764 Nov  7 01:53 Test5.zip

!zipinfo -1 Test5.zip

doy2019190052709_aid0001.tif
__MACOSX/._doy2019190052709_aid0001.tif
mymatrix.mat
__MACOSX/._mymatrix.mat

D=unzip('Test5.zip')

D = 1×4 cell array
    {'doy2019190052709_aid0001.tif'}    {'__MACOSX/._doy2019190052709_aid0001.tif'}    {'mymatrix.mat'}    {'__MACOSX/._mymatrix.mat'}

ls -la

total 148
drwxr-xr-x 3 mluser worker  4096 Nov  7 01:53 .
drwx------ 1 mluser worker  4096 Nov  7 01:52 ..
-rw-r--r-- 1 mluser worker 40220 Jul 19 19:45 doy2019190052709_aid0001.tif
drwxr-xr-x 2 mluser worker  4096 Nov  7 01:53 __MACOSX
-rw-r--r-- 1 mluser worker 27157 Nov  6 16:00 mymatrix.mat
-rw-r--r-- 1 mluser worker 68764 Nov  7 01:53 Test5.zip

Answer 2

[Sulaymon Eshkabilov]

If understood your question correctly, this is how to get it solved:

D=unzip('Test5.zip');
%%
D1 = imread(D{1,1});
DS = load(D{1,3}).image_median;
subplot(211);imagesc(D1); title('doy2019190052709_aid0001.tif')
subplot(212); imagesc(DS); title('mymatrix.mat')

[X, Y] = find(DS(:,:,1)>0.9);
D1(X,Y,:)=DS(X,Y,:);
%%
figure
imagesc(D1); title('Updated TIFF file')
saveas(gcf, 'UPDATED.tiff') % Save the updated image data

[Rafat]

Hi,

Many thanks!!

The idea behiend copy 'matrix.mat' to 'doy2019190052709_aid0001.tif', is that the tif have geografical information that I want to use it as template.

In your code there is no any more geografical information, so I cannot use it in qgis, ther is no lat/long any more.

Since I am biggner so maybe, I am asking no logic qustion, or may be impossible to do what I think, sorry.

[Rafat]

@Sulaymon Eshkabilov,

If you please what is the meaning of the following comand line? is it 2 or 3?

DS = load(D{1,3}).image_median;

Cheers

[Sulaymon Eshkabilov]

read variable called image_median from the mat file.

[Rafat]

@Sulaymon Eshkabilov Thank you so much. Cheers

[Sulaymon Eshkabilov]

Most welcome. Glad to shed some light :)