ROI displacement to ROI e11 strain

[Patcurran]

hi everyone,

working on master branch v 0.93.6

i wanted to calculate the average stain in the loading direction for each region of interest (ROI) for each step. Working with Dongchen we got the code below. please could someone double check it make sense. :)

Displacement = DicMap.xd.reshape(np.shape(DicMap))
Position = DicMap.xc.reshape(np.shape(DicMap))

strain_list=[]
crop = 80
row = 20
displacement = np.mean(Displacement[crop:-crop,crop:row+crop])- np.mean(Displacement[crop:-crop,-crop-row:-crop])
position = np.mean(Position[crop:-crop,crop:row+crop])- np.mean(Position[crop:-crop,-crop-row:-crop])
strain_list.append(abs(displacement)/abs(position))
strain_list

plt.imshow(Displacement[:,:],cmap='viridis')
#plt.imshow(Displacement[crop:-crop,crop:-crop])
# Define rectangle patches
rect1 = Rectangle((crop, crop), row, Displacement.shape[0] - 2*crop, linewidth=1, edgecolor='r', facecolor='none')  # Rectangle 1
rect2 = Rectangle((Displacement.shape[1] - crop - row, crop), row, Displacement.shape[0] - 2*crop, linewidth=1, edgecolor='r', facecolor='none')  # Rectangle 2

# Add rectangles to the plot
plt.gca().add_patch(rect1)
plt.gca().add_patch(rect2)
plt.colorbar() 
plt.show()

[Patcurran]

strain_list=[]
crop = 80
row = 20
for n in range(1,20):
    if n <10:
        DicMap = hrdic.Map(DicFilePath, "B0000{}.txt".format(n))
        Displacement = DicMap.xd.reshape(np.shape(DicMap))
        Position = DicMap.xc.reshape(np.shape(DicMap))
        print(n)
        displacement = np.mean(Displacement[crop:-crop,crop:row+crop])- np.mean(Displacement[crop:-crop,-crop-row:-crop])
        position = np.mean(Position[crop:-crop,crop:row+crop])- np.mean(Position[crop:-crop,-crop-row:-crop])
        strain_list.append(abs(displacement)/abs(position))    
    else:
        DicMap = hrdic.Map(DicFilePath, "B000{}.txt".format(n))
        Displacement = DicMap.xd.reshape(np.shape(DicMap))
        Position = DicMap.xc.reshape(np.shape(DicMap))
        print(n)
        displacement = np.mean(Displacement[crop:-crop,crop:row+crop])- np.mean(Displacement[crop:-crop,-crop-row:-crop])
        position = np.mean(Position[crop:-crop,crop:row+crop])- np.mean(Position[crop:-crop,-crop-row:-crop])
        strain_list.append(abs(displacement)/abs(position))    
strain_list