亚洲国产人成在线观看,亚洲乱亚洲乱妇24p,无码欧精品亚洲日韩一区夜夜嗨

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,有空间字�D�the_geom,但是是基于SRS4326的，我想�?00913d的投影，所以就学要��d��一个新的空间字�D��ƈ转换投媄坐标。做法如下：

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)

上面的表�l�构是添加和转换后的。添加空间字�D늚�同时�Q�也为空间字�D�|��加了索引�?br />
操作前提�Q�postgis/postgres能支�?00913�Q�如果不支持的话�Q�就学要先添�?00913投媄支持�?br />

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

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

在调试mapfile的时候，往往不知道怎么才能看到mapserver的执行细节，所以，能查看日至不�׃ؓ一个好的办法�?br /> 在mapserver中��用日志的�Ҏ��如下:

(1)在MAP中添加输出日��x��?br />
CONFIG "MS_ERRORFILE" "/tmp/map_error.log"

(2)在LAYER里面打开DEBUG
也就是设�|�DEBUG的��别�?br />
DEBUG 5
�Ҏ��mapserver的参考文档，可以讄��ON|OFF|0|1|2|3|4|5.
调试阶段�Q�设�|�成5�Q�可以得到最多的输出信息�?br />
(3)实时查看日志
tail -f /tmp/map_error.log

(4)错误调试
比如�Q�我��d��一个错误的标签AAAA,�pȝ��输出日志如下�Q?br /> [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)
。。。。。。。。。。。。�?br />
�q�样�Q�你��可以根据错误信息很�Ҏ��的找到错误�ƈ更正错误�?br />

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

mapserver-�l�道路添加label

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

Mapfile�?需要��用如下配�|?

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

LABELITEM "pathname"

(2)在MAP中添加字体引�?br />
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

同上一��?a id="viewpost1_TitleUrl" class="postTitle2" href="../../zhyiwww/archive/2009/08/19/291809.html">mapfile中��用子查询实现道�\分��?br /> 在mapfile中，也可以用FILTER来数据的�q��o�Q�通过此方法也可以实现道�\的分�U��?br />
DATA "the_geom from rbeijing as subquery using unique gid USING SRID=4326"
FILTER "pathclass='3'"

pathclass是表rbeijing中的字段的名�U��?是其��|��表示道�\的一�U��别�?br />

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

ogr2ogr

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

在linux下面�Q�默认的是没有ogr2ogr命��o的，如果想用�Q�就要自己安装�?br />
在ubuntu下面�Q�可以用
aptitude install gdal-bin -y
安装完成��可以用ogr2ogr命��o了�?br />

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

mapfile中��用子查询实现道�\分��

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

查看道�\分��如下�Q?br />
mapdb=# select distinct pathclass from rbeijing;
pathclass
-----------
0
15
2
3
4
(5 rows)

使用子查询实现道路分�U?br />
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 发表评论

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

(1)主要在mapscript(beijing.map)文�g里面��d��如下配置�Q?br />         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
                }
            );

效果�Q�请参看上篇文章

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

让你的mapserver支持AGG

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

�pȝ��q�_��Q�在ubuntu上面
(1)aptitude install mapserver libagg-dev

(2)在map文�g里面��d��输出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
                }
            );

前后的效果比�?br /> 没有使用AGG支持的效�?br />

使用AGG支持的效�?br />

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

shp2pgsql的SRID参数

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

#shp2pgsql Rbeijing Rbeijing > rbeijing.sql
导出文�g如下�Q?br />     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
导出文�g如下�Q?br />     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字段的��g��因�ؓSRID的不同而不同�?br />

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 发表评论

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

我的错误日志如下�Q?br />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

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

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

Postgis的安装和配置

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

把下面的代码存�ؓ install_postgis.sh
然后执行�?br />必须在postgres用户下执行�?br />

######################################################################
## 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 发表评论

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

debian 下安装mapserv5.2.2的简单方�?br />[1]安装gd
参�?br />[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]中遇到缺��包的错误，那么�Q�再用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

�W�[4]步出现的问题�Q?br /> /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
安装后再�~�译�Q�就可以了�?br />

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

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

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

# 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数据库的部分�Q�和取数据的部分

[3]road.html

讉K��http://localhost:8888/postgis/road.html�Q�可以看到地囄��果如�?

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

mapserver - loadLayer()

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

在windows下面配置mapserver�Q�用shp做数据源�Q��L��出现如下的错误，
loadLayer(): Unknown identifier. Parsing error near (:):(line 46)
��试多次也没有找到结果，最后，扑ֈ�了�?br />主要是在map的文�?br />我的错误的配�|�是:

# 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

正确的配�|�是�Q?br />
# 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

�q�个��是正确的配�|�，仅仅是引��L��差别�?br />

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

