BlogJava-zhyiwww -随笔分类-gis

添加空间字段的另一种方法

zhyiwww — Thu, 22 Oct 2009 07:46:00 GMT

zydb=# select addGeometryColumn('public','pbeijing','the_geom_900913',900913,'POINT',2);
                       addgeometrycolumn
---------------------------------------------------------------
public.pbeijing.the_geom_900913 SRID:900913 TYPE:POINT DIMS:2

(1 row)

zydb=# \d pbeijing
                                      Table "public.pbeijing"
     Column      |          Type          |                       Modifiers
-----------------+------------------------+--------------------------------------------------------
gid             | integer                | not null default nextval('pbeijing_gid_seq'::regclass)
mapid           | character varying(8)   |
kind            | character varying(4)   |
name            | character varying(100) |
py              | character varying(250) |
exonym          | character varying(250) |
adminname       | character varying(50) |
dzname          | character varying(50) |
zipcode         | character varying(6)   |
address         | character varying(100) |
telephone       | character varying(15) |
admincode       | character varying(6)   |
dzcode          | character varying(6)   |
isindex         | boolean                |
isother         | boolean                |
the_geom        | geometry               |
x               | double precision       |
y               | double precision       |
the_geom_900913 | geometry               |
Indexes:
    "pbeijing_pkey" PRIMARY KEY, btree (gid)
    "pbeijing_the_geom_gist" gist (the_geom)
Check constraints:
    "enforce_dims_the_geom" CHECK (ndims(the_geom) = 2)
    "enforce_dims_the_geom_900913" CHECK (ndims(the_geom_900913) = 2)
    "enforce_geotype_the_geom" CHECK (geometrytype(the_geom) = 'POINT'::text OR the_geom IS NULL)
    "enforce_geotype_the_geom_900913" CHECK (geometrytype(the_geom_900913) = 'POINT'::text OR the_geom_900913 IS NULL)
    "enforce_srid_the_geom" CHECK (srid(the_geom) = 4326)
    "enforce_srid_the_geom_900913" CHECK (srid(the_geom_900913) = 900913)

zydb=# update pbeijing set the_geom_900913=transform(the_geom,900913);
UPDATE 34159

zydb=# select the_geom ,the_geom_900913 from pbeijing limit 10;
                      the_geom                      |                  the_geom_900913
----------------------------------------------------+----------------------------------------------------
0101000020E6100000927A4FE5B41F5D408D2782380FCF4340 | 010100002031BF0D00B4A0714C22BC684182BB38C5C9595241
0101000020E610000098A59D9ACB125D409949D40B3EED4340 | 010100002031BF0D00FBE9FE1C2BB16841EF87D46C1F7B5241
0101000020E6100000B857E6ADBA085D404AF086342ADE4340 | 010100002031BF0D009FC3C3879EA868417558A1E4746A5241
0101000020E6100000041DAD6A49135D404A09C1AA7AF94340 | 010100002031BF0D0075C44FF795B16841D7D4B67FAB885241
0101000020E610000017838769DF175D4083BF5FCC962E4440 | 010100002031BF0D00E581730A7BB5684137D06310B0C35241
0101000020E6100000ABAE4335251C5D401E705D3123F24340 | 010100002031BF0D004A7A81001CB96841E7D1D22E8A805241
0101000020E6100000E083D72E6D135D409CDB847B650A4440 | 010100002031BF0D00575CA557B4B168419AA92F116E9B5241
0101000020E610000064B2B8FFC8165D40C7D8092FC1F14340 | 010100002031BF0D004B7295958EB46841215637B21D805241
0101000020E6100000F60A0BEE071B5D40F7729F1C05EE4340 | 010100002031BF0D009F0D11B729B86841BFAFDEA9FB7B5241
0101000020E6100000BD8DCD8E54155D40A6EECA2E18F44340 | 010100002031BF0D0029501D4552B36841F88C1DCEB4825241
(10 rows)

zydb=#

zhyiwww 2009-10-22 15:46 发表评论

给空间数据表添加空间字段并实现投影转换

zhyiwww — Wed, 21 Oct 2009 08:45:00 GMT

我有一个表rbeijing,有空间字段the_geom,但是是基于SRS4326的，我想用900913d的投影，所以就学要添加一个新的空间字段并转换投影坐标。做法如下：

mapdb=# alter table rbeijing add column the_geom_900913 geometry;
mapdb=# create index index_rbeijing_the_geom_900913 on rbeijing using gist(the_geom_900913);
mapdb=# create index index_rbeijing_the_geom on rbeijing gist(the_geom);
mapdb=# update table rbeijing set the_geom_900913=st_tranform(the_geom,900913);

mapdb=# update rbeijing set the_geom_900913 = st_transform(the_geom,900913);

mapdb=# \d rbeijing
                                      Table "public.rbeijing"
     Column      |          Type          |                       Modifiers
