# Difference: JeanDuffyPythonCode (1 vs. 2)

#### Revision 225 Oct 2013 - JeanDuffy

Line: 1 to 1

 META TOPICPARENT name="JeanDuffyBScProject"
-- JeanDuffy - 25 Oct 2013
Deleted:
<
<

### Extraction of useful information from a fit file

#Change directory

Line: 65 to 64
plt.figure(1), plt.imshow(d2)
>
>
y=d.sum(1)

import numpy as np

x=np.arange(0,len(y),1)

xmean=(x[600:]*y[600:]).sum()/y[600:].sum()

print xmean

plt.figure(2),plt.plot(x,y)

plt.show()

#### Revision 125 Oct 2013 - JeanDuffy

Line: 1 to 1
>
>
 META TOPICPARENT name="JeanDuffyBScProject"
-- JeanDuffy - 25 Oct 2013

### Extraction of useful information from a fit file

#Change directory

import os

os.chdir("C:/Users/Jean/Desktop/Python/Bsc/Fits")

# Import fits from astropy.io

from astropy.io import fits

# Open the fit file "Atlas_pleione_1sec.fit" which

# opens it as an HDU (Header Data Unit) list and set to f.

f=fits.open("Atlas_pleione_1sec.fit")

# Set i to be an HDU object.

i=f[0]

# Set d to be an array of data.

d=i.data

# Imports matplotlib.pyplot module

import matplotlib.pyplot as plt

# Plots the array

plt.figure(0),plt.imshow(d)

# Creates a figure and shows the plot

# Not needed in the kernel (for some reason)

# Gets all the Header information from the fit file

# These can be run in the kernel

print d.mean()

print d.std()

print d.sum()

print d.sum(0)

print d.sum(1)

# Create a sub area

d2=d[700:800,850:950]

# Plot the subarea on a new figure

plt.figure(1), plt.imshow(d2)

plt.show()

Copyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding RHUL Physics Department TWiki? Send feedback