Dynamics are inherent characteristics of batch processes, which may be not only within a batch, but also from batch to batch. Two-dimensional dynamic principal component analysis (2-D-DPCA) method [Lu, N., Yao, Y., Gao, F., 2005. Two-dimensional dynamic PCA for batch process monitoring. A.I.Ch.E. Journal 51, 3300–3304] can model both kinds of batch dynamics, but may lead to the inclusion of large number of lagged variables and make the contribution plot difficult to read. To solve this problem, subspace identification technique is combined with 2-D-DPCA in this paper. The state space model of a 2-D batch process can be identified with canonical variate analysis (CVA) method based on the auto-determined support region (ROS). In 2-D-DPCA modeling, the utilization of state variables instead of lagged process variables reduces the number of variables and provides a clearer contribution plot for fault diagnosis.