Commit f2d2242a authored by Olga Stamati's avatar Olga Stamati
Browse files

sort fig. 3-4-5

parent 3dcc707b
Pipeline #78319 passed with stage
in 1 minute and 25 seconds
......@@ -3,12 +3,11 @@ import matplotlib.pyplot as plt
import numpy
import tifffile
import scipy.ndimage
from skimage.transform import hough_line, hough_line_peaks
from matplotlib import cm
from matplotlib.colors import LogNorm
import progressbar
plt.style.use('./tools/radioSphere.mplstyle')
plt.style.use('../../tools/radioSphere.mplstyle')
verbose = False
use_cache = False
......@@ -36,7 +35,8 @@ if use_cache:
else:
to_plot = []
# for radialOffsetMM in numpy.linspace(0,9.5,50):
for projectedAspectRatio in progressbar.progressbar(numpy.linspace(1,3,20)):
#for projectedAspectRatio in progressbar.progressbar(numpy.linspace(1,3,20)):
for projectedAspectRatio in (numpy.linspace(1,3,20):
alpha = numpy.arccos(1./projectedAspectRatio)
radialOffsetMM = numpy.tan(alpha)*sourceObjectDistMM
# alpha = numpy.arctan(radialOffsetMM/sourceObjectDistMM) # angle of particle centroid from beam centre
......@@ -131,4 +131,4 @@ plt.xlabel('Projected particle aspect ratio (-)')
plt.ylabel('Error in position detection\n(radii)')
plt.subplots_adjust(left=0.19,bottom=0.21,top=0.99,right=0.99)
plt.savefig('./figures/aspectRatioSensitivity.pdf')
plt.savefig('aspectRatioSensitivity.pdf')
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment