In the book of Stein, on page 9:
"Let us remark that these additional properties easily lead to the following conclusions among others. First note that $\mu (B) >0$ for any ball $B$, which is a consequence of the doubling condition, (iii), and the fact that $\mu( \mathbb{R}^n) >0$. It then follows from (iv) that for any locally integrable $f$ and any $\delta>0$, the mean value:
$$A_\delta(f(x)) = \frac{1}{\mu(B(x,\delta))}\int_{B(x,\delta)} f(y) d\mu(y)$$
is continuous as a function of $x$.
"
The second conclusion of $A_\delta(f(X))$ is what I don't know how to conclude. Do you have any hints?
Here are the conditions listed:
(i) $B(x,\delta)\cap B(y,\delta) \ne \emptyset $ implies $B(y,\delta)\subset B(x, c_1 \delta)$.
(ii) $\mu(B(x,c_1\delta))\le c_2 \mu(B(x,\delta))$ where $c_1,c_2>1$ are constants.(this is the doubling condition).
(iii) $\bigcap_\delta \overline{B}(x,\delta) = \{ x\} $, $\bigcup_\delta B(x,\delta)=\mathbb{R}^n$.
(iv) For each open set $U$ and each $\delta>0$ the function $x\mapsto \mu(\{ B(x,\delta)\cap U\} ) $ is continuous.
As you pointed out to Felipeh, it's important to remember that the bounded open (with respect to Euclidean topology on $\mathbb{R}^{n}$) sets $B(x,\delta)$ are a priori "balls", not actual balls define by a metric or some weaker quasidistance function. Perhaps, you have already shown this, but I don't think it is that obvious that $\mu$ is finite on the balls $B(x,\delta)$ and inner regular, which Stein uses in the proof of the weak type $(1,1)$ inequality for the associated maximal operator. One can show this, but I thought it was a nontrivial exercise.
For the benefit of those users who don't have access to the book, here's a relevant passage on pg. 8 at the beginning of the section describing what these sets are:
The OP has the hypotheses in the question.
Regarding your question, note that condition (iv) implies that
Now fix $x\in\mathbb{R}^{n}$ and $\delta>0$, and suppose that
$$\mu(B(x,\delta)\triangle B(y,\delta))\not\rightarrow 0, \quad y\rightarrow x$$ Then there exists some $c>0$ and a sequence $y_{k}\rightarrow y$, such that
$$\mu(B(x,\delta)\setminus B(y_{k},\delta))+\mu(B(y_{k},\delta)\setminus B(x,\delta))=\mu(B(x,\delta)\triangle B(y_{k},\delta))\geq c, \quad\forall k$$
Given $0<\epsilon<\min\left\{c/2,\mu(B(x,\delta))/2\right\}$, there exists an open neighborhood $U$ about $x$, such that $$y\in U\Longrightarrow \left|\mu(B(x,\delta)-\mu(B(y,\delta))\right|<\epsilon$$ Since the function $y\mapsto \mu(B(x,\delta)\cap B(y,\delta))$ is also continuous, replacing $U$ by a smaller open neighborhood of $x$, we have that $$y\in U\Longrightarrow \mu(B(x,\delta)\setminus B(y,\delta))=\mu(B(x,\delta))-\mu(B(x,\delta)\cap B(y,\delta))<\epsilon,$$ Without loss of generality, we may assume that $y_{k}\in U$ for all $k$, therefore $$\mu(B(y_{k},\delta)\setminus B(x,\delta))\geq c/2,\quad\forall k$$ and therefore \begin{align*} \mu(B(x,\delta))&=\lim_{k\rightarrow\infty}\mu(B(y_{k},\delta))=\mu(B(y_{k},\delta)\setminus B(x,\delta))+\mu(B(y_{k},\delta)\cap B(x,\delta))\\ &\geq c/2+\lim_{k\rightarrow\infty}\mu(B(y_{k},\delta)\cap B(x,\delta))\\ &=c/2+\mu(B(x,\delta)), \end{align*} which is a contradiction. Summarizing,
We'll need the following lemma.
Proof. By local integrability, the function $f\chi_{K}\in L_{\mu}^{1}(\mathbb{R}^{n})$. Let $g\in L_{\mu}^{1}(\mathbb{R}^{n})$ be a bounded function such that $\left\|f-g\right\|_{1}<\epsilon$. Then \begin{align*} \int_{E}\left|f(y)\right|d\mu(y)&=\int_{E}\left|[f\chi_{K}(y)-g(y)]+g(y)\right|d\mu(y)\\ &\leq\left\|f\chi_{K}-g\right\|_{L_{\mu}^{1}}+\int_{E}\left|g(y)\right|d\mu(y)\\ &\leq\epsilon/2+\left\|g\right\|_{L_{\mu}^{\infty}}\left|E\right| \end{align*} Taking $\gamma<\epsilon/(2\left\|g\right\|_{\infty})$ completes the proof. $\Box$
In order to apply the lemma, we need to show that for all $y$ closed to $x$, the union $B(x,\delta)\cup B(y,\delta)$ are contained in some bounded measurable set, which may depend on $x$. But this is trivial, as $U=B(x,\delta)$ is an open neighborhood about $x$, and for all $y\in U$, $B(x,\delta)\cap B(y,\delta)\neq\emptyset$, whence by condition (i), $$B(x,\delta)\cup B(y,\delta)\subset B(x,c_{1}\delta),$$ which is a bounded, measurable set. As shown above, there exists an open neighborhood $U$ about $x$ such that $$y\in U\Rightarrow\mu(B(x,\delta)\triangle B(y,\delta))<\gamma\Rightarrow\int_{B(x,\delta)\triangle B(y,\delta)}\left|f(z)\right|d\mu(z)<\epsilon,$$ given $\epsilon>0$. It follows immediately that $$\lim_{y\rightarrow x}\int_{B(y,\delta)}\left|f(z)\right|d\mu(z)=\int_{B(x,\delta)}\left|f(z)\right|d\mu(z)$$
Putting the above result together with $\mu(B(y,\delta))\rightarrow\mu(B(x,\delta))$, as $y\rightarrow x$, completes the proof.