Systems and methods for model-based optimization of spinal cord stimulation electrodes and devices are disclosed. According to an aspect a method includes providing a patient-specific electroanatomical model including the spine, spinal cord, and a map of target neural elements and non-target neural elements. The method also includes using model electrodes to stimulate the target neural elements. Further, the method includes determining differences in activation thresholds between the target neural elements and the non-target neural elements in a plurality of different configurations of the model electrodes. The method also includes generating an optimal spinal cord stimulation electrode configuration based on the determined differences in activation thresholds.