|
#!/usr/local/bin/python |
|
# encoding: utf-8 |
|
""" |
|
ned-galaxy-semi-major-axis-vs-mag.py |
|
=============================== |
|
:Summary: |
|
Plot Galaxy Semi-Major Axis vs Galaxy Magnitude |
|
|
|
:Author: |
|
David Young |
|
|
|
:Date Created: |
|
February 14, 2017 |
|
|
|
Usage: |
|
ned-galaxy-semi-major-axis-vs-mag |
|
|
|
Options: |
|
-h, --help show this help message |
|
-v, --version show version |
|
-s, --settings the settings file |
|
""" |
|
## USER DEFINED VARIABLES - ONLY PLACE THIS SCRIPT NEEDS MODIFIED ## |
|
|
|
|
|
################# GLOBAL IMPORTS #################### |
|
import sys |
|
import os |
|
os.environ['TERM'] = 'vt100' |
|
import readline |
|
import math |
|
import re |
|
import glob |
|
import docopt |
|
from fundamentals import tools, times |
|
|
|
|
|
def main(arguments=None): |
|
""" |
|
The main function used when ``gsc-bright-star-associations.py`` is run as a single script from the cl, or when installed as a cl command |
|
""" |
|
|
|
# setup the command-line util settings |
|
su = tools( |
|
arguments=arguments, |
|
docString=__doc__, |
|
logLevel="DEBUG", |
|
options_first=False, |
|
projectName=False |
|
) |
|
arguments, settings, log, dbConn = su.setup() |
|
|
|
csvFile = os.path.dirname(__file__) + \ |
|
"/ned-galaxy-semi-major-diameter-vs-mag.csv" |
|
|
|
import codecs |
|
pathToReadFile = csvFile |
|
try: |
|
log.debug("attempting to open the file %s" % (pathToReadFile,)) |
|
readFile = codecs.open(pathToReadFile, encoding='utf-8', mode='r') |
|
thisData = readFile.read() |
|
readFile.close() |
|
except IOError, e: |
|
message = 'could not open the file %s' % (pathToReadFile,) |
|
log.critical(message) |
|
raise IOError(message) |
|
readFile.close() |
|
|
|
lines = thisData.split("\n")[1:] |
|
|
|
mags = [] |
|
axes = [] |
|
|
|
reMag = re.compile(r'\d+\.\d+') |
|
|
|
for l in lines: |
|
l = l.split(",") |
|
if len(l) < 2: |
|
continue |
|
match = reMag.match(l[1]) |
|
if match: |
|
mags.append(float(match.group(0))) |
|
axes.append(float(l[0]) * 30.) |
|
|
|
import numpy as np |
|
import matplotlib.pyplot as plt |
|
|
|
line = plt.figure() |
|
|
|
x = axes |
|
y = mags |
|
|
|
x = np.array(x) |
|
y = np.array(y) |
|
|
|
ax = plt.gca() |
|
ax.invert_yaxis() |
|
ax.set_xscale("log", nonposy='clip') |
|
|
|
ax.set_ylim(30, 10) |
|
ax.set_xlim(1, 1000) |
|
|
|
ax.set_title("NED Galaxy Mags vs Semi-Major Axis Radii") |
|
ax.set_xlabel("Semi Major Axis Radius (Arcsec)") |
|
ax.set_ylabel("Magnitude") |
|
|
|
plt.plot(x, y, "o", ms=1) |
|
|
|
xLine = [1, 10, 100, 1000] |
|
abline_values = [-6 * math.log10(i) + 25.3 for i in xLine] |
|
plt.plot(xLine, abline_values, 'b') |
|
|
|
# # draw vertical line from (70,100) to (70, 250) |
|
# plt.plot([70, 70], [100, 250], 'k-', lw=2) |
|
|
|
# # draw diagonal line from (70, 90) to (90, 200) |
|
# plt.plot([70, 90], [90, 200], 'k-') |
|
|
|
plt.show() |
|
return |
|
|
|
if __name__ == '__main__': |
|
main() |
