The aim of this work is to provide a characterization for the Besov spaces \(B^{r,q}_p(\mathbb {R}^n)\) via the continuous shearlet transform in higher dimensions with isotropic dilations to show the regularity of weak solutions \(u \in W^{m,p}(\mathbb {R}^n)\) under the equation \(Qu =f\) , so that for a given \(f \in B^{r,q}_p(\mathbb {R}^n)\) , \(u \in B^{r+m,q}_p(\mathbb {R}^n)\) . In this case, \(1< p,q < \infty \) , \(0< r < 1\) , and where \(Q = \sum _{\beta \le m}c_{\beta } \partial ^{\beta }\) is a partial differential operator of pure order m with positive constant coefficients \(c_{\beta }\) .