1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
function ConvertHDF5_VTK(hdf_file, vtk_prefix, varargin)
% ConvertHDF5_VTK(hdf_file, vtk_prefix, varargin)
%
% Convert openEMS field data stored in the given hdf5 file to a vtk file.
%
% arguments:
% hdf_file: source hdf5 file
% vtk_prefix: output vtk files prefix
%
% optional arguments:
% 'TD_Dump': activate dump for time-domain data (default is off)
% 'FD_Dump': activate dump for frequency-domain data (default is on)
% 'NumPhase': number of phase to dump frequency domain data animation
% (default is 36 --> 10°)
% 'FieldName': field name written to vtk, e.g. 'E-Field'
% 'weight': field weighting
%
% for more optional aguments have a look at ReadHDF5Dump
%
% example:
% % read time-domian data from hdf5, perform dft and dump as vtk
% ConvertHDF5_VTK('Et.h5','Ef','NumPhase',18,'Frequency',1e9)
%
% openEMS matlab interface
% -----------------------
% author: Thorsten Liebig
%
% See also ReadHDF5Dump Dump2VTK
do_FD_dump = 1;
do_TD_dump = 0;
phase_N = 36;
weight = 1;
fieldname = 'unknown';
for n=1:2:numel(varargin)
if (strcmp(varargin{n},'TD_Dump')==1);
do_TD_dump = varargin{n+1};
elseif (strcmp(varargin{n},'FD_Dump')==1);
do_FD_dump = varargin{n+1};
elseif (strcmp(varargin{n},'NumPhase')==1);
phase_N = varargin{n+1};
elseif (strcmp(varargin{n},'FieldName')==1);
fieldname = varargin{n+1};
elseif (strcmp(varargin{n},'weight')==1);
weight = varargin{n+1};
end
end
[field mesh] = ReadHDF5Dump(hdf_file, varargin{:});
if ((do_TD_dump==0) && (do_FD_dump==0))
warning('openEMS:ConvertHDF5_VTK','FD and TD dump disabled, nothing to be done...');
end
if (do_FD_dump)
if (~isfield(field,'FD'))
warning('openEMS:ConvertHDF5_VTK','no FD data found skipping frequency domian vtk dump...');
else
%set weighting
if (numel(weight)~=numel(field.FD.frequency))
FD_weight = ones(size(field.FD.frequency))*weight(1);
else
FD_weight = weight;
end
if (field.FD.DataType==1) % dump complex value FD data
ph = linspace(0,360,phase_N+1);
ph = ph(1:end-1);
for n = 1:numel(field.FD.frequency)
for p = ph
filename = [vtk_prefix '_' num2str(field.FD.frequency(n)) '_' num2str(p,'%03d') '.vtk' ];
Dump2VTK(filename, real(FD_weight(n)*field.FD.values{n}*exp(1j*p*pi/180)), mesh, fieldname, varargin{:});
end
filename = [vtk_prefix '_' num2str(field.FD.frequency(n)) '_abs.vtk' ];
Dump2VTK(filename, abs(FD_weight(n)*field.FD.values{n}), mesh, fieldname, varargin{:});
filename = [vtk_prefix '_' num2str(field.FD.frequency(n)) '_ang.vtk' ];
Dump2VTK(filename, angle(FD_weight(n)*field.FD.values{n}), mesh, fieldname, varargin{:});
end
else % dump real value FD data
for n = 1:numel(field.FD.frequency)
filename = [vtk_prefix '_' num2str(field.FD.frequency(n)) '.vtk' ];
Dump2VTK(filename, real(FD_weight(n)*field.FD.values{n}), mesh, fieldname, varargin{:});
end
end
end
end
if (do_TD_dump)
if (~isfield(field,'TD'))
warning('openEMS:ConvertHDF5_VTK','no TD data found skipping time domian vtk dump...');
else
disp('dumping time domain data...')
acc = ['%0' int2str(ceil(log10(numel(field.TD.time)+1))) 'd'];
for n = 1:numel(field.TD.time)
filename = [vtk_prefix '_TD_' num2str(n,acc) '.vtk' ];
Dump2VTK(filename, abs(weight(1))*field.TD.values{n}, mesh, fieldname, varargin{:});
end
end
end
|
