gamit-工程处理流程－详细

gamit-工程处理流程－详细.txt

工程说明：工程名为：db05。数据来自中国大陆7个IGS站数据，其中kunm、urum、lhaz三个作为已知站，其他4站有：wuhn,shao,bjfs,guao,计算日期从2009年5月1日到2009年5月4日。年纪日（121～124），以下为具体操作流程：

1、在工程文件夹db05下，首先建立tables文件夹的链接，命令如下。完成后一定要检查utl. pole. soltbl. luntab. nutabl.等文件，主要是处理的数据日期是否在这些文件有效范围内，否则到网站http://sopac.ucsd.edu/processing/gamit/ 下载
 sh_setup -yr 2009

2、在工程文件夹下建立各种工作必须目录，命令如下：
 mkdir rinex igs brdc

3、在rinex文件夹下，运行下面命令，下载4天的IGS的观测文件，主要是中国大陆的IGS站，共用7站，其中kunm、urum、lhaz3站作为已知站
 sh_get_rinex -archive sopac -yr 2009 -doy 121 -ndays 4 -sites kunm urum lhaz khaj wuhn shao bjfs guao
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/rinex/YYYY/DDD
ftpcmd ftp -inv garner.ucsd.edu
wgetsite garner.ucsd.edu
wlogin anonymous zhao@linux
########################################
directory = /pub/rinex/2009/121
Try 1: Getting rinex data for day 121 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1210.09d.Z urum1210.09d.Z lhaz1210.09d.Z khaj1210.09d.Z wuhn1210.09d.Z shao1210.09d.Z bjfs1210.09d.Z guao1210.09d.Z
directory = /pub/rinex/2009/122
Try 1: Getting rinex data for day 122 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1220.09d.Z urum1220.09d.Z lhaz1220.09d.Z khaj1220.09d.Z wuhn1220.09d.Z shao1220.09d.Z bjfs1220.09d.Z guao1220.09d.Z
directory = /pub/rinex/2009/123
Try 1: Getting rinex data for day 123 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1220.09d.Z kunm1230.09d.Z urum1230.09d.Z lhaz1230.09d.Z khaj1230.09d.Z wuhn1230.09d.Z shao1230.09d.Z bjfs1230.09d.Z guao1230.09d.Z
directory = /pub/rinex/2009/124
Try 1: Getting rinex data for day 124 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1240.09d.Z urum1240.09d.Z lhaz1240.09d.Z khaj1240.09d.Z wuhn1240.09d.Z shao1240.09d.Z bjfs1240.09d.Z guao1240.09d.Z


4、在SP3文件夹下运行如下命令，获取4天精密星历文件。
 sh_get_orbits -archive sopac -yr 2009 -doy 121 -ndays 4 -makeg no
ftp_info = sh_get_ftp_info -archive sopac -type sp3
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/products/GPSW
ftpcmd ftp -inv garner.ucsd.edu
########################################
Getting requested orbit files

5、在rinex文件夹下运行如下命令，获取4天的广播星历，这个星历包含所有卫星，文件以brdc开头，而不是测站出观测到的广播星历。
 sh_get_nav -archive sopac -yr 2009 -doy 121 -ndays 4
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/rinex/YYYY/DDD
ftpcmd ftp -inv garner.ucsd.edu
########################################
Getting auto or brdc files for 2009 121 for 4 days from garner.ucsd.edu

6、准备lfile.文件。这里分为两部分，已知约束点坐标（IGS）从tables目录下的itrf00.apr文件中拷贝，求解的未知点利用观测o文件手动生成（脚本如下），
制作批处理脚本genLfile(注意修改年份，这里是09o)放在rinex文件夹下，
#!/bin/bash
for file in *09o
do sh_rx2apr -site $file
done
cat *.apr >>./lfile.
rm *.apr

把上面文件保存在rinex文件夹下，执行方法如下：
（1）、修改可执行行：chmod +x genLfile
  执行：./genLfile
 或者利用如下命令：
（2）、sh genLfile

