How to read NanoDCAL output with Python?

vincentm

It is sometimes inconvenient to parse through output files with Matlab. Is it possible to do so with another platform, like Python?

vincentm

MAT-files are but HDF5 files. As such, one may read them with the h5py module as follows

import scipy.io as sio
import h5py

filename = "TotalEnergy.mat"
f = h5py.File(filename, "r")
print(f.keys())
print(f["data"][0])

which produces

<KeysViewHDF5 ['data']>
[-59.08137326]

vincentm

Another way would be to load the data in MATLAB and write it to a JSON file as follows

clear
load('NanodcalObject.mat');
s = struct(object);
f = fieldnames(s);
nf = numel(f);
for ii = 1:nf
   try
      jsonencode(s.(f{ii}));
   catch
      s = rmfield(s, f{ii});
      fprintf('remove field %s\n', f{ii});
   end
end
txt = jsonencode(s);
fid = fopen('NanodcalObject.json','w');
fprintf(fid, '%s', txt);
fclose(fid);

and then read it with Python

import json

filename = "NanodcalObject.json"
with open(filename, "r") as f:
    d = json.load(f)

print(d.keys())