def bestnet(goodnet):
goodnet.write(write('data.sdss.org'))
goodnet.write(write('api.sdss.org'))
config._check_access()
config.access = 'collab'
config.setRelease('MPL-6')
yield goodnet
def setUp(self):
super(TestGalaxyPage, self).setUp()
self.blue = 'galaxy_page'
config.setRelease('MPL-5')
self.mode = config.mode
self.release = config.release
self.params = {'release': self.release}
# set up cube and expected values
self.cube = Cube(plateifu=self.plateifu, mode=self.mode)
# NSA params for 8485-1901
self.cols = [
'z', 'elpetro_logmass', 'sersic_n', 'elpetro_absmag_i',
'elpetro_absmag_g_r', 'elpetro_th50_r', 'elpetro_absmag_u_r',
'elpetro_absmag_i_z', 'elpetro_ba', 'elpetro_phi',
'elpetro_mtol_i', 'elpetro_th90_r'
]
self.exp_nsa_plotcols = {
'z': 0.0407447,
'elpetro_logmass': 9.565475912843823,
'sersic_n': 3.29617,
'elpetro_absmag_i': -19.1125469207764,
'elpetro_absmag_g_r': 0.5961246490479013,
'elpetro_th50_r': 1.33067,
'elpetro_absmag_u_r': 1.7617149353028019,
'elpetro_absmag_i_z': 0.20068740844720168,
'elpetro_ba': 0.87454,
'elpetro_phi': 154.873,
'elpetro_mtol_i': 1.30610692501068,
'elpetro_th90_r': 3.6882
}
def set_the_config(release):
"""Set config release without parametrizing.
Using ``set_release`` combined with ``galaxy`` double parametrizes!"""
config.access = 'collab'
config.setRelease(release)
set_sasurl(loc='local')
config.login()
config._traceback = None
def test_tree(self, bestnet, release):
assert config.access == 'collab'
assert 'mangawork' in os.environ['MANGA_SPECTRO_REDUX']
assert 'MPL' in config.release
config.setRelease(release)
if 'mpl' in release:
assert 'mangawork' in os.environ['MANGA_SPECTRO_REDUX']
else:
assert release in os.environ['MANGA_SPECTRO_REDUX']
assert config.release == release.upper()
def _setrun_query(self, limit=10):
config.setRelease("MPL-4")
self.filter = 'nsa.z < 0.1'
self._set_remote(limit=limit)
r = self.q.run()
r.sort('z')
self.assertEqual(1213, r.totalcount)
self.assertEqual(limit, r.count)
self.assertEqual(limit, len(r.results))
self.assertEqual(limit, r.limit)
self.assertEqual(limit, r.chunk)
return r
def check_request_for_release(request):
''' Check Flask request for release information and set it in config
Parameters:
request (Flask.Request):
Flask Request object
'''
# get release from request header
release = request.headers.get('Release', None)
# if no release, get release variable from Flask request
if not release:
release = request.environ.get('MARVIN_RELEASE', None)
# if release then set change it in the config but only when the session is empty
if release:
config.setRelease(release.upper())
def test_sasurl(self, bestnet, release):
assert 'api.sdss.org/marvin' in config.sasurl
config.setRelease(release)
sasurl = 'dr15.sdss.org/api/marvin' if 'DR' in release else 'api.sdss.org/marvin'
assert sasurl in config.sasurl
def test_invalid_release(self):
with pytest.raises(MarvinError) as cm:
config.setRelease('MPL-1')
assert 'trying to set an invalid release version.' in str(cm.value)
def test_drpall(self, bestnet, release):
assert 'mangawork' in config.drpall
config.setRelease(release)
if config.drpall:
word = 'mangawork' if 'mpl' in release else release
assert word in config.drpall
def vacs():
config.setRelease('DR15')
v = VACs()
yield v
v = None
def test_sasurl(self, bestnet, release):
assert 'sas.sdss.org' in config.sasurl
config.setRelease(release)
sasurl = 'dr15.sdss.org' if 'DR' in release else 'sas.sdss.org'
assert sasurl in config.sasurl
#!/usr/bin/env python
import numpy as np
import matplotlib.pyplot as plt
from marvin.tools.maps import Maps
from marvin.utils.general.general import downloadList
from marvin import config
from adtools import *
import time
if __name__ == '__main__':
start = time.time()
config.setRelease('MPL-5')
spx = readall()
#spx = loadfits('SPX-GAU-MILESHC-coposite_1.00Re.fits')
nbins = 5
stepsize = .25
bins = np.arange(0, stepsize * (nbins + 1), stepsize)
alldata = []
#for i in range(len(spx)):
for i in [3]:
print(len(spx[i]))
plateifu = []
d4000 = np.zeros(len(spx[i]))
res = np.zeros_like(d4000)
for j in range(len(spx[i])):
plateifu += [
str(spx[i]['plate'][j]) + '-' + str(spx[i]['ifudesign'][j])
]
def prepare_marvin():
from marvin import config
config.access = "collab"
config.mode = "auto"
config.setRelease("MPL-10")
config.login()
def set_the_config(release):
"""Set config release without parametrizing.
Using ``set_release`` combined with ``galaxy`` double parametrizes!"""
config.setRelease(release)
config._traceback = None
import matplotlib.pyplot as plt
import numpy as np
import marvin
from marvin import config, marvindb
import math
from bresenham import bresenham
import numpy.ma as ma
from scipy.optimize import curve_fit
from scipy import asarray as ar, exp
import cv2
get_ipython().run_line_magic('matplotlib', 'inline')
# Make sure you have collaborator access - if not go set it up following these instructinos:
# https://sdss-marvin.readthedocs.io/en/stable/installation.html
config.access = 'collab'
# Choose the data release you would like to use (could also use MPL)
config.setRelease('MPL-8')
print(config.access)
# I had to re-log in when I first ran this code:
config.login(refresh=True)
print(config.token)
#Strong bars:
#8444-12703
#8984-12704
#8135-6103
#9196-12701
#9886-12705
#Weak bars (some of these might not have enough signal to do anything with their velocity fields. Let me know if so and I’ll find some others):
#8987-12702
#Defining class to do linear regession using torch NN
class linearRegression(torch.nn.Module):
def __init__(self, inputSize, outputSize):
super(linearRegression, self).__init__()
self.linear = torch.nn.Linear(inputSize, outputSize)
def forward(self, x):
out = self.linear(x)
return out
plt.ion()
#set config attributes and turn on global downloads of Marvin data
config.setRelease('DR15')
config.mode = 'local'
config.download = True
########################################################### Question 1- Marvin ############################################
#Importing All MaNGA Data from DPRall Schema
data = pd.read_csv('CompleteTable.csv')
# my_cube1=Cube('7957-12703')
# central_spaxel1=my_cube1.getSpaxel(0,0)
# map1=my_cube1.getMaps()
# image1=Image('7957-12703')
# central_spaxel1.flux.plot()
# # plt.show()
def test_convert_tool_no_modelcube(self):
config.setRelease('MPL-4')
with self.assertRaises(MarvinError) as cm:
self._convertTool('modelcube')
errmsg = "ModelCube requires at least dapver='2.0.2'"
self.assertIn(errmsg, str(cm.exception))
def setUp(self):
super(TestIndexPage, self).setUp()
self.blue = 'index_page'
config.setRelease('MPL-5')
self.release = config.release
