clc;clear;close?all;
pointdata=importdata（‘CSite4_orig.txt‘）;??
%坐標最小值
minx=min（pointdata（:1））;
miny=min（pointdata（:2））;
minz=min（pointdata（:3））;
%?平移至原點
pointx=pointdata（:1）-minx;
pointy=pointdata（:2）-miny;
pointz=pointdata（:3）-minz;
%確定步長
step=max（pointy）/2048
%遍歷點云
X=min（pointx）:step:max（pointx）;
Y=min（pointy）:step:max（pointy）;
%生成曲面并平滑
[X1Y1]=meshgrid（XY）;
Z=griddata（pointxpointypointzX1Y1）;
[mn]=size（Z）;
TZ=zeros（mn）;
for?i=1:m
????for?j=1:n
????????TZ（ij）=Z（in-j+1）;
????end
end
maxTZ=max（max（TZ））;
minTZ=min（min（TZ））;
%進行灰度重新分配
I=uint8（TZ/（maxTZ-minTZ）*255）;
figure
imshow（I）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-26?22:43??點云柵格化\
?????文件????20722400??2001-12-16?22:55??點云柵格化\CSite4_orig.txt
?????文件?????????670??2018-06-26?16:54??點云柵格化\Rasterized.m