Compare commits
4 Commits
9eb24b9be0
...
61d30e9365
| Author | SHA1 | Date | |
|---|---|---|---|
| 61d30e9365 | |||
| 4612dd5a6e | |||
| 5510e4086b | |||
| 0c6d3629e1 |
41
kisli.py
41
kisli.py
@@ -138,13 +138,15 @@ def done():
|
|||||||
|
|
||||||
def calc_vector(px, py, point_specific_matrix):
|
def calc_vector(px, py, point_specific_matrix):
|
||||||
vect = [0, 0]
|
vect = [0, 0]
|
||||||
|
print("spcific point = {}, {}".format(px, py))
|
||||||
for x in range(9):
|
for x in range(9):
|
||||||
for y in range(10):
|
for y in range(10):
|
||||||
if(px != x or py != y):
|
if(px != x or py != y):
|
||||||
diff = [px-x, py-y] #Abstandsvektor
|
diff = [x-px, y-py] #Abstandsvektor
|
||||||
mag = np.sqrt(diff[0]**2 + diff[1]**2) #betrag
|
mag = np.sqrt(diff[0]**2 + diff[1]**2) #betrag
|
||||||
vect[0] += point_specific_matrix[x][y]*diff[0] / mag
|
vect[0] += point_specific_matrix[x][y]*diff[0] / mag
|
||||||
vect[1] += point_specific_matrix[x][y]*diff[1] / mag
|
vect[1] += point_specific_matrix[x][y]*diff[1] / mag
|
||||||
|
print("{}: vect += {} / {} == {}".format([x,y],diff, mag, vect))
|
||||||
return vect
|
return vect
|
||||||
|
|
||||||
""" ==============================================================================================================
|
""" ==============================================================================================================
|
||||||
@@ -181,14 +183,40 @@ if args.fromfile:
|
|||||||
vectormap = np.zeros(shape=(90, 2))
|
vectormap = np.zeros(shape=(90, 2))
|
||||||
i = 0
|
i = 0
|
||||||
for row in full_matrix:
|
for row in full_matrix:
|
||||||
|
print(i)
|
||||||
point_specific_matrix=get_mapped(row)
|
point_specific_matrix=get_mapped(row)
|
||||||
vectormap[i] = calc_vector(0, 0, point_specific_matrix)
|
vectormap[i] = calc_vector(channels.map[i][0], channels.map[i][1], point_specific_matrix)
|
||||||
|
i+=1
|
||||||
|
print()
|
||||||
|
|
||||||
|
|
||||||
|
if args.view == "N":
|
||||||
|
neighview = np.zeros(shape=(30,32))
|
||||||
|
|
||||||
|
i = 0
|
||||||
|
for row in full_matrix:
|
||||||
|
point_specific_matrix=get_mapped(row)
|
||||||
|
|
||||||
|
#coordinates of specific point
|
||||||
|
y_pos = channels.map[i][1] - 1
|
||||||
|
x_pos = channels.map[i][0] - 1
|
||||||
|
|
||||||
|
if(x_pos+1 < 9):
|
||||||
|
neighview[(x_pos*3)+2][(y_pos*3)+1] = point_specific_matrix[x_pos+1][y_pos]
|
||||||
|
|
||||||
|
if(x_pos-1 >= 0):
|
||||||
|
neighview[(x_pos*3)+0][(y_pos*3)+1] = point_specific_matrix[x_pos-1][y_pos]
|
||||||
|
|
||||||
|
if(y_pos+1 < 10):
|
||||||
|
neighview[(x_pos*3)+1][(y_pos*3)+2] = point_specific_matrix[x_pos][y_pos+1]
|
||||||
|
|
||||||
|
if(y_pos-1 >= 0):
|
||||||
|
neighview[(x_pos*3)+1][(y_pos*3)+0] = point_specific_matrix[x_pos][y_pos-1]
|
||||||
|
|
||||||
i+=1
|
i+=1
|
||||||
|
|
||||||
mapped_vectormap = get_vector_mapped(vectormap)
|
viewer = liveview.ResistanceView(maximum=4000)
|
||||||
mapped_vectormap /= max(mapped_vectormap.min(), mapped_vectormap.max(), key=abs)
|
viewer.updateView(neighview)
|
||||||
print(mapped_vectormap)
|
|
||||||
viewer = liveview.VectorView(mapped_vectormap)
|
|
||||||
|
|
||||||
else:
|
else:
|
||||||
s = socket.socket() # Establish a TCP/IP socket object
|
s = socket.socket() # Establish a TCP/IP socket object
|
||||||
@@ -221,6 +249,7 @@ if args.view == "M":
|
|||||||
viewer = liveview.ResistanceView(maximum=full_matrix.max())
|
viewer = liveview.ResistanceView(maximum=full_matrix.max())
|
||||||
viewer.updateView(get_mapped(np.mean(full_matrix, axis = 0)))
|
viewer.updateView(get_mapped(np.mean(full_matrix, axis = 0)))
|
||||||
|
|
||||||
|
|
||||||
if args.tofile:
|
if args.tofile:
|
||||||
np.savetxt(args.tofile, full_matrix, delimiter="\t")
|
np.savetxt(args.tofile, full_matrix, delimiter="\t")
|
||||||
|
|
||||||
|
|||||||
@@ -12,7 +12,7 @@ class ResistanceView:
|
|||||||
self.fig.canvas.flush_events()
|
self.fig.canvas.flush_events()
|
||||||
|
|
||||||
class VectorView:
|
class VectorView:
|
||||||
def __init__(self, vectormap):
|
def __init__(self, vectormap, title=""):
|
||||||
# Vector origin location
|
# Vector origin location
|
||||||
X = np.repeat(range(9), 10)
|
X = np.repeat(range(9), 10)
|
||||||
Y = np.resize(range(10), 10*9)
|
Y = np.resize(range(10), 10*9)
|
||||||
@@ -23,11 +23,11 @@ class VectorView:
|
|||||||
|
|
||||||
# Creating plot
|
# Creating plot
|
||||||
pyplot.quiver(X, Y, U, V, color='b', units='xy', scale=1)
|
pyplot.quiver(X, Y, U, V, color='b', units='xy', scale=1)
|
||||||
pyplot.title('Single Vector')
|
pyplot.title(title)
|
||||||
|
|
||||||
## x-lim and y-lim
|
## x-lim and y-lim
|
||||||
#pyplot.xlim(-2, 5)
|
pyplot.xlim(-1, 9)
|
||||||
#pyplot.ylim(-2, 2.5)
|
pyplot.ylim(-2, 11)
|
||||||
|
|
||||||
# Show plot with grid
|
# Show plot with grid
|
||||||
pyplot.grid()
|
pyplot.grid()
|
||||||
|
|||||||
Reference in New Issue
Block a user