In this work, a method to compute the radio interference (RI) lateral profiles generated by corona discharge in high voltage direct current (HVDC) transmission lines is presented. The method is based on a transmission line model that considers the skin effect, through the concept of complex penetration depth, in the conductors and in the ground plane. The attenuation constants are determined from the line parameters and the bipolar system is decoupled by using modal decomposition theory. As application cases, +/- 500 and +/- 600 kV bipolar transmission lines were analyzed. Afterwards, parametric sweeps of five variables that affect the RI levels are presented. Both the RI and the maximum electric field were calculated as a function of sub-conductor radius, bundle spacing, and the number of sub-conductors in the bundle. Additionally, the RI levels were also calculated as a function of the soil resistivity, and the RIV (radio interference voltage) frequency. Following this, vector optimization was applied to minimize the RI levels produced by the HVDC lines and differences between the designs with nominal and optimal values are discussed.