App下載

Python怎么操作Matlab的Mat格式數(shù)據(jù)

猿友 2021-07-29 10:52:53 瀏覽數(shù) (3109)
反饋

作為老牌的科學(xué)計(jì)算常用語言,matlab稱霸了這個(gè)行業(yè)(就算現(xiàn)在python風(fēng)頭正勁,matlab也還是能保住其地位)。作為后來者的python想要在科學(xué)計(jì)算這方面分一杯羹,就必須要支持matlab的很多操作,比如說操作matlab的mat格式數(shù)據(jù)。那么python怎么操作mat格式數(shù)據(jù)呢?接下來這篇文章告訴你!

1. 非 matlab v7.3 files 讀寫

import scipy.io as sio
import numpy
# matFile 讀取
matFile = 'matlabdata.mat'
datas = sio.loadmat(matFile)
# 加載 matFile 內(nèi)的數(shù)據(jù)
# 假設(shè) mat 內(nèi)保存的變量為 matlabdata
matlabdata = datas['matlabdata']

# matFile 寫入
save_matFile = 'save_matlabdata.mat'
save_matlabdata = np.array([1,2,3,4,5])
sio.savemat(save_matFile, {'array':save_matlabdata})

2. matlab v7.3 files 讀取

如果 matlab 保存 data 時(shí),采用的是 ‘-v7.3',scipy.io.loadmat函數(shù)加載數(shù)據(jù)會(huì)出現(xiàn)錯(cuò)誤:

File "/usr/local/lib/python2.7/dist-packages/scipy/io/matlab/mio.py", line 64, in mat_reader_factory
raise NotImplementedError('Please use HDF reader for matlab v7.3 files')
NotImplementedError: Please use HDF reader for matlab v7.3 files

可以采用:

import h5py
with h5py.File('matlabdata.mat', 'r') as f:
    f.keys() # matlabdata.mat 中的變量名
datas = h5py.File('matlabdata.mat')['matlabdata'].value

補(bǔ)充:【Matlab/Python】Matlab和Python之間的數(shù)據(jù)傳輸

很多時(shí)候,我們需要把matlab里的數(shù)據(jù)保存下來,然后用python來處理。

方法一(.mat格式)

直接將matlab的數(shù)據(jù)存儲(chǔ)成.mat格式,然后在python中利用scipy.io中的loadmat函數(shù)來讀取

import scipy.io     as sio
data = sio.loadmat(...)

方法二(.h5格式)

當(dāng)要存儲(chǔ)的.mat文件比較大時(shí),matlab中需要用save -v7.3才能存儲(chǔ)。但是利用方法一在python中讀取時(shí),會(huì)不支持。可以用以下方法讀取

with h5py.File("mydata.mat") as f:
    data = f["mydata"][:]

可以正確讀取數(shù)據(jù),但是數(shù)組維度會(huì)倒過來,即本來是(2,3,4,5),讀出來會(huì)是(5,4,3,2)

但是只有.mat數(shù)據(jù)用save -v7.3保存時(shí),才能用此方法讀取,否則,應(yīng)用方法一讀取

所以可以在matlab中將數(shù)據(jù)保存成.h5格式,統(tǒng)一按照方法二來讀取

matlab中存儲(chǔ).h5格式用如下方法

h5create('data.h5','/data',[2,3,4,5]);
data = rand(2,3,4,5)
h5write('data.h5','/data',data)

但是用python讀取時(shí),矩陣維度仍然會(huì)倒過來

以上就是python操作mat格式數(shù)據(jù)的全部內(nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持W3Cschool。



0 人點(diǎn)贊