Plotter2D.py @ 7fb69f26

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload

Last change on this file since 7fb69f26 was 0aca693, checked in by Jae Cho <jhjcho@…>, 12 years ago
Added showing "DataInfor?" from data panel as well as from plot
Property mode set to `100644`
File size: 28.8 KB

Rev	Line
[1bf33c1]	1
[d955bf19]	2	################################################################################
	3	#This software was developed by the University of Tennessee as part of the
	4	#Distributed Data Analysis of Neutron Scattering Experiments (DANSE)
	5	#project funded by the US National Science Foundation.
	6	#
	7	#See the license text in license.txt
	8	#
	9	#copyright 2008, University of Tennessee
	10	################################################################################
[1bf33c1]	11
	12
	13	import wx
[32c0841]	14	import sys
[ed8ad21]	15	import os
[32c0841]	16	import math
[003fa4e]	17	import numpy
[0d9dae8]	18	import pylab
[1bf33c1]	19	import danse.common.plottools
	20	from danse.common.plottools.PlotPanel import PlotPanel
[4ac8556]	21	from danse.common.plottools.plottables import Graph
[49a470d]	22	from danse.common.plottools.TextDialog import TextDialog
[55a0dc1]	23	from sans.guiframe.events import EVT_NEW_PLOT
	24	from sans.guiframe.events import EVT_SLICER_PARS
	25	from sans.guiframe.events import StatusEvent
	26	from sans.guiframe.events import NewPlotEvent
[a07e72f]	27	from sans.guiframe.events import PanelOnFocusEvent
[55a0dc1]	28	from sans.guiframe.events import SlicerEvent
[0d9dae8]	29	from sans.guiframe.utils import PanelMenu
[1bf33c1]	30	from binder import BindArtist
	31	from Plotter1D import ModelPanel1D
[32c0841]	32	from danse.common.plottools.toolbar import NavigationToolBar
[4ac8556]	33	from sans.guiframe.dataFitting import Data1D
[49a470d]	34	from matplotlib.font_manager import FontProperties
[32c0841]	35	(InternalEvent, EVT_INTERNAL) = wx.lib.newevent.NewEvent()
[0d9dae8]	36
	37	DEFAULT_QMAX = 0.05
[1bf33c1]	38	DEFAULT_QSTEP = 0.001
	39	DEFAULT_BEAM = 0.005
[ef0c170]	40	BIN_WIDTH = 1.0
[32c0841]	41
[d955bf19]	42
[ac8671e]	43	class NavigationToolBar2D(NavigationToolBar):
[d955bf19]	44	"""
	45	"""
[ac8671e]	46	def __init__(self, canvas, parent=None):
	47	NavigationToolBar.__init__(self, canvas=canvas, parent=parent)
	48
	49	def delete_option(self):
	50	"""
[d955bf19]	51	remove default toolbar item
[ac8671e]	52	"""
	53	#delete reset button
	54	self.DeleteToolByPos(0)
	55	#delete dragging
[88ca6db]	56	self.DeleteToolByPos(2)
[ac8671e]	57	#delete unwanted button that configures subplot parameters
	58	self.DeleteToolByPos(4)
	59
	60	def add_option(self):
	61	"""
[d955bf19]	62	add item to the toolbar
[ac8671e]	63	"""
[dc51a7f]	64	#add button
	65	id_context = wx.NewId()
	66	context_tip = 'Graph Menu'
	67	context = wx.ArtProvider.GetBitmap(wx.ART_LIST_VIEW, wx.ART_TOOLBAR)
	68	self.InsertSimpleTool(0, id_context, context,
	69	context_tip, context_tip)
	70	wx.EVT_TOOL(self, id_context, self.parent.onToolContextMenu)
	71	self.InsertSeparator(1)
[ac8671e]	72	#add print button
	73	id_print = wx.NewId()
	74	print_bmp = wx.ArtProvider.GetBitmap(wx.ART_PRINT, wx.ART_TOOLBAR)
	75	self.AddSimpleTool(id_print, print_bmp,
	76	'Print', 'Activate printing')
	77	wx.EVT_TOOL(self, id_print, self.on_print)
	78
	79
	80	class ModelPanel2D(ModelPanel1D):
[1bf33c1]	81	"""
[d955bf19]	82	Plot panel for use with the GUI manager
[1bf33c1]	83	"""
	84
	85	## Internal name for the AUI manager
	86	window_name = "plotpanel"
	87	## Title to appear on top of the window
	88	window_caption = "Plot Panel"
	89	## Flag to tell the GUI manager that this panel is not
	90	# tied to any perspective
	91	ALWAYS_ON = True
	92	## Group ID
	93	group_id = None
	94
	95
[32c0841]	96	def __init__(self, parent, id=-1, data2d=None, color = None,
	97	dpi=None, style=wx.NO_FULL_REPAINT_ON_RESIZE, **kwargs):
[1bf33c1]	98	"""
[d955bf19]	99	Initialize the panel
[1bf33c1]	100	"""
[32c0841]	101	ModelPanel1D.__init__(self, parent, id=id, style=style, **kwargs)
[1bf33c1]	102
	103	## Reference to the parent window
	104	self.parent = parent
[6c0568b]	105	## Dictionary containing Plottables
[1bf33c1]	106	self.plots = {}
[6c0568b]	107	## Save reference of the current plotted
	108	self.data2D = data2d
[1bf33c1]	109	## Unique ID (from gui_manager)
	110	self.uid = None
	111	## Action IDs for internal call-backs
	112	self.action_ids = {}
[6c0568b]	113	## Create Artist and bind it
[1bf33c1]	114	self.connect = BindArtist(self.subplot.figure)
[6c0568b]	115	## Beam stop
[1bf33c1]	116	self.beamstop_radius = DEFAULT_BEAM
[6c0568b]	117	## to set the order of lines drawn first.
[f15ed33]	118	self.slicer_z = 5
[6c0568b]	119	## Reference to the current slicer
[1bf33c1]	120	self.slicer = None
[6c0568b]	121	## event to send slicer info
[d468daa]	122	self.Bind(EVT_INTERNAL, self._onEVT_INTERNAL)
[1bf33c1]	123
[6c0568b]	124	self.axes_frozen = False
	125	## panel that contains result from slicer motion (ex: Boxsum info)
[32c0841]	126	self.panel_slicer = None
[49a470d]	127	self.title_label = None
	128	self.title_font = None
	129	self.title_color = 'black'
[1bf33c1]	130	## Graph
	131	self.graph = Graph()
	132	self.graph.xaxis("\\rm{Q}", 'A^{-1}')
[32c0841]	133	self.graph.yaxis("\\rm{Intensity} ", "cm^{-1}")
[1bf33c1]	134	self.graph.render(self)
[8dfdd20]	135	## store default value of zmin and zmax
	136	self.default_zmin_ctl = self.zmin_2D
	137	self.default_zmax_ctl = self.zmax_2D
[ac8671e]	138
[a07e72f]	139	def onLeftDown(self, event):
	140	"""
	141	left button down and ready to drag
	142
	143	"""
	144	# Check that the LEFT button was pressed
[53cf669]	145	PlotPanel.onLeftDown(self, event)
[e85e2bc]	146	ax = event.inaxes
	147	if ax != None:
	148	# data coordinate position
	149	pos_x = "%8.3g"% event.xdata
	150	pos_y = "%8.3g"% event.ydata
	151	position = "x: %s y: %s" % (pos_x, pos_y)
	152	wx.PostEvent(self.parent, StatusEvent(status=position))
[a07e72f]	153	self.plottable_selected(self.data2D.id)
	154	self._manager.set_panel_on_focus(self)
	155	wx.PostEvent(self.parent, PanelOnFocusEvent(panel=self))
	156
[ac8671e]	157	def add_toolbar(self):
	158	"""
[d955bf19]	159	add toolbar
[ac8671e]	160	"""
	161	self.enable_toolbar = True
[32c0841]	162	self.toolbar = NavigationToolBar2D(parent=self, canvas=self.canvas)
[ac8671e]	163	self.toolbar.Realize()
	164	# On Windows platform, default window size is incorrect, so set
	165	# toolbar width to figure width.
	166	tw, th = self.toolbar.GetSizeTuple()
	167	fw, fh = self.canvas.GetSizeTuple()
	168	self.toolbar.SetSize(wx.Size(fw, th))
[32c0841]	169	self.sizer.Add(self.toolbar, 0, wx.LEFT\|wx.EXPAND)
[ac8671e]	170	# update the axes menu on the toolbar
	171	self.toolbar.update()
	172
[a07e72f]	173	def plot_data(self, data):
[1bf33c1]	174	"""
[d955bf19]	175	Data is ready to be displayed
	176
	177	:TODO: this name should be changed to something more appropriate
	178	Don't forget that changing this name will mean changing code
	179	in plotting.py
	180
	181	:param event: data event
[1bf33c1]	182	"""
[e4a703a]	183	xlo = None
	184	xhi = None
	185	ylo = None
	186	yhi = None
[b837e1f]	187	if data.__class__.__name__ == 'Data1D':
	188	return
[6c0568b]	189	## Update self.data2d with the current plot
[857d00f]	190	self.data2D = data
[a07e72f]	191	if data.id in self.plots.keys():
	192	#replace
[e4a703a]	193	xlo, xhi = self.subplot.get_xlim()
	194	ylo, yhi = self.subplot.get_ylim()
[a07e72f]	195	self.graph.replace(data)
	196	self.plots[data.id] = data
[ab8f936]	197	else:
[a07e72f]	198	self.plots[data.id] = data
	199	self.graph.add(self.plots[data.id])
[ac9a5f6]	200	# update qmax with the new xmax of data plotted
[a07e72f]	201	self.qmax = data.xmax
[4b91fd1]	202
[32c0841]	203	self.slicer = None
[1bf33c1]	204	# Check axis labels
	205	#TODO: Should re-factor this
[6c0568b]	206	## render the graph with its new content
[7fff5cd]	207
	208	#data2D: put 'Pixel (Number)' for axis title and unit in case of having no detector info and none in _units
[49a470d]	209	if len(data.detector) < 1:
[a07e72f]	210	if len(data._xunit) < 1 and len(data._yunit) < 1:
	211	data._xaxis = '\\rm{x}'
	212	data._yaxis = '\\rm{y}'
	213	data._xunit = 'pixel'
	214	data._yunit = 'pixel'
[857d00f]	215
	216	# graph properties
[a07e72f]	217	self.graph.xaxis(data._xaxis, data._xunit)
	218	self.graph.yaxis(data._yaxis, data._yunit)
[36288ca]	219	if self._is_changed_legend_label:
	220	data.label = self.title_label
[49a470d]	221
[fe48fcc]	222	if data.label == None:
	223	data.label = data.name
[49a470d]	224
	225	if not self.title_font:
	226	self.graph.title(data.label)
	227	self.graph.render(self)
[857d00f]	228	# Set the axis labels on subplot
	229	self._set_axis_labels()
[49a470d]	230	self.draw_plot()
	231	else:
	232	self.graph.render(self)
	233	self.draw_plot()
	234	self.subplot.set_title(label=data.label,
	235	fontproperties=self.title_font,
	236	color=self.title_color)
	237	self.subplot.figure.canvas.draw_idle()
	238
[857d00f]	239
[8dfdd20]	240	## store default value of zmin and zmax
	241	self.default_zmin_ctl = self.zmin_2D
	242	self.default_zmax_ctl = self.zmax_2D
[7022fdc]	243	# Check if zoomed
	244	toolbar_zoomed = self.toolbar.GetToolEnabled(self.toolbar._NTB2_BACK)
	245	if not self.is_zoomed and not toolbar_zoomed:
	246	return
	247
[e4a703a]	248	# Recover the x,y limits
	249	if (xlo and xhi and ylo and yhi) != None:
	250	if (xlo > data.xmin and xhi < data.xmax and\
	251	ylo > data.ymin and yhi < data.ymax):
	252	self.subplot.set_xlim((xlo, xhi))
[fe48fcc]	253	self.subplot.set_ylim((ylo, yhi))
	254	else:
	255	self.toolbar.update()
[7022fdc]	256	self._is_zoomed = False
[53cf669]	257
	258	# Update Graph menu and help string
	259	pos = self.parent._window_menu.FindItem(self.window_caption)
	260	helpString = 'Show/Hide Graph: '
[d74fbff]	261	helpString += (' ' + str(data.label) +';')
[53cf669]	262	self.parent._window_menu.SetHelpString(pos, helpString)
[a4964b8e]	263
[857d00f]	264	def _set_axis_labels(self):
	265	"""
	266	Set axis labels
	267	"""
	268	data = self.data2D
	269	# control axis labels from the panel itself
	270	yname, yunits = data.get_yaxis()
	271	if self.yaxis_label != None:
	272	yname = self.yaxis_label
	273	yunits = self.yaxis_unit
	274	else:
	275	self.yaxis_label = yname
	276	self.yaxis_unit = yunits
	277	xname, xunits = data.get_xaxis()
	278	if self.xaxis_label != None:
	279	xname = self.xaxis_label
	280	xunits = self.xaxis_unit
	281	else:
	282	self.xaxis_label = xname
	283	self.xaxis_unit = xunits
[2a09236]	284	self.xaxis(xname, xunits, self.xaxis_font,
	285	self.xaxis_color, self.xaxis_tick)
	286	self.yaxis(yname, yunits, self.yaxis_font,
	287	self.yaxis_color, self.yaxis_tick)
[857d00f]	288
[1bf33c1]	289	def onContextMenu(self, event):
	290	"""
[d955bf19]	291	2D plot context menu
	292
	293	:param event: wx context event
[15550f4]	294
[d955bf19]	295	"""
[1bf33c1]	296	slicerpop = PanelMenu()
	297	slicerpop.set_plots(self.plots)
	298	slicerpop.set_graph(self.graph)
[9a585d0]	299
	300	id = wx.NewId()
[6d727ae]	301	slicerpop.Append(id, '&Save Image')
[9a585d0]	302	wx.EVT_MENU(self, id, self.onSaveImage)
	303
	304	id = wx.NewId()
[6d727ae]	305	slicerpop.Append(id,'&Print Image', 'Print image')
[9a585d0]	306	wx.EVT_MENU(self, id, self.onPrint)
	307
[1ce365f8]	308	id = wx.NewId()
[6d727ae]	309	slicerpop.Append(id,'&Print Preview', 'Print preview')
[1ce365f8]	310	wx.EVT_MENU(self, id, self.onPrinterPreview)
[003fa4e]	311
	312	id = wx.NewId()
	313	slicerpop.Append(id, '&Copy to Clipboard', 'Copy to the clipboard')
	314	wx.EVT_MENU(self, id, self.OnCopyFigureMenu)
[c5a769e]	315	slicerpop.AppendSeparator()
[ed8ad21]	316	# saving data
	317	plot = self.data2D
	318	id = wx.NewId()
	319	name = plot.name
[0aca693]	320	slicerpop.Append(id, "&Data Info" )
	321	wx.EVT_MENU(self, id, self._onDataShow)
	322
	323	id = wx.NewId()
	324	name = plot.name
[c17760d]	325	slicerpop.Append(id, "&Save as a File (DAT)" )
[ed8ad21]	326	self.action_ids[str(id)] = plot
	327	wx.EVT_MENU(self, id, self._onSave)
	328
[9a585d0]	329	slicerpop.AppendSeparator()
[7fff5cd]	330	if len(self.data2D.detector) == 1:
[0e13148]	331
[a07e72f]	332	item_list = self.parent.get_context_menu(self)
[6d727ae]	333	if (not item_list == None) and (not len(item_list) == 0) and\
	334	self.data2D.name.split(" ")[0] != 'Residuals':
	335	# The line above; Not for trunk
[32c0841]	336	for item in item_list:
	337	try:
	338	id = wx.NewId()
	339	slicerpop.Append(id, item[0], item[1])
	340	wx.EVT_MENU(self, id, item[2])
	341	except:
	342	msg = "ModelPanel1D.onContextMenu: "
	343	msg += "bad menu item %s"%sys.exc_value
	344	wx.PostEvent(self.parent, StatusEvent(status=msg))
	345	pass
	346	slicerpop.AppendSeparator()
[0e13148]	347
[15550f4]	348	id = wx.NewId()
[c17760d]	349	slicerpop.Append(id, '&Perform Circular Average')
[6d727ae]	350	wx.EVT_MENU(self, id, self.onCircular) \
	351	# For Masked Data
	352	if not plot.mask.all():
	353	id = wx.NewId()
[c17760d]	354	slicerpop.Append(id, '&Masked Circular Average')
[6d727ae]	355	wx.EVT_MENU(self, id, self.onMaskedCircular)
[15550f4]	356	id = wx.NewId()
[c17760d]	357	slicerpop.Append(id, '&Sector [Q View]')
[15550f4]	358	wx.EVT_MENU(self, id, self.onSectorQ)
	359	id = wx.NewId()
[c17760d]	360	slicerpop.Append(id, '&Annulus [Phi View ]')
[15550f4]	361	wx.EVT_MENU(self, id, self.onSectorPhi)
[92c2345]	362	id = wx.NewId()
[15550f4]	363	slicerpop.Append(id, '&Box Sum')
	364	wx.EVT_MENU(self, id, self.onBoxSum)
	365	id = wx.NewId()
[c17760d]	366	slicerpop.Append(id, '&Box Averaging in Qx')
[15550f4]	367	wx.EVT_MENU(self, id, self.onBoxavgX)
	368	id = wx.NewId()
[c17760d]	369	slicerpop.Append(id, '&Box Averaging in Qy')
[15550f4]	370	wx.EVT_MENU(self, id, self.onBoxavgY)
[32c0841]	371	if self.slicer != None:
[eba08f1a]	372	id = wx.NewId()
[c17760d]	373	slicerpop.Append(id, '&Clear Slicer')
[15550f4]	374	wx.EVT_MENU(self, id, self.onClearSlicer)
[32c0841]	375	if self.slicer.__class__.__name__ != "BoxSum":
[15550f4]	376	id = wx.NewId()
	377	slicerpop.Append(id, '&Edit Slicer Parameters')
	378	wx.EVT_MENU(self, id, self._onEditSlicer)
	379	slicerpop.AppendSeparator()
[fe48fcc]	380
	381	id = wx.NewId()
[dc51a7f]	382	slicerpop.Append(id, '&Edit Graph Label', 'Edit Graph Label')
[fe48fcc]	383	wx.EVT_MENU(self, id, self.onEditLabels)
	384	slicerpop.AppendSeparator()
	385
