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