Mercurial > repos > imgteam > 2d_local_threshold

comparison auto_local_threshold.py @ 0:a20c14eb5f98 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/BMCV/galaxy-image-analysis/tools/2d_auto_local_threshold/ commit f3903c908786b9b0ea5a46e9cc35ee025770ecda
author imgteam
date Thu, 18 Jul 2019 09:22:07 -0400
parents
children c90b91f4a07b
comparison
equal deleted inserted replaced
-1:000000000000 0:a20c14eb5f98
1 import argparse
2 import sys
3 import skimage.io
4 import skimage.filters
5 import skimage.util
6
7 threshOptions = {
8 'gaussian': lambda img_raw, bz: skimage.filters.threshold_local(img_raw, bz, method='gaussian'),
9 'mean': lambda img_raw, bz: skimage.filters.threshold_local(img_raw, bz, method='mean'),
10 'median': lambda img_raw, bz: skimage.filters.threshold_local(img_raw, bz, method='median')
11 }
12
13 if __name__ == "__main__":
14 parser = argparse.ArgumentParser(description='Segment Foci')
15 parser.add_argument('input_file', type=argparse.FileType('r'), default=sys.stdin, help='input file')
16 parser.add_argument('out_file', type=argparse.FileType('w'), default=sys.stdin, help='out file (TIFF)')
17 parser.add_argument('block_size', type=int, default=5, help='Odd size of pixel neighborhood which is used to calculate the threshold value')
18 parser.add_argument('thresh_type', choices=threshOptions.keys(), help='thresholding method')
19 parser.add_argument('dark_background', default=True, type=bool, help='True if background is dark')
20 args = parser.parse_args()
21
22 img_in = skimage.io.imread(args.input_file.name)
23 thresh = threshOptions[args.thresh_type](img_in, args.block_size)
24
25 if args.dark_background:
26 res = img_in > thresh
27 else:
28 res = img_in <= thresh
29
30 res = skimage.util.img_as_uint(res)
31 skimage.io.imsave(args.out_file.name, res, plugin="tifffile")