Objective. To investigate different dosimetric aspects of 90Y-IsoPet™ intratumoral therapy in canine soft tissue sarcomas, model the spatial spread of the gel post-injection, evaluate absorbed dose to clinical target volumes, and assess dose distributions and treatment efficacy. Approach. Six canine cases treated with 90Y-IsoPet™ for soft tissue sarcoma at the Veterinary Health Center, University of Missouri are analyzed in this retrospective study. The dogs received intratumoral IsoPet™ injections, following a grid pattern to achieve a near-uniform dose distribution in the clinical target volume. Two dosimetry methods were performed retrospectively using the Monte Carlo toolkit OpenTOPAS: imaging-based dosimetry obtained from post-injection PET/CT scans, and stylized phantom-based dosimetry modeled from the planned injection points to the gross tumor volume. For the latter, a Gaussian parameter with variable sigma was introduced to reflect the spatial spread of IsoPet™. The two methods were compared using dose-volume histograms (DVHs) and dose homogeneity, allowing an approximation of the closest sigma for the spatial spread of the gel post-injection. In addition, we compared Monte Carlo-based dosimetry with voxel S-value (VSV)-based dosimetry to investigate the dosimetric differences. Main results. Imaging-based dosimetry showed differences between Monte Carlo and VSV calculations in tumor high-density areas with higher self-absorption. Stylized phantom-based dosimetry indicated a more homogeneous target dose with increasing sigma. The sigma approximation of the 90Y-IsoPet™ post-injection gel spread resulted in a median sigma of approximately 0.44 mm across all cases to reproduce the dose heterogeneity observed in Monte Carlo calculations. Significance. The results indicate that dose modeling based on planned injection points can serve as a first-order approximation for the delivered dose in 90Y-IsoPet™ therapy for canine soft tissue sarcomas. The dosimetry evaluation highlights the non-uniformity of absorbed doses despite the gel spread, emphasizing the importance of considering tumor dose heterogeneity in treatment evaluation. Our findings suggest that using Monte Carlo for dose calculation seems more suitable for this type of tumor where high-density areas might play an important role in dosimetry.