Moeck P, DeStefano P
Nano-Crystallography Group, Department of Physics, Portland State University, P.O. Box 751, Portland, OR 97201 USA.
Adv Struct Chem Imaging. 2018;4(1):5. doi: 10.1186/s40679-018-0051-z. Epub 2018 Mar 28.
Three different algorithms, as implemented in three different computer programs, were put to the task of extracting direct space lattice parameters from four sets of synthetic images that were per design more or less periodic in two dimensions (2D). One of the test images in each set was per design free of noise and, therefore, genuinely 2D periodic so that it adhered perfectly to the constraints of a Bravais lattice type, Laue class, and plane symmetry group. Gaussian noise with a mean of zero and standard deviations of 10 and 50% of the maximal pixel intensity was added to the individual pixels of the noise-free images individually to create two more images and thereby complete the sets. The added noise broke the strict translation and site/point symmetries of the noise-free images of the four test sets so that all symmetries that existed per design turned into pseudo-symmetries of the second kind. Moreover, motif and translation-based pseudo-symmetries of the first kind, a.k.a. genuine pseudo-symmetries, and a metric specialization were present per design in the majority of the noise-free test images already. With the extraction of the lattice parameters from the images of the synthetic test sets, we assessed the robustness of the algorithms' performances in the presence of both Gaussian noise and pre-designed pseudo-symmetries. By applying three different computer programs to the same image sets, we also tested the reliability of the programs with respect to subsequent geometric inferences such as Bravais lattice type assignments. Partly due to per design existing pseudo-symmetries of the first kind, the lattice parameters that the utilized computer programs extracted in their default settings disagreed for some of the test images even in the absence of noise, i.e., in the absence of pseudo-symmetries of the second kind, for any reasonable error estimates. For the noisy images, the disagreement of the lattice parameter extraction results from the algorithms was typically more pronounced. Non-default settings and re-interpretations/re-calculations on the basis of program outputs allowed for a reduction (but not a complete elimination) of the differences in the geometric feature extraction results of the three tested algorithms. Our lattice parameter extraction results are, thus, an illustration of Kenichi Kanatani's dictum that no extraction algorithm for geometric features from images leads to results because they are all aiming at an intrinsically impossible task in all real-world applications (Kanatani in Syst Comput Jpn 35:1-9, 2004). Since 2D-Bravais lattice type assignments are the natural end result of lattice parameter extractions from more or less 2D-periodic images, there is also a section in this paper that describes the intertwined metric relations/holohedral plane and point group symmetry hierarchy of the five translation symmetry types of the Euclidean plane. Because there is no definitive lattice parameter extraction algorithm, the outputs of computer programs that implemented such algorithms are also not definitive. Definitive assignments of higher symmetric Bravais lattice types to real-world images should, therefore, not be made on the basis of the numerical values of extracted lattice parameters and their error bars. Such assignments require (at the current state of affairs) arbitrarily set thresholds and are, therefore, always so that they cannot claim objective definitiveness. This is the essence of Kenichi Kanatani's comments on the vast majority of computerized attempts to extract symmetries and other hierarchical geometric features from noisy images (Kanatani in IEEE Trans Pattern Anal Mach Intell 19:246-247, 1997). All there should be instead for noisy and/or genuinely pseudo-symmetric images are rankings of the relative likelihoods of classifications into higher symmetric Bravais lattice types, Laue classes, and plane symmetry groups.
三种不同的算法,分别在三个不同的计算机程序中实现,被用于从四组合成图像中提取直接空间晶格参数。这些合成图像在二维(2D)上按设计或多或少具有周期性。每组中的一张测试图像按设计无噪声,因此是真正的二维周期性图像,完全符合布拉菲晶格类型、劳厄类和平面对称群的约束。将均值为零、标准差分别为最大像素强度的10%和50%的高斯噪声分别添加到无噪声图像的各个像素上,以创建另外两张图像,从而完成图像集。添加的噪声打破了四个测试集无噪声图像的严格平移和位点/点对称性,使得按设计存在的所有对称性都变成了第二类伪对称性。此外,在大多数无噪声测试图像中,按设计已经存在第一类基于图案和平移的伪对称性,即真正的伪对称性,以及一种度量特化。通过从合成测试集的图像中提取晶格参数,我们评估了算法在存在高斯噪声和预先设计的伪对称性情况下性能的稳健性。通过将三个不同的计算机程序应用于相同的图像集,我们还测试了程序在后续几何推断(如布拉菲晶格类型分配)方面的可靠性。部分由于按设计存在的第一类伪对称性,即使在没有噪声(即没有第二类伪对称性)的情况下,对于任何合理的误差估计,所使用的计算机程序在其默认设置下提取的晶格参数对于某些测试图像也不一致。对于有噪声的图像,算法提取的晶格参数结果的不一致通常更为明显。基于程序输出的非默认设置和重新解释/重新计算,使得三个测试算法在几何特征提取结果上的差异有所减少(但没有完全消除)。因此,我们的晶格参数提取结果说明了金谷谦一的格言:从图像中提取几何特征的任何算法都不会得出确定的结果,因为在所有实际应用中它们都旨在完成一项本质上不可能的任务(金谷谦一,《系统计算日本》35:1 - 9,2004)。由于二维布拉菲晶格类型分配是从或多或少二维周期性图像中提取晶格参数的自然最终结果,本文还有一部分描述了欧几里得平面五种平移对称类型相互交织的度量关系/全对称平面和点群对称层次。由于没有确定的晶格参数提取算法,实现此类算法的计算机程序的输出也不是确定的。因此,不应基于提取的晶格参数的数值及其误差条对现实世界图像进行更高对称布拉菲晶格类型的确定分配。此类分配需要(在当前情况下)任意设置阈值,因此总是主观的,无法声称具有客观确定性。这就是金谷谦一对绝大多数从有噪声图像中提取对称性和其他层次几何特征的计算机化尝试的评论的本质(金谷谦一,《IEEE模式分析与机器智能汇刊》19:246 - 247,1997)。对于有噪声和/或真正伪对称的图像,应该做的只是对更高对称布拉菲晶格类型、劳厄类和平面对称群分类的相对可能性进行排序。