With the increasing adoption of building integrated photovoltaics (BIPV), concerns arise about potential glare. While recommended criteria to assess glare exist, it is challenging to apply these in the spatial and temporal domains and communicate the complex data to planning authorities and clients. This paper presents a new computational workflow using annual daylight simulation, material modelling using bi-directional scattering distribution functions (BSDFs) and image-based postprocessing to obtain 3-dimensional renderings of cumulative annual irradiance and glare duration on the built environment. The annual daylight simulation considers relevant sun positions in high temporal resolution (15-minute timesteps) and measured BSDFs to model different PV materials. The postprocessing includes a relative irradiance visualisation comparing the impact of a proposed PV proportional to a reference material. It also includes a new spatio-temporal workflow to assess the glare duration based on recommended thresholds. This workflow is demonstrated with a case study of a proposed PV roof for a church, assessing the glare potential of two different PV materials. The visualisations indicate glare durations well below the thresholds with satinated PVs, and in noncritical zones outside observer positions with standard PVs. Thus the proposed PV roof does not cause any disturbing glare.