区域填充代码-MATLAB.doc

三、区域填充 %Fill_Test.m文件，这是主文件 clear all; clc; num GetFigure ; ET Create_ET num ; %输出桶 size ET ; for i 1:length ET T ET i .AET; if ~isempty T ET i .ystart; length T ; for k 1:length T T k ; end end end Point Fill ET ; figure ; scatter Point :,1 ,Point :,2 ,.,r ; %GetFigure.m文件，获取图形数据文件,数据源为程序自带的*.xls文件 function Figure GetFigure %从指定文件中读取图形数据矩阵 %图形数据矩阵的各列含义如下： %点编号 点坐标x 点坐标y 线段起点编号 线段终点编号 % [FileName,PathName,FilterIndex] uigetfile *.xls,select an office file ; [num,text,head] xlsread strcat PathName,FileName ; Figure num; end %Create_ET.m文件，创建边表文件 function ET Create_ET num %num共五列数据，每一列的数据含义如下： %点编号 点坐标x 点坐标y 线段起点编号 线段终点编号 %采用矩阵存储相同深度的有效边，ET i .AET即指向该矩阵 %构造桶,先确定桶深 Ymin min num :,3 ; Ymax max num :,3 ; Ymin Ymin 1 ;Ymax Ymax 1 ; %构造了空桶ET,同ystart字段为起点的y值,next为此层桶内的有效边矩阵 for i 1:Ymax-Ymin+1 % ET i .x -1; % ET i .ystart i+Ymin-1; % ET i .k1 -100; % ET i .AET []; ET i struct x,-1,ystart,i+Ymin-1,k1,-100,AET,[] ; end %构造有效边AET for i 1:length num :,4 %取出第i条线段的端点 p [num num i,4 ,[2,3] ;num num i,5 ,[2,3] ]; % break; %做两个端点的差 if p 1,2 p 2,2 t p 1,: ; p 1,: p 2,: ; p 2,: t; end %开始计算斜率倒数 d p 2,: -p 1,: ; if d 2 0 k1 Inf; else k1 d 1 /d 2 ; end %斜率为0的边不输入桶中 if k1 Inf continue; end %构造有效边, %x字段为边起始点横坐标,ymax为终点纵坐标,k1为斜率的倒数,next为具有相同起点纵坐标的下一条有效边 % AET.x p 1,1 ; % AET.ymax p 2,2 ; % AET.k1 k1; AET struct x,p 1,1 ,ymax,p 2,2 ,k1,k1 ; %定位有效边深度 Tindex p 1,2 -Ymin+1; T ET Tindex .AET; j 1;%j表示新节点的插入位置 if isempty T %若当前矩阵为空,则指定结构体矩阵的各个字段 T 1 .x AET.x; T 1 .ymax AET.ymax; T 1 .k1 AET.k1; else while j length T if T j .x AET.x||T j .x AET.xT j .k1 AET.k1 j j+1; else break; end end %将后序点后移 for k length T :-1:j T