[0e13148]	386	id = wx.NewId()
[5637362]	387	slicerpop.Append(id, '&Edit Y Axis Label')
[857d00f]	388	wx.EVT_MENU(self, id, self._on_yaxis_label)
	389	id = wx.NewId()
[5637362]	390	slicerpop.Append(id, '&Edit X Axis Label')
[857d00f]	391	wx.EVT_MENU(self, id, self._on_xaxis_label)
[5637362]	392	id = wx.NewId()
	393	slicerpop.Append(id, '&Toggle Grid On/Off', 'Toggle Grid On/Off')
	394	wx.EVT_MENU(self, id, self.onGridOnOff)
[857d00f]	395	slicerpop.AppendSeparator()
	396
	397	id = wx.NewId()
[52f3c98]	398	slicerpop.Append(id, '&2D Color Map')
[32c0841]	399	wx.EVT_MENU(self, id, self._onEditDetector)
[857d00f]	400	slicerpop.AppendSeparator()
	401
[1bf33c1]	402	id = wx.NewId()
[c17760d]	403	slicerpop.Append(id, '&Toggle Linear/Log Scale')
[1bf33c1]	404	wx.EVT_MENU(self, id, self._onToggleScale)
[53cf669]	405
	406
	407	slicerpop.AppendSeparator()
	408	id = wx.NewId()
	409	slicerpop.Append(id, '&Window Title')
	410	wx.EVT_MENU(self, id, self.onChangeCaption)
	411
[dc51a7f]	412	try:
	413	pos_evt = event.GetPosition()
	414	pos = self.ScreenToClient(pos_evt)
	415	except:
	416	pos_x, pos_y = self.toolbar.GetPositionTuple()
	417	pos = (pos_x, pos_y + 5)
[1bf33c1]	418	self.PopupMenu(slicerpop, pos)
[37c36d9]	419
[fe48fcc]	420
	421	def onEditLabels(self, event):
	422	"""
	423	Edit legend label
	424	"""
	425	selected_plot = self.plots[self.graph.selected_plottable]
	426	label = selected_plot.label
[49a470d]	427	dial = TextDialog(None, -1, 'Change Label', label)
[fe48fcc]	428	if dial.ShowModal() == wx.ID_OK:
[49a470d]	429	try:
	430	FONT = FontProperties()
	431	newlabel = dial.getText()
	432	font = FONT.copy()
	433	font.set_size(dial.getSize())
	434	font.set_family(dial.getFamily())
	435	font.set_style(dial.getStyle())
	436	font.set_weight(dial.getWeight())
	437	colour = dial.getColor()
	438	if len(newlabel) > 0:
[53cf669]	439	# update Label
[49a470d]	440	selected_plot.label = newlabel
[53cf669]	441	self.graph.title(newlabel)
[49a470d]	442	self.title_label = selected_plot.label
	443	self.title_font = font
	444	self.title_color = colour
	445	## render the graph
	446	self.subplot.set_title(label=self.title_label,
	447	fontproperties=self.title_font,
	448	color=self.title_color)
[36288ca]	449	self._is_changed_legend_label = True
[49a470d]	450	self.subplot.figure.canvas.draw_idle()
	451	except:
	452	if self.parent != None:
	453	from sans.guiframe.events import StatusEvent
	454	msg= "Add Text: Error. Check your property values..."
	455	wx.PostEvent(self.parent, StatusEvent(status = msg ))
	456	else:
	457	raise
[fe48fcc]	458	dial.Destroy()
[53cf669]	459
	460	# Update Graph menu and help string
	461	if self.title_label != None:
	462	pos = self.parent._window_menu.FindItem(self.window_caption)
	463	helpString = 'Show/Hide Graph: '
[e2160ab]	464	helpString += (' ' + str(self.title_label) +';')
[53cf669]	465	self.parent._window_menu.SetHelpString(pos, helpString)
[49a470d]	466
[fe48fcc]	467
[ea290ee]	468	def _onEditDetector(self, event):
[6d920cd]	469	"""
[d955bf19]	470	Allow to view and edits detector parameters
	471
	472	:param event: wx.menu event
	473
[6d920cd]	474	"""
[ea290ee]	475	import detector_dialog
[8dfdd20]	476	dialog = detector_dialog.DetectorDialog(self, -1,base=self.parent,
	477	reset_zmin_ctl =self.default_zmin_ctl,
	478	reset_zmax_ctl = self.default_zmax_ctl,cmap=self.cmap)
[6c0568b]	479	## info of current detector and data2D
[ea290ee]	480	xnpts = len(self.data2D.x_bins)
	481	ynpts = len(self.data2D.y_bins)
	482	xmax = max(self.data2D.xmin, self.data2D.xmax)
	483	ymax = max(self.data2D.ymin, self.data2D.ymax)