-----------------+------------------------+--------------------------------------------------------
gid             | integer                | not null default nextval('rbeijing_gid_seq'::regclass)
mapid           | character varying(8)   |
id              | character varying(13) |
kind_num        | character varying(2)   |
kind            | character varying(23) |
width           | character varying(3)   |
direction       | character varying(150) |
const_st        | character varying(1)   |
pathname        | character varying(40) |
pathpy          | character varying(250) |
snodeid         | character varying(13) |
enodeid         | character varying(13) |
pathclass       | character varying(2)   |
pathno          | character varying(13) |
the_geom        | geometry               |
the_geom_900913 | geometry               |
Indexes:
    "rbeijing_pkey" PRIMARY KEY, btree (gid)
    "index_rbeijing_the_geom" gist (the_geom)
    "index_rbeijing_the_geom_900913" gist (the_geom_900913)
Check constraints:
    "enforce_dims_the_geom" CHECK (ndims(the_geom) = 2)
    "enforce_geotype_the_geom" CHECK (geometrytype(the_geom) = 'MULTILINESTRING'::text OR the_geom IS NULL)
    "enforce_srid_the_geom" CHECK (srid(the_geom) = 4326)

上面的表结构是添加和转换后的。添加空间字段的同时，也为空间字段添加了索引。

操作前提：postgis/postgres能支持900913，如果不支持的话，就学要先添加900913投影支持。

zhyiwww 2009-10-21 16:45 发表评论

mapserver中使用mapfile的调试日志

zhyiwww — Wed, 02 Sep 2009 09:15:00 GMT

在调试mapfile的时候，往往不知道怎么才能看到mapserver的执行细节，所以，能查看日至不失为一个好的办法。
在mapserver中使用日志的方法如下:

(1)在MAP中添加输出日至文件

CONFIG "MS_ERRORFILE" "/tmp/map_error.log"

(2)在LAYER里面打开DEBUG
也就是设置DEBUG的级别。

DEBUG 5
根据mapserver的参考文档，可以设置ON|OFF|0|1|2|3|4|5.
调试阶段，设置成5，可以得到最多的输出信息。

(3)实时查看日志
tail -f /tmp/map_error.log

(4)错误调试
比如，我添加一个错误的标签AAAA,系统输出日志如下：
[Wed Sep 2 16:59:58 2009].311416 msPOSTGISLayerClose -- layerinfo is NULL
[Wed Sep 2 17:12:25 2009].871161 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].888432 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].902453 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].907025 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].942961 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].972666 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].977157 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].998068 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:25 2009].998070 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].19756 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].43084 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].51924 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].68473 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].73118 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
[Wed Sep 2 17:12:26 2009].82840 loadLayer(): Unknown identifier. Parsing error near (AAAA):(line 149)
。。。。。。。。。。。。。

这样，你就可以根据错误信息很容易的找到错误并更正错误。

zhyiwww 2009-09-02 17:15 发表评论

mapserver-给道路添加label

zhyiwww — Tue, 25 Aug 2009 11:54:00 GMT

Mapfile中,需要使用如下配置:

(1)在LAYER中添加LABEL所用的数据字段

LABELITEM "pathname"

(2)在MAP中添加字体引用

FONTSET "fonts/fonts-lite.list"

(3)在CLASS中添加LABEL设置

LABEL

FONT "msyh"
ENCODING GB2312
TYPE TRUETYPE
COLOR 250 20 20
SIZE 8
POSITION uc
ANGLE 25
PARTIALS FALSE
BUFFER 60
END

效果如下图:

zhyiwww 2009-08-25 19:54 发表评论

mapfile中使用FILTER实现道路分级

zhyiwww — Mon, 24 Aug 2009 04:01:00 GMT

同上一篇mapfile中使用子查询实现道路分级。
在mapfile中，也可以用FILTER来数据的过滤，通过此方法也可以实现道路的分级。

DATA "the_geom from rbeijing as subquery using unique gid USING SRID=4326"
FILTER "pathclass='3'"

pathclass是表rbeijing中的字段的名称。3是其值，表示道路的一种级别。

zhyiwww 2009-08-24 12:01 发表评论

ogr2ogr

zhyiwww — Thu, 20 Aug 2009 03:11:00 GMT

在linux下面，默认的是没有ogr2ogr命令的，如果想用，就要自己安装。

在ubuntu下面，可以用
aptitude install gdal-bin -y
安装完成就可以用ogr2ogr命令了。

zhyiwww 2009-08-20 11:11 发表评论

mapfile中使用子查询实现道路分级

zhyiwww — Wed, 19 Aug 2009 09:51:00 GMT

查看道路分级如下：

mapdb=# select distinct pathclass from rbeijing;
pathclass
-----------
0
15
2
3
4
(5 rows)

使用子查询实现道路分级

CONNECTIONTYPE postgis
PROCESSING "CLOST_CONNECTION=DEFER"
CONNECTION "user=map password=postgres dbname=mapdb host=192.168.1.11 port=5432"
DATA "the_geom from (select * from rbeijing where pathclass='3') as subquery using unique gid USING SRID=4326"

分级前的查询代码
CONNECTIONTYPE postgis
PROCESSING "CLOST_CONNECTION=DEFER"
CONNECTION "user=zhangyi password=postgres dbname=zydb host=192.168.10.166 port=5432"
DATA "the_geom from rbeijing using unique gid USING SRID=4326"
#DATA "the_geom from (select * from rbeijing where pathclass='3') as subquery using unique gid USING SRID=4326"

分级前的效果

zhyiwww 2009-08-19 17:51 发表评论

让TileCache支持AGG-反锯齿

zhyiwww — Tue, 18 Aug 2009 09:22:00 GMT

(1)主要在mapscript(beijing.map)文件里面添加如下配置：
        OUTPUTFORMAT
                NAME AGG
                DRIVER AGG/PNG
                IMAGEMODE RGB
                MIMETYPE IMAGE/PNG
                EXTENSION PNG
                FORMATOPTION TRANSPARENT=ON
                FORMATOPTION INTERLANCE=OFF
        END

(2) 在OpenLayer里面使用AGG
           var serverURL = "/cgi-bin/tilecache.cgi";

            var layer_obj = new OpenLayers.Layer.WMS(
                "beijing",
                serverURL,
                {
                    layers: 'beijing',
                    map: '/wwwroot/shp/bj.map',
                    format: 'AGG'
                },
                {
                    reproject: false,
                    numZoomLevels: 20
                }
            );

效果，请参看上篇文章

zhyiwww 2009-08-18 17:22 发表评论

让你的mapserver支持AGG

zhyiwww — Tue, 18 Aug 2009 07:07:00 GMT

系统平台，在ubuntu上面
(1)aptitude install mapserver libagg-dev

(2)在map文件里面添加输出AGG支持

OUTPUTFORMAT

NAME'AGG'

DRIVERAGG/PNG

IMAGEMODERGB

END

(3)在OpenLayers里面使用AGG

        var layer_obj = new OpenLayers.Layer.WMS(
                "beijing",
                serverURL,
                {
                    layers: 'rbeijing',
                    map: '/wwwroot/shp/bj.map',
                    format: 'AGG'
                },
                {
                    reproject: false,
                    numZoomLevels: 20,
                    gutter: 15,
                    buffer: 0
                }
            );

前后的效果比较
没有使用AGG支持的效果

使用AGG支持的效果

zhyiwww 2009-08-18 15:07 发表评论

shp2pgsql的SRID参数

zhyiwww — Fri, 10 Jul 2009 08:31:00 GMT

#shp2pgsql Rbeijing Rbeijing > rbeijing.sql
导出文件如下：
     1 BEGIN;
     2 CREATE TABLE "rbeijing" (gid serial PRIMARY KEY,
     3 "mapid" varchar(8),
     4 "id" varchar(13),
     5 "kind_num" varchar(2),
     6 "kind" varchar(23),
     7 "width" varchar(3),
     8 "direction" varchar(150),
     9 "const_st" varchar(1),
    10 "pathname" varchar(40),
    11 "pathpy" varchar(250),
    12 "snodeid" varchar(13),
    13 "enodeid" varchar(13),
    14 "pathclass" varchar(2),
    15 "pathno" varchar(13));
    16 SELECT AddGeometryColumn('','rbeijing','the_geom','-1','MULTILINESTRING',2);

7 INSERT INTO "rbeijing" ("mapid","id","kind_num","kind","width","direction","const_st","pathname","pathpy","snodeid","enodeid","pathclass","pathno",th       e_geom) VALUES ('595525','59552500001','1','0401','55','0','1',NULL,NULL,'59552500003','59552500136','2','0','010500000001000000010200000007000000BD5       296218EED5C40D844662E70C7434069C4CC3E8FED5C40103FFF3D78C743405DFC6D4F90ED5C40897956D28AC743405DFC6D4F90ED5C405D6DC5FEB2C74340A374E95F92ED5C40A587A1D5       C9C743400806103E94ED5C40AB9509BFD4C74340DCD6169E97ED5C405CAD1397E3C74340');

#shp2pgsql -s 4326 Rbeijing Rbeijing > rbeijing.sql
导出文件如下：
     1 BEGIN;
     2 CREATE TABLE "rbeijing" (gid serial PRIMARY KEY,
     3 "mapid" varchar(8),
     4 "id" varchar(13),
     5 "kind_num" varchar(2),
     6 "kind" varchar(23),
     7 "width" varchar(3),
     8 "direction" varchar(150),
     9 "const_st" varchar(1),
    10 "pathname" varchar(40),
    11 "pathpy" varchar(250),
    12 "snodeid" varchar(13),
    13 "enodeid" varchar(13),
    14 "pathclass" varchar(2),
    15 "pathno" varchar(13));
    16 SELECT AddGeometryColumn('','rbeijing','the_geom','4326','MULTILINESTRING',2);

INSERT INTO "rbeijing" ("mapid","id","kind_num","kind","width","direction","const_st","pathname","pathpy","snodeid","enodeid","pathclass","pathno",the       _geom) VALUES ('595525','59552500001','1','0401','55','0','1',NULL,NULL,'59552500003','59552500136','2','0','SRID=4326;0105000000010000000102000000070       00000BD5296218EED5C40D844662E70C7434069C4CC3E8FED5C40103FFF3D78C743405DFC6D4F90ED5C40897956D28AC743405DFC6D4F90ED5C405D6DC5FEB2C74340A374E95F92ED5C40A       587A1D5C9C743400806103E94ED5C40AB9509BFD4C74340DCD6169E97ED5C405CAD1397E3C74340');

