ifp.startlogger(filename='/media/julian/Daten/neuraldata/cremi_2016/develop/160927_determine_shortest_paths/160929_pthsovrdist_pow10/160929.log', type='a') # # Cropping # if crop: # ifp.crop([10, 200, 200], [110, 712, 712]) # ifp.write() # ifp.logging('ifp.get_image = {}', ifp.get_image('labels')[0, 0, 0]) # ifp.logging('ifp.amax = {}\n', ifp.amax('labels')) ifp.logging('keys() = {}', ifp.get_data().keys()) # For multiple labels ifp.logging('Starting label iterator for {} labels', len(ifp.anytask_rtrn(np.unique, ids='labels'))) c = 0 for lblo in ifp.label_bounds_iterator('labels', 'curlabel', ids='disttransf', targetids='curdisttransf', maskvalue=0, value=0): ifp.logging('------------\nCurrent label {} in iteration {}', lblo['label'], c) ifp.logging('Bounding box = {}', lblo['bounds']) local_maxima_found = find_local_maxima(ifp) ifp.logging('Local maxima found: {}', local_maxima_found) if local_maxima_found: if ifp.amax('locmax') != 0:
# Seed the randomize function random.seed(1) # # Choose one label # lbl = random.choice(labels) # ifp.logging('Choice: {}', lbl) # # ifp.astype(np.uint32, ids='labels') # (grag, rag) = ifp.anytask_rtrn(graphs.gridRegionAdjacencyGraph, ignoreLabel=0, ids='labels') # ifp.logging('RAG: {}', rag) # # ifp.logging('Node ids: {}', rag.nodeIds()) ifp.astype(np.uint32, ids='labels') (grag, rag) = ifp.anytask_rtrn(graphs.gridRegionAdjacencyGraph, ignoreLabel=0, ids='labels') edge_ids = rag.edgeIds() ifp.logging('Edge ids: {}', edge_ids) # Type 1: # Select edges by size (smallest edges) ifp.logging('Number of edgeLengths = {}', len(rag.edgeLengths())) edgelen_ids = dict(zip(edge_ids, rag.edgeLengths())) ifp.logging('edgelen_ids = {}', edgelen_ids) sorted_edgelens = np.sort(rag.edgeLengths()) # smallest_merge_lens = sorted_edgelens[0:numberbysize] ifp.logging('Lengths selected for merging: {}', smallest_merge_lens) # smallest_merge_ids = [] for x in smallest_merge_lens:
type='a') # # Cropping # if crop: # ifp.crop([10, 200, 200], [110, 712, 712]) # ifp.write() # ifp.logging('ifp.get_image = {}', ifp.get_image('labels')[0, 0, 0]) # ifp.logging('ifp.amax = {}\n', ifp.amax('labels')) ifp.logging('keys() = {}', ifp.get_data().keys()) # For multiple labels ifp.logging('Starting label iterator for {} labels', len(ifp.anytask_rtrn(np.unique, ids='labels'))) c = 0 for lblo in ifp.label_bounds_iterator('labels', 'curlabel', ids='disttransf', targetids='curdisttransf', maskvalue=0, value=0): ifp.logging('------------\nCurrent label {} in iteration {}', lblo['label'], c) ifp.logging('Bounding box = {}', lblo['bounds']) local_maxima_found = find_local_maxima(ifp)