从CrysAlisPro衍射图片获取二维散射数据

casjxm

单晶衍射测试的获得的数据精度较高，一些有缺陷或无序的的晶体也会产生明显的散射信息，这些散射信息可用于缺陷结构的解析，但与一维散射数据（如PDF等）不同的是，单晶衍射测试获得的数据为二维照片或重构出来的三维ewald倒空间数据，下面介绍如何从理学CrysAlisPro（CAP）软件产生的二维照片中读取散射数据。
CAP有个功能是获得二维Unwarping Precesssion images，包含有完整的bragg与散射信息，默认是保存在unwarp文件夹中img格式的文件中，img文件为非标准的图片格式，直接读取比较麻烦，可以通过CAP界面左下角的照片转换成通用的jpg格式，然后就可以用通用代码处理了。
注意img或jpg中的数据是直角坐标系表示的强度信息，因此需要通过取向矩阵转换成倒格子晶体学坐标系下的强度数据。
晶体学坐标与直角坐标的转换参考：https://www.matstr.com/forum.php ... 1&extra=#pid880


读取jpg文件中的散射强度（matlab）：
jpg包含有3个通道（RGB）的散射数据，但每个RGB值是一样的，实际上只需要一个通道的数据就可以了。
读取：
img = imread('test.jpg');

size(img)
img(1,1)

转换成灰度图像：
img_gray = rgb2gray(img);

可视化验证：

img = imread('test.jpg');
subplot(1,2,1);
imshow(img);
title('原始图像');
if ndims(img) == 3

    img_gray = rgb2gray(img);
else
    img_gray = img;
end
subplot(1,2,2);
imshow(img_gray);
title('灰度化结果');