function [vec_seg_vw,vec_roi_vw]=volume2vec_3w(vw,seg_roi)
vec_roi_vw=vw(:);
seg_roi_vec=seg_roi(:);
Nm=find(seg_roi_vec==1);
vec_seg_vw=vec_roi_vw(Nm);