%%
clear all
clc
%
fidA=fopen('mulu.txt');
list =textscan(fidA, '%s');
[nn,mm]=size(list{1,1});
%读入数据
t=1;
a={};
A=[];
for j=1:nn
disp(j)
aName=list{1,1}{j,1};
fidB=fopen(aName);
a=textread(aName,'%s');
[NN,MM]=size(a);
A(1,2)=str2num(a{44,1}); %间隔
A(1,3)=str2num(a{39,1}); %点数
[row,col]=find(strcmp(a,'CSMNC'));
kk=row+1;
if kk>0
else
[row,col]=find(strcmp(a,'XJB'));
kk=row+1;
end
if a{31,1}=='EW'
for i=kk:NN %EW
A(i-kk+2,1)=str2num(a{i,1});
end
elseif a{31,1}=='NS'
for i=kk:NN
A(i-kk+2,2)=str2num(a{i,1}); %NS
end
else
for i=kk:NN
A(i-kk+2,3)=str2num(a{i,1}); %UD
end
end
[m,n]=size(A);
if rem(j/3,1)==0
%**************************文件的输出******************************
%输出文件
aName01=strcat(a{1,1}(1:6),a{2,1}(2:end),'.dat');
name{t,1}=aName01;
t=t+1;
b=tabulate(name);
[m1,n1]=size(b);
for k=1:m1
if aName01==b{k,1}
if b{k,2}>1
ll=b{k,2};
aName00=strcat(aName01(1:18),'_',num2str(ll),'.dat');
else
aName00=aName01;
end
end
end
% aName02=strcat(a{1,1}(1:6),a{2,1}(2:end),'.csv');
%*************CSV输出*****************
% csvwrite(aName02,A);
%*************fprintf输出*****************
fid=fopen(aName00,'wt');
fprintf(fid,'%12.2f\t',A(1,1));
fprintf(fid,'%12.4f\t',A(1,2));
fprintf(fid,'%12.0f\n',A(1,3));
for k1=2:1:m
for k2=1:1:n
if k2==n
fprintf(fid,'%12.6f\n',A(k1,k2));
else
fprintf(fid,'%12.6f\t',A(k1,k2));
end
end
end
fclose(fid);
a={};
A=[];
%*************dlmwrite输出*****************
% dlmwrite(aName02,A,'delimiter','\n','precision','%12.6f');
end
fclose(fidB);
end