Background: The MYCN copy number category is closely related to the prognosis of neuroblastoma (NB). Therefore, this study aimed to assess the predictive ability of Methods: A retrospective analysis was performed on 104 pediatric patients with NB that had been confirmed by pathology. To develop the Bio-omics model (B-model), which incorporated clinical and biological aspects, PET/CT radiographic features, PET quantitative parameters, and significant features with multivariable stepwise logistic regression were preserved. Important radiomics features were identified through least absolute shrinkage and selection operator (LASSO) and univariable analysis. On the basis of radiomics features obtained from PET and CT scans, the radiomics model (R-model) was developed. The significant bio-omics and radiomics features were combined to establish a Multi-omics model (M-model). The above 3 models were established to differentiate MYCN wild from MYCN gain and MYCN amplification (MNA). The calibration curve and receiver operating characteristic (ROC) curve analyses were performed to verify the prediction performance. Post hoc analysis was conducted to compare whether the constructed M-model can distinguish MYCN gain from MNA. Results: The M-model showed excellent predictive performance in differentiating MYCN wild from MYCN gain and MNA, which was better than that of the B-model and R-model [area under the curve (AUC) 0.83, 95% confidence interval (CI): 0.74-0.92 Conclusions: The M-model model based on bio-omics and radiomics features is an effective tool to distinguish MYCN copy number category in pediatric patients with NB.