Changeset 4325


Ignore:
Timestamp:
Feb 22, 2020 11:14:33 AM (20 months ago)
Author:
vondreele
Message:

small RMCProfile fixes
New auto spot mask version - works

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIimage.py

    r4324 r4325  
    14781478    for it,TTh in enumerate(TThs):
    14791479        band.mask = ma.masked_outside(TA,TTh,TTh+dtth).mask+tam
    1480         pcmax = np.percentile(band,prob)
     1480        pcmax = np.percentile(band.compressed(),prob)
    14811481        mband = ma.masked_greater(band,pcmax)
    1482         mean = ma.median(mband)
     1482        mean = ma.mean(mband)
    14831483        std = ma.std(mband)
    14841484        anom = ma.masked_greater((band-mean)/std,esdMul)
    14851485        mask ^= (anom.mask^band.mask)
    14861486        if not dlg is None:
    1487             dlg.Raise()
    14881487            GoOn = dlg.Update(it,newmsg='Processed 2-theta rings = %d'%(it))
    14891488            if not GoOn[0]:
  • trunk/GSASIIimgGUI.py

    r4324 r4325  
    19471947    spotSizer.Add(wx.StaticText(G2frame.dataWindow,label=' Select n*sigma spot rejection (n=1-10): '),0,WACV)
    19481948    spotSizer.Add(G2G.ValidatedTxtCtrl(G2frame.dataWindow,loc=data['SpotMask'],
    1949         key='esdMul',min=1,max=10,size=(40,25)),0,WACV)
     1949        key='esdMul',min=1.,max=10.,size=(40,25)),0,WACV)
    19501950    numPix = 0
    19511951    if data['SpotMask']['spotMask'] is not None:
  • trunk/GSASIIphsGUI.py

    r4323 r4325  
    53795379                    OutFile.close()
    53805380                    print('RMCProfile file %s read'%(pName+item))
    5381                 else:
    5382                     print('RMCProfile file %s not found'%(pName+item))
     5381#                else:
     5382#                    print('RMCProfile file %s not found'%(pName+item))
    53835383#total result plots
    53845384            Labels = {'_PDF1.csv':[r'$\mathsf{R,\AA}$','G(R)','RMCP G(R) for '],
     
    55025502                        numx,numy = odfData[:2]
    55035503                        G2plt.Plot3dXYZ(G2frame,int(numx),int(numy),odfData[2:],
    5504                             newPlot=False,Title='Bond %s-%s'%(bond[0],bond[1]),Centro=True)       
     5504                            newPlot=False,Title='Bond %s-%s'%(bond[0],bond[1]),Centro=True) 
     5505                    OutFile.close()
    55055506       
    55065507           
  • trunk/GSASIIplot.py

    r4319 r4325  
    52335233    if np.any(X) and np.any(Y) and np.any(Z):
    52345234        np.seterr(all='ignore')
    5235         Plot.plot_surface(X,Y,Z,rstride=1,cstride=1,color='g',linewidth=1)
     5235        try:
     5236            Plot.plot_surface(X,Y,Z,rstride=1,cstride=1,color='g',linewidth=1)
     5237        except:
     5238            pass
    52365239        xyzlim = np.array([Plot.get_xlim3d(),Plot.get_ylim3d(),Plot.get_zlim3d()]).T
    52375240        XYZlim = [min(xyzlim[0]),max(xyzlim[1])]
Note: See TracChangeset for help on using the changeset viewer.