From 4878b181fdb68fe2f51cb2b400976af58627897b Mon Sep 17 00:00:00 2001 From: rskoupy Date: Fri, 8 Dec 2023 19:01:04 +0100 Subject: [PATCH] After ER meeting --- .../ptychoScopy-checkpoint.ipynb | 360 +++++++++++++----- calibrations.xlsx | Bin 33397 -> 33439 bytes ptychoScopy.ipynb | 360 +++++++++++++----- 3 files changed, 546 insertions(+), 174 deletions(-) diff --git a/.ipynb_checkpoints/ptychoScopy-checkpoint.ipynb b/.ipynb_checkpoints/ptychoScopy-checkpoint.ipynb index 268e78d..1456a12 100644 --- a/.ipynb_checkpoints/ptychoScopy-checkpoint.ipynb +++ b/.ipynb_checkpoints/ptychoScopy-checkpoint.ipynb @@ -2,19 +2,10 @@ "cells": [ { "cell_type": "code", - "execution_count": 1, + "execution_count": 79, "id": "7425242d-3c91-4c1e-a424-08625a38ee7a", "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "C:\\Users\\skoupy_r\\Anaconda3\\lib\\site-packages\\scipy\\__init__.py:155: UserWarning: A NumPy version >=1.18.5 and <1.25.0 is required for this version of SciPy (detected version 1.26.1\n", - " warnings.warn(f\"A NumPy version >={np_minversion} and <{np_maxversion}\"\n" - ] - } - ], + "outputs": [], "source": [ "### Initial packages import ###################################################\n", "import numpy as np\n", @@ -25,7 +16,7 @@ "from pandas import read_excel\n", "from IPython.display import display\n", "from plotly.subplots import make_subplots\n", - "from ipywidgets import interactive_output, HBox, VBox, Layout, Label, Valid, ToggleButtons, RadioButtons, Dropdown, IntSlider, Checkbox, Button, ButtonStyle, GridBox\n", + "from ipywidgets import interactive_output, HBox, VBox, Layout, Label, Valid, ToggleButtons, RadioButtons, Dropdown, IntSlider, Checkbox, Button, ButtonStyle, GridBox,IntSlider\n", "\n", "# print(\"Numpy version \", np.__version__)\n", "# print(\"Pandas version \", pd.__version__)\n", @@ -300,31 +291,37 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 389, "id": "8055b802-cf83-4250-aea9-54e3e6b73db0", "metadata": {}, "outputs": [], "source": [ "### CONTROLS ##########################################\n", "\n", - "ali = dict(style = {'description_width': '60px','button_width': '33px'}, button_style='',disabled=False)\n", - "ali2 = dict(style = {'description_width': '60px','button_width': '90px'}, button_style='',disabled=False)\n", - "header = Button(description='STEM control panel',layout=Layout(width='auto', grid_area='header'), style=ButtonStyle(button_color='silver'))\n", + "ali = dict(style = {'description_width': '60px','button_width': '35px', 'font_weight': 'bold'}, button_style='',disabled = False)\n", + "ali2 = dict(style = {'description_width': '60px','button_width': '90px', 'font_weight': 'bold'}, button_style='',disabled = False)\n", "\n", - "beam_set = Button(description='Electron beam settings',layout=Layout(width='auto', grid_area='beam_set'), style=ButtonStyle(button_color='lightsteelblue')) \n", + "\n", + "\n", + "\n", + "\n", + "\n", + "header = Button(description='STEM control panel',layout=Layout(width='auto', grid_area='header'), style=ButtonStyle(button_color='silver', font_weight= 'bold',))\n", + "\n", + "beam_set = Button(description='Electron beam settings',layout=Layout(width='auto', grid_area='beam_set'), style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold')) \n", "beam = ToggleButtons(options=opt.energies(), value=200,description='E (keV)', layout=Layout(width='auto', grid_area='beam_set1'), **ali)\n", "aperture = ToggleButtons(options=opt.apertures(), value='30um', description='Aperture', layout=Layout(width='auto', grid_area='beam_set2'), **ali)\n", "probe = ToggleButtons(options=opt.probes(), value='8C', description='Probe', layout=Layout(width='auto', grid_area='beam_set3'), **ali)\n", "defocus = IntSlider(description='Δf (nm)', value=0, min=0, max=1000, step=1, continuous_update=False, layout=Layout(width='auto', grid_area='beam_set4'), **ali)\n", "\n", - "scanning_set = Button(description='Scanning parameters',layout=Layout(width='auto', grid_area='scanning_set'),style=ButtonStyle(button_color='lightsteelblue'))\n", + "scanning_set = Button(description='Scanning parameters',layout=Layout(width='auto', grid_area='scanning_set'),style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold'))\n", "mag = ToggleButtons(options=opt.magnifications(),value=3,description='Mag. Mx', layout=Layout(width='auto', grid_area='scanning_set1'), **ali)\n", "matrix = ToggleButtons(options=opt.mappings(), value=2048, description='Matrix', layout=Layout(width='auto', grid_area='scanning_set2'), **ali)\n", "dwell_time = ToggleButtons(options=opt.dwelltimes(), value=10, description='t (μs)', layout=Layout(width='auto', grid_area='scanning_set3'), **ali)\n", "\n", "\n", "\n", - "camera_set = Button(description='Detection',layout=Layout(width='auto', grid_area='camera_set'),style=ButtonStyle(button_color='lightsteelblue'))\n", + "camera_set = Button(description='Detection',layout=Layout(width='auto', grid_area='camera_set'),style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold'))\n", "cl = ToggleButtons(options=opt.cameralengths(), value=12, description='CL (cm)', layout=Layout(width='auto', grid_area='cl_set1'), **ali)\n", "camera = ToggleButtons(options=opt.detectors(), description='Detector', layout=Layout(width='auto', grid_area='camera_set1'), **ali2)\n", "restriction = ToggleButtons(options=[('.',False), ('..',True)], value = False, description='PAAR',icons = ['times','check'], layout=Layout(width='auto', grid_area='camera_set2'), **ali)\n", @@ -348,7 +345,7 @@ "footer = Button(description='Footer',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='silver'))\n", "\n", "\n", - "controls = GridBox(children=[header, beam_set, beam, aperture, probe, defocus, scanning_set, mag,matrix,dwell_time, cl,camera_set, camera, restriction, binning, beam_res, aperture_res,aperture_res2,fov_res, dwell_time_res,footer],\n", + "controls = GridBox(children=[header, beam_set, beam, aperture, probe, defocus, scanning_set, mag,matrix,dwell_time, cl,camera_set, camera, restriction, binning, beam_res, aperture_res,aperture_res2,fov_res, dwell_time_res],\n", " layout=Layout(\n", " width='22%',\n", " grid_template_rows='auto auto auto',\n", @@ -375,49 +372,56 @@ "\n", "# SMALL CONTROLS\n", "\n", - "descr = '50px'\n", - "butt = '30px'\n", + "descr = '80px'\n", + "butt = '60px'\n", "\n", - "header = Button(description='Graph controls',layout=Layout(width='auto', grid_area='header'),style=ButtonStyle(button_color='whitesmoke'))\n", + "header = Button(description='Graph controls',layout=Layout(width='auto', grid_area='header'),style=ButtonStyle(button_color='whitesmoke', font_weight= 'bold'))\n", "\n", - "method = ToggleButtons(options=[('SSB','direct'), ('Iter','iterative')], description='Method', tooltips=['SSB', 'PIE, MLc, DM'],\n", - " layout=Layout(width='auto', grid_area='methods'), style = {'description_width': descr,'button_width': butt}, button_style='',disabled=False,)\n", + "method = ToggleButtons(options=[('SSB','direct'), ('Iterative','iterative')], description='Method', tooltips=['SSB', 'PIE, MLc, DM'],\n", + " layout=Layout(width='auto', grid_area='methods'), style = {'description_width': descr,'button_width': butt, 'font_weight': 'bold'}, button_style='',disabled=False,)\n", "\n", - "graph_size = ToggleButtons(options=[('XXS', 50),('XS', 70),('S', 90), ('M', 100), ('L', 120)], value=100, description='Graphs size',\n", + "graph_size = IntSlider(value=100,min=10,max=200,step=5, description='Graphs size',\n", " layout=Layout(width='auto', grid_area='gr_size'), style = {'description_width': descr,'button_width': butt}, button_style='',disabled=False,)\n", "\n", "ctf_xaxis = RadioButtons(options=['α','mrad', 'A'], description='CTF x-axis:',\n", " layout=Layout(width='auto', grid_area='ctxaxis'), style = {'description_width': descr,'button_width': butt}, disabled=False,)\n", "\n", "cl_check_laa = Checkbox(value=True, description='Low anlge approximation', indent=False,\n", - " layout=Layout(width='auto', grid_area='check1'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Low angle approximation check - scattering angle < 10 deg\n", - "cl_check_pix = Checkbox(value=True, description='Recon. pixel < step size',indent=False,\n", - " layout=Layout(width='auto', grid_area='check2'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Reconstructed pixel size < scanning step size \n", + " layout=Layout(width='auto', grid_area='check1'), style = {'description_width': descr}, disabled=False,) # Low angle approximation check - scattering angle < 10 deg\n", + "cl_check_pix = Checkbox(value=True, description='Ptycho pixel < step size',indent=False,\n", + " layout=Layout(width='auto', grid_area='check2'), style = {'description_width': descr}, disabled=False,) # Reconstructed pixel size < scanning step size \n", "cl_check_def = Checkbox(value=True, description='Beam fits in probe window', indent=False,\n", - " layout=Layout(width='auto', grid_area='check3'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Is probe window big enough to accomodate defocused probe?\n", + " layout=Layout(width='auto', grid_area='check3'), style = {'description_width': descr}, disabled=False,) # Is probe window big enough to accomodate defocused probe?\n", "cl_check_1bf = Checkbox(value=True, description='Detector cover < 1α', indent=False,\n", - " layout=Layout(width='auto', grid_area='check4'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Detector cover is lower than 1 α\n", + " layout=Layout(width='auto', grid_area='check4'), style = {'description_width': descr}, disabled=False,) # Detector cover is lower than 1 α\n", "cl_check_nbf = Checkbox(value=True, description='Detector cover > 6α', indent=False,\n", - " layout=Layout(width='auto', grid_area='check5'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Detector cover is higher than 6 α\n", + " layout=Layout(width='auto', grid_area='check5'), style = {'description_width': descr}, disabled=False,) # Detector cover is higher than 6 α\n", "\n", - "footer = Button(description='Interactive results',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='whitesmoke'))\n", "\n", - "small_controls = GridBox(children=[header, method, graph_size, ctf_xaxis, cl_check_laa,cl_check_pix,cl_check_def, cl_check_1bf, cl_check_nbf, footer],\n", + "chart = ToggleButtons(options=[('STEM','micr'), ('Charts','dedic'),('Both','both')],value = 'both', description='Show',\n", + " layout=Layout(width='auto', grid_area='charts'), style = {'description_width': descr,'button_width': butt, 'font_weight': 'bold'}, button_style='',disabled=False,)\n", + "\n", + "\n", + "footer = Button(description='Interactive results',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='whitesmoke', font_weight= 'bold'))\n", + "\n", + "small_controls = GridBox(children=[header, method, chart, graph_size, ctf_xaxis, cl_check_laa,cl_check_pix,cl_check_def, cl_check_1bf, cl_check_nbf, footer],\n", " layout=Layout(\n", " width='100%',\n", " grid_template_rows='auto auto auto',\n", - " grid_template_columns='40% 20% 20% 20%',\n", + " grid_template_columns='50% 20% 30%',\n", " grid_template_areas='''\n", - " \"header header header header \"\n", - " \"methods ctxaxis check1 check2\"\n", - " \"gr_size ctxaxis check3 check5\"\n", - " \". ctxaxis check4 . \" \n", - " \"footer footer footer footer\"\n", + " \"header header header \"\n", + " \"methods ctxaxis check1 \"\n", + " \"charts ctxaxis check2 \"\n", + " \"gr_size ctxaxis check3 \"\n", + " \". . check4 \"\n", + " \". . check5 \"\n", + " \"footer footer footer\"\n", " ''')\n", " )\n", "\n", "\n", - "def ptycho_interact(beam, aperture, aperture_res, aperture_res2, probe, cl, matrix, defocus, mag, camera, binning, dwell_time, restriction, method, ctf_xaxis, cl_check_laa, cl_check_pix, cl_check_def, cl_check_1bf, cl_check_nbf, graph_size):\n", + "def ptycho_interact(beam, aperture, aperture_res, aperture_res2, probe, cl, matrix, defocus, mag, camera, binning, dwell_time, restriction, method, ctf_xaxis, cl_check_laa, cl_check_pix, cl_check_def, cl_check_1bf, cl_check_nbf, graph_size, chart):\n", " \n", " ### SINGLE SETTING PARAMETERS ###\n", " wavelength = pty.get_wavelength(beam)*1e12\n", @@ -464,59 +468,90 @@ " pctf = pty.get_ssb_ctf()\n", " apertury = opt.apertures()\n", " \n", + "\n", + "\n", + "\n", + "\n", + "\n", " \n", " ### SAMPLE PLANE ##########################################\n", " points = 5\n", " yyy = np.append(np.linspace(0,overlap,100), 200)\n", - " wid = np.linspace(0,points-1,points) * step_size_corr\n", + " wid = np.array([-2, -1, 0, 1, 2]) * step_size_corr\n", " \n", " fig5 = make_subplots(specs=[[{\"secondary_y\": True}]])\n", - " fig5.update_xaxes(title_text=\"Sample plane (nm)\", range=[0, points*step_size_corr], showgrid=False, zeroline=False)\n", - " fig5.update_yaxes(range=[0, points*step_size_corr],showgrid=False, showticklabels=False)\n", + " fig5.update_xaxes(title_text=\"Sample plane (nm)\", range=[-4*step_size_corr, 4*step_size_corr], showgrid=False, zeroline=False)\n", + " fig5.update_yaxes(range=[-4*step_size_corr, 4*step_size_corr],showgrid=False, showticklabels=False)\n", " for x in range(len(wid)):\n", " fig5.add_trace(go.Scatter(showlegend=False, x=wid, y=np.ones(len(wid))*wid[x], marker_color='blue'),secondary_y=False)\n", " for x in range(len(wid)): \n", " for y in range(len(wid)): \n", " fig5.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", opacity=0.1, fillcolor=\"#FF7F7F\",name=False, x0=wid[x]-beam_diameter/2, y0=wid[y]-beam_diameter/2, x1=wid[x]+beam_diameter/2, y1=wid[y]+beam_diameter/2, line_color=\"red\",secondary_y=False)\n", "\n", - " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='blue', name='Step '+str(np.round(10*step_size,3))+' Å ('+str(np.round(10*step_size_corr,3))+')'))\n", + " fig5.add_annotation(x=-1*step_size_corr, y=-2*step_size_corr, ax=-2*step_size_corr, ay=-2*step_size_corr, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=-2*step_size_corr, y=-1*step_size_corr, ax=-2*step_size_corr, ay=-2*step_size_corr, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=-1.25*step_size_corr, y=-1.5*step_size_corr, text=str(np.round(10*step_size_corr,3))+' Å', showarrow=False, yshift=0)\n", + " \n", + " fig5.add_annotation(x=2*step_size_corr-0.75*beam_diameter/2, y=-2*step_size_corr-0.75*beam_diameter/2, ax=2*step_size_corr, ay=-2*step_size_corr,xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=2*step_size_corr+0.75*beam_diameter/2, y=-2*step_size_corr+0.75*beam_diameter/2, ax=2*step_size_corr, ay=-2*step_size_corr,xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=2.25*step_size_corr, y=-2*step_size_corr, text=str(np.round(10*beam_diameter,2))+' Å', showarrow=False, yshift=-15)\n", + " \n", + " \n", " if method == 'direct':\n", " xx = 2*wavelength/(np.sin(semi_angle_corr*omega/1000))/100\n", " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='green', name='Rec step '+str(np.round(xx[-1]/3,2))+' Å')) \n", " \n", - " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='#FF7F7F', name='Beam ⌀ '+ str(np.round(10*beam_diameter,2))+\" Å\"))\n", " fig5.update_layout(legend=dict(orientation=\"v\",yanchor=\"bottom\",y=0.85,xanchor=\"right\",x=0.9))\n", - " fig5.add_trace(go.Scatter(showlegend=False, x = points*step_size_corr*np.ones(len(yyy)), y=yyy, opacity=1,mode='markers',marker_symbol = \"triangle-down\",marker=dict(size=20, color=yyy, colorscale='rainbow_r', showscale=False)),secondary_y=True) \n", - " fig5.update_yaxes(title_text=\"Overlap %\", range=[0, 100],showgrid=False, showticklabels=True, tickvals = [0, 20, 40,50, 60,70, 80, 90, 100], secondary_y=True)\n", - " fig5.update_layout(title={'text': \"Sample\",'y':0.95, 'x':0.30,'xanchor': 'left','yanchor': 'top'})\n", + " fig5.add_trace(go.Scatter(showlegend=False, x = -4*step_size_corr*np.ones(len(yyy)), y=yyy, opacity=1,mode='markers',marker_symbol = \"triangle-down\",marker=dict(size=20, color=yyy, colorscale='rainbow_r', showscale=False)),secondary_y=True) \n", + " fig5.update_yaxes(range=[0, 100], showgrid=False, showticklabels=False, secondary_y=True)\n", + " fig5.update_yaxes(title_text=\"Overlap %\", showgrid=False, showticklabels=True, tickvals = [-4*step_size_corr, -2*step_size_corr, 0, 2*step_size_corr, 4*step_size_corr],ticktext = ['0', '25', '50', '75', '100'], secondary_y=False)\n", + " fig5.update_layout(title={'text': \"Sample\",'y': 0.975, 'x': 0.30,'xanchor': 'left','yanchor': 'top'})\n", " fig5.update_layout(width=4.5*graph_size, height=3.0*graph_size, margin =dict(l=1.3*graph_size, r=0.3*graph_size, t=0.1*graph_size, b=0))\n", "\n", " \n", + "\n", + " \n", + " \n", + " \n", + " \n", " \n", " \n", " \n", " ### DETECTOR PLANE ##########################################\n", - " fig6 = go.Figure() \n", - " fig6.update_xaxes(title_text=\"Pixels\", range=[1, num_pixels], zeroline=False)\n", - " fig6.update_yaxes(range=[1, num_pixels], showticklabels=False,)\n", - " for x in np.linspace(1,int(np.round(covered_alfas)),int(np.round(covered_alfas))):\n", - " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", x0=num_pixels/2-x*beam_diameter_pix, y0=num_pixels/2-x*beam_diameter_pix, x1=num_pixels/2+x*beam_diameter_pix, y1=num_pixels/2+x*beam_diameter_pix,line_color=\"gray\")\n", - " fig6.add_trace(go.Scatter(showlegend=False,x=[num_pixels/2],y=[num_pixels/2-x*beam_diameter_pix],mode=\"lines+text\", text=str(int(2*x))+\"α\",textposition=\"top center\"))\n", + " fig6 = make_subplots(specs=[[{\"secondary_y\": True}]])\n", " \n", - " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", fillcolor=\"#FF7F7F\", x0=num_pixels/2-beam_diameter_pix/2, y0=num_pixels/2-beam_diameter_pix/2, x1=num_pixels/2+beam_diameter_pix/2, y1=num_pixels/2+beam_diameter_pix/2,line_color=\"#FF7F7F\")\n", - " fig6.add_annotation(x=num_pixels/2-beam_diameter_pix/2, y=num_pixels/2, ax=num_pixels/2+beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=3,arrowsize=1,arrowwidth=1,arrowcolor='black')\n", - " fig6.add_annotation(x=num_pixels/2+beam_diameter_pix/2,y=num_pixels/2, ax=num_pixels/2-beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y', text='', showarrow=True,arrowhead=3,arrowsize=1,arrowwidth=1,arrowcolor='black')\n", + " for x in np.linspace(1,int(np.round(covered_alfas)),int(np.round(covered_alfas))):\n", + " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", x0=num_pixels/2-x*beam_diameter_pix, y0=num_pixels/2-x*beam_diameter_pix, x1=num_pixels/2+x*beam_diameter_pix, y1=num_pixels/2+x*beam_diameter_pix, line=dict(color=\"gray\",width=1,dash=\"dot\",),secondary_y=False)\n", + " fig6.add_trace(go.Scatter(showlegend=False,x=[num_pixels/2],y=[num_pixels/2+x*beam_diameter_pix],mode=\"lines+text\", text=str(int(2*x))+\"α\",textposition=\"top center\"),secondary_y=False)\n", + " \n", + " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\",opacity=0.5, fillcolor=\"#FF7F7F\", x0=num_pixels/2-beam_diameter_pix/2, y0=num_pixels/2-beam_diameter_pix/2, x1=num_pixels/2+beam_diameter_pix/2, y1=num_pixels/2+beam_diameter_pix/2,line_color=\"#FF7F7F\",secondary_y=False)\n", + " fig6.add_annotation(x=num_pixels/2-beam_diameter_pix/2, y=num_pixels/2, ax=num_pixels/2+beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig6.add_annotation(x=num_pixels/2+beam_diameter_pix/2,y=num_pixels/2, ax=num_pixels/2-beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y', text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", " fig6.add_annotation(x=num_pixels/2, y=num_pixels/2, text=str(np.round(beam_diameter_pix,1))+\" pix\", showarrow=False, yshift=10)\n", - " fig6.add_trace(go.Scatter(x=[1,1.5,1,1.5], y=[1,1,1.5,1.5], name='Pix ∠ ' +str(np.round(pixel_angle,2))+' mrad'))\n", - " fig6.add_trace(go.Scatter(x=[1,1.5,1,1.5], y=[1,1,1.5,1.5], name='Pixel size ' +str(np.round(size_pixel,2))+' um'))\n", - " fig6.update_layout(legend=dict(orientation=\"v\",yanchor=\"bottom\",y=0.80,xanchor=\"right\",x=1))\n", - " fig6.update_layout(title={'text': \"Detector\",'y':0.95, 'x':0.11,'xanchor': 'left','yanchor': 'top'})\n", - " fig6.update_layout(xaxis = dict(tickmode = 'array', tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels])) \n", - " fig6.update_layout(width=3.5*graph_size, height=3*graph_size, margin = dict(l=0.35*graph_size, r=0.65*graph_size, t=0.3*graph_size, b=0))\n", + " fig6.add_annotation(x=num_pixels/2, y=num_pixels/2, text=\"1 pix = \"+str(np.round(pixel_angle,2))+\" mrad\", showarrow=False, yshift=-10)\n", + " fig6.add_annotation(x=num_pixels/2, y=num_pixels/10, text=\"Real pixel size = \"+str(np.round(size_pixel,2))+\" μm\", showarrow=False, yshift=0)\n", + " fig6.add_trace(go.Scatter(showlegend=True,x=[-1], y=[-1], marker_color='red', name = \"BF disk\"),secondary_y=False)\n", + " \n", + " fig6.add_trace(go.Scatter(showlegend=False,x=[-1], y=[-1]),secondary_y=True)\n", + " \n", + " fig6.update_layout(legend=dict(orientation=\"h\",yanchor=\"bottom\",y=1.00,xanchor=\"right\",x=0.9))\n", + " fig6.update_layout(title={'text': \"Detector\",'y':0.975, 'x':0.11,'xanchor': 'left','yanchor': 'top'})\n", + " \n", + " fig6.update_xaxes(title_text=\"Pixels\", range=[1, num_pixels], zeroline=False, tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels]) \n", + " fig6.update_yaxes(range=[1, num_pixels], showticklabels=False, showgrid=False, secondary_y=False)\n", + " fig6.update_yaxes(title_text=\"Angular cover (mrad)\", range=[1, num_pixels], showticklabels=True, tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels], ticktext = ['', np.round(detector_cover/4,1), np.round(detector_cover/2,1), np.round(3*detector_cover/4,1), np.round(detector_cover,1)], secondary_y=True) \n", + " \n", + " fig6.update_layout(width=3.5*graph_size, height=3*graph_size, margin = dict(l=0.35*graph_size, r=0.48*graph_size, t=0.3*graph_size, b=0))\n", "\n", " \n", " \n", " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " focusdepth = wavelength/(semi_angle_corr/1000)**2/1000\n", " ssb_cl, kolik = pty.cl4ssb(detector_cover_a_all)\n", " ### MICROSCOPE SCHEME ##########################################\n", @@ -654,20 +689,27 @@ " fig4.update_layout(xaxis = dict(tickmode = 'array', tickvals = opt.cameralengths()))\n", " fig4.update_layout(width=8.0*graph_size, height=2.25*graph_size, margin =dict(l=1.3*graph_size, r=0.25*graph_size, t=0.3*graph_size, b=0.3*graph_size))\n", "\n", - " ### SHOWING ###\n", - " right_column = VBox([go.FigureWidget(fig),HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", - " left_column = VBox([go.FigureWidget(fig8)])\n", - " \n", - "\n", + " \n", " \n", + " ### SHOWING ###\n", + " if chart == 'micr':\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column]) \n", " \n", - " total = HBox([left_column, right_column])\n", - " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px')\n", - " \n", + " elif chart == 'dedic':\n", + " right_column = VBox([go.FigureWidget(fig) ,HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", + " total = HBox([right_column])\n", + " \n", + " else: \n", + " right_column = VBox([go.FigureWidget(fig) ,HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column, right_column])\n", + " \n", + " \n", + " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px') \n", " display(total)\n", - " \n", - " \n", - " \n", + " \n", + "\n", " \n", " \n", " if method == 'iterative': \n", @@ -781,7 +823,7 @@ " \n", " dictionary[i] = ToggleButtons(options=data_tuple, button_style=button_style, icons = icons, **align6)\n", " \n", - " legend = ToggleButtons(options=['Ptycho pix. (Å)','Det. cover (mrad)','Det. cover (α)', 'Max. def. (nm)','Selected by CL ','Nominal CL (cm)',], \n", + " legend = ToggleButtons(options=['Ptycho pixel (Å)','Det. cover (mrad)','Det. cover (α)', 'Max. def. (nm)','Selected by CL ','Nominal CL (cm)',], \n", " tooltips=['Dependent on: detector, beam energy, camera length',\n", " 'Dependent on: detector, camera length',\n", " 'Dependent on: detector, camera length, probe semi-angle',\n", @@ -796,22 +838,32 @@ " \n", " \n", " \n", + " \n", " ### SHOWING ###\n", - " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", - " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", - " left_column = VBox([go.FigureWidget(fig8)]) \n", - " \n", - " \n", - " total = HBox([left_column, right_column])\n", - " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px')\n", - " \n", + " if chart == 'micr':\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column]) \n", + " \n", + " elif chart == 'dedic':\n", + " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", + " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", + " total = HBox([right_column])\n", + " \n", + " else: \n", + " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", + " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", + " left_column = VBox([go.FigureWidget(fig8)]) \n", + " total = HBox([left_column, right_column]) \n", + " \n", + " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px') \n", " display(total)\n", + "\n", " \n", " return \n", " \n", "gui = interactive_output(ptycho_interact, {\"beam\": beam, \"aperture\": aperture,\"aperture_res\": aperture_res,\"aperture_res2\": aperture_res2, \"probe\": probe, \"cl\": cl, \"matrix\": matrix, \"defocus\": defocus, \"mag\": mag, \"camera\": camera,\n", " \"binning\": binning, \"dwell_time\": dwell_time, \"restriction\": restriction, \"method\": method, \"ctf_xaxis\": ctf_xaxis,\n", - " \"cl_check_laa\":cl_check_laa, \"cl_check_pix\": cl_check_pix, \"cl_check_def\": cl_check_def, \"cl_check_1bf\": cl_check_1bf, \"cl_check_nbf\": cl_check_nbf, \"graph_size\": graph_size}) \n", + " \"cl_check_laa\":cl_check_laa, \"cl_check_pix\": cl_check_pix, \"cl_check_def\": cl_check_def, \"cl_check_1bf\": cl_check_1bf, \"cl_check_nbf\": cl_check_nbf, \"graph_size\": graph_size,\"chart\": chart}) \n", "\n" ] }, @@ -845,14 +897,14 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 390, "id": "7937f054-fcd0-4e67-a20f-7696f5903a94", "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "47642ef0d8b94bceb11f878293997c33", + "model_id": "6f84c5e6c29142a9bb9db08cd48aefbd", "version_major": 2, "version_minor": 0 }, @@ -865,6 +917,7 @@ } ], "source": [ + "\n", "display(HBox([controls, VBox([small_controls, gui]) ]))" ] }, @@ -888,6 +941,139 @@ "# import ipywidgets\n", "# print(\"Ipywidgets version \", ipywidgets.__version__)" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "33727cb4-c3d7-4134-be89-f10b9922fbda", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 75, + "id": "1cb358ba-42c4-420a-a0a6-0de425e670a6", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Ipywidgets version 7.6.5\n" + ] + } + ], + "source": [ + "import ipywidgets\n", + "print(\"Ipywidgets version \", ipywidgets.__version__)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a52ff861-8243-4d39-b7da-1ee8fb0ee494", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "61970cec-8678-4edf-99f6-4f832d2aad8d", + "metadata": {}, + "outputs": [], + "source": [ + "https://ipywidgets.readthedocs.io/en/latest/user_install.html" + ] + }, + { + "cell_type": "code", + "execution_count": 386, + "id": "a4628def-4c43-456e-a93b-95952e2a3dea", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "" + ], + "text/plain": [ + "" + ] + }, + "execution_count": 386, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "from collections import defaultdict\n", + "from IPython.display import HTML\n", + "import ipywidgets\n", + "from pprint import pprint\n", + "\n", + "reverse_lut = defaultdict(set)\n", + "styles = set()\n", + "for export_name in dir(ipywidgets.widgets):\n", + " export = getattr(ipywidgets.widgets, export_name)\n", + " try:\n", + " if issubclass(export, ipywidgets.Widget) and 'style' in export.class_trait_names():\n", + " reverse_lut[export.style.klass.__name__].add(export.__name__)\n", + " styles.add(export.style.klass)\n", + " except TypeError:\n", + " pass\n", + "\n", + "html = '\"\n", + "HTML(html)" + ] + }, + { + "cell_type": "code", + "execution_count": 74, + "id": "54c1add7-f5ec-43f0-a447-e22dd03926a7", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "
AttributeToggleButtonsStyleButtonStyleDescriptionStyleProgressStyleSliderStyle
description_width✓✓✓✓
button_width✓
font_weight✓✓
bar_color✓
handle_color✓
button_color✓
" + ], + "text/plain": [ + "" + ] + }, + "execution_count": 74, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "attributes = defaultdict(set)\n", + "base_traits = set(ipywidgets.Style.class_trait_names())\n", + "\n", + "for s in styles:\n", + " for t in s.class_trait_names():\n", + " if not t.startswith(\"_\") and t not in base_traits:\n", + " attributes[s.__name__].add(t)\n", + "all_attributes = set().union(*attributes.values())\n", + "\n", + "html = '\\n'\n", + "html = f\"{html}{ ''.join(f'' for s in attributes.keys()) }\"\n", + "for a in all_attributes:\n", + " html = f\"\"\"{html}{ ''.join(f'' for attribs in attributes.values()) }\"\"\"\n", + "html += \"
Attribute{s}
{a}{\"✓\" if a in attribs else \"\"}
\"\n", + "HTML(html)" + ] } ], "metadata": { diff --git a/calibrations.xlsx b/calibrations.xlsx index db7ef42ab14f7cbabd0732d0dc783a0672c62a0d..10025df9cb36bf46f1519252a505c2f808e5d39a 100644 GIT binary patch delta 4376 zcmZWscQ~BQ7GJ%F=yjD4L``&}6KxYBdKcltCwlbq>OI;UQNofS1ZzdNI#Hv1h`M@B zq7xRYTsQZ*-@VVh=b!n_oH}RbnK|b;5;rdjw-QSTIB!5=D)1nXa1scF5(0tvy%h0v z^K`a#b8{B``VgKtmlWaQYYzm_wOQ2>l#iwAScT_Ot+%X zS!#GpOl4a7E_(Sp)5pfy`+gV96m=G3vG{_eqG<`yksvKuiZxFSpLv}iDq6TJ$dbxIBP8`(nk>iO z$Vp8ay6ed;KQ`hwOsCDbllY0l!`l1#faF>OMDL)K(Lg-Pa$FgzI!%$;;aS^d9si_;dQ`bq>VRHJ!zA9TZ(&NAP-Z>wi#XX8tsEjk$-@ zJ?9)OMP(6oP1(Wdc}>u!d^dmL;0nlBP%72ilt^xfu>wwcJX2Zj44sTBRp&f6-jbc=v9H}-Jk@K=0<_z1fU`x|%k-j?Z zJ}%vov&6-tkphR$0!w4at9bz`JAs~}5~TYpK07+o=sjk6m{?~L!gL<~$s3XM=@z1$ z+(nAKQ+Y@9?ELU5na3*61{o}(oldf0+v8CaUH{{aAW=i)3q4~MqLyM?A|GFyS4tlw zRJevzW_5);(+eSLw zcDQ~4PaO5#WhNyj%-_vj;(-r*XsgFa8EHvgUTz%SB7;EaL1lhs@D@7}fKEbl+>YK8 zIED)DLW2f7;;Xf*U!$Z9YvnCzf}4rN?*P2&bnx|0_HLL)F!>N@mljLS~@oUT`QMOXUgw#T_6^rk&|x=Jt_qF zYMXJow>5qLT6v6kZA%dY`~h=u2sP~Dx$8vz=tXI!V}I7u%s@s|1wo9)P#^)~p(?rN z#(e-CJ8b0`zZjuExP30XX679`O?+?8DKt&#V8sdjXO|@PmO>BHYRBjjXaf@w?ZU9! zz1uKqH02s;7chKQcfpja*MlrXpYL^kVUnn1xPn+e?MFY0KN1iC{KJaeB3~o#UpGA- zO(^9C4k>43;UCA{Sw&Xwf1N6G&ssCh3`qoEidxzCi(kU%kF=Lau%)P<&t~0lC7d{>h<3LZSubuDQNB z<)1gF%n6`f(f7%t_lUNk`18G0xrLUA(%qFGocE&Z&=qu9Nm1XIv^**JPIe2nSxNC3 z!UaMwk}k*ufLQOPTBnx1_3+oTiQJ}{2Rbc61f`i{k~I|X!!XGBv(c3O5- zdJ7TEnU6Z%e;z9RRLuQ% zr2~1-S1CJPbt+GrEJpiF3*CM!GcP|(Pz_*93Iu-GT4|NL55`NEz|JKE;Jq{h^hd2> zAzIVtzFyR)hXk~-5WX}$YOI&jjx4jn^0Ro&MHV$~24&YJ7CytVJ(+B*VvBZ80scd$ z<~;UewsZK^vL=DVo`GK-S|g>5Cx4$-40w%rW`GUc#|@+{0)v=8b@Hn%Y1Smh+zL$G zi?*?+%l@7__s5TyK?$^@XEw)sT|FJa?q9~pll&y>i6VXh@N(lTA_(OA`sNqDzJ7zA zah;QR6g`ROT4n>rgY-R`(pFGV~?fpS9;X{BWOW8Vule&_`ZGhM{Hh%@(bj^1?zDe13AqAYBoULQ_ zbfu<7J9f@iHxZ|K=Evht-pPs`4(%s9&;dc$t7og&5sfFPsngxYo(zc@tJOGvTvJf7 zl*InYK?zTq#CilX!tRMT(jNeV;?Hjf2ZDrSDzciCOdMQx;e#mzRq*weehgbxEA63? zOB8E^QD5Ne%@$!C=)+mkQ$`HNW8GwN-y%_9M62AU+s{$gdbg`}q1gA1^Ns=LAQ0Xg zSuDpXR2cCq1xHkcr;pr}W6M+xzEC7r|EIu6A1Y<6WkX^_WTXOA17?q|Q?m`)UQRV6 zv0C5NdTIG{=Acw)T?FYHX{{I>e71hFyY)4r*nA_HuReTnP=3uPPpoO=hin*WOP_Bt z0QbJUo>{|2@|bw-9SpK6=SQsVyqpbx^tNn;l+m~t#=D}SCbxoYHs&{#xWjkjzx zv^69-Bb#3YwY9m+UB5e6aqEFSsDNUxswnf0N`CyRO&jc;G&f8eH2DJo9@?lDXZ#07~YT}qF-vIf9go{{S9$Xb;|NNc^? zCe*v_7uRBbdwrwhg>`>#jJ#!kuYgX!s+lIy4^BM0$j2#6RZ5f9v1?thez%E@NxnLd zN%dvEnF7N}Ne-HVjV5R}fBXf9kGx2%+(SP-m)*wq$S&zvy!B7g2MBXD&eEYSpeGU< z`FWEFE~yI*Br}!t`{SF}F7@ze2tewGSg9c~jw`;!;}dF>=<{@IoJL$wDLgZA5~ zS(vSqQP?IL3tAXpb*Kr?YfQ@R4f`6Q7zvMcRH}B+j`?08?F6OpS@1>9#X-FO;lDBt z{3qpi{JoB?q;Xs)j1DJzgygkP*5ubsRmrf0J@KFE58n>2RR{YWTi_p)@} zRTe-_b2=$>fc-4TuV02NsFmTp%_yo>?6+N5VH#_ccAL$%}84;MXVtZo&j&R z`(xPB7*?&mmtbG`HHo47frLNlq~XP#?Uc!y2&Cmrz;qgVnQEVrcc9THFPh)%y-EI` zi2pladz6^*tb%!CdW*ME5mKLRliN#F%zq}($CN5NVhRm+LGPaVTvW~dVzf{^e6dsd zqRpW`JPFevqHr|P{pUfsh8p11!c-Aju+)sdVCCgX3bj7wtu6nlDeuJbx-d>0q#xtg ztcXTq@V=RZ2l*Er4BI?acj=07>j*F6RD^-lBkPxB2= z^Nml7hg527bA`*J3t;weQ4GA&yF}Z%xOtIfumakoZ2CTU7r4 z@c-7}O3ghx{{C?Luru3I<_aCrA;&Z5)49utp4Y$%x)1Ui;jtvTSsw?v*(9VI{`bBX zSn58G08ta}zhZ~Wa1=SDz1KtJh^XK+wsYfgjbhyI8y}v=InmEBO#hdx^TJ0NBtadMG>sz(X!j_++STTs@XB3Ls8M&?XQ#rmk+sB_v zxsD6N#M5}j9-C>Fab)dI(xOs?QQDtCS5a1y>gjWFhAO23UU1|AJIUsn3cX`XA?X`X zTNFxRZ_$T4l=2D-joIPD9)X&o97H!OJ8=eKpu0K)WmF9*wdr@~%bY1s530gKaP$Ey z@f7MFv%}Gwf9kUkt(b|ZuKK}2H+4qx8{srli!E+P)?|IqMGQ*l%k zH>GJX-cm6mh!bpf-rSw?^xp+J#i4`^KD#G=zahyDl$WqK8(_%Lb6*U4E zL|Mpcz3R^jxa@ESOwz$8VvHmwf(zFd0@K$x)wSYYDg-TRnHI2WAzyJEz`iP~0xw-& z=C+hx5Ly4c(PSZK|MvOblO}`??`dX^S6hykmcqjl{&W1T0W)2yC&cke31)8t#2z8b z{N>bBk1SwT6U-k=2x9)|-gzQSPE#A@V9rw6iqAlh=JPo7-fVyyV~@^DWseN}7>X5C zT|CykhUI0$^uevAHz^Ox+X1hB`CjOKNj$3DlxO6c)%AvyDY&Wxq&2>oJhVWdyvh0% z>-TiBV*X%fm16xlSRJi~hdaB44CA|xTG{2vREPQ(t$t=FhO5tfo2WzMEl7C;-L28K zjb)gRudNyc4UGeU!3ja+B$g|mr#yY~f9%kHsj{&0`w^$cGV2OQ^75>i3qfZ%sq@X6l~ zkL2%mALI|#!RZ72gEesYz>Z)6*1sDDWe9}vrUdy9EwO=0x9C7<2p#L+_2&Nxaq<5a zdJ!UmYa;m@+DHxRyXC)lXO_3g|IL9wXsvD(b#WQ2!O>8C+}S^F25wM0OawRA v^*7856Tz_oN5Z&ReqSFP2*e2&0-^ekfg7X>-U?U1=>{K%a}jm~{gd?%hBaL7 delta 4346 zcmZWtXH=6*yQNBR(pv!OphOS^^#}q=Po#sW)PNvJkxmG&gf88eURA2n1R{!blrBY( zUP1{1(pw^gayfeMJ?nhm{Mu`twfFPP?3pzsH;imQlg)E@+1Aa?S;p^)OR0&oTb>fzapT|#)fZpdAITI7{HoNGI!@x*MRa4{#CfViT{aMGbxevX~ zN0+u!V40jXo3BH|;$N5Z8Pxf#WlDvnMh22Rlve%x z(iaB()w1R2Ux}}Ty;9J<_Kf9rRo?KK;?g-)c3+j^?wjtlnsSJvtGr~+3YOb~{sIV3 z#tFUs4IkAjJ5v^BE%pbYp7-A9naDkevQeYdK7347=mLwQV`i`Dd0EK4mLXF`Z_;#Q z%>Szm*?R(TMb$5ulKY+7vtge>VAPa>AZTxGYp-@vE-CKr)2znK7^8&Ws^_swlv3zF z{Mnp|DidtI_ukgK7PWW<#tY=m06`hYciO1?9~F$XOP1K0x;_$*`;v-3Wy<+dsTalQ zCO>oJD@!1>K*GkJax$YarVWO(FCZg8C zxc%vOSfV6(%KBRR3vX{sNS0J4r-g{iw3z}zO@DHC<S5~tA(numoi z$1uQy+DJQ9Y$sd45T_i4(a*J#M~AJZRulp-#z%*Htido|%2&N6l02(eCQZBLk*b#N zP~AM_4$YB}s$8Sd8U?4WNPq76c%yk?dIR{!{|NJ|2#h4XAOTzoc1RchhKf zBIOOk+5q46Y#SwdJ}0>Y{LV={WQNb$?bISk%IIrb(t;l)bcp62#YeQ!H)A!bh9X>IEsQW*OVl#A1@p(6AEAkhUiwO0*S{H*(gWKqj1$9|bu=6b7|6JHMy+^<;|J{k z^9LBZ8^~-0X1DX*^_v?f>|1oVa}%XozU?MAw+t=pZR@N#DxpwO<&?yD)n%)Af;j6b z(d*ZfQ|e zN9d1KDe|Q%>trki3ne9ZZztUn)|1TiRt1h9oU_x(T7q>YJH33V>AvTRKN#|ITWcRC z$q{_zL0Q(|wtLy@^|{{ZW(+Vet7l*O1ER8i-K&h!60dC$=E26G-fORZH*{*I^|Fa5 zGZ<3OG>@#4wL6vmpNFZ-{h;{^>r|v0wb}Nfd%h2$#rW5 zsiMtde4nPQbW|i^*hZ9$4oF{oHtHQ+zuT(a*%szCeG21B3Fq*snrsFXjgKl{nETkY zxxK@MyqlvOIhM(Ev%fbmx~BUHv4G&CAJ3w6s(C=MtE3x(la3zx=APKTfP-ekbfa!e zdEdEcOqpGMBx^691w)eG)E%qymE8{9r=Sg`VzWs)lGRY5DW63y-3qf~38?d2?@oiI zo;;%yW_yaV5dRXEdc)lMS?l613(qM*JK%d)IOjBByEfc@b5|fGNq@U+C*I+d0rv?n zk!3^3M(ynQF~MKD8Tc_j;y2+6qYgZOsg69dae*pup>nckaFk0dge6+6i(1xpwf-CQ zRQ&L5K(h)-^@1RrhW7gBTgd)1)5rdpe&vtp51l73vxRR`5t-Bol`i5=oOVS?$yGD4 zRj#AY1NowtKZ8lo_7h^nBL8Ny z4-So<&8wFA>^(z!aZWM}{M2D}VK$i>w0>2#B`u+(^S8F2?GsIXg7Td_=-{3o#(h|A zydxr$jpg`|5~3DtP-i=SE6FutLhnaGstDh>elNWtps72*z>&YQ70Z;Ry1}iuh3*DaObxH9RPWHa9IB2@b8`Tx(lu4(%ei}AF!#Rs1j}$)e5y>3 z#XMt2d1Of8>-c-y*Ru_3E+WZ~lim?d%e40vx+h@t_0~CTED*NHW%^c+i54?`TJ0;C zy-m^y^6~VK7on?(CsQ|m)d?%e2bQ#!OatV#uH=x^Yk2lpgYQ6x_!b#|R1_Ep>0vV! zBGY_r|ChzGM9}h5L#aCNa2F~?rjVe;a$lhI#G5t?7s9M)%RH~mm%M`18F_RyG|)_v~Qcy{Cp~dNc_(8Z$1O$L{uHRa8)i_>Lx(KxBP!Pedw!Iyr7T*~rt85je?1_gg8S z>|l=$7JGI#u>#P>qR%};MbP`(Dj$A!MAvxx;M#WL#ghXAer@lq%=K7$HVjjoREi}i zjF+`X`m5aou_aLmm`8oz|e+-fGA9F;i} zo&DaZ0XRVB)bZ_vViFh!)z5WOrwd*|1ea>E_-7dTyAH#1Dt84Rx+zrn2cdQv-}IYk>(jqF@R^{t`c zSUMBLd@~l{UOG6f?s8+DuwI|BTjR4S{=B!IS{-u|f9*p5sP?kBo4{vZ{M*sZxoz7m z7(_V&|9UK6EHL-8Ki}J9&!Yz@h5TSl0`TubS7nvnLq0rybn{%_F3>Q%MvUG(*yx_@ za3j)n>^5;jPaNvm7A8Cisi0t- zc`%DL%G$V0&E9!G?U{E=Q#&oRyI6!F>(XoHVZ6QFXU4oQ6tBMF%N8J{x`e!sK+EPr z#H`SNE`jIN%F66Vrt1`&NA=~%t{^gu3a}I^(VmQoL0M(WZ?m#zG~dj{y}rIRDDgd_ zZT@13THpFXMU*r)-h^kl`9P$gxN*1hOd9v!(podh4dSlP?!O1?J>7x24bcy{*qy9x zy_k|Q-7lVNS$duZ{79SENA5TFL_Hh>vFUXQo=BUEHAqE&1C;my>L}r_|DO4-qM{er zx7oPtyd7_7t*8LE2(%rBdqWuAs<~3v75$h!puJdv{`F16rG6RQRF2NPlFovNeq|hR zfuLa7eYbn`NI6IgGI+~Ei)uxP91*IQz+H7~Brg(|yxy1PM08n%3VFQhKkbj(tDBsf~FFWvmWjse;y=c|EF=9k+2?VNlJ2vTGQ4rw|91C*Ab=J{5$XaK?`UTmkSqV$#OetT#ee zA^lJ9K_sF33&!aZI>!=`>b{EZ4x^S&kp!wlfw8c;i@k&C$9<1)eJBeeX4pk$~ep z4BfQ+T+o4uc+C16iAn#bZe&#t4)7)mM5naJdBXNnJFV#t{aYrZTN=gIJX+Ve>ZUab zGFw-%9&Xkli6)oc2=C_y%fApC)U(hk5v)~of#tH6*le=!n~|TOX9tpr+idNAcOyOW zjOwoBCr4PN6bNMgo9W51<9$-gyLTw$+RSIA$Nc!2_P0yMclZgV*9Jz$pr8#Fi> z&qk5nHgfG+BBce?j1qL0=R04yF9$~iJ9K&z6&WW@@Jk!npr#x@wMbnHba<^|abD;nh`qv3m9<*vIMrJLCTeog3zeY}97_S*JgL8=!FmlY`rm14 z;pleJL!GlE{O2~rDm5mkWfES>^i4y-;ZKlkBQ|zD0E-!6NS1j(#!2aZ#dE)mc3aZFcQY-mx)=Zld zP@tPUgmz=*$!Vf;6^VT*XrRF%`+v!Zxp4`U5_TqVDJ{`A7l!l0zyI5qPxjkeyFkkc zQ}v9zryAwP_nkDIYExf%Qo!z~&SY?R4qs%ZL}M7vRoZzl0pW8@C|ox!VGY#ltF`K@ z70DK`=MF?y-YJl~EsQ49uRk9);jrx8K)tDN&2NNGPDax!SbzJgsZpmE?q@$cZp;%S z-Fz-5-=)<%+^27a!mh_xFg>>y14X_{LIo^yw7Ugbip`qRnt zUGlGdh3B~wihj*|;;g_Ha(3`Onz4G$|1mku|9=92#wa-H7cdW{NhZPoHY>A(q(SVw z|2!s0NVv}KB>%K2P&P=K@1J|;f5xSxev8BgLCK(3f5V|Oh*UmIaK}0>MzyO7L#5 qF`1X^?~bms4n^=FSoqK40EN)(q&%P^TA4Hx^g#=sOAY+H?*9R&a9U^p diff --git a/ptychoScopy.ipynb b/ptychoScopy.ipynb index 268e78d..1456a12 100644 --- a/ptychoScopy.ipynb +++ b/ptychoScopy.ipynb @@ -2,19 +2,10 @@ "cells": [ { "cell_type": "code", - "execution_count": 1, + "execution_count": 79, "id": "7425242d-3c91-4c1e-a424-08625a38ee7a", "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "C:\\Users\\skoupy_r\\Anaconda3\\lib\\site-packages\\scipy\\__init__.py:155: UserWarning: A NumPy version >=1.18.5 and <1.25.0 is required for this version of SciPy (detected version 1.26.1\n", - " warnings.warn(f\"A NumPy version >={np_minversion} and <{np_maxversion}\"\n" - ] - } - ], + "outputs": [], "source": [ "### Initial packages import ###################################################\n", "import numpy as np\n", @@ -25,7 +16,7 @@ "from pandas import read_excel\n", "from IPython.display import display\n", "from plotly.subplots import make_subplots\n", - "from ipywidgets import interactive_output, HBox, VBox, Layout, Label, Valid, ToggleButtons, RadioButtons, Dropdown, IntSlider, Checkbox, Button, ButtonStyle, GridBox\n", + "from ipywidgets import interactive_output, HBox, VBox, Layout, Label, Valid, ToggleButtons, RadioButtons, Dropdown, IntSlider, Checkbox, Button, ButtonStyle, GridBox,IntSlider\n", "\n", "# print(\"Numpy version \", np.__version__)\n", "# print(\"Pandas version \", pd.__version__)\n", @@ -300,31 +291,37 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 389, "id": "8055b802-cf83-4250-aea9-54e3e6b73db0", "metadata": {}, "outputs": [], "source": [ "### CONTROLS ##########################################\n", "\n", - "ali = dict(style = {'description_width': '60px','button_width': '33px'}, button_style='',disabled=False)\n", - "ali2 = dict(style = {'description_width': '60px','button_width': '90px'}, button_style='',disabled=False)\n", - "header = Button(description='STEM control panel',layout=Layout(width='auto', grid_area='header'), style=ButtonStyle(button_color='silver'))\n", + "ali = dict(style = {'description_width': '60px','button_width': '35px', 'font_weight': 'bold'}, button_style='',disabled = False)\n", + "ali2 = dict(style = {'description_width': '60px','button_width': '90px', 'font_weight': 'bold'}, button_style='',disabled = False)\n", "\n", - "beam_set = Button(description='Electron beam settings',layout=Layout(width='auto', grid_area='beam_set'), style=ButtonStyle(button_color='lightsteelblue')) \n", + "\n", + "\n", + "\n", + "\n", + "\n", + "header = Button(description='STEM control panel',layout=Layout(width='auto', grid_area='header'), style=ButtonStyle(button_color='silver', font_weight= 'bold',))\n", + "\n", + "beam_set = Button(description='Electron beam settings',layout=Layout(width='auto', grid_area='beam_set'), style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold')) \n", "beam = ToggleButtons(options=opt.energies(), value=200,description='E (keV)', layout=Layout(width='auto', grid_area='beam_set1'), **ali)\n", "aperture = ToggleButtons(options=opt.apertures(), value='30um', description='Aperture', layout=Layout(width='auto', grid_area='beam_set2'), **ali)\n", "probe = ToggleButtons(options=opt.probes(), value='8C', description='Probe', layout=Layout(width='auto', grid_area='beam_set3'), **ali)\n", "defocus = IntSlider(description='Δf (nm)', value=0, min=0, max=1000, step=1, continuous_update=False, layout=Layout(width='auto', grid_area='beam_set4'), **ali)\n", "\n", - "scanning_set = Button(description='Scanning parameters',layout=Layout(width='auto', grid_area='scanning_set'),style=ButtonStyle(button_color='lightsteelblue'))\n", + "scanning_set = Button(description='Scanning parameters',layout=Layout(width='auto', grid_area='scanning_set'),style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold'))\n", "mag = ToggleButtons(options=opt.magnifications(),value=3,description='Mag. Mx', layout=Layout(width='auto', grid_area='scanning_set1'), **ali)\n", "matrix = ToggleButtons(options=opt.mappings(), value=2048, description='Matrix', layout=Layout(width='auto', grid_area='scanning_set2'), **ali)\n", "dwell_time = ToggleButtons(options=opt.dwelltimes(), value=10, description='t (μs)', layout=Layout(width='auto', grid_area='scanning_set3'), **ali)\n", "\n", "\n", "\n", - "camera_set = Button(description='Detection',layout=Layout(width='auto', grid_area='camera_set'),style=ButtonStyle(button_color='lightsteelblue'))\n", + "camera_set = Button(description='Detection',layout=Layout(width='auto', grid_area='camera_set'),style=ButtonStyle(button_color='lightsteelblue', font_weight= 'bold'))\n", "cl = ToggleButtons(options=opt.cameralengths(), value=12, description='CL (cm)', layout=Layout(width='auto', grid_area='cl_set1'), **ali)\n", "camera = ToggleButtons(options=opt.detectors(), description='Detector', layout=Layout(width='auto', grid_area='camera_set1'), **ali2)\n", "restriction = ToggleButtons(options=[('.',False), ('..',True)], value = False, description='PAAR',icons = ['times','check'], layout=Layout(width='auto', grid_area='camera_set2'), **ali)\n", @@ -348,7 +345,7 @@ "footer = Button(description='Footer',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='silver'))\n", "\n", "\n", - "controls = GridBox(children=[header, beam_set, beam, aperture, probe, defocus, scanning_set, mag,matrix,dwell_time, cl,camera_set, camera, restriction, binning, beam_res, aperture_res,aperture_res2,fov_res, dwell_time_res,footer],\n", + "controls = GridBox(children=[header, beam_set, beam, aperture, probe, defocus, scanning_set, mag,matrix,dwell_time, cl,camera_set, camera, restriction, binning, beam_res, aperture_res,aperture_res2,fov_res, dwell_time_res],\n", " layout=Layout(\n", " width='22%',\n", " grid_template_rows='auto auto auto',\n", @@ -375,49 +372,56 @@ "\n", "# SMALL CONTROLS\n", "\n", - "descr = '50px'\n", - "butt = '30px'\n", + "descr = '80px'\n", + "butt = '60px'\n", "\n", - "header = Button(description='Graph controls',layout=Layout(width='auto', grid_area='header'),style=ButtonStyle(button_color='whitesmoke'))\n", + "header = Button(description='Graph controls',layout=Layout(width='auto', grid_area='header'),style=ButtonStyle(button_color='whitesmoke', font_weight= 'bold'))\n", "\n", - "method = ToggleButtons(options=[('SSB','direct'), ('Iter','iterative')], description='Method', tooltips=['SSB', 'PIE, MLc, DM'],\n", - " layout=Layout(width='auto', grid_area='methods'), style = {'description_width': descr,'button_width': butt}, button_style='',disabled=False,)\n", + "method = ToggleButtons(options=[('SSB','direct'), ('Iterative','iterative')], description='Method', tooltips=['SSB', 'PIE, MLc, DM'],\n", + " layout=Layout(width='auto', grid_area='methods'), style = {'description_width': descr,'button_width': butt, 'font_weight': 'bold'}, button_style='',disabled=False,)\n", "\n", - "graph_size = ToggleButtons(options=[('XXS', 50),('XS', 70),('S', 90), ('M', 100), ('L', 120)], value=100, description='Graphs size',\n", + "graph_size = IntSlider(value=100,min=10,max=200,step=5, description='Graphs size',\n", " layout=Layout(width='auto', grid_area='gr_size'), style = {'description_width': descr,'button_width': butt}, button_style='',disabled=False,)\n", "\n", "ctf_xaxis = RadioButtons(options=['α','mrad', 'A'], description='CTF x-axis:',\n", " layout=Layout(width='auto', grid_area='ctxaxis'), style = {'description_width': descr,'button_width': butt}, disabled=False,)\n", "\n", "cl_check_laa = Checkbox(value=True, description='Low anlge approximation', indent=False,\n", - " layout=Layout(width='auto', grid_area='check1'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Low angle approximation check - scattering angle < 10 deg\n", - "cl_check_pix = Checkbox(value=True, description='Recon. pixel < step size',indent=False,\n", - " layout=Layout(width='auto', grid_area='check2'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Reconstructed pixel size < scanning step size \n", + " layout=Layout(width='auto', grid_area='check1'), style = {'description_width': descr}, disabled=False,) # Low angle approximation check - scattering angle < 10 deg\n", + "cl_check_pix = Checkbox(value=True, description='Ptycho pixel < step size',indent=False,\n", + " layout=Layout(width='auto', grid_area='check2'), style = {'description_width': descr}, disabled=False,) # Reconstructed pixel size < scanning step size \n", "cl_check_def = Checkbox(value=True, description='Beam fits in probe window', indent=False,\n", - " layout=Layout(width='auto', grid_area='check3'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Is probe window big enough to accomodate defocused probe?\n", + " layout=Layout(width='auto', grid_area='check3'), style = {'description_width': descr}, disabled=False,) # Is probe window big enough to accomodate defocused probe?\n", "cl_check_1bf = Checkbox(value=True, description='Detector cover < 1α', indent=False,\n", - " layout=Layout(width='auto', grid_area='check4'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Detector cover is lower than 1 α\n", + " layout=Layout(width='auto', grid_area='check4'), style = {'description_width': descr}, disabled=False,) # Detector cover is lower than 1 α\n", "cl_check_nbf = Checkbox(value=True, description='Detector cover > 6α', indent=False,\n", - " layout=Layout(width='auto', grid_area='check5'), style = {'description_width': descr,'button_width': butt}, disabled=False,) # Detector cover is higher than 6 α\n", + " layout=Layout(width='auto', grid_area='check5'), style = {'description_width': descr}, disabled=False,) # Detector cover is higher than 6 α\n", "\n", - "footer = Button(description='Interactive results',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='whitesmoke'))\n", "\n", - "small_controls = GridBox(children=[header, method, graph_size, ctf_xaxis, cl_check_laa,cl_check_pix,cl_check_def, cl_check_1bf, cl_check_nbf, footer],\n", + "chart = ToggleButtons(options=[('STEM','micr'), ('Charts','dedic'),('Both','both')],value = 'both', description='Show',\n", + " layout=Layout(width='auto', grid_area='charts'), style = {'description_width': descr,'button_width': butt, 'font_weight': 'bold'}, button_style='',disabled=False,)\n", + "\n", + "\n", + "footer = Button(description='Interactive results',layout=Layout(width='auto', grid_area='footer'),style=ButtonStyle(button_color='whitesmoke', font_weight= 'bold'))\n", + "\n", + "small_controls = GridBox(children=[header, method, chart, graph_size, ctf_xaxis, cl_check_laa,cl_check_pix,cl_check_def, cl_check_1bf, cl_check_nbf, footer],\n", " layout=Layout(\n", " width='100%',\n", " grid_template_rows='auto auto auto',\n", - " grid_template_columns='40% 20% 20% 20%',\n", + " grid_template_columns='50% 20% 30%',\n", " grid_template_areas='''\n", - " \"header header header header \"\n", - " \"methods ctxaxis check1 check2\"\n", - " \"gr_size ctxaxis check3 check5\"\n", - " \". ctxaxis check4 . \" \n", - " \"footer footer footer footer\"\n", + " \"header header header \"\n", + " \"methods ctxaxis check1 \"\n", + " \"charts ctxaxis check2 \"\n", + " \"gr_size ctxaxis check3 \"\n", + " \". . check4 \"\n", + " \". . check5 \"\n", + " \"footer footer footer\"\n", " ''')\n", " )\n", "\n", "\n", - "def ptycho_interact(beam, aperture, aperture_res, aperture_res2, probe, cl, matrix, defocus, mag, camera, binning, dwell_time, restriction, method, ctf_xaxis, cl_check_laa, cl_check_pix, cl_check_def, cl_check_1bf, cl_check_nbf, graph_size):\n", + "def ptycho_interact(beam, aperture, aperture_res, aperture_res2, probe, cl, matrix, defocus, mag, camera, binning, dwell_time, restriction, method, ctf_xaxis, cl_check_laa, cl_check_pix, cl_check_def, cl_check_1bf, cl_check_nbf, graph_size, chart):\n", " \n", " ### SINGLE SETTING PARAMETERS ###\n", " wavelength = pty.get_wavelength(beam)*1e12\n", @@ -464,59 +468,90 @@ " pctf = pty.get_ssb_ctf()\n", " apertury = opt.apertures()\n", " \n", + "\n", + "\n", + "\n", + "\n", + "\n", " \n", " ### SAMPLE PLANE ##########################################\n", " points = 5\n", " yyy = np.append(np.linspace(0,overlap,100), 200)\n", - " wid = np.linspace(0,points-1,points) * step_size_corr\n", + " wid = np.array([-2, -1, 0, 1, 2]) * step_size_corr\n", " \n", " fig5 = make_subplots(specs=[[{\"secondary_y\": True}]])\n", - " fig5.update_xaxes(title_text=\"Sample plane (nm)\", range=[0, points*step_size_corr], showgrid=False, zeroline=False)\n", - " fig5.update_yaxes(range=[0, points*step_size_corr],showgrid=False, showticklabels=False)\n", + " fig5.update_xaxes(title_text=\"Sample plane (nm)\", range=[-4*step_size_corr, 4*step_size_corr], showgrid=False, zeroline=False)\n", + " fig5.update_yaxes(range=[-4*step_size_corr, 4*step_size_corr],showgrid=False, showticklabels=False)\n", " for x in range(len(wid)):\n", " fig5.add_trace(go.Scatter(showlegend=False, x=wid, y=np.ones(len(wid))*wid[x], marker_color='blue'),secondary_y=False)\n", " for x in range(len(wid)): \n", " for y in range(len(wid)): \n", " fig5.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", opacity=0.1, fillcolor=\"#FF7F7F\",name=False, x0=wid[x]-beam_diameter/2, y0=wid[y]-beam_diameter/2, x1=wid[x]+beam_diameter/2, y1=wid[y]+beam_diameter/2, line_color=\"red\",secondary_y=False)\n", "\n", - " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='blue', name='Step '+str(np.round(10*step_size,3))+' Å ('+str(np.round(10*step_size_corr,3))+')'))\n", + " fig5.add_annotation(x=-1*step_size_corr, y=-2*step_size_corr, ax=-2*step_size_corr, ay=-2*step_size_corr, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=-2*step_size_corr, y=-1*step_size_corr, ax=-2*step_size_corr, ay=-2*step_size_corr, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=-1.25*step_size_corr, y=-1.5*step_size_corr, text=str(np.round(10*step_size_corr,3))+' Å', showarrow=False, yshift=0)\n", + " \n", + " fig5.add_annotation(x=2*step_size_corr-0.75*beam_diameter/2, y=-2*step_size_corr-0.75*beam_diameter/2, ax=2*step_size_corr, ay=-2*step_size_corr,xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=2*step_size_corr+0.75*beam_diameter/2, y=-2*step_size_corr+0.75*beam_diameter/2, ax=2*step_size_corr, ay=-2*step_size_corr,xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig5.add_annotation(x=2.25*step_size_corr, y=-2*step_size_corr, text=str(np.round(10*beam_diameter,2))+' Å', showarrow=False, yshift=-15)\n", + " \n", + " \n", " if method == 'direct':\n", " xx = 2*wavelength/(np.sin(semi_angle_corr*omega/1000))/100\n", " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='green', name='Rec step '+str(np.round(xx[-1]/3,2))+' Å')) \n", " \n", - " fig5.add_trace(go.Scatter(showlegend=True, x=[-1], y=[-1],marker_size=12, marker_color='#FF7F7F', name='Beam ⌀ '+ str(np.round(10*beam_diameter,2))+\" Å\"))\n", " fig5.update_layout(legend=dict(orientation=\"v\",yanchor=\"bottom\",y=0.85,xanchor=\"right\",x=0.9))\n", - " fig5.add_trace(go.Scatter(showlegend=False, x = points*step_size_corr*np.ones(len(yyy)), y=yyy, opacity=1,mode='markers',marker_symbol = \"triangle-down\",marker=dict(size=20, color=yyy, colorscale='rainbow_r', showscale=False)),secondary_y=True) \n", - " fig5.update_yaxes(title_text=\"Overlap %\", range=[0, 100],showgrid=False, showticklabels=True, tickvals = [0, 20, 40,50, 60,70, 80, 90, 100], secondary_y=True)\n", - " fig5.update_layout(title={'text': \"Sample\",'y':0.95, 'x':0.30,'xanchor': 'left','yanchor': 'top'})\n", + " fig5.add_trace(go.Scatter(showlegend=False, x = -4*step_size_corr*np.ones(len(yyy)), y=yyy, opacity=1,mode='markers',marker_symbol = \"triangle-down\",marker=dict(size=20, color=yyy, colorscale='rainbow_r', showscale=False)),secondary_y=True) \n", + " fig5.update_yaxes(range=[0, 100], showgrid=False, showticklabels=False, secondary_y=True)\n", + " fig5.update_yaxes(title_text=\"Overlap %\", showgrid=False, showticklabels=True, tickvals = [-4*step_size_corr, -2*step_size_corr, 0, 2*step_size_corr, 4*step_size_corr],ticktext = ['0', '25', '50', '75', '100'], secondary_y=False)\n", + " fig5.update_layout(title={'text': \"Sample\",'y': 0.975, 'x': 0.30,'xanchor': 'left','yanchor': 'top'})\n", " fig5.update_layout(width=4.5*graph_size, height=3.0*graph_size, margin =dict(l=1.3*graph_size, r=0.3*graph_size, t=0.1*graph_size, b=0))\n", "\n", " \n", + "\n", + " \n", + " \n", + " \n", + " \n", " \n", " \n", " \n", " ### DETECTOR PLANE ##########################################\n", - " fig6 = go.Figure() \n", - " fig6.update_xaxes(title_text=\"Pixels\", range=[1, num_pixels], zeroline=False)\n", - " fig6.update_yaxes(range=[1, num_pixels], showticklabels=False,)\n", - " for x in np.linspace(1,int(np.round(covered_alfas)),int(np.round(covered_alfas))):\n", - " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", x0=num_pixels/2-x*beam_diameter_pix, y0=num_pixels/2-x*beam_diameter_pix, x1=num_pixels/2+x*beam_diameter_pix, y1=num_pixels/2+x*beam_diameter_pix,line_color=\"gray\")\n", - " fig6.add_trace(go.Scatter(showlegend=False,x=[num_pixels/2],y=[num_pixels/2-x*beam_diameter_pix],mode=\"lines+text\", text=str(int(2*x))+\"α\",textposition=\"top center\"))\n", + " fig6 = make_subplots(specs=[[{\"secondary_y\": True}]])\n", " \n", - " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", fillcolor=\"#FF7F7F\", x0=num_pixels/2-beam_diameter_pix/2, y0=num_pixels/2-beam_diameter_pix/2, x1=num_pixels/2+beam_diameter_pix/2, y1=num_pixels/2+beam_diameter_pix/2,line_color=\"#FF7F7F\")\n", - " fig6.add_annotation(x=num_pixels/2-beam_diameter_pix/2, y=num_pixels/2, ax=num_pixels/2+beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=3,arrowsize=1,arrowwidth=1,arrowcolor='black')\n", - " fig6.add_annotation(x=num_pixels/2+beam_diameter_pix/2,y=num_pixels/2, ax=num_pixels/2-beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y', text='', showarrow=True,arrowhead=3,arrowsize=1,arrowwidth=1,arrowcolor='black')\n", + " for x in np.linspace(1,int(np.round(covered_alfas)),int(np.round(covered_alfas))):\n", + " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\", x0=num_pixels/2-x*beam_diameter_pix, y0=num_pixels/2-x*beam_diameter_pix, x1=num_pixels/2+x*beam_diameter_pix, y1=num_pixels/2+x*beam_diameter_pix, line=dict(color=\"gray\",width=1,dash=\"dot\",),secondary_y=False)\n", + " fig6.add_trace(go.Scatter(showlegend=False,x=[num_pixels/2],y=[num_pixels/2+x*beam_diameter_pix],mode=\"lines+text\", text=str(int(2*x))+\"α\",textposition=\"top center\"),secondary_y=False)\n", + " \n", + " fig6.add_shape(type=\"circle\",xref=\"x\", yref=\"y\",opacity=0.5, fillcolor=\"#FF7F7F\", x0=num_pixels/2-beam_diameter_pix/2, y0=num_pixels/2-beam_diameter_pix/2, x1=num_pixels/2+beam_diameter_pix/2, y1=num_pixels/2+beam_diameter_pix/2,line_color=\"#FF7F7F\",secondary_y=False)\n", + " fig6.add_annotation(x=num_pixels/2-beam_diameter_pix/2, y=num_pixels/2, ax=num_pixels/2+beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y',text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", + " fig6.add_annotation(x=num_pixels/2+beam_diameter_pix/2,y=num_pixels/2, ax=num_pixels/2-beam_diameter_pix/2, ay=num_pixels/2, xref='x',yref='y',axref='x',ayref='y', text='', showarrow=True,arrowhead=2,arrowsize=2,arrowwidth=1,arrowcolor='black')\n", " fig6.add_annotation(x=num_pixels/2, y=num_pixels/2, text=str(np.round(beam_diameter_pix,1))+\" pix\", showarrow=False, yshift=10)\n", - " fig6.add_trace(go.Scatter(x=[1,1.5,1,1.5], y=[1,1,1.5,1.5], name='Pix ∠ ' +str(np.round(pixel_angle,2))+' mrad'))\n", - " fig6.add_trace(go.Scatter(x=[1,1.5,1,1.5], y=[1,1,1.5,1.5], name='Pixel size ' +str(np.round(size_pixel,2))+' um'))\n", - " fig6.update_layout(legend=dict(orientation=\"v\",yanchor=\"bottom\",y=0.80,xanchor=\"right\",x=1))\n", - " fig6.update_layout(title={'text': \"Detector\",'y':0.95, 'x':0.11,'xanchor': 'left','yanchor': 'top'})\n", - " fig6.update_layout(xaxis = dict(tickmode = 'array', tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels])) \n", - " fig6.update_layout(width=3.5*graph_size, height=3*graph_size, margin = dict(l=0.35*graph_size, r=0.65*graph_size, t=0.3*graph_size, b=0))\n", + " fig6.add_annotation(x=num_pixels/2, y=num_pixels/2, text=\"1 pix = \"+str(np.round(pixel_angle,2))+\" mrad\", showarrow=False, yshift=-10)\n", + " fig6.add_annotation(x=num_pixels/2, y=num_pixels/10, text=\"Real pixel size = \"+str(np.round(size_pixel,2))+\" μm\", showarrow=False, yshift=0)\n", + " fig6.add_trace(go.Scatter(showlegend=True,x=[-1], y=[-1], marker_color='red', name = \"BF disk\"),secondary_y=False)\n", + " \n", + " fig6.add_trace(go.Scatter(showlegend=False,x=[-1], y=[-1]),secondary_y=True)\n", + " \n", + " fig6.update_layout(legend=dict(orientation=\"h\",yanchor=\"bottom\",y=1.00,xanchor=\"right\",x=0.9))\n", + " fig6.update_layout(title={'text': \"Detector\",'y':0.975, 'x':0.11,'xanchor': 'left','yanchor': 'top'})\n", + " \n", + " fig6.update_xaxes(title_text=\"Pixels\", range=[1, num_pixels], zeroline=False, tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels]) \n", + " fig6.update_yaxes(range=[1, num_pixels], showticklabels=False, showgrid=False, secondary_y=False)\n", + " fig6.update_yaxes(title_text=\"Angular cover (mrad)\", range=[1, num_pixels], showticklabels=True, tickvals = [1, num_pixels/4, num_pixels/2, 3*num_pixels/4, num_pixels], ticktext = ['', np.round(detector_cover/4,1), np.round(detector_cover/2,1), np.round(3*detector_cover/4,1), np.round(detector_cover,1)], secondary_y=True) \n", + " \n", + " fig6.update_layout(width=3.5*graph_size, height=3*graph_size, margin = dict(l=0.35*graph_size, r=0.48*graph_size, t=0.3*graph_size, b=0))\n", "\n", " \n", " \n", " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " focusdepth = wavelength/(semi_angle_corr/1000)**2/1000\n", " ssb_cl, kolik = pty.cl4ssb(detector_cover_a_all)\n", " ### MICROSCOPE SCHEME ##########################################\n", @@ -654,20 +689,27 @@ " fig4.update_layout(xaxis = dict(tickmode = 'array', tickvals = opt.cameralengths()))\n", " fig4.update_layout(width=8.0*graph_size, height=2.25*graph_size, margin =dict(l=1.3*graph_size, r=0.25*graph_size, t=0.3*graph_size, b=0.3*graph_size))\n", "\n", - " ### SHOWING ###\n", - " right_column = VBox([go.FigureWidget(fig),HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", - " left_column = VBox([go.FigureWidget(fig8)])\n", - " \n", - "\n", + " \n", " \n", + " ### SHOWING ###\n", + " if chart == 'micr':\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column]) \n", " \n", - " total = HBox([left_column, right_column])\n", - " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px')\n", - " \n", + " elif chart == 'dedic':\n", + " right_column = VBox([go.FigureWidget(fig) ,HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", + " total = HBox([right_column])\n", + " \n", + " else: \n", + " right_column = VBox([go.FigureWidget(fig) ,HBox([go.FigureWidget(fig5),go.FigureWidget(fig6)]),go.FigureWidget(fig4)])\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column, right_column])\n", + " \n", + " \n", + " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px') \n", " display(total)\n", - " \n", - " \n", - " \n", + " \n", + "\n", " \n", " \n", " if method == 'iterative': \n", @@ -781,7 +823,7 @@ " \n", " dictionary[i] = ToggleButtons(options=data_tuple, button_style=button_style, icons = icons, **align6)\n", " \n", - " legend = ToggleButtons(options=['Ptycho pix. (Å)','Det. cover (mrad)','Det. cover (α)', 'Max. def. (nm)','Selected by CL ','Nominal CL (cm)',], \n", + " legend = ToggleButtons(options=['Ptycho pixel (Å)','Det. cover (mrad)','Det. cover (α)', 'Max. def. (nm)','Selected by CL ','Nominal CL (cm)',], \n", " tooltips=['Dependent on: detector, beam energy, camera length',\n", " 'Dependent on: detector, camera length',\n", " 'Dependent on: detector, camera length, probe semi-angle',\n", @@ -796,22 +838,32 @@ " \n", " \n", " \n", + " \n", " ### SHOWING ###\n", - " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", - " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", - " left_column = VBox([go.FigureWidget(fig8)]) \n", - " \n", - " \n", - " total = HBox([left_column, right_column])\n", - " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px')\n", - " \n", + " if chart == 'micr':\n", + " left_column = VBox([go.FigureWidget(fig8)])\n", + " total = HBox([left_column]) \n", + " \n", + " elif chart == 'dedic':\n", + " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", + " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", + " total = HBox([right_column])\n", + " \n", + " else: \n", + " sample_overlap = HBox([go.FigureWidget(fig5), go.FigureWidget(fig7), go.FigureWidget(fig6)]) \n", + " right_column = VBox([HBox([go.FigureWidget(fig),checks]), sample_overlap, cltab]) \n", + " left_column = VBox([go.FigureWidget(fig8)]) \n", + " total = HBox([left_column, right_column]) \n", + " \n", + " total.layout = Layout(border='solid 0px gray',margin='10px 10px 10px 40px', padding='10px 10px 10px 10px') \n", " display(total)\n", + "\n", " \n", " return \n", " \n", "gui = interactive_output(ptycho_interact, {\"beam\": beam, \"aperture\": aperture,\"aperture_res\": aperture_res,\"aperture_res2\": aperture_res2, \"probe\": probe, \"cl\": cl, \"matrix\": matrix, \"defocus\": defocus, \"mag\": mag, \"camera\": camera,\n", " \"binning\": binning, \"dwell_time\": dwell_time, \"restriction\": restriction, \"method\": method, \"ctf_xaxis\": ctf_xaxis,\n", - " \"cl_check_laa\":cl_check_laa, \"cl_check_pix\": cl_check_pix, \"cl_check_def\": cl_check_def, \"cl_check_1bf\": cl_check_1bf, \"cl_check_nbf\": cl_check_nbf, \"graph_size\": graph_size}) \n", + " \"cl_check_laa\":cl_check_laa, \"cl_check_pix\": cl_check_pix, \"cl_check_def\": cl_check_def, \"cl_check_1bf\": cl_check_1bf, \"cl_check_nbf\": cl_check_nbf, \"graph_size\": graph_size,\"chart\": chart}) \n", "\n" ] }, @@ -845,14 +897,14 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 390, "id": "7937f054-fcd0-4e67-a20f-7696f5903a94", "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "47642ef0d8b94bceb11f878293997c33", + "model_id": "6f84c5e6c29142a9bb9db08cd48aefbd", "version_major": 2, "version_minor": 0 }, @@ -865,6 +917,7 @@ } ], "source": [ + "\n", "display(HBox([controls, VBox([small_controls, gui]) ]))" ] }, @@ -888,6 +941,139 @@ "# import ipywidgets\n", "# print(\"Ipywidgets version \", ipywidgets.__version__)" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "33727cb4-c3d7-4134-be89-f10b9922fbda", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 75, + "id": "1cb358ba-42c4-420a-a0a6-0de425e670a6", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Ipywidgets version 7.6.5\n" + ] + } + ], + "source": [ + "import ipywidgets\n", + "print(\"Ipywidgets version \", ipywidgets.__version__)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a52ff861-8243-4d39-b7da-1ee8fb0ee494", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "61970cec-8678-4edf-99f6-4f832d2aad8d", + "metadata": {}, + "outputs": [], + "source": [ + "https://ipywidgets.readthedocs.io/en/latest/user_install.html" + ] + }, + { + "cell_type": "code", + "execution_count": 386, + "id": "a4628def-4c43-456e-a93b-95952e2a3dea", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
    \n", + "
  • DescriptionStyle: BoundedFloatText, BoundedIntText, Checkbox, ColorPicker, Combobox, DatePicker, Dropdown, FloatText, HTML, HTMLMath, IntText, Label, Password, Play, RadioButtons, Select, SelectMultiple, SelectionRangeSlider, SelectionSlider, Text, Textarea, ToggleButton, Valid
    • \n", + "
  • ButtonStyle: Button, FileUpload
    • \n", + "
  • SliderStyle: FloatLogSlider, FloatRangeSlider, FloatSlider, IntRangeSlider, IntSlider
    • \n", + "
  • ProgressStyle: FloatProgress, IntProgress
    • \n", + "
  • ToggleButtonsStyle: ToggleButtons
" + ], + "text/plain": [ + "" + ] + }, + "execution_count": 386, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "from collections import defaultdict\n", + "from IPython.display import HTML\n", + "import ipywidgets\n", + "from pprint import pprint\n", + "\n", + "reverse_lut = defaultdict(set)\n", + "styles = set()\n", + "for export_name in dir(ipywidgets.widgets):\n", + " export = getattr(ipywidgets.widgets, export_name)\n", + " try:\n", + " if issubclass(export, ipywidgets.Widget) and 'style' in export.class_trait_names():\n", + " reverse_lut[export.style.klass.__name__].add(export.__name__)\n", + " styles.add(export.style.klass)\n", + " except TypeError:\n", + " pass\n", + "\n", + "html = '
    '\n", + "for style, widgets in reverse_lut.items():\n", + " html = f\"{html}\\n
  • {style}: {', '.join(sorted(widgets))}
    • \"\n", + "html += \"
\"\n", + "HTML(html)" + ] + }, + { + "cell_type": "code", + "execution_count": 74, + "id": "54c1add7-f5ec-43f0-a447-e22dd03926a7", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "\n", + "
AttributeToggleButtonsStyleButtonStyleDescriptionStyleProgressStyleSliderStyle
description_width✓✓✓✓
button_width✓
font_weight✓✓
bar_color✓
handle_color✓
button_color✓
" + ], + "text/plain": [ + "" + ] + }, + "execution_count": 74, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "attributes = defaultdict(set)\n", + "base_traits = set(ipywidgets.Style.class_trait_names())\n", + "\n", + "for s in styles:\n", + " for t in s.class_trait_names():\n", + " if not t.startswith(\"_\") and t not in base_traits:\n", + " attributes[s.__name__].add(t)\n", + "all_attributes = set().union(*attributes.values())\n", + "\n", + "html = '\\n'\n", + "html = f\"{html}{ ''.join(f'' for s in attributes.keys()) }\"\n", + "for a in all_attributes:\n", + " html = f\"\"\"{html}{ ''.join(f'' for attribs in attributes.values()) }\"\"\"\n", + "html += \"
Attribute{s}
{a}{\"✓\" if a in attribs else \"\"}
\"\n", + "HTML(html)" + ] } ], "metadata": {