首先準備乙個mask，reslice到與fmri相同的解析度。注意這個mask裡面不是roi的地方值為零，而且每個roi體素的值都是這個roi對應的編號。比如我用的power 264 roi，roi值最大是264。
然後用matlab寫**。把fmri從四維轉成二維去做會比較方便。
fmri =
spm_vol
('ica_denoised_sub-600001103037_task-rest_bold.nii');
fdata =
spm_read_vols
(fmri)
;mask =
spm_vol
('rroi_264_mni_1mm.nii');
mdata =
spm_read_vols
(mask)
;d =
size
(fdata,4)
;time =
124;
if d == time
fdatatran =
reshape
(fdata,
,d);
mdata =
reshape
(mdata,
,1);
% 將mask轉換為乙個列向量，當成索引表
slide264 =
zeros
(264
,time)
;for j =1:
264
roi_index =
find
(mdata==j)
;slide264
(j,:)=
mean
(fdatatran
(roi_index,:)
);end
end
金融文字時間提取
瀏覽github發現乙個有趣的庫，在裡面發現了乙個非常非常有用的庫 from coconlp.extractor import extractor ex extractor times ex.extract time text print times time記錄一下。順便感謝開源大佬！以及 尋人微...
根據id提取fasta序列
bioperl讀入寫出fasta 要求根據序列id，從fasta檔案中提取目標序列並輸出 資料序列id fasta檔案 思路以序列id為鍵，構建雜湊 用bioperl讀入fasta，獲得序列id 如果id存在於雜湊中，輸出序列 die perl 0 unless ar 3 0程式名 use bio ...
CBHG模組 提取序列特徵
cbhg模組如下圖所示。首次提出在goggle的一篇文章 tacotron towards end to end speech synthesis 回到cbhg模組，該模組善於提取序列特徵。下面為模組步驟 1.輸入序列，先經過k個1 d卷積，第k個卷積核 filter 通道為k，這些卷積核可以對當前...