Математика за MKMetersPerMapPointAtLatitude
Я пытаюсь преобразовать код отображения Apple в Java. Я почти полностью конвертировал его, за исключением нескольких вызовов MKMetersPerMapPointAtLatitude
У меня очень близкое решение... но это не точно, и я не уверен, почему нет. Любые идеи?
#import <Foundation/Foundation.h>
#import <Math.h>
@import MapKit;
#define MERCATOR_OFFSET 268435456.0 / 2.0
#define MERCATOR_RADIUS (MERCATOR_OFFSET/M_PI)
#define WGS84_RADIUS 6378137.0
#define POINTS_PER_METER (MERCATOR_RADIUS / WGS84_RADIUS)
double MyMetersPerMapPointAtLatitude(double latitude) {
return 1.0 / (POINTS_PER_METER / cos(latitude * M_PI / 180.0));
}
int main(int argc, const char * argv[]) {
@autoreleasepool {
double latitude = 33.861315;
for (int i = 0; i < 100; i++) {
double a = MKMetersPerMapPointAtLatitude(latitude);
double b = MyMetersPerMapPointAtLatitude(latitude);
NSLog(@"%f %f", a, b);
latitude += .1;
}
}
return 0;
}
Распечатывает следующие
2015-05-19 09:13:00.334 Test[92619:5369062] 0.123522 0.123969
2015-05-19 09:13:00.335 Test[92619:5369062] 0.123379 0.123824
2015-05-19 09:13:00.335 Test[92619:5369062] 0.123236 0.123678
2015-05-19 09:13:00.335 Test[92619:5369062] 0.123092 0.123532
2015-05-19 09:13:00.335 Test[92619:5369062] 0.122948 0.123386
2015-05-19 09:13:00.335 Test[92619:5369062] 0.122804 0.123239
2015-05-19 09:13:00.335 Test[92619:5369062] 0.122659 0.123092
...etc
Ответы
Ответ 1
Для начала мы можем изменить вашу функцию, чтобы сделать ее более читаемой:
double MyMetersPerMapPointAtLatitude(double latitude) {
return cos(latitude * M_PI / 180.0) / POINTS_PER_METER;
}
Теперь, проблема, как указывает Томми, заключается в том, что вы не учитываете сглаживание Земли. Вы можете сделать это с помощью
double f = 1/298.257223563; // WGS84 flattening
double MyMetersPerMapPointAtLatitude(double latitude) {
return (1-f) * cos(latitude * M_PI / 180.0) / POINTS_PER_METER;
}
Это приведет к ошибке до той точки, где я бы привязал ее к округлению и усечению.