任何人都可以将此60位值转换为地理坐标吗？

问题描述：

我有一个包含位置信息的（DB2）数据库表，其中一列是CHARACTER（16）并包含一个单一的十六进制数字。我有一个程序将这些显示为地图上的点，但无法访问其来源。改变数字会移动地图上的点 - 我只是不知道算法。任何人都可以将此60位值转换为地理坐标吗？

一些例子（编辑以包括更多）：

04867C279DE2D6EC -32.063657° 115.7658683° 
04867C27C030085E -32.0633982° 115.7649085° 
04867C27C230A5FE -32.0633846° 115.7653336°

是否有这样做的任何标准的方法？这是一些DB2约定还是什么？关于如何解决这个问题的任何想法？

谢谢！

还有一组将使这个问题交代。 – msw 2010-06-28 06:35:14

已添加。但是，我不确定那些坐标是否完全准确;可能会有一些四舍五入的情况发生。 – 2010-06-28 06:50:37

答

这可能是morton codes。他们用于包装坐标的一个维。

生成莫顿号码是很容易。你所做的只是将x和y坐标转换为二进制。然后“交错”这些比特以获得莫顿数。先走哪个并不重要，但你必须保持一致。的

是的！就是这个！这些值正在被成功分割，并且当一个被修改并且被重新组合时，地图上的点仅在一个维度上移动。所以，现在我只需要弄清楚他们指的是什么网格。 – 2010-07-01 05:59:00

答

也许它是一个UTM网格格式？

http://www.uwgb.edu/dutchs/usefuldata/utmformulas.htm

有很多在UTM变量当然（例如所使用的基准）。

可能是，尤其是MGA是UTM投影。我会进一步调查。 – 2010-06-28 06:52:18

答

好了，只是想解开各种方式的十六进制值...

#! /usr/bin/python 
import struct 
import binascii 

a = '04867C279DE2D6EC' 
b = '04867C27C030085E' 
c = '04867C27C230A5FE' 

formats = ['2I', '2i', '2f', 'd', '4h', '4H'] 
formats += ['>'+item for item in formats] 

for fmt in formats: 
    print fmt, '-->' 
    for item in [a,b,c]: 
     coords = struct.unpack(fmt, binascii.unhexlify(item)) 
     print ' ', coords

取得了一些想法......

2I --> 
    (662472196, 3973505693) 
    (662472196, 1577595072) 
    (662472196, 4272238786) 
2i --> 
    (662472196, -321461603) 
    (662472196, 1577595072) 
    (662472196, -22728510) 
2f --> 
    (3.5044675291578432e-15, -2.07824221089183e+27) 
    (3.5044675291578432e-15, 2.4533886735682109e+18) 
    (3.5044675291578432e-15, -1.0978789217059195e+38) 
d --> 
    (-1.9722947342913136e+216,) 
    (9.4395557694675488e+144,) 
    (-1.135288151092706e+302,) 
4h --> 
    (-31228, 10108, -7523, -4906) 
    (-31228, 10108, 12480, 24072) 
    (-31228, 10108, 12482, -347) 
4H --> 
    (34308, 10108, 58013, 60630) 
    (34308, 10108, 12480, 24072) 
    (34308, 10108, 12482, 65189) 
>2I --> 
    (75922471, 2648889068) 
    (75922471, 3224373342) 
    (75922471, 3257968126) 
>2i --> 
    (75922471, -1646078228) 
    (75922471, -1070593954) 
    (75922471, -1036999170) 
>2f --> 
    (3.1617264522911893e-36, -6.0043925910101893e-21) 
    (3.1617264522911893e-36, -2.7505106925964355) 
    (3.1617264522911893e-36, -44.162101745605469) 
>d --> 
    (7.3832340678903009e-287,) 
    (7.3832347392384709e-287,) 
    (7.383234778429458e-287,) 
>4h --> 
    (1158, 31783, -25118, -10516) 
    (1158, 31783, -16336, 2142) 
    (1158, 31783, -15824, -23042) 
>4H --> 
    (1158, 31783, 40418, 55020) 
    (1158, 31783, 49200, 2142) 
    (1158, 31783, 49712, 42494)

拆开它作为一个大端无符号32位整数（ > 2I）看起来有点像它可能被投射某种坐标......

(-32.063657, 115.7658683) --> (75922471, 2648889068) 
(-32.0633982, 115.7649085) --> (75922471, 3224373342) 
(-32.0633846, 115.7653336) --> (75922471, 3257968126)

如果是在UTM，它的可能是UTM Zone 50S，基于lat，long ...我们不知道基准，但是这不应该超过几百米的差距。

有关区域50S：

(lat, long) --> (Easting, Northing) 
(-32.063657, 115.7658683) --> (383506.31320936838, 6451842.2821839228) 
(-32.0633982, 115.7649085) --> (383415.3800562254, 6451869.9348384682) 
(-32.0633846, 115.7653336) --> (383455.49221963808, 6451871.9016738012)

（使用OSR的（可怕unpythonic）Python包装......）

from osgeo import osr 
def latlong2utm(lat, long): 
    epsg_wgs84 = 4326 
    epsg_utm50S = 32750 
    inproj = osr.SpatialReference() 
    inproj.ImportFromEPSG(epsg_wgs84) 
    outproj = osr.SpatialReference() 
    outproj.ImportFromEPSG(epsg_utm50S) 
    transform = osr.CoordinateTransformation(inproj, outproj) 
    x,y,_ = transform.TransformPoint(long, lat) 
    return x, y

不幸的是，它看起来并不像UTM区50S，反正.. 。和区域49不看任何好转...

所以，没有太大的帮助，但我想我会张贴此帮助谁可能试图在同一路径其他人......（我也试过拆包IBM转而采用o的东西f IEEE漂浮...没有运气那里要么...）也许我只是愚蠢的，但我没有想法。希望这有助于一点，但...

使用OSR +1。尽管Python包装是......奇怪的...我仍然发现它是交互式坐标系变换的一个令人难以置信的强大工具。 – 2010-06-29 21:10:54

任何人都可以将此60位值转换为地理坐标吗？

相关推荐