Lusin approximation for functions of bounded variation

We prove a Lusin approximation of functions of bounded variation. If $f$ is a function of bounded variation on an open set $\mathrm{\Omega }\subset X$, where $X=(X,d,\mu )$ is a given complete doubling metric measure space supporting a $1$-Poincaré inequality, then for every $\epsilon >0$, there exist a function $f_{\epsilon }$ on $\mathrm{\Omega }$ and an open set $U_{\epsilon }\subset \mathrm{\Omega }$ such that the following properties hold true: \begin{enumerate} \item ${\mathrm{C}\mathrm{a}\mathrm{p}}_1(U_{\epsilon })<\epsilon$; \item $\mathrm{\parallel }f-f_{\epsilon }{\mathrm{\parallel }}_{\text{BV}(\mathrm{\Omega })}<\epsilon$; \item $f^{\vee }\equiv f_{\epsilon }^{\vee }$ and $f^{\wedge }\equiv f_{\epsilon }^{\wedge }$ on $\mathrm{\Omega }\setminus U_{\epsilon }$; \item $f_{\epsilon }^{\vee }$ is upper semicontinuous on $\mathrm{\Omega }$, and $f_{\epsilon }^{\wedge }$ is lower semicontinuous on $\mathrm{\Omega }$. \end{enumerate} If the space $X$ is unbounded, then such an approximating function $f_{\epsilon }$ can be constructed with the additional property that the uniform limit at infinity of both $f_{\epsilon }^{\vee }$ and $f_{\epsilon }^{\wedge }$ is $0$. Moreover, when $X={\mathbb{R}}^d$, we show that the non-centered maximal function of $f_{\epsilon }$ is continuous in $\mathrm{\Omega }$.