Here is a simple jupyter notebook scriptto find out how to idemtify spacially expressive genes in glioblastoma data.
The dataset is freely available from 10X Genomics.
Packages used : scanpy (for pre-processing) , leidenalg (for clustering), SpatialDE (for identification of spatially varied genes).