在postgis使用google 900913投媄

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

在postgis中是使用google epsg:900913 投媄
�q�_��Q�debian
[1]�~�辑/opt/proj/share/proj/epsg文�g
注意�Q?br />我安装的proj�?opt/proj目录�Q�所以，��q��?opt/proj/share/proj/epsg文�g
如果你用aptitude install proj�Q�那么就�~�辑 /usr/share/proj/epsg文�g

在文件的�l�尾��d��

# 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 <>

�Q?�Q�在postgis数据库里面，��d��如下记录

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文�g做�ؓ数据源在Mapserver发布WMS服务

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

�pȝ��
Linux+MapServer-5.2.1+shpfile
�Q?�Q�安装apache2.2

�Q?�Q�安装mapserver-5.2.1
把mapserv复制�?{apache_home}/cgi-bin

�Q?�Q�用qgis导出mapscript文�g

# 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

�Q?�Q�在OpenLayers中调用地图文�?br />

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 发表评论

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]导入�I�间数据
http://www.tkk7.com/zhyiwww/archive/2008/11/12/240106.html
[4]创徏�I�间索引
     �q�一步是最关键也是最重要的步骤�?/font>
     CREATE INDEX idx_beijng_polyline_geom ON beijng_polyline(geom) INDEXTYPE IS MDSYS.SPATIAL_INDEX;

[5]创徏命名�I�间

[6]创徏Datastore for oracle

查看是否配置成功

[7]配置FeatureType
   选择刚才配置好的oracle数据库连�?br />

�~�辑FeatureType内容

[8]查看地图

如果配置成功�Q�就可以看到下面的地�?br />

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

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

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

[1]在mapinfo中打开shp文�g�Q��ƈ保存成tab文�g
[2]打开easyloader工具�Q?br />打开eaayloader

[3]配置数据源，到oracle spatial的连�?br />
注意�Q�用戯��有��够的权限�Q�能够创建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;
   // 创徏用户�Q�锁定状�?br />   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能够识别的数据格式，利用下面的工�?br />shp2sdo -o c:\tmp\beijing\beijing_polyline c:\oracle_sql\beijng_polyline -s 4326

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

命��o成功执行后，会生成三个文�?br />c:\oracle_sql\beijng_polyline.sql
c:\oracle_sql\beijng_polyline.dat
c:\oracle_sql\beijng_polyline.ctl

[2]修改文�g
原始的beijng_polyline.ctl文�g�?br />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�Q�这个地�Ҏ��讄��~�码�Q?br /> INFILE beijng_polyline.dat�Q�这个地方很重要�Q�一定要把�\径去掉，否则的话��L��有错误）
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]执行导入
在命令行下，改变目录到数据所在的目录�Q�如
c:\oracle_sql\
然后执行
sqlldr zy/password@orcl beijing_polyline
��可以了�Q�如果成功，��׃��出现如图所�C�的�l�果

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

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

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

要实��C��postgress中导出空间数据成shp格式的数据，可以使用如下的命�?br />
pgsql2shp.exe -f beijing_polyline -h localhost -u postgres -P xxxxx mapdb public
.beijing_polyline

-h 数据库服务器
-u 用户�?br />-p 密码
-f 参数�Q�可以指定导出后的文件名�U?br />mapdb 是数据库名称
public.beijing_polyline是数据表的名�U?br />
执行�Q�就可以把数据导出成Shpae格式的文件了�?br />
不过�Q�要实现�q�一步，此数据库需要时�I�间数据库，也就说需要在创徏数据库的时候创建成一个postgis的数据库�?br />

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

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,如下图所�C�：

[5]配置maplab,mapbrowser,mapedit

[6]查看配置�Q�如果都是绿色的“√”，��p��明��正确�?br />

ok,如果你看到这个，说明配置成功了�?br />最可能不好配置的就是mapscript的配�|�了�?br />

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

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

[1]创徏普通烦�?br />   CRREATE INDEX idx_beijing_polyline_mapid ON beijing_polyline(mapid);
[2]创徏�I�间数据索引
   CREATE INDEX idx_beijing_polyline_geom ON beijing_polyline USING GIST(the_geom);
   需要先安装postgis支持�?br />

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 发表评论

如何取得�I�间两条新的交点

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

�q�段�Q�在学习geotools�Q�去实现一些空间的�q�算�?br />下面的一�D늨�序，实现了求得空间两条线的交点的�E�序�Q�仅仅是一个demo,所以，��固定了两个�U�段的坐标，然后��L��其交炏V�?br />需要导入geotools下的包，其实主要是用的是JTS�?br />

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();
    }

}

�E�序�q�行�Q�返回的�l�果是：
POINT (116.34822145369408 39.87716939623774)
如果没有交点�Q�那么返回的�l�果会告诉你是empty.

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