导出三维面数据.docxVIP

导出三维面数据所示为同三维面(3DFACE)组成的企鹅模型，每个小的面片为一个三角形，要求：（1） 将每个三角形的坐标读出并保存到文本文件中，文本文件中的一行为一个三角形3个点的坐标（一共9个数据），坐标数值保留到小数点后3位；（2） 相同颜色的面连续输出，不同颜色的面不能够交叉输出。图0-1 小企鹅解题思路将小企鹅所在模型空间设为当前，通过引用Thisdrawing.ModelSpace.Item(i)可遍历几何元素。本题中感兴趣的是3dface,其颜色属性color为Autocad颜色索引表中的某个值(0~255)。为保证连续输出相同颜色面，且不同颜色面不可交叉，需要对模型中的各3dface的颜色排序。为使编程简化，将问题解决过程分成3部分：1. 导出模型中所有3dface的有益信息这些信息包括：ModelSpace.Item变量号i，3dface颜色索引值，各顶点三维坐标值。此过程为下面的Matlab分析做准备，导出的数据保存在文件3dface.dat中。2. 颜色排序测试根据1中导出的数据，测试使用快速排序法对颜色排序的效果（准确性）。此过程可积累相关算法经验，以便在VBA中实现。图1-1 MATLAB重回小企鹅图1-2颜色排序效果3. 在VBA中的最终实现由过程2验证排序算法实现要求2的准确性，并且积累了快速排序算法的实现经验，此时重新返回VBA编制满足要求的程序。图1-3 程序过程总图程序1. 初始程序Public Sub read3dface() Dim ent As AcadEntity Dim x As Double, y As Double, z As Double Dim i As Integer, j As Integer, N As Integer Dim K As Integer, thisColor As Integer Dim Coord() As Double 将实际路径改变到当前文件夹下 ChDrive ThisDrawing.Path ChDir ThisDrawing.Path 输出必要数据，以便使用Matlab测试 Open 3dface.dat For Output As #1 遍历所有ModelSpace中几何体 N = ThisDrawing.ModelSpace.Count - 1 For i = 0 To N Set ent = ThisDrawing.ModelSpace.Item(i) If TypeOf ent Is Acad3DFace Then 输出顺序为 index color coordinates thisColor = ent.color coords = ent.Coordinates Print #1, i, thisColor, K = UBound(coords) For j = 0 To K - 1 Print #1, coords(j), Next j Print #1, coords(K) End If Next i Close End Sub 2. 颜色排序测试A 主程序% 测试来自cad的网格数据，绘制，并排序保存 %% 导入数据 clear;clc % 小企鹅顶点数据 load my3dface; % cad颜色索引与RGB对应表 load cadcolormap; % 三角形数目 N = length(face3d.index); %% 绘制小企鹅 figure(1) hold on axis equal for i = 1:N colorspec = colormap.rgb(face3d.color(i),:); X = face3d.coords(i,1:3:7); Y = face3d.coords(i,2:3:8); Z = face3d.coords(i,3:3:9); fill3(X,Y,Z,colorspec); end %% 排序，相同颜色合并 colors = face3d.color; index = face3d.index; [colors,index] = quick_sort(colors,index,1,length(index)); % 绘制颜色序列 figure(2) subplot(2,1,1) drawColorOrder(face3d.color,colormap,3dface原来的顺序); subplot(2,1,2) drawColorOrder(colors,colormap,3dface整理后的顺序); %% 导出文件 fname = myface3d.dat; fid = fopen(fname,w); for i=1:length(index) fprintf(fid,