  Most of the literature in the \emph{bulk reconstruction program} in
holography focuses on recovering local bulk operators propagating on a
quasilocal bulk geometry and the knowledge of the bulk geometry is always
assumed or guessed. The fundamental problem of the bulk reconstruction program,
which is \emph{recovering the bulk background geometry (metric)} from the
boundary CFT state is still outstanding. In this work, we formulate a recipe to
extract the bulk metric itself from the boundary state, specifically, the
modular Hamiltonian information of spherical subregions in the boundary. Our
recipe exploits the recent construction of Kabat and Lifschytz
\cite{Kabat:2017mun} to first compute the bulk two point function of scalar
fields directly in the CFT without knowledge of the bulk metric or the
equations of motion, and then to take a large scaling dimension limit (WKB) to
extract the geodesic distance between two close points in the bulk i.e. the
metric. As a proof of principle, we consider three dimensional bulk and
selected CFT states such as the vacuum and the thermofield double states. We
show that they indeed reproduce the pure AdS and the regions outside the
Rindler wedge and the BTZ black hole \emph{up to a rigid conformal factor}.
Since our approach does not rely on symmetry properties of the CFT state, it
can be applied to reconstruct asymptotically AdS geometries dual to arbitrary
general CFT states provided the modular Hamiltonian is available. We discuss
several obvious extensions to the case of higher spacetime dimensions as well
as some future applications, in particular, for constructing metric beyond the
causal wedge of a boundary region. In the process, we also extend the
construction of \cite{Kabat:2017mun} to incorporate the first order
perturbative locality for AdS scalars.
