Changeset 2613
- Timestamp:
- Jan 6, 2017 9:34:29 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIimage.py
r2604 r2613 30 30 import GSASIIpwd as G2pwd 31 31 import GSASIIspc as G2spc 32 import GSASIImath as G2mth 32 33 33 34 # trig functions in degrees … … 1141 1142 cent[1,1] = mag 1142 1143 msk = np.array(Image[ind[0]-1:ind[0]+2,ind[1]-1:ind[1]+2]) 1143 msk = m sk-cent1144 if mag > 1. 5*np.mean(msk):1144 msk = ma.array(msk,mask=(msk==mag)) 1145 if mag > 1.33*ma.mean(msk): 1145 1146 jndx.append([ind[1]+.5,ind[0]+.5]) 1147 print 'Spots found: ',len(jndx) 1146 1148 if len(jndx) > 100: 1147 1149 txt = 'Found: %d. Too many spots found; are rings spotty?'%(len(jndx)) … … 1150 1152 peaks = jndx*pixelSize/1000. 1151 1153 tth = GetTth(peaks.T[0],peaks.T[1],Controls) 1152 histo = np.histogram(tth,2500) 1154 Peakarray = np.vstack((tth,peaks.T)).T 1155 Peakarray = np.array(G2mth.sortArray(Peakarray,0)) #now in 2theta order 1153 1156 Points = np.ones((peaks.shape[0],3)) 1154 Points[:,:2] = peaks1157 Points[:,:2] = Peakarray[:,1:] 1155 1158 Masks['Points'] = Points 1156 1159 return None
Note: See TracChangeset
for help on using the changeset viewer.