datos = fits.open(\'/home/citlali/Documentos/Servicio/Lista.fits\')\ndata = datos[1].data\n\n\n#Linea [SIII] 9532\nMask_1 = data[\'flux_[SIII]9531.1_Re_fit\'] / data[\'e_flux_[SIII]9531.1_Re_fit\'] > 5\nnewdata1 = data[Mask_1]\ndat_flux = newdata1[\'flux_[SIII]9069.0_Re_fit\']\ndat_eflux = newdata1[\'e_flux_[SIII]9069.0_Re_fit\']\nMask_2 = dat_flux / dat_eflux > 5\nnewdata2 = newdata1[Mask_2]\n\n\nH1_alpha = newdata1[\'log_NII_Ha_Re\']\nH1_beta = newdata1[\'log_OIII_Hb_Re\']\nH2_alpha = newdata2[\'log_NII_Ha_Re\']\nH2_beta = newdata2[\'log_OIII_Hb_Re\']\n\n\nM = H1_alpha < -0.9\nnewx = H1_alpha[M]\nnewy = H1_beta[M] \nex = newx \ney = newy \n#print("Elementos de SIII [9532]: ", len(newx))\nm = H2_alpha < -0.9\nnewxm = H2_alpha[m]\nnewym = H2_beta[m] \n#print("Elementos de SIII [9069]: ", len(newxm))\n\nsm = heapq.nsmallest(3000, zip(newx, newy)) # zip them to sort together\nnewx, newy …
Run Code Online (Sandbox Code Playgroud)