| 387 |
387 |
return R_2, pc
|
| 388 |
388 |
|
| 389 |
389 |
|
| 390 |
|
def analyze_curve1(name, y, pixelsize, outfilename, graph_name='graph.png'):
|
|
390 |
def circle(R, c):
|
|
391 |
phi=np.linspace(np.pi/6., 5*np.pi/6., 51)
|
|
392 |
x=c[0]+R*np.cos(phi)
|
|
393 |
y=c[1]+R*np.sin(phi)
|
|
394 |
return x,y
|
|
395 |
|
|
396 |
|
|
397 |
|
|
398 |
def analyze_curve1(name, y, pixelsize, outfilename, title, graph_name='graph.png'):
|
| 391 |
399 |
# isotropic pixels
|
| 392 |
400 |
y[:,0] = y[:,0]*pixelsize[0]
|
| 393 |
401 |
y[:,1] = y[:,1]*pixelsize[1]
|
| ... | ... | |
| 418 |
426 |
middle_point = yy[int(len(yy)/2)]
|
| 419 |
427 |
# compute curvature radius by fitting the sepal section to a circle
|
| 420 |
428 |
yyy=y[first:last]
|
| 421 |
|
R, pc = compute_curvature_radius(yyy[:,0],yyy[:,1])
|
|
429 |
R, pcR = compute_curvature_radius(yyy[:,0],yyy[:,1])
|
| 422 |
430 |
print("R=",R," um")
|
| 423 |
431 |
# cutting curve in pieces in order to estimate curvature radius on portions of the curve
|
| 424 |
432 |
slength=len(yyy)
|
| ... | ... | |
| 428 |
436 |
y423=yyy[int(slength/4):3*int(slength/4)]
|
| 429 |
437 |
y44=yyy[3*int(slength/4):slength]
|
| 430 |
438 |
# and computing the curvature radius for each portion
|
| 431 |
|
R21, pc = compute_curvature_radius(y21[:,0],y21[:,1])
|
| 432 |
|
R22, pc = compute_curvature_radius(y22[:,0],y22[:,1])
|
| 433 |
|
R41, pc = compute_curvature_radius(y41[:,0],y41[:,1])
|
| 434 |
|
R423, pc = compute_curvature_radius(y423[:,0],y423[:,1])
|
| 435 |
|
R44, pc = compute_curvature_radius(y44[:,0],y44[:,1])
|
| 436 |
|
fig = plt.figure()
|
| 437 |
|
plt.plot(yy[:,0], yy[:,1], 'bo')
|
| 438 |
|
plt.plot(y[:,0], y[:,1], 'r')
|
| 439 |
|
plt.plot(y[first][0], y[first][1], 'ro')
|
| 440 |
|
plt.plot(middle_point[0], middle_point[1], 'ro')
|
| 441 |
|
plt.plot(y[last][0], y[last][1], 'ro')
|
| 442 |
|
plt.plot(yy[height_index][0], yy[height_index][1], 'go')
|
| 443 |
|
plt.plot(pc[0], pc[1], 'bx')
|
|
439 |
R21, pc21 = compute_curvature_radius(y21[:,0],y21[:,1])
|
|
440 |
R22, pc22 = compute_curvature_radius(y22[:,0],y22[:,1])
|
|
441 |
R41, pc41 = compute_curvature_radius(y41[:,0],y41[:,1])
|
|
442 |
R423, pc423 = compute_curvature_radius(y423[:,0],y423[:,1])
|
|
443 |
R44, pc44 = compute_curvature_radius(y44[:,0],y44[:,1])
|
|
444 |
fig = plt.figure(figsize=(15, 5))
|
|
445 |
plt.plot(yy[:,0], yy[:,1], color='lightgreen', linewidth=10, label='length='+"{:.0f}".format(length)+" um")
|
|
446 |
#plt.plot(y[:,0], y[:,1], '--')
|
|
447 |
#plt.plot(y[first][0], y[first][1], 'ro')
|
|
448 |
#plt.plot(middle_point[0], middle_point[1], 'ro')
|
|
449 |
#plt.plot(y[last][0], y[last][1], 'ro')
|
|
450 |
#plt.plot(yy[height_index][0], yy[height_index][1], 'go')
|
|
451 |
#plt.plot(pcR[0], pcR[1], 'ro', color='royalblue')
|
|
452 |
#plt.plot(pc423[0], pc423[1], 'ro', color='darkorange')
|
|
453 |
xR,yR = circle(R, pcR)
|
|
454 |
plt.plot(xR, yR, "--",color='royalblue', linewidth=3, label='R='+"{:.0f}".format(R)+" um")
|
|
455 |
xR,yR = circle(R423, pc423)
|
|
456 |
plt.plot(xR, yR, "--", color='tomato', linewidth=3, label='R423='+"{:.0f}".format(R423)+" um")
|
| 444 |
457 |
plt.gca().set_aspect('equal', adjustable='box')
|
|
458 |
plt.xlabel('x (um)')
|
|
459 |
plt.ylabel('z (um)')
|
|
460 |
plt.legend()
|
|
461 |
plt.grid(True)
|
|
462 |
plt.title(title)
|
| 445 |
463 |
#plt.show()
|
| 446 |
464 |
fig.savefig(graph_name)
|
| 447 |
465 |
#plt.close()
|