帮我写一个python函数，任意输入一个经纬坐标，输出距离这个经纬坐标在200km到800km范围距离内的所有经纬坐标，要求输出的坐标都是0.25的整数倍数

以下是一个实现该功能的Python函数：

from math import radians, cos, sin, asin, sqrt

def get_nearby_coordinates(lat, lon):
    R = 6371  # 地球半径，单位为千米
    d = 200  # 距离范围下限，单位为千米
    max_d = 800  # 距离范围上限，单位为千米
    step = 0.25  # 坐标间隔，单位为度

    lat, lon = radians(lat), radians(lon)

    result = []

    for i in range(int(d / step), int(max_d / step) + 1):
        d_lat = i * step / R
        for j in range(0, 360, int(step / 0.25)):
            d_lon = i * step / (R * cos(lat))
            lat_new = asin(sin(lat) * cos(d_lat) + cos(lat) * sin(d_lat) * cos(radians(j)))
            lon_new = lon + asin(sin(radians(j)) * sin(d_lon) / cos(lat_new))
            if abs(lat_new) < radians(90) and abs(lon_new) < radians(180):
                result.append((round(lat_new / radians(step)) * step, round(lon_new / radians(step)) * step))

    return result

该函数的参数为经度和纬度，单位为度。函数返回一个列表，包含距离输入坐标在200km到800km范围内的所有坐标，每个坐标为一个元组，包含经度和纬度，单位为度，坐标值为0.25的整数倍数。