sravan/EECS-442-WN19-DIC-Project
1
0
Fork 0
mirror of https://github.com/balajsra/EECS-442-WN19-DIC-Project.git synced 2025-09-02 17:03:15 +00:00
Code Issues Packages Projects Releases Wiki Activity

Debugging Displacement Tracking

Browse Source 
* Implemented gradient, but not sure if it is correct.
* Debugging displacement tracking
This commit is contained in:
Sravan Balaji
2019-04-26 11:30:55 -04:00
parent e8c0d1bc71
commit 82eb35b28d
1 changed files with 25 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
src/__main__.py
View File

@@ -27,19 +27,19 @@ def find_displacement(match_method):
plt.figure(1)
reference = images[8]
compare_img = images[96]
reference = images[7]
compare_img = images[8]
plt.imshow(reference, cmap="gray", vmin=0, vmax=255)
subset_size = 5
subset_spacing = 20
subset_spacing = 30
search_size = 5
x_range = range(650, 2080, subset_spacing)
y_range = range(120, 500, subset_spacing)
im_data = ImageData(len(y_range), len(x_range))
im_data = image_data.ImageData(len(y_range), len(x_range))
disp_mag = np.zeros((len(y_range), len(x_range)))
@@ -74,9 +74,27 @@ def find_displacement(match_method):
im_data.displacement[i, j, :] = np.array([dx, dy])
disp_mag[i, j] = np.sqrt((dx ** 2) + (dy ** 2))
plt.quiver(X=x_range, Y=y_range,
U=im_data.displacement[:, :, 0], V=im_data.displacement[:, :, 1],
C=im_data.disp_mag, cmap=plt.cm.jet)
# Strain in x-direction
im_data.strain[:, :, 0] = np.gradient(
im_data.displacement[:, :, 0], # x displacements
im_data.location[:, 0, 1], # y (row) positions
im_data.location[0, :, 0])[1] # x (col) positions
im_data.strain[:, :, 1] = np.gradient(
im_data.displacement[:, :, 1], # y displacements
im_data.location[:, 0, 1], # y (row) positions
im_data.location[0, :, 0])[0] # x (col) positions
# plt.quiver(X=im_data.location[0, :, 0], Y=im_data.location[:, 0, 1],
# U=im_data.strain[:, :, 0], V=im_data.strain[:, :, 1],
# C=strain_mag, cmap=plt.cm.jet)
plt.quiver(im_data.location[:, :, 0], # x coordinates of arrow locations
im_data.location[:, :, 1], # y coordinates of arrow locations
im_data.displacement[:, :, 0], # x components of arrow vectors
im_data.displacement[:, :, 1], # y components of arrow vectors
disp_mag, # arrow color (vector magnitude)
cmap=plt.cm.jet) # color map (jet)
plt.show()