def __init__(self, working_dir, view=None, imgData=None, maskData=None):
self.working_dir = working_dir
if view == None:
self.view = MaskView()
else:
self.view = view
if imgData == None:
self.img_data = ImgData()
else:
self.img_data = imgData
if maskData == None:
self.mask_data = MaskData()
else:
self.mask_data = maskData
self.view.img_view.mouse_left_clicked.connect(self.process_click)
self.state = None
self.clicks = 0
self.create_signals()
self.rect = None
self.circle = None
self.polygon = None
self.point = None
self.raise_window()
def setUp(self):
self.img_data = ImgData()
self.img_data.load('Data/Mg2SiO4_ambient_001.tif')
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data.load('Data/calibration.poni')
self.mask_data = MaskData()
self.mask_data.load_mask('Data/test.mask')
self.spectrum_data = SpectrumData()
class CombinedDataTest(unittest.TestCase):
def setUp(self):
self.img_data = ImgData()
self.img_data.load('Data/Mg2SiO4_ambient_001.tif')
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data.load('Data/calibration.poni')
self.mask_data = MaskData()
self.mask_data.mask_ellipse(500, 500, 100, 100)
self.spectrum_data = SpectrumData()
def test_dependencies(self):
tth1, int1 = self.calibration_data.integrate_1d()
self.img_data.load_next_file()
self.assertEqual(os.path.abspath(self.img_data.filename),
os.path.abspath('Data/Mg2SiO4_ambient_002.tif'))
tth2, int2 = self.calibration_data.integrate_1d()
self.assertFalse(np.array_equal(int1, int2))
plt.figure(1)
plt.plot(tth1, int1)
plt.plot(tth2, int2)
plt.savefig('Results/dependencies1.png')
tth3, int3 = self.calibration_data.integrate_1d(mask=self.mask_data.get_mask())
self.assertFalse(np.array_equal(int2, int3))
plt.figure(2)
plt.plot(tth2, int2)
plt.plot(tth3, int3)
plt.savefig('Results/dependencies2.png')
tth4, int4 = self.calibration_data.integrate_1d(polarization_factor=0.90, mask=None)
plt.figure(3)
plt.plot(tth2, int2)
plt.plot(tth4, int4)
plt.savefig('Results/dependencies3.png')
tth5, int5 = self.calibration_data.integrate_1d(polarization_factor=.5, mask=None)
plt.figure(4)
plt.plot(tth4, int4)
plt.plot(tth5, int5)
plt.savefig('Results/dependencies4.png')
def test_automatism(self):
def integrate_and_set_spectrum():
tth, I = self.calibration_data.integrate_1d()
self.spectrum_data.set_spectrum(tth, I, self.img_data.filename)
self.img_data.subscribe(integrate_and_set_spectrum)
y1 = self.spectrum_data.spectrum.data[1]
self.img_data.load_next_file()
y2 = self.spectrum_data.spectrum.data[1]
self.assertFalse(np.array_equal(y1, y2))
class CombinedDataTest(unittest.TestCase):
def setUp(self):
self.img_data = ImgData()
self.img_data.load('Data/Mg2SiO4_ambient_001.tif')
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data.load('Data/calibration.poni')
self.mask_data = MaskData()
self.mask_data.load_mask('Data/test.mask')
self.spectrum_data = SpectrumData()
def test_dependencies(self):
tth1, int1 = self.calibration_data.integrate_1d()
self.img_data.load_next()
tth2, int2 = self.calibration_data.integrate_1d()
self.assertFalse(np.array_equal(int1, int2))
plt.figure(1)
plt.plot(tth1, int1)
plt.plot(tth2, int2)
plt.savefig('Results/dependencies1.jpg')
tth3, int3 = self.calibration_data.integrate_1d(mask=self.mask_data.get_mask())
self.assertFalse(np.array_equal(int2, int3))
plt.figure(2)
plt.plot(tth2, int2)
plt.plot(tth3, int3)
plt.savefig('Results/dependencies2.jpg')
tth4, int4 = self.calibration_data.integrate_1d(polarization_factor=0.90, mask=None)
plt.figure(3)
plt.plot(tth2, int2)
plt.plot(tth4, int4)
plt.savefig('Results/dependencies3.jpg')
tth5, int5 = self.calibration_data.integrate_1d(polarization_factor=.5, mask=None)
plt.figure(4)
plt.plot(tth4, int4)
plt.plot(tth5, int5)
plt.savefig('Results/dependencies4.jpg')
def test_automatism(self):
def integrate_and_set_spectrum():
tth, I = self.calibration_data.integrate_1d()
self.spectrum_data.set_spectrum(tth, I, self.img_data.filename)
self.img_data.subscribe(integrate_and_set_spectrum)
y1 = self.spectrum_data.spectrum.data[1]
self.img_data.load_next()
y2 = self.spectrum_data.spectrum.data[1]
self.assertFalse(np.array_equal(y1, y2))
def __init__(self, app):
self.splash_img = QtGui.QPixmap("UiFiles/splash.png")
self.splash_screen = QtGui.QSplashScreen(
self.splash_img, QtCore.Qt.WindowStaysOnTopHint)
self.splash_screen.show()
app.processEvents()
time.sleep(1)
app.processEvents()
self.view = MainView()
#create data
self.img_data = ImgData()
self.calibration_data = CalibrationData(self.img_data)
self.mask_data = MaskData()
self.spectrum_data = SpectrumData()
self.phase_data = PhaseData()
self.load_directories()
#create controller
self.calibration_controller = CalibrationController(
self.working_dir, self.view.calibration_widget, self.img_data,
self.mask_data, self.calibration_data)
self.mask_controller = MaskController(self.working_dir,
self.view.mask_widget,
self.img_data, self.mask_data)
self.integration_controller = IntegrationController(
self.working_dir, self.view.integration_widget, self.img_data,
self.mask_data, self.calibration_data, self.spectrum_data,
self.phase_data)
self.create_signals()
self.set_title()
self.raise_window()
self.splash_screen.finish(self.view)
def __init__(self, working_dir, view=None, imgData=None, maskData=None):
self.working_dir = working_dir
if view == None:
self.view = MaskView()
else:
self.view = view
if imgData == None:
self.img_data = ImgData()
else:
self.img_data = imgData
if maskData == None:
self.mask_data = MaskData()
else:
self.mask_data = maskData
self.view.img_view.add_left_click_observer(self.process_click)
self.state = None
self.clicks = 0
self.create_signals()
self.rect = None
self.circle = None
self.polygon = None
self.point = None
self.raise_window()
def __init__(self):
self.view = MainView()
#create data
self.img_data = ImgData()
self.calibration_data = CalibrationData(self.img_data)
self.mask_data = MaskData()
self.spectrum_data = SpectrumData()
self.phase_data = PhaseData()
self.settings_directory = os.path.join(os.path.expanduser("~"), '.Dioptas')
self.working_directories = {'calibration': '', 'mask': '', 'image': '', 'spectrum': '', 'overlay': '',
'phase': ''}
self.load_settings()
#create controller
self.calibration_controller = CalibrationController(self.working_directories,
self.view.calibration_widget,
self.img_data,
self.mask_data,
self.calibration_data)
self.mask_controller = MaskController(self.working_directories,
self.view.mask_widget,
self.img_data,
self.mask_data)
self.integration_controller = IntegrationController(self.working_directories,
self.view.integration_widget,
self.img_data,
self.mask_data,
self.calibration_data,
self.spectrum_data,
self.phase_data)
self.create_signals()
self.set_title()
self.raise_window(self.view)
def __init__(self):
self.view = MainView()
#create data
self.img_data = ImgData()
self.calibration_data = CalibrationData(self.img_data)
self.mask_data = MaskData()
self.spectrum_data = SpectrumData()
self.phase_data = PhaseData()
self.load_directories()
#create controller
self.calibration_controller = CalibrationController(
self.working_dir, self.view.calibration_widget, self.img_data,
self.calibration_data)
self.mask_controller = MaskController(self.working_dir,
self.view.mask_widget,
self.img_data, self.mask_data)
self.integration_controller = IntegrationController(
self.working_dir, self.view.integration_widget, self.img_data,
self.mask_data, self.calibration_data, self.spectrum_data,
self.phase_data)
self.create_signals()
self.raise_window()
def setUp(self):
self.img_data = ImgData()
self.img_data.load('Data/Mg2SiO4_ambient_001.tif')
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data.load('Data/calibration.poni')
self.mask_data = MaskData()
self.mask_data.mask_ellipse(500, 500, 100, 100)
self.spectrum_data = SpectrumData()
def setUp(self):
self.app = QtGui.QApplication(sys.argv)
self.img_data = ImgData()
self.mask_data = MaskData()
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data._calibrants_working_dir = 'Data/Calibrants'
self.calibration_view = CalibrationView()
self.working_dir = {}
self.calibration_controller = CalibrationController(working_dir=self.working_dir,
img_data=self.img_data,
mask_data=self.mask_data,
view=self.calibration_view,
calibration_data=self.calibration_data)
def __init__(self,
working_dir,
view=None,
img_data=None,
mask_data=None,
calibration_data=None,
spectrum_data=None,
phase_data=None):
self.working_dir = working_dir
if view == None:
self.view = IntegrationView()
else:
self.view = view
if img_data == None:
self.img_data = ImgData()
else:
self.img_data = img_data
if mask_data == None:
self.mask_data = MaskData()
else:
self.mask_data = mask_data
if calibration_data == None:
self.calibration_data = CalibrationData(self.img_data)
else:
self.calibration_data = calibration_data
if spectrum_data == None:
self.spectrum_data = SpectrumData()
else:
self.spectrum_data = spectrum_data
if phase_data == None:
self.phase_data = PhaseData()
else:
self.phase_data = phase_data
self.create_sub_controller()
self.view.setWindowState(self.view.windowState()
& ~QtCore.Qt.WindowMinimized
| QtCore.Qt.WindowActive)
self.view.activateWindow()
self.view.raise_()
class MaskController(object):
def __init__(self, working_dir, view=None, imgData=None, maskData=None):
self.working_dir = working_dir
if view == None:
self.view = MaskView()
else:
self.view = view
if imgData == None:
self.img_data = ImgData()
else:
self.img_data = imgData
if maskData == None:
self.mask_data = MaskData()
else:
self.mask_data = maskData
self.view.img_view.mouse_left_clicked.connect(self.process_click)
self.state = None
self.clicks = 0
self.create_signals()
self.rect = None
self.circle = None
self.polygon = None
self.point = None
self.raise_window()
def raise_window(self):
self.view.show()
self.view.setWindowState(self.view.windowState()
& ~QtCore.Qt.WindowMinimized
| QtCore.Qt.WindowActive)
self.view.activateWindow()
self.view.raise_()
def connect_click_function(self, emitter, function):
self.view.connect(emitter, QtCore.SIGNAL('clicked()'), function)
def create_signals(self):
self.connect_click_function(self.view.circle_btn,
self.activate_circle_btn)
self.connect_click_function(self.view.rectangle_btn,
self.activate_rectangle_btn)
self.connect_click_function(self.view.polygon_btn,
self.activate_polygon_btn)
self.connect_click_function(self.view.point_btn,
self.activate_point_btn)
self.connect_click_function(self.view.undo_btn, self.undo_btn_click)
self.connect_click_function(self.view.redo_btn, self.redo_btn_click)
self.connect_click_function(self.view.below_thresh_btn,
self.below_thresh_btn_click)
self.connect_click_function(self.view.above_thresh_btn,
self.above_thresh_btn_click)
self.connect_click_function(self.view.cosmic_btn,
self.cosmic_btn_click)
self.connect_click_function(self.view.invert_mask_btn,
self.invert_mask_btn_click)
self.connect_click_function(self.view.clear_mask_btn,
self.clear_mask_btn_click)
self.connect_click_function(self.view.save_mask_btn,
self.save_mask_btn_click)
self.connect_click_function(self.view.load_mask_btn,
self.load_mask_btn_click)
self.connect_click_function(self.view.add_mask_btn,
self.add_mask_btn_click)
self.connect_click_function(self.view.mask_rb, self.mask_rb_click)
self.connect_click_function(self.view.unmask_rb, self.unmask_rb_click)
self.connect_click_function(self.view.fill_rb, self.fill_rb_click)
self.connect_click_function(self.view.transparent_rb,
self.transparent_rb_click)
self.view.connect(self.view.point_size_sb,
QtCore.SIGNAL('valueChanged(int)'),
self.set_point_size)
self.view.img_view.mouse_moved.connect(self.show_img_mouse_position)
self.view.keyPressEvent = self.key_press_event
def uncheck_all_btn(self, except_btn=None):
btns = [self.view.circle_btn, self.view.rectangle_btn, self.view.polygon_btn, \
self.view.point_btn]
for btn in btns:
if btn is not except_btn:
if btn.isChecked():
btn.toggle()
# if not except_btn.isChecked() and except_btn is not None:
# except_btn.toggle()
shapes = [self.rect, self.circle, self.polygon]
for shape in shapes:
if shape is not None:
self.view.img_view.img_view_box.removeItem(shape)
self.view.img_view.mouse_moved.disconnect(shape.set_size)
try:
self.view.img_view.mouse_moved.disconnect(self.point.set_position)
self.view.img_view.img_view_box.removeItem(self.point)
self.point = None
except AttributeError:
pass
def activate_circle_btn(self):
if self.view.circle_btn.isChecked():
self.state = 'circle'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.circle_btn)
else:
print 'hmm'
self.state = None
self.clicks = 0
self.uncheck_all_btn()
def activate_rectangle_btn(self):
if self.view.rectangle_btn.isChecked():
self.state = 'rectangle'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.rectangle_btn)
else:
self.state = None
self.uncheck_all_btn()
def activate_polygon_btn(self):
if self.view.polygon_btn.isChecked():
self.state = 'polygon'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.polygon_btn)
else:
self.state = None
self.uncheck_all_btn()
def activate_point_btn(self):
if self.view.point_btn.isChecked():
self.state = 'point'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.point_btn)
self.point = self.view.img_view.draw_point(
self.view.point_size_sb.value())
self.view.img_view.mouse_moved.connect(self.point.set_position)
else:
self.state = 'None'
self.uncheck_all_btn()
def undo_btn_click(self):
self.mask_data.undo()
self.view.img_view.plot_mask(self.mask_data.get_img())
def redo_btn_click(self):
self.mask_data.redo()
self.view.img_view.plot_mask(self.mask_data.get_img())
def plot_image(self):
self.view.img_view.plot_image(self.img_data.get_img_data(), False)
self.view.img_view.auto_range()
def process_click(self, x, y):
if self.state == 'circle':
self.draw_circle(x, y)
elif self.state == 'rectangle':
self.draw_rectangle(x, y)
elif self.state == 'point':
self.draw_point(x, y)
elif self.state == 'polygon':
self.draw_polygon(x, y)
def draw_circle(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.circle = self.view.img_view.draw_circle(x, y)
self.view.img_view.mouse_moved.connect(self.circle.set_size)
elif self.clicks == 1:
self.clicks = 0
self.mask_data.mask_QGraphicsEllipseItem(self.circle)
self.view.img_view.img_view_box.removeItem(self.circle)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.mouse_moved.disconnect(self.circle.set_size)
self.circle = None
def draw_rectangle(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.rect = self.view.img_view.draw_rectangle(x, y)
self.view.img_view.mouse_moved.connect(self.rect.set_size)
elif self.clicks == 1:
self.clicks = 0
self.mask_data.mask_QGraphicsRectItem(self.rect)
self.view.img_view.img_view_box.removeItem(self.rect)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.mouse_moved.disconnect(self.rect.set_size)
self.rect = None
def draw_point(self, x, y):
radius = self.view.point_size_sb.value()
self.mask_data.mask_ellipse(y, x, radius, radius)
self.view.img_view.plot_mask(self.mask_data.get_img())
def set_point_size(self, radius):
try:
self.point.set_radius(radius)
except AttributeError:
pass
def draw_polygon(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.polygon = self.view.img_view.draw_polygon(x, y)
self.view.img_view.mouse_moved.connect(self.polygon.set_size)
self.view.img_view.mouse_left_double_clicked.connect(
self.finish_polygon)
elif self.clicks == 1:
self.polygon.set_size(x, y)
self.polygon.add_point(x, y)
def finish_polygon(self, x, y):
self.view.img_view.mouse_moved.disconnect(self.polygon.set_size)
self.view.img_view.mouse_left_double_clicked.disconnect(
self.finish_polygon)
self.polygon.add_point(y, x)
self.clicks = 0
self.mask_data.mask_QGraphicsPolygonItem(self.polygon)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.img_view_box.removeItem(self.polygon)
self.polygon = None
def below_thresh_btn_click(self):
thresh = np.float64(self.view.below_thresh_txt.text())
self.mask_data.mask_below_threshold(self.img_data.get_img_data(),
thresh)
self.view.img_view.plot_mask(self.mask_data.get_img())
def above_thresh_btn_click(self):
thresh = np.float64(self.view.above_thresh_txt.text())
self.mask_data.mask_above_threshold(self.img_data.get_img_data(),
thresh)
self.view.img_view.plot_mask(self.mask_data.get_img())
def invert_mask_btn_click(self):
self.mask_data.invert_mask()
self.view.img_view.plot_mask(self.mask_data.get_img())
def clear_mask_btn_click(self):
self.mask_data.clear_mask()
self.view.img_view.plot_mask(self.mask_data.get_img())
def cosmic_btn_click(self):
self.mask_data.remove_cosmic(self.img_data.get_img_data())
self.view.img_view.plot_mask(self.mask_data.get_img())
def save_mask_btn_click(self, filename=None):
if filename is None:
filename = str(
QtGui.QFileDialog.getSaveFileName(
self.view,
caption="Save mask data",
directory=self.working_dir['mask'],
filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
np.savetxt(filename, self.mask_data.get_mask(), fmt="%d")
def load_mask_btn_click(self, filename=None):
if filename is None:
filename = str(
QtGui.QFileDialog.getOpenFileName(
self.view,
caption="Load mask data",
directory=self.working_dir['mask'],
filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
mask_data = np.loadtxt(filename)
if self.img_data.get_img_data().shape == mask_data.shape:
self.mask_data.set_mask(np.loadtxt(filename))
self.plot_mask()
else:
QtGui.QMessageBox.critical(
self.view, 'Error',
'Image data and mask data in selected file do not have '
'the same shape. Mask could not be loaded.')
def add_mask_btn_click(self, filename=None):
if filename is None:
filename = str(
QtGui.QFileDialog.getOpenFileName(
self.view,
caption="Add mask data",
directory=self.working_dir['mask'],
filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
mask_data = np.loadtxt(filename)
if self.mask_data.get_mask().shape == mask_data.shape:
self.mask_data.add_mask(np.loadtxt(filename))
self.plot_mask()
else:
QtGui.QMessageBox.critical(
self.view, 'Error',
'Image data and mask data in selected file do not have '
'the same shape. Mask could not be added.')
def plot_mask(self):
self.view.img_view.plot_mask(self.mask_data.get_mask())
def key_press_event(self, ev):
if self.state == "point":
if ev.text() == 'q':
self.view.point_size_sb.setValue(
self.view.point_size_sb.value() + 1)
if ev.text() == 'w':
self.view.point_size_sb.setValue(
self.view.point_size_sb.value() - 1)
if ev.modifiers() == QtCore.Qt.ControlModifier:
if ev.key() == 90: #for pressing z
self.undo_btn_click()
elif ev.key() == 89: #for pressing y
self.redo_btn_click()
elif ev.key() == 83: #for pressing s
self.save_mask_btn_click()
elif ev.key == 79: #for pressing o
self.load_mask_btn_click()
elif ev.key == 65:
self.add_mask_btn_click()
def mask_rb_click(self):
self.mask_data.set_mode(True)
def unmask_rb_click(self):
self.mask_data.set_mode(False)
def fill_rb_click(self):
self.view.img_view.set_color([255, 0, 0, 255])
self.plot_mask()
#
def transparent_rb_click(self):
self.view.img_view.set_color([255, 0, 0, 100])
self.plot_mask()
def show_img_mouse_position(self, x, y):
try:
if x > 0 and y > 0:
str = "x: %8.1f y: %8.1f I: %6.f" % (
x, y, self.view.img_view.img_data.T[np.floor(x),
np.floor(y)])
else:
str = "x: %.1f y: %.1f" % (x, y)
except (IndexError, AttributeError):
str = "x: %.1f y: %.1f" % (x, y)
self.view.pos_lbl.setText(str)
class MaskController(object):
def __init__(self, working_dir, view=None, imgData=None, maskData=None):
self.working_dir = working_dir
if view == None:
self.view = MaskView()
else:
self.view = view
if imgData == None:
self.img_data = ImgData()
else:
self.img_data = imgData
if maskData == None:
self.mask_data = MaskData()
else:
self.mask_data = maskData
self.view.img_view.add_left_click_observer(self.process_click)
self.state = None
self.clicks = 0
self.create_signals()
self.rect = None
self.circle = None
self.polygon = None
self.point = None
self.raise_window()
def raise_window(self):
self.view.show()
self.view.setWindowState(self.view.windowState() & ~QtCore.Qt.WindowMinimized | QtCore.Qt.WindowActive)
self.view.activateWindow()
self.view.raise_()
def connect_click_function(self, emitter, function):
self.view.connect(emitter, QtCore.SIGNAL('clicked()'), function)
def create_signals(self):
self.connect_click_function(self.view.circle_btn, self.activate_circle_btn)
self.connect_click_function(self.view.rectangle_btn, self.activate_rectangle_btn)
self.connect_click_function(self.view.polygon_btn, self.activate_polygon_btn)
self.connect_click_function(self.view.point_btn, self.activate_point_btn)
self.connect_click_function(self.view.undo_btn, self.undo_btn_click)
self.connect_click_function(self.view.redo_btn, self.redo_btn_click)
self.connect_click_function(self.view.below_thresh_btn, self.below_thresh_btn_click)
self.connect_click_function(self.view.above_thresh_btn, self.above_thresh_btn_click)
self.connect_click_function(self.view.cosmic_btn, self.cosmic_btn_click)
self.connect_click_function(self.view.invert_mask_btn, self.invert_mask_btn_click)
self.connect_click_function(self.view.clear_mask_btn, self.clear_mask_btn_click)
self.connect_click_function(self.view.save_mask_btn, self.save_mask_btn_click)
self.connect_click_function(self.view.load_mask_btn, self.load_mask_btn_click)
self.connect_click_function(self.view.add_mask_btn, self.add_mask_btn_click)
self.connect_click_function(self.view.mask_rb, self.mask_rb_click)
self.connect_click_function(self.view.unmask_rb, self.unmask_rb_click)
self.connect_click_function(self.view.fill_rb, self.fill_rb_click)
self.connect_click_function(self.view.transparent_rb, self.transparent_rb_click)
self.view.connect(self.view.point_size_sb, QtCore.SIGNAL('valueChanged(int)'), self.set_point_size)
self.view.keyPressEvent = self.key_press_event
def uncheck_all_btn(self, except_btn=None):
btns = [self.view.circle_btn, self.view.rectangle_btn, self.view.polygon_btn, \
self.view.point_btn]
for btn in btns:
if btn is not except_btn:
if btn.isChecked():
btn.toggle()
# if not except_btn.isChecked() and except_btn is not None:
# except_btn.toggle()
shapes = [self.rect, self.circle, self.polygon]
for shape in shapes:
if shape is not None:
self.view.img_view.img_view_box.removeItem(shape)
self.view.img_view.del_mouse_move_observer(shape.set_size)
try:
self.view.img_view.del_mouse_move_observer(self.point.set_position)
self.view.img_view.img_view_box.removeItem(self.point)
self.point = None
except AttributeError:
pass
def activate_circle_btn(self):
if self.view.circle_btn.isChecked():
self.state = 'circle'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.circle_btn)
else:
print 'hmm'
self.state = None
self.clicks = 0
self.uncheck_all_btn()
def activate_rectangle_btn(self):
if self.view.rectangle_btn.isChecked():
self.state = 'rectangle'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.rectangle_btn)
else:
self.state = None
self.uncheck_all_btn()
def activate_polygon_btn(self):
if self.view.polygon_btn.isChecked():
self.state = 'polygon'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.polygon_btn)
else:
self.state = None
self.uncheck_all_btn()
def activate_point_btn(self):
if self.view.point_btn.isChecked():
self.state = 'point'
self.clicks = 0
self.uncheck_all_btn(except_btn=self.view.point_btn)
self.point = self.view.img_view.draw_point(self.view.point_size_sb.value())
self.view.img_view.add_mouse_move_observer(self.point.set_position)
else:
self.state = 'None'
self.uncheck_all_btn()
def undo_btn_click(self):
self.mask_data.undo()
self.view.img_view.plot_mask(self.mask_data.get_img())
def redo_btn_click(self):
self.mask_data.redo()
self.view.img_view.plot_mask(self.mask_data.get_img())
def plot_image(self):
self.view.img_view.plot_image(self.img_data.get_img_data(), False)
self.view.img_view.auto_range()
def process_click(self, x, y):
if self.state == 'circle':
self.draw_circle(x, y)
elif self.state == 'rectangle':
self.draw_rectangle(x, y)
elif self.state == 'point':
self.draw_point(x, y)
elif self.state == 'polygon':
self.draw_polygon(x, y)
def draw_circle(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.circle = self.view.img_view.draw_circle(x, y)
self.view.img_view.add_mouse_move_observer(self.circle.set_size)
elif self.clicks == 1:
self.clicks = 0
self.mask_data.mask_QGraphicsEllipseItem(self.circle)
self.view.img_view.img_view_box.removeItem(self.circle)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.del_mouse_move_observer(self.circle.set_size)
self.circle = None
def draw_rectangle(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.rect = self.view.img_view.draw_rectangle(x, y)
self.view.img_view.add_mouse_move_observer(self.rect.set_size)
elif self.clicks == 1:
self.clicks = 0
self.mask_data.mask_QGraphicsRectItem(self.rect)
self.view.img_view.img_view_box.removeItem(self.rect)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.del_mouse_move_observer(self.rect.set_size)
self.rect = None
def draw_point(self, x, y):
radius = self.view.point_size_sb.value()
self.mask_data.mask_ellipse(y, x, radius, radius)
self.view.img_view.plot_mask(self.mask_data.get_img())
def set_point_size(self, radius):
try:
self.point.set_radius(radius)
except AttributeError:
pass
def draw_polygon(self, x, y):
if self.clicks == 0:
self.clicks += 1
self.polygon = self.view.img_view.draw_polygon(x, y)
self.view.img_view.add_mouse_move_observer(self.polygon.set_size)
self.view.img_view.add_left_double_click_observer(self.finish_polygon)
elif self.clicks == 1:
self.polygon.set_size(x, y)
self.polygon.add_point(x, y)
def finish_polygon(self, x, y):
self.view.img_view.del_mouse_move_observer(self.polygon.set_size)
self.view.img_view.del_left_double_click_observer(self.finish_polygon)
self.polygon.add_point(y, x)
self.clicks = 0
self.mask_data.mask_QGraphicsPolygonItem(self.polygon)
self.view.img_view.plot_mask(self.mask_data.get_img())
self.view.img_view.img_view_box.removeItem(self.polygon)
self.polygon = None
def below_thresh_btn_click(self):
thresh = np.float64(self.view.below_thresh_txt.text())
self.mask_data.mask_below_threshold(self.img_data.get_img_data(), thresh)
self.view.img_view.plot_mask(self.mask_data.get_img())
def above_thresh_btn_click(self):
thresh = np.float64(self.view.above_thresh_txt.text())
self.mask_data.mask_above_threshold(self.img_data.get_img_data(), thresh)
self.view.img_view.plot_mask(self.mask_data.get_img())
def invert_mask_btn_click(self):
self.mask_data.invert_mask()
self.view.img_view.plot_mask(self.mask_data.get_img())
def clear_mask_btn_click(self):
self.mask_data.clear_mask()
self.view.img_view.plot_mask(self.mask_data.get_img())
def cosmic_btn_click(self):
self.mask_data.remove_cosmic(self.img_data.get_img_data())
self.view.img_view.plot_mask(self.mask_data.get_img())
def save_mask_btn_click(self, filename=None):
if filename is None:
filename = str(QtGui.QFileDialog.getSaveFileName(self.view, caption="Save mask data",
directory=self.working_dir['mask'], filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
np.savetxt(filename, self.mask_data.get_mask(), fmt="%d")
def load_mask_btn_click(self, filename=None):
if filename is None:
filename = str(QtGui.QFileDialog.getOpenFileName(self.view, caption="Load mask data",
directory=self.working_dir['mask'], filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
mask_data = np.loadtxt(filename)
if self.img_data.get_img_data().shape == mask_data.shape:
self.mask_data.set_mask(np.loadtxt(filename))
self.plot_mask()
else:
QtGui.QMessageBox.critical(self.view, 'Error', 'Image data and mask data in selected file do not have '
'the same shape. Mask could not be loaded.')
def add_mask_btn_click(self, filename=None):
if filename is None:
filename = str(QtGui.QFileDialog.getOpenFileName(self.view, caption="Add mask data",
directory=self.working_dir['mask'], filter='*.mask'))
if filename is not '':
self.working_dir['mask'] = os.path.dirname(filename)
mask_data = np.loadtxt(filename)
if self.mask_data.get_mask().shape == mask_data.shape:
self.mask_data.add_mask(np.loadtxt(filename))
self.plot_mask()
else:
QtGui.QMessageBox.critical(self.view, 'Error', 'Image data and mask data in selected file do not have '
'the same shape. Mask could not be added.')
def plot_mask(self):
self.view.img_view.plot_mask(self.mask_data.get_mask())
def key_press_event(self, ev):
if self.state == "point":
if ev.text() == 'q':
self.view.point_size_sb.setValue(self.view.point_size_sb.value() + 1)
if ev.text() == 'w':
self.view.point_size_sb.setValue(self.view.point_size_sb.value() - 1)
if ev.modifiers() == QtCore.Qt.ControlModifier:
if ev.key() == 90: #for pressing z
self.undo_btn_click()
elif ev.key() == 89: #for pressing y
self.redo_btn_click()
elif ev.key() == 83: #for pressing s
self.save_mask_btn_click()
elif ev.key == 79: #for pressing o
self.load_mask_btn_click()
elif ev.key == 65:
self.add_mask_btn_click()
def mask_rb_click(self):
self.mask_data.set_mode(True)
def unmask_rb_click(self):
self.mask_data.set_mode(False)
def fill_rb_click(self):
self.view.img_view.set_color([255, 0, 0, 255])
self.plot_mask()
#
def transparent_rb_click(self):
self.view.img_view.set_color([255, 0, 0, 100])
self.plot_mask()
class CalibrationControllerTest(unittest.TestCase):
def setUp(self):
self.app = QtGui.QApplication(sys.argv)
self.img_data = ImgData()
self.mask_data = MaskData()
self.calibration_data = CalibrationData(self.img_data)
self.calibration_data._calibrants_working_dir = 'Data/Calibrants'
self.calibration_view = CalibrationView()
self.working_dir = {}
self.calibration_controller = CalibrationController(working_dir=self.working_dir,
img_data=self.img_data,
mask_data=self.mask_data,
view=self.calibration_view,
calibration_data=self.calibration_data)
def tearDown(self):
del self.app
def test_automatic_calibration1(self):
self.calibration_controller.load_img('Data/LaB6_p49_40keV_006.tif')
self.calibration_controller.search_peaks(1179.6, 1129.4)
self.calibration_controller.search_peaks(1268.5, 1119.8)
self.calibration_controller.view.sv_wavelength_txt.setText('0.31')
self.calibration_controller.view.sv_distance_txt.setText('200')
self.calibration_controller.view.sv_pixel_width_txt.setText('79')
self.calibration_controller.view.sv_pixel_height_txt.setText('79')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(7)
self.calibration_controller.calibrate()
self.calibration_controller.view.cake_view.set_vertical_line_pos(1419.8, 653.4)
def test_automatic_calibration2(self):
self.calibration_controller.load_img('Data/LaB6_WOS_30keV_005.tif')
self.calibration_controller.search_peaks(1245.2, 1919.3)
self.calibration_controller.search_peaks(1334.0, 1823.7)
self.calibration_controller.view.sv_wavelength_txt.setText('0.3344')
self.calibration_controller.view.sv_distance_txt.setText('500')
self.calibration_controller.view.sv_pixel_width_txt.setText('200')
self.calibration_controller.view.sv_pixel_height_txt.setText('200')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(7)
self.calibration_controller.calibrate()
self.calibration_controller.view.cake_view.set_vertical_line_pos(206.5, 171.6)
def test_automatic_calibration3(self):
self.calibration_controller.load_img('Data/CeO2_Oct24_2014_001_0000.tif')
QTest.mouseClick(self.calibration_view.automatic_peak_num_inc_cb, QtCore.Qt.LeftButton)
self.assertFalse(self.calibration_view.automatic_peak_num_inc_cb.isChecked())
self.calibration_controller.search_peaks(517.664434674, 647.529865592)
self.calibration_controller.search_peaks(667.380513299, 525.252854758)
self.calibration_controller.search_peaks(671.110095329, 473.571503774)
self.calibration_controller.search_peaks(592.788872703, 350.495296791)
self.calibration_controller.search_peaks(387.395462348, 390.987901686)
self.calibration_controller.search_peaks(367.94835605, 554.290314848)
self.calibration_controller.view.sv_wavelength_txt.setText('0.406626')
self.calibration_controller.view.sv_distance_txt.setText('200')
self.calibration_controller.view.sv_pixel_width_txt.setText('172')
self.calibration_controller.view.sv_pixel_height_txt.setText('172')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(4)
start_values = self.calibration_view.get_start_values()
self.assertAlmostEqual(start_values['wavelength'], 0.406626e-10)
self.assertAlmostEqual(start_values['pixel_height'], 172e-6)
self.assertAlmostEqual(start_values['pixel_width'], 172e-6)
self.calibration_controller.load_calibrant()
self.assertAlmostEqual(self.calibration_data.calibrant.wavelength, 0.406626e-10)
QTest.mouseClick(self.calibration_view.integrate_btn, QtCore.Qt.LeftButton)
calibration_parameter = self.calibration_data.get_calibration_parameter()[0]
self.assertAlmostEqual(calibration_parameter['dist'], .2086, places=4)
def test_automatic_calibration_with_supersampling(self):
self.calibration_controller.load_img('Data/LaB6_p49_40keV_006.tif')
self.calibration_controller.search_peaks(1179.6, 1129.4)
self.calibration_controller.search_peaks(1268.5, 1119.8)
self.calibration_controller.view.sv_wavelength_txt.setText('0.31')
self.calibration_controller.view.sv_distance_txt.setText('200')
self.calibration_controller.view.sv_pixel_width_txt.setText('79')
self.calibration_controller.view.sv_pixel_height_txt.setText('79')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(7)
self.img_data.set_supersampling(2)
self.calibration_data.set_supersampling(2)
self.calibration_controller.calibrate()
def test_automatic_calibration_with_supersampling_and_mask(self):
self.calibration_controller.load_img('Data/LaB6_p49_40keV_006.tif')
self.calibration_controller.search_peaks(1179.6, 1129.4)
self.calibration_controller.search_peaks(1268.5, 1119.8)
self.calibration_controller.view.sv_wavelength_txt.setText('0.31')
self.calibration_controller.view.sv_distance_txt.setText('200')
self.calibration_controller.view.sv_pixel_width_txt.setText('79')
self.calibration_controller.view.sv_pixel_height_txt.setText('79')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(7)
self.img_data.set_supersampling(2)
self.mask_data.mask_below_threshold(self.img_data._img_data, 1)
self.mask_data.set_supersampling(2)
self.calibration_data.set_supersampling(2)
self.calibration_controller.view.use_mask_cb.setChecked(True)
self.calibration_controller.calibrate()
def test_calibrating_one_image_size_and_loading_different_image_size(self):
self.calibration_controller.load_img('Data/LaB6_p49_40keV_006.tif')
self.calibration_controller.search_peaks(1179.6, 1129.4)
self.calibration_controller.search_peaks(1268.5, 1119.8)
self.calibration_controller.view.sv_wavelength_txt.setText('0.31')
self.calibration_controller.view.sv_distance_txt.setText('200')
self.calibration_controller.view.sv_pixel_width_txt.setText('79')
self.calibration_controller.view.sv_pixel_height_txt.setText('79')
self.calibration_controller.view.calibrant_cb.setCurrentIndex(7)
self.calibration_controller.view.options_automatic_refinement_cb.setChecked(False)
self.calibration_controller.view.use_mask_cb.setChecked(True)
self.calibration_controller.calibrate()
self.calibration_data.integrate_1d()
self.calibration_data.integrate_2d()
self.calibration_controller.load_img('Data/CeO2_Pilatus1M.tif')
self.calibration_data.integrate_1d()
self.calibration_data.integrate_2d()
def test_loading_and_saving_of_calibration_files(self):
self.calibration_controller.load_calibration('Data/calibration.poni')
self.calibration_controller.save_calibration('Data/calibration.poni')
class MainController(object):
"""
Creates a the main controller for Dioptas. Loads all the data objects and connects them with the other controllers
"""
def __init__(self):
self.view = MainView()
#create data
self.img_data = ImgData()
self.calibration_data = CalibrationData(self.img_data)
self.mask_data = MaskData()
self.spectrum_data = SpectrumData()
self.phase_data = PhaseData()
self.settings_directory = os.path.join(os.path.expanduser("~"), '.Dioptas')
self.working_directories = {'calibration': '', 'mask': '', 'image': '', 'spectrum': '', 'overlay': '',
'phase': ''}
self.load_settings()
#create controller
self.calibration_controller = CalibrationController(self.working_directories,
self.view.calibration_widget,
self.img_data,
self.mask_data,
self.calibration_data)
self.mask_controller = MaskController(self.working_directories,
self.view.mask_widget,
self.img_data,
self.mask_data)
self.integration_controller = IntegrationController(self.working_directories,
self.view.integration_widget,
self.img_data,
self.mask_data,
self.calibration_data,
self.spectrum_data,
self.phase_data)
self.create_signals()
self.set_title()
self.raise_window(self.view)
@staticmethod
def raise_window(widget):
widget.show()
widget.setWindowState(widget.windowState() & ~QtCore.Qt.WindowMinimized | QtCore.Qt.WindowActive)
widget.activateWindow()
widget.raise_()
def create_signals(self):
self.view.tabWidget.currentChanged.connect(self.tab_changed)
self.view.closeEvent = self.close_event
self.img_data.subscribe(self.set_title)
self.spectrum_data.subscribe(self.set_title)
def tab_changed(self, ind):
if ind == 2:
self.mask_data.set_supersampling()
self.integration_controller.image_controller.plot_mask()
self.integration_controller.view.calibration_lbl.setText(self.calibration_data.calibration_name)
auto_scale_previous = self.integration_controller.image_controller._auto_scale
self.integration_controller.image_controller._auto_scale = False
self.integration_controller.spectrum_controller.image_changed()
self.integration_controller.image_controller.update_img()
self.integration_controller.image_controller._auto_scale = auto_scale_previous
elif ind == 1:
self.mask_controller.plot_mask()
self.mask_controller.plot_image()
elif ind == 0:
self.calibration_controller.plot_mask()
try:
self.calibration_controller.update_calibration_parameter_in_view()
except (TypeError, AttributeError):
pass
def set_title(self):
img_filename = os.path.basename(self.img_data.filename)
spec_filename = os.path.basename(self.spectrum_data.spectrum_filename)
calibration_name = self.calibration_data.calibration_name
str = 'Dioptas ' + __VERSION__
if img_filename is '' and spec_filename is '':
self.view.setWindowTitle(str + u' - © 2014 C. Prescher')
self.view.integration_widget.img_frame.setWindowTitle(str + u' - © 2014 C. Prescher')
return
if img_filename is not '' or spec_filename is not '':
str += ' - ['
if img_filename is not '':
str += img_filename
elif img_filename is '' and spec_filename is not '':
str += spec_filename
if not img_filename == spec_filename:
str += ', ' + spec_filename
if calibration_name is not None:
str += ', calibration: ' + calibration_name
str += ']'
str += u' - © 2014 C. Prescher'
self.view.setWindowTitle(str)
self.view.integration_widget.img_frame.setWindowTitle(str)
def load_settings(self):
if os.path.exists(self.settings_directory):
self.load_directories()
self.load_xml_settings()
def load_directories(self):
working_directories_path = os.path.join(self.settings_directory, 'working_directories.csv')
if os.path.exists(working_directories_path):
reader = csv.reader(open(working_directories_path, 'r'))
self.working_directories = dict(x for x in reader)
def load_xml_settings(self):
xml_settings_path = os.path.join(self.settings_directory, "settings.xml")
if os.path.exists(xml_settings_path):
tree = ET.parse(xml_settings_path)
root = tree.getroot()
filenames = root.find("filenames")
calibration_path=filenames.find("calibration").text
if os.path.exists(str(calibration_path)):
self.calibration_data.load(calibration_path)
def save_settings(self):
if not os.path.exists(self.settings_directory):
os.mkdir(self.settings_directory)
self.save_directories()
self.save_xml_settings()
def save_directories(self):
working_directories_path = os.path.join(self.settings_directory, 'working_directories.csv')
writer = csv.writer(open(working_directories_path, 'w'))
for key, value in list(self.working_directories.items()):
writer.writerow([key, value])
writer.writerow([key, value])
def save_xml_settings(self):
root = ET.Element("DioptasSettings")
filenames = ET.SubElement(root, "filenames")
calibration_filename = ET.SubElement(filenames, "calibration")
calibration_filename.text = self.calibration_data.filename
tree = ET.ElementTree(root)
tree.write(os.path.join(self.settings_directory, "settings.xml"))
def close_event(self, _):
self.save_settings()
QtGui.QApplication.closeAllWindows()
QtGui.QApplication.quit()