[32c0841]	484	qmax = math.sqrt(math.pow(xmax, 2) + math.pow(ymax, 2))
[ea290ee]	485	beam = self.data2D.xmin
[6c0568b]	486	## set dialog window content
[ea290ee]	487	dialog.setContent(xnpts=xnpts,ynpts=ynpts,qmax=qmax,
	488	beam=self.data2D.xmin,
	489	zmin = self.zmin_2D,
	490	zmax = self.zmax_2D)
	491	if dialog.ShowModal() == wx.ID_OK:
	492	evt = dialog.getContent()
	493	self.zmin_2D = evt.zmin
	494	self.zmax_2D = evt.zmax
[32c0841]	495	self.cmap = evt.cmap
[ea290ee]	496	dialog.Destroy()
[6c0568b]	497	## Redraw the current image
[32c0841]	498	self.image(data=self.data2D.data,
[20b6760]	499	qx_data=self.data2D.qx_data,
	500	qy_data=self.data2D.qy_data,
[ea290ee]	501	xmin= self.data2D.xmin,
	502	xmax= self.data2D.xmax,
	503	ymin= self.data2D.ymin,
	504	ymax= self.data2D.ymax,
	505	zmin= self.zmin_2D,
	506	zmax= self.zmax_2D,
[8dfdd20]	507	cmap= self.cmap,
[32c0841]	508	color=0, symbol=0, label=self.data2D.name)
[ea290ee]	509	self.subplot.figure.canvas.draw_idle()
[1bf33c1]	510
	511	def freeze_axes(self):
[d955bf19]	512	"""
	513	"""
[1bf33c1]	514	self.axes_frozen = True
	515
	516	def thaw_axes(self):
[d955bf19]	517	"""
	518	"""
[1bf33c1]	519	self.axes_frozen = False
	520
	521	def onMouseMotion(self,event):
[d955bf19]	522	"""
	523	"""
[1bf33c1]	524	pass
[d955bf19]	525
[1bf33c1]	526	def onWheel(self, event):
[d955bf19]	527	"""
	528	"""
[6c0568b]	529	pass
	530
[1bf33c1]	531	def update(self, draw=True):
	532	"""
[d955bf19]	533	Respond to changes in the model by recalculating the
	534	profiles and resetting the widgets.
[1bf33c1]	535	"""
[6d727ae]	536	self.draw_plot()
[1bf33c1]	537
	538	def _getEmptySlicerEvent(self):
[6c0568b]	539	"""
[d955bf19]	540	create an empty slicervent
[6c0568b]	541	"""
[32c0841]	542	return SlicerEvent(type=None, params=None, obj_class=None)
[d955bf19]	543
[1bf33c1]	544	def _onEVT_INTERNAL(self, event):
	545	"""
[d955bf19]	546	Draw the slicer
	547
	548	:param event: wx.lib.newevent (SlicerEvent) containing slicer
[6c0568b]	549	parameter
[d955bf19]	550
[1bf33c1]	551	"""
	552	self._setSlicer(event.slicer)
	553
	554	def _setSlicer(self, slicer):
[6c0568b]	555	"""
[d955bf19]	556	Clear the previous slicer and create a new one.Post an internal
	557	event.
	558
	559	:param slicer: slicer class to create
	560
[6c0568b]	561	"""
	562	## Clear current slicer
[1bf33c1]	563	if not self.slicer == None:
	564	self.slicer.clear()
[6c0568b]	565	## Create a new slicer
[1bf33c1]	566	self.slicer_z += 1
	567	self.slicer = slicer(self, self.subplot, zorder=self.slicer_z)
[240c805]	568	self.subplot.set_ylim(self.data2D.ymin, self.data2D.ymax)
	569	self.subplot.set_xlim(self.data2D.xmin, self.data2D.xmax)
[6c0568b]	570	## Draw slicer
[1bf33c1]	571	self.update()
	572	self.slicer.update()
[32c0841]	573	msg = "Plotter2D._setSlicer %s"%self.slicer.__class__.__name__
	574	wx.PostEvent(self.parent, StatusEvent(status=msg))
[1bf33c1]	575	# Post slicer event
	576	event = self._getEmptySlicerEvent()
	577	event.type = self.slicer.__class__.__name__
	578	event.obj_class = self.slicer.__class__
	579	event.params = self.slicer.get_params()
[d468daa]	580	wx.PostEvent(self, event)
[6d727ae]	581
	582	def onMaskedCircular(self, event):
	583	"""
	584	perform circular averaging on Data2D with mask if it exists
	585
	586	:param event: wx.menu event
	587
	588	"""
	589	self.onCircular(event, True)
	590
	591	def onCircular(self, event, ismask=False):
