Las点云数据读取代码

function outfile = LASRead （infilenameoutfilenamenFields）
% LASREAD reads in files in LAS 1.1 format and outputs comma delimited text files
%
% INPUT
% infilename:   input file name in LAS 1.1 format 
%               （for example ‘myinfile.las‘） 
% outfilename:  output file name in text format 
%               （for example ‘myoutfile.txt‘）
% nFields:      default value of 1 outputs X Y and Z coordinates of the 
%               point - [X Y Z]. 
%               A value of 2 gives Intensity as an additional attribute - [X Y Z I].
%               A value of 3 gives the Return number and the Number of returns 
%               in addition to the above - [X Y Z I R N].                
%           
% OUTPUT
% outfile:      the output matrix
% 
% EXAMPLE
% A = LASRead （‘infile.las‘ ‘outfile.txt‘ 3）
%
% Cici Alexander
% September 2008 （updated 26.09.2008）

% Open the file
fid =fopen（infilename）;

% Check whether the file is valid
if fid == -1
    error（‘Error opening file‘）
end

% Check whether the LAS format is 1.1
% fseek（fid 24 ‘bof‘）;
% VersionMajor = fread（fid1‘uchar‘）;
% VersionMinor = fread（fid1‘uchar‘）;
% if VersionMajor ~= 1 || VersionMinor ~= 1
%     error（‘LAS format is not 1.1‘）
% end

% Read in the offset to point data
fseek（fid 96 ‘bof‘）;
OffsetToPointData = fread（fid1‘uint32‘）;

% Read in the scale factors and offsets required to calculate the coordinates
fseek（fid 131 ‘bof‘）;
XScaleFactor = fread（fid1‘double‘）;
YScaleFactor = fread（fid1‘double‘）;
ZScaleFactor = fread（fid1‘double‘）;
XOffset = fread（fid1‘double‘）;
YOffset = fread（fid1‘double‘）;
ZOffset = fread（fid1‘double‘）;

% The number of bytes from the begin