|
1.INTRODUCTIONIn recent years, the application of CT equipment in medical and industrial fields has become more and more extensive, and CT quality control has attracted much attention as a key link to ensure the performance and image quality of CT equipment. CT equipment quality control mainly refers to the use of professional testing modules and equipment to test the technical indexes and performance parameters of the CT equipment, in order to assess the quality of the CT equipment and the practicability of whether it meets the clinical needs, and to ensure that the CT equipment meets the medical needs. Quality control of CT is the process of evaluating the quality and utility of CT equipment to ensure that it meets clinical needs. Quality control of CT is an important part of CT imaging. It not only concerns the accuracy of medical diagnosis, but also directly affects the treatment effect and life safety of patients. Under the premise of safe radiation dose, CT quality control can ensure that the equipment provides sufficiently clear and accurate diagnostic information to help doctors make accurate judgment and treatment decisions. At the same time, effective quality control can also minimize the potential radiation risk to patients and reduce the damage caused by radiation to the patient’s body[1]. Since the world’s first commercial CT was applied to the clinic in 1972, experts and scholars have recognized the importance of CT quality control. Since then, Europe, the United States and other developed countries have invested a large amount of research resources in this field, and have continuously optimized the techniques and methods of quality control. For example, Germany has achieved remarkable results in the accuracy of detection equipment, the United States is in a leading position in the research of radiation dose control, and Japan has made outstanding contributions to the improvement of image clarity and resolution. China’s CT equipment quality control work started late, and the technology is relatively behind Germany, the United States, Japan and other developed countries. However, with the introduction of the first CT equipment in 1979, China’s research and practice in related fields has gradually increased. In terms of testing technology, from the initial dependence on foreign technology, to today’s independent research and development of a series of high-precision and high-efficiency testing methods; testing module has also developed from a simple copy of the development of independent innovation to meet the needs of different types of CT equipment; supervision and management, from the initial lack of norms to the establishment of a comprehensive regulatory system to ensure that the quality control of CT equipment can be effectively carried out. supervision and management, from the initial lack of standardization to the establishment of a perfect regulatory system to ensure the effective implementation of CT equipment quality control. The quality control of CT equipment in China has developed rapidly, and the detection technology, detection model, supervision and management have been continuously improved and perfected [4]. At home and abroad, many studies have been devoted to CT quality control automation, and certain results have been achieved, but there are still some limitations. In foreign countries, the relevant research mainly focuses on the use of advanced technology and equipment to improve the automation level of CT quality control. For example, some studies have adopted high-precision detection modules and instruments to obtain more accurate detection data; and some studies have been devoted to the development of intelligent detection systems that can automatically identify and analyze problems in CT images. For example, in 2020, Polaris Imaging (NSI) introduced the robotiX fully automated X-ray CT scanning system [5], which integrates industrial robotics with standard digital imaging and computed tomography X-ray equipment to improve the automation capabilities and flexibility of industry CT systems. However, these methods tend to be costly and may be limited by a variety of factors in practical applications. The current means of CT image detection are manually detected by human beings, which is time-consuming and laborious and may have subjective errors. To address this deficiency, this study makes full use of the powerful image processing function of MATLAB and combines it with APP Designer to design an independent APP to realize the automation of CT quality control. Compared with the existing automatic methods, this method has the following features and advantages:
2.CODE IMPLEMENTATION OF MATLAB-BASED DETECTION AUTOMATIONAt present, CT quality control process in China is mainly based on the implementation of two standards, respectively, “JJG961-2017 medical diagnostic spiral computed tomography device (CT) X-ray radiation source verification procedures” [7] and “GB17589-2011 X-ray computed tomography device quality assurance testing specifications” [8]. According to China’s current “Radiological Diagnosis and Treatment Management Regulations” regulatory requirements, according to the standard GB17589-2011 [8], the quality control of CT equipment includes three aspects: acceptance testing before use, annual condition testing and daily stability testing. 2.1Testing programThe testing items include diagnostic bed positioning accuracy, positioning light accuracy, scanning frame inclination accuracy, reconstruction layer thickness deviation, CT-weighted dose index, consistency between images, CT number, uniformity, noise, CT number linearity, spatial resolving power (rate), contrast resolving power (rate) [8]. In this paper, several items of CT number, noise, and layer thickness of water are automated.
2.2Basis of testingThe detection module used in this paper is CatPhan500 phantom (shown in Figure 1), CatPhanTM500 includes four detection modules: CTP401: layer thickness, CT number linearity and contrast scale, CTP528: high contrast resolution, CTP515: low contrast resolution, CTP486: field uniformity and noise[11]. In this paper, we specifically utilize CTP401 and CTP486, and the detection indexes are CT number, noise and layer thickness of water, respectively. CT quality control is detected when the CT number of the water in the homogeneous body water model (Figure 2), by the CT number of the formula can be obtained, the CT number of water should be 0, so the detected performance requirements: the CT number of the water should be within ± 4 (acceptance testing); ± 6 range (state testing). 2.3Code implementation of the CT number of water and the noise sectionFirstly, CT images are read and preprocessed, CT image data are usually saved in DICOM file format by medical devices or imaging software. The CT image data (DICOM format) was read using the dicomread function (Digital Imaging and Communications in Medicine) [6] and grayscaled using imshow (I,[]). In this paper, obtaining ROI region link belongs to a semi-automatic process, first of all, use plot and poly2mask function to convert ROI region to region mask to realize the extraction of ROI region. When using the plot function to draw a circle, the center of the circular region is obtained by calculating the center of mass, i.e., the coordinates of the center of mass are found by setting a threshold and performing a binary loop on the image data to find the center of mass, thus realizing the center of mass, and then the center of mass is used as the center of the circle to successfully draw a 500 pixel-point circle. The ROI region in this paper refers to the circular region of 500 pixel points in the center of the image required in the CT number of water and noise detection project, as well as the small window region where the line segments need to be cut out for convenient processing in layer thickness detection. 2.4Code implementation of the layer thickness detection partLayer thickness detection is to measure the length of the metal diagonal line in the image, according to the requirements of the half-value width and height, you need to find the highest CT number of the line L2, as well as the average CT number of the background L1. In this paper, we use the imbinarize function for the grey-scale image binarization, and then use the bwlabeln function for the connectivity component to label, to find the background region and diagonal area, and then calculate L1 and L2. The automatic setting of the window position is L=(L1+L2)/2, and the window width is adjusted to the minimum. It can make the line segments to be measured can make the line segments to be measured (the four line segments in Fig. 3) partially visible (as in Fig. 3), omitting the step of conventional manual detection to find the optimal window position, and also more accurately through the pixel points to calculate the length directly [12]. Marking the connected components is to find each connected region in the binary map and mark it with some kind of notation, so that the number and location of connected regions can be determined. In the layer thickness detection part, it is necessary to find the position of the line segment to be measured, which can be regarded as a connected component, and the rest of the interfering scatter is another region, by this way, the line segment and the interfering scatter can be separated in the matrix, which is the basis for further point selection and calculation. In the layer thickness detection part of this thesis, the bwlabeln function is applied to a part of the intercepted area to create a labeling matrix using 4-connected objects. The effect of using this is shown in Figure 4. Labeling after the acquisition of the coordinates of the target region, as shown in Figure 4, connected region labeling after the line segment region and other regions in the matrix of different numbers, this step needs to do is to find the number of regions and get the endpoint coordinates, the next step can be used to calculate the distance only with these coordinates. In this paper, in order to achieve this purpose, we first use histc function to find out the number of line segment region represented by the number, screening out the most consecutive occurrences of the number of times, and then use find function with the value of the same time restrictions to find the line number of the largest and smallest of the two points as well as to get the coordinates. 3.INSPECTION SOFTWARE DESIGN BASED ON APP DESIGNER3.1Brief description of Matlab APP DesignerAPP Designer is a MATLAB application design platform that does not require users to have a background in software design and high-end computer design, but only need to drag and drop visual components to realize the graphical user interface (GUI) design layout, followed by controls, layouts, callback functions, and other features to quickly build function and data-driven application program interface, you can output a simple APP. 3.2Sharing data within APP DesignerEach callback of APP Designer is relatively independent, the program and data behind each button are not connected, and in the program design of this paper, there is a need to call the data of other callback modules, for example, the detection process requires the pictures in the coordinate system, and the measurement value of CT number is required for the calculation of noise. Using attributes is the best way to share data within an App. If you want to share data that needs to be accessed by an intermediate result or multiple callbacks, you can define public or private attributes to store the data. This thesis uses the addition of private properties to enable data sharing between individual callbacks. 3.3General framework design of CT quality control softwareThe CT quality control inspection software is divided into two modules containing three inspection items, and the specific design flow is shown in Figure 5. The specific APP design interface is shown in Figure 6, which mainly consists of the coordinate area and the three detection project buttons and their data output text boxes, in addition to placing the button to select the image. In the use phase, more emphasis is placed on the output of the results so the amount of process is reduced to show a more concise interface and simpler operation. When measuring the CT number and noise of water just click the button to automatically select the ROI and output the result. For layer thickness detection, you only need to interactively click on the desired area to get the vertex coordinates to intercept the image for the next step of calculation. 3.4Main function modules and realization
4.RESULT AND DISCUSSIONIn this experiment, the measurement results of manual testing and the developed CT quality control software were compared, and the results are shown in Table 1. Table 1.Comparison of measurement results between the two methods
As can be seen from Table 1, the results of the detection software in this thesis are very similar to the manual detection results, and the results are informative, which can solve the error problems in manual detection, and compared to the inevitable errors in the manual detection process, such as inaccurate selection of the ROI range in the process of CT number detection, and accidental errors brought about by the measurement of the length of the layer thickness in the process of layer thickness detection, the automated detection software can solve all these problems, and the faster and more efficient. Therefore, this study can be used for daily CT quality control inspection. 5.CONCLUSIONIn this thesis, with the standard of CT quality control, the powerful image processing ability of MATLAB software and the environmental support of APP Designer, we have designed the convenient and simple CT quality control inspection software, meanwhile, this thesis has designed all the inspection items into a whole software through page design and callback routine, which is hoped to be used for CT quality control inspection. Meanwhile, it can be improved in the following aspects in the future project research:
ACKNOWLEDGMENTSThis research was funded by The University Basic research project of The Educational Department of Liaoning Province (Project JYTMS20231373), Natural Science Foundation of Liaoning Province (Project 2022-KF-12-03), Rolling support of Shenyang Pharmaceutical University Young and Middle-aged Teachers’ Career Development Support Plan (Project ZQN2019019), Shenyang Pharmaceutical University 2024 Innovation and Entrepreneurship Training Program for College Students (Project Song Chenkai) REFERENCESFang Kun Chen YJ.,
“Discussion on quality control management of CT equipment[J],”
China Medical Devices Information, 27
(13), 169
–171
(2021). Google Scholar
McDougald, W.A., Mannheim, J.G., Understanding the importance of quality control and quality assurance in preclinical PET/CT imaging, 9
(1), 77
(2022). Google Scholar
Li N.,
“CT image quality control testing and influencing factors[J],”
China Medical Engineering, 31
(6), 47
–49
(2023). Google Scholar
TAN Xianjian, YAO Guoqing, LI Yiyong, MI Yongwei, LI Tao, ZHENG Xiaoxi, DING Chaofei,
“Analysis of the current status of medical CT quality control[J],”
China Medical Equipment, 14
(2), 106
–109
(2017). Google Scholar
Chatterjee, D., Shen, T.C., Mukherjee, P. et al, Automated detection of incidental abdominal aortic aneurysms on computed tomography, 49
(2), 642
–650
(2024). Google Scholar
QIN Yongchun, XU Xiaosan,
“A preliminary study on linear automatic analysis of CT number in performance testing of X-ray computed tomography devices[J],”
Medical and Health Equipment, 44
(5), 50
–53
(2023). Google Scholar
National Ionizing Radiation Measurement Technical Committee, JJG 961-2017 X-ray radiation source for medical diagnostic spiral computed tomography device (CT) [S], China Metrology Press, Beijing
(2017). Google Scholar
Ministry of Health of the People’s Republic of China, Specification for quality assurance testing of X-ray computed tomography devices: GB 17589-2011 [S], China Standard Press, Beijing
(2012). Google Scholar
Liao X, Su LN, Ma YZ,
“Quality Control Test and Analysis of 105 CT in Beijing [J],”
China Medical Equipment, 39
(5), 37
–4098
(2024). Google Scholar
Wang Shaohua,
“CT machine metrology detection and image quality control[J],”
China Metrology,
(6), 100
–102
(2016). Google Scholar
Huang Haixiao, Bao Lingyan,
“Image data SPSS analysis of the effect of different modalities on CT image quality control[J],”
Equipment Management and Maintenance,
(4), 27
–29
(2023). Google Scholar
YAN Hongbo, WANG Mengwei, XU Yang, ZHAO Pengbo, LIU Yu,
“Design of medical image processing simulation system based on Matlab GUI[J],”
Computer Applications, 40
(S1), 226
–230
(2020). Google Scholar
|