如何将纬度标准化为-90至90?

时间:2018-12-23 09:21:22

标签: python normalize piexif

使用piexif时,我试图将十进制和经度的DMS转换为十进制度,但是对于某些图像,我的纬度值为184.62583333333333,该值超出[-90,90]范围。

查看下面的代码,

exif_dict = piexif.load('images/DJI_0026.JPG')
long = 0
latt = 0
value = exif_dict['GPS']
if value:
    lat = value[2]
    lon = value[4]
    for i in range(3):
        if i == 1:
            latt += lat[i][0]/60.0
        elif i == 2:
            latt += lat[i][0]/3600.0
        else:
            latt += lat[i][0]
    for i in range(3):
        if i == 1:
            long += lon[i][0]/60.0
        elif i == 2:
            long += lon[i][0]/3600.0
        else:
            long += lon[i][0]
print(latt, long)

value = {0:(2,3,0,0),1:b'N',2:(((19,1),(8,1),(595773,10000)),3:b 'E',4:((73,1),(0,1),(131775,10000)),5:0,6:(70989,1000)}

我关心的是纬度和经度,它们存储在键2和4的值中。

纬度= 19 + 8 / 60.0 + 595773 / 3600.0

经度= 73 + 0 / 60.0 + 131775 / 3600.0

这就是输出。

输出:184.62583333333333 109.60416666666666

请让我知道如何规范[-90,90]范围内的纬度。

2 个答案:

答案 0 :(得分:2)

piexif返回的GPS坐标数据具有以下格式:

exif_data = {0: (2, 3, 0, 0),
             # Latitude: b'N' or b'S'
             1: b'N', 
             # deg, min, sec as (numerator,denominator) of rationals
             2: ((19, 1), (8, 1), (595773, 10000)),
             # Longitude: b'E' or b'W'
             3: b'E', 
             4: ((73, 1), (0, 1), (131775, 10000)), 
             5: 0, 6: (70989, 1000)}

纬度和经度是朝N或S(分别为E或W)给出的正的,有理的值。

我们需要将DMS的正值转换为十进制,然后根据方向给它们正确的符号:

def convert_DMS_tuple(tup):
    d, m, s = [t[0]/t[1] for t in tup]
    degrees = d + m/60 + s/3600
    return degrees

def EXIF_to_lat_long(exif_data):
    # Latitude is counted positive towards N
    lat_sign = 1 if exif_data[1] == b'N' else -1
    lat = lat_sign*convert_DMS_tuple(exif_data[2])

    # Longitude is counted positive towards E
    long_sign = 1 if exif_data[3] == b'E' else -1
    long = long_sign*convert_DMS_tuple(exif_data[4])

    return lat, long

EXIF_to_lat_long(exif_data)
# (19.149882583333333, 73.00366041666666)

您可以这样使用:

exif_dict = piexif.load('images/DJI_0026.JPG')
value = exif_dict['GPS']
if value:
   print(EXIF_to_lat_long(value))

答案 1 :(得分:1)

从给出的示例值中,很明显应该将val [0]除以val [1]。

例如595773/10000 = 59,这很有意义。

顺便说一句,您不需要for循环。它们使您的代码超出了必要的时间。

value[1] == b'N'可能是一个原因。您当前的代码虽然不会对其进行评估。这意味着您的代码仅适用于地球表面的1/4。

请参阅https://sno.phy.queensu.ca/~phil/exiftool/TagNames/GPS.html以获得良好的概述。您的代码也应该能够解释南方和西方。