将地质图（shp文件）叠加到天地图上

1.环境依赖

python2.7+pyshp模块
下载python2.7

2.配置环境

安装时可以直接add path。
在cmd中输入python,出现以下界面代表配置成功。
手动配置环境变量教程
网上有很多教程可以去搜搜。

环境配置好后，在cmd中输入python -m pip install pyshp -i http://mirrors.aliyun.com/pypi/simple --trusted-host mirrors.aliyun.com安装pyshp模块。
至此环境就配置好了。

3.主要代码

首先到开发目录开启服务器python -m CGIHTTPServer 8889，在浏览器中输入localhost:8889。
在你的开发目录新建index.html文件

<!DOCTYPE html>
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8"/>
<meta name="keywords" content="天地图"/>
<title>天地图－地图API－范例－经纬度直投地图</title>
<script type="text/javascript" src="http://api.tianditu.gov.cn/api?v=4.0&tk=您的密钥"></script>
<style type="text/css">body,html{width:100%;height:100%;margin:0;font-family:"Microsoft YaHei"}#mapDiv{width:100%;height:400px}input,b,p{margin-left:5px;font-size:14px}</style>
<script type="text/javascript" src="http://localhost:8889/cgi-bin/readShp.py"></script>
</head>
<body onLoad="onLoad()">
<div id="mapDiv"></div>
</body>
</html>

在你的开发目录新建一个文件cgi-bin，写的py文件必须放在这个目录。
readShp.py

#!/usr/bin/env python
# -*- coding:utf-8 -*-
# @Author: wardseptember
print "Content-Type:text/javascript"
printimport shapefile
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
sf = shapefile.Reader("./cgi-bin/shapefile/aotao.shp")
shapes = sf.shapes()first_str = '''var map;var zoom = 9;map = new T.Map('mapDiv', {projection: 'EPSG:32650'});map.centerAndZoom(new T.LngLat(112.858577, 35.497694), zoom);'''
name_str=''''''
for num in range(len(shapes)):name_str=name_str+'''points_autao%d=[];'''% (num)first_str=first_str+name_strpolyline_str = ''''''
for shp in range(len(shapes)):shap = shapes[shp]#一个面文件，里面可能有很多ploygon,shapes[shp]就是一个，一共有len(shapes)个ploygonfor i in range(len(shap.points)):#读每个ploygon中点坐标polyline_str=polyline_str+'''points_autao%d.push(new T.LngLat(%f, %f));'''% (shp,shap.points[i][0],shap.points[i][1])finally_str=''
for num in range(len(shapes)):finally_str=finally_str+'''var polygon_autao%d = new T.Polygon(points_autao%d,{color: "blue", weight: 3, opacity: 0.5, fillColor: "red", fillOpacity: 0.5});map.addOverLay(polygon_autao%d);'''% (num,num,num)print first_str+polyline_str+finally_str
#必须要挨个读取ploygon，不能一次读取，不然画出来的图形就乱了。一次压入所有点，就依点画，明显不对。

读其他ploygon图层写法

#!/usr/bin/env python
# -*- coding:utf-8 -*-
# @Author: wardseptemberprint "Content-Type:text/javascript"
printimport shapefile
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
sf = shapefile.Reader("./cgi-bin/shapefile/sandie.shp")
shapes = sf.shapes()first_str = '''
'''
name_str=''''''
for num in range(len(shapes)):name_str=name_str+'''points_sd%d=[];'''% (num)first_str=first_str+name_strpolyline_str = ''''''
for shp in range(len(shapes)):shap = shapes[shp]for i in range(len(shap.points)):polyline_str=polyline_str+'''points_sd%d.push(new T.LngLat(%f, %f));'''% (shp,shap.points[i][0],shap.points[i][1])finally_str=''
for num in range(len(shapes)):finally_str=finally_str+'''var polygon_sd%d = new T.Polygon(points_sd%d,{color: "blue", weight: 3, opacity: 0.5, fillColor: "#37FA3F", fillOpacity: 0.5});map.addOverLay(polygon_sd%d);'''% (num,num,num)print first_str+polyline_str+finally_str

读线文件

