天天看點

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

想要悄悄的擷取某人的位置,隻需通過她拍攝的照片,就能可以實作。此方法需要用到識别圖檔中繼資料的庫。這裡操作是在虛拟機ubuntu上,使用python最好是在3.6及以上。

pip3 install exifread           

1. 粗擷取

擷取拍攝照片的經緯度,時間,裝置資訊。這裡使用的照片需是拍攝的原圖,圖檔路徑本文是放在了與運作程式相同的檔案夾。實作代碼如下:

import exifread

# 擷取照片位置
class Image_Location():

    def __init__(self, image_path):
        self.img_path = image_path

    def lati_long_date(self, data):
        """
        對經緯度進行處理,保留6位小數
        """

        # 删除左右括号,最後以逗号分隔為一個清單
        data_list_tmp = str(data) .replace('[', '').replace(']', '').split(',')

        # 循環取出每個元素,删除元素兩邊的空格,得到一個新清單
        data_list = [date.strip() for date in data_list_tmp]

        # 替換秒的值
        data_tmp = data_list[-1].split('/')

        # 秒的值
        data_sec = int(data_tmp[0]) / int(data_tmp[1]) / 3600

        # 替換分的值
        data_tmp = data_list[-2]

        # 分的值
        data_minute = int(data_tmp) / 60

        # 度的值
        data_degree = int(data_list[0])

        # 由于高德API隻能識别到小數點後的6位
        # 需要轉換為浮點數,并保留為6位小數
        result = "%.6f" % (data_degree + data_minute + data_sec)
        return float(result)

    def get_image_exif(self):
        """
        擷取圖檔的屬性:緯度,經度,拍攝時間等
        """

        # 通過exifread擷取圖檔的屬性
        img_read = open(self.img_path, 'rb')
        img_exif = exifread.process_file(img_read)

        # 讀取到的屬性
        if img_exif:
            for exif in img_exif:

                # 緯度
                if exif == "GPS GPSLatitude":
                    latitude = img_exif["GPS GPSLatitude"]

                # 緯度的方向
                elif exif == "GPS GPSLatitudeRef":
                    latitude_direction = img_exif["GPS GPSLatitudeRef"]

                # 經度
                elif exif == "GPS GPSLongitude":
                    longitude = img_exif["GPS GPSLongitude"]

                # 經度方向
                elif exif == "GPS GPSLongitudeRef":
                    longitude_direction = img_exif["GPS GPSLongitudeRef"]

                # 拍攝時間
                elif exif == "EXIF DateTimeDigitized":
                    take_time = img_exif["EXIF DateTimeDigitized"]

                # 拍攝裝置
                elif exif == "Image Make":
                    take_equipment = img_exif["Image Make"]

                # 拍攝型号
                elif exif == "Image Model":
                    take_model = img_exif["Image Model"]

            # 對擷取的經緯度資訊進一步處理
            if latitude and longitude:
                latitude_date = self.lati_long_date(latitude)
                longitude_date = self.lati_long_date(longitude)
                print("緯度為:", latitude_date)
                print("經度為:", longitude_date)
            print("拍攝時間為:", take_time)
            print("拍攝裝置為: ", take_equipment)
            print("拍攝型号為:", take_model)

        else:
            print("圖像資訊為空,可能是上傳的不是原圖")


if __name__ == '__main__':
    location = Image_Location('li.jpg')
    location.get_image_exif()           

運作結果:

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

可通過GPS查詢網址:

http://www.gpsspg.com/maps.htm

,查詢具體的位置,如下圖:

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

查詢的位置與實際位置有一點誤差,不要擔心下面有更準确,更智能的方法。

2.精擷取

這種方法需要先在高德開放平台上申請一個 Web 服務的應用,獲得一個秘鑰Key 用于逆地理編碼 API。

進入高德開放平台,選擇【開發支援】,【Web服務】

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

選擇【擷取Key】,就可以看到擷取Key的步驟,這裡需要注冊一個高德的賬号。

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

選擇進入【控制台】,就可以建立一個應用,按照上述步驟就可以擷取Key了。

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

擷取Key界面。

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊

整個操作分為三步,首先是擷取圖檔的經緯度,然後對經緯度進行資料矯正,最後調用高德逆地理編碼API擷取具體位置。

因為GPS擷取的經緯度和高德地圖的坐标存在一定的誤差,這裡需要把坐标轉換為【火星坐标系】。

轉換代碼如下:

# transfer.py
# -*- coding: utf-8 -*-
import json
import math

x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 長半軸
ee = 0.00669342162296594323  # 扁率


def wgs84togcj02(lng, lat):
    """
    WGS84轉GCJ02(火星坐标系)
    :param lng:WGS84坐标系的經度
    :param lat:WGS84坐标系的緯度
    :return:
    """
    if out_of_china(lng, lat):  # 判斷是否在國内
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [mglng, mglat]


def gcj02towgs84(lng, lat):
    """
    GCJ02(火星坐标系)轉GPS84
    :param lng:火星坐标系的經度
    :param lat:火星坐标系緯度
    :return:
    """
    if out_of_china(lng, lat):
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [lng * 2 - mglng, lat * 2 - mglat]


