Changeset 1189


Ignore:
Timestamp:
Jan 10, 2014 9:29:52 AM (8 years ago)
Author:
vondreele
Message:

change ring search defaults for stress/strainCalc
implement a fit all images for stress strain
image fit is on delta-d (was delta d2)

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIgrid.py

    r1184 r1189  
    9797
    9898[ wxID_MASKCOPY, wxID_MASKSAVE, wxID_MASKLOAD,wxID_NEWMASKSPOT,wxID_NEWMASKARC,wxID_NEWMASKRING,
    99   wxID_NEWMASKFRAME, wxID_NEWMASKPOLY,
     99    wxID_NEWMASKFRAME, wxID_NEWMASKPOLY,
    100100] = [wx.NewId() for item in range(8)]
    101101
    102 
    103102[ wxID_STRSTACOPY, wxID_STRSTAFIT, wxID_STRSTASAVE, wxID_STRSTALOAD,wxID_APPENDDZERO,
    104 ] = [wx.NewId() for item in range(5)]
     103    wxID_STRSTAALLFIT,
     104] = [wx.NewId() for item in range(6)]
    105105
    106106[ wxID_BACKCOPY,wxID_LIMITCOPY,wxID_SAMPLECOPY, wxID_BACKFLAGCOPY, wxID_SAMPLEFLAGCOPY,
     
    24242424        self.StrStaEdit.Append(help='Fit stress/strain data',
    24252425            id=wxID_STRSTAFIT, kind=wx.ITEM_NORMAL,text='Fit stress/strain')
     2426        self.StrStaEdit.Append(help='Fit stress/strain data for all images',
     2427            id=wxID_STRSTAALLFIT, kind=wx.ITEM_NORMAL,text='All image fit')
    24262428        self.StrStaEdit.Append(help='Copy stress/strain data to other images',
    24272429            id=wxID_STRSTACOPY, kind=wx.ITEM_NORMAL,text='Copy stress/strain')
  • trunk/GSASIIimage.py

    r1188 r1189  
    952952        V = 1.+np.sum(np.sum(E*calcFij(90.,phi,azm,th).T/1.e6,axis=2),axis=1)
    953953        dspc = dset*V
    954         return dspo**2-dspc**2
     954        return dspo-dspc
    955955       
    956956    names = ['e11','e12','e22']
  • trunk/GSASIIimgGUI.py

    r1188 r1189  
    13221322   
    13231323    def OnAppendDzero(event):
    1324         data['d-zero'].append({'Dset':1.0,'Dcalc':0.0,'pixLimit':10,'cutoff':10.0,
     1324        data['d-zero'].append({'Dset':1.0,'Dcalc':0.0,'pixLimit':10,'cutoff':1.0,
    13251325            'ImxyObs':[[],[]],'ImtaObs':[[],[]],'ImtaCalc':[[],[]],'Emat':[1.0,1.0,1.0]})
    13261326        UpdateStressStrain(G2frame,data)
     
    13921392
    13931393    def OnFitStrSta(event):
    1394         Masks = G2frame.PatternTree.GetItemPyData(
    1395             G2gd.GetPatternTreeItemId(G2frame,G2frame.Image, 'Masks'))
    13961394        Controls = G2frame.PatternTree.GetItemPyData(
    13971395            G2gd.GetPatternTreeItemId(G2frame,G2frame.Image, 'Image Controls'))
     
    14011399        G2plt.PlotExposedImage(G2frame,event=event)
    14021400        G2plt.PlotStrain(G2frame,data,newPlot=False)
     1401       
     1402    def OnAllFitStrSta(event):
     1403        TextList = [[False,'All IMG',0]]
     1404        Names = []
     1405        if G2frame.PatternTree.GetCount():
     1406            id, cookie = G2frame.PatternTree.GetFirstChild(G2frame.root)
     1407            while id:
     1408                name = G2frame.PatternTree.GetItemText(id)
     1409                Names.append(name)
     1410                if 'IMG' in name:
     1411                    TextList.append([False,name,id])
     1412                id, cookie = G2frame.PatternTree.GetNextChild(G2frame.root, cookie)
     1413            if len(TextList) == 1:
     1414                G2frame.ErrorDialog('Nothing to fit','There must some "IMG" patterns')
     1415                return
     1416            dlg = G2frame.CopyDialog(G2frame,'Stress/Strain fitting','Select images to fit:',TextList)
     1417            try:
     1418                if dlg.ShowModal() == wx.ID_OK:
     1419                    result = dlg.GetData()
     1420                    if result[0][0]:                    #the 'All IMG' is True
     1421                        result = TextList[1:]
     1422                        for item in result: item[0] = True
     1423                    for item in result:
     1424                        ifFit,name,id = item
     1425                        if ifFit:
     1426                            Controls = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,id, 'Image Controls'))
     1427                            StaCtrls = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,id, 'Stress/Strain'))
     1428                            id = G2gd.GetPatternTreeItemId(G2frame, G2frame.root, name)
     1429                            Npix,imagefile = G2frame.PatternTree.GetItemPyData(id)
     1430                            image = G2IO.GetImageData(G2frame,imagefile,True)
     1431                            G2img.FitStrSta(image,StaCtrls,Controls)
     1432                            G2plt.PlotStrain(G2frame,StaCtrls,newPlot=True)
     1433            finally:
     1434                dlg.Destroy()
     1435            print 'All images fitted'
    14031436       
    14041437    def SamSizer():
     
    14791512                value = min(10.0,max(0.5,float(Obj.GetValue())))
    14801513            except ValueError:
    1481                 value = 10.0
     1514                value = 1.0
    14821515            Obj.SetValue("%.1f"%(value))
    14831516            data['d-zero'][Indx[Obj.GetId()]]['cutoff'] = value
     
    15461579    G2frame.dataFrame.Bind(wx.EVT_MENU, OnAppendDzero, id=G2gd.wxID_APPENDDZERO)
    15471580    G2frame.dataFrame.Bind(wx.EVT_MENU, OnFitStrSta, id=G2gd.wxID_STRSTAFIT)
     1581    G2frame.dataFrame.Bind(wx.EVT_MENU, OnAllFitStrSta, id=G2gd.wxID_STRSTAALLFIT)
    15481582    G2frame.dataFrame.Bind(wx.EVT_MENU, OnCopyStrSta, id=G2gd.wxID_STRSTACOPY)
    15491583    G2frame.dataFrame.Bind(wx.EVT_MENU, OnLoadStrSta, id=G2gd.wxID_STRSTALOAD)
  • trunk/GSASIIplot.py

    r1188 r1189  
    23872387                return
    23882388            dsp = G2img.GetDsp(Xpos,Ypos,Data)
    2389             StrSta['d-zero'].append({'Dset':dsp,'Dcalc':0.0,'pixLimit':10,'cutoff':10.0,
     2389            StrSta['d-zero'].append({'Dset':dsp,'Dcalc':0.0,'pixLimit':10,'cutoff':1.0,
    23902390                'ImxyObs':[[],[]],'ImtaObs':[[],[]],'ImtaCalc':[[],[]],'Emat':[1.0,1.0,1.0]})
    23912391            R,r = G2img.MakeStrStaRing(StrSta['d-zero'][-1],G2frame.ImageZ,Data)
Note: See TracChangeset for help on using the changeset viewer.