The ground state of macroscopic samples of magnetically ordered materials is a domain state because of magnetostatic energy or entropy, yet we have limited experimental means for imaging the bulk domain structure and the magnetization process directly. The common methods available reveal the domains at the surface or in electron- or x-ray transparent lamellae, not those in the bulk. The magnetization curve just reflects the vector sum of the moments of all the domains in the sample, but magnetostriction curves are more informative. They are strongly influenced by the domain structure in the unmagnetized state and its evolution during the magnetization process in an applied field. Here we report a method of determining the bulk domain struct...