geojson合并python的简单介绍

本文目录一览：

• 1、Geojson笔记二：geojson-python-util
• 2、python 合并两个json文件
• 3、怎么样用Python实现地理编码
• 4、如何生成geojson格式的tiles
• 5、有哪些 GIS+Python 的开发经验值得分享

Geojson笔记二：geojson-python-util

随着geojson变得越来越流行，需要给geojson一些具体的数据操作。我在写这库的工程中也是学习的过程，希望也能给你的GIS学习一点点帮助。

将脚本 geojson_utils.py 复制到你的文件夹中，然后输入：

验证两条线是否相交并求出交点。 原理

判断点是否在多边形内部。原理大致如下：

判断点是否在多个多边形内

通过一个中心点和半径获得一个圆形的图形

取出矩形的中心点

两点的互联网球面距离

(参考 )，注意这个是球型距离，非大地陀球距离

判断点线面是否在某点半价内

求多边形面积

多边形中心点

adapted from

通过起点，距离和角度来计算终点

python 合并两个json文件

先用json.load解析各自的文件内容,

然后用第一个接触出来的的列表，extend方法第二个列表，然后两个列表就合二为一了。

怎么样用Python实现地理编码

工具原料：地理编码工具geocoder用Python实现地理编码方法如下：一、看一下它都支持哪些公司的地理编码服务:二、安装pipinstallgeocoder三、地理编码importgeocoderg=geocoder.google(“1403WashingtonAve,NewOrleans,LA70130″)g=geocoder.arcgis(u”北京市海淀区上地十街10号”)g.latlng输出为[29.9287839,-90.08421849999999]也可以查看完整的geojsong.geojson输出为{‘bbox’:[-90.0855674802915,29.9274349197085,-90.0828695197085,29.9301328802915],’geometry’:{‘coordinates’:[-90.08421849999999,29.9287839],’type’:’Point’},’properties’:{‘accuracy’:u’ROOFTOP’,’address’:u’1403WashingtonAve,NewOrleans,LA70130,USA’,’bbox’:[-90.0855674802915,29.9274349197085,-90.0828695197085,29.9301328802915],’city’:u’NewOrleans’,’confidence’:9,’country’:u’US’,’county’:u’OrleansParish’,’encoding’:’utf-8′,’housenumber’:u’1403′,’lat’:29.9287839,’lng’:-90.08421849999999,’location’:’1403WashingtonAve,NewOrleans,LA70130′,’neighborhood’:u’GardenDistrict’,’ok’:True,’place’:u’ChIJGyFHWc2lIIYRYSoneaXAUiw’,’postal’:u’70130′,’provider’:’google’,’quality’:u’street_address’,’state’:u’LA’,’status’:’OK’,’status_code’:200,’street’:u’WashingtonAve’},’type’:’Feature’}直接用Google尝试查询中文地址时失败g=geocoder.google(u”北京市海淀区上地十街10号”)g.ok输出为False用百度应该没问题，不过我没有申请相应的key。切换到arcgis，能够成功编码g=geocoder.arcgis(u”北京市海淀区上地十街10号”)g.latlng输出为[40.050934,116.30079]逆地理编码g=geocoder.google([29.9287839,-90.08421849999999],method=’reverse’)printg.addressprintg.cityprintg.stateprintg.country输出为1403WashingtonAve,NewOrleans,LA70115,USANewOrleansLAUS换成中国的地址g=geocoder.google([40.050934,116.30079],method=’reverse’)printg.addressprintg.cityprintg.stateprintg.country输出为BaiDuDaSha,HaidianQu,BeijingShi,China,100193BeijingBeijingShiCN用arcgis的服务试试g=geocoder.arcgis([40.050934,116.30079],method=’reverse’)printg.addressprintg.cityprintg.stateprintg.country输出为None北京市北京市CHNGoogle转换成的是英文，但地址比较全。arcgis虽然是中文，但是详细的地址居然输出为了None，这有个X用。其他geocoder的功能不止于此，它还可以查IP（包括自己的）。g=geocoder.ip(‘199.7.157.0’)printg.latlngprintg.cityg=geocoder.ip(‘me’)printg.latlngprintg.city输出为[43.6934,-79.4857]Toronto[51.05,13.75]Dresden查询一个城市的空间包围盒g=geocoder.arcgis(u”山东”)g.bbox输出为{‘northeast’:[38.976997,121.976998],’southwest’:[33.022997,116.022998]}

如何生成geojson格式的tiles

转载自 需下载的工具有：python：python-2.6

Python Setup tools ：setuptools-0.6c11.win32-py2.6

GDAL：GDAL-1.9.2.win32-py2.6

TileStache：TileStache-1.42.1

1.首先安装python2.6 并配置环境变量

2.安装Setup tools （方面安装其他软件）和GDAL

3.使用Setup tools安装TileStache, ModestMaps, PIL(TileStache 依赖于这两项)

easy_install tilestache

easy_install PIL

easy_install ModestMaps

4.在下载的TileStache中的文件夹中配置cfg文件

{

“cache”:

{

“name”: “Disk”,

“path”: “/tmp/stache”,

“umask”: “0000”

},

“layers”:

{

“postgis”:

{

“provider”: {“name”: “vector”, “driver”: “PostgreSQL”,

“parameters”: {“dbname”: “postgis”, “user”: “postgres”,”password”:”postgres”,

“table”: “mainroad”}}

}

}

}

}

5.然后用tilestache-seed进行切片

6.结果

后面会推出矢量切片在OL中的应用。OpenLayers中加载矢量切片，可有效解决大数据量的问题

有哪些 GIS+Python 的开发经验值得分享

GISer入门指南电子书第一季 分享版.pptx免费下载

链接:

提取码: wttn 

ArcGIS入门级教程,详细介绍了ArcGIS产品族及成员产品 1、GIS 的概念和需求 2、什么是ArcGIS 3、Geodatabase中的GIS数据概念 4、桌面GIS产品:ArcView, ArcEditor和ArcInfo等