def main(par): verb = par.bool('verb', False) inp = par.string('inp') # ifile.rsf out_basenm = os.path.basename(inp).rstrip(ivlad.ext) + '_split' outdir = par.string('outdir', (out_basenm + ivlad.ext)) # Default is ifile_split.rsf nthick = par.int('nthick', 1) # slice thickness ivlad.chk_param_limit(nthick, 'nthick', 1, '>') lastdim = str(ivlad.ndims(inp)) n = int(ivlad.getout('sfget', ['parform=n', 'n' + lastdim], inp)) ivlad.chk_param_limit(nthick, 'nthick', n, '<') if not os.path.isdir(outdir): os.mkdir(outdir) # Slices will be called 01.rsf, etc. Avoid datapath aliasing: dpath = 'datapath=$DATAPATH/' + out_basenm + '_' nslices_whole = int(math.floor(n / nthick)) f = 0 for i in range(nslices_whole): i_str = ivlad.add_zeros(i, nslices_whole) i_slc = os.path.join(outdir, i_str + '_stdout' + ivlad.ext) cmd = '<%s sfwindow f%s=%d ' % (inp, lastdim, f) cmd += 'n%s=%d %s > %s' % (lastdim, nthick, dpath, i_slc) ivlad.exe(cmd, verb) f += nthick ist = n - nslices_whole * nthick # Incomplete Slice Thickness if ist > 0: i_str = ivlad.add_zeros(i + 1, nslices_whole) i_slc = os.path.join(outdir, i_str + '_stdout' + ivlad.ext) cmd = '<%s sfwindow f%s=%d ' % (inp, lastdim, f) cmd += 'n%s=%d > %s' % (lastdim, ist, i_slc) ivlad.exe(cmd, verb) return ivlad.unix_success
def main(par): verb = par.bool('verb', False) inp = par.string('inp') # ifile.rsf out_basenm = os.path.basename(inp).rstrip(ivlad.ext) + '_split' outdir = par.string('outdir',(out_basenm+ivlad.ext)) # Default is ifile_split.rsf nthick = par.int('nthick', 1) # slice thickness ivlad.chk_param_limit(nthick, 'nthick', 1, '>') lastdim = str(ivlad.ndims(inp)) n = int(ivlad.getout('sfget',['parform=n','n'+lastdim], inp)) ivlad.chk_param_limit(nthick, 'nthick', n, '<') if not os.path.isdir(outdir): os.mkdir(outdir) # Slices will be called 01.rsf, etc. Avoid datapath aliasing: dpath = 'datapath=$DATAPATH/' + out_basenm + '_' nslices_whole = int(math.floor(n/nthick)) f = 0 for i in range(nslices_whole): i_str = ivlad.add_zeros(i,nslices_whole) i_slc = os.path.join(outdir, i_str + '_stdout' + ivlad.ext) cmd = '<%s sfwindow f%s=%d ' % (inp, lastdim,f) cmd += 'n%s=%d %s > %s' % (lastdim, nthick, dpath, i_slc) ivlad.exe(cmd, verb) f += nthick ist = n - nslices_whole * nthick # Incomplete Slice Thickness if ist > 0: i_str = ivlad.add_zeros(i+1,nslices_whole) i_slc = os.path.join(outdir, i_str + '_stdout' + ivlad.ext) cmd = '<%s sfwindow f%s=%d ' % (inp, lastdim, f) cmd += 'n%s=%d > %s' % (lastdim, ist, i_slc) ivlad.exe(cmd, verb) return ivlad.unix_success
def main(par): inp = par.string('inp') # input file out = par.string('out') # output file verb = par.bool('verb', False) # if y, print system commands, outputs ivlad.chk_file_dims(inp, 2) n1 = ivlad.getout('sfget',['parform=n','n1'], inp, verb) n2 = ivlad.getout('sfget',['parform=n','n2'], inp, verb) n1 = int(n1) n2 = int(n2) h = par.int('h', 768) # output height w = par.int('w', 1024) # output width has_h = par.string('h') has_w = par.string('w') if has_h and has_w: # both w and h were read. Check for sanity: ivlad.chk_param_limit(w, 'w') ivlad.chk_param_limit(h, 'h') if (h,w) == (n1,n2): ivlad.msg('Change h or w if you want out != inp') sf.cp(inp, out, verb) return ivlad.unix_success # Transform h and w to pixels, if they are not # No default value for par.string -- Quirk of rsf, replicated in rsfbak unit = par.string('unit', 'px') # unit of h and w. Can be: px, mm, cm, in ivlad.chk_par_in_list(unit,['mm','cm','in','px']) if unit != 'px': ppi = par.int('ppi') # outp. resolution (px/in). Necessary when unit!=px ivlad.chk_param_limit(ppi, 'ppi') # Transform w and h to px if unit == 'in': scale = 1 elif unit == 'mm': scale = 254 elif unit == 'cm': scale = 25.4 w *= ppi / float(scale) h *= ppi / float(scale) # Don't worry, we'll convert to int after prar block del scale # Now h and w are in pixels # If prar=y, then h and/or w define a bounding box. # Find the dimensions of the image inside this box prar = par.bool('prar', True) # if y, PReserve Aspect Ratio of input if prar: # preserve aspect ratio if has_h and not has_w: w = n2 * float(h) / n1 elif has_w and not has_h: h = n1 * float(w) / n2 else: # Full bounding box specified hscale = float(h) / n1 wscale = float(w) / n2 if hscale < wscale: w = n2 * hscale else: h = n1 * wscale h = int(h) w = int(w) assert h > 1 assert w > 1 h = ivlad.valswitch(h, n1, None) w = ivlad.valswitch(w, n2, None) tmp = tempfile.mktemp(dir=rsf.path.datapath()) # Interpolation and, if needed, bandpass if h != None: d1 = ivlad.getout('sfget', ['parform=n','d1'], inp, verb) d1 = float(d1) * (n1-1)/float(h-1) if h < n1: ready_for_remap_1 = tmp + '1' sf.bandpass(inp, ready_for_remap_1, fhi=0.5/d1, verb=verb) rem2del_junk1 = True else: ready_for_remap_1 = inp rem2del_junk1 = False if w: out_remap1 = tmp + '2' rem2del_junk2 = True else: out_remap1 = out rem2del_junk2 = False sf.remap1(ready_for_remap_1, out_remap1, n1=h, d1=d1, verb=verb) if rem2del_junk1: sf.rm(ready_for_remap_1, verb) else: # no action on axis 1 out_remap1 = inp rem2del_junk2 = False if w != None: d2 = ivlad.getout('sfget', ['parform=n','d2'], inp, verb) d2 = float(d2) * (n2-1)/float(w-1) out_transp1 = tmp + '3' sf.transp(out_remap1, out_transp1, verb=verb) if rem2del_junk2: sf.rm(out_remap1, verb) if w < n2: ready_for_remap_2 = tmp + '4' sf.bandpass(out_transp1, ready_for_remap_2, fhi=0.5/d2, verb=verb) rem2del_junk4 = True else: ready_for_remap_2 = out_transp1 rem2del_junk4 = False ready_for_transp2 = tmp + '5' sf.remap1(ready_for_remap_2,ready_for_transp2,n1=w,d1=d2,verb=verb) sf.rm(out_transp1, verb) if rem2del_junk4: sf.rm(ready_for_remap_2, verb) sf.transp(ready_for_transp2, out, verb=verb) sf.rm(ready_for_transp2, verb) return ivlad.unix_success
def main(par): inp = par.string('inp') # input file out = par.string('out') # output file verb = par.bool('verb', False) # if y, print system commands, outputs ivlad.chk_file_dims(inp, 2) n1 = ivlad.getout('sfget', ['parform=n', 'n1'], inp, verb) n2 = ivlad.getout('sfget', ['parform=n', 'n2'], inp, verb) n1 = int(n1) n2 = int(n2) h = par.int('h', 768) # output height w = par.int('w', 1024) # output width has_h = par.string('h') has_w = par.string('w') if has_h and has_w: # both w and h were read. Check for sanity: ivlad.chk_param_limit(w, 'w') ivlad.chk_param_limit(h, 'h') if (h, w) == (n1, n2): ivlad.msg('Change h or w if you want out != inp') sf.cp(inp, out, verb) return ivlad.unix_success # Transform h and w to pixels, if they are not # No default value for par.string -- Quirk of rsf, replicated in rsfbak unit = par.string('unit', 'px') # unit of h and w. Can be: px, mm, cm, in ivlad.chk_par_in_list(unit, ['mm', 'cm', 'in', 'px']) if unit != 'px': ppi = par.int( 'ppi') # outp. resolution (px/in). Necessary when unit!=px ivlad.chk_param_limit(ppi, 'ppi') # Transform w and h to px if unit == 'in': scale = 1 elif unit == 'mm': scale = 254 elif unit == 'cm': scale = 25.4 w *= ppi / float(scale) h *= ppi / float(scale) # Don't worry, we'll convert to int after prar block del scale # Now h and w are in pixels # If prar=y, then h and/or w define a bounding box. # Find the dimensions of the image inside this box prar = par.bool('prar', True) # if y, PReserve Aspect Ratio of input if prar: # preserve aspect ratio if has_h and not has_w: w = n2 * float(h) / n1 elif has_w and not has_h: h = n1 * float(w) / n2 else: # Full bounding box specified hscale = float(h) / n1 wscale = float(w) / n2 if hscale < wscale: w = n2 * hscale else: h = n1 * wscale h = int(h) w = int(w) assert h > 1 assert w > 1 h = ivlad.valswitch(h, n1, None) w = ivlad.valswitch(w, n2, None) tmp = tempfile.mktemp(dir=rsf.path.datapath()) # Interpolation and, if needed, bandpass if h != None: d1 = ivlad.getout('sfget', ['parform=n', 'd1'], inp, verb) d1 = float(d1) * (n1 - 1) / float(h - 1) if h < n1: ready_for_remap_1 = tmp + '1' sf.bandpass(inp, ready_for_remap_1, fhi=0.5 / d1, verb=verb) rem2del_junk1 = True else: ready_for_remap_1 = inp rem2del_junk1 = False if w: out_remap1 = tmp + '2' rem2del_junk2 = True else: out_remap1 = out rem2del_junk2 = False sf.remap1(ready_for_remap_1, out_remap1, n1=h, d1=d1, verb=verb) if rem2del_junk1: sf.rm(ready_for_remap_1, verb) else: # no action on axis 1 out_remap1 = inp rem2del_junk2 = False if w != None: d2 = ivlad.getout('sfget', ['parform=n', 'd2'], inp, verb) d2 = float(d2) * (n2 - 1) / float(w - 1) out_transp1 = tmp + '3' sf.transp(out_remap1, out_transp1, verb=verb) if rem2del_junk2: sf.rm(out_remap1, verb) if w < n2: ready_for_remap_2 = tmp + '4' sf.bandpass(out_transp1, ready_for_remap_2, fhi=0.5 / d2, verb=verb) rem2del_junk4 = True else: ready_for_remap_2 = out_transp1 rem2del_junk4 = False ready_for_transp2 = tmp + '5' sf.remap1(ready_for_remap_2, ready_for_transp2, n1=w, d1=d2, verb=verb) sf.rm(out_transp1, verb) if rem2del_junk4: sf.rm(ready_for_remap_2, verb) sf.transp(ready_for_transp2, out, verb=verb) sf.rm(ready_for_transp2, verb) return ivlad.unix_success