We survey recent results on inverse problems for geodesic X-ray transforms and other linear and non-linear geometric inverse problems for Riemannian metrics, connections and Higgs fields defined on manifolds with boundary.Johann Radon and his contemporaries formulated several integral geometric problems, not only in linear but also in non-linear settings [42,123]. Such problems, namely travel-time tomography and boundary rigidity as later formulated in [72,62], are concerned with recovering a Riemannian metric from the shortest length between any two boundary points. Such problems and their cousins (described below), now make the field of integral geometry, or how to reconstruct geometric features of a manifold from integral functionals defined over that manifold.Nowadays this field forms the basis of several non-invasive approaches to imaging internal properties of materials: seismology [42,123], or how to reconstruct the density inside the Earth from first arrival times of seismic wavefronts; medical imaging since the development of X-ray Computerized Tomography [75,119,25]; Single-Photon Emission Computerized Tomography using the attenuated X-ray transform [76,74,78]; vector tomography in helio-seismology [51,94,52]; ocean imaging [73]; X-ray diffraction strain tomography [59,19] and tomography in elastic media [103, Ch. 7][108]; neutron imaging, as applied to the imaging of vertebrate remains [102] and shales [11]. Non-linear integral geometric problems also continue to find new applications: recently, Neutron Spin Tomography [99] as a means to measure magnetic fields in materials, has arised as a novel method which can be of use in electrical engineering, superconductivity, etc. The transform to invert in this case is a non-linear operator, the so-called "non-abelian X-ray transform" of the magnetic field, see Problem 3 below.Recent breakthroughs have fuelled the field, exploiting a combination of old and new methods. Examples of such methods are: the systematic use of analysis on the unit sphere bundle combining energy methods (also coined "Pestov identities"), initiated by Mukhometov [71] and generalized in [91,103], and harmonic analysis on the tangent fibers [15,83,86]; in dimensions three and higher, the discovery in [120] that the existence of a foliation of the domain by strictly convex hyperfsurfaces, local or global, yields a powerful and robust approach to integral geometric inversions [116,118,125, 127], via a successful use of Melrose's scattering calculus [60]; the systematic use of analytic microlocal analysis to produce 'generic' results, implying the unique identifiability of unknown parameters in an open and dense subset of all cases [112,45, 128]; finally, recent results in the context of Anosov flows, leading to positive results for certain geometries with trapped sets [33,34,38].This review article aims at giving an overview of the arsenal of these methods, and to describe to what extent they help coping with various geometric settings, whose complexity is mainly governed by two features ...