带调色板TIFF到Jpeg的转化Python脚本

发现东西存成文件放到电脑里边超级难找，找不到就成了破烂，说不定下次就被Format掉了。最近赶51前赶快把机器里不积累的小东西发到BLOG上来存档。
下边是很早转DRG数据时写的Python脚本都快在电脑里边沉下去了。

 1  #  功能：用于将调色板GEOTIFF转化为JPEG并获取元数据.py
 2  #  作者：贠建明
 3  #  日期：2007-1-22
 4 
 5  import  os,sys,string
 6 
 7  #  遍历目录获取TIF文件名
 8  def  GetTiffFileName(dir):
 9      filenames  =  []
10      
11      dirlist = os.listdir(dir)
12       for  subdir  in  dirlist:
13          subdir  =  os.path.join(dir, subdir)
14           if  os.path.isdir(subdir):
15              subfilelist = os.listdir(subdir)
16               for  filename  in  subfilelist:
17                   if  filename[ - 4 :].lower() == " .tif " :
18                      filenames.append(os.path.join(subdir, filename))
19           else :
20              filenames.append(subdir)
21 
22      filelist = file( " c:/filelist.txt " , ' w ' )
23       for  filename  in  filenames:
24          filelist.write(filename  +   " \n " )
25 
26      filelist.close()        
27       return  filenames    
28 
29                 
30  #  获取元数据到文件
31  def  GetMetadata(filenames):
32      metainfo  =  file( ' c:/metainfo.txt ' , ' w ' )
33       for  filename  in  filenames:
34          os.system( " gdalinfo > C:/info.txt  "   +  filename)   
35          info  =  file( ' C:/info.txt ' , ' r ' )
36           while  True:
37              line  =  info.readline()
38               if  len(line)  ==  0:
39                   break
40              
41               if  line[: 4 ]  ==   ' Size ' :
42                   # 文件地址
43                  metainfo.write( ' \n ' )
44                  i  =  filename.rfind( ' \\ ' ) + 1
45                  path  =   " c:\\JPG\\ "   +  filename[i: - 4 ]  +   " .jpg "               
46                  metainfo.write(path)
47                  metainfo.write( ' , ' )
48                   # 分幅号
49                  noOrigin  =  filename[ - 17 : - 7 ]
50                  intNO  =  ord(noOrigin[0]) - ord( ' A ' )  +   1   
51                  mapNo  =  str(intNO) + noOrigin[ 1 : 3 ] + noOrigin[ 4 :]
52                  metainfo.write(mapNo)
53                  metainfo.write( ' , ' )
54                   #  列行号
55                  metainfo.write(line[ 8 :len(line) - 1 ])
56               if  line[: 10 ]  ==   ' Upper Left ' :
57                  metainfo.write( ' , ' )
58                  metainfo.write(line[ 14 :len(line) - 3 ])
59               if  line[: 11 ]  ==   ' Lower Right ' :
60                  metainfo.write( ' , ' )
61                  metainfo.write(line[ 14 :len(line) - 3 ])
62                  metainfo.write( ' , ' )
63                   #  投影带号
64                  strip  =  string.atoi(filename[ - 16 : - 14 ])  # 'H:\DRG\I48E009021\I48E009021DRG.tif'[-16:-14]='48'
65                  strip  =  strip - 30
66                  metainfo.write(str(strip))
67 
68          info.close()    
69 
70      metainfo.close();
71      
72       for  filename  in  filenames:
73          os.system( " gdalinfo >> C:/AllInfo.txt  "   +  filename)    
74 
75 
76  #  进行TIFF到Jpeg的转化
77  def  GeoTiff2Jpeg(filenames):
78       for  filename  in  filenames:
79          os.system( " call C:\FWTools1.0.5\setfw.bat & pct2rgb  " + filename + "  c:/rgb.tif " )
80          i  =  filename.rfind( ' \\ ' ) + 1
81          os.system( " gdal_translate c:/rgb.tif -of JPEG c:/JPG/ "   +  filename[i: - 4 ]  +   " .jpg " )
82 
83      
84  if   __name__ == " __main__ " :
85      files  =  GetTiffFileName( ' I:\\5万实验数据\\DRG ' )
86       # GetMetadata(files)
87      GeoTiff2Jpeg(files)
88