We discuss the magnetization M_m in the m-th column of the zig-zag layered 2D Ising model on a half-plane using Kadanoff-Ceva fermions and orthogonal polynomials techniques. Our main result gives an explicit representation of M_m via m×m Hankel determinants constructed from the spectral measure of a certain Jacobi matrix which encodes the interaction parameters between the columns. We also illustrate our approach by giving short proofs of the classical Kaufman-Onsager-Yang and McCoy-Wu theorems in the homogeneous setup and expressing M_m as a Toeplitz+Hankel determinant for the homogeneous sub-critical model in presence of a boundary magnetic field.