まず、質問のリンク先から入力データ(画像2枚とピクセルラベルイメージデータ2枚の計4枚)を抜き出し、これらをイメージデータストア、ピクセルラベルデータストアとして読み込み、「深層学習を使用したセマンティックセグメンテーションでの画像のセグメント後」の例を作成しました。
mkdir labels % ピクセルラベルイメージデータ格納用のフォルダを作成（画像はカレントフォルダ)
movefile 0016E5_08430_L.png labels % ピクセルラベルイメージデータをlabelsフォルダに移動
movefile Seq05VD_f00840_L.png labels % ピクセルラベルイメージデータをlabelsフォルダに移動
cmap = [128 128 128;128 0 0;192 192 192;128 64 128;60 40 222;128 128 0;192 128 128;64 64 128;64 0 128;64 64 0;0 128 192] ./ 255;
classes = ["Sky" "Building" "Pole" "Road" "Pavement" "Tree" "SignSymbol" "Fence" "Car" "Pedestrian" "Bicyclist"];
labelIDs = {[128 128 128] [000 128 064; 128 000 000; 064 192 000; 064 000 064; 192 000 128] [192 192 128; 000 000 064] [128 064 128; 128 000 192; 192 000 064] [000 000 192; 064 192 128; 128 128 192] [128 128 000; 192 192 000] [192 128 128; 128 128 064; 000 064 064] [064 064 128] [064 000 128; 064 128 192; 192 128 192; 192 064 128; 128 064 064] [064 064 000; 192 128 064; 064 000 192; 064 128 064] [000 128 192; 192 000 192]};
imds = imageDatastore(pwd); % 添付した'Seq05VD_f00840.png'と'0016E5_08430.png'を読む
pxds = pixelLabelDatastore([pwd '/labels'],classes,labelIDs); % 添付した'Seq05VD_f00840_L.png'と'0016E5_08430_L.png'を読む
質問のリンク先における「データセット統計の解析」の内容を、データセットを分割し画像1枚分毎に実行しました。countEachLabel関数により各ラベル毎の画像内に占める割合を算出しました。
n = numpartitions(pxds);
for i = 1:n
subimds = partition(imds,n,i); % データセットを分割し1枚ずつ取り出す
subpxds = partition(pxds,n,i); % データセットを分割し1枚ずつ取り出す
I = histeq(readimage(subimds,1)); % 元画像の読み込みとコントラスト加工
C = readimage(subpxds,1); % ピクセルラベルイメージデータの読み込み
B = labeloverlay(I,C,'ColorMap',cmap); % 重ねる
figure; imshow(B); colormap(gca,cmap); % 表示する
c = colorbar('peer', gca); c.TickLabels = classes; numClasses = size(cmap,1); c.Ticks = 1/(numClasses*2):1/numClasses:1; c.TickLength = 0;
% ===================================================================================
% ここから下が、countEachLabel関数により各ラベル毎の画像内に占める割合を算出する方法です。
tbl = countEachLabel(subpxds) % 各ラベル毎の画素数及び画像全体の画素数を出力する
frequency = tbl.PixelCount/sum(tbl.PixelCount); % 各ラベル毎の画像内に占める割合を計算
figure; bar(1:numel(classes),frequency); % 表示する
xticks(1:numel(classes)); xticklabels(tbl.Name); xtickangle(45); ylabel('Frequency');
end
