生成ply文件

function write_to_ply (name,re_s1)
[b1 ,b2]=size(re_s1);
%fid=fopen('SmallCar-transfer.ply','w');
fid=fopen(name,'w');
fprintf(fid,'ply\n');
fprintf(fid,'format ascii 1.0\n');
fprintf(fid,'element vertex ');
fprintf(fid,'%5d\n',b1);
fprintf(fid,'property float x\n');
fprintf(fid,'property float y\n');
fprintf(fid,'property float z\n');
fprintf(fid,'property uint16 reflectance\n');
% fprintf(fid,'property float nx\n');
% fprintf(fid,'property float ny\n');
% fprintf(fid,'property float nz\n');
fprintf(fid,'property uchar red\n');
fprintf(fid,'property uchar green\n');
fprintf(fid,'property uchar blue\n');
fprintf(fid,'end_header\n');
%re_s1(:,4:6)=round(re_s1(:,4:6));

for i=1:b1
    for j=1:b2
        if(j<4)
            fprintf(fid,'%12.4f',re_s1(i,j));
        else
            fprintf(fid,'%5u',re_s1(i,j));  
        end
    end
    if(i ~= b1)
        fprintf(fid,'\n');
    end
end
fclose(fid);
end