Changeset 2731 for trunk/GSASIIimage.py

Timestamp:

Mar 1, 2017 9:08:17 AM (8 years ago)

Author:

vondreele

Message:

add save of stress/strain ring intensities as MRD in txt file
calculate stress/strain ring intensities from 5x5 pixel blocks (was 3x3)
show sig(MRD) instead of var(MRD) in ring intensities

File:

• trunk/GSASIIimage.py (modified) (3 diffs)

trunk/GSASIIimage.py ¶

@@ -214 +214 @@
     ypix = int(Ypix)
     if not w:
-        ZMax = np.sum(image[ypix-1:ypix+1,xpix-1:xpix+1])
+        ZMax = np.sum(image[ypix-2:ypix+2,xpix-2:xpix+2])
         return xpix,ypix,ZMax,0.0001
     if (w2 < xpix < sizex-w2) and (w2 < ypix < sizey-w2) and image[ypix,xpix]:
@@ -1030 +1030 @@
             ringint = np.array([float(Image[int(x*scalex),int(y*scaley)]) for y,x in np.array(ringxy)[:,:2]])
             ringint /= np.mean(ringint)
-            print ' %s %.3f %s %.3f'%('d-spacing',ring['Dcalc'],'var(MRD):',np.var(ringint))
+            print ' %s %.3f %s %.3f'%('d-spacing',ring['Dcalc'],'sig(MRD):',np.sqrt(np.var(ringint)))
             RingsAI.append(np.array(zip(ringazm,ringint)).T)
-#            GSASIIpath.IPyBreak()
     return RingsAI
     
@@ -1056 +1055 @@
         else:                       #compression
             ring['Dcalc'] = dmin+3.*(dmax-dmin)/4.
-#        ring['Dcalc'] = np.mean(ring['ImtaCalc'][0])
 
 def calcFij(omg,phi,azm,th):