Построение смещений на карте

[TryingToSolve]

Пытаюсь выполнить задачу построения карты используя подложку tile сервера OpenStreetMap на которую нанесены: опорные точки, точки профилей 1 и 2 соединенные линиями различных цветов. Для каждой из точек профиля должны быть указаны направление (az в градусах) и величина смещения (value) в виде стрелки. Исходные считываемые данные приложены в .zip.

При выполнении кода не отображаются направление (az в градусах) и величинf смещения (value) в виде стрелок. Изменения значений scale_factor не помогают. Подскажите, пожалуйста, что же не так? Никак не могу найти ошибку. Заранее спасибо за ответ!

import numpy as np
import matplotlib.pyplot as plt
from cartopy import crs as ccrs
import shapefile
from cartopy.io.img_tiles import OSM
from matplotlib.patches import FancyArrowPatch

shpfile = shapefile.Reader("/User/Desktop/tdata/tdata")

fields = [x[0] for x in shpfile.fields[1:]]

points = {'base': [], 'profile1': [], 'profile2': []}
for item in shpfile.iterShapeRecords():
kind = item.record[fields.index('kind')]
position = item.shape.points[0]
points[kind].append((position, dict(zip(fields, item.record))))

tiles_provider = OSM()

plt.figure(figsize=(16, 16))

ax = plt.axes(projection=tiles_provider.crs)

ax.set_extent([shpfile.bbox[0]-0.1, shpfile.bbox[2]+0.1,
shpfile.bbox[1]-0.1, shpfile.bbox[3]+0.1], ccrs.PlateCarree())

ax.add_image(tiles_provider, 11)

def draw_arrow(ax, point, angle, length, scale_factor=1, **kwargs):
angle_rad = np.deg2rad(angle)
dx = scale_factor * length * np.cos(angle_rad)
dy = scale_factor * length * np.sin(angle_rad)
point_proj = ccrs.PlateCarree().transform_point(point[0], point[1], ax.projection)
arrow = FancyArrowPatch(point_proj, (point_proj[0] + dx, point_proj[1] + dy), connectionstyle="arc3,rad=.5", **kwargs)
ax.add_patch(arrow)

for point, data in points['base']:
ax.plot(point[0], point[1], marker='o', markersize=5, color='red', transform=ccrs.PlateCarree())

for i, (point, data) in enumerate(points['profile1']):
if i > 0:
prev_point = points['profile1'][i-1][0]
ax.plot([prev_point[0], point[0]], [prev_point[1], point[1]], color='blue', transform=ccrs.PlateCarree())
draw_arrow(ax, point, data['az'], data['value'], scale_factor=10, color='blue', transform=ax.projection)

for i, (point, data) in enumerate(points['profile2']):
if i > 0:
prev_point = points['profile2'][i-1][0]
ax.plot([prev_point[0], point[0]], [prev_point[1], point[1]], color='green', transform=ccrs.PlateCarree())
draw_arrow(ax, point, data['az'], data['value'], scale_factor=10, color='green', transform=ax.projection)

plt.show()

[BDA]

arrow = FancyArrowPatch(point, (point[0] + dx, point[1] + dy), connectionstyle="arc3,rad=.5", **kwargs)
...
draw_arrow(ax, point, data['az'], data['value'], scale_factor=0.002, color='blue', transform=ccrs.PlateCarree())

Вот так рисуются какие-то "загогулины"

[TryingToSolve]

BDA, да, но почему-то значения направления в градусах все-равно не отображаются, хотя все прописано

[BDA]

Почему не отображаются? Чуть-чуть арки наклонены, но в файле как раз значения углов не превышают 12 градусов:

blue 0
blue 10
blue -5
blue 12
blue 7
blue -3
blue 0
blue -10
blue -2
blue 4
green -7
green 3
green 0
green -1
green -4
green -6