Advertisement
Guest User

code

a guest
Feb 23rd, 2024
26
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.65 KB | None | 0 0
  1. df = gaia_v2 #change to desired dataframe
  2. twomass_matches = []
  3. for i in range(0, len(df)):
  4. result = Vizier(columns=["**"], catalog="II/246").query_region(SkyCoord(ra=df['ra'][i], dec=df['dec'][i],
  5. unit=(u.deg, u.deg)), radius='20s')
  6. if len(result) == 0:
  7. values = [np.nan for i in range(0, 62)]
  8. twomass_matches.append(values)
  9. else:
  10. result = result[0]
  11.  
  12. if len(result) == 1:
  13. values = [result[colname][0] for colname in result.columns]
  14. twomass_matches.append(values)
  15. else:
  16. diffs = [(((df['ra'][1]-result['RAJ2000'][j])**2)+(df['dec'][i]-result['DEJ2000'][j])**2) for j in range(0, len(result))]
  17. index = diffs.index(min(diffs))
  18. values = [result[colname][index] for colname in result.columns]
  19. twomass_matches.append(values)
  20.  
  21. colnames = result.columns
  22. for i in range(0, len(colnames)):
  23. values = [twomass_matches[j][i] for j in range(0, len(twomass_matches))]
  24. df[colnames[i]] = values
  25. gaia_2mass = df
  26.  
  27. gaia_2mass.to_csv("2MASS_data") #insert desired filename to download as csv
  28.  
  29.  
  30. print (gaia_2mass)
  31.  
  32.  
  33.  
  34. #full error is
  35. AttributeError Traceback (most recent call last)
  36. Cell In[53], line 21
  37. 18 values = [result[colname][index] for colname in result.columns]
  38. 19 twomass_matches.append(values)
  39. ---> 21 colnames = result.columns
  40. 22 for i in range(0, len(colnames)):
  41. 23 values = [twomass_matches[j][i] for j in range(0, len(twomass_matches))]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement