Question

4 MARKS QUEStION 3 Background You are part of a team working for the United Nations Environment Programme (UNEP) to investigate the deforestation process in Borneo. You are provided six images of the forest area in Borneo from 1950-2020 which comprise of historical and projection data. Forests are represented as green pixels and deforested areas as yellow pixels. Q3a In the Q3a.m file, use the imread() function to read the images. For each year (1950, 1985, 2000, 2005, 2010, 2020), calculate the percentage forested area (PFA) which is defined as: forest area green pixels PFA 10096-greenandyellow pixe >< 100% green and Write the year, forested area, total land area and PFA to a file named "BorneoForestData.txt". The file should look like the following Forest Area T tal Land Area (%) Year 1950 1985 PEA Hint: Green (forest) in grayscale 75 to 115 (inclusive) Light blue (water) in grayscale 240 Black (text) in grayscale 0 All other colours in grayscale represent land You should still have five figure windows by the end of this task. Q3b In the Q3b.m file, plot the PFA against the year as blue circles in a new figure. You believe that the data can be appropriately represented by either an exponential model or a 2nd order polynomial model. Find the coefficients for each model based on the form of the equations below: Exponential model: PFA-aePt 2nd order polynomial model: PFA at2 bt c where a, B, a, b, c are coefficients, and t represents the year Plot both fitted models for years 1950 to 2100 (yearly increments) on the newly created figure. Print the equations of the fitted models to the command window using the exponential specifier for fprintf You should have six figure windows by the end of this task.

Answer 1

%Read the image
I1950 = imread('1950.png');
I1985 = imread('1985.png');
I2000 = imread('2000.png');
I2005 = imread('2005.png');
I2010 = imread('2010.png');
I2020 = imread('2020.png');

%Display the image
subplot(3, 2, 1);
imshow(I1950);
subplot(3, 2, 2);
imshow(I1985);
subplot(3, 2, 3);
imshow(I2000);
subplot(3, 2, 4);
imshow(I2005);
subplot(3, 2, 5);
imshow(I2010);
subplot(3, 2, 6);
imshow(I2020);

%Identify Green and Yellow Pixels for image 1950
countG = 0;
countY = 0;
[x, y, z] = size(I1950);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I1950(i, j, k)>=75 && I1950(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I1950(i, j, k)~=0 && I1950(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fileID = fopen('BorneoForestData.txt', 'w');
fprintf(fileID, 'Year Forest Area Total Land Area PFA ');
fprintf(fileID, '%d %d %d %f ', 1950, countG, countG+countY, PFA);

%Identify Green and Yellow Pixels for image 1985
countG = 0;
countY = 0;
[x, y, z] = size(I1985);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I1985(i, j, k)>=75 && I1985(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I1985(i, j, k)~=0 && I1985(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fprintf(fileID, '%d %d %d %f ', 1985, countG, countG+countY, PFA);

%Identify Green and Yellow Pixels for image 2000
countG = 0;
countY = 0;
[x, y, z] = size(I2000);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I2000(i, j, k)>=75 && I2000(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I2000(i, j, k)~=0 && I2000(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fprintf(fileID, '%d %d %d %f ', 2000, countG, countG+countY, PFA);

%Identify Green and Yellow Pixels for image 2005
countG = 0;
countY = 0;
[x, y, z] = size(I2005);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I2005(i, j, k)>=75 && I2005(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I2005(i, j, k)~=0 && I2005(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fprintf(fileID, '%d %d %d %f ', 2005, countG, countG+countY, PFA);

%Identify Green and Yellow Pixels for image 2010
countG = 0;
countY = 0;
[x, y, z] = size(I2010);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I2010(i, j, k)>=75 && I2010(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I2010(i, j, k)~=0 && I2010(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fprintf(fileID, '%d %d %d %f ', 2010, countG, countG+countY, PFA);

%Identify Green and Yellow Pixels for image 2020
countG = 0;
countY = 0;
[x, y, z] = size(I2020);
for i = 1:x
for j = 1:y
for k = 1:z
%Check for Green Pixels
if I2020(i, j, k)>=75 && I2020(i, j, k)<=115
countG = countG+1;
  
%Check for Yellow Pixels
elseif I2020(i, j, k)~=0 && I2020(i, j, k)~=240
countY = countY+1;
end
end
end
end

%Calculate PFA
PFA = countG*100/(countG+countY);

%Write to file
fprintf(fileID, '%d %d %d %f ', 2020, countG, countG+countY, PFA);

Screenshot of BorneoForestData.txt

MATLAB R2017a HOME PLOTS EDITOR lSearch Documentation Log In Find Files 凶Compare ▼ 연 Run Section comment 9e骗3, Go To New Open Save Breakpoints Run Run and Advance Advance Run and Print ▼ Find Indent Time FILE Program FilesトMATLAB R2017aトbin Command Window > O3a > type BorneoForestData.txt Year Forest Area Total Landd Area PEA 1950 119365 1985 2000 64610 2005 2010 46135 2020 36679 1066012 1067328 1066839 1065192 1066176 4.327147 1065497 11.197341 6.817398 6.056209 5.293600 72764 56387 3.442431 02:16 OType here to search x ENG 03.01-2019 TO

Screenshot of Image: