То, что вы хотите сделать, на самом деле должно работать.Имейте в виду, однако, что если север находится сверху, горизонтальная ось (x) - это ДЛИНА, а вертикальная ось (y) - это РЕШЕНИЕ (я думаю, вы могли бы запутать это).
Если вы параметризовалиline as lat = func (long), у вас возникнут проблемы с вертикальными линиями (то есть теми, которые идут точно с севера на юг), поскольку широта меняется, а долгота фиксирована.
Поэтому я бы лучше использовал другую параметризацию:
long(alpha) = long_1 + alpha * (long_2 - long_1)
lat(alpha) = lat_1 + alpha * (lat_2 - lat_1)
и варьируются alpha
от 0 до 1.
Это не совсем точно совпадает с большим кругом (кратчайший путь на сфере), но чем меньшерегион, на который вы смотрите, тем менее заметна будет разница (как отмечали другие авторы).