예제 #1
0
def plotGrid(plt):
    plt.figure(figsize=(18, 11.25))
    # plotUtil.setGridSize(-2, 33, 22, 52)
    # plotUtil.setGridSize(9, 16, 36, 43)
    plotUtil.setGridSize(9.3, 12, 40.2, 42.4)
    plt.xlabel('RA')
    plt.ylabel('DEC')
    plotUtil.plotGrid(plt, 0.1)
예제 #2
0
def plotGalaxy(plt):
    plt.plot(common.M31_RA, common.M31_DEC, '*', color='r')
    plt.plot(common.M32_RA, common.M32_DEC, '*', color='r')
    plt.plot(common.M33_RA, common.M33_DEC, '*', color='r')

    plotUtil.plotEllipseWithParam(plt, common.M31_RA, common.M31_DEC, common.M31_Width, common.M31_height,
                                  common.M31_Angle, (0, 0, 0, 0), (0, 0, 0, 0.4))
    plotUtil.plotEllipseWithParam(plt, common.M32_RA, common.M32_DEC, common.M32_Width, common.M32_height,
                                  common.M31_Angle, (0, 0, 0, 0), (0, 0, 0, 0.4))
예제 #3
0
def plotOver(lamostTable, sdssTable, mu):
    dic_lamost = getCoordinateDic(lamostTable)
    dic_sdss = getCoordinateDic(sdssTable)

    dic_over = {}

    for ra in dic_lamost.keys():
        for dec in dic_lamost[ra].keys():
            lcount = dic_lamost[ra][dec]
            scount = util.getDic2Value(dic_sdss, ra, dec)

            if scount is None:
                util.setItemToDic2(dic_over, ra, dec, 10)
            else:
                util.setItemToDic2(dic_over, ra, dec, float(lcount) / float(scount))

    for ra in dic_sdss.keys():
        for dec in dic_sdss[ra].keys():
            if not util.isInDic2(dic_lamost, ra, dec):
                util.setItemToDic2(dic_over, ra, dec, 11)

    for ra in dic_over.keys():
        for dec in dic_over[ra].keys():
            alpha = dic_over[ra][dec]

            if alpha == 10:
                plotUtil.plotRectangleWithParam(plt, ra, dec, 1, 1, "red")
            elif alpha == 11:
                plotUtil.plotRectangleWithParam(plt, ra, dec, 1, 1, "blue")

            else:
                newAlpha = alpha * mu
                if newAlpha > 1:
                    newAlpha = 1

                plotUtil.plotRectangleWithParam(plt, ra, dec, 1, 1, "black", newAlpha)

    plt.title("Lamost/SDSS Covered")
    CommonPlot.plotGrid(plt)
    CommonPlot.plotGalaxy(plt)

    plt.savefig('/media/zdwdong/my/fits/Cover/Lamost_OVER_SDSS_' + str(mu) + '.png')
    plt.close()
예제 #4
0
x_mean = []
y_mean = []
x_3std1 = []
y_3std1 = []
x_3std2 = []
y_3std2 = []

for key in dic.keys():
    value = dic[key]
    std = np.std(value)
    value_new = value / std

    dic[key] = value_new
    mean = np.mean(value_new)

    plotUtil.plotText(plt, key, 8.5, round(std, 2))
    plotUtil.plotText(plt, key, 9, round(mean, 1))

    count_value = len(value)
    x = []
    for i in range(count_value):
        x.append(key)
    plt.plot(x, value_new, '.')

    x_mean.append(key)
    y_mean.append(mean)

    x_3std1.append(key)
    x_3std2.append(key)

    y_3std1.append(mean + 3)
예제 #5
0
# coding=utf-8
__author__ = 'zdwdong'

"""
画 teff/logg散点图
"""

import CodeLib.Database.MySQLHepler as sqlHelper
import matplotlib.pyplot as plt
import CodeLib.Util.PlotUtil as plotUtil

sql_select = "select teff,logg from dr3.M31_ELL where teff>0 and logg>0"

data = sqlHelper.getFetchAll(sql_select)

teff = []
logg = []

for dataRow in data:
    teff.append(float(dataRow[0]))
    logg.append(float(dataRow[1]))

plt.title("Lamost Logg_teff")
plt.plot(teff, logg, '.')

plt.xlabel('teff')
plt.ylabel('logg')
plotUtil.plotGrid(plt)
plotUtil.invertX(plt)
plt.show()