Convert Lat/Longs to X/Y Co-ordinates

by

You need more information than just a single lat/lon pair to be able to do this.

At this stage, the information you have provided is missing two things:

  • how large an area does your image cover (in terms of lat/lon)? Based on what you’ve provided, I don’t know if the image shows an area a metre wide or a kilometre wide.
  • what spot on your image does your reference coordinate (-37.803134, 145.132377) refer to? Is it one of the corners? In the middle somewhere?

I’m also going to assume that your image is aligned north/south – for example, it doesn’t have north pointing towards the top left corner. That would tend to complicate things.

The easiest approach is to figure out exactly what lat/lon coordinates correspond to the (0, 0) pixel and the (1017, 915) pixel. Then you can figure out the pixel corresponding to a given lat/lon coordinate through interpolation.

To briefly outline that process, imagine that your (-37.803134, 145.132377) lat/lon corresponds to your (0, 0) pixel, and that you’ve discovered that your (1017, 915) pixel corresponds to the lat/lon (-37.798917, 145.138535). Assuming the usual convention with pixel (0, 0) in the bottom-left corner, this means north is up in the image.

Then, if you are interested in the target coordinate (-37.801465, 145.134984), you can figure out the corresponding number of pixels up the image it is as follows:

pixelY = ((targetLat - minLat) / (maxLat - minLat)) * (maxYPixel - minYPixel)
       = ((-37.801465 - -37.803134) / (-37.798917 - -37.803134)) * (915 - 0)
       = 362.138

That is, the corresponding pixel is 362 pixels from the bottom of the image. You can then do exactly the same for the horizontal pixel placement, but using longitudes and X pixels instead.

The part ((targetLat - minLat) / (maxLat - minLat)) works out how far you are between the two reference coordinates, and gives 0 to indicate that you are at the first one, 1 to indicate the second one, and numbers in between to indicate locations in between. So for example, it would produce 0.25 to indicate that you are 25% of the way north between the two reference coordinates. The last bit converts that into the equivalent pixels.

HTH!

EDIT Ok, based on your comment I can be a little more specific. Given that you seem to be using the top left corner as your primary reference point, I’ll use the following definitions:

minLat = -37.803134
maxLat = -37.806232
MAP_HEIGHT = 916

Then, if we use an example coordinate of (-37.804465, 145.134984), the Y coordinate of the corresponding pixel, relative to the top-left corner, is:

pixelY = ((targetLat - minLat) / (maxLat - minLat)) * (MAP_HEIGHT - 1)
       = ((-37.804465 - -37.803134) / (-37.806232 - -37.803134)) * 915
       = 393.11

Therefore, the corresponding pixel is 393 pixels down from the top. I’ll let you work out the horizontal equivalent for yourself – it’s basically the same. NOTE The -1 with the MAP_HEIGHT is because if you start at zero, the maximum pixel number is 915, not 916.

EDIT: One thing I’d like to take the opportunity to point out is that this is an approximation. In reality, there isn’t a simple linear relationship between latitude and longitude coordinates and other forms of Cartesian coordinates for a number of reasons including the projections that are used while making maps, and the fact that the Earth is not a perfect sphere. In small areas this approximation is close enough as to make no significant difference, but on larger scales discrepancies may become evident. Depending on your needs, YMMV. (My thanks to uray, whose answer below reminded me that this is the case).

More Related Contents:

  1. What is the best way to create a sparse array in C++?
  2. iterator adapter to iterate just the values in a map?
  3. Why do i have a random result on this template?
  4. Resolve build errors due to circular dependency amongst classes
  5. Why isn’t sizeof for a struct equal to the sum of sizeof of each member?
  6. How do I iterate over cin line by line in C++?
  7. Testing stream.good() or !stream.eof() reads last line twice [duplicate]
  8. Move capture in lambda
  9. Non-const reference bound to temporary, Visual Studio bug?
  10. Is there a limit on number of open files in Windows
  11. Why does this function push RAX to the stack as the first operation?
  12. Using stdlib’s rand() from multiple threads
  13. Transform the modelMatrix
  14. std::stoi doesn’t exist in g++ 4.6.1 on MinGW
  15. Why does C need “struct” keyword and not C++?
  16. Vectors in Arduino
  17. Properly print utf8 characters in windows console
  18. Constexpr Math Functions
  19. How can I set the baud rate to 307,200 on Linux?
  20. Accessing private members [closed]
  21. array vs vector vs list
  22. CMake is not finding Boost
  23. Semicolon after class declaration braces
  24. GCC linker can’t find standard library?
  25. std::function and std::bind: what are they, and when should they be used?
  26. constexpr vs. static const: Which one to prefer?
  27. Performance difference of “if if” vs “if else if”
  28. Difference between A* pA = new A; and A* pA = new A();
  29. Local type as template arguments in C++ [duplicate]
  30. In lambda functions syntax, what purpose does a ‘capture list’ serve?

Leave a Comment