粗体部分是最主要的区别点。geom字段的值会因为SRID的不同而不同。

zhyiwww 2009-07-10 16:31 发表评论

让postgis支持900913投影

zhyiwww — Thu, 09 Jul 2009 09:57:00 GMT

# insert into spatial_ref_sys VALUES (900913, 'EPSG', 900913, '', '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs');


gisdb=# select astext(transform(geomfromtext('POINT(116.78787 39.29983)',4326),900913));
                 astext                  
-----------------------------------------
 POINT(13000766.219231 4764711.10113251)
(1 row)

zydb=# select astext(transform(geomfromtext('POINT(116.78787 39.29983)',4326),32650));
                  astext                  
------------------------------------------
 POINT(481709.059137956 4350071.34479768)
(1 row)

gisdb=#

zhyiwww 2009-07-09 17:57 发表评论

tilecache一个异常问题

zhyiwww — Sun, 28 Jun 2009 15:52:00 GMT

我的错误日志如下：
root@XXX:/etc/rc2.d# tilecache_http_server
Listening on port 8080
Traceback (most recent call last):
File "/usr/lib/python2.6/wsgiref/handlers.py", line 93, in run
    self.result = application(self.environ, self.start_response)
File "/var/lib/python-support/python2.6/TileCache/Service.py", line 328, in wsgiApp
    return wsgiHandler(environ, start_response, theService)
File "/var/lib/python-support/python2.6/TileCache/Service.py", line 240, in wsgiHandler
    from paste.request import parse_formvars
ImportError: No module named paste.request
localhost - - [28/Jun/2009 23:47:49] "GET /tilecache/tilecache.py/1.0.0/NASA_JPL_WMS/1/3/1.png HTTP/1.1" 500 59
Traceback (most recent call last):
File "/usr/lib/python2.6/wsgiref/handlers.py", line 93, in run
    self.result = application(self.environ, self.start_response)
File "/var/lib/python-support/python2.6/TileCache/Service.py", line 328, in wsgiApp
    return wsgiHandler(environ, start_response, theService)
File "/var/lib/python-support/python2.6/TileCache/Service.py", line 240, in wsgiHandler
    from paste.request import parse_formvars
ImportError: No module named paste.request

错误的关键就是上面的黑体部分。
解决的方法：
# aptitude install python-paste -y
就可以解决问题了。

zhyiwww 2009-06-28 23:52 发表评论

Postgis的安装和配置

zhyiwww — Fri, 26 Jun 2009 08:55:00 GMT

把下面的代码存为 install_postgis.sh
然后执行。
必须在postgres用户下执行。

######################################################################
## install postgis on postgresql 8.3
##
## NOTICE:
## su to postgres user,then execute sh install_postgis.sh
##
######################################################################

dropdb postgis_template

createdb -T template1 postgis_template

createlang -d postgis_template plpgsql

cd /usr/share/postgresql-8.3-postgis

psql -d postgis_template -f lwpostgis.sql
psql -d postgis_template -f spatial_ref_sys.sql

psql -d postgis_template -c "GRANT ALL ON geometry_columns TO PUBLIC"
psql -d postgis_template -c "GRANT ALL ON spatial_ref_sys TO PUBLIC"

psql -d postgis_template -c "VACUUM FREEZE"

psql -l

zhyiwww 2009-06-26 16:55 发表评论

debian下安装mapserver的简单方法

zhyiwww — Fri, 03 Apr 2009 06:53:00 GMT

debian 下安装mapserv5.2.2的简单方法
[1]安装gd
参考
[2]安装所需的包
aptitude install gdal-bin libagg-dev libgeos-dev libgdal1-dev libxml2-dev -y

[3]安装mapserver

./configure --enable-debug --with-freetype --with-zlib --with-png --with-jpeg --with-xpm --with-gd=/opt/gd --with-pdf --with-agg --with-eppl --with-proj --with-threads --with-geos --with-ogr --with-gdal --with-tiff --with-postgis --with-ming --with-wfs --with-wcs --with-wmsclient --with-wfsclient --with-sos

[4]make
[5]make instll

如果在[3][4]中遇到缺少包的错误，那么，再用aptitude install安装包就可以了.

zhyiwww 2009-04-03 14:53 发表评论

install GD-2.0.35 on debian 5.0

zhyiwww — Fri, 03 Apr 2009 05:38:00 GMT

[1]下载 GD
     wget http://www.libgd.org/releases/gd-2.0.35.tar.gz
[2]解压
     tar vxfz gd-2.0.35
[3]./configure --prefix=/opt/gd
[4]make
[5]make install

