A method of performing dose calculations for ion radiotherapy compensating for tissue in which the density within a voxel may be inhomogeneous by approximating a portion of the voxel as an air cavity. For each dose voxel, the voxel is inscribed in a three- dimensional grid comprising a number of cells and the propagation of ions through the voxel is calculated based on the cell pattern in the at least one cell overlapping the voxel. Preferably, the voxel is inscribed in the three-dimensional grid in such a way that it overlaps at least one cell fully. Each cell comprises a first portion representing a first density corresponding to a density of a tissue and a second portion representing a second density corresponding to the density of air, the first and second portions forming a cell pattern.L'invention porte sur un procédé de réalisation de calculs de dose destiné à la compensation de radiothérapie par ions pour un tissu dans lequel la densité à l'intérieur d'un voxel peut être inhomogène par approximation d'une partie du voxel en tant que cavité d'air. Pour chaque voxel de dose, le voxel est inscrit dans une grille tridimensionnelle comprenant un nombre de cellules et la propagation d'ions à travers le voxel est calculée sur la base du motif de cellules dans la ou les cellules recouvrant le voxel. De préférence, le voxel est inscrit dans la grille tridimensionnelle de telle sorte qu'il recouvre entièrement au moins une cellule. Chaque cellule comprend une première partie représentant une première densité correspondant à une densité d'un tissu et une seconde partie représentant une seconde densité correspondant à la densité d'air, les première et seconde parties formant un motif de cellule.