openEMS/matlab/ConvertHDF5_VTK.m

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 augments 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