第[4]步出现的问题：
/bin/sh ./libtool --mode=link gcc -g -O2 -o gdparttopng gdparttopng.o
./libgd.la -lm
gcc -g -O2 -o .libs/gdparttopng gdparttopng.o ./.libs/libgd.so -lm
-Wl,--rpath -Wl,/usr/local/lib
gdparttopng.o(.text+0xac): In function `main':
/root/gd-2.0.21/gdparttopng.c:42: undefined reference to `gdImageCreateFromGd2Part'
collect2: ld returned 1 exit status
make[2]: *** [gdparttopng] Error 1
make[2]: Leaving directory `/root/gd-2.0.21'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/root/gd-2.0.21'
make: *** [all] Error 2
解决的方法：
     aptitude install libgd2-xpm libgd2-xpm-dev
安装后再编译，就可以了。

zhyiwww 2009-04-03 13:38 发表评论

Mapserver配置postgis数据源

zhyiwww — Fri, 20 Feb 2009 10:19:00 GMT

[1]把shp数据导入到postgres(postgis)
[2]road.map文件内容如下

# Map file generated by QGIS version 0.7.4
# Edit this file to customize for your interface
# Not all sections are complete. See comments for details.
NAME road
STATUS ON

# Map image size. Change size as desired
SIZE 256 256
#
# Start of web interface definition. Only the TEMPLATE parameter
# must be specified to display a map. See Mapserver documentation
#
WEB
# HEADER
# TEMPLATE
# FOOTER
MINSCALE 1
MAXSCALE 8
# Set IMAGEPATH to the path where mapserver should
# write its output
IMAGEPATH 'd:\\tmp\\'
# Set IMAGEURL to the url that points to IMAGEPATH
#IMAGEURL '/map_output/'
END

PROJECTION
    "init=epsg:4326"
END

# Extent based on full extent of QGIS view
EXTENT 115.275 39.2204 117.475 40.9462

UNITS METERS

IMAGECOLOR 255 255 255
IMAGETYPE gif
# Projection definition
# Projections are not currenlty supported. If desired, add your own
# projection information based on Mapserver documentation.
#
LAYER
NAME road
TYPE LINE
STATUS ON

CONNECTIONTYPE POSTGIS
CONNECTION "host=localhost dbname=mapdb user=postgres password=postgres port=5432"
DATA "the_geom from beijing_polyline using unique gid using srid=4326"

#DATA "d:\\wwwroot\\shp\\Rbeijing.shp"

CLASS
    NAME "ROAD"
    # TEMPLATE
    OUTLINECOLOR 112 238 238
END
END
END # Map File

蓝色的部分就是连接postgis数据库的部分，和取数据的部分

[3]road.html

访问http://localhost:8888/postgis/road.html，可以看到地图结果如下:

zhyiwww 2009-02-20 18:19 发表评论

mapserver - loadLayer()

zhyiwww — Fri, 20 Feb 2009 09:40:00 GMT

在windows下面配置mapserver，用shp做数据源，总是出现如下的错误，
loadLayer(): Unknown identifier. Parsing error near (:):(line 46)
尝试多次也没有找到结果，最后，找到了。
主要是在map的文件
我的错误的配置是:

# Map file generated by QGIS version 0.7.4
# Edit this file to customize for your interface
# Not all sections are complete. See comments for details.
NAME road
STATUS ON

# Map image size. Change size as desired
SIZE 256 256
#
# Start of web interface definition. Only the TEMPLATE parameter
# must be specified to display a map. See Mapserver documentation
#
WEB
# HEADER
# TEMPLATE
# FOOTER
MINSCALE 1
MAXSCALE 8
# Set IMAGEPATH to the path where mapserver should
# write its output
IMAGEPATH 'd:\\tmp\\'
# Set IMAGEURL to the url that points to IMAGEPATH
#IMAGEURL '/map_output/'
END

PROJECTION
    "init=epsg:4326"
END

# Extent based on full extent of QGIS view
EXTENT 115.275 39.2204 117.475 40.9462

UNITS METERS

IMAGECOLOR 255 255 255
IMAGETYPE gif
# Projection definition
# Projections are not currenlty supported. If desired, add your own
# projection information based on Mapserver documentation.
#
LAYER
NAME road
TYPE LINE
STATUS ON
DATA d:\\wwwroot\\shp\\Rbeijing.shp
CLASS
    NAME "ROAD"
    # TEMPLATE
    OUTLINECOLOR 112 238 238
END
END
END # Map File

正确的配置是：

# Map file generated by QGIS version 0.7.4
# Edit this file to customize for your interface
# Not all sections are complete. See comments for details.
NAME road
STATUS ON

# Map image size. Change size as desired
SIZE 256 256
#
# Start of web interface definition. Only the TEMPLATE parameter
# must be specified to display a map. See Mapserver documentation
#
WEB
# HEADER
# TEMPLATE
# FOOTER
MINSCALE 1
MAXSCALE 8
# Set IMAGEPATH to the path where mapserver should
# write its output
IMAGEPATH 'd:\\tmp\\'
# Set IMAGEURL to the url that points to IMAGEPATH
#IMAGEURL '/map_output/'
END

PROJECTION
    "init=epsg:4326"
END

# Extent based on full extent of QGIS view
EXTENT 115.275 39.2204 117.475 40.9462

UNITS METERS

IMAGECOLOR 255 255 255
IMAGETYPE gif
# Projection definition
# Projections are not currenlty supported. If desired, add your own
# projection information based on Mapserver documentation.
#
LAYER
NAME road
TYPE LINE
STATUS ON
DATA "d:\\wwwroot\\shp\\Rbeijing.shp"
CLASS
    NAME "ROAD"
    # TEMPLATE
    OUTLINECOLOR 112 238 238
END
END
END # Map File

这个就是正确的配置，仅仅是引号的差别。

zhyiwww 2009-02-20 17:40 发表评论

在postgis使用google 900913投影

zhyiwww — Thu, 19 Feb 2009 13:36:00 GMT

在postgis中是使用google epsg:900913 投影
平台：debian
[1]编辑/opt/proj/share/proj/epsg文件
注意：
我安装的proj在/opt/proj目录，所以，就编辑/opt/proj/share/proj/epsg文件
如果你用aptitude install proj，那么就编辑 /usr/share/proj/epsg文件

在文件的结尾添加

# Google Mercator
<900913> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs <>

［2］在postgis数据库里面，添加如下记录

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, srtext, proj4text) values (900913 ,'EPSG',900913,'GEOGCS["WGS 84", DATUM["World Geodetic System
1984", SPHEROID["WGS 84", 6378137.0, 298.257223563,AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], NIT["degree",0.017453292519943295], AXIS["Longitude", EAST], AXIS["Latitude", NORTH],AUTHORITY["EPSG","4326"]], PROJECTION["Mercator_1SP"],PARAMETER["semi_minor", 6378137.0],
PARAMETER["latitude_of_origin",0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor",1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0],UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH],AUTHORITY["EPSG","900913"]] |','+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m
+nadgrids=@null +no_defs');

zhyiwww 2009-02-19 21:36 发表评论

用shape文件做为数据源在Mapserver发布WMS服务

zhyiwww — Fri, 06 Feb 2009 08:19:00 GMT

系统
Linux+MapServer-5.2.1+shpfile
［1］安装apache2.2

［2］安装mapserver-5.2.1
把mapserv复制到${apache_home}/cgi-bin

［3］用qgis导出mapscript文件

# Map file generated by QGIS version 0.7.4
# Edit this file to customize for your interface
# Not all sections are complete. See comments for details.
MAP
NAME bbeijing_polyline
STATUS ON

# Map image size. Change size as desired
SIZE 800 600
#
# Start of web interface definition. Only the TEMPLATE parameter
# must be specified to display a map. See Mapserver documentation
#
WEB
# HEADER
# TEMPLATE
# FOOTER
MINSCALE 1
MAXSCALE 13
# Set IMAGEPATH to the path where mapserver should
# write its output
IMAGEPATH '/tmp/'
# Set IMAGEURL to the url that points to IMAGEPATH
#IMAGEURL '/map_output/'
END

PROJECTION
    "init=epsg:4326"
END

# Extent based on full extent of QGIS view
EXTENT 115.275 39.2204 117.475 40.9462

UNITS METERS

IMAGECOLOR 255 255 255
IMAGETYPE gif
# Projection definition
# Projections are not currenlty supported. If desired, add your own
# projection information based on Mapserver documentation.
#
LAYER
NAME bbeijing_polyline
TYPE LINE
STATUS ON
DATA /opt/wwwroot/map_shp/Bbeijing_polyline.shp
CLASS
    NAME "Bbeijing_polyline"
    # TEMPLATE
    OUTLINECOLOR 112 238 238
END
END
END # Map File

［4］在OpenLayers中调用地图文件

zhyiwww 2009-02-06 16:19 发表评论

在debian上安装mapserver5.2.1

zhyiwww — Thu, 05 Feb 2009 10:39:00 GMT

[1]configure

./configure --prefix=/opt/mapserver-5.2.1 --with-gd=/opt/gd-2.0.35/ --with-wmsclient --with-wms --with-png --with-jpeg --with-agg --with-proj=/opt/proj-4.6.0/ --with-gdal=/opt/gdal-1.6.0/bin/gdal-config --with-ogr=/opt/gdal-1.6.0/bin/gdal-config --with-freetype=/opt/freetype/ --with-curl-config=/opt/curl-7.19.3/bin/curl-config --with-geos=/opt/geos-3.0.3/bin/geos-config --with-wfs --with-wfsclient --with-wcs --with-sos --with-xml2-config=/opt/libxml/bin/xml2-config --with-tiff --with-postgis=/opt/PostgreSQL/8.3/bin/pg_config

[2] make

[3]make install

zhyiwww 2009-02-05 18:39 发表评论

GeoServer中配置Oracle Spatial数据库

zhyiwww — Fri, 14 Nov 2008 06:55:00 GMT

[1]安装GeoServer1.7

download http://geoserver.org/display/GEOS/Download
install http://geoserver.org/display/GEOSDOC/1.1+Install+GeoServer

[2]安装Oracle spatial plugin for Geoserver

download oracle plugin for geoserver from http://sourceforge.net/project/downloading.php?groupname=geoserver&filename=geoserver-1.6.5-oracle-plugin.zip&use_mirror=nchc
config
1. Copy gt2-oracle jar to [GEOSERVER_HOME]/server/geoserver/WEB-INF/lib/
2. copy ojdbc14.jar to [GEOSERVER_HOME]/lib
3. restart GeoServer

[3]导入空间数据
http://www.blogjava.net/zhyiwww/archive/2008/11/12/240106.html
[4]创建空间索引
     这一步是最关键也是最重要的步骤。
     CREATE INDEX idx_beijng_polyline_geom ON beijng_polyline(geom) INDEXTYPE IS MDSYS.SPATIAL_INDEX;

[5]创建命名空间

[6]创建Datastore for oracle

查看是否配置成功

[7]配置FeatureType
   选择刚才配置好的oracle数据库连接

编辑FeatureType内容

[8]查看地图

如果配置成功，就可以看到下面的地图

zhyiwww 2008-11-14 14:55 发表评论

通过mapinfo的easyloader把shp数据导入到oracle

zhyiwww — Wed, 12 Nov 2008 08:03:00 GMT

[1]在mapinfo中打开shp文件，并保存成tab文件
[2]打开easyloader工具，
打开eaayloader

[3]配置数据源，到oracle spatial的连接

注意：用户要有足够的权限，能够创建Map Catalog

[4]导入数据

zhyiwww 2008-11-12 16:03 发表评论

oracle用户管理

zhyiwww — Wed, 12 Nov 2008 06:34:00 GMT

[1]创建删除用户
   create user usera identified by password;
   drop user usera;
   // 创建用户，锁定状态
   create user userb identified by password account lock;

zhyiwww 2008-11-12 14:34 发表评论

把shp格式的数据导入到oracle

zhyiwww — Wed, 12 Nov 2008 06:06:00 GMT

[1]利用shp2sdo把shp格式的数据转换成oracle能够识别的数据格式，利用下面的工具
shp2sdo -o c:\tmp\beijing\beijing_polyline c:\oracle_sql\beijng_polyline -s 4326

c:\tmp\beijing\beijing_polyline 是你的shp文件
c:\oracle_sql\beijng_polyline 是你要生成的目标文件

命令成功执行后，会生成三个文件
c:\oracle_sql\beijng_polyline.sql
c:\oracle_sql\beijng_polyline.dat
c:\oracle_sql\beijng_polyline.ctl

[2]修改文件
原始的beijng_polyline.ctl文件是
LOAD DATA CHARACTERSET UTF8
INFILE c:\oracle_sql\beijng_polyline.dat
TRUNCATE
CONTINUEIF NEXT(1:1) = '#'
INTO TABLE BEIJNG_POLYLINE
FIELDS TERMINATED BY '|'
TRAILING NULLCOLS (
   MAPID     NULLIF MAPID = BLANKS,
   ID     NULLIF ID = BLANKS,
   KIND     NULLIF KIND = BLANKS,
   GEOM COLUMN OBJECT
   (
     SDO_GTYPE       INTEGER EXTERNAL,
     SDO_SRID        INTEGER EXTERNAL,
     SDO_ELEM_INFO   VARRAY TERMINATED BY '|/'
       (X            FLOAT EXTERNAL),
     SDO_ORDINATES   VARRAY TERMINATED BY '|/'
       (X            FLOAT EXTERNAL)
   )
)
我修改后的文件是

LOAD DATA CHARACTERSET UTF8（这个地方是设置编码）
INFILE beijng_polyline.dat（这个地方很重要，一定要把路径去掉，否则的话总是有错误）
TRUNCATE
CONTINUEIF NEXT(1:1) = '#'
INTO TABLE BEIJNG_POLYLINE
FIELDS TERMINATED BY '|'
TRAILING NULLCOLS (
   MAPID     NULLIF MAPID = BLANKS,
   ID     NULLIF ID = BLANKS,
   KIND     NULLIF KIND = BLANKS,
   GEOM COLUMN OBJECT
   (
     SDO_GTYPE       INTEGER EXTERNAL,
     SDO_SRID        INTEGER EXTERNAL,
     SDO_ELEM_INFO   VARRAY TERMINATED BY '|/'
       (X            FLOAT EXTERNAL),
     SDO_ORDINATES   VARRAY TERMINATED BY '|/'
       (X            FLOAT EXTERNAL)
   )
)

[3]执行导入
在命令行下，改变目录到数据所在的目录，如
c:\oracle_sql\
然后执行
sqlldr zy/password@orcl beijing_polyline
就可以了，如果成功，就会出现如图所示的结果

zhyiwww 2008-11-12 14:06 发表评论

如何从postgis数据导出成shp格式

zhyiwww — Wed, 12 Nov 2008 05:32:00 GMT

要实现从postgress中导出空间数据成shp格式的数据，可以使用如下的命令

pgsql2shp.exe -f beijing_polyline -h localhost -u postgres -P xxxxx mapdb public
.beijing_polyline

-h 数据库服务器
-u 用户名
-p 密码
-f 参数，可以指定导出后的文件名称
mapdb 是数据库名称
public.beijing_polyline是数据表的名称

执行，就可以把数据导出成Shpae格式的文件了。

不过，要实现这一步，此数据库需要时空间数据库，也就说需要在创建数据库的时候创建成一个postgis的数据库。

zhyiwww 2008-11-12 13:32 发表评论

windows环境下MapServer和MapLab的配置

zhyiwww — Tue, 11 Nov 2008 10:07:00 GMT

[1]下载安装mapserver
下载http://maptools.org/ms4w/index.phtml?page=downloads.html
安装到 d:\MS4W\MS4W
[2]下载maplab,http://www.maptools.org/maplab/
[3]解压maplab-2.2.1,复制到d:\ms4w\ms4w\apache\htdocs\
[4]在浏览器打开http://localhost/maplab-2.2.1/htdocs/
可以看到

[4]点击配置maplab,如下图所示：

[5]配置maplab,mapbrowser,mapedit

[6]查看配置，如果都是绿色的“√”，就说明使正确的

ok,如果你看到这个，说明配置成功了。
最可能不好配置的就是mapscript的配置了。

zhyiwww 2008-11-11 18:07 发表评论

在postgres上创建索引

zhyiwww — Fri, 07 Nov 2008 05:32:00 GMT

[1]创建普通索引
   CRREATE INDEX idx_beijing_polyline_mapid ON beijing_polyline(mapid);
[2]创建空间数据索引
   CREATE INDEX idx_beijing_polyline_geom ON beijing_polyline USING GIST(the_geom);
   需要先安装postgis支持。

zhyiwww 2008-11-07 13:32 发表评论

how to install Mapserver5 on Fedora9

zhyiwww — Sat, 25 Oct 2008 14:47:00 GMT

http://blog.londatiga.net/it/installing-mapserver-and-php-mapscript-on-linux/

http://mapserver.gis.umn.edu/docs/howto/freebsd62mapserver410

zhyiwww 2008-10-25 22:47 发表评论

DoubleBits.cpp:51: error: ‘memcpy’ was not declared in this scope

zhyiwww — Sat, 25 Oct 2008 13:27:00 GMT

xian xiang :

DoubleBits.cpp: In static member function 'static double
geos::index::quadtree::DoubleBits::powerOf2(int)':
DoubleBits.cpp:51: error: 'memcpy' was not declared in this scope
DoubleBits.cpp: In constructor
'geos::index::quadtree::DoubleBits::DoubleBits(double)':
DoubleBits.cpp:94: error: 'memcpy' was not declared in this scope
make[3]: *** [DoubleBits.lo] Error 1
make[3]: *** Waiting for unfinished jobs....
 i686-pc-linux-gnu-g++ -DHAVE_CONFIG_H -I. -I. -I../../../source/headers
-I../../../source/headers/geos -I../../../source/headers -O2 -march=i686
-mtune=athlon-xp -msse -mmmx -pipe -ggdb3 -DGEOS_INLINE -Wall -ansi -pedantic
-Wno-long-long -MT IntervalSize.lo -MD -MP -MF .deps/IntervalSize.Tpo -c
IntervalSize.cpp -o IntervalSize.o >/dev/null 2>&1
make[3]: Leaving directory
`/var/tmp/portage/sci-libs/geos-3.0.0/work/geos-3.0.0/source/index/quadtree'
make[2]: *** [all-recursive] Error 1
make[2]: Leaving directory
`/var/tmp/portage/sci-libs/geos-3.0.0/work/geos-3.0.0/source/index'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory
`/var/tmp/portage/sci-libs/geos-3.0.0/work/geos-3.0.0/source'
make: *** [all-recursive] Error 1

