ich habe mit einem Numpy Array ein Voxel Gitter erstellt, in meinem array bedeutet Array[100][100][40] = 20: Der Voxel an der Stelle x=100 y=100 z=40 besitzt den Wert 20.
So nun möchte ich Messdaten in diesem array abspeichern, ich bekomme Messdate in dieser Form: x=10 y=10 z=10 r=3 s= 77, d.h. alle Voxel auf mit einem Radius 3 auf der Halbkugel um den Punkt (x,y,z) sollen den Wert 77 addiert bekommen. Jedoch kenn ich keine wirklich geeignete Funktion, mein Lösungsansatz wäre zurzeit:
Code: Alles auswählen
def Koordinaten(x, y, z, r, s):
a = (100, 100, 100)
koord = np.zeros(a)
phi = 0
eta = 0
while phi <= 2pi:
while eta <= pi:
x = r * sin eta * cos phi
y = r * sin eta * sin phi
z = t * cos eta
a[x][y][z] = s
eta + = 0.1
phi += 0.1