Changeset 4593 for trunk


Ignore:
Timestamp:
Oct 13, 2020 2:16:15 PM (14 months ago)
Author:
vondreele
Message:

fix image plotting issues with 2-theta limits, arc & ring masks & azimuth limits with new GetDetectorXY code

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIimage.py

    r4589 r4593  
    507507    xyz -= np.array([0.,0.,dist])                 #translate back
    508508    xyz = np.inner(xyz,iMN)
    509     return xyz[:2]+cent
     509    return np.squeeze(xyz)[:2]+cent
    510510       
    511511def GetDetectorXY2(dsp,azm,data):
     
    638638    distsq = data['distance']**2
    639639    G = ((dx-x0x)**2+(dy-x0y)**2+distsq)/distsq       #for geometric correction = 1/cos(2theta)^2 if tilt=0.
    640     return np.array([tth,azm,G,dsp])
     640    return tth,azm,G,dsp
    641641   
    642642def GetTth(x,y,data):
  • trunk/GSASIIplot.py

    r4589 r4593  
    68226822    Azm = np.linspace(*aR)
    68236823    for azm in Azm:
    6824         XY = G2img.GetDetectorXY(Dsp(angI,wave),azm,Data)
     6824        XY = G2img.GetDetectorXY(Dsp(np.squeeze(angI),wave),azm,Data)
    68256825        if XY is not None:
    68266826            xy1.append(XY)      #what about hyperbola
    6827         XY = G2img.GetDetectorXY(Dsp(angO,wave),azm,Data)
     6827        XY = G2img.GetDetectorXY(Dsp(np.squeeze(angO),wave),azm,Data)
    68286828        if XY is not None:
    68296829            xy2.append(XY)      #what about hyperbola
     
    72487248                    azm[0] += off
    72497249                    azm[1] -= off
    7250                 (x1,y1),(x2,y2) = ComputeArc(angI,angO,wave,*azm)
     7250                (x1,y1),(x2,y2) = ComputeArc(np.squeeze(angI),np.squeeze(angO),wave,*azm)
    72517251                pI,pO,pL,pU = G2frame.arcList[pick.itemNumber]
    72527252                pI.set_data((x2,y2))
     
    77177717                xyO = []
    77187718                arcxO = []
    7719                 arcxI = []
    77207719                for azm in Azm:
    77217720                    xy = G2img.GetDetectorXY(dspO,azm,Data)
     
    77267725                    arcxO,arcyO = xyO.T               
    77277726                    Plot.plot(arcxO,arcyO,picker=3,label='Otth')
    7728             if ellO and ellI and len(arcxO) and len(arcxI):
     7727            if ellO and ellI and len(arcxO):
    77297728                Plot.plot([arcxI[0],arcxO[0]],[arcyI[0],arcyO[0]],picker=3,label='Lazm')
    77307729                Plot.plot([arcxI[-1],arcxO[-1]],[arcyI[-1],arcyO[-1]],picker=3,label='Uazm')
     
    77347733                    cake += delAzm/2.
    77357734                ind = np.searchsorted(Azm,cake)
    7736                 if len(arcxO) and len(arcxI):
     7735                if len(arcxO):
    77377736                    Plot.plot([arcxI[ind],arcxO[ind]],[arcyI[ind],arcyO[ind]],color='k',dashes=(5,5))
    77387737        if 'linescan' in Data and Data['linescan'][0] and G2frame.GPXtree.GetItemText(G2frame.PickId) in ['Image Controls',]:
     
    77947793        for iarc,arc in enumerate(Masks['Arcs']):      # drawing arc masks
    77957794            if arc:
    7796                 tth,azm,thick = arc           
     7795                tth,azm,thick = arc
     7796                azm = np.squeeze(azm)
    77977797                wave = Data['wavelength']
    77987798                (x1,y1),(x2,y2) = ComputeArc(tth-thick/2.,tth+thick/2.,wave,azm[0],azm[1])
Note: See TracChangeset for help on using the changeset viewer.