In this paper, we propose a novel numerical technique to 2D multi-term time and space fractional Bloch–Torrey equations defined on an irregular convex domain. First, we consider the problem with space integral Laplacian operator. We present the semi-discrete and fully-discrete schemes by using the L1 formula on a temporal graded mesh and an unstructured-mesh Galerkin finite element method (FEM) based on the matrix transfer technique (MTT). Moreover, the unconditional stability and convergence of the schemes are discussed and theoretically proved. Then we extend the technique and develop a numerical scheme for the case of space fractional Laplacian operator. In addition, the recent L1+ formula is utilized to improving the temporal accuracy. Finally, three numerical examples in rectangular, elliptical and human brain-like irregular domains are given to demonstrate the accuracy and efficiency of the proposed numerical scheme.