Month: August 2010

前兩天被這個東西給陰到了，只怪我沒搞清楚，這在 Learning python 2nd 裡的 21.1 裡寫的很清楚。不過我想，從Java/C++/C#跳槽過來的人應該都會搞糊塗吧~

以C++舉例，下面的code：
[cpp]class Object {
private:
int id;
};
[/cpp]

到了 python，應該都會很自然的寫成：
[python]# 例A
class Object:
id=0
[/python]

看起來沒錯，對吧~
不過，事實上，應該要寫成這樣才對：
[python]# 例 B
class Object:
def init(self):
self.id=0
[/python]

換句話說，例A翻譯成C++的話，其實是類似 static 的用法，也就是說，id 會變成一個共用的 static 變數：
[cpp]class Object {
private:
static int id;
};
[/cpp]

要存GPS的相關tag到圖片時，因為某些tag要求存Rational或GPSCoordinate，我不知道該怎麼轉換，所以花了不少功夫查。我基本上是參考 這個老先生的code，他主要是使用 GPS tracker，所以會要讀取 tracker 裡的紀錄，然後以時間跟照片做比對，再把 GPS 位置寫進圖片裡。

可是遇到兩個問題：

1. 他的 code 是用 surd 處理 Rational，我找不到 surd 了。
2. 算出來的度、分、秒不準確。

首先是解決Rational問題，python 2.6 已經有內置的 module – Fraction 可以解決了。

def R(f):
"""R(float) - get a Rational number for a float"""
from fractions import Fraction
from pyexiv2 import Rational
r = Fraction.from_float( float(f) )
return Rational( r.numerator, r.denominator )

再來就是轉換的問題，從 Google gear 那邊拿到的經緯度，要轉換成度分秒，我後來是參考 上河文化網站-座標轉換程式與 Google Maps 地圖定位裡的 javascript 來改寫，這才得到正確的度、分、秒：

def d(angle):
"""d(any) - get degrees from a number :eg d(33.41) = 33"""
return math.floor(angle)
def m(angle):
"""m(any) - get minutes from a number :eg d(33.41) = 24"""
return math.floor( (angle-d(angle))*60 )
def s(angle):
"""s(any) - get seconds from a number :eg s(33.41) = 36"""
return round((angle-d(angle)-m(angle)/60)*3600)
def degree2coordinate( lat, lon ):
"""Convert degree to coordinate string."""
latR    = 'N'
lonR    = 'E'
if lat  < 0:
lat = -lat
latR= 'S'
if lon  < 0:
lon = -lon
lonR= 'W'
slat = "%03d,%02d,%02d%s" % (d(lat),m(lat),s(lat),latR )
slon = "%03d,%02d,%02d%s" % (d(lon),m(lon),s(lon),lonR )
return (pyexiv2.GPSCoordinate.from_string( slat ), pyexiv2.GPSCoordinate.from_string( slon ) )

最後，記一下，Tag裡要存的資料型態：

• Xmp.exif.GPSLatitude、Xmp.exif.GPSLongitude 裡存的是 pyexiv2.GPSCoordinate
• Xmp.exif.GPSMapDatum、Exif.GPSInfo.GPSMapDatum裡存的是字串，我是用 WGS-84
• Exif.GPSInfo.GPSLatitude、Exif.GPSInfo.GPSLongitude 裡存的是度、分、秒的tuple，要注意，度、分、秒都要轉成 Rational
• Exif.GPSInfo.GPSLatitudeRef、Exif.GPSInfo.GPSLongitudeRef 裡存的是度、分、秒的表示字串，如：125.3,67.5,12E

在之後打算要寫入 geotagging 到圖片的時候，才發現 python XMP toolkit 的不足，於是又踏上尋覓的路程。後來找到 pyexiv2，這個 library 更好用，不僅可以讀，也能寫，用法也不會太複雜。唯一的缺點是，ubuntu 10.04 包的是 1.x 版，好在，有人包了最新的版本在ppa裡。
安裝：

sudo apt-add-repository ppa:pyexiv2-developers/ppa
sudo apt-get install python-pyexiv2

讀取跟寫入，基本上參考官方的tutorial與API 文件就行了。Cell ID 資訊在這幾個tag裡：

• Xmp.cell.mcc
• Xmp.cell.mnc
• Xmp.cell.lac
• Xmp.cell.cellid

去年就看到這篇文章：Cell ID Photo PC Tool 0.3，想試試看了。因為是c#寫的，照理來說，Linux 上的 mono 應該可以用才對，哪知道卻發現他用了 System.Windows.Media.Imaging 這個 mono 沒實作的類別庫。
沒辦法，於是就想說用 python 寫一個好了，當作練習。
首先遇到的問題就是 cell id 怎麼取得，原本以為是在 exif 裏面，用了 PIL 去撈，卻沒撈到。胡亂找了一陣，才知道這是放在 XMP metadata 裏面，得用 XML toolkit sdk 來撈才行，python 也是有 binding：Python XMP toolkit，只是得另外裝，ubuntu 沒有包進去。

是故，得先裝 libexempi3 (sudo apt-get install libexempi3)，再下載 Python XMP toolkit，解開後，執行 sudo python setup.py install 來安裝。這樣就大功告成了。

程式很簡單，簡單到爆：

#!/usr/bin/env python
import sys
from libxmp import file_to_dict
CELL_NS="http://developer.sonyericsson.com/cell/1.0/"
if __name__ == '__main__':
xmp = file_to_dict( sys.argv[1] )
for prop_name, prop_value, prop_options in xmp[CELL_NS]:
print prop_name, prop_value

就這樣就可以把 Cell ID 相關的資訊都撈出來了。