比如以图像中某点为中心,取到中心距离<n的点集
for循环来计算距离比较耗时
x = 1:rows; //rows为图像行数
y = 1:cols; //cols为图像列数
[xx,yy] = meshgrid(y,x);
dist = (yy-row).^2 + (xx-col).^2; //(row, col)为指定为中心点
dist = sqrt(dist);
inds = find(dist < n); //到中心点距离< n个像素的index(按列扫描顺序)