In this study, tubular pinewood (Pinus Sylvestris L.) specimens are tested and shear strain measurements are performed by applying torsion in z-direction. Strain gauge measurements are performed for the maximum shear stresses which develop on the tubular specimen, along the radial r (rin, rout), circumferential φ (φin, φout) and z directions, in a point-wise manner. The data is gathered and examined for the determination of the local variations of empirical shear modulus functions on transversely isotropic surfaces of the specimens. The coordinate dependent shear modulus functions of Gzφ(r), Gzφ(φ), Gzφ(z) are derived for Gzφ(r, φ, z) as the function of r, φ and z, respectively, by analyzing the gathered data. It is proposed to represent the shear modulus functions, Gzφ(φ) and Gzφ(z) with the parabolic polynomials, and, to represent the shear modulus function Gzφ(r) with a linear equation.