Exemple #1
0
  def tst_centroid_points3d(self):
    from dials.algorithms.image.centroid import centroid_points
    from scitbx import matrix

    centroid = centroid_points(self.pixels3d.as_1d(), self.points3d.as_1d())

    assert(abs(self.gold3d - matrix.col(centroid.mean())) < self.EPS)
    assert(abs(self.gold3dvar - matrix.col(centroid.variance())) < self.EPS)
    assert(abs(self.gold3dubvar - matrix.col(centroid.unbiased_variance())) < self.EPS)
    print 'OK'
Exemple #2
0
    def test_centroid_points3d(self):
        from dials.algorithms.image.centroid import centroid_points
        from scitbx import matrix

        centroid = centroid_points(self.pixels3d.as_1d(), self.points3d.as_1d())

        assert abs(self.gold3d - matrix.col(centroid.mean())) < self.EPS
        assert abs(self.gold3dvar - matrix.col(centroid.variance())) < self.EPS
        assert (
            abs(self.gold3dubvar - matrix.col(centroid.unbiased_variance())) < self.EPS
        )