77 77 | stokes, axis, mask, minpts, multifit, spxtype, spxest, |
78 78 | spxfix, div, spxsol, spxerr, |
79 79 | model, residual, wantreturn, |
80 80 | stretch, logresults, logfile, append, |
81 81 | sigma, outsigma |
82 82 | ): |
83 83 | casalog.origin('spxfit') |
84 84 | myia = image() |
85 85 | retval = None |
86 86 | try: |
87 - | if type(imagename) == list and len(imagename) > 1: |
88 - | myia = myia.imageconcat(outfile="", infiles=imagename, axis=axis, relax=True) |
87 + | if type(imagename) == list and len(imagename) > 1: |
88 + | myia = myia.imageconcat(outfile="", infiles=imagename, axis=axis, relax=True) |
89 + | else: |
90 + | if type(imagename) == list and len(imagename) == 1: |
91 + | imagename = imagename[0] |
92 + | if (not myia.open(imagename)): |
93 + | raise RuntimeError("Cannot create image analysis tool using " + str(imagename)) |
94 + | sigmacopy = sigma |
95 + | if type(sigma) == list and type(sigma) == str: |
96 + | if len(sigma) == 1: |
97 + | sigmacopy = sigma[0] |
90 - | if type(imagename) == list and len(imagename) == 1: |
91 - | imagename = imagename[0] |
92 - | if (not myia.open(imagename)): |
93 - | raise Exception("Cannot create image analysis tool using " + str(imagename)) |
94 - | sigmacopy = sigma |
95 - | if type(sigma) == list and type(sigma) == str: |
96 - | if len(sigma) == 1: |
97 - | sigmacopy = sigma[0] |
98 - | else: |
99 - | sigia = myia.imageconcat(outfile="", infiles=sigma, axis=axis, relax=True) |
100 - | sigmacopy = sigia.getchunk() |
101 - | retval = myia.fitprofile( |
102 - | box=box, region=region, chans=chans, |
103 - | stokes=stokes, axis=axis, mask=mask, |
104 - | minpts=minpts, ngauss=0, multifit=multifit, |
105 - | spxtype=spxtype, spxest=spxest, spxfix=spxfix, |
106 - | div=div, model=model, residual=residual, |
107 - | stretch=stretch, logresults=logresults, |
108 - | spxsol=spxsol, spxerr=spxerr, logfile=logfile, |
109 - | append=append, |
110 - | sigma=sigmacopy, outsigma=outsigma |
111 - | ) |
112 - | except Exception as instance: |
113 - | casalog.post( str( '*** Error ***') + str(instance), 'SEVERE') |
114 - | retval = None |
115 - | myia.done() |
116 - | if (wantreturn): |
117 - | return retval |
118 - | else: |
119 - | if (retval): |
120 - | del retval |
121 - | return None |
122 - | |
123 - | |
124 - | |
99 + | sigia = myia.imageconcat(outfile="", infiles=sigma, axis=axis, relax=True) |
100 + | sigmacopy = sigia.getchunk() |
101 + | retval = myia.fitprofile( |
102 + | box=box, region=region, chans=chans, |
103 + | stokes=stokes, axis=axis, mask=mask, |
104 + | minpts=minpts, ngauss=0, multifit=multifit, |
105 + | spxtype=spxtype, spxest=spxest, spxfix=spxfix, |
106 + | div=div, model=model, residual=residual, |
107 + | stretch=stretch, logresults=logresults, |
108 + | spxsol=spxsol, spxerr=spxerr, logfile=logfile, |
109 + | append=append, |
110 + | sigma=sigmacopy, outsigma=outsigma |
111 + | ) |
112 + | finally: |
113 + | myia.done() |
114 + | if (wantreturn): |
115 + | return retval |
116 + | else: |
117 + | if (retval): |
118 + | del retval |
119 + | return |