jie jue :

http://bugs.gentoo.org/show_bug.cgi?id=217834

http://bugs.gentoo.org/attachment.cgi?id=149838

zhyiwww 2008-10-25 21:27 发表评论

如何取得空间两条新的交点

zhyiwww — Fri, 06 Jun 2008 01:47:00 GMT

近段，在学习geotools，去实现一些空间的运算。
下面的一段程序，实现了求得空间两条线的交点的程序，仅仅是一个demo,所以，就固定了两个线段的坐标，然后去求其交点。
需要导入geotools下的包，其实主要是用的是JTS。

package com.cybersoftchina.jts.demo;

import java.util.ArrayList;
import java.util.List;

import org.geotools.geometry.jts.FactoryFinder;

import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.Geometry;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.LineString;
import com.vividsolutions.jts.io.ParseException;
import com.vividsolutions.jts.io.WKTReader;

public class LineCross {

    public void getCrossPoints() {


        Coordinate[] coordinates = new Coordinate[] {
        new Coordinate(116.29440, 39.90131),
        new Coordinate(116.41457, 39.84741) };
        Geometry g1 = new GeometryFactory().createLineString(coordinates);

        Coordinate[] coordinates1 = new Coordinate[] {
        new Coordinate(116.42418, 39.9136),
        new Coordinate(116.26968, 39.8395) };
        Geometry g2 = new GeometryFactory().createLineString(coordinates1);

        Geometry g = g1.intersection(g2);
        System.out.println(g);


    }

    public void lineCross(){



    }

    public static void main(String[] args) {
        LineCross lc = new LineCross();
        lc.getCrossPoints();
    }

}

程序运行，返回的结果是：
POINT (116.34822145369408 39.87716939623774)
如果没有交点，那么返回的结果会告诉你是empty.

zhyiwww 2008-06-06 09:47 发表评论