保存当前进度
下次运行可以从上次进度运行
可以更新县区代码
可以更新时间

最终的代码如下：

#!/bin/sh
# cd /opt/to_reback/data/landuse
#set -x
d1=./tmpdata1
d2=./tmpdata2
YEAR="80年代末 1995年 2000年"
for a in $YEAR; do
for b in $a/*; do #province, 目录判断
if [ -d $b ]; then
rm -fdr $d1
rm -fdr $d2
for c in $b/ld*; do # areacode, 县级区域
bb=`echo $b | sed "s/\///g"` # 去除/特殊字符，否则sed罢工
areacode=`echo $c | sed "s/\///g" | sed "s/$bb//g" | sed "s/ld//g"`
#echo $c
if [ `cat work.log | grep "$c"` ]; then
echo $c
echo "already done!"
else
ogr2ogr $d1 $c && ogr2ogr $d2 $d1/PAL.shp -f "ESRI Shapefile" -select area,perimeter,r_area,ln71,ln72,ln73,ln74,aftln,tm20,tm30,tm40,tm50,tm60,afttm,ph41,ph52,ph53,net || exit
shp2pgsql -a -g the_geom -i $d2/PAL.shp landuse | psql -d westdc || exit
psql -d westdc -c "update landuse set areacode='$areacode' where areacode=''" || exit
echo $c >> work.log
# echo "here not done"
echo $c
#exit
fi
done
#vacuumdb westdc
fi
done
psql -d westdc -c "update landuse set year='$a' where year=''" || exit
done
#set +x

#!/bin/sh
# cd /opt/to_reback/data/landuse
tmpdata=./tmpdata
YEAR="80年代末 1995年 2000年"
for a in $YEAR; do
for b in $a/*; do #province, 目录判断
if [ -d $b ]; then
rm -fdr $tmpdata
for c in $b/ld*; do # areacode, 县级区域
bb=`echo $b | sed "s/\///g"` # 去除/特殊字符，否则sed罢工
areacode=`echo $c | sed "s/\///g" | sed "s/$bb//g" | sed "s/ld//g"` #主要时间都花费在这儿了
ogr2ogr $tmpdata $c
ogr2ogr $tmpdata/poly.shp $tmpdata/PAL.shp -select area,perimeter,r_area,ln71,ln72,ln73,ln74,aftln,tm20,tm30,tm40,tm50,tm60,afttm,ph41,ph52,ph53,net
shp2pgsql -a -g the_geom -i $tmpdata/poly.shp landuse | psql -d westdc
psql -d westdc -c "update landuse set areacode='$areacode' where areacode=''"
psql -d westdc -c "update landuse set year='$a' where year=''"
# exit
done
fi
done
done

最后再添加一个超级用户：
create user xxx with password 'passwd' createuser createdb

pg_dump -f dumpfile.sql yourdatabase
dropdb yourdatabase
createdb yourdatabase
createlang plpgsql yourdatabase
psql -f lwpostgis.sql -d yourdatabase
psql -f dumpfile.sql -d yourdatabase
vacuumdb -z yourdatabase

（这样好像有点安全问题，但若本机只有自己能用，安全性也还过得去吧）

1、安装

安装很简单，首先要安装postgresql，注意其有两个帐号，一是系统帐号，二是数据库管理帐号。安装程序可以新建用户，要提到一点，目前的postgresql数据库不支持远程安装，就是说，不能通过远程桌面进行安装（这个限制，不爽）。windows版本还包括了pgadmin3，图形化的管理界面。
然后安装postgis，需要前面安装postgresql的管理员帐号，同时其还要安装一个postgis数据库。
默认情况下，可执行文件都安装到postgresql的bin目录，几个SQL文件安装在postgresql的share/contrib目录下。

2、使用

对于中文用户来讲，首先要考虑数据库的编码问题，可以使用unicode,utf8,euc_cn等编码，默认的是unicode编码。
第一次使用unicode编码时移植linux下的应用，遇到了很多非常奇怪的问题。改用euc_cn编码，问题就少多了。

建库：
createdb -U postgres -E EUC_CN db
（linux环境下本机默认都不需要认证密码，但在windows系统下，默认是需要验证密码的）
alter database db set client_encode="EUC_CN"
（这个要进入到psql环境中才能运行的，不确定是否必须，在linux环境下不需要）
搭建postgis环境：
creatlang -U postgres plpgsql db
(这个在windows下默认好象是不需要处理的）
引入postgis支持：
psql -U postgres -d db -f lwpostgis.sql
psql -U postgres -d db -f spatial_ref_sys.sql
生成矢量数据的sql文件：
shp2pgsql xx.shp xx_table > xx.sql
导入到数据库：
psql -U postgres -d db -f xx.sql

3、开始应用

到此，您的基本的postgis数据库就已经搭建好了，剩下的工作就是看您怎么去最大的压榨postgis的功能了。