第11讲：三维点云处理

Halcom

三维点云处理：百度网盘链接：https://pan.baidu.com/s/1YG6tUTWXCy4zGQhkbqYThA
具体链接在halcom.cn论坛，联系人QQ：3283892722
该论坛是一个学习交流平台，我会逐一的和大家分享学习。
欢迎大家录制视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1

主程序如下：

1. %% 人头部
   
2. clc,clear,close all
   
3. filename = 'head1.cgo';
   
4. [data,count]= cgo_XYZ(filename);
   
5. x = data(:,1);
   
6. y = data(:,2);
   
7. z = data(:,3);
   
8. figure('color',[1,1,1])
   
9. plot3(x,y,z,'.')
   
10. view([45.5 8]);
    
11. xlabel('x');ylabel('y');zlabel('z');

复制代码

cgo_XYZ函数如下：

1. function [data,count]= cgo_XYZ(filename)
   
2. % 输入：
   
3. %      filename:输入的.cgo文件名
   
4. % 输出：
   
5. %      data:获取的XYZ三维坐标信息数据
   
6. %      count:数据长度  
   
7. %% 读取.cgo数据，XYZ三维坐标信息
   
8. % fid = fopen('y.cgo');
   
9. fid = fopen(filename);
   
10. tline = fgetl(fid);
    
11. count = 1;
    
12. data=[];
    
13. while ischar(tline)
    
14.     str1 = tline;
    
15.     kk=[];
    
16.     for i=1:length(str1)
    
17.         if strcmp(str1(i),' ')
    
18.             kk=[kk,i];
    
19.         end
    
20.     end
    
21.     if length(kk)>=2
    
22.         str11 = str1(1:kk(1)-1);               str11 = str2num(str11);     
    
23.         str12 = str1(kk(1)+1:kk(2)-1);         str12 = str2num(str12);  
    
24.         str13 = str1(kk(2)+1:length(str1));    str13 = str2num(str13);  
    
25.         data = [data;str11,str12,str13];
    
26.     end
    
27.     tline = fgetl(fid);
    
28.     count = count+1;
    
29. end

复制代码

得到的图像如下：