def estimateNonUniformity(self):
        # TODO: base things on the input file size (and allow for overwritten defaults)
        impFields = []
        for i in range(len(self.inputs)):
            inputName   = self.inputFilenames[i]
            outFileBase = fh.removeBaseAndExtension(inputName) + "_nu_estimate.imp"
            outFileDir  = self.inputs[i].tmpDir
            outFile     = fh.createBaseName(outFileDir, outFileBase)
            impFields.append(outFile)
            cmd  = ["nu_estimate", "-clobber"]
            cmd += ["-distance", "8"]
            cmd += ["-iterations", "100"]
            cmd += ["-stop", "0.0001"]
            cmd += ["-fwhm", "0.15"]
            cmd += ["-shrink", "4"]
            cmd += ["-lambda", "5.0e-02"]
            if(self.masks):
                mask = self.masks[i]
                cmd += ["-mask", InputFile(mask)]
            cmd += [InputFile(inputName)]
            cmd += [OutputFile(outFile)]
            nu_estimate = CmdStage(cmd)
            nu_estimate.colour = "red"
            nu_estimate.setLogFile(LogFile(fh.logFromFile(self.inputs[i].logDir, outFile)))
            self.p.addStage(nu_estimate)
        self.impFields = impFields
more info... create the imp non uniformity estimation
dataset/ETHPy150Open Mouse-Imaging-Centre/pydpiper/atoms_and_modules/LSQ6.py/NonUniformityCorrection.estimateNonUniformity