基于pyexiv2的土工格栅

Question

我正在使用pyexiv2模块对JPEG进行地理定位，使用我在另一个中找到的代码，所以答案(请参阅：使用Python对jpeg图像进行地理标记的最佳方法是什么？)，我有一个关于GPSTag值的问题。

答案中给出的代码有以下几行：

exiv_image["Exif.Image.GPSTag"] = 654
exiv_image["Exif.GPSInfo.GPSMapDatum"] = "WGS-84"
exiv_image["Exif.GPSInfo.GPSVersionID"] = '2 0 0 0'

我已经查看了Exiv2文档，并找到了对GPSTag、GPSMapDatum和GPSVersionID的描述，但仍然对GPSTag的价值感到困惑。

从文件上说：

指向GPS信息IFD的指针。与Exif IFD一样，GPS信息IFD的互操作性结构没有图像数据。

这个描述并没有真正解释如何确定使用什么值，而且我还没有找到更好的GPSTag在线描述。

所以我的问题是：

1. 给定一个新的图像，如何确定Exif.Image.GPSTag的值？
2. 为什么代码示例使用的值为654 (这可能由问题一来回答)？

谢谢你的帮助。

Answer 1

使用pyexiv2进行地理标记照片的最佳方法肯定是使用我的程序，GottenGeography ;-)

但是说真的，如果您想从pyexiv2访问GPS数据，代码如下所示：

    GPS = 'Exif.GPSInfo.GPS'
    try:
        self.latitude = dms_to_decimal(
            *self.exif[GPS + 'Latitude'].value +
            [self.exif[GPS + 'LatitudeRef'].value]
        )
        self.longitude = dms_to_decimal(
            *self.exif[GPS + 'Longitude'].value +
            [self.exif[GPS + 'LongitudeRef'].value]
        )
    except KeyError:
        pass
    try:
        self.altitude = float(self.exif[GPS + 'Altitude'].value)
        if int(self.exif[GPS + 'AltitudeRef'].value) > 0:
            self.altitude *= -1
    except KeyError:
        pass

文字看起来是这样的：

    self.exif[GPS + 'AltitudeRef']  = '0' if self.altitude >= 0 else '1'
    self.exif[GPS + 'Altitude']     = Fraction(self.altitude)
    self.exif[GPS + 'Latitude']     = decimal_to_dms(self.latitude)
    self.exif[GPS + 'LatitudeRef']  = 'N' if self.latitude >= 0 else 'S'
    self.exif[GPS + 'Longitude']    = decimal_to_dms(self.longitude)
    self.exif[GPS + 'LongitudeRef'] = 'E' if self.longitude >= 0 else 'W'
    self.exif[GPS + 'MapDatum']     = 'WGS-84'

具备这些支助职能：

class Fraction(fractions.Fraction):
    """Only create Fractions from floats.

    >>> Fraction(0.3)
    Fraction(3, 10)
    >>> Fraction(1.1)
    Fraction(11, 10)
    """

    def __new__(cls, value, ignore=None):
        """Should be compatible with Python 2.6, though untested."""
        return fractions.Fraction.from_float(value).limit_denominator(99999)

def dms_to_decimal(degrees, minutes, seconds, sign=' '):
    """Convert degrees, minutes, seconds into decimal degrees.

    >>> dms_to_decimal(10, 10, 10)
    10.169444444444444
    >>> dms_to_decimal(8, 9, 10, 'S')
    -8.152777777777779
    """
    return (-1 if sign[0] in 'SWsw' else 1) * (
        float(degrees)        +
        float(minutes) / 60   +
        float(seconds) / 3600
    )


def decimal_to_dms(decimal):
    """Convert decimal degrees into degrees, minutes, seconds.

    >>> decimal_to_dms(50.445891)
    [Fraction(50, 1), Fraction(26, 1), Fraction(113019, 2500)]
    >>> decimal_to_dms(-125.976893)
    [Fraction(125, 1), Fraction(58, 1), Fraction(92037, 2500)]
    """
    remainder, degrees = math.modf(abs(decimal))
    remainder, minutes = math.modf(remainder * 60)
    return [Fraction(n) for n in (degrees, minutes, remainder * 60)]

尽管我目前正在研究一种替代pyexiv2的方法，它使用GObject内省来更直接地访问exiv2库，称为GExiv2，并希望得到一些反馈。gexiv2和pyexiv2都是围绕同一个exiv2库的包装器，但区别在于pyexiv2是一个非常大的项目，有大量的胶水，只在python中工作，处于放弃的边缘；而gexiv2是轻巧灵活的，可以从任何编程语言中访问，并且由于Shotwell的使用而得到了很好的维护。

希望这能有所帮助！

* pyexiv2's author, Olivier Tilloy, has asked me for help with maintainership as he no longer has much time

Answer 2

我的版本有点冗长..。

from fractions import Fraction
import pyexiv2

try:
    metadata = pyexiv2.metadata.ImageMetadata(image_file)
    metadata.read();
    thumb = metadata.exif_thumbnail

    try:
        latitude = metadata.__getitem__("Exif.GPSInfo.GPSLatitude")
        latitudeRef = metadata.__getitem__("Exif.GPSInfo.GPSLatitudeRef")
        longitude = metadata.__getitem__("Exif.GPSInfo.GPSLongitude")
        longitudeRef = metadata.__getitem__("Exif.GPSInfo.GPSLongitudeRef")

        latitude = str(latitude).split("=")[1][1:-1].split(" ");
        latitude = map(lambda f: str(float(Fraction(f))), latitude)
        latitude = latitude[0] + u"\u00b0" + latitude[1] + "'" + latitude[2] + '"' + " " + str(latitudeRef).split("=")[1][1:-1]

        longitude = str(longitude).split("=")[1][1:-1].split(" ");
        longitude = map(lambda f: str(float(Fraction(f))), longitude)
        longitude = longitude[0] + u"\u00b0" + longitude[1] + "'" + longitude[2] + '"' + " " + str(longitudeRef).split("=")[1][1:-1]

        latitude_value = dms_to_decimal(*metadata.__getitem__("Exif.GPSInfo.GPSLatitude").value + [metadata.__getitem__("Exif.GPSInfo.GPSLatitudeRef").value]);
        longitude_value = dms_to_decimal(*metadata.__getitem__("Exif.GPSInfo.GPSLongitude").value + [metadata.__getitem__("Exif.GPSInfo.GPSLongitudeRef").value]);

        print "--- GPS ---"
        print "Coordinates: " + latitude + ", " + longitude
        print "Coordinates: " + str(latitude_value) + ", " + str(longitude_value)
        print "--- GPS ---"
    except Exception, e:
        print "No GPS Information!"
        #print e

    # Check for thumbnail
    if(thumb.data == ""):
        print "No thumbnail!"
except Exception, e:
    print "Error processing image..."
    print e;