11import numpy as np
22from klib .baseio import *
33from scipy .ndimage import filters as ndfilter
4- from klib .glib .DrawSimulationSWCModel import simulate3DTreeModel , save_swc
4+ from klib .glib .DrawSimulationSWCModel import simulate3DTreeModel_dendrite , simulate3DTreeModel_axon , save_swc
55import copy
66import cv2 as cv
77import multiprocessing as mp
88from skimage import morphology
99import time
1010
11+
1112# init parameter
12- data_type = np .uint8
13+ data_type = np .uint16
1314internal_feature = True
1415external_feature = True
15- image_number = 4
16+ image_number = 6
1617
1718# size of training data
1819size_x = 64
1920size_y = 64
20- channel = 16
21- steps = 8
21+ channel = 32
22+
2223
2324# data dir
24- data_image_dir = 'data/sim_img/'
25- data_label_dir = 'data/sim_label/'
26- data_swc_dir = 'data/sim_swc/'
25+ data_image_dir = 'data/simulator_data/ sim_img/'
26+ data_label_dir = 'data/simulator_data/ sim_label/'
27+ data_swc_dir = 'data/simulator_data/ sim_swc/'
2728
28- train_data_dir = 'data/data_divide/'
29- train_label_dir = 'data/label_divide/'
29+ train_data_dir = 'data/simulator_data/ data_divide/'
30+ train_label_dir = 'data/simulator_data/ label_divide/'
3031
3132def main (img_num ):
3233 np .random .seed ()
@@ -40,7 +41,8 @@ def main(img_num):
4041 MAX_BOX_WIDTH = [MAX_BOX_Z ,MAX_BOX_X ,MAX_BOX_Y ]
4142
4243
43- img_sim , label_sim , swc_data = simulate3DTreeModel (MAX_BOX_WIDTH , internal_feature , external_feature , data_type = data_type )
44+ img_sim , label_sim , swc_data = simulate3DTreeModel_dendrite (MAX_BOX_WIDTH , internal_feature , external_feature , data_type = data_type )
45+ # img_sim, label_sim, swc_data = simulate3DTreeModel_axon(MAX_BOX_WIDTH, internal_feature, external_feature, data_type=data_type)
4446 SHAPE = img_sim .shape
4547
4648 img_new = np .zeros ([int (SHAPE [0 ] // 2 ), int (SHAPE [1 ]), int (SHAPE [2 ])], dtype = data_type )
@@ -103,7 +105,7 @@ def main(img_num):
103105
104106if __name__ == '__main__' :
105107 # multiprocessing
106- cpu_core_num = 4
108+ cpu_core_num = 5
107109 pool = mp .Pool (processes = cpu_core_num ) # we set cpu core is 4
108110 pool .map (main , range (0 , image_number ))
109111
0 commit comments