Changeset 89


Ignore:
Timestamp:
Jun 18, 2010 9:22:13 AM (12 years ago)
Author:
vondreel
Message:

make arc masks bigger

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIimage.py

    r88 r89  
    469469    tax,tay = np.dsplit(TA,2)    #azimuth, 2-theta
    470470    for tth,thick in rings:
    471         tam = ma.mask_or(tam.flatten(),ma.getmask(ma.masked_inside(tay.flatten(),tth-thick/2.,tth+thick/2.)))
     471        tam = ma.mask_or(tam.flatten(),ma.getmask(ma.masked_inside(tay.flatten(),max(0.01,tth-thick/2.),tth+thick/2.)))
    472472    for tth,azm,thick in arcs:
    473         tam = ma.mask_or(tam.flatten(),ma.getmask(ma.masked_inside(tay.flatten(),tth-thick/2.,tth+thick/2.))* \
     473        tam = ma.mask_or(tam.flatten(),ma.getmask(ma.masked_inside(tay.flatten(),max(0.01,tth-thick/2.),tth+thick/2.))* \
    474474            ma.getmask(ma.masked_inside(tax.flatten(),azm[0],azm[1])))
    475475    taz = ma.masked_greater(ma.masked_less(image,Zlim[0]),Zlim[1]).flatten()
  • trunk/GSASIIplot.py

    r88 r89  
    845845        wave = Data['wavelength']
    846846        x1,y1 = np.hsplit(np.array(G2img.makeIdealRing(G2img.GetEllipse(Dsp(tth+thick/2.,wave),Data),azm)),2)
    847         x2,y2 = np.hsplit(np.array(G2img.makeIdealRing(G2img.GetEllipse(Dsp(tth-thick/2.,wave),Data),azm)),2)
     847        x2,y2 = np.hsplit(np.array(G2img.makeIdealRing(G2img.GetEllipse(Dsp(max(0.01,tth-thick/2.),wave),Data),azm)),2)
    848848        self.arcList.append([Plot.plot(x1,y1,'r',picker=3),iarc,'o'])           
    849849        self.arcList.append([Plot.plot(x2,y2,'r',picker=3),iarc,'i'])
Note: See TracChangeset for help on using the changeset viewer.