In [17]:

    

folder = '/usr/users/spardy/data/Gals/'
m101_file = 'frame-u-003712-3-0187.fits'
ngc3906_file = 'NGC3906.stellar.fits'

import numpy as np
import matplotlib.pyplot as plt
import astropy.io
from astropy.io import fits
%matplotlib inline


    

In [3]:

    

m101_img = fits.open(folder+m101_file)


    

In [4]:

    

m101_img.info()


    

    




Filename: /usr/users/spardy/data/Gals/frame-u-003712-3-0187.fits
No.    Name         Type      Cards   Dimensions   Format
0    PRIMARY     PrimaryHDU      96   (2048, 1489)   float32   
1                ImageHDU         6   (2048,)      float32   
2                BinTableHDU     27   1R x 3C      [49152E, 2048E, 1489E]   
3                BinTableHDU     79   1R x 31C     [J, 3A, J, A, D, D, 2J, J, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, D, E, E]   

In [11]:

    

myarray = m101_img[0].data


    

In [12]:

    

myarray.shape


    

    
Out[12]:





(1489, 2048)

In [13]:

    

print myarray[1,4]


    

    




0.0927734

In [14]:

    

myarray.mean()


    

    
Out[14]:





0.047707528

In [15]:

    

myarray.max()


    

    
Out[15]:





155.5

In [16]:

    

myarray.min()


    

    
Out[16]:





-0.20556641

In [21]:

    

myarray.std()


    

    
Out[21]:





0.24568631

In [30]:

    

arraylog10 = np.log10(myarray)

arraylog10[np.isnan(arraylog10)] = 0.0

print arraylog10.min()
print arraylog10.max()
print arraylog10.mean()
print arraylog10.std()


    

    




-2.70337
2.19173
-1.03553
0.714922






    




/usr/users/spardy/.local/lib/python2.7/site-packages/IPython/kernel/__main__.py:1: RuntimeWarning: invalid value encountered in log10
  if __name__ == '__main__':

In [34]:

    

plt.imshow(arraylog10, vmin=0, vmax=1)


    

    
Out[34]:





<matplotlib.image.AxesImage at 0x9448c350>






    

In [35]:

    

ngc3906img = fits.open(folder+ngc3906_file)


    

In [36]:

    

ngc3906img.info()


    

    




Filename: /usr/users/spardy/data/Gals/NGC3906.stellar.fits
No.    Name         Type      Cards   Dimensions   Format
0    PRIMARY     PrimaryHDU     577   (569, 632)   float32   

In [37]:

    

ngcarray = ngc3906img[0].data


    

In [40]:

    

ngcarraylog10 = np.log10(ngcarray)


    

    




/usr/users/spardy/.local/lib/python2.7/site-packages/IPython/kernel/__main__.py:1: RuntimeWarning: divide by zero encountered in log10
  if __name__ == '__main__':
/usr/users/spardy/.local/lib/python2.7/site-packages/IPython/kernel/__main__.py:1: RuntimeWarning: invalid value encountered in log10
  if __name__ == '__main__':

In [41]:

    

plt.imshow(ngcarraylog10)


    

    
Out[41]:





<matplotlib.image.AxesImage at 0x960f0650>






    

In [ ]: