This repository has been archived by the owner on Jul 27, 2019. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 1
/
simple_iterator.py
55 lines (46 loc) · 1.71 KB
/
simple_iterator.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
import numpy as np
def template1(img)
# returns the list of coordinates where template1 is discovered
return list_coordinates
#iterate in three directions and mark 1 if template 1 is found
def new_god(list_axial,list_saggital,list_coronal):
length=len(list_axial)
m=np.zeros((length,length,length),dtype=np.int8)
fidu_coordinates = []
for x in range(length):
for y in range(length):
for z in range(length):
if (y,x) in template1(list_axial[z]) or (length-z,y) in template1(list_saggital[x]) or (length-z,x) in template1(list_coronal[y]):
m[x,y,z]=1
#refining over axial images
for i in range(length):
for j in range(length):
k = 0
while k < length:
if m[i,j,k] == 1:
fidu_coordinates.append(i,j,k)
k += 1
while(m[i,j,k] == 1):
m[i,j,k] = 0
k += 1
#refining over saggital images
for i in range(length):
for j in range(length):
k=0
if m[k,i,j]==1:
fidu_coordinates.append(k,i,j)
k+=1
while(m[k,i,j] == 1):
m[k,i,j] = 0
k += 1
#refining over coronal images
for i in range(length):
for j in range(length):
k=0
if m[j,k,i]==1:
fidu_coordinates.append(j,k,i)
k+=1
while(m[j,k,i] == 1):
m[j,k,i] = 0
k += 1
return fidu_coordinates