#!/usr/bin/env python
# -*- coding:utf-8 -*-
# @Author: wardseptember
print "Content-Type:text/javascript"
printimport shapefile
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
sf = shapefile.Reader("./cgi-bin/shapefile/JCborder.shp")
shapes = sf.shapes()
first_str = '''points = [];'''
polyline_str = ''''''
for shp in range(len(shapes)):shap = shapes[shp]for i in range(len(shap.points)):polyline_str=polyline_str+'''points.push(new T.LngLat(%f, %f));'''% (shap.points[i][0],shap.points[i][1])
finally_str='''var line = new T.Polyline(points);map.addOverLay(line);
'''print first_str+polyline_str

可以发现读线读点读面写法类似。
下面是批量读点的代码

#!/usr/bin/env python
# -*- coding:utf-8 -*-
# @Author: wardseptember
print "Content-Type:text/javascript"
printimport shapefile
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
sf = shapefile.Reader("./cgi-bin/shapefile/JCborder.shp")
shapes = sf.shapes()polyline_str = ''''''
lng_num=[0.0]
lat_num=[0.0]
total_points=1
for shp in range(len(shapes)):shap = shapes[shp]for i in range(len(shap.points)):total_points=total_points+1lng_num.append(shap.points[i][0])lat_num.append(shap.points[i][1])for j in range(total_points):polyline_str=polyline_str+'''var point_%d=new T.LngLat(%f, %f);var marker_%d = new T.Marker(point_%d);map.addOverLay(marker_%d);marker_%d.disableDragging();  '''% (j,lng_num[j],lat_num[j],j,j,j,j)polyline_str=polyline_str+'map.removeOverLay(marker_0)'
print polyline_str

将地质图（shp文件）叠加到天地图上相关推荐

shp文件显示 c语言,上传并在地图中显示Shp文件
前段时间参与了一个项目,客户有一个功能需求是上传SHP文件并在地图上显示,然后在此基础上做缓冲区处理.经过对比测试,最终选择了shapefile.js工具,在此做个记录. shapfe.js能够将Es ...
shp文件纯前端的上传、解析、编辑、下载
本文主要讲述一种体量较小的shp文件纯前端的上传.解析.编辑.下载的技术流程,适用于要素量少的shp文件修改操作. 准备工作下载一下几个包,详细用法请见结尾参考. npm install file- ...
上传并在地图中显示Shp文件
前段时间参与了一个项目,客户有一个功能需求是上传SHP文件并在地图上显示,然后在此基础上做缓冲区处理.经过对比测试,最终选择了shapefile.js工具,在此做个记录. shapfe.js能够将Es ...
arcgis导出shp文件_RegionManager GIS导出shp文件编码说明
在<RegionManager GIS数据如何上报到国家水土保持重点工程项目管理系统>一文介绍了RM GIS软件如何导出shp文件的相关操作,以及如何定义投影等相关操作.从RM GIS软件 ...
将shp文件转化为osm文件，并导入到sumo中建立路网
前言最近在做的一个课题要用到sumo仿真,正好我们手头也有arcgis的shp文件,比osm上面的要简洁准确不少,在座的过程中遇到不少问题,这里做一下记录要注意的地方有几个: shp转osm的工具 ...
osm服务器显示乱码,怎样获得osm上的行政区划shp文件
2018.4.18 :刚刚发现我这里的wifi(电信的宽带)打不开下面的其中几个网站,需要使用网络工具才行. 但是用手机流量(联通)是可以直接打开的. 2018.1.25 :方法二的网站又恢复了. ...
在Vs2017上集成osgearth3.2和qt5.9，并加载shp文件。
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档文章目录前言一.QT5.9在vs上部署 1.1 在Vs中下载插件 1.2 配置QT 二.OsgEarth3.2环境配置. 三.在Q ...
10显示不支持导入dwg文件_Global Mapper中80坐标系高程DEM与kml文件叠加实例
概述在GIS行业中,有很多经常用到的软件,比如CAD.ArcGIS.Erdas和ENVI等大型专业软件,也会常用到像Global Mapper这样的"迷你"软件,这里,我们以80 ...
tif 高程_Global Mapper中80坐标系高程DEM与kml文件叠加实例
概述在GIS行业中,有很多经常用到的软件,比如CAD.ArcGIS.Erdas和ENVI等大型专业软件,也会常用到像Global Mapper这样的"迷你"软件,这里,我们以80 ...

将地质图（shp文件）叠加到天地图上

1.环境依赖

2.配置环境

3.主要代码

将地质图（shp文件）叠加到天地图上相关推荐

最新文章

热门文章