wflynny

import re
from pathlib import Path
from xml.etree import ElementTree as ET

import numpy as np
from skimage import exposure
from skimage.io import imread, imsave, imshow
from skimage.filters import threshold_otsu, gaussian

wflynny

from matplotlib.text import Annotation
from matplotlib.legend_handler import HandlerBase
from matplotlib import patheffects
from matplotlib.colors import CSS4_COLORS, hex2color

class AnnotationHandler(HandlerBase):
    def create_artists(self, legend, artist, xdescent, ydescent,
                       width, height, fontsize, trans):
        a = Annotation(
            artist.get_text(), [width/2, height/2],

wflynny

###################
# chord diagram
import matplotlib.pyplot as plt
from matplotlib.path import Path as mplPath
import matplotlib.patches as patches

import numpy as np
import pandas as pd

LW = 0.3

wflynny

Introduction

This process has gotten simpler with 10X Genomics Loupe Browser 5.0.0+. I try to keep a cache of all 10X downloads internal to JAX here (must be on the VPN):

http://singlecell-software.s3-far.jax.org/list.html

where you can download the latest version of Loupe Browser for your OS. Otherwise, head over to the 10X site:

https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest#loupe

wflynny

import scanpy as sc
import pandas as pd
import GEOparse

def gse_to_adata(gse):
    X = gse.pivot_samples("VALUE")
    key = list(gse.gpls.keys())[0]
    annotated = X.index.to_frame(index=False).merge(
        gse.gpls[key].table[["ID", "Gene Symbol"]], 
        left_on="ID_REF", right_on="ID"