概要

MATLABのデータをnumpyで利用したかったので、scipyで読み込んでみた。

Octaveでデータの作成

MATLABが手元になかったので、Octaveでm.matを作成。

$ octave -q
octave:1> m = zeros(8, 8)
m =

   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0

octave:2> m(3:6, 3:6) = 1
m =

   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0
   0   0   1   1   1   1   0   0
   0   0   1   1   1   1   0   0
   0   0   1   1   1   1   0   0
   0   0   1   1   1   1   0   0
   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0

octave:3> save data.mat m -mat
octave:4> save data.txt m

ちなみに -mat オプションを付けない場合はテキストデータなので、コンソールで読める。

$ cat data.txt
# Created by Octave 4.2.1, Fri Dec 08 01:32:35 2017 JST <user@hostname>
# name: m
# type: matrix
# rows: 8
# columns: 8
 0 0 0 0 0 0 0 0
 0 0 0 0 0 0 0 0
 0 0 1 1 1 1 0 0
 0 0 1 1 1 1 0 0
 0 0 1 1 1 1 0 0
 0 0 1 1 1 1 0 0
 0 0 0 0 0 0 0 0
 0 0 0 0 0 0 0 0

scipyでMATLAB形式のデータを読み込み

scipy.ioを使用しdata.matを読みこむ。

$ ipython
Python 3.6.1 (v3.6.1:69c0db5050, Mar 21 2017, 01:21:04)
Type 'copyright', 'credits' or 'license' for more information
IPython 6.2.1 -- An enhanced Interactive Python. Type '?' for help.

In [1]: import scipy.io

In [2]: data = scipy.io.loadmat("data.mat")

In [3]: type(data["m"])
Out[3]: numpy.ndarray

In [4]: data["m"]
Out[4]:
array([[ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  1.,  1.,  1.,  1.,  0.,  0.],
       [ 0.,  0.,  1.,  1.,  1.,  1.,  0.,  0.],
       [ 0.,  0.,  1.,  1.,  1.,  1.,  0.,  0.],
       [ 0.,  0.,  1.,  1.,  1.,  1.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.]])

↑のようにして読み込むことができた。

scipyでMATLAB形式のデータを書き込み

逆に書き込むこともできる。 行列データmの中身をすべて-1したものをdata-1.matとして書き込む。

In [5]: data["m"] = data["m"] - 1

In [6]: data["m"]
Out[6]:
array([[-1., -1., -1., -1., -1., -1., -1., -1.],
       [-1., -1., -1., -1., -1., -1., -1., -1.],
       [-1., -1.,  0.,  0.,  0.,  0., -1., -1.],
       [-1., -1.,  0.,  0.,  0.,  0., -1., -1.],
       [-1., -1.,  0.,  0.,  0.,  0., -1., -1.],
       [-1., -1.,  0.,  0.,  0.,  0., -1., -1.],
       [-1., -1., -1., -1., -1., -1., -1., -1.],
       [-1., -1., -1., -1., -1., -1., -1., -1.]])

In [7]: scipy.io.savemat("data-1.mat", data)

OctaveでMATLAB形式のデータを読み込み

上記でscipyによって作成したdata-1.matをOctaveで読み込む。

$ octave -q
octave:1> load data-1.mat m -mat
octave:2> m
m =

  -1  -1  -1  -1  -1  -1  -1  -1
  -1  -1  -1  -1  -1  -1  -1  -1
  -1  -1   0   0   0   0  -1  -1
  -1  -1   0   0   0   0  -1  -1
  -1  -1   0   0   0   0  -1  -1
  -1  -1   0   0   0   0  -1  -1
  -1  -1  -1  -1  -1  -1  -1  -1
  -1  -1  -1  -1  -1  -1  -1  -1

読み込めた。 めでたしめでたし。

参考

• https://docs.scipy.org/doc/scipy/reference/tutorial/io.html