---
abstract: 'We present a calculation of the electromagnetic form factor of the pion in $N_f=2+1$ flavor lattice QCD. Calculations are made on the PACS-CS gauge field configurations generated using Iwasaki gauge action and Wilson-clover quark action on a $32^3\times64$ lattice volume with the lattice spacing estimated as $a=0.0907(13)$ fm at the physical point. Measurements of the form factor are made using the technique of partially twisted boundary condition to reach small momentum transfer as well as periodic boundary condition with integer momenta. Additional improvements including random wall source techniques and a judicious choice of momenta carried by the incoming and outgoing quarks are employed for error reduction. Analyzing the form factor data for the pion mass at $M_\pi \approx 411$ MeV and 296 MeV, we find that the NNLO SU(2) chiral perturbation theory fit yields $\left< r^2\right>=0.441 \pm 0.046~{\rm fm}^2$ for the pion charge radius at the physical pion mass. Albeit the error is quite large, this is consistent with the experimental value of $0.452\pm 0.011~{\rm fm}^2$. Below $M_\pi\approx 300$ MeV, we find that statistical fluctuations in the pion two- and three-point functions become too large to extract statistically meaningful averages on a $32^3$ spatial volume. We carry out 