bim = bigimage('tumor_091R.tif');

h = bigimageshow(bim);
hROI = drawcircle(gca,'Radius',470,'Position',[1477 2284]);

maskLevel = 3;

ref = bim.SpatialReferencing(maskLevel);
pixelExtent = [ref.PixelExtentInWorldX,ref.PixelExtentInWorldY];

bmask = bigimage(ref,1,'logical');

for cStart = 1:bmask.BlockSize(2):ref.ImageSize(2)
    for rStart = 1:bmask.BlockSize(1):ref.ImageSize(1)

        % Get the center of top left pixel of this block in world units.
        xyStart = [cStart,rStart].*pixelExtent;

        % Get the block size. The |'BlockSize'| property represents the
        % size as a 2-element vector of the form [row,column]. Switch the
        % order of the elements so that the block size is represented as
        % [x,y].
        bsize = bmask.BlockSize
        bsize = fliplr(bsize);

        % Get the center of the bottom right pixel of this block in world
        % units.
        xyEnd = ([cStart,rStart] + (bsize-1)).*pixelExtent;

        % Use the |meshgrid| function to determine the (x,y) coordinates in
        % world units for all pixels in this block.
        [xgrid, ygrid] = meshgrid(xyStart(1):pixelExtent(1):xyEnd(1),...
            xyStart(2):pixelExtent(2):xyEnd(2));

        % Determine which pixels have coordinates inside the ROI.
        roiPositions = hROI.Vertices;
        blockMask = inpolygon(xgrid,ygrid, ...
            roiPositions(:,1),roiPositions(:,2));

        % Set the pixel values of the block.
        setBlock(bmask,1,xyStart,blockMask);
    end
end

figure
bigimageshow(bmask)