制作好的lfile. 如下替换掉tables文件夹中的相应文件：

 BJFS_GPS  -2148743.784   4426641.236   4044655.935  -0.0444   0.0141  -0.0013  1997.0  21601M001 BEIJING              .014  .019  .017  .0040 .0055 .0049
 GUAO_GPS   228380.790000     4631962.420000     4367041.080000   0.0000 0.0000 0.0000 2009.34
 KUNM_GPS  -1281255.473   5640746.079   2682880.117  -0.0318  -0.0024  -0.0203  1997.0  21609M001 KUNMING              .010  .022  .012  .0037 .0083 .0043
 LHAS_GPS   -106937.669   5549269.591   3139215.762  -0.0462  -0.0076   0.0121  1997.0  21613M001 LHASA                .001  .003  .002  .0005 .0011 .0007
 SHAO_GPS   -2831733.583000     4675665.958000     3275369.410000   0.0000 0.0000 0.0000 2009.34
 URUM_GPS    193030.873   4606851.324   4393311.421  -0.0306  -0.0055  -0.0003  1997.0  21612M001 URUMQI               .003  .006  .006  .0009 .0022 .0020
 WUHN_GPS  -2267749.162   5009154.325   3221290.762  -0.0325  -0.0077  -0.0119  1997.0  21602M001 WUHAN                .002  .003  .002  .0007 .0012 .0009

7、准备sites.defaults 如下：（kunm,urum,lhaz做为已知点）:xstinfo表示不更新station.info文件，glrepu表示重复解，glreps表示作为参考框架
   BJFS_GPS  db05 ftprnx           glrepu  
   GUAO_GPS  db05 ftprnx           glrepu  
   KUNM_GPS  db05 ftprnx  xstinfo glrepu glreps
   LHAS_GPS  db05 ftprnx  xstinfo glrepu glreps 
   SHAO_GPS  db05 ftprnx           glrepu  
   URUM_GPS  db05 ftprnx  xstinfo glrepu  glreps
   WUHN_GPS  db05 ftprnx           glrepu

8、制作station.info文件，利用在rinex文件夹下用命令手动生成，命令如下：
 sh_upd_stnfo -files *.09o -expt db05
最终如下：
 BJFS  BJFS              2009 121  0  0 30  2009 121 23 59 30   0.0460  DHARP   0.0000   0.0000  ASHZ12   9.20  ATDMCB   mstinf2: bjfs1210.09o
 BJFS  BJFS              2009 122  0  0 30  2009 122 23 59 30   0.0460  DHARP   0.0000   0.0000  ASHZ12   9.20  ATDMCB   mstinf2: bjfs1220.09o
 BJFS  BJFS              2009 123  0  0 30  2009 123 23 59 30   0.0460  DHARP   0.0000   0.0000  ASHZ12   9.20  ATDMCB   mstinf2: bjfs1230.09o
 BJFS  BJFS              2009 124  0  0 30  2009 124 23 59 30   0.0460  DHARP   0.0000   0.0000  ASHZ12   9.20  ATDMCB   mstinf2: bjfs1240.09o
.
 GUAO  GUAO              2009 121  0  0  0  2009 121 23 59 30   0.0800  DHARP   0.0000   0.0000  ATUZ12   9.93  ATDMRB   mstinf2: guao1210.09o
 GUAO  GUAO              2009 122  0  0  0  2009 122 23 59 30   0.0800  DHARP   0.0000   0.0000  ATUZ12   9.93  ATDMRB   mstinf2: guao1220.09o
 GUAO  GUAO              2009 123  0  0  0  2009 123 23 59 30   0.0800  DHARP   0.0000   0.0000  ATUZ12   9.93  ATDMRB   mstinf2: guao1230.09o
 GUAO  GUAO              2009 124  0  0  0  2009 124 23 59 30   0.0800  DHARP   0.0000   0.0000  ATUZ12   9.93  ATDMRB   mstinf2: guao1240.09o
.
 KUNM  KUNM              2009 121  0  0  0  2009 121 23 59 30   0.0793  DHARP   0.0000   0.0000  TR8000   0.00  TRBROG  mstinf2: kunm1210.09o?
 KUNM  KUNM              2009 122  0  0  0  2009 122 23 59 30   0.0793  DHARP   0.0000   0.0000  TR8000   0.00  TRBROG  mstinf2: kunm1220.09o?
 KUNM  KUNM              2009 123  0  0  0  2009 123 23 59 30   0.0793  DHARP   0.0000   0.0000  TR8000   0.00  TRBROG  mstinf2: kunm1230.09o?
 KUNM  KUNM              2009 124  0  0  0  2009 124 23 59 30   0.0793  DHARP   0.0000   0.0000  TR8000   0.00  TRBROG  mstinf2: kunm1240.09o?
.
 LHAZ  LHAZ              2009 121  0 38  0  2009 121 23 59 30   0.1330  DHARP   0.0000   0.0000  TRMSST   0.00  ATGGDB mstinf2: lhaz1210.09o??
 LHAZ  LHAZ              2009 122  0 38 30  2009 122 23 59 30   0.1330  DHARP   0.0000   0.0000  TRMSST   0.00  ATGGDB mstinf2: lhaz1220.09o??
.LHAZ  LHAZ              2009 123  0 38 30  2009 123 23 59 30   0.1330  DHARP   0.0000   0.0000  TRMSST   0.00  ATGGDB mstinf2: lhaz1230.09o??
 LHAZ  LHAZ              2009 124  0 38 30  2009 124 23 59 30   0.1330  DHARP   0.0000   0.0000  TRMSST   0.00  ATGGDB mstinf2: lhaz1240.09o??

 SHAO  SHAO              2009 121  0  0  0  2009 121 23 59 30   0.0814  DHARP   0.0000   0.0000  ATUZ12   0.00  TRBROG  mstinf2: shao1210.09o?
 SHAO  SHAO              2009 122  0  0  0  2009 122 23 59 30   0.0814  DHARP   0.0000   0.0000  ATUZ12   0.00  TRBROG  mstinf2: shao1220.09o?
 SHAO  SHAO              2009 123  0  0  0  2009 123 23 59 30   0.0814  DHARP   0.0000   0.0000  ATUZ12   0.00  TRBROG  mstinf2: shao1230.09o?
 SHAO  SHAO              2009 124  0  0  0  2009 124 23 59 30   0.0814  DHARP   0.0000   0.0000  ATUZ12   0.00  TRBROG  mstinf2: shao1240.09o?
.
 URUM  URUM              2009 121  0  0  0  2009 121 23 59 30   0.0460  DHARP   0.0000   0.0000  TRMSST   3.30  TPSC3D  mstinf2: urum1210.09o?
 URUM  URUM              2009 122  0  0  0  2009 122 23 59 30   0.0460  DHARP   0.0000   0.0000  TRMSST   3.30  TPSC3D  mstinf2: urum1220.09o?
 URUM  URUM              2009 123  0  0  0  2009 123 23 59 30   0.0460  DHARP   0.0000   0.0000  TRMSST   3.30  TPSC3D  mstinf2: urum1230.09o?
 URUM  URUM              2009 124  0  0  0  2009 124 23 59 30   0.0460  DHARP   0.0000   0.0000  TRMSST   3.30  TPSC3D  mstinf2: urum1240.09o?
.
 WUHN  WUHN              2009 121  0  0 30  2009 121 23 59 30   2.3610  DHARP  -0.0094  -0.0022  ASHZ12   9.20  ATDMCE   mstinf2: wuhn1210.09o
 WUHN  WUHN              2009 122  0  0 30  2009 122 23 59 30   2.3610  DHARP  -0.0094  -0.0022  ASHZ12   9.20  ATDMCE   mstinf2: wuhn1220.09o
 WUHN  WUHN              2009 123  0  0 30  2009 123 23 59 30   2.3610  DHARP  -0.0094  -0.0022  ASHZ12   9.20  ATDMCE   mstinf2: wuhn1230.09o
 WUHN  WUHN              2009 124  0  0 30  2009 124 23 59 30   2.3610  DHARP  -0.0094  -0.0022  ASHZ12   9.20  ATDMCE   mstinf2: wuhn1240.09o

 9、利用gamit批处理模式解算基线，在工程目录下执行如下命令：
sh_gamit -d 2009 121 122 123 124 -expt db05 -noftp > & ! sh_gamit_db05.log
可用gedit打开 sh_gamit_db05.log查看进行到了那里，是否出了错误等。
运行完毕后查看各年纪日文件夹下的sh_gamit_***.summary、qdb05a.***等文件，确定是否解算成功、精度符合规范等。此工程 Postfit nrms均在0.2附近，说明解算成功。
其中一个总结文件sh_gamit_121.summary如下：

 Input options -d 2009 121 122 123 124 -expt db05 -noftp

Processing 2009 121 GPS week 1529 5 Using node: linux-zhao Started at: 09_06_01_16:09:35

Processing directory: /home/zhao/project/db05/121
Disk Usage:  82428.7  Free  17969.6 Mbyte. Used 83%
 
Number of stations used 7 Total xfiles 7
X-files below minxf (300KB) threshold: 法访问 法访问 法访问 法访问 法访问 法访问 法访问
Sites excluded by xsite command
 
Postfit RMS rms, total and by satellite
RMS  IT Site   All  01  02  03  04  05  06  07  08  09  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32
RMS  16 ALL    8.0  69  66 113  69  60  84  85  76  73  66  88  89  78  87  66  82  70  68  77  95  74  82  94  73  77  69  73  78  76  82  88  90
Best and Worst two sites:
RMS  16 GUAO   6.9   7   5   8   7   6   8   6   8   6   6   8   9   8   6   6   6   6   5   6   8   7   7   8   6   5   6   7   6   8   6   7   8
RMS  16 SHAO   6.9   0   6   7   6   0   6   7   7   7   7   6   8   8   8   5   8   7   7   6   8   5   6   6   8   7   8   5   7   7   8   7   8
RMS  16 WUHN   8.6   0   6  10   8   0   8   8   9  10   8   7  10   9   9   7   9   8   8   8  11   7   9   9   8   7   8   8   8   9  10  11   9
RMS  16 KUNM  10.8   8   9  22   7   6   7  15   8   7   6  14   8   6  12   8  11   8   6  13  15   8   9  16  10  10   7  12  11   6  12  10  12

Double difference statistics
 Prefit nrms:  0.15731E+01    Postfit nrms: 0.19834E+00
 Prefit nrms:  0.15688E+01    Postfit nrms: 0.20417E+00
 Prefit nrms:  0.15731E+01    Postfit nrms: 0.18983E+00
 Prefit nrms:  0.15688E+01    Postfit nrms: 0.19448E+00
 Number of double differences:  27955

 Phase ambiguities (Total  WL-fixed   NL-fixed): 180 171 151  AUTCLN (Total Expected)  180 180
 Phase ambiguities WL fixed  95.0% NL fixed  83.9%

Processing 2009 121 GPS week 1529 5 Using node: linux-zhao Finished at: 09_06_01_16:14:55

基线解算就此完成。平差可以采用其他软件（如武大的PowerAdj）或采用gamit相应的globk。本工程采用globk平差。过程如下：

10、利用如下脚本生成要求的的文件夹链接，在工程目录下执行这个脚本，方法同前：
for file in 1??   # this is for current directory
   do
      ln -s $file 2009_$file
   done

11、运行如下命令，实现批处理平差模式：
 sh_glred -s 2009 121 2009 124 -yrext -local db05 -opt H G E
因为上面命令没有指定工程名，因此采用了默认的expr工程名。运行结束后查看gsoln文件夹下的sum.expr文件，详细结果在globk_expt_09***.org等4个文件中。
其中sum.expr 如下：

 ENSUM: GLOBK NEU component summary: includes the following solutions
 File tmp.ensum_input:  Combination of ALL networks
  Stations    enu  S#   #     mean length   sig   wrms  nrms     slope    sig    wrms nrms   dur   mean    line #
                                  (m)       (m)   (mm)         (mm/yr)  (mm/yr)  (mm)       (yrs)  (yrs)
