BlogJava-zhyiwww -随笔分类-Postgishttp://www.blogjava.net/zhyiwww/category/37784.html用平实的笔，记录编程路上的点点滴滴………zh-cnTue, 25 Aug 2009 19:48:53 GMTTue, 25 Aug 2009 19:48:53 GMT60mapserver-给道路添加labelhttp://www.blogjava.net/zhyiwww/archive/2009/08/25/292550.htmlzhyiwwwzhyiwwwTue, 25 Aug 2009 11:54:00 GMThttp://www.blogjava.net/zhyiwww/archive/2009/08/25/292550.htmlhttp://www.blogjava.net/zhyiwww/comments/292550.htmlhttp://www.blogjava.net/zhyiwww/archive/2009/08/25/292550.html#Feedback0http://www.blogjava.net/zhyiwww/comments/commentRss/292550.htmlhttp://www.blogjava.net/zhyiwww/services/trackbacks/292550.html 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 发表评论
]]>shp2pgsql的SRID参数http://www.blogjava.net/zhyiwww/archive/2009/07/10/286275.htmlzhyiwwwzhyiwwwFri, 10 Jul 2009 08:31:00 GMThttp://www.blogjava.net/zhyiwww/archive/2009/07/10/286275.htmlhttp://www.blogjava.net/zhyiwww/comments/286275.htmlhttp://www.blogjava.net/zhyiwww/archive/2009/07/10/286275.html#Feedback0http://www.blogjava.net/zhyiwww/comments/commentRss/286275.htmlhttp://www.blogjava.net/zhyiwww/services/trackbacks/286275.html导出文件如下：
     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学习http://www.blogjava.net/zhyiwww/archive/2009/07/10/286270.htmlzhyiwwwzhyiwwwFri, 10 Jul 2009 07:49:00 GMThttp://www.blogjava.net/zhyiwww/archive/2009/07/10/286270.htmlhttp://www.blogjava.net/zhyiwww/comments/286270.htmlhttp://www.blogjava.net/zhyiwww/archive/2009/07/10/286270.html#Feedback0http://www.blogjava.net/zhyiwww/comments/commentRss/286270.htmlhttp://www.blogjava.net/zhyiwww/services/trackbacks/286270.html#select AddGeometryColumn('','gtest','geom',4326,'LINESTRING',2);


# INSERT INTO gtest (ID, NAME, GEOM)
VALUES (
   1,
   ’First Geometry’,
   GeomFromText(’LINESTRING(2 3,4 5,6 5,7 8)’, -1)
);


# select astext(geom) as geom from gtest;
              geom              
---------------------------------
 LINESTRING(2 3,4 5)
 LINESTRING(2 3,4 5,5 6)
 LINESTRING(2 3,4 5,5 6,3 5)
 LINESTRING(2 3,4 5,5 6,3 5,1 6)
(4 rows)

# create index idx_geom on gtest using GIST(geom);
CREATE INDEX

# select asbinary(geom) as geom from gtest;
                                                                                                                                                                  geom                                                                                                                                                                 
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 \001\002\000\000\000\002\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@
 \001\002\000\000\000\003\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@
 \001\002\000\000\000\004\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@\000\000\000\000\000\000\010@\000\000\000\000\000\000\024@
 \001\002\000\000\000\005\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@\000\000\000\000\000\000\010@\000\000\000\000\000\000\024@\000\000\000\000\000\000\360?\000\000\000\000\000\000\030@
(4 rows)

# select GeomFromText('POINT(116.8839 36.88299)',900913);
                    geomfromtext                   
----------------------------------------------------
 010100002031BF0D00CD3B4ED191385D40F758FAD005714240
(1 row)


# select asEWKB(geom) as geom from gtest;
                                                                                                                                                                        geom                                                                                                                                                                        
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 \001\002\000\000 \346\020\000\000\002\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@
 \001\002\000\000 \346\020\000\000\003\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@
 \001\002\000\000 \346\020\000\000\004\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@\000\000\000\000\000\000\010@\000\000\000\000\000\000\024@
 \001\002\000\000 \346\020\000\000\005\000\000\000\000\000\000\000\000\000\000@\000\000\000\000\000\000\010@\000\000\000\000\000\000\020@\000\000\000\000\000\000\024@\000\000\000\000\000\000\024@\000\000\000\000\000\000\030@\000\000\000\000\000\000\010@\000\000\000\000\000\000\024@\000\000\000\000\000\000\360?\000\000\000\000\000\000\030@
