def radialMaxLevel(self, rmin, rmax, Nbin): "Compute the maximum (absolute) level for the image and return the array of radius (arcsec) and values" dr = (rmax - rmin) / Nbin radius = np.arange(Nbin) * dr + rmin radMaxLevel = np.zeros(Nbin) ## loop over the pixels ia.open(self.beamImage) for i in range(self.Nx): for j in range(self.Ny): xx = (i - self.refx) * self.dx yy = (j - self.refy) * self.dy rr = sqrt(xx * xx + yy * yy) indexR = floor((rr - rmin) / dr) if indexR >= 0 and indexR < Nbin: val = ia.pixelvalue([i, j])['value']['value'] if abs(val) > abs(radMaxLevel[indexR]): radMaxLevel[indexR] = val ia.close() return (radius, radMaxLevel)
def radialMaxLevel(self,rmin,rmax,Nbin): "Compute the maximum (absolute) level for the image and return the array of radius (arcsec) and values" dr=(rmax-rmin)/Nbin radius=np.arange(Nbin)*dr+rmin radMaxLevel=np.zeros(Nbin) ## loop over the pixels ia.open(self.beamImage) for i in range(self.Nx): for j in range(self.Ny): xx=(i-self.refx)*self.dx yy=(j-self.refy)*self.dy rr=sqrt(xx*xx+yy*yy) indexR=floor((rr-rmin)/dr) if indexR >=0 and indexR < Nbin: val=ia.pixelvalue([i,j])['value']['value'] if abs(val) > abs(radMaxLevel[indexR]): radMaxLevel[indexR]=val ia.close() return(radius,radMaxLevel)