Single crystals of Magnesium have been deformed in plane strain compression from room temperature to 450 °C. The deformed crystals have been analysed by EBSD to identify the operative twin variants and parent rotations for comparison with the expected systems. In parallel with these experiments, we have developed a crystal plasticity model based on Taylor principles and the Schmid law, to compute the activated systems together with the twin reorientations and the lattice rotations for large strains. The CRSS values are determined by an iterative procedure that best correlates the experimental data with the numerical simulations for each orientation. From the measured stress/strain data the CRSS for twinning and some slip systems are given over this wide temperature range.