[1bf33c1]	592	"""
[d955bf19]	593	perform circular averaging on Data2D
	594
	595	:param event: wx.menu event
	596
[1bf33c1]	597	"""
[003fa4e]	598	# Find the best number of bins
	599	npt = math.sqrt(len(self.data2D.data[numpy.isfinite(self.data2D.data)]))
	600	npt = math.floor(npt)
[fe857e2]	601	from sans.dataloader.manipulations import CircularAverage
[6c0568b]	602	## compute the maximum radius of data2D
[32c0841]	603	self.qmax = max(math.fabs(self.data2D.xmax),
	604	math.fabs(self.data2D.xmin))
	605	self.ymax = max(math.fabs(self.data2D.ymax),
	606	math.fabs(self.data2D.ymin))
	607	self.radius = math.sqrt(math.pow(self.qmax, 2)+ math.pow(self.ymax, 2))
[6c0568b]	608	##Compute beam width
[003fa4e]	609	bin_width = (self.qmax + self.qmax)/npt
[6c0568b]	610	## Create data1D circular average of data2D
[32c0841]	611	Circle = CircularAverage(r_min=0, r_max=self.radius,
	612	bin_width=bin_width)
[6d727ae]	613	circ = Circle(self.data2D, ismask=ismask)
[1bf33c1]	614	from sans.guiframe.dataFitting import Data1D
[32c0841]	615	if hasattr(circ, "dxl"):
	616	dxl = circ.dxl
[1bf33c1]	617	else:
[32c0841]	618	dxl = None
	619	if hasattr(circ, "dxw"):
	620	dxw = circ.dxw
[1bf33c1]	621	else:
[32c0841]	622	dxw = None
[003fa4e]	623
	624	new_plot = Data1D(x=circ.x, y=circ.y, dy=circ.dy, dx=circ.dx)
[32c0841]	625	new_plot.dxl = dxl
	626	new_plot.dxw = dxw
	627	new_plot.name = "Circ avg " + self.data2D.name
	628	new_plot.source = self.data2D.source
	629	#new_plot.info = self.data2D.info
[1bf33c1]	630	new_plot.interactive = True
[32c0841]	631	new_plot.detector = self.data2D.detector
[6d727ae]	632
[6c0568b]	633	## If the data file does not tell us what the axes are, just assume...
[32c0841]	634	new_plot.xaxis("\\rm{Q}", "A^{-1}")
[2ca51f44]	635	if hasattr(self.data2D, "scale") and \
	636	self.data2D.scale == 'linear':
	637	new_plot.ytransform = 'y'
	638	new_plot.yaxis("\\rm{Residuals} ", "normalized")
[6d727ae]	639	else:
	640	new_plot.yaxis("\\rm{Intensity} ", "cm^{-1}")
	641
[32c0841]	642	new_plot.group_id = "Circ avg " + self.data2D.name
	643	new_plot.id = "Circ avg " + self.data2D.name
	644	new_plot.is_data = True
[13382fc7]	645	self.parent.update_theory(data_id=self.data2D, \
	646	theory=new_plot)
[32c0841]	647	wx.PostEvent(self.parent,
	648	NewPlotEvent(plot=new_plot, title=new_plot.name))
[ef0c170]	649
[1bf33c1]	650	def _onEditSlicer(self, event):
[6c0568b]	651	"""
[d955bf19]	652	Is available only when a slicer is drawn.Create a dialog
	653	window where the user can enter value to reset slicer
	654	parameters.
	655
	656	:param event: wx.menu event
	657
[6c0568b]	658	"""
[32c0841]	659	if self.slicer != None:
[1bf33c1]	660	from SlicerParameters import SlicerParameterPanel
[4f8a00c]	661	dialog = SlicerParameterPanel(self, -1, "Slicer Parameters")
[1bf33c1]	662	dialog.set_slicer(self.slicer.__class__.__name__,
	663	self.slicer.get_params())
	664	if dialog.ShowModal() == wx.ID_OK:
	665	dialog.Destroy()
	666
	667	def onSectorQ(self, event):
	668	"""
[d955bf19]	669	Perform sector averaging on Q and draw sector slicer
[1bf33c1]	670	"""
[ef0c170]	671	from SectorSlicer import SectorInteractor
[1bf33c1]	672	self.onClearSlicer(event)
[32c0841]	673	wx.PostEvent(self, InternalEvent(slicer=SectorInteractor))
[1bf33c1]	674
	675	def onSectorPhi(self, event):
	676	"""
[d955bf19]	677	Perform sector averaging on Phi and draw annulus slicer
[1bf33c1]	678	"""
[ef0c170]	679	from AnnulusSlicer import AnnulusInteractor
[1bf33c1]	680	self.onClearSlicer(event)
[32c0841]	681	wx.PostEvent(self, InternalEvent(slicer=AnnulusInteractor))
[1bf33c1]	682
[d955bf19]	683	def onBoxSum(self, event):
	684	"""
	685	"""