(4 rows)

# select asEWKT(geom) as geom from gtest;
                   geom                   
-------------------------------------------
 SRID=4326;LINESTRING(2 3,4 5)
 SRID=4326;LINESTRING(2 3,4 5,5 6)
 SRID=4326;LINESTRING(2 3,4 5,5 6,3 5)
 SRID=4326;LINESTRING(2 3,4 5,5 6,3 5,1 6)
(4 rows)

# select GeomFromEWKT('SRID=4326;POINT(119.8782 36.29288)');
                    geomfromewkt                   
----------------------------------------------------
 0101000020E610000039D6C56D34F85D4084D382177D254240
(1 row)

# select GeomFromEWKT('SRID=900913;POINT(119.8782 36.29288)');
                    geomfromewkt                   
----------------------------------------------------
 010100002031BF0D0039D6C56D34F85D4084D382177D254240
(1 row)


# select 'SRID=4326;POINT(119.28983 36.99289)'::geometry;
                      geometry                     
----------------------------------------------------
 0101000020E610000094D920938CD25D404243FF04177F4240
(1 row)

# select 'SRID=900913;POINT(119.28983 36.99289)'::geometry;
                      geometry                     
----------------------------------------------------
 010100002031BF0D0094D920938CD25D404243FF04177F4240
(1 row)






zhyiwww 2009-07-10 15:49 发表评论
]]>让postgis支持900913投影http://www.blogjava.net/zhyiwww/archive/2009/07/09/286125.htmlzhyiwwwzhyiwwwThu, 09 Jul 2009 09:57:00 GMThttp://www.blogjava.net/zhyiwww/archive/2009/07/09/286125.htmlhttp://www.blogjava.net/zhyiwww/comments/286125.htmlhttp://www.blogjava.net/zhyiwww/archive/2009/07/09/286125.html#Feedback0http://www.blogjava.net/zhyiwww/comments/commentRss/286125.htmlhttp://www.blogjava.net/zhyiwww/services/trackbacks/286125.html# 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 发表评论
]]>Mapserver配置postgis数据源http://www.blogjava.net/zhyiwww/archive/2009/02/20/255879.htmlzhyiwwwzhyiwwwFri, 20 Feb 2009 10:19:00 GMThttp://www.blogjava.net/zhyiwww/archive/2009/02/20/255879.htmlhttp://www.blogjava.net/zhyiwww/comments/255879.htmlhttp://www.blogjava.net/zhyiwww/archive/2009/02/20/255879.html#Feedback0http://www.blogjava.net/zhyiwww/comments/commentRss/255879.htmlhttp://www.blogjava.net/zhyiwww/services/trackbacks/255879.html[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
<html xmlns="http://www.w3.org/1999/xhtml">
  <head>
    <style type="text/css">
        ##mapdiv {
            width: 800px;
            height: 475px;
            border: 1px solid black;
        }
    </style>
    <script src="/OpenLayers/OpenLayers.js"></script>

    <script type="text/javascript">
        //Google has 20 scales
        //Resolutions always start big to small, descending order
        //numZoomLevels=20 (corresponds to Googles 20 zoom levels)
        //map.setCenter(center, 19); (Zoom Level starts at 0 and ends in 19)

        window.onload = function() {

            map = new OpenLayers.Map('mapdiv');

            var serverURL = "/cgi-bin/mapserv";

            var layer_obj = new OpenLayers.Layer.WMS(
                "road",
                serverURL,
                {  
                    layers: 'road',
            map : 'd:\\wwwroot\\postgis\\road.map',
                    format: 'gif'
                },
                {
                    //reproject: false,
                    'numZoomLevels': 20,
                    gutter: 15,
                    buffer: 0
                }
            );
            layer_obj.setIsBaseLayer(true);  
            layer_obj.setVisibility(true);  
            map.addLayer(layer_obj);

            map.addControl(new OpenLayers.Control.PanZoomBar());
            map.addControl(new OpenLayers.Control.MousePosition());
            map.addControl(new OpenLayers.Control.LayerSwitcher());
            map.addControl(new OpenLayers.Control.Scale());
            var center = new OpenLayers.LonLat(116.39166, 39.90625);
            map.setCenter(center, 8);
        }
         
    </script>
  </head>
  <body>
    <div id="mapdiv"></div>
  </body>
</html>


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

mapserver_postgis.PNG






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