convert Rijksdriehoek coordinates into Latitude and longitude

example_convert_rd_to_gps.py

from pyproj import Proj, transform

def convert_rd_to_gps(eastings, northings):
    # Define the Rijksdriehoek coordinate system
    rd_proj = Proj(proj="sterea", ellps="bessel", lat_0=52.15616055555555, lon_0=5.3876388888889, x_0=155000, y_0=463000, k=1.0008035)

    # Define the WGS84 coordinate system (GPS)
    wgs84_proj = Proj(proj="latlong", datum="WGS84")

    # Convert RD coordinates to GPS coordinates
    longitude, latitude = transform(rd_proj, wgs84_proj, eastings, northings)

    return latitude, longitude

# Rijksdriehoek coordinates
eastings = 111750
northings = 557250

# Call the function and get GPS coordinates
latitude, longitude = convert_rd_to_gps(eastings, northings)

# Print the individual coordinates
print(f"Latitude: {latitude}")  # Print the latitude
print(f"Longitude: {longitude}")  # Print the longitude