def transformlat(lng, lat):
    ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \
        0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng))
    ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lat * pi) + 40.0 *
            math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
    ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *
            math.sin(lat * pi / 30.0)) * 2.0 / 3.0
    return ret


def transformlng(lng, lat):
    ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \
        0.1 * lng * lat + 0.1 * math.sqrt(math.fabs(lng))
    ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *
            math.sin(2.0 * lng * pi)) * 2.0 / 3.0
    ret += (20.0 * math.sin(lng * pi) + 40.0 *
            math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
    ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *
            math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
    return ret


def out_of_china(lng, lat):
    """
    判斷是否在國内,不在國内不做偏移
    :param lng:
    :param lat:
    :return:
    """
    if lng < 72.004 or lng > 137.8347:
        return True
    if lat < 0.8293 or lat > 55.8271:
        return True
    return False
           

擷取照片資訊代碼,這裡需要将Key,換成你在高德開放平台擷取的Key:

import os
import exifread
from decimal import Decimal
from transfer import *
import requests
import json
import datetime

#  pip3 install exifread
class Location(object):
    def __init__(self, image_path):
        self.img_path = image_path
        self.api_key = "将你的Key粘貼到這裡"
        self.url_get_position = 'https://restapi.amap.com/v3/geocode/regeo?key={}&location={}'

    def run(self):
        coordinate, time = self.__get_image_ability()
        print(f'拍攝圖檔的緯度為:{coordinate}')
        print(f'拍攝照片的時間為:{time}')
        if not coordinate:
            return

        # 根據經度和緯度,擷取到詳細位址
        address = self.__get_address(coordinate)

        # 檢驗坐标值
        # https://lbs.amap.com/console/show/picker
        print(f'拍攝圖檔的位址為:{address}')

    def __get_address(self, location):
        """
        根據坐标得到詳細位址
        :param location: 經緯度值
        :return:
        """
        resp = requests.get(
            self.url_get_position.format(self.api_key, location))
        location_data = json.loads(resp.text)
        address = location_data.get('regeocode').get('formatted_address')
        return address

    def __format_lati_long_data(self, data):
        """
        對經度和緯度資料做處理,保留6位小數
        :param data: 原始經度和緯度值
        :return:
        """
        # 删除左右括号和空格
        data_list_tmp = str(data).replace('[', '').replace(']', '').split(',')
        data_list = [data.strip() for data in data_list_tmp]

        # 替換秒的值
        data_tmp = data_list[-1].split('/')

        # 秒的值
        data_sec = int(data_tmp[0]) / int(data_tmp[1]) / 3600

        # 替換分的值
        data_tmp = data_list[-2]

        # 分的值
        data_minute = int(data_tmp) / 60

        # 度的值
        data_degree = int(data_list[0])

        # 由于高德API隻能識别到小數點後的6位
        # 需要轉換為浮點數,并保留為6位小數
        result = "%.6f" % (data_degree + data_minute + data_sec)
        return float(result)

    def __get_image_ability(self):
        """
        擷取圖檔的屬性值,包含:經緯度、拍攝時間等
        :param picture_name:
        :return:
        """
        # 利用exifread庫,讀取圖檔的屬性
        img_exif = exifread.process_file(open(self.img_path, 'rb'))

        # 能夠讀取到屬性
        if img_exif:
            # 緯度數
            latitude_gps = img_exif['GPS GPSLatitude']

            # N,S 南北緯方向
            latitude_direction = img_exif['GPS GPSLatitudeRef']

            # 經度數
            longitude_gps = img_exif['GPS GPSLongitude']

            # E,W 東西經方向
            longitude_direction = img_exif['GPS GPSLongitudeRef']

            # 拍攝時間
            take_time = img_exif['EXIF DateTimeOriginal']

            # 拍攝裝置
            take_equipment = img_exif["Image Make"]

            # 拍攝型号
            take_model = img_exif["Image Model"]

            # 緯度、經度、拍攝時間
            if latitude_gps and longitude_gps and take_time and take_equipment and take_model:

                # 對緯度、經度值原始值作進一步的處理
                latitude = self.__format_lati_long_data(latitude_gps)
                longitude = self.__format_lati_long_data(longitude_gps)

                print(f'拍攝的裝置型号為:{take_equipment} {take_model}')
                # print(f'{longitude},{latitude}')
                # 注意:由于gps擷取的坐标在國内高德等主流地圖上逆編碼不夠精确,這裡需要轉換為火星坐标系
                location = wgs84togcj02(longitude, latitude)
                return f'{location[0]},{location[1]}', take_time
            else:
                print(f'擷取的圖檔資料屬性不完整')
                return ''
        else:
            print('抱歉,圖檔不是原圖,沒法擷取到圖檔屬性。')
            return ''

if __name__ == '__main__':
    # 她發過來的圖檔【原圖】
    location = Location('li.jpg')

    # 找到她的地理位置
    location.run()
           

執行代碼就可以得到,你想知道的了,一張照片就可以知道對方的位置。

python 通過圖檔(原圖)精确擷取圖檔拍攝的位置,時間,裝置等資訊