Nanostructures, such a quantum dots or nanoparticles, made of three-dimensional topological insulators (3DTIs) have been recently attracting increasing interest, especially for their optical properties. In this paper we calculate the energy spectrum, the surface states and the dipole matrix elements for optical transitions with in-plane polarization of 3DTI nanocylinders of finite height L L and radius R R . We first derive an effective 2D Hamiltonian by exploiting the cylindrical symmetry of the problem. We develop two approaches: the first one is an exact numerical tight-binding model obtained by discretising the Hamiltonian. The second one, which allows us to obtain analytical results, is an approximated model based on a large- R R expansion and on an effective boundary condition to account for the finite height of the nanocylinder. We find that the agreement between the two models, as far as eigenenergies and eigenfunctions are concerned, is excellent for the lowest absolute value of the longitudinal component of the angular momentum. Finally, we derive analytical expressions for the dipole matrix elements by first considering the lateral surface alone and the bases alone, and then for the whole nanocylinder. In particular, we focus on the two limiting cases of tall and squat nanocylinders. The latter case is compared with the numerical results finding a good agreement.