Metastable states of homogeneous 2D and 3D Ising models are studied under free boundary conditions. The states are defined in terms of weak and strict local minima of the total interaction energy. The morphology of these minima is characterized locally and globally on square and cubic grids. Furthermore, in the 2D case, transition from any spin configuration that is not a strict minimum to a strict minimum is possible via non-energy-increasing single flips.