Let $A \subset [0,1]^n$ with $A$ measurable and such that $\mathcal{L}^n (A)= \delta >0$, and consider a partition of $[0,1]^n$ in $\epsilon$-cubes (i.e. cubes of side $\epsilon)$.
For $\epsilon \to 0$ (say, $\epsilon= 1/n$ with $n \in \mathbb{N}$ so that $1/\epsilon$ is always an integer) we should have that, at the limit, it holds the Lebesgue differentiation theorem, so you have that points in $A$ have density 1 and those in $A^c$ have density 0. What I mean is that given every (in fact $\mathcal{L}^n$ almost every, but whatever) point $x$ in $[0,1]^n$ I think that if you take the sequence of the $\epsilon$ -cubes $Q_\epsilon (x)$ such that they contain $x$, one should have that
$$
\frac{| A \cap Q_\epsilon (x)|}{\epsilon ^n} \to I_{A} (x),
$$ with $I_A$ the indicator function of $A$.
In fact the Lebesgue theorem is about shrinking balls centered at the point $x$ you consider, but at the limit I think this shouldn't matter if you consider finer and finer partitions.
What I would like to know is if there is a sort of quantitative version of this result for fixed $\epsilon$. What I mean is, given a partition $\{Q_\epsilon (i)\}_{i=1, \dots, \epsilon^{-n}}$ , having some result which gives you a sort of balance between how big $f_{A, \epsilon } (i) =\frac{| A \cap Q_\epsilon (i)|}{\epsilon ^n}$ is and in how many cubes you have this.
So for example, if I fix $\epsilon$ we can't have $ f_{A, \epsilon} (i)= 1$ for more than $\delta/ \epsilon ^n$ indices, but it's easy to give an example of an $A$ such that $f_{A, \epsilon} (i)= \delta$ for every $i$.
I was thinking that since the convergence from the Lebesgue theorem is pointwise I can play a bit with Egoroff's theorem or something like that, but maybe these kind of estimates are already considered and have a name in the literature (although I couldn't find anything).