from __future__ import division
from __future__ import print_function
from scitbx.array_family import flex
from dials.model.data import Reflection, ReflectionList
from dials.algorithms import shoebox
rl = ReflectionList()
r1 = Reflection()
r1.shoebox = (10, 20, 10, 20, 10, 20)
r2 = Reflection()
r2.shoebox = (15, 25, 15, 25, 15, 25)
r3 = Reflection()
r3.shoebox = (20, 30, 20, 30, 20, 30)
rl.append(r1)
rl.append(r2)
rl.append(r3)
overlapping = shoebox.find_overlapping(rl)
for e in overlapping.edges():
print("Edge: ", overlapping.edge_vertices(e))
for v in overlapping.vertices():
print("Vertex: ", v, " => ", [a for a in overlapping.adjacent_vertices(v)])
from __future__ import division
from scitbx.array_family import flex
from dials.model.data import Reflection, ReflectionList
from dials.algorithms import shoebox
rl = ReflectionList()
r1 = Reflection()
r1.shoebox = (10, 20, 10, 20, 10, 20)
r2 = Reflection()
r2.shoebox = (15, 25, 15, 25, 15, 25)
r3 = Reflection()
r3.shoebox = (20, 30, 20, 30, 20, 30)
rl.append(r1)
rl.append(r2)
rl.append(r3)
overlapping = shoebox.find_overlapping(rl)
for e in overlapping.edges():
print "Edge: ", overlapping.edge_vertices(e)
for v in overlapping.vertices():
print "Vertex: ", v, " => ", [a for a in overlapping.adjacent_vertices(v)]
data2d = flex.double(flex.grid(1, 3, 3),15)
data2d[0, 1, 1] = 50
for row in range(3):
for col in range(3):
data2d[0,row, col] += row * 2
data2d[0,row, col] += col * 2
mask2d = flex.int(flex.grid(1, 3, 3),3)
mask2d[0, 1, 1] = 5
background2d = flex.double(flex.grid(1, 3, 3),0)
from dials.model.data import Reflection, ReflectionList
from scitbx.array_family import flex
r = Reflection()
r.shoebox = data2d
r.shoebox_mask = mask2d
r.shoebox_background = background2d
rlist = ReflectionList()
rlist.append(r)
#from dials.algorithms.background.flat_background_subtractor \
# import layering_and_background_avg
#layering_and_background_avg(rlist)
#from dials.algorithms.background.curved_background_subtractor \
# import layering_and_background_modl
#layering_and_background_modl(rlist)
from dials.algorithms.background.inclined_background_subtractor \
data2d[0, 1, 1] = 50
for row in range(3):
for col in range(3):
data2d[0, row, col] += row * 2
data2d[0, row, col] += col * 2
mask2d = flex.int(flex.grid(1, 3, 3), 3)
mask2d[0, 1, 1] = 5
background2d = flex.double(flex.grid(1, 3, 3), 0)
from dials.model.data import Reflection, ReflectionList
from scitbx.array_family import flex
r = Reflection()
r.shoebox = data2d
r.shoebox_mask = mask2d
r.shoebox_background = background2d
rlist = ReflectionList()
rlist.append(r)
# from dials.algorithms.background.flat_background_subtractor \
# import layering_and_background_avg
# layering_and_background_avg(rlist)
# from dials.algorithms.background.curved_background_subtractor \
# import layering_and_background_modl
# layering_and_background_modl(rlist)
from dials.algorithms.background.inclined_background_subtractor import (
