画像内の各色毎の割合算出方法について
显示 更早的评论
深層学習を使用したセマンティックセグメンテーションでの画像のセグメント後に、各ラベル毎の画像内に占める割合を算出したいのですが方法がわかりません。解決策がありましたらご回答の程よろしくお願い致します。
采纳的回答
まず、質問のリンク先から入力データ(画像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
更多回答（1 个）
2 个评论
Atsushi Ueno
2021-11-28
私はMATLAB Answersのlive editorで実行する事だけを考えた為、その点回答を端折ってしまいました。countEachLabel関数の対象は(イメージ)データストアです。画像1枚分のデータストアに対しては1枚分の結果を表示します。しかしそれでは本質問の主旨に対する回答になりませんよね。
質問の趣旨(複数の画像が含まれるデータストアの処理に関して、画像1枚だけの結果を表示したい)に対しては、データストアを分割する必要があります。その点後ほど1番目の回答に追記致します。
Atsushi Ueno
2021-11-29
そうですね。下記のように色々と疑問点があります。
skyClassIds = uint8([1]); これはSkyのラベルIDが1ではないからなのでしょうか⇒そう思います
skyPixels = ismember(y(:),skyClassIds); 画素値とクラスIDを比較する根拠は？
validPixels = (y~=0); ← 単に縦×横で良く、画素値=0を除く必要は無いのでは？
numskyPixels = sum(skyPixels(:)) ← カラー画像なので3倍の値になる
numValidPixels = sum(validPixels(:)) ← カラー画像なので3倍の値になる
せっかくImage Processing Toolboxを使ってるのですから、回答に挙げたcountEachLabel関数を使うのはいかがでしょうか。
类别
在 帮助中心 和 File Exchange 中查找有关 セマンティック セグメンテーション 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
