A novel numerical transfer matrix approach to obtain the critical behavior of Ising models in dimensions higher than two is described. The approach augments the physical transfer matrices with other matrices which do not necessarily have a straightforward physical interpretation. These additional matrices have approximately the same form as the physical transfer matrices, and for high dimensions may provide the only way of performing finite size scaling using transfer matrix calculations. Results are presented for the Ising model in three, four, and five dimensions.