目 標(biāo) 場(chǎng) 景
有時(shí)候女朋友一個(gè)人在外面玩耍，問她在哪個(gè)地方，就是不告訴我。但是，你又很想知道女朋友的「位置」，這該如何是好？
其實(shí)你可以這樣套路女朋友，假裝自己在家很無(wú)聊，可以幫她修圖，讓她微信發(fā)原圖給你，拿到「微信原圖」后，你就可以利用 Python 快速獲取到女友的具體位置了。
準(zhǔn) 備 工 作
首先，在虛擬環(huán)境中安裝識(shí)別圖片元數(shù)據(jù)的庫(kù)。
pip3 install exifread
然后，進(jìn)入高德開放平臺(tái)，申請(qǐng)一個(gè) Web 服務(wù)的應(yīng)用，獲取到一個(gè) 「Key」用于逆地理編碼 API。
編 寫 腳 本
整個(gè)操作分為 3 步驟，分別是獲取圖片的經(jīng)度和緯度、對(duì)經(jīng)度和緯度進(jìn)行數(shù)據(jù)矯正、調(diào)用高德逆地理編碼 API 獲取具體位置。
第 1 步，獲取圖片的「經(jīng)度和緯度」。
使用 exifread 庫(kù)可以直接讀取圖片文件，獲取到圖片的元數(shù)據(jù)，包含經(jīng)度、緯度、南北緯方向、東西經(jīng)方向和拍攝時(shí)間。
# 使用 exifread 獲取圖片的元數(shù)據(jù)
img_exif = exifread.process_file(open(self.img_path, 'rb'))
# 能夠讀取到屬性
if img_exif:
# 緯度數(shù)
latitude_gps = img_exif['GPS GPSLatitude']
# N,S 南北緯方向
latitude_direction = img_exif['GPS GPSLatitudeRef']
# 經(jīng)度數(shù)
longitude_gps = img_exif['GPS GPSLongitude']
# E,W 東西經(jīng)方向
longitude_direction = img_exif['GPS GPSLongitudeRef']
# 拍攝時(shí)間
take_time = img_exif['EXIF DateTimeOriginal']
如果元數(shù)據(jù)存在，然后判斷拍攝時(shí)間是否合理。如果拍攝時(shí)間不在今天，那只能很遺憾地通知你，你的女朋友在向你撒謊「撒謊」。
def judge_time_met(self, take_time):
"""
判斷拍攝時(shí)間是否是在今天
:param take_time:
:return:
"""
# 拍攝時(shí)間
format_time = str(take_time).split(" ")[0].replace(":", "-")
# 當(dāng)天日期
today = str(datetime.date.today())
if format_time == today:
return False
else:
return True
if is_lie:
print('很遺憾的通知你，你的女朋友在撒謊?。?！')
return
如果女友沒有撒謊，那么可以進(jìn)行第 2 步的操作。
因?yàn)橥ㄟ^ GPS 獲取的經(jīng)度、緯度和高德地圖的坐標(biāo)存在一定的誤差，這里需要把坐標(biāo)轉(zhuǎn)換為「火星坐標(biāo)系」。
x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626 # π
a = 6378245.0 # 長(zhǎng)半軸
ee = 0.00669342162296594323 # 扁率
def wgs84togcj02(lng, lat):
"""
WGS84轉(zhuǎn)GCJ02(火星坐標(biāo)系)
:param lng:WGS84坐標(biāo)系的經(jīng)度
:param lat:WGS84坐標(biāo)系的緯度
:return:
"""
if out_of_china(lng, lat): # 判斷是否在國(guó)內(nèi)
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]
另外需要注意的是，接口中經(jīng)度、緯度參數(shù)只能識(shí)別小數(shù)點(diǎn)后 6 位，需要對(duì)經(jīng)緯度中的度、分、秒做一定的數(shù)據(jù)處理，然后再進(jìn)行四舍五入。
def __format_lati_long_data(self, data):
"""
對(duì)經(jīng)度和緯度數(shù)據(jù)做處理，保留6位小數(shù)
:param data: 原始經(jīng)度和緯度值
:return:
"""
# 刪除左右括號(hào)和空格
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只能識(shí)別到小數(shù)點(diǎn)后的6位
# 需要轉(zhuǎn)換為浮點(diǎn)數(shù)，并保留為6位小數(shù)
result = "%.6f" % (data_degree + data_minute + data_sec)
return float(result)
第 3 步，調(diào)用高德的反地理編碼 API，傳入申請(qǐng)的應(yīng)用 Key，就能拿到女朋友的詳細(xì)地址。
def __get_address(self, location):
"""
根據(jù)坐標(biāo)得到詳細(xì)地址
:param location: 經(jīng)緯度值
: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
結(jié) 果 結(jié) 論
確保圖片是原圖的基礎(chǔ)上，可以快速幫你判斷女朋友是否在撒謊；如果女朋友沒有撒謊，就返回女朋友具體的位置。
本文首發(fā)于公眾號(hào)「 AirPython 」，關(guān)注公眾號(hào)后，回復(fù)「 圖片定位 」即可獲得。
如果你覺得文章還不錯(cuò)，請(qǐng)大家點(diǎn)贊分享下。你的肯定是我最大的鼓勵(lì)和支持。
推薦閱讀：
