Journal of Robotics and Control (JRC) Volume 6. Issue 4, 2025 ISSN: 2715-5072 DOI: 10. 18196/jrc. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Phichitphon Chotikunnan 1. Rawiphon Chotikunnan 2. Tasawan Puttasakul 3*. Wanida Khotakham 4. Pariwat Imura 5. Jaroonrut Prinyakupt 6. Nuntachai Thongpance 7. Anuchart Srisiriwat 8* 1, 2, 3, 4, 5, 6, 7 College of Biomedical Engineering. Rangsit University. Pathum Thani. Thailand Department of Electrical Engineering. Pathumwan Institute of Technology. Bangkok. Thailand Email: 1 phichitphon. c@rsu. th, 2 rawiphon. c@rsu. th, 3 tasawan. p@rsu. th, 4 wanida. k@rsu. i@rsu. th, 6 jaroonrut. p@rsu. th, 7 nuntachai. t@rsu. th, 8 anuchart@pit. *Corresponding Author AbstractAiThis study offers a systematic approach to improving the reconstruction of three-dimensional anatomical models from CT imaging data. The main difficulty tackled is the maintenance of internal bone features during denoising, essential for producing clinically relevant models. A nonlinear filtering strategy was implemented, utilizing a 3y3 median filter alongside manual refinement to eliminate salt-and-pepper noise while preserving anatomical information. The study presents a reproducible image-processing pipeline that improves structural clarity and enables material-efficient 3D printing while preserving internal bone integrity. A publicly available dataset including 813 anonymized chest CT scans . y512 pixels, 16-bit grayscal. from Zenodo was employed. Preprocessing included grayscale normalization, brightness adjustment, and the application of median filters with kernel sizes from 3y3 to 9y9, followed by artifact removal using FlashPrint software before STL conversion. The 3y3 median filter achieved an excellent balance between noise reduction and anatomical clarity, outperforming mean filtering and larger kernels in maintaining edge detail. Although statistical evaluation was not conducted, visual analysis validated an 18. percent decrease in print time and a 17. 88 percent reduction in filament consumption. The technology exhibited actual efficacy in generating high-quality anatomical models. Future endeavors will incorporate automated segmentation and sophisticated denoising methodologies to enhance applicability in surgical simulation, clinical education, and personalized healthcare KeywordsAiMedian Filtering. Medical Image Processing. STL File Generation. DICOM. INTRODUCTION Advancements in three-dimensional . D) printing have revolutionized healthcare by facilitating the creation of patient-specific organ models for surgery planning, medical teaching, and patient communication . , . Computed tomography (CT) is crucial for obtaining high-resolution anatomical information . CT images often exhibit saltand-pepper noise, a form of random impulse noise marked by intermittent high-contrast pixels. If not treated, such noise can hide intricate anatomical boundaries, including vascular structures and organ interfaces, thereby diminishing the accuracy of 3D reconstruction for therapeutic and educational applications . Conventional denoising methods may unintentionally obscure edges or eliminate minor tissue variations essential for diagnostic accuracy. Recent work on advanced preprocessing methods and convolutional neural networks has explored approaches to improve image clarity while maintaining structural integrity . Numerous studies have demonstrated that effective image filtering is crucial for dependable picture-based applications in agriculture . , medical diagnostics . , facial recognition . , robot-assisted perception and object segmentation in robotics . , and healthcare . Deep learning methodologies, including convolutional neural networks (CNN. and ResNet-based architectures, are being widely utilized and have shown effectiveness in chest disease diagnosis, skin lesion classification, and semantic segmentation within both medical and industrial domains . Nonetheless, preprocessing remains crucial, as noise reduction enhances feature clarity and improves the reliability of downstream tasks, including diagnosis and 3D reconstruction . In surgical planning and patientspecific anatomical modeling, minimizing noise is essential to avoid misinterpretation of critical structures and to improve the 3D printing process . , . This study proposes the implementation of a median filter, a recognized non-linear filtering method that mitigates saltand-pepper noise while preserving critical anatomical Median filtering offers a straightforward and effective method that preserves sufficient anatomical fidelity in contrast to more intricate techniques such as bilateral filtering or anisotropic diffusion, which frequently necessitate parameter adjustment and increased computational resources. This work is unique in minimizing internal structural loss usually associated with high-threshold segmentation and in efficiently removing noise artifacts that may compromise the 3D model. Particularly in osseous tissue, this approach helps to maintain internal density, thereby allowing precise physical The organ or region of interest is delineated using thresholding and an erasing function in anatomical modeling software during post-processing. This study meticulously assesses and compares several median filter kernel sizes . y3, 5y5, 7y7, and 9y. to Journal Web site: http://journal. id/index. php/jrc Journal Email: jrc@umy. Journal of Robotics and Control (JRC) ISSN: 2715-5072 determine the best balance between noise reduction and feature preservation . Median filtering is prioritized in this work due to its simplicity and consistent performance in retaining edge clarity with minimal parameter adjustment, despite the availability of more advanced denoising methods. The resulting 3D model is refined and printable, with potential utility in supporting medical professionals in diagnosis, surgical planning, and patient communication . The main aim was to guarantee compatibility with medical imaging software and filtration tools. The growing accessibility of sophisticated imageprocessing methods and open-source imaging datasets has enabled wider applicability across other fields. This research concentrates on medical imaging and anatomical modeling, although certain publications in . emphasize applications in agriculture and robotics. The extensive applications demonstrate the adaptability of filtering procedures, although the current study emphasizes the clinical significance of generating rapid, accurate, and cost-effective anatomical models from CT data. The main contribution of this study is a reproducible and efficient CT image filtering and modeling framework that reduces salt-and-pepper noise, preserves internal bone structure, and enables the generation of anatomically detailed 3D-printed models that may potentially support future applications in clinical planning, education, and simulation. The structure of this research is as follows. Section II offers the research technique, including segmentation of CT images, thresholding, and filtering techniques. Section i reviews the experimental results and assesses the efficiency of the proposed approach in terms of noise attenuation and 3D reconstruction quality. The study concludes in Section IV, which also provides future research directions. II. RESEARCH METHOD This work introduces a comprehensive image-processing pipeline that is designed to improve the quality of CT scans for the purpose of three-dimensional organ modeling and medical 3D printing. The process depicted in Fig. 1 involves a number of phases, such as dataset acquisition, preprocessing, filtering, segmentation, model refinement, and the fabrication of physical 3D models. A comparative assessment across four median filter kernel sizes . y3, 5y5, 7y7, and 9y. was conducted, focusing on visual evaluation of edge preservation and noise suppression. Median filtering using a 3y3 kernel size is highlighted for reducing salt-andpepper noise while keeping important anatomical details needed for 3D reconstruction. The procedure starts with dataset acquisition from a publicly available Zenodo repository . The dataset consists of 813 anonymized DICOM slices, each with a size of 512y512 pixels and a 16-bit grayscale depth, obtained by helical CT scanning. These images were selected to represent clinically relevant spatial resolution and contrast fidelity. It is important to note that pixel-based evaluation was adopted, as the DICOM files lacked physical metadata such as PixelSpacing and SliceThickness. In the preprocessing step using MATLAB R2024b. Hounsfield Unit (HU) values were made consistent, and the brightness was adjusted to ensure uniform contrast across all Fig. Flowchart of the proposed 3D medical modeling pipeline using median filtering, covering all stages from DICOM acquisition to STL-based 3D printing preparation During the median filtering phase, four kernel sizes, 3y3, 5y5, 7y7, and 9y9, were utilized with MATLAB's medfilt2 No quantitative metrics such as PSNR or SSIM were applied in this stage. however, the filtered outputs were qualitatively compared for noise suppression, anatomical Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 border clarity, and visual print quality. The 3y3 kernel exhibited the most advantageous compromise. The filtered images preserved anatomical border clarity and removed high-frequency noise, thereby offering a reliable foundation for further segmentation and modeling tasks. After filtering, thresholding and segmentation were used to highlight anatomical features such as bones and soft tissues based on specific HU ranges determined through experimentation . or example, >200 HU for bon. Binary masks were created to define pertinent anatomical areas, serving as the foundation for volume depiction. Segmentation and Thresholding The DCM_01. rar archive . , which is accessible at . ttps://zenodo. org/records/15244881/files/DCM_01. rar?do wnload=. , was accessed to obtain a publicly available CT chest dataset for preliminary preprocessing. MATLAB R2024b was employed to accomplish this task. This collection consists of 813 anonymized DICOM segments in a standard 16-bit grayscale format, each with a spatial resolution of 512y512 pixels. The data were acquired using a helical CT scan technique and have been completely de-identified in accordance with international privacy standards. The DICOM series was imported into MATLAB for The Hounsfield Unit (HU) values of anatomical areas of interest, primarily osseous structures, were used to delineate them using intensity-based Non-relevant regions, such as soft tissues and air, were excluded, and voxels exceeding 200 HU were identified as bone tissue. In the absence of metadata fields such as PixelSpacing and SliceThickness, all measurements and filtering evaluations were conducted at the pixel level without spatial calibration. To improve anatomical contrast, the HU thresholds were determined by clinical imaging standards and visual inspection. The binary templates that were generated were employed for the purpose of volumetric reconstruction, segmentation, and additional denoising. No annotated ground-truth segmentation masks were available for benchmarking, so qualitative validation relied on visual consistency with anatomical references derived from standard CT atlases and prior literature. No clinical expert validation or inter-observer assessment was performed, which is acknowledged as a limitation and an area for future improvement. The 3D compilation of the original DICOM segments is illustrated in Fig. 2, which facilitates 3D rendering and multiplanar visualization . xial, sagittal, and corona. These outputs are well-suited for research, instructional purposes, and preliminary prototyping in biomedical engineering Fig. Workflow of CT image processing, starting from stacked DICOM slices . , through multiplanar reformatting . , to final 3D volume rendering . Segmentation and Thresholding The quality of medical imaging is significantly enhanced by noise reduction, particularly in CT scans that are used for 3D reconstruction and anatomical modeling. Salt-and-pepper noise is a notable occurrence, characterized by the arbitrary arrangement of brilliant and dark pixels. This noise has the potential to substantially alter anatomical boundaries and reduce the accuracy of segmentation, potentially affecting the reliability of downstream applications including surgical planning, diagnosis, and 3D printing. Median filtering is the primary preprocessing technique employed in this investigation to address this matter. It is a nonlinear technique that is highly regarded for its ability to reduce impulse noise while preserving critical edge details. This method is especially advantageous in medical applications that necessitate precise organ boundary Unlike linear filters such as the mean filter, which may blur structural features, the median filter replaces each pixel with the median value of its neighborhood, thereby supporting the retention of fine anatomical structures. The primary objective of this investigation is to enhance the precision of CT image data obtained from a modified DICOM dataset by employing median filtering. To mitigate noise artifacts before key operations, the filtering process is performed prior to thresholding and 3D model generation. This section outlines the justification for applying median filtering and presents the framework adopted to assess its The subsequent sections explain . why different median filter window sizes were chosen and how they were selected based on visual fidelity and processing considerations and . the steps taken using MATLAB, including preparation steps and the method used to measure how well noise was reduced, how long processing took, and how well anatomical details were kept. Choice of Median Filter An often-occurring artifact in medical imaging, especially in CT scans, is salt-and-pepper noise. It appears as haphazardly distributed high-contrast pixels that significantly reduce image quality and compromise the accuracy of subsequent image segmentation and 3D reconstruction. Such noise introduces abrupt intensity variations that are challenging to remove using traditional linear filters. The median filter, a nonlinear spatial filter that is effective in reducing impulsive noise while preserving the structural integrity of anatomical details, is employed in this study to resolve this issue. The median filter replaces the intensity of each pixel with the median value of its neighbors, in contrast to linear filters like the Gaussian or average filter, which frequently blur boundaries. By refraining from averaging, which would otherwise obfuscate fine details, this replacement mechanism maintains sharp anatomical This method is especially beneficial for the precise preservation of organ boundaries, which is crucial for clinical applications. The median filtering process is quantitatively expressed by the following . Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 yci. cu, y. = ycoyceyccycnycaycu. cn, y. cn, y. OO yui. cu, y. } . where yce. cn, y. represents the original intensity values inside the neighborhood yui. cu, y. , and yci. cu, y. indicates the output intensity following filtering. This formulation effectively removes outlier values while preserving surrounding anatomical consistency. Four different window sizes . y3, 5y5, 7y7, and 9y. were evaluated in this work. Each size offers a different balance between detail retention and noise-lowering Smaller windows retain more anatomical detail but may be less effective at eliminating extensive noise, whereas larger windows suppress noise more thoroughly but may reduce structural fidelity. The decision to assess several kernel sizes was motivated by the need to find the best tradeoff between fidelity and efficiency. As shown in Fig. 3, the 3y3 median filter ranks the intensity levels of every pixel within its local 3y3 neighborhood window . arked in re. In the example, the center pixel . ighlighted in gree. and its surrounding neighbors are used. The window's sorted pixel values allow one to choose a median value to replace the center pixel. Applying this method across the image results in a filtered output that maintains important structures and lowers impulse The right panel in Fig. 3 shows clear edge preservation together with noise artifact reduction, thereby demonstrating the pragmatic efficiency of median filtering. Fig. Example of 3y3 median filtering Implementation This work used 813 axial CT images in DICOM format, each with a 16-bit grayscale depth and 512y512 pixels. Obtained from the CT chest collection AuDCM_01,Ay which was developed to comply with clinical imaging guidelines and is publicly accessible and anonymized, these images were used for analysis. For preprocessing, each segment was imported in MATLAB R2024b using the functions AudicominfoAy and AudicomreadAy. Bit-depth normalization and grayscale scaling were used to standardize intensity representation across all Using four different square kernel sizes . y3, 5y5, 7y7, and 9y. , the aim was to reduce salt-and-pepper noise across all images. This was accomplished using the built-in function This filtering made comparative evaluation of anatomical detail retention against noise suppression possible. The syntax employed in this procedure is delineated in . yceycnycoycyceycyceyccyaycoycayciyce = ycoyceyccyceycnycoyc2. cuycycnyciycnycuycaycoyaycoycayciyce, . c y. ) . The variable yc denotes the width of the square kernel and belongs to the set . ,5,7,. Every image underwent filtering separately to ensure consistency and preserve localized Archived and then evaluated using FlashPrint 5 tools, the filtered sections served for visual inspection and 3D model For each kernel size, key performance parameters were recorded, including approximate processing time per segment, edge fidelity, and visual noise reduction. These assessments helped to identify the 3y3 kernel arrangement that offers the best desired balance between anatomical integrity and filtering performance. Quantitative metrics were not computed, which is acknowledged as a limitation and will be addressed in future investigations. Mean Filtering for Noise Reduction The mean filter is a fundamental linear filtering technique that is frequently employed in the field of image processing to attenuate noise and normalize images. It operates by replacing each pixel with the average value of its neighboring pixels within a predetermined window size. This technique is effective in reducing Gaussian noise or minor intensity variations by smoothing random fluctuations in pixel The mathematical expression of the mean filter's output at pixel location . cu, y. is as follows . cu, y. = Oc. cn,ycOOyui. cu,y. ) yce. cn, y. The intensity value of the input image at coordinates . cn, y. is denoted by yce. cn, y. , while the neighborhood window centered around the pixel at . cu, y. is defined by yui. cu, y. The variable ycA represents the total number of pixels in the window, such as 9 for a 3y3 kernel. The average of all pixel values within the window is used to calculate the output intensity yci. cu, y. This method smooths local intensity variations but may blur structural boundaries. The computational simplicity and ease of implementation of the mean filter make it advantageous in general image processing tasks. Nevertheless, its primary drawback is its tendency to blur edges and reduce structural clarity. The uniform treatment of all pixels within the window frequently mitigates abrupt transitions, such as edges or complex anatomical structures, which are essential in medical imaging applications such as CT and MRI. The mean filter is less suitable for tasks that require structural accuracy, such as organ segmentation or 3D model generation, due to its tendency to degrade anatomical boundaries, even though it improves image uniformity in medical applications. In such contexts, median filtering is the preferable method because it effectively eliminates impulsive noise, including salt-and-pepper noise, while preserving In conclusion, the mean filter is appropriate for generalpurpose denoising, but its application in high-precision medical imaging requires careful evaluation, especially when the preservation of anatomical detail is critical. Model Refinement, 3D Reconstruction, and Printing Workflow The erase tool in FlashPrint 5 was utilized to manually eradicate any remaining noise artifacts, such as floating Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 voxels, speckle noise, or broken parts, following the median filtering and segmentation procedure. This step helped improve model clarity and reduce non-anatomical fragments, enhancing the suitability of the model for fabrication using FDM printing. The erase function effectively eliminated falsepositive areas that were introduced by thresholding or aggressive filtration while preserving the main anatomical Subsequently, a comprehensive 3D organ model was produced by volumetrically stacking the filtered and corrected image slices. The model was converted into a surface geometry using MATLAB and FlashPrint 5, and the resulting file was exported in STL format for slicing applications. order to prevent printing failures and maintain the structure of the model, the mesh structure of each STL file was visually examined, and obvious artifacts were manually removed. The desktop FDM 3D printer was utilized, and FlashPrint 5 was the slicing software employed in the construction stage. PLA filament was used for printing, and its specifications included a layer height of 0. 1 to 0. 2 mm, an infill density of 15 percent to 25 percent, and a nozzle temperature of 210 degrees Celsius. These printing parameters were selected to achieve stable fabrication and recognizable anatomical form. The computational framework for image processing and reconstruction efficiently processed the 813-slice dataset on a workstation outfitted with an Intel Core i7 CPU, 32 GB of RAM, and 1 TB of SSD storage, utilizing MATLAB R2024b. A consistent pipeline was implemented for converting CT data affected by noise into 3D-printable anatomical representations through filtering, manual refinement, model reconstruction, and physical printing. The STL file generated by MATLAB and prepared in FlashPrint 5 prior to 3D printing is illustrated in Fig. The reconstructed rib and pelvic skeleton model from CT data is depicted in the left panel, while the same model with automatically generated support structures is depicted in the right panel. This example demonstrates the connection between digital processing and fabrication steps in anatomical modeling workflows. Fig. STL file of the rib and pelvic skeleton created in MATLAB and configured in FlashPrint 5 with automatically built supports for 3D printing i. RESULTS AND DISCUSSION The median filter's efficacy in noise reduction and its influence on the reconstruction of 3D models were assessed using a CT dataset featuring 813 DICOM segments . y512 pixels, 16-bit grayscal. The study evaluated the efficacy of four distinct median filter window sizes, 3y3, 5y5, 7y7, and 9y9, in suppressing salt-and-pepper noise while maintaining anatomical fidelity. The evaluation was conducted through visual inspection rather than quantitative image quality metrics due to the lack of ground truth annotations and the unavailability of standardized segmentation references. The primary emphasis was placed on visual clarity and practical feasibility for medical 3D printing rather than statistical image analysis, which is acknowledged as a methodological limitation in this experimental study. Comparative Examination of Filtering Methods for CT Spine Visualization The initial part of the image processing technique involved selecting an area of interest (ROI) from a standard axial CT This step was intended to evaluate the preliminary effectiveness of the pipeline in separating anatomical features before the application of further segmentation and 3D reconstruction techniques. Fig. shows the raw CT slice, highlighting different body parts such as bone, soft tissue, and air spaces, along with a basic grayscale intensity pattern. A basic thresholding method was applied to identify clinically significant regions. The threshold range was manually adjusted through visual inspection until skeletal structures, including ribs and spinal segments, were sufficiently separated. Fig. shows the initial selection, with green highlights indicating the areas identified based on this visual segmentation. The segmentation results were improved through filtering and noise reduction after thresholding, as illustrated in Fig. This resulted in a more precise and focused Soft tissue artifacts and background noise were minimal, while key features and skeletal structures persisted. This refinement step contributed to improved visual clarity and supported STL export for subsequent 3D printing. Fig. CT image preprocessing workflow: . original, . segmented ROI, . final binary mask The simplicity and adaptability of this threshold-based method are its greatest assets. Specific threshold values were not predefined. however, dynamic adjustments were implemented for each scan to improve visibility and anatomical precision. The results indicate that when combined with visual confirmation and domain knowledge, fundamental segmentation techniques can yield dependable results. In general, this segmentation procedure is a critical preprocessing stage that converts raw CT data into refined binary masks. It guarantees that the final model contains only pertinent anatomical structures and reinforces the stability of downstream tasks in MATLAB and 3D reconstruction Comparative Examination of Filtering Methods for CT Spine Visualization This study assessed the effectiveness of several filtering techniques in enhancing CT image quality for 3D reconstruction through visual and comparative evaluations of multiple image preprocessing methods. The focus was on Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 preserving anatomical integrity, particularly the spinal structure, while effectively reducing salt-and-pepper noise. were significantly lost, especially in areas with faint edge Fig. 6 shows sagittal CT scans of the spine that have been processed through several stages. The original image Fig. 6(A) contained a considerable amount of high-contrast noise, but raw anatomical features were retained. A thresholding process was used in Fig. 6(B) to separate skeletal components and eliminate soft tissue. In order to improve the contrast between bone features and background and enable more accurate segmentation, the threshold range was manually adjusted through visual inspection. A 3y3 mean filter was used in Fig. 6(C), which produced a smoother output but blurred small bony structures and reduced edge definition. On the other hand, the 3y3 median filter in (D) significantly reduced high-frequency noise while preserving edge sharpness. As the kernel size was increased to 5y5, 7y7, and 9y9 Fig. 6(E) to (G), over-smoothing caused an increasing loss of anatomical detail. Fig. Close-up view of CT slice filtering results: (A) original, (B) thresholded . , (C) mean 3y3, (D) median 3y3, (E) median 5y5, (F) median 7y7, (G) median 9y9 TABLE I. QUALITATIVE COMPARISON OF FILTERING PERFORMANCE ON CT IMAGES Filter Type Mean 3y3 Median 3y3 Median 5y5 Median 7y7 Median 9y9 Noise Removal Edge Clarity Anatomical Detail Visual Quality Moderate Low Blurred Fair High High Preserved Excellent Very High Moderate Slightly Good Very High Low Flattened Fair Very High Very Low Flattened Poor Due to the absence of annotated ground truth, the evaluation was based on qualitative visual inspection rather than statistical validation. Metrics such as PSNR. SSIM, or edge preservation scores were not applied. The emphasis was on practical outcomes relevant to 3D printing and segmentation workflows, despite the limitations in numerical Fig. Sample 2D CT image processing: (A) original, (B) thresholded . , (C) mean 3y3, (D) median 3y3, (E) median 5y5, (F) median 7y7, (G) median 9y9 These results are in line with earlier studies that show median filters preserve important edge information in medical imaging better than mean filters, especially for bone This work is unique in that it confirms that the 3y3 kernel offers a good trade-off between fidelity and filtering efficiency by applying these observations to a real-world 3D printing scenario. These differences are more noticeable in Fig. 7, which shows a magnified view of the lumbar spine region. While larger window sizes Fig. 7(F) and (G) made it harder to see vertebral boundaries, the 3y3 median filter reduced noise while maintaining structural clarity. The mean filter Fig. 7(C) resulted in transition blurring and disrupted the continuity between adjacent vertebrae. The findings support that the 3y3 median filter is the most suitable option for preprocessing CT data for anatomical It minimizes the need for further correction while supporting diagnostic precision, especially when threshold values are carefully adjusted for each dataset. Larger filters may be acceptable in low-detail applications but are not recommended for tasks involving skeletal segmentation. The comparative outcomes of the filtering techniques are summarized in Table I. The best visual quality was shown by the 3y3 median filter, which successfully reduced noise while preserving anatomical features. Stronger noise reduction was achieved with larger kernel sizes, but fine skeletal features Evaluation and Comparison of Median Filter Effectiveness The performance of median filters applied to CT images derived from DICOM data with varying kernel sizes was investigated in this study. These images enabled the Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 construction of precise three-dimensional medical models, where anatomical fidelity is crucial. window sizes, on the other hand, frequently resulted in oversmoothing and the loss of clinically meaningful detail. Applying each median filter size to the initial CT slices was the first step in the evaluation. To evaluate the degree of noise reduction and anatomical contour preservation, the results were visually examined. As demonstrated in Fig. 8, the 3y3 median filter successfully removed small noise while preserving anatomical structure and edge sharpness. Significant smoothing was introduced by larger filters, like 7y7 and 9y9 . ee Fig. 8(E)), which reduced boundary clarity and eliminated fine anatomical information. To examine filter performance in localized regions. Fig. 9 presents an enlarged pixel cluster containing 64 pixels. In this context, the 3y3 median filter demonstrated good noise reduction while maintaining pixel structure, especially in regions with delicate intensity changes or complex textures. Since there were no annotated ground truth data available for quantitative comparison, this evaluation was done solely One known drawback is that no statistical measures like SSIM or Dice scores were calculated. Instead of objectively benchmarking image quality, the focus was on the usefulness of filtered images for 3D reconstruction and These outcomes are consistent with prior research findings, which reported that median filtering is more effective than linear filters in preserving edge features in medical imaging, especially for osseous structures. While those prior studies emphasized algorithmic evaluation, this work applies such findings to a physical 3D modeling workflow and demonstrates that the 3y3 kernel offers a practical and efficient balance between detail preservation and noise suppression. The results support the conclusion that the 3y3 median filter is the most appropriate choice for CT preprocessing in anatomical modeling. Particularly when threshold values are customized for specific datasets, it minimizes the need for post-processing while maintaining diagnostic integrity. Larger filters are generally discouraged for tasks involving bone segmentation, though they might be useful in situations where high detail is not required. Fig. Median filtering results for a CT region with size 66 y 55 pixels. (A) Original image, (B) Median filter 3y3, (C) Median filter 5y5, (D) Median filter 7y7, (E) Median filter 9y9 Fig. Example of a small pixel cluster consisting of 64 pixels. (A) Original image, (B) Median filter 3y3, (C) Median filter 5y5 A comparison with the mean filter was conducted to highlight the relative advantages of median filtering. shown in Fig. 10, four images are displayed in sequence: the original, the thresholded version, and the outputs from the mean and median filters. Compared to the mean filter, the 3y3 median filter showed superior edge retention and noise More intricate anatomical shapes were the subject of additional comparisons. For example, a region with irregular contours and curvature is shown in Fig. While the mean filter blurred internal boundaries, the 3y3 median filter produced more distinct structural edges. The observed benefit of median filtering is further supported by Fig. 12 and Fig. 13, which show this behavior in smaller segmented regions with various geometric properties. The collective visual evidence from Fig. 8 to Fig. indicates that the 3y3 median filter is the most effective in reducing noise while preserving anatomical accuracy. Larger Fig. Comparison of image filtering results for a region with size 66 y 55 (A) Original image, (B) Thresholded image, (C) Mean filter 3y3, (D) Median filter 3y3 Fig. Comparison of image filtering results for a complex region with size 40 y 48 pixels. (A) Original image, (B) Thresholded image . , (C) Mean filter 3y3, (D) Median filter 3y3 Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 Fig. Example of segmented region with size 12 y 30 pixels. (A) Original image, (B) Thresholded image . , (C) Mean filter 3y3, (D) Median filter 3y3 Fig. Example of segmented region with size 21 y 28 pixels. (A) Original image, (B) Thresholded image . , (C) Mean filter 3y3, (D) Median filter 3y3. Assessment of 3D Model Precision Post Noise Mitigation This section assesses the influence of median filtering on the precision of 3D reconstructions obtained from DICOMformatted CT scan data. Four preprocessing strategies were analyzed, including original pictures, thresholding without filtering, a 3y3 mean filter, and a 3y3 median filter. The reconstructed models were examined in three dimensions to assess structural clarity and residual noise artifacts. Fig. 3D Bone Model Comparison Using Different Filtering Techniques: (A) Original Image, (B) Thresholding . , (C) Mean Filter 3y3, (D) Median Filter 3y3 Fig. 14 demonstrates that thresholding alone produced considerable surface abnormalities, notably fine voxel clusters linked to salt-and-pepper noise. The mean filter, although generating smoother surfaces, induced blurring effects that obscured fine anatomical contours such as spinal outlines and rib curvatures, thereby reducing visibility. In contrast, the application of the 3y3 median filter preserved skeletal boundaries while effectively minimizing random noise, resulting in structurally coherent models. These conclusions are further substantiated by the broadened viewpoint in Fig. The raw and thresholded datasets clearly displayed irregularities and discontinuities, whereas the median-filtered reconstruction exhibited smoother and more continuous anatomical outlines. The mean filter compromised the visual clarity of critical skeletal features, which may affect diagnostic interpretability. The results collectively offer persuasive qualitative evidence that the 3y3 median filter effectively balances noise reduction with the preservation of anatomical detail. This observation aligns with earlier research highlighting the strength of median filtering in edge preservation, especially in bone-focused imaging contexts. Although no quantitative metrics were calculated, the visual consistency and structural integrity of the output validate the practical use of the 3y3 median filter in clinical modeling and medical prototyping. Fig. Zoomed-in Rib Cage Region Comparison: (A) Original Image, (B) Thresholding . , (C) Mean Filter 3y3, (D) Median Filter 3y3 Assessment of 3D Model Readiness through Slicing Simulation and Material Estimation This investigation assessed the preparedness of CTderived anatomical models for 3D printing with FlashPrint The workflow encompassed slicing simulation, support generation, and material consumption estimation. Three preprocessing procedures were evaluated, including thresholding without filtering, the use of a 3y3 average filter followed by thresholding, and the use of a 3y3 median filter followed by thresholding. The results illustrated in Fig. 16 through Fig. 18 and Table II indicate that thresholding alone produced the maximum number of support structures and the largest filament usage. Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 amounting to 103. 67 meters. The total print duration was almost 61 hours and 3 minutes. The pronounced voxel-level noise seen in the unfiltered model led to suboptimal material use and extended fabrication time. The application of an average filter prior to thresholding significantly diminished noise artifacts, which consequently enhanced the surface smoothness of the resultant model. This method, however, also resulted in edge blurring, reducing anatomical clarity. Filament usage was minimized to 90. meters, and the printing time was shortened to 40 hours and 25 minutes. The utilization of a 3y3 median filter yielded the most structurally distinct and print-optimized model. This arrangement maintained essential anatomical contours while reducing superfluous supports. Consequently, filament consumption was diminished to 85. 13 meters, with the projected printing duration being 50 hours and 1 minute. shown in Table II, the model processed with a 3y3 median filter necessitated the least filament and produced the most balanced print duration. This configuration resulted in an 07 percent reduction in time and a 17. 88 percent drop in material usage compared to the unfiltered baseline. These findings underscore the practical advantage of incorporating median filtering into the preprocessing pipeline for medical 3D printing. The results confirm that a 3y3 median filter improves printing efficiency while better maintaining anatomical characteristics compared to average filtering or unfiltered This method corresponds with the overarching goal of enhancing computational models for subsequent physical reconstruction in biomedical applications. Fig. 3D model Preview from thresholded CT data without filtering Fig. 3D model preview with average filter 3y3 and thresholding Fig. 3D model preview with median filter 3y3 and thresholding TABLE II. COMPARISON OF 3D PRINTING TIME AND FILAMENT CONSUMPTION BASED ON PREPROCESSING METHODS Method Print Time Filament Used . Time Saved . Time Saved (%) Filament Saved . Filament Saved (%) 61h 3min Avg Filter 3y3 Thresholding 40h 25min Median Filter 3y3 Thresholding 50h 1min Thresholding Only IV. CONCLUSION This study presents an effective technique to improve the accuracy of 3D organ models obtained from CT medical images by utilizing a 3y3 median filter alongside further manual modification. The suggested method efficiently diminishes salt-and-pepper noise while maintaining anatomical features crucial for precise 3D reconstructions. The 3y3 median filter attained the most equitable compromise between noise reduction and structural integrity, enhancing the surface quality of models while maximizing material efficiency for 3D printing. These advantages indicate the method's potential utility in medical environments, such as anatomical modeling for clinical teaching. However, the current study does not include clinical validation or crossdataset performance, so its broader applicability remains Nonetheless, additional confirmation is required prior to incorporation into diagnostic or surgical planning The limitations of this study include the lack of quantitative image quality assessments, such as the Dice coefficient or SSIM, and reliance on a single dataset that lacks anatomical variety. The segmentation pipeline depends on manually adjusted thresholds and deletion steps, which may affect repeatability and scalability. Furthermore, no physical measurement or expert-based validation of printed outputs was conducted to confirm anatomical accuracy. Future research may address these constraints by incorporating automated. AI-driven segmentation methods, such as U-Net or GANs, to improve accuracy and applicability. Extending experiments to heterogeneous datasets and incorporating reproducibility controls is also essential. Investigating comparative benchmarking alongside advanced denoising methods and the physical validation of printed models is crucial to enhance therapeutic significance. Phichitphon Chotikunnan. Noise-Reduced 3D Organ Modeling from CT Images Using Median Filtering for Anatomical Preservation in Medical 3D Printing Journal of Robotics and Control (JRC) ISSN: 2715-5072 ACKNOWLEDGMENT The researcher would like to thank the Research Institute. Academic Services Center, and College of Biomedical Engineering. Rangsit University for the grant of research funding to the research team. The authors would like to convey their profound appreciation to Assoc. Prof. Dr. Manas Sangworasil and Prof. Dr. Takenobu Matsuura for their foundational advice and oversight during the first phases of this research effort. Their mentorship offered essential guidance for the effective completion of the project. Furthermore, it is confirmed that the project has been reviewed by the Ethics Review Board of Rangsit University, with reference number RSUERB2025-007, which certifies that the research does not involve human subjects. Moreover. AI-driven methods (QuillBot Premiu. were utilized for grammatical verification, paraphrasing, and linguistic augmentation to ensure the accuracy and clarity of the text. All CT DICOM images included in this investigation were anonymized prior to analysis. The anonymization process involved removing any personal health information (PHI) and anything related to the patient's identity, hospital, study details, device serial numbers, and private DICOM tags. The experiment was conducted in accordance with established rules for ethical data exchange and the protection of human REFERENCES