SHAO_GPS-       N  1    4    3461996.3402 0.0355  18.3  0.26    1418.45 11634.00  21.5  0.30  0.01 2009.33       1-      7
SHAO_GPS-       E  1    4   11552738.9544 1.2501 225.6  0.09   19993.08 ******** 262.4  0.10  0.01 2009.33      12-     18
SHAO_GPS-       U  1    4         21.5172 0.4081 202.6  0.25   ******** ******** 238.5  0.29  0.01 2009.33      23-     29
BJFS_GPS-       N  1    4    4409209.3915 0.1631  52.8  0.16    5482.05 53351.20  60.1  0.18  0.01 2009.33      34-     40
BJFS_GPS-       E  1    4    9939262.4549 1.2346 195.4  0.08   18376.46 ******** 225.8  0.09  0.01 2009.33      45-     51
BJFS_GPS-       U  1    4         86.9263 0.4167 169.2  0.20   ******** ******** 189.5  0.23  0.01 2009.33      56-     62
WUHN_GPS-       N  1    4    3398768.0109 0.1700  23.1  0.07    2624.29 55538.72  25.9  0.08  0.01 2009.33      67-     73
WUHN_GPS-       E  1    4   10965014.4971 1.2472 198.7  0.08   18470.75 ******** 229.8  0.09  0.01 2009.33      78-     84
WUHN_GPS-       U  1    4         25.3512 0.4111 235.7  0.29   ******** ******** 276.5  0.34  0.01 2009.33      89-     95
KUNM_GPS-       N  1    4    2786275.3981 0.4423  17.2  0.02    3132.77 ********  16.1  0.02  0.01 2009.33     100-    106
KUNM_GPS-       E  1    4   10368693.2207 1.2557 144.6  0.06   14968.96 ******** 164.8  0.07  0.01 2009.33     111-    117
KUNM_GPS-       U  1    4       1985.7652 0.4130 310.3  0.38   ******** ******** 367.3  0.44  0.01 2009.33     122-    128
LHAZ_GPS-       N  1    4    3301439.1381 0.7038  59.8  0.04    5547.11 ********  69.2  0.05  0.01 2009.33     133-    139
LHAZ_GPS-       E  1    4    8813208.7511 1.2231  86.0  0.04    8353.36 ********  99.0  0.04  0.01 2009.33     144-    150
LHAZ_GPS-       U  1    4       3624.1771 0.4132 305.0  0.37   ******** ******** 353.8  0.43  0.01 2009.33     155-    161
URUM_GPS-       N  1    4    4876678.6982 0.7840 124.8  0.08   12312.75 ******** 143.2  0.09  0.01 2009.33     166-    172
URUM_GPS-       E  1    4    7037402.0994 1.1504  98.7  0.04    4520.03 ******** 119.3  0.05  0.01 2009.33     177-    183
URUM_GPS-       U  1    4        858.4000 0.4196 196.7  0.23   ******** ******** 206.9  0.25  0.01 2009.33     188-    194
GUAO_GPS-       N  1    4    4839181.6651 0.7923 124.1  0.08   12314.34 ******** 142.3  0.09  0.01 2009.33     199-    205
GUAO_GPS-       E  1    4    7042905.1263 1.1501  96.0  0.04    4482.59 ******** 115.9  0.05  0.01 2009.33     210-    216
GUAO_GPS-       U  1    4       2028.2399 0.4193 198.1  0.24   ******** ******** 207.0  0.25  0.01 2009.33     221-    227

对比各站的IGS坐标，可见符合要求！

至此整个工程的解算过程就完成了，当然这里只是示范性的解算，有很多参数可以在process.defaults、sestbl.、sittbl.中具体设置，而不采用默认值。
gamit\globk软件虽然采用可以批处理模式，使得相对简化了处理的复杂性，但只用在不断的反复解算中才能掌握各种处理的技巧、问题解决的方法！

使用 Scrapbook作为个人资料管理系统

