def getRMS(p,ra,dec):
e = p.split('.')
field = e[0]+'_'+e[1]+'_'+e[2]
name = field+'_850.sdf'
scuba = p+'_850.fits'
# print name,'\n',scuba
# ra = p[1]
# dec = p[2]
# if field != 'G001_36+20_96': #and field != 'G005_93-01_00':
# continue
# if field != 'G005_93-01_00':
# continue
center_deg = (ra_to_deg(ra),dec_to_deg(dec))
ra0,dec0 = center_deg
print ra,' ',dec
print center_deg
scubaFile_3D = scubaFits_in+scuba
f = fits.open(scubaFile_3D)#[0]
sc = FixFits(f) # sc is same type as f
fData = sc[0].data
fig = aplpy.FITSFigure(sc)
I,J = np.indices(fData.shape)
RA,D = fig.pixel2world(J,I)
dist = distance_on_sphere_deg(ra0,dec0,RA,D)
m = np.nonzero((dist<0.02)&(np.isfinite(fData)))
rms = np.std(fData[m])
rmsVal.write(p[0]+' '+str(rms)+'\n')
fig.close()
return rms
####
f = fits.open(filename)[0]
fData = f.data
######################
# fig = aplpy.FITSFigure(filename)
# J,I = np.indices(f.data.shape)
# RA,D = fig.pixel2world(J,I)
# fig.close()
RA, D = pixel_to_wcs(filename)
######################
# Calculate mean
dist = distance_on_sphere_deg(ra0, dec0, RA, D)
m = np.nonzero((dist < 0.02) & (np.isfinite(f.data)))
# m = np.nonzero((dist<0.02)&(np.isfinite(fData)))
# m = np.nonzero(np.isfinite(f.data))
# m = np.nonzero(dist<0.02)
# if scuba == 'G039.74+01.98_850.fits' and filename.endswith('_160_L2.5_J_res_1.fits'):
# print f.data[m]
# print f.data.shape
# print f.data[m].shape
e = np.mean(f.data[m])
# print e