Daily Python Blogger

通常Escape sequence都是表示特殊字元，如$, ?, ^...，所以在python中，要用\$, \?, \^，分別來表示，但是backquote本身該如何表示呢，\\即可，所以在regexp中，要寫成r'\\'，雖然己經用了raw表示法，但仍然還需再加一個backquote。

為了要找出最適合的善歌，用以下的程式來列舉出所有的mp3，然後再找出適合的歌曲。

from sys import stdin
from re import findall
data=stdin.read()
Files=findall(r'\\+\\(.*).mp3, data)
for file in Files:
     if not Temp.count(file): Temp.append(file)
for file in Temp: print file

# posted by samuel @ 12:50 AM 0 comments

LaTeX可以非常方便地，將TeX語法轉成數學圖形，於是模仿TextoGif的作法，且結合MoinMoin下的latexparser的圖形製作方式，自己改成了Tex2Png，主要是可以從stdin讀入latex的數學指令，存成Temp.tex，然後再裁切成png的圖檔，改寫的主要的步驟如下:

Tex.Tmpl內容如下:

\begin{eqnarray*}
$TexData
\end{eqnarray*}

1. 如stdin讀入eqnarray的指令，將中間的字串換成要顯示的數學指令。
2. 進入/tmp目錄夾
3. 使用latex -interaction=batchmode Temp.tex的指令，作出Temp.dvi的暫存檔。
4. dvips -E Temp.dvi -o Temp.ps，作出Temp.ps檔
5. convert -crop 0x0 -density 120x120 Temp.ps Temp.png，作出png的圖檔，同時也裁切圖檔，避免圖檔空白的地方太多

現在要再深入的是

1. dvips的指令
2. gs的指令: 因為textogif是使用gs -dNOPAUSE -r -sQutputFile=- -sDEVICE=pbmraw Temp.ps來達成

# posted by samuel @ 6:38 PM 0 comments

想要用python來處理影像，PIL(Python Imaging Library)，是蠻不錯的選擇。

PIL來自http://effbot.org/downloads是在Python社群中，專門處理圖片的Library，它有2種release版本，一為source，另一為binary，最新的版本分別如下：

• Source
• Binary

如果下載source，安裝後在$Imaging-xx.yy/Docs的目錄中，可以看到豐富的說明文件，及範例程式，文件內容包含了：系統需求、如何去編譯、安裝...。

1. 系統需求

2. 編譯及安裝

a. 先作出libImaging/ImConfig.h
[samuel@IBM Imaging-1.1.4]$ cd libImaging
    ./configure  ===> Makefile
    make         ===> ImConfig.h
b. 安裝
    [samuel@IBM Imaging-1.1.4]$sudo python setup.py install
c. 測試
    [samuel@sl304 ~]$ python
    Python 2.4 (#9, Jan 22 2005, 13:05:42)
    [GCC 3.2.3 20030502 (Red Hat Linux 3.2.3-42)] on linux2
    Type "help", "copyright", "credits" or "license" for more information.
    >>> import PIL
    >>>
如果沒有錯誤訊息，表示安裝成功。

3. 範例

1. 簡單測試
   >>> import Image
   >>> im = Image.open("lena.ppm")
   >>> print im.format, im.size, im.mode
   PPM (512, 512) RGB
   >>> im.show() #using 'xv', handy for debugging and testing
2. 開檔
    import Image
    im=Image.open('Conficius.gif')
    im.crop((360,150,700,745)).save('Conficius.png')
3. Crop Images
    box = (100, 100, 400, 400)      #cropping
    region = im.crop(box)
    region = region.transpose(Image.ROTATE_180)     #box processed and pasted back
    im.paste(region, box)
4. Rolling an image

    def roll(image, delta):
            "Roll an image sideways"
            xsize, ysize = image.size
            delta = delta % xsize
            if delta == 0: return image
            part1 = image.crop((0, 0, delta, ysize))
            part2 = image.crop((delta, 0, xsize, ysize))
            image.paste(part2, (0, 0, xsize-delta, ysize))
            image.paste(part1, (xsize-delta, 0, xsize, ysize))
            return image
5. Geometrical Transforms
    out = im.resize((128, 128))
    out = im.rotate(45) # degrees counter-clockwise
    out = im.transpose(Image.FLIP_LEFT_RIGHT)
    out = im.transpose(Image.FLIP_TOP_BOTTOM)
    out = im.transpose(Image.ROTATE_90)
    out = im.transpose(Image.ROTATE_180)
    out = im.transpose(Image.ROTATE_270)

# posted by samuel @ 4:59 PM 0 comments

Howto to Talk with GoogleTalk in Python

• GoogleTalk
• GTalk-dev

import xmpp
login = 'Your.Login' # @gmail.com
pwd   = 'YourPassword'
cnx = xmpp.Client('gmail.com')
cnx.connect( server=('talk.google.com',5223) )
cnx.auth(login,pwd, 'botty')
cnx.send( xmpp.Message( "YourFriend@gmail.com" ,"Hello World form Python" ) )

# posted by samuel @ 4:58 PM 0 comments

如果你有個list的結構如下:

L=['index.html', 'preface.html', 'd0e34.html', 'd0e37.html', 'd0e44.html', 'consul
ting.html', 'consulting.html', 'preface.html', 'preface.html', 'Audience.html',
'ChangesThirdEd.html', 'Acknowledgement.html', 'ToolsPart.html', 'intro.html', '
intro.html', 'OnlineResources.html', 'XSLprocessors.html', 'XSLprocessors.html',
'FOprocessors.html', 'Portability.html', 'ToolsSetup.html', 'ToolsSetup.html',
'ToolsSetup.html', 'ToolsSetup.html', 'ToolsSetup.html', 'InstallStylesheets.htm
l', 'InstallingAProcessor.html']

也就是說，有些值是duplicated，這時該怎何將此list的資料結構purified?

這可以應用在存取網路文件資料的範例中，例如在http://www.sagehill.net/docbookxsl，上頭有很多有用的文章，但它是以html的方式存在，並不是單一的檔案(如pdf的格式)，這對於如果想要讀取文件，但總是需要到同一網站來進行online方式的閱讀，最好的方式是將之全數下載，但對於python，要用何方式，來作到呢?

利用urllib2中的urlopen可以幫我們作到這一點。

#!/usr/bin/env python
from urllib2 import urlopen
from re import findall
from os import system, stat

def Purified(List):
        Temp=[]
        for value in List:
                if not Temp.count(value): Temp.append(value)
        return Temp

def RetDocxsl():
      try:
              SageFlag=stat('sagehill')
              if not SageFlag[6]:
                      data=urlopen('http://www.sagehill.net/docbookxsl/').read()
                      Webpages=findall('(.*?.html)',data)
                      Webpages=Purified(Webpages)
                      for page in Webpages:
                              system('wget --no-clobber http://63.249.72.176/docbookxsl/%(page)s'%vars())
                              #system('wget --no-clobber http://www.sagehill.net/docbookxsl/%(page)s'%vars())
              else:
                      data=open('sagehill').read()
                      Webpages=findall('(.*?.html)',data)
                      for page in Webpages:
                              system('wget --no-clobber http://63.249.72.176/docbookxsl/%(page)s'%vars())

      except:
              data=urlopen('http://www.sagehill.net/docbookxsl/').read()
              open('sagehill'),write(data)
              Webpages=findall('(.*?.html)',data)
              for page in Webpages:
                      system('wget --no-clobber http://63.249.72.176/docbookxsl/%(page)s'%vars())
if __name__=='__main__':
      RetDocxsl()

以上的程式碼，是利用urlopen將網址內容(index.html)讀取，並將文件內容存取成檔案(檔名為sagehill，方便下次來進行存取)，然後用re來parse所需要的文件網址位置。

利用wget來捉取檔案(--no-clobber選項的目的是，如果系統內有重複的檔案，就不要捉取具有相同檔名的文件，wget會略過，以利捉取檔案的速度)。

若是不用--no-clobber的選項，就要讓List中的檔案項目單一化(purified)，這一點也可以化作讀者的作業習題?!

此時可以利用getopt的方法，讓網址由getopt的選項，傳入去程式裡，使得程式本身更有彈性，加上利用urlparse，使得網址可以方便地找出來。這方面的便利性，可以作為程式改良的思考方向。

延伸問題: 如果是在windows中，沒有wget這個方便的util，如何用python幫忙將網頁資料取回??

>>> urlparse.urlparse('http://www.sagehill.net/docbookxsl/')
('http', 'www.sagehill.net', '/docbookxsl/', '', '', '')

# posted by samuel @ 4:57 PM 0 comments

如果在文件檔中，有一些是屬於character reference的文字如下:

<gentext key="Abstract"       text="&# x6458; &# x8981;"/>

#!/usr/bin/env python
# -*- coding: utf-8 -*-
import os,sys,re
from optparse import OptionParser

class XmlParser:
   def __init__(self, XmlFile):
       self.XmlData=open(XmlFile).read().decode('utf8')
   def PrintOut(self):
       from xml.dom.minidom import parseString
       self.XmlObj= parseString(self.XmlData)
       self.XmlStr=self.XmlObj.toxml(encoding='utf8')
       print self.XmlStr

if __name__=='__main__':
   usage='%prog [option]'
   parser=OptionParser(usage=usage,version='%prog, Reversion 0.1')
   parser.add_option('-X','--Xml',help='-X file.xml, parse XML')
   options,args=parser.parse_args()
   if len(sys.argv)==1: parser.print_help();raise SystemExit
   if options.Xml:
       xmlParser=XmlParser(options.Xml)
       xmlParser.PrintOut()

利用xml模組中的minidom，可以方便地將character reference文字轉換，google中相關的網站連結資料很豐富:

• unicode
• character reference
• characters
• 在emacs中轉換character reference
• vim中轉換
• Ascii vs Latin-1

# posted by samuel @ 4:53 PM 0 comments

對於python可以運用在教學情境裡，可以用一個很有趣的例子，來作這方面的應用，「時鐘追趕問題」：如果時間在二點到三點鐘之間，分針與時針什麼時候重合？

「重合」表示2針會成為一直線，這代表會有2次重疊的時刻，一次是短針和長針同方向重合在一起，另一次則是2針反方向而成為一直線的狀態。

現在想要更一般性地來表示這個問題，如果時刻是在h時m分的時候，請問短針、長針重合的時刻為何？

我們將證明，如果時間為h時m分，和時間為h時整點，其重合的時間是一致的。

證明：假設x為重合的時間，ExtraTime為額外需要的時間，如果題目是要2針夾角成Theta度，則換算成為Theta/6分；因為分針的速度為時針的12倍，所以下列式子A成立。

x-m=t
t=t/12+ExtraTime+(5h+m/12)-m   <=== 式子A
x=(60h+12*ExtraTime)/11 

ExtraTime為0的表示同方向重合，為30則表示反方向重合的狀態，剛好與時針在整點時，所算出的結果是一致的。

import sys
from optparse import OptionParser
from math import modf
def SynClockTime(*args):
  Hour, Min=args
  MinAngle, HourAngle=Min*6, Hour*30+Min/2.
  if HourAngle>MinAngle:      #MinAngle over the HourAngle
      for ExtraMin in (30, 60):
          AlignMin=(60*Hour+ExtraMin*12)/11.
          Sec, Min=modf(AlignMin)
          Sec, Min=Sec*60, int(Min)
          if Min<60: Min%=60; Hour+=1
          print 'The next alignment is= %(Hour)s:%(Min)s:%(Sec)s'%vars()
  else:
      for ExtraMin in (0, 30):
          AlignMin=(60*Hour+ExtraMin*12)/11.
          Sec, Min=modf(AlignMin)
          Sec, Min=Sec*60, int(Min)
          if Min>60: Min%=60; Hour+=1
          print 'The next alignment is= %(Hour)s:%(Min)s:%(Sec)s'%vars()

if __name__=='__main__':
  usage='%prog [option]'
  parser=OptionParser(usage=usage,version='%prog, Reversion 0.1')
  parser.add_option('-T','--Time',nargs=2,type='int',help='-T Hour(24-format) Min')
  options,args=parser.parse_args()
  if len(sys.argv)==1:parser.print_help();raise SystemExit
  if options.Time: SynClockTime(*options.Time)

程式碼分為2個情形，如果時針所走的角度比較大，則在同一整點時刻，會有2次重合的時間，如果分針所走的角度較大，則在同一整點及下一整點時刻，分別都有一次的重合時刻。

這個問題，可以幫我們理解，如何用python來表示，追趕問題中所需的數學式子。

延伸問題：
是否可以算至重合時間精確到秒?亦即為找出"hour時min分sec秒"，2針重合。

# posted by samuel @ 4:51 PM 0 comments

Teaching Math with Python

biblio

# posted by samuel @ 4:47 PM 0 comments

fac(n) = (n==0) ? 1 : n * fac(n-1)
fac(x) = (int(x)==0) ? 1.0 : int(x) * fac(int(x)-1.0)
stirling(x) = sqrt(2*pi*x) * x**x * exp(-x)
set xrange [1:15]
set yrange [1:1e+10]
set log y
set sample 15
plot stirling(x) notitle with lines,   fact(x)     notitle with points
#fact(n) ~ sqrt(2pi*x)x^x*exp(-x)

http://t16web.lanl.gov/Kawano/gnuplot/index-e.html

http://members.aol.com/demalion2/index.html http://www.blender3d.org/cms/Home.2.0.html

http://www.geom.uiuc.edu/~banchoff/projects.html torus

http://mathworld.wolfram.com/Torus.html

http://www.mmit.stc.sh.cn/telecenter/CnHisScience/yuanzhou.htm

http://www.mathland.idv.tw/ 昌爸工作坊

http://zh.wikipedia.org/wiki/%E5%9C%93%E5%91%A8%E7%8E%87 割圓術、圓周率

lanl Kawano

gnuplot是用來繪圖專用的軟體，它符合組合學習的方式，所以就來練習一下。

參考連結:

• 誰來上機
• gnuplot介紹

如果用攝氏和華氏的溫度轉換來學習，應該是很好的應用範例。

F=C*1.8+32，可以列出下列的數據對:

-49 -56.20
-42 -43.60
-35 -31.00
-28 -18.40
-21 -5.80
-14 6.80
-7 19.40
 0 32.00
 7 44.60
14 57.20
21 69.80
28 82.40
35 95.00
42 107.60

然後透過python令gnuplot來畫圖。

同樣的方式，也可以來畫2次函式圖，如繪出3x^2+9x-83的函式圖，比較可以直覺地教授2次函 式的特性如:頂點、向上vs向下、與x軸的交點、面積...>等等。

# posted by samuel @ 4:46 PM 0 comments