Hi, assuming your image variable is Im; I've created this small code, which provides a ThisSector image (B&W) and ThisSectorIDs (indexes to sector positions) which you can use to select part of your image to crop (or just multiply the image by the sector) to obtain black everywhere instead where the sector is one. Hope it helps. The imshow() etc. is just for you to see the produced sectors. The code allows you to interactively draw the circle around the area of interest. If you already have the circle defined, just comment the section where the circle is created interactively
[Nrows,Ncols,Ncolors]=size(Im);
hFig=imshow(Im);
%% Comment this section if x0,y0,r are already selected
hCircle=drawcircle;
x0=hCircle.Center(1);
y0=hCircle.Center(2);
r=hCircle.Radius;
%% if you already have x0,y0,r, then use
hCircle=drawcircle('Center',[x0,y0],'Radius',r);
%%
CircleMask=createMask(hCircle);
n=6;
tet=linspace(-pi,pi,n+1);
xi=r*cos(tet)+x0;
yi=r*sin(tet)+y0;
for k=1:numel(xi)-1
if (k<0.5*(numel(xi)))
CornerX1=xi(k);
CornerY1=1;
CornerX2=xi(k+1);
CornerY2=1;
else
CornerX1=xi(k);
CornerY1=Nrows;
CornerX2=xi(k+1);
CornerY2=Nrows;
end
hPoly=drawpolygon('Position',[x0 y0; xi(k) yi(k); CornerX1 CornerY1; CornerX2 CornerY2; xi(k+1) yi(k+1);x0 y0]);
ThisMask=createMask(hPoly);
ThisSector=(ThisMask & CircleMask);
ThisSectorIDs=find(ThisSector);
figure(k);
imshow(ThisSector);
end
