Octaveで画像のマスク処理を行う

%%%%%画像マスク化
pkg load image;
img = rgb2gray(imread(IMAGEURL)); %%グレイスケール化
BWs = img > 20; %% 20より上の濃さのみ残す
BW2 = bwareaopen(BWs, 300); %% 面積300より小さいのを消したものをマスクとして設定
img(BW2 == 0) = 0; %%マスクをグレイケールイメージに適用
imshow(ROI); %%表示