Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 1 | /* |
| 2 | * Copyright (C) 2012 The Android Open Source Project |
| 3 | * |
| 4 | * Licensed under the Apache License, Version 2.0 (the "License"); |
| 5 | * you may not use this file except in compliance with the License. |
| 6 | * You may obtain a copy of the License at |
| 7 | * |
| 8 | * http://www.apache.org/licenses/LICENSE-2.0 |
| 9 | * |
| 10 | * Unless required by applicable law or agreed to in writing, software |
| 11 | * distributed under the License is distributed on an "AS IS" BASIS, |
| 12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. |
| 13 | * See the License for the specific language governing permissions and |
| 14 | * limitations under the License. |
| 15 | */ |
| 16 | |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 17 | #ifndef LATINIME_GEOMETRY_UTILS_H |
| 18 | #define LATINIME_GEOMETRY_UTILS_H |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 19 | |
| 20 | #include <cmath> |
| 21 | |
| 22 | #define MAX_DISTANCE 10000000 |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 23 | #define MAX_PATHS 2 |
| 24 | |
| 25 | #define DEBUG_DECODER false |
| 26 | |
Ken Wakasa | bcec82d | 2012-08-12 11:10:48 +0900 | [diff] [blame] | 27 | #define M_PI_F 3.14159265f |
| 28 | |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 29 | namespace latinime { |
| 30 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 31 | static inline float squareFloat(float x) { |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 32 | return x * x; |
| 33 | } |
| 34 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 35 | static inline float getNormalizedSquaredDistanceFloat(float x1, float y1, float x2, float y2, |
| 36 | float scale) { |
| 37 | return squareFloat((x1 - x2) / scale) + squareFloat((y1 - y2) / scale); |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 38 | } |
| 39 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 40 | static inline float getSquaredDistanceFloat(float x1, float y1, float x2, float y2) { |
| 41 | return squareFloat(x1 - x2) + squareFloat(y1 - y2); |
| 42 | } |
| 43 | |
| 44 | static inline float getDistanceFloat(float x1, float y1, float x2, float y2) { |
| 45 | return hypotf(x1 - x2, y1 - y2); |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 46 | } |
| 47 | |
Ken Wakasa | bcec82d | 2012-08-12 11:10:48 +0900 | [diff] [blame] | 48 | static inline int getDistanceInt(int x1, int y1, int x2, int y2) { |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 49 | return static_cast<int>(getDistanceFloat(static_cast<float>(x1), static_cast<float>(y1), |
| 50 | static_cast<float>(x2), static_cast<float>(y2))); |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 51 | } |
| 52 | |
| 53 | static inline float getAngle(int x1, int y1, int x2, int y2) { |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 54 | const int dx = x1 - x2; |
| 55 | const int dy = y1 - y2; |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 56 | if (dx == 0 && dy == 0) return 0; |
| 57 | return atan2f(static_cast<float>(dy), static_cast<float>(dx)); |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 58 | } |
| 59 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 60 | static inline float getAngleDiff(float a1, float a2) { |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 61 | const float diff = fabsf(a1 - a2); |
Ken Wakasa | bcec82d | 2012-08-12 11:10:48 +0900 | [diff] [blame] | 62 | if (diff > M_PI_F) { |
| 63 | return 2.0f * M_PI_F - diff; |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 64 | } |
| 65 | return diff; |
| 66 | } |
| 67 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 68 | // static float pointToLineSegSquaredDistanceFloat( |
Ken Wakasa | bcec82d | 2012-08-12 11:10:48 +0900 | [diff] [blame] | 69 | // float x, float y, float x1, float y1, float x2, float y2) { |
| 70 | // float A = x - x1; |
| 71 | // float B = y - y1; |
| 72 | // float C = x2 - x1; |
| 73 | // float D = y2 - y1; |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 74 | // return fabsf(A * D - C * B) / sqrtf(C * C + D * D); |
Ken Wakasa | bcec82d | 2012-08-12 11:10:48 +0900 | [diff] [blame] | 75 | // } |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 76 | |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 77 | static inline float pointToLineSegSquaredDistanceFloat( |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 78 | float x, float y, float x1, float y1, float x2, float y2) { |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 79 | const float ray1x = x - x1; |
| 80 | const float ray1y = y - y1; |
| 81 | const float ray2x = x2 - x1; |
| 82 | const float ray2y = y2 - y1; |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 83 | |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 84 | const float dotProduct = ray1x * ray2x + ray1y * ray2y; |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 85 | const float lineLengthSqr = squareFloat(ray2x) + squareFloat(ray2y); |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 86 | const float projectionLengthSqr = dotProduct / lineLengthSqr; |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 87 | |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 88 | float projectionX; |
| 89 | float projectionY; |
| 90 | if (projectionLengthSqr < 0.0f) { |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 91 | projectionX = x1; |
| 92 | projectionY = y1; |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 93 | } else if (projectionLengthSqr > 1.0f) { |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 94 | projectionX = x2; |
| 95 | projectionY = y2; |
| 96 | } else { |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 97 | projectionX = x1 + projectionLengthSqr * ray2x; |
| 98 | projectionY = y1 + projectionLengthSqr * ray2y; |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 99 | } |
Ken Wakasa | 0fb9508 | 2012-08-13 17:31:46 +0900 | [diff] [blame^] | 100 | return getSquaredDistanceFloat(x, y, projectionX, projectionY); |
Satoshi Kataoka | 6b4a1d7 | 2012-08-10 15:42:56 +0900 | [diff] [blame] | 101 | } |
| 102 | } // namespace latinime |
Ken Wakasa | 507113a | 2012-08-13 11:36:41 +0900 | [diff] [blame] | 103 | #endif // LATINIME_GEOMETRY_UTILS_H |