We study the homogenized energy densities of periodic ferromagnetic Ising systems. We prove that, for finite range interactions, the homogenized energy density, identifying the effective limit, is crystalline, i.e. its Wulff crystal is a polytope, for which we can (exponentially) bound the number of vertices. This is achieved by deriving a dual representation of the energy density through a finite cell formula. This formula also allows easy numerical computations: we show a few experiments where we compute periodic patterns which minimize the anisotropy of the surface tension.where (z) + denotes the positive part of z ∈ R, c i,j : L × L → [0, +∞) are T -periodic, that is c i+T z,j+T z = c i,j for all i, j ∈ L and z ∈ Z d and satisfy the following decay assumption (iii) (Decay of interactions) For all i ∈ L there holds j∈L c i,j |i − j| < +∞ .