This paper proposes parallel implementation of the image informative segments extraction method. The images are segmented in the spatial spectrum domain. Median energy in each selected segment is viewed as an area. For time saving purpose parallel implementation was developed for the areas calculation phase. The developed software implementation was tested on the high performance multicore computing system. Currently computer processing of medical diagnostic images is one of the vital research tools and a way to improve efficiency of early detection of various diseases. Change of the body fluids composition is one of the information-bearing health condition areas. Metabolic change that occurs due to pathological conditions affects the fluid composition; there are numerous changes in the molecular composition of tissues and body fluids. Converting the fluid from one phase state to another is one of the ways to detect such changes. Crystallization is one of the most convenient methods to change the fluid phase. Crystal properties modification is caused by changed physical and chemical properties of a body fluid. The investigation of these properties is the crucial problem of crystal analysis [1]. In medicine, studied crystallograms are the structures formed by salt crystallization caused by body fluid drying. In clinical practice the crystallogram analysis is based on their images. It is not always possible to visually identify changes in such key crystallogram parameters as predominant bar direction, bar density etc., which contribute to major pathologic signs. Quantified analysis and objectivity are among computer analysis advantages. The information contained in the image is structurally excessive. It is known that if the parallel bars of certain direction were predominant on the original image, the bars of the same direction would dominate the Fourier transform of an original image. This property can be used to analyze crystallograms [2, 3] and other images of branched structure. The developed method based on discriminatory analysis algorithm is applied to generate the informative areas set, which is used in this paper to identify the characteristics of the initial crystallogram images. In the article, we propose a parallel implementation of the method to speed up computations.
The methods based on the discriminative analysis algorithm proved to be a good solution to form up new problem-specified
areas [
The discriminative analysis is used to eliminate the correlation between the areas and consequently to reduce the size of the areas set. The usage of this algorithm allows, on the one hand, to maintain the informativeness of the feature set for classification and, on the other hand, to reduce the number of areas to apply less complicated classification methods and to reduce the classification error value.
An individual separability criterion was calculated for each area based on equation: J=tr((T)-1B), where T = B + W, B is between-group scattering matrix, W is intragroup scattering matrix.
The informative areas set is further generated in the following manner: The areas are ranged in the order of decreasing of individual separability criteria values. The initial areas set consists of a feature with the largest criteria. Classification is carried out. Then the area with the next value of separability criterion is added to the set. Classification based on the new set of areas is carried out. Repeat item 3 until all areas are included into the set.
The informative set of areas is the one with classification results yielding the minimum error value. The classification error that defines the number of cases with the classifier acquiring incorrect value is calculated from the equation: ε = (m/n) · 100%, – where m is the number of classification errors, n is the total number of images tested.
The entire algorithm of informative areas generation may be presented as a diagram shown in Figure 1. The analysis of the algorithm structure showed that the first stage of areas computing has the maximum computational complexity.
As mentioned above most of the algorithm run time is spent on computing the areas. This stage includes pre-processing of the learning sample, generation of spatial spectrum for each learning image and calculation of the area values. The number of calculations can be reduced if one take into account that the spatial spectrum image is symmetrical relative to its center, and for this reason the areas calculation may and must be performed only on the half of the spatial spectrum image. In order to speed up the algorithm run time this paper will apply task division. This paper will not use an MPI technology but will apply task distribution by threads, that is why the way of splitting the image into tasks is not important. In case of sample pre-processing and spatial spectrum image generation a single image is sent to each thread. Then at the areas calculation stage each thread reveives a separate element - an image segment calculated on the basis of the proposed segmentation. The next task is sent to the thread as soon as it completes handling of the element.
Application of this segmentation pattern allows substantial run-time saving at the first stage of the algorithm. During this resource-intensive stage such method of tasks segmentation by threads has permitted to achieve the three-fold acceleration when four threads were used. The curve in Figure 2 shows the relation between acceleration and the number of areas.
After classification for all possible segmentations, areas were selected using the developed area selection algorithm. Figure 3 shows the relation between classification error and the number of areas taken in the descending order of their separability criterion in case of classification with 2-nearest neighbors methods to split into 4 sectors and 8 rings r s ebm fo tcer
o uN s 5. Conclusion
This paper presented the method to generate a set of informative local areas of a spatial spectrum in order to classify medical crystallogram images, and an option of parallel implementation of such method. In addition, the research included experimental testing of the developed software implementation, which demonstrated that parallel algorithm implementation provided almost three-fold acceleration at the areas calculation stage. The next step is to implement parallel computing at the smart area analysis stage in order to improve the informative area set calculation speed.
This work was partially supported by the Ministry of education and science of the Russian Federation in the framework of the implementation of the Program of increasing the competitiveness of SSAU among the world’s leading scientific and educational centers for 2013-2020 years; by the Russian Foundation for Basic Research grants (# 15-29-03823, # 15-29-07077, # 16-41
630761; # 16-29-11698, # 17-01-00972); by the ONIT RAS program # 6 “Bioinformatics, modern information technologies and mathematical methods in medicine” 2017