def getCoordinateCountMap(data):
"""
get dictionary of ra and dec
ra - dec - count
"""
dic_radec = {}
for row in data:
ra = int(row[0])
dec = int(row[1])
count = int(row[2])
util.setItemToDic2(dic_radec, ra, dec, count)
return dic_radec
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()
