0001
0002 function [roadPoints, costMap, Green, Water] = extractRoadsCoverWaterPx(A)
0003
0004
0005
0006
0007
0008
0009
0010
0011 water = 1;
0012 road = 1;
0013 urabanArea = 1;
0014 grass = 1;
0015
0016
0017 roadPoints = zeros(size(A(:,:,1)));
0018 costMap= ones(size(A(:,:,1)));
0019 Green = zeros(size(A(:,:,1)));
0020 Water = zeros(size(A(:,:,1)));
0021
0022 for i = 1:length(A(:,1,1))
0023 for j=1:length(A(1,:,1))
0024 costMap(i,j) = urabanArea;
0025 [col] = findColor(A(round(i),round(j),:));
0026 if (col(2) == 1)
0027 costMap(i,j) = grass;
0028 Green(i,j) = 1;
0029 end;
0030 if (col(3) == 1)
0031 costMap(i,j) = water;
0032 Water(i,j) = 1;
0033 end;
0034 if (col(5) == 1)
0035 roadPoints(i,j) = 1;
0036 costMap(i,j) = road;
0037 end;
0038
0039 if costMap(i,j) == 0
0040 costMap(i,j) = grass;
0041 end;
0042
0043
0044
0045 col = zeros(size(col));
0046 end
0047 end
0048
0049
0050
0051