我现在用Scrapbook（一 个Firefox的插件）作为我的个人资料管理系统，顺便还存储一点个人效率系统的资料。这里的资料主要是指每天获得的零散的信息（互联网上的阅读等）以 及自己的随手笔记。关于电子书等其他资料和文件系统的管理，以后我会专门撰文叙述。
GTD中提到建立一个简单而高效的个人资料管理系统是具有重要作用的。实际上，我看完GTD这本书就觉得资料管理系统实际上跟个人效率系统是两套不 同的体系。个人效率（时间管理）系统侧重于让我们更有效率的完成工作。而资料管理系统则是为了“有效组织和管理资料，做到随时需要，随时获得”（古道语，虽然我觉得古道把 事情搞得有点麻烦^^）。
于是当时Google了一番，发现善用佳软上有一个系 列“寻找最好用的笔 记软件”，比较的很全面。于是按图索骥，对十来个软件进行了试用，最后选了MyBase作为我的资料管理系统，并花了一两天时间把 以前存在Scrapbook里的零散资料转到了MyBase里。在Windows下，MyBase是不错的笔记软件了，树形组织十分方便，还提供了IE和 Firefox的抓取插件。除了启动慢一点，Firefox抓网页乱码几率很高，搜索基本上没用以外，已经算比较完美了。我也用的很happy。
但是几个月前，我实在受不了越来越慢的Windows，把整个硬盘格了装成了Ubuntu，这下问题来了，MyBase没有提供任何Linux的支持。 试用了Ubuntu下的几个笔记软件，比如Basket、Jarnal什 么的，也没发现顺手的工具。正巧那几天发现了Scrapbook的几个插件，感觉比较兴奋，于是重新启用Scrapbook。
Scrapbook是Firefox的一个插件，用来保存网页以便离线的时候也可以查阅。目前的版本是1.3.3.2，已经基本上可以胜任一个笔记软件的需求。再加上几个插件，可以让你用的很 High。另外，很重要的一点是，Scrapbook是Firefox的插件，跨平台，而且它所使用的技术、代码、数据格式等都是开放的（Mybase的 数据导出可折腾死我了）。在你觉得某个功能不爽的时候，你甚至可以自己动手来完善它LoL
首先来看一下需要安装的插件，如下图。首先Scrapbook本身，然后是Integrated Search和AutoSave。 安装完Scrapbook之后，就可以整合到Firefox的右键菜单里面。具体的设置各位可以自己琢磨。当然，也有热心用户为Scrapbook专门写 了用 户指南。

安装完之后的Scrapbook是什么样子呢？

所有的文件夹都是自己建立，你想怎么组织你的知识管理系统都可以，关键就是要用起来舒服。比如说我就不喜欢清单系统里面的东西太多，于是把每周要回 顾的内容和的清单放到了这里面，每周进行回顾和改进。
图的最上面那个输入框是Integrated Search，它可以让你搜索你的整个Scrapbook资料系统，而且速度很快，这比起Google Desktop这种笨东西来，要好的多了。
分割线是我自己加的，下面是AutoSave保存的网页。因为每天看的东西很多，而我又不是随时记得把网页保存到资料系统里去。AutoSave就 是把每一个我看过的页面保存下来，够暴力吧？但是“一力破十慧”，它还就是好用，我在看东西的时候再也不担心会丢掉有用的内容了。而且没有网络连接的时 候，我想起某个看过的网页也可以拿出来翻翻。什么？担心占的磁盘越来越大？我们这不是在实践GTD么？把“清理资料系统”加入到你的每周回顾里面不就得 了？我已经这么干了四、五周，基本上每次清理的时间不会超过15分钟吧。
需要注意的是，Windows用户的Scrapbook的默认目录是在Application Data下面的mozilla目录，大家还是把它移到数据集中存放的地点吧（Alt+K打开Scrapbook面板，工具->组织）。要不然如果某 次重装忘了备份，就只有泪流满面了。当然，也可以使用上面图中的BackupHelper。不过我在Ubuntu下，主文件夹还是蛮安全的。
好啦，总结一下，Firefox+Scrapbook+Integrated Search+AutoSave＝轻便好用一站式个人资料系统，自由开放跨平台，而且将会越来越快（因为Firefox 3会有大的性能提升）。如果你担心你的资料管理系统某天崩掉或者某天你心血来潮要切换到另一个平台，那么试一下Scrapbook吧，应该不会让你失望。
你也使用Scrapbook么？如果你有什么独门秘技，不要藏私哦，欢迎留言与我们分享。
http://www.productivelife.cn/2008/04/personal-km-system.html

坐标转换程序（xyz－>geo）||convertc(gamit附带实用程序)

