In [1]:
%pylab inline

from __future__ import (division, print_function)

import os
import sys
import copy
import fnmatch
import warnings
import collections

import numpy as np
import scipy
try:
    from scipy.stats import scoreatpercentile
except:
    scoreatpercentile = False
from scipy.interpolate import interp1d
import cPickle as pickle

# Astropy
from astropy.io import fits
from astropy    import units as u
from astropy.stats import sigma_clip
from astropy.table import Table, Column
from astropy.utils.console import ProgressBar

# AstroML
from astroML.plotting import hist

# Matplotlib related
import matplotlib.pyplot as plt
from matplotlib.patches import Ellipse
from matplotlib.ticker import NullFormatter
from matplotlib.ticker import MaxNLocator
# Matplotlib default settings
rcdef = plt.rcParams.copy()
pylab.rcParams['figure.figsize'] = 12, 10
pylab.rcParams['xtick.major.size'] = 8.0
pylab.rcParams['xtick.major.width'] = 2.5
pylab.rcParams['xtick.minor.size'] = 4.0
pylab.rcParams['xtick.minor.width'] = 2.5
pylab.rcParams['ytick.major.size'] = 8.0
pylab.rcParams['ytick.major.width'] = 2.5
pylab.rcParams['ytick.minor.size'] = 4.0
pylab.rcParams['ytick.minor.width'] = 2.5

# Personal
import hscUtils as hUtil
import galSBP


Populating the interactive namespace from numpy and matplotlib

In [2]:
def readProfile(ellFile):
    """ Load the pickle format 1-D profile. """
    if os.path.isfile(ellFile):
        return pickle.load(open(ellFile, 'rb'))
    else:
        warnings.warn("!!! Can not find the Ellipse Output at %s" % 
                      ellFile )
        return None

In [3]:
location = '/Users/songhuang/work/hscs/hsc_cutout/test/127/HSC-I/default'
testFile = 'redBCG_127_HSC-I_full_psf_ellip_3.pkl'
testEll = os.path.join(location, testFile)
testSbp = readProfile

In [6]:
np.max(testSbp['sma'])


Out[6]:
21.435890000000001

In [42]:
rad = (testSbp['sma'] ** 0.25)
sbp = - 2.5 * np.log10(testSbp['intens'])


/usr/local/lib/python2.7/site-packages/ipykernel/__main__.py:2: RuntimeWarning: invalid value encountered in log10
  from ipykernel import kernelapp as app

In [49]:
rsma_common = np.arange(0.4, 4.1, 0.01)
intrpFunc = interp1d(rad, sbp, kind='slinear', bounds_error=False)

In [50]:
sbpInterp = intrpFunc(rsma_common)

In [51]:
plt.plot(rad, sbp, c='r')
plt.plot(rsma_common, sbpInterp, c='b')


Out[51]:
[<matplotlib.lines.Line2D at 0x10f296050>]

In [52]:
sbpInterp


Out[52]:
array([  3.07062774,   3.07266096,   3.07469419,   3.07672742,
         3.07876064,   3.08079387,   3.0828271 ,   3.08486033,
         3.08689355,   3.08892678,   3.09096001,   3.09299323,
         3.09502646,   3.09705969,   3.09909291,   3.10112614,
         3.10315937,   3.1051926 ,   3.10722582,   3.10925905,
         3.11129228,   3.1133255 ,   3.11535873,   3.11739196,
         3.11942518,   3.12145841,   3.12349164,   3.12552487,
         3.12755809,   3.12959132,   3.13162455,   3.13365777,
         3.135691  ,   3.13772423,   3.13975745,   3.14179068,
         3.14382391,   3.14585714,   3.14789036,   3.14992359,
         3.15195682,   3.15399004,   3.15602327,   3.1580565 ,
         3.16008972,   3.1624086 ,   3.17131779,   3.18022699,
         3.18989827,   3.19957501,   3.21000985,   3.22053381,
         3.23175771,   3.24321973,   3.25525078,   3.26775557,
         3.28060032,   3.2942687 ,   3.30793709,   3.32303859,
         3.33816295,   3.35572613,   3.37477907,   3.39498237,
         3.41977951,   3.44457666,   3.47290784,   3.50260394,
         3.53323831,   3.56806919,   3.60290007,   3.64127963,
         3.68219885,   3.72311808,   3.771716  ,   3.82097901,
         3.87216889,   3.93138885,   3.99060881,   4.05270595,
         4.11910102,   4.18549609,   4.24989237,   4.3125938 ,
         4.37529523,   4.44020973,   4.50652538,   4.57284103,
         4.65919457,   4.7582058 ,   4.85721703,   4.96182706,
         5.07118554,   5.18054401,   5.28513128,   5.38241835,
         5.47970541,   5.57901739,   5.68797266,   5.79692792,
         5.90588318,   6.02595185,   6.1476205 ,   6.26928915,
         6.38707935,   6.50091589,   6.61475243,   6.7288877 ,
         6.85131481,   6.97374192,   7.09616903,   7.21423767,
         7.32778614,   7.44133461,   7.55488308,   7.66797746,
         7.78099889,   7.89402032,   8.00606877,   8.11236687,
         8.21866497,   8.32496308,   8.43167439,   8.53919415,
         8.64671392,   8.75423368,   8.8607974 ,   8.96614425,
         9.0714911 ,   9.17683794,   9.27909173,   9.37753003,
         9.47596833,   9.57440663,   9.66798629,   9.75289757,
         9.83780884,   9.92272012,  10.00756452,  10.09208713,
        10.17660974,  10.26113235,  10.34565495,  10.45142957,
        10.5599937 ,  10.66855783,  10.77712195,  10.86768155,
        10.93968692,  11.0116923 ,  11.08369767,  11.15570305,
        11.22215323,  11.28857604,  11.35499886,  11.42142168,
        11.51213872,  11.64097206,  11.7698054 ,  11.89863874,
        12.02747207,  12.11858429,  12.19131577,  12.26404726,
        12.33677874,  12.40951022,  12.73110843,  13.09950083,
        13.46789322,  13.83628562,  14.20467802,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan,          nan,          nan,
                nan,          nan])

In [ ]: