def read_strong_i_from_xds_ascii(xds_ascii_in): tmp = XDS_ASCII(xds_ascii_in, i_only=True).i_obs(anomalous_flag=False) sel = tmp.sigmas() > 0 sel &= tmp.data() / tmp.sigmas() > 2 sel &= tmp.d_spacings() > 3 if sel.count(True) < 10: return None tmp = tmp.select(sel) merge = tmp.merge_equivalents(use_internal_variance=False) return merge.array()
def read_strong_i_from_xds_ascii(xds_ascii_in): tmp = XDS_ASCII(xds_ascii_in, i_only=True).i_obs(anomalous_flag=False) sel = tmp.sigmas() > 0 sel &= tmp.data()/tmp.sigmas() > 2 sel &= tmp.d_spacings() > 3 if sel.count(True) < 10: return None tmp = tmp.select(sel) merge = tmp.merge_equivalents(use_internal_variance=False) return merge.array()