1. python讀取圖片exif屬性中的GPS信息

• 智能手機(jī)或平板如果在拍照時(shí)開啟定位服務(wù)，照片中就會(huì)記錄拍照位置信息和拍攝時(shí)間。如果將原始照片直接發(fā)送發(fā)布到網(wǎng)上，無(wú)意中就泄漏了自己的位置信息，有惡意企圖的人可能會(huì)通過(guò)照片分析出你的家庭住址和工作單位。
• python有很多工具庫(kù)可以解析圖片的exif元數(shù)據(jù)信息，筆者喜歡使用exifread這個(gè)庫(kù)。

2. 完整代碼（需要使用python3）

• 圖片exif屬性中的經(jīng)緯度是“度，分，秒”的形式，如：[25, 34, 5927/100]，需要轉(zhuǎn)換成類似“16.439683,39.941649999”這樣的形式。
• 除經(jīng)緯度外，還要讀取GPS元數(shù)據(jù)中的東西半球和南北半球標(biāo)識(shí)。
• 經(jīng)緯度轉(zhuǎn)換成地理位置需要調(diào)用地圖服務(wù)接口，本程序代碼采用百度的逆地理編碼服務(wù)API接口，調(diào)用接口需要申請(qǐng)服務(wù)密鑰。
• 筆者自己利用手機(jī)拍照測(cè)試，程序分析后得到的位置信息基本沒有誤差。

