Multiferroic magnetoelectric materials, which offer the possibility of manipulating the magnetic order by a low energy consuming electric field, have been attracting considerable attention during the last decade. However, materials that exhibit multiferroic properties above room temperature, and are thus suitable for future device applications, are very rare. First-principles calculations based on density functional theory have become of high relevance in exploring and predicting novel multiferroic materials. Likewise, the experimental validation of the magnetoelectric phenomena both at the atomic and macroscopic scales is crucial to improve our insight into the microscopic mechanisms underlying the predicted phenomena.Here, in this context, we study the effect of doping and epitaxial strain on the structure and the ferroelectric properties of the promising Aurivillius-phase compounds [1]. In particular, we investigate the potential multiferroic behavior of Bi 5 FeTi 3 O 15 as predicted by ab-initio calculations [2,3]. Its crystal structure consists of 4 perovskite layers stacked periodically along the [001] direction and separated by fluorite-like (Bi 2 O 2 ) 2+ layers. This phase has been predicted to have a large response of in-plane polarization under biaxial strain, and a preferential occupation of the inner (outer) perovskite site with the Fe 3+ for tensile (compressive) strain.The experimental characterization of this system is very challenging, dealing with the requirement of very high spatial resolution (< 1Å). In fact, while the detection of the Fe substitutional species in preferential lattice sites requires analytical investigations with atomic resolution, the estimation of the local polarization at the atomic scale is nevertheless more demanding, being the typical atomic displacements involved in the polar distortions in the 10-50 pm range [4,5]. To perform the abovementioned investigations at sub-Ångstrom resolution we use probe-corrected scanning transmission electron microscopy (STEM). Specifically, we exploit high-angle annular dark-field (HAADF) and annular bright-field (ABF) images to get access to both heavy and light atomic columns.In order to investigate the microscopic properties of strained and doped Aurivillius phases we have implemented an analytical algorithm able to retrieve the polarization at the atomic scale from the structural distortions, i.e. from atomic or charge's center of mass displacements. In addition, other distortions like tilting of the oxygens octahedrons, tetragonal deformation and metal-oxygen-metal bonding angle bending, that play a pivotal role in determining the ferroelectric/ferromagnetic ordering are measured by our algorithm. Combining the HAADF and ABF signals, the accurate positions of all the different atomic species in the structure are precisely measured fitting the peaks of the images with 2D asymmetric Gaussians [6]. The polarization at the atomic scale can therefore be computed as the displacement of the B cation with respect to the center of mass of th...