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).