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