sc = read_rds("4_harmony_with_covidcontrol_res1.0/sc.NormByLocationRep.Harmony.final.rds")
## Registered S3 method overwritten by 'spatstat.geom':
##   method     from
##   print.boxx cli

Summary

dim(sc)
## [1]  50780 173243