A computer-implemented inverse rendering method is provided. The method comprises: computing an adjoint image by differentiating an objective function that evaluates the quality of a rendered image, image elements of the adjoint image encoding the sensitivity of spatially corresponding image elements of the rendered image with respect to the objective function; sampling the adjoint image at a respective sample position to determine a respective adjoint radiance value associated with the respective sample position; emitting the respective adjoint radiance value into a scene model characterised by scene parameters; determining an interaction location of a respective incident adjoint radiance value with a surface and/or volume of the scene model; determining a respective incident radiance value or an approximation thereof at the interaction location, the respective incident radiance value expressing the amount of incident radiance that travels from the surrounding environment towards the interaction location; and updating a scene parameter gradient by using at least the interaction location and the respective incident adjoint radiance value and the respective incident radiance value or its approximation at the interaction location.