“样本match”的版本间差异

2023年7月17日 (一) 13:12的最新版本

ID1=list(zip(Brick,BOID))
ID2=list(zip(Brick2,BOID2))
xy,x_ind,y_ind=np.intersect1d(ID1,ID2,assume_unique=True,return_indices=True)

参见[1]

@@ 第1行： / 第1行： @@
+==python==
-==一维数组的match==
+ ID1=list(zip(Brick,BOID))
-[[astro_lib]]中的match和match2
+ ID2=list(zip(Brick2,BOID2))
-*match,（数组中不能有重复）
+ xy,x_ind,y_ind=np.intersect1d(ID1,ID2,assume_unique=True,return_indices=True)
-CALLING SEQUENCE:
+参见[https://numpy.org/doc/stable/reference/generated/numpy.intersect1d.html]
- match, a, b, suba, subb, [ COUNT =, /SORT, EPSILON =  ]
-''超大整数的match，比如SDSS的objID会出问题''
+==[[idl中数组match]]==
-*match2可以用在某一样本是另外一个样本的子样本，这样可以完全匹配
- match2,a,b,suba,subb
-如果a完全是b的子集，那么b[suba]=a
-'''该程序可用于找出子样本的补集'''
- a=[3,2,1]
- b=indgen(5)
- match2,a,b,suba,subb
- print,b[where(subb eq -1)]
-===多维数据的match===
-:matchnd
-==天球坐标的match==
-*[[idlutils]]中的spherematch
-:spherematch,Ra1,Dec1,Ra2,Dec2,matchlength,ID1,ID2,dis
-:'''注意matchlength的单位是度'''
-:我自己根据spherematch写了一个spherematch_kth,可以自动搜索在某个距离范围内，第kth个match
-:[[idlutils]]里面有djs_ANGLE_NMATCH 和djs_ANGLE_MATCH程序，功能类似。
-*[[sdssidl]]中有一个CLOSE_MATCH_RADEC