Matlab二值图像栅格化和圆域范围框定

本文以2017年华为杯数学建模A题中震区高程数据，实现对海拔高于3000m的震区做阈值截断，将海拔大于等于3000m的地区赋予新值255，将海拔小于3000m地区赋予新值0，从而获得震区3000m海拔高度下的二值图像。震区图和3000m海拔高度下的二值图像分别如图1，图2 所示：

··················································图1 震区三维高程地形图

··················································图2 震区3000m海拔高度下的二值图像

1.给二值图像添加网格线(辅助线)或叫栅格化:

2.在1的基础之上对圆域范围框定:

代码:

% 重要参数设置
r = 69;
c = [252,323]
max(c(1)-r,1)
min(c(1)+r,n)
max(c(2)-r,1)
min(c(2)+r,m)
% 二值图像中显示出障碍区域大小
I=imread(&