[x,y] = ind2sub(size(masked(:,:,1)), conncomps{1}.PixelIdxList{1});
figure(1)
subplot(1,2,1)
imshow(masked(:,:,1));
subplot(1,2,2)
imshow(masked(:,:,1));
hold on
plot(y,x,'rx')
plot(round(F{1,1}(4)),round(F{1,1}(5)), 'bx')

map1 = cell(size(F,1),1);