#include #include #include #include my_bool maptile_for_point_init(UDF_INIT *initid, UDF_ARGS *args, char *message) { if ( args->arg_count != 3 || args->arg_type[0] != INT_RESULT || args->arg_type[1] != INT_RESULT || args->arg_type[2] != INT_RESULT ) { strcpy( message, "Your maptile_for_point arguments are bogus!" ); return 1; } return 0; } void maptile_for_point_deinit(UDF_INIT *initid) { return; } long long maptile_for_point(UDF_INIT *initid, UDF_ARGS *args, char *is_null, char *error) { double lat = *(long long *)args->args[0] / 10000000.0; double lon = *(long long *)args->args[1] / 10000000.0; long long zoom = *(long long *)args->args[2]; double scale = pow(2, zoom); double r_per_d = M_PI / 180; unsigned int x; unsigned int y; x = floor((lon + 180.0) * scale / 360.0); y = floor((1 - log(tan(lat * r_per_d) + 1.0 / cos(lat * r_per_d)) / M_PI) * scale / 2.0); return (x << zoom) | y; }