v3d_import_in.m

matlab 三维显示。所有代码都有。可以显示三维网格。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [x,y,z,v] = v3d_import_in(filename)
%
% Prof. Spitzers IN-Datei importieren
% --------------------------------------------------------------------
%
% Format der Datei:
% 
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
fid = fopen(filename,'rt');

% 1 Zeile einlesen -> erstmal uninteressant
fgetl(fid);

% Zahlen an n鋍hster Stelle der Variable anzsrc und anzrec zuweisen
anzsrc=fscanf(fid,'%4f',1);
anzrec=fscanf(fid,'%4f',1);
fgetl(fid);

% Zahl an n鋍hster Stelle der Variable anzx zuweisen
% usw.
anzx=fscanf(fid,'%4f',1);
anzy=fscanf(fid,'%4f',1);
anzz=fscanf(fid,'%4f',1);
fgetl(fid);

% Quellen, Kommentar und Empf鋘ger einlesen -> uninteressant
for i=1:anzsrc
    fgetl(fid);
end
fgetl(fid);
for i=1:anzrec
    fgetl(fid);
end

% Blockkoordinaten einlesen
for i=1:anzx
    xx(i)=fscanf(fid,'%8f',1);
end
for i=1:anzy
    yy(i)=fscanf(fid,'%8f',1);
end
for i=1:anzz
    zz(i)=fscanf(fid,'%8f',1);
end
fgetl(fid);

% n鋍hste Zeile komplett einlesen, auseinander nehmen und res zuweisen
res=str2num(fgetl(fid));

% einzeln Z-Layer einlesen
w=waitbar(0,'Lade...');
for k=1:length(zz)-1
    waitbar(k/(length(zz)-1));
    for j=1:length(yy)-1
        for i=1:length(xx)-1
            c=fscanf(fid,'%c',1);
%           switch c
%               case {'1','2','3','4','5','6','7','8','9'}
                    c=str2num(c);      
%           end
            vv(j,i,k)=res(c);
        end
        c=fscanf(fid,'%c',1);
    end
    fgetl(fid);
end
close(w); 

% Datei schliessen
fclose(fid);

[x,y,z,v]=v3d_import_block(xx,yy,zz,vv);