We compute the molecular line emission of massive protostellar disks by solving the equation of radiative transfer through the cores and disks produced by the recent radiation-hydrodynamic simulations of Krumholz, Klein, & McKee. We find that in several representative lines the disks show brightness temperatures of hundreds of Kelvin over velocity channels ~10 km s^-1 wide, extending over regions hundreds of AU in size. We process the computed intensities to model the performance of next-generation radio and submillimeter telescopes. Our calculations show that observations using facilities such as the EVLA and ALMA should be able to detect massive protostellar disks and measure their rotation curves, at least in the nearest massive star-forming regions. They should also detect significant sub-structure and non-axisymmetry in the disks, and in some cases may be able to detect star-disk velocity offsets of a few km s^-1, both of which are the result of strong gravitational instability in massive disks. We use our simulations to explore the strengths and weaknesses of different observational techniques, and we also discuss how observations of massive protostellar disks may be used to distinguish between alternative models of massive star formation.