Exemplo n.º 1
0
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()
Exemplo n.º 2
0
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()