[7ab9241]	686	from boxSum import BoxSum
	687	self.onClearSlicer(event)
[54cc36a]	688	self.slicer_z += 1
	689	self.slicer = BoxSum(self, self.subplot, zorder=self.slicer_z)
	690	self.subplot.set_ylim(self.data2D.ymin, self.data2D.ymax)
	691	self.subplot.set_xlim(self.data2D.xmin, self.data2D.xmax)
	692	self.update()
	693	self.slicer.update()
[6c0568b]	694	## Value used to initially set the slicer panel
	695	type = self.slicer.__class__.__name__
	696	params = self.slicer.get_params()
	697	## Create a new panel to display results of summation of Data2D
[8a7a21b]	698	from slicerpanel import SlicerPanel
[32c0841]	699	new_panel = SlicerPanel(parent=self.parent, id=-1,
	700	base=self, type=type,
	701	params=params, style=wx.RAISED_BORDER)
	702
	703	new_panel.window_caption = self.slicer.__class__.__name__ + " " + \
	704	str(self.data2D.name)
	705	new_panel.window_name = self.slicer.__class__.__name__+ " " + \
	706	str(self.data2D.name)
[6c0568b]	707	## Store a reference of the new created panel
[32c0841]	708	self.panel_slicer = new_panel
[6c0568b]	709	## save the window_caption of the new panel in the current slicer
[32c0841]	710	self.slicer.set_panel_name(name=new_panel.window_caption)
[6c0568b]	711	## post slicer panel to guiframe to display it
[55a0dc1]	712	from sans.guiframe.events import SlicerPanelEvent
[32c0841]	713	wx.PostEvent(self.parent, SlicerPanelEvent(panel=self.panel_slicer,
	714	main_panel=self))
[6c0568b]	715
[8a7a21b]	716	def onBoxavgX(self,event):
[6c0568b]	717	"""
[d955bf19]	718	Perform 2D data averaging on Qx
	719	Create a new slicer .
	720
	721	:param event: wx.menu event
[6c0568b]	722	"""
[8a7a21b]	723	from boxSlicer import BoxInteractorX
[38224f10]	724	self.onClearSlicer(event)
[32c0841]	725	wx.PostEvent(self, InternalEvent(slicer=BoxInteractorX))
[d468daa]	726
[8a7a21b]	727	def onBoxavgY(self,event):
[6c0568b]	728	"""
[d955bf19]	729	Perform 2D data averaging on Qy
	730	Create a new slicer .
	731
	732	:param event: wx.menu event
	733
[6c0568b]	734	"""
[8a7a21b]	735	from boxSlicer import BoxInteractorY
	736	self.onClearSlicer(event)
[32c0841]	737	wx.PostEvent(self, InternalEvent(slicer=BoxInteractorY))
[6c0568b]	738
[1bf33c1]	739	def onClearSlicer(self, event):
	740	"""
[d955bf19]	741	Clear the slicer on the plot
[1bf33c1]	742	"""
[32c0841]	743	if not self.slicer == None:
[1bf33c1]	744	self.slicer.clear()
	745	self.subplot.figure.canvas.draw()
	746	self.slicer = None
	747	# Post slicer None event
	748	event = self._getEmptySlicerEvent()
[d468daa]	749	wx.PostEvent(self, event)
[ed8ad21]	750
	751	def _onSave(self, evt):
	752	"""
	753	Save a data set to a dat(text) file
	754
	755	:param evt: Menu event
	756
	757	"""
	758	id = str(evt.GetId())
[c553b18]	759	if self.parent != None:
	760	self._default_save_location = self.parent._default_save_location
[3858e0f]	761	default_name = self.plots[self.graph.selected_plottable].label
	762	if default_name.count('.') > 0:
	763	default_name = default_name.split('.')[0]
	764	default_name += "_out"
[ed8ad21]	765	if id in self.action_ids:
	766	path = None
[176fbf1]	767	self.parent.save_data2d(self.data2D, default_name)
[0aca693]	768
	769	def _onDataShow(self, evt):
	770	"""
	771	Show the data set in text
	772
	773	:param evt: Menu event
	774
	775	"""
	776	menu = evt.GetEventObject()
	777	id = evt.GetId()
	778	self.set_selected_from_menu(menu, id)
	779	data = self.plots[self.graph.selected_plottable]
	780	default_name = data.label
	781	if default_name.count('.') > 0:
	782	default_name = default_name.split('.')[0]
	783	#default_name += "_out"
	784	if self.parent != None:
	785	self.parent.show_data2d(data, default_name)
[176fbf1]	786

Note: See TracBrowser for help on using the repository browser.

SasView

source: sasview/sansguiframe/src/sans/guiframe/local_perspectives/plotting/Plotter2D.py @ 7fb69f26

Download in other formats: