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

xml@python

Resource Description Framework (RDF)

Primer: Getting into RDF & Semantic Web using N3

xmltex

# posted by samuel @ 4:44 PM 0 comments

import win32com.client
strComputer = "."
objWMIService = win32com.client.Dispatch("WbemScripting.SWbemLocator")
objSWbemServices = objWMIService.ConnectServer(strComputer,"root\cimv2")
colItems = objSWbemServices.ExecQuery("Select * from CIM_VideoControllerResolution")
for objItem in colItems:
    print "Caption: ", objItem.Caption
    print "Description: ", objItem.Description
    print "Horizontal Resolution: ", objItem.HorizontalResolution
    print "Max Refresh Rate: ", objItem.MaxRefreshRate
    print "Min Refresh Rate: ", objItem.MinRefreshRate
    print "Number Of Colors: ", objItem.NumberOfColors
    print "Refresh Rate: ", objItem.RefreshRate
    print "Scan Mode: ", objItem.ScanMode
    print "Setting ID: ", objItem.SettingID
    print "Vertical Resolution: ", objItem.VerticalResolution


import win32com.client
strComputer = "."
objWMIService = win32com.client.Dispatch("WbemScripting.SWbemLocator")
objSWbemServices = objWMIService.ConnectServer(strComputer,"root\cimv2")
colItems = objSWbemServices.ExecQuery("Select * from Win32_DisplayConfiguration")
for objItem in colItems:
    print "Bits Per Pel: ", objItem.BitsPerPel
    print "Caption: ", objItem.Caption
    print "Description: ", objItem.Description
    print "Device Name: ", objItem.DeviceName
    print "Display Flags: ", objItem.DisplayFlags
    print "Display Frequency: ", objItem.DisplayFrequency
    print "Dither Type: ", objItem.DitherType
    print "Driver Version: ", objItem.DriverVersion
    print "ICM Intent: ", objItem.ICMIntent
    print "ICM Method: ", objItem.ICMMethod
    print "Log Pixels: ", objItem.LogPixels
    print "Pels Height: ", objItem.PelsHeight
    print "Pels Width: ", objItem.PelsWidth
    print "Setting ID: ", objItem.SettingID
    print "Specification Version: ", objItem.SpecificationVersion


import win32com.client
strComputer = "."
objWMIService = win32com.client.Dispatch("WbemScripting.SWbemLocator")
objSWbemServices = objWMIService.ConnectServer(strComputer,"root\cimv2")
colItems = objSWbemServices.ExecQuery("Select * from Win32_DesktopMonitor")
for objItem in colItems:
    print "Availability: ", objItem.Availability
    print "Bandwidth: ", objItem.Bandwidth
    print "Caption: ", objItem.Caption
    print "Config Manager Error Code: ", objItem.ConfigManagerErrorCode
    print "Config Manager User Config: ", objItem.ConfigManagerUserConfig
    print "Creation Class Name: ", objItem.CreationClassName
    print "Description: ", objItem.Description
    print "Device ID: ", objItem.DeviceID
    print "Display Type: ", objItem.DisplayType
    print "Error Cleared: ", objItem.ErrorCleared
    print "Error Description: ", objItem.ErrorDescription
    print "Install Date: ", objItem.InstallDate
    print "Is Locked: ", objItem.IsLocked
    print "Last Error Code: ", objItem.LastErrorCode
    print "Monitor Manufacturer: ", objItem.MonitorManufacturer
    print "Monitor Type: ", objItem.MonitorType
    print "Name: ", objItem.Name
    print "Pixels Per X Logical Inch: ", objItem.PixelsPerXLogicalInch
    print "Pixels Per Y Logical Inch: ", objItem.PixelsPerYLogicalInch
    print "PNP Device ID: ", objItem.PNPDeviceID
    z = objItem.PowerManagementCapabilities
    if z is None:
        a = 1
    else:
        for x in z:
            print "Power Management Capabilities: ", x
    print "Power Management Supported: ", objItem.PowerManagementSupported
    print "Screen Height: ", objItem.ScreenHeight
    print "Screen Width: ", objItem.ScreenWidth
    print "Status: ", objItem.Status
    print "Status Info: ", objItem.StatusInfo
    print "System Creation Class Name: ", objItem.SystemCreationClassName
    print "System Name: ", objItem.SystemName


import win32com.client
strComputer = "."
objWMIService = win32com.client.Dispatch("WbemScripting.SWbemLocator")
objSWbemServices = objWMIService.ConnectServer(strComputer,"root\cimv2")
colItems = objSWbemServices.ExecQuery("Select * from Win32_DisplayControllerConfiguration")
for objItem in colItems:
    print "Bits Per Pixel: ", objItem.BitsPerPixel
    print "Caption: ", objItem.Caption
    print "Color Planes: ", objItem.ColorPlanes
    print "Description: ", objItem.Description
    print "Device Entries In A Color Table: ", objItem.DeviceEntriesInAColorTable
    print "Device Specific Pens: ", objItem.DeviceSpecificPens
    print "Horizontal Resolution: ", objItem.HorizontalResolution
    print "Name: ", objItem.Name
    print "Refresh Rate: ", objItem.RefreshRate
    print "Reserved System Palette Entries: ", objItem.ReservedSystemPaletteEntries
    print "Setting ID: ", objItem.SettingID
    print "System Palette Entries: ", objItem.SystemPaletteEntries
    print "Vertical Resolution: ", objItem.VerticalResolution
    print "Video Mode: ", objItem.VideoMode


import win32com.client
strComputer = "."
objWMIService = win32com.client.Dispatch("WbemScripting.SWbemLocator")
objSWbemServices = objWMIService.ConnectServer(strComputer,"root\cimv2")
colItems = objSWbemServices.ExecQuery("Select * from Win32_VideoController")
for objItem in colItems:
    z = objItem.AcceleratorCapabilities
    if z is None:
        a = 1
    else:
        for x in z:
            print "Accelerator Capabilities: ", x
    print "Adapter Compatibility: ", objItem.AdapterCompatibility
    print "Adapter DAC Type: ", objItem.AdapterDACType
    print "Adapter RAM: ", objItem.AdapterRAM
    print "Availability: ", objItem.Availability
    z = objItem.CapabilityDescriptions
    if z is None:
        a = 1
    else:
        for x in z:
            print "Capability Descriptions: ", x
    print "Caption: ", objItem.Caption
    print "Color Table Entries: ", objItem.ColorTableEntries
    print "Config Manager Error Code: ", objItem.ConfigManagerErrorCode
    print "Config Manager User Config: ", objItem.ConfigManagerUserConfig
    print "Creation Class Name: ", objItem.CreationClassName
    print "Current Bits Per Pixel: ", objItem.CurrentBitsPerPixel
    print "Current Horizontal Resolution: ", objItem.CurrentHorizontalResolution
    print "Current Number Of Colors: ", objItem.CurrentNumberOfColors
    print "Current Number Of Columns: ", objItem.CurrentNumberOfColumns
    print "Current Number Of Rows: ", objItem.CurrentNumberOfRows
    print "Current Refresh Rate: ", objItem.CurrentRefreshRate
    print "Current Scan Mode: ", objItem.CurrentScanMode
    print "Current Vertical Resolution: ", objItem.CurrentVerticalResolution
    print "Description: ", objItem.Description
    print "Device ID: ", objItem.DeviceID
    print "Device Specific Pens: ", objItem.DeviceSpecificPens
    print "Dither Type: ", objItem.DitherType
    print "Driver Date: ", objItem.DriverDate
    print "Driver Version: ", objItem.DriverVersion
    print "Error Cleared: ", objItem.ErrorCleared
    print "Error Description: ", objItem.ErrorDescription
    print "ICM Intent: ", objItem.ICMIntent
    print "ICM Method: ", objItem.ICMMethod
    print "Inf Filename: ", objItem.InfFilename
    print "Inf Section: ", objItem.InfSection
    print "Install Date: ", objItem.InstallDate
    print "Installed Display Drivers: ", objItem.InstalledDisplayDrivers
    print "Last Error Code: ", objItem.LastErrorCode
    print "Max Memory Supported: ", objItem.MaxMemorySupported
    print "Max Number Controlled: ", objItem.MaxNumberControlled
    print "Max Refresh Rate: ", objItem.MaxRefreshRate
    print "Min Refresh Rate: ", objItem.MinRefreshRate
    print "Monochrome: ", objItem.Monochrome
    print "Name: ", objItem.Name
    print "Number Of Color Planes: ", objItem.NumberOfColorPlanes
    print "Number Of Video Pages: ", objItem.NumberOfVideoPages
    print "PNP Device ID: ", objItem.PNPDeviceID
    z = objItem.PowerManagementCapabilities
    if z is None:
        a = 1
    else:
        for x in z:
            print "Power Management Capabilities: ", x
    print "Power Management Supported: ", objItem.PowerManagementSupported
    print "Protocol Supported: ", objItem.ProtocolSupported
    print "Reserved System Palette Entries: ", objItem.ReservedSystemPaletteEntries
    print "Specification Version: ", objItem.SpecificationVersion
    print "Status: ", objItem.Status
    print "Status Info: ", objItem.StatusInfo
    print "System Creation Class Name: ", objItem.SystemCreationClassName
    print "System Name: ", objItem.SystemName
    print "System Palette Entries: ", objItem.SystemPaletteEntries
    print "Time Of Last Reset: ", objItem.TimeOfLastReset
    print "Video Architecture: ", objItem.VideoArchitecture
    print "Video Memory Type: ", objItem.VideoMemoryType
    print "Video Mode: ", objItem.VideoMode
    print "Video Mode Description: ", objItem.VideoModeDescription
    print "Video Processor: ", objItem.VideoProcessor

# posted by samuel @ 4:42 PM 0 comments

Python & Webservices

1. tempeststrings
2. pywebsvcs
3. schema

pymedia

http://mail.python.org/pipermail/python-announce-list/2004-June/003160.html

http://pymedia.org/
http://pymedia.org/tut/index.html

cvs -d:pserver:anonymous@cvs.sourceforge.net:/cvsroot/pymedia login
cvs -z3 -d:pserver:anonymous@cvs.sourceforge.net:/cvsroot/pymedia co -P pymedia

PyMedia library is a Python module for wav, mp3, ogg, avi, divx, dvd, cdda etc files manipulations. It allows you to parse, demutiplex, multiplex, decode and encode all supported formats. Also it features direct interface to a DVD-ROM drive allowing you to extract Audio disks, DVD and VCD movies. Sound interface is making it a snap to play sound on a selected audio card. Resample and reformat your sound before you play it. Features page can help you to see what else you can do with PyMedia.

PyMedia is very small in size and very flexible at the same time. See tutorial for numerous list of examples. It allows you to create your own mutimedia applications in a matter of minutes by even a newbie in programming. Python language is choosen because of simple semantics, complete and reach set of features.

There is a PyCar application which takes advantage of the PyMedia library extensivelly and serve in my car as a powerful car media center.

secure directory

GPS written in Python

Tkinter Usage

# posted by samuel @ 4:40 PM 0 comments

VPython Link

1. Bounce
2. gas simulation tutorial

Python@Physics&Science -- Py@Psi

一直想要將物理與Programming Technique連結，但是無法順利地將夢想實現，總算在cyberworld中，有這樣想法的人還是存在的，我想將這部份要發展的software，在sourceforge上開發專案的Codename叫作Py@Psi。

Py@Psi名稱即取自python的諧音，同時也是取自數學的重要常數pi的發音，另外Physics&Science的諧音為Psi，Psi是在量子物理中很重要的函數代號，此2者合起來叫作Py@Psi，表示在Sourceforge上開發的專案，是專注在將物理與科學上的現象，呈現出圖象化，現階段先來介紹vpython。

Visual Python

為了要2D/3D的繪圖加速，vpython使用boost來將python的程式，用C++包裝起來(Wrapping)，且使用POV-RAY的繪圖、OpenGL的繪圖方式，來將物理定律呈現。

CVS snapshot

cvs -d:pserver:anonymous@cvs.sourceforge.net:/cvsroot/visualpython login
(just Enter)
cvs -z3 -d:pserver:anonymous@cvs.sourceforge.net:/cvsroot/visualpython co vpython

但是在安裝CVS snapshot，無法找到configure，因為那是要自己作出來的，所以就安裝下列的套件，試試看可否安裝成功。

autotools-dev libtool autoconf automake1.4 m4

最快的方式，就是安裝python-visual，同時也會安裝了下列的套件

gtkglarea5 libboost-python1.33.0 python-visual

Visual Python links

1. VRML Gallery of EM
2. vpython
3. mathworld.wolfram/LogPlot
4. vpython/contributed
5. Esalgado's vpython/

Other Links

1. Matter&Interaction
2. Ruth Chabay's POV-Ray
3. cT archives

# posted by samuel @ 4:31 PM 0 comments