zhao@zhao-laptop:~$ convertc

 CONVERTC: Converts coordinate types and files
 
 Runstring:
 
  convertc
       OR  
  convertc
       OR
  converc
 
  where
 
  is the name of the input file.
   Type determined by convertc; supported:
    GLOBK apr file (new GAMIT l-file), XYZ
    GLOBK velocity file (lon/lat ENU velocity
    GLOBK glist file (lon/lat)  
  X Y Z (m) or lat lon ht (deg, m)
  Xdot Ydot Zdot OR Ndot Edot Udot (m/yr)
  is the name of the output file
  is the type of output file or coords:
    LFILE GAMIT oldstyle, spherical deg/min/sec
    APR   GLOBK apr file (new GAMIT l-file), XYZ
    VEL   GLOBK velocity file (lon/lat ENU velocity
    KML   Google Earth input (lat/lon/ht)
    GEO   Geodetic coords (lat/lon/ht), full precision
  If coords input on command-line, output is to screen
 
  Examples:
   convertc emed.apr emed.kml KML
   convertc 918129.451 -4346071.255  4561977.839 GEO
   convertc  43.36 -112.567 27.2  XYZ
   convertc  43.36 247.433  27.2  LFILE
zhao@zhao-laptop:~$ convertc 388042.6123 -740382.4073 6302001.9040  geo
Invalid screen output type geo
zhao@zhao-laptop:~$ convertc 388042.6123 -740382.4073 6302001.9040  GEO
   82.49429  -62.34047   78.2
zhao@zhao-laptop:~$

统计软件教程---主成分分析

统计软件教程

这是北京大学概率统计系《统计软件》课程的讲义， 根据MS Word版本转换而成。现有七章，以后还可能 有新内容。

欢迎批评与建议。

作者 李东风
Email: ldf@math.pku.edu.cn

目录

• 前 言
• SAS 初阶
• SAS 语言与数据管理
• SAS 过程初步
• SAS 的基本统计分析功能
• SAS 多元统计分析
• S 语言介绍
• Gauss 系统介绍

前言

统计学是研究如何收集数据、分析数据并进行推断的学科。 统计学的应用必然要涉及数据的收集、存贮、整理， 以及各种统计方法的实现， 这些都要靠统计软件的帮助来完成。

从电子计算机出现至今，统计软件已经有了长足的发展。 一方面经典的统计方法都已被实现到统计软件中， 另一方面，帮助统计学家实现新的统计方法的软件也极大地 推动了新的统计计算方法的研究与开发。 国内统计软件的应用这些年来也有很大进步， 引进了国外的SAS、SPSS、Splus、Gauss等统计软件并出版了 一些有关书籍。然而，现在介绍统计软件的书多数是着重讲软件本身的功能， 没有考虑到初学者循序渐进的要求， 可供自学的教材性质的书比较少， 所以本书试图填补这一空白，以自学教材的形式讲述， 不强调全面介绍软件功能，而是让读者能够通过阅读本书掌握统计软件的基本用法， 然后可以进一步阅读更详细的资料以达到熟练使用统计软件完成统计应用的目的。

本书主要讲SAS系统和Splus系统。SAS系统是一个集大型数据库管理、 统计分析、报表图形、信息系统开发等多种强大功能为一体的大型软件系统， 是国际公认的优秀统计分析软件，但初学有些困难， 本书力图采取循序渐进的讲法使读者轻易掌握SAS的要点。 Splus是一种面向对象的统计编程与统计分析系统，它具有强大的探索性数据分析功能， 实现了很多最新的统计方法，而且用Splus编制统计计算程序既容易学会又容易使用， 可以极大地降低统计计算程序的开发时间，方便了统计学家研究新的统计方法和进行统计计算。

本书可用作大学统计软件课程的教材，也可以由需要使用统计软件的读者作为自学教程。 阅读本书需要读者有基本的数理统计知识， 另外，读者最好已经学过计算机程序设计， 本书讲的是专用语言，所以不适合作为一般的程序设计入门教材。 本书的着重点是用统计软件进行数据管理、统计分析和统计计算编程， 适用于统计软件的初学者，对于用统计软件开发商业化应用程序所需知识涉及较少， 希望进行统计应用开发的读者可以在阅读本书后进一步学习有关知识。