Mercurial > repos > imgteam > 2d_auto_threshold

comparison auto_threshold.py @ 8:699a5e9146b3 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/tree/master/tools/2d_auto_threshold/ commit c86a1b93cb7732f7331a981d13465653cc1a2790
author imgteam
date Wed, 24 Apr 2024 08:11:33 +0000
parents e5c8e7e72373
children 50fa6150e340
comparison
equal deleted inserted replaced
7:e5c8e7e72373 8:699a5e9146b3
5 See file LICENSE for detail or copy at https://opensource.org/licenses/MIT 5 See file LICENSE for detail or copy at https://opensource.org/licenses/MIT
6 """ 6 """
7 7
8 import argparse 8 import argparse
9 9
10 import giatools.io
10 import numpy as np 11 import numpy as np
11 import skimage.filters 12 import skimage.filters
12 import skimage.io
13 import skimage.util 13 import skimage.util
14 import tifffile 14 import tifffile
15
15 16
16 th_methods = { 17 th_methods = {
17 'manual': lambda thres, **kwargs: thres, 18 'manual': lambda thres, **kwargs: thres,
18 19
19 'otsu': lambda img_raw, **kwargs: skimage.filters.threshold_otsu(img_raw), 20 'otsu': lambda img_raw, **kwargs: skimage.filters.threshold_otsu(img_raw),
26 'loc_mean': lambda img_raw, bz, **kwargs: skimage.filters.threshold_local(img_raw, bz, method='mean') 27 'loc_mean': lambda img_raw, bz, **kwargs: skimage.filters.threshold_local(img_raw, bz, method='mean')
27 } 28 }
28 29
29 30
30 def do_thresholding(in_fn, out_fn, th_method, block_size, offset, threshold, invert_output=False): 31 def do_thresholding(in_fn, out_fn, th_method, block_size, offset, threshold, invert_output=False):
31 img = skimage.io.imread(in_fn) 32 img = giatools.io.imread(in_fn)
32 img = np.squeeze(img) 33 img = np.squeeze(img)
33 assert img.ndim == 2 34 assert img.ndim == 2
34 35
35 th = offset + th_methods[th_method](img_raw=img, bz=block_size, thres=threshold) 36 th = offset + th_methods[th_method](img_raw=img, bz=block_size, thres=threshold)
36 res = img > th 37 res = img > th