matlab中读取一个影像矩阵,以矩阵中的四个点为角点,画一个四边形,给这个四边形内部的点赋值为1
dem=imread(dem_name);
x = [66 10597 10987 425];
y = [3904 7 685 4693];
temp_mode=zeros(size(dem));
for i=1:size(dem,2)
disp(i/size(dem,2))
temp=inpolygon(ones(1,5043)*i, [1:1:5043], x,y);
temp_mode(:,i)=temp;
end
dem_new=dem.*temp_mode;