import matplotlib.pyplot as plt
from geopy.distance import geodesic
import numpy as np
from pyproj import Transformer

from shapely.geometry import Point
from shapely.ops import transform as shapely_transform



GEO_TO_CART = Transformer.from_crs("EPSG:4326", "EPSG:25832", always_xy=True)

CART_TO_GEO = Transformer.from_crs("EPSG:25832", "EPSG:4326", always_xy=True)


def geo_to_cart_shape(polygon):
    new_poly = shapely_transform(GEO_TO_CART.transform, polygon)
    return new_poly


def cart_to_geo_shape(polygon):
    new_poly = shapely_transform(CART_TO_GEO.transform, polygon)
    return new_poly




center = (23.252628, 107.531652)  # (lat, lon)
radius_m = 28986.88  
cart_circle = geo_to_cart_shape(Point((center[1],center[0]))).buffer(radius_m)
geo_circle = cart_to_geo_shape(cart_circle)


num_points = 360
circle_lats = []
correct_circle_lons = []

for angle in np.linspace(0, 360, num_points):
    # geodesic use (lat, lon)
    destination = geodesic(meters=radius_m).destination(center, angle)
    circle_lats.append(destination.latitude)
    correct_circle_lons.append(destination.longitude)

points = list(geo_circle.boundary.coords)

buffer_circle_lons = [pt[0] for pt in points]
buffer_circle_lats = [pt[1] for pt in points]

plt.figure(figsize=(10, 10))
plt.plot(correct_circle_lons, circle_lats, "b-", linewidth=2)
plt.plot(buffer_circle_lons, buffer_circle_lats, "r-o", markersize=4, linewidth=2)
plt.scatter(center[1], center[0], c="blue", s=100, marker="*")
plt.xlabel("long")
plt.ylabel("lat")
plt.legend()
plt.grid(True)
plt.axis("equal")
plt.show()

cart_circle = geo_to_cart_shape(Point((center[1],center[0]))).buffer(radius_m)
geo_circle = cart_to_geo_shape(cart_circle)

for angle in np.linspace(0, 360, num_points):
    # geodesic use (lat, lon)
    destination = geodesic(meters=radius_m).destination(center, angle)
    circle_lats.append(destination.latitude)
    correct_circle_lons.append(destination.longitude)

cart_circle = geo_to_cart_shape(Point((center[1],center[0]))).buffer(radius_m)
geo_circle = cart_to_geo_shape(cart_circle)