def aggregate_and_plot(self): global SIC_ORIG global SIC_ROOT global SIC_PROCESSED global SIC_RESULTS global SIC_DATA_PICKLE mskpath = join(SIC_ROOT, SIC_PROCESSED) global d try: d except NameError: d = pickle.load(file(join(SIC_ROOT, SIC_PROCESSED, SIC_DATA_PICKLE))) o2n = d["o2n"] spots = aggregate_spots(o2n, mskpath) d["spots"] = spots pickle.dump(d, file(join(SIC_ROOT, SIC_RESULTS, SIC_DATA_PICKLE), "w")) histogram_intensities(spots, mskpath) scatterplot_intensities(spots, mskpath) spots_per_cell_distribution(spots, mskpath) pl.show()
def aggregate_and_plot(self): global SIC_ORIG global SIC_ROOT global SIC_PROCESSED global SIC_RESULTS global SIC_DATA_PICKLE mskpath = join(SIC_ROOT, SIC_PROCESSED) global d try: d except NameError: d = pickle.load( file(join(SIC_ROOT, SIC_PROCESSED, SIC_DATA_PICKLE))) o2n = d["o2n"] spots = aggregate_spots(o2n, mskpath) d["spots"] = spots pickle.dump(d, file(join(SIC_ROOT, SIC_RESULTS, SIC_DATA_PICKLE), "w")) histogram_intensities(spots, mskpath) scatterplot_intensities(spots, mskpath) spots_per_cell_distribution(spots, mskpath) pl.show()
def make_plots(spots, d): pf.histogram_intensities(spots) pf.scatterplot_intensities(spots) pf.spots_per_cell_distribution(spots) #pf.plot_time2ratio_between_one_dot_number_and_cell_number(d) pl.show()