#!/bin/python #coding:utf8 import os import exifread import re import sys import requests import json__author__ = 'DaDaLuLa'#************************************************************************ #代碼功能： # # 1.讀取所有圖片文件的exif信息 # # 2.提取圖片中的經(jīng)緯度，將度、分、秒轉(zhuǎn)換為小數(shù)形式 # # 3.利用百度地圖API接口將經(jīng)緯度轉(zhuǎn)換成地址形式 # #************************************************************************#遍歷文件夾及子文件夾中的所有圖片,逐個(gè)文件讀取exif信息 def get_pic_GPS(pic_dir): items = os.listdir(pic_dir) for item in items: path = os.path.join(pic_dir, item) if os.path.isdir(path): get_pic_GPS(path) else: imageread(path)#將經(jīng)緯度轉(zhuǎn)換為小數(shù)形式 def convert_to_decimal(*gps):#度if '/' in gps[0]:deg = gps[0].split('/')if deg[0] == '0' or deg[1] == '0':gps_d = 0else:gps_d = float(deg[0]) / float(deg[1])else:gps_d = float(gps[0])#分if '/' in gps[1]:minu = gps[1].split('/')if minu[0] == '0' or minu[1] == '0':gps_m = 0else:gps_m = (float(minu[0]) / float(minu[1])) / 60else:gps_m = float(gps[1]) / 60#秒if '/' in gps[2]:sec = gps[2].split('/')if sec[0] == '0' or sec[1] == '0':gps_s = 0else:gps_s = (float(sec[0]) / float(sec[1])) / 3600else:gps_s = float(gps[2]) / 3600decimal_gps = gps_d + gps_m + gps_s#如果是南半球或是西半球if gps[3] == 'W' or gps[3] == 'S' or gps[3] == "83" or gps[3] == "87":return str(decimal_gps * -1)else:return str(decimal_gps)#讀取圖片的經(jīng)緯度和拍攝時(shí)間 def imageread(path):f = open(path,'rb')GPS = {}Data = ""try:tags = exifread.process_file(f)except:return '''for tag in tags: print(tag,":",tags[tag])'''#南北半球標(biāo)識(shí)if 'GPS GPSLatitudeRef' in tags:GPS['GPSLatitudeRef'] = str(tags['GPS GPSLatitudeRef'])else:GPS['GPSLatitudeRef'] = 'N' #缺省設(shè)置為北半球#東西半球標(biāo)識(shí)if 'GPS GPSLongitudeRef'in tags:GPS['GPSLongitudeRef'] = str(tags['GPS GPSLongitudeRef'])else:GPS['GPSLongitudeRef'] = 'E' #缺省設(shè)置為東半球#海拔高度標(biāo)識(shí)if 'GPS GPSAltitudeRef' in tags:GPS['GPSAltitudeRef'] = str(tags['GPS GPSAltitudeRef'])#獲取緯度if 'GPS GPSLatitude' in tags:lat = str(tags['GPS GPSLatitude'])#處理無(wú)效值if lat == '[0, 0, 0]' or lat == '[0/0, 0/0, 0/0]':returndeg, minu, sec = [x.replace(' ', '') for x in lat[1:-1].split(',')]#將緯度轉(zhuǎn)換為小數(shù)形式GPS['GPSLatitude'] = convert_to_decimal(deg, minu, sec,GPS['GPSLatitudeRef'])#獲取經(jīng)度if 'GPS GPSLongitude' in tags:lng = str(tags['GPS GPSLongitude'])#處理無(wú)效值if lng == '[0, 0, 0]' or lng == '[0/0, 0/0, 0/0]':returndeg, minu, sec = [x.replace(' ', '') for x in lng[1:-1].split(',')]#將經(jīng)度轉(zhuǎn)換為小數(shù)形式GPS['GPSLongitude'] = convert_to_decimal(deg, minu, sec,GPS['GPSLongitudeRef'])#對(duì)特殊的經(jīng)緯度格式進(jìn)行處理#獲取海拔高度if 'GPS GPSAltitude' in tags:GPS['GPSAltitude'] = str(tags["GPS GPSAltitude"])#獲取圖片拍攝時(shí)間if 'Image DateTime' in tags:GPS["DateTime"] = str(tags["Image DateTime"])elif "EXIF DateTimeOriginal" in tags:GPS["DateTime"] = str(tags["EXIF DateTimeOriginal"])if 'GPSLatitude' in GPS:#將經(jīng)緯度轉(zhuǎn)換為地址convert_gps_to_address(GPS)#利用百度全球逆地理編碼服務(wù)（Geocoder）Web API接口服務(wù)將經(jīng)緯轉(zhuǎn)換為位置信息 def convert_gps_to_address(GPS):secret_key = 'xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx'#百度密鑰lat, lng = GPS['GPSLatitude'], GPS['GPSLongitude']#注意coordtype為wgs84ll(GPS經(jīng)緯度),否則定位會(huì)出現(xiàn)偏差baidu_map_api = "http://api.map.baidu.com/geocoder/v2/?coordtype=wgs84ll&location={0},{1}&output=json&pois=0&ak={2}".format(lat,lng,secret_key)content = requests.get(baidu_map_api).textgps_address = json.loads(content)#結(jié)構(gòu)化的地址formatted_address = gps_address["result"]["formatted_address"]#國(guó)家（若需訪問(wèn)境外POI，需申請(qǐng)逆地理編碼境外POI服務(wù)權(quán)限）country = gps_address["result"]["addressComponent"]["country"]#省province = gps_address["result"]["addressComponent"]["province"]#城市city = gps_address["result"]["addressComponent"]["city"]#區(qū)district = gps_address["result"]["addressComponent"]["district"]#語(yǔ)義化地址描述sematic_description = gps_address["result"]["sematic_description"]#將轉(zhuǎn)換后的信息寫入文件 with open("gps_address.csv","a+") as csv:csv.write(GPS["DateTime"] + "|" + formatted_address + "|" + country + "|" + province + "|" + city + "|" + district + "|" + sematic_description + "\n") if __name__ == "__main__":get_pic_GPS("./photo/")

• 從網(wǎng)站下載多張照片，測(cè)試運(yùn)行結(jié)果如下：
  2017:07:27 22:02:52|天津市河?xùn)|區(qū)長(zhǎng)征路8號(hào)|中國(guó)|天津市|天津市|河?xùn)|區(qū)|水利家園東79米
  2018:10:28 17:53:32|廣東省廣州市天河區(qū)員村二橫路|中國(guó)|廣東省|廣州市|天河區(qū)|程界子富新村內(nèi)
  2017:08:30 05:52:54|山東省濟(jì)南市濟(jì)陽(yáng)縣經(jīng)四路|中國(guó)|山東省|濟(jì)南市|濟(jì)陽(yáng)縣|嘉景苑南

3. 某色情網(wǎng)站照片GPS批量提取

• 筆者無(wú)意中發(fā)現(xiàn)，互聯(lián)網(wǎng)上傳播的一些色情圖片保留了完整的GPS等元數(shù)據(jù)信息。為了分析色情圖片的拍攝位置，從某網(wǎng)站爬取了34萬(wàn)張少兒不宜的圖片。利用上文中的程序進(jìn)行分析，共獲取了3000余張圖片包含的經(jīng)緯度信息，利用百度地圖API接口轉(zhuǎn)換成對(duì)應(yīng)的精確地理位置信息。
  photo.jpg

• 經(jīng)程序分析后獲取的部分?jǐn)?shù)據(jù)如下（多數(shù)位置信息能夠精確所在的小區(qū)和酒店,以下結(jié)果對(duì)具體地址加*號(hào)處理）：
  2017:10:16 15:13:43|四川省成都市雙流縣華陽(yáng)大道4段-200|中國(guó)|四川省|成都市|雙流縣|成都栢頓*大酒店內(nèi)
  2017:10:16 15:43:39|四川省成都市雙流縣華陽(yáng)大道4段-200|中國(guó)|四川省|成都市|雙流縣|泊爾*酒店
  2017:10:24 19:41:37|四川省成都市雙流縣華陽(yáng)大道4段-200|中國(guó)|四川省|成都市|雙流縣|泊爾*酒店
  2017:08:31 22:04:58|山東省濰坊市諸城市|中國(guó)|山東省|濰坊市|諸城市|萬(wàn)*家園內(nèi)
  2017:08:31 22:05:20|山東省濰坊市諸城市|中國(guó)|山東省|濰坊市|諸城市|萬(wàn)*家園內(nèi)
  2017:08:31 22:04:29|山東省濰坊市諸城市|中國(guó)|山東省|濰坊市|諸城市|萬(wàn)*家園內(nèi)
  2017:08:31 22:07:36|山東省濰坊市諸城市|中國(guó)|山東省|濰坊市|諸城市|萬(wàn)*家園內(nèi)
  2018:01:24 23:21:23|河南省許昌市長(zhǎng)葛市鐵東路|中國(guó)|河南省|許昌市|長(zhǎng)葛市|幸*小區(qū)
  2018:01:24 23:21:20|河南省許昌市長(zhǎng)葛市鐵東路|中國(guó)|河南省|許昌市|長(zhǎng)葛市|幸*小區(qū)
  018:08:31 22:56:52|廣東省深圳市寶安區(qū)航城大道|中國(guó)|廣東省|深圳市|寶安區(qū)|南航明*花園內(nèi)
  2018:11:14 21:43:42|廣東省深圳市寶安區(qū)航城大道|中國(guó)|廣東省|深圳市|寶安區(qū)|南航明*花園內(nèi)
  2018:06:21 21:43:45|遼寧省沈陽(yáng)市于洪區(qū)中央大街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|沈陽(yáng)世紀(jì)高爾*俱樂(lè)部?jī)?nèi)
  2018:08:23 19:04:38|遼寧省沈陽(yáng)市大東區(qū)滂江街|中國(guó)|遼寧省|沈陽(yáng)市|大東區(qū)|龍之*暢園內(nèi)
  2018:10:24 23:16:13|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水墨*青內(nèi)
  2018:10:24 23:12:40|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水墨*青內(nèi)
  2018:10:24 23:34:06|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水墨*青內(nèi)
  2018:10:24 23:12:30|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水墨*青內(nèi)
  2017:10:10 07:37:31|廣東省深圳市寶安區(qū)澳桂路|中國(guó)|廣東省|深圳市|寶安區(qū)|尚品*居內(nèi)
  2017:10:10 07:37:36|廣東省深圳市寶安區(qū)澳桂路|中國(guó)|廣東省|深圳市|寶安區(qū)|尚品*居內(nèi)
  2017:10:21 20:28:06|江蘇省徐州市云龍區(qū)大工巷1號(hào)樓2樓|中國(guó)|江蘇省|徐州市|云龍區(qū)|建*小區(qū)內(nèi)
  2017:10:21 20:22:25|江蘇省徐州市云龍區(qū)大工巷1號(hào)樓2樓|中國(guó)|江蘇省|徐州市|云龍區(qū)|建*小區(qū)內(nèi)
  2017:11:07 16:02:51|廣東省深圳市寶安區(qū)澳桂路|中國(guó)|廣東省|深圳市|寶安區(qū)|尚品*居內(nèi)
  2017:11:07 16:02:50|廣東省深圳市寶安區(qū)澳桂路|中國(guó)|廣東省|深圳市|寶安區(qū)|尚品*居內(nèi)
  2018:04:04 15:40:21|上海市寶山區(qū)春雷路372號(hào)|中國(guó)|上海市|上海市|寶山區(qū)|上海中*醫(yī)院內(nèi)
  2018:04:04 15:40:14|上海市寶山區(qū)春雷路372號(hào)|中國(guó)|上海市|上海市|寶山區(qū)|上海中*醫(yī)院內(nèi)
  2017:10:16 22:42:54|浙江省溫州市永嘉縣|中國(guó)|浙江省|溫州市|永嘉縣|西后村附近37米
  2018:11:14 21:44:13|廣東省深圳市寶安區(qū)寶安大道5005-9|中國(guó)|廣東省|深圳市|寶安區(qū)|匯*居內(nèi)
  2013:05:13 00:20:17|四川省成都市成華區(qū)建設(shè)北路2段-4號(hào)|中國(guó)|四川省|成都市|成華區(qū)|電子科技大學(xué)(沙河校區(qū))內(nèi)
  2017:10:30 11:47:33|廣東省汕頭市龍湖區(qū)長(zhǎng)平路東段|中國(guó)|廣東省|汕頭市|龍湖區(qū)|汕頭龍光喜來(lái)*酒店內(nèi),龍光世紀(jì)大廈南59米
  2017:10:30 11:35:11|廣東省汕頭市龍湖區(qū)金環(huán)南路|中國(guó)|廣東省|汕頭市|龍湖區(qū)|汕頭龍光喜來(lái)*酒店內(nèi),龍光世紀(jì)大廈內(nèi)0米
  2017:07:28 23:11:27|山東省濟(jì)南市濟(jì)陽(yáng)縣緯一路|中國(guó)|山東省|濟(jì)南市|濟(jì)陽(yáng)縣|雅*園-四區(qū)內(nèi)
  2017:07:29 05:20:51|山東省濟(jì)南市濟(jì)陽(yáng)縣緯一路|中國(guó)|山東省|濟(jì)南市|濟(jì)陽(yáng)縣|雅*園-四區(qū)內(nèi)
  2017:07:28 23:12:03|山東省濟(jì)南市濟(jì)陽(yáng)縣緯一路|中國(guó)|山東省|濟(jì)南市|濟(jì)陽(yáng)縣|雅*園-四區(qū)內(nèi)
  2018:01:09 13:14:27|河北省滄州市任丘市燕山南道|中國(guó)|河北省|滄州市|任丘市|源*美璟商業(yè)廣場(chǎng)內(nèi),源*大酒店內(nèi)
  2017:12:31 11:28:04|浙江省杭州市江干區(qū)秋濤北路38號(hào)|中國(guó)|浙江省|杭州市|江干區(qū)|浙江大學(xué)醫(yī)學(xué)院附屬邵逸夫醫(yī)院（慶春院區(qū)）西南298米
  2017:12:31 11:24:09|浙江省杭州市江干區(qū)秋濤北路52號(hào)3023|中國(guó)|浙江省|杭州市|江干區(qū)|杭州錢*精品酒店內(nèi)
  2018:01:13 10:47:13|浙江省杭州市上城區(qū)浣紗路18號(hào)|中國(guó)|浙江省|杭州市|上城區(qū)|浙江煙草大樓東54米
  2017:12:31 11:30:59|浙江省杭州市江干區(qū)秋濤北路52號(hào)3023|中國(guó)|浙江省|杭州市|江干區(qū)|杭州錢*精品酒店內(nèi)
  2018:06:21 21:43:45|遼寧省沈陽(yáng)市于洪區(qū)中央大街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|沈陽(yáng)世紀(jì)高*俱樂(lè)部?jī)?nèi)
  2018:08:23 19:04:38|遼寧省沈陽(yáng)市大東區(qū)滂江街|中國(guó)|遼寧省|沈陽(yáng)市|大東區(qū)|龍之*暢園內(nèi)
  2018:10:24 23:16:13|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水*丹青內(nèi)
  2018:10:24 23:12:40|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水*丹青內(nèi)
  2018:10:24 23:34:06|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水*丹青內(nèi)
  2018:10:24 23:12:30|遼寧省沈陽(yáng)市于洪區(qū)十一號(hào)街|中國(guó)|遼寧省|沈陽(yáng)市|于洪區(qū)|水*丹青內(nèi)
  2017:08:07 21:47:14|江蘇省蘇州市張家港市小河壩西路|中國(guó)|江蘇省|蘇州市|張家港市|尚*國(guó)際內(nèi)
  2017:05:01 23:17:23|江蘇省蘇州市張家港市小河壩西路|中國(guó)|江蘇省|蘇州市|張家港市|尚*國(guó)際內(nèi)
  2017:05:01 22:56:49|江蘇省蘇州市張家港市小河壩西路|中國(guó)|江蘇省|蘇州市|張家港市|尚*國(guó)際內(nèi)
  2017:05:01 23:13:07|江蘇省蘇州市張家港市小河壩西路|中國(guó)|江蘇省|蘇州市|張家港市|尚*國(guó)際內(nèi)
  2017:09:15 09:29:44|陜西省西安市未央?yún)^(qū)浐灞二路|中國(guó)|陜西省|西安市|未央?yún)^(qū)|滹*小區(qū)西69米
  2017:09:15 09:25:36|陜西省西安市未央?yún)^(qū)浐灞二路|中國(guó)|陜西省|西安市|未央?yún)^(qū)|泘*小學(xué)內(nèi)
  2017:09:15 09:27:12|陜西省西安市未央?yún)^(qū)浐灞二路|中國(guó)|陜西省|西安市|未央?yún)^(qū)|滹*錦繡-北區(qū)內(nèi)
  2017:09:15 09:25:38|陜西省西安市未央?yún)^(qū)浐灞二路|中國(guó)|陜西省|西安市|未央?yún)^(qū)|泘沱小學(xué)內(nèi)
  2018:08:23 23:06:26|福建省南平市建陽(yáng)市黃溪路93號(hào)|中國(guó)|福建省|南平市|建陽(yáng)市|金*大酒店(上水南路店)內(nèi)
  2017:09:15 21:53:39|四川省成都市成華區(qū)和錦路|中國(guó)|四川省|成都市|成華區(qū)|成都*畔生活酒店(成都火車東站魅力店)
  2017:06:30 23:22:51|福建省福州市閩侯縣廣賢路|中國(guó)|福建省|福州市|閩侯縣|福建華南*職業(yè)學(xué)院(旗山校區(qū))內(nèi)
  2017:12:04 00:32:06|福建省福州市倉(cāng)山區(qū)永南路|中國(guó)|福建省|福州市|倉(cāng)山區(qū)|領(lǐng)*新城內(nèi)
  2017:08:06 10:07:41|四川省成都市新都區(qū)詹家灣路|中國(guó)|四川省|成都市|新都區(qū)|潤(rùn)*花園內(nèi)

4. 拍攝位置分析

• 拍攝地分析

利用python的科學(xué)庫(kù)numpy、數(shù)據(jù)分析庫(kù)panda、繪圖庫(kù)matplotlib和數(shù)據(jù)可視化庫(kù)pyecharts,對(duì)獲取的不良圖片拍攝位置數(shù)據(jù)進(jìn)行處理，分析不良圖片拍攝地分布情況。

本文進(jìn)對(duì)數(shù)據(jù)進(jìn)行粗略分析，展示結(jié)果僅作一般性觀察分析，不代表筆者任何觀點(diǎn)和傾向。

#代碼在notebook中執(zhí)行 %matplotlib inline import pandas as pd import matplotlib.pyplot as plt import numpy as np plt.rcParams['font.sans-serif'] = ['Simhei'] df = pd.read_csv('gps_address.csv',sep='|') fig = plt.figure(figsize=(15,20)) data = df["province"].value_counts() num_of_province = data.values plt.barh(range(len(num_of_province),0,-1),num_of_province,height=0.7,color='steelblue',alpha=0.8) plt.title("各省分布情況",fontsize=20) plt.xlabel("數(shù)量",fontsize=15) plt.ylim(0,len(num_of_province)+1) plt.yticks(range(len(num_of_province),0,-1),data.index,fontsize=15) for x,y in enumerate(np.sort(num_of_province)):plt.text(y + 0.5,x + 0.9,'%s' % y,fontsize=18 ) plt.show()

?

province.png

fig = plt.figure(figsize=(15,40)) data = df["city"].value_counts() num_of_city = data.values plt.barh(range(len(num_of_city),0,-1),num_of_city,height=0.8,color='steelblue',alpha=0.8) plt.title("城市分布情況",fontsize=20) plt.xlabel("數(shù)量",fontsize=15) plt.yticks(range(len(num_of_city),0,-1),data.index,fontsize=15) plt.ylim(0,len(num_of_city)+2) for x,y in enumerate(np.sort(num_of_city)):plt.text(y +0.2,x + 0.8,'%s' % y,fontsize=18 ) plt.show()

?

city.png

• 各省分布情況地圖展示

from pyecharts import Map, Geo import re provinces = list(df["province"].value_counts().index) provinces = [re.sub("壯族自治區(qū)|自治區(qū)|省|市","",x) for x in provinces] print(provinces)pro_values = list(df["province"].value_counts().values) city = list(df["province"].value_counts().index) city_values = list(df["province"].value_counts().values)geo = Geo("各省分布情況", "", title_color="#fff",title_pos="center", width=1000,height=600, background_color='#404a59') geo.add("", provinces, pro_values, visual_range=[0, 200], maptype='china',visual_text_color="#fff",symbol_size=10, is_visualmap=True)geo

?

各省分布情況.png

• 城市分布情況地圖展示

city = list(df["city"].value_counts().index) city= [re.sub("白族自治州|市","",x) for x in city] city_values = list(df["city"].value_counts().values) geo = Geo("城市分布情況", "", title_color="#fff",title_pos="center", width=1000,height=600, background_color='#404a59') geo.add("", city, city_values, visual_range=[0, 200], maptype='china',visual_text_color="#fff", symbol_size=10, is_visualmap=True) geo

?

城市分布情況.png

總結(jié)

以上是生活随笔為你收集整理的某少儿不宜网站图片拍摄位置分析,Python批量读取图片GPS位置！的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。