osm数据下载 python_osm数据下载 python_批量下载osm的分区域的osm文件

# coding: utf-8

#!/usr/bin/python

# #### http://www.openstreetmap.org 的分区后的osm数据批量下载工具(需要python的requests和BeautifulSoup包)。

# 从 http://download.geofabrik.de/ 提取网页中下载连接，然后写入osm-from-geofabrik.txt。

# 然后可以批量下载：** wget -c -i osm-from-geofabrik.txt ** ，或者载入其它的下载工具。

#

# 或者使用crontab定时下载,编写一个geofabrik-dl.txt(内容如下,指定晚上11:00,每周一到周五执行)。

# 然后执行：crontab crontab.txt。

# ```

# 0 23 * * 1-5 wget -c -N -i osm-from-geofabrik.txt

# 0 6 * * 1-5 killall wget

# ```

# 原始数据下载,参考：** http://planet.openstreetmap.org/planet/ **,  最新版本：http://planet.openstreetmap.org/planet/2016/planet-160411.osm.bz2

#

# ##### by openthings@163.com,http://my.oschina.net/u/2306127/blog

from bs4 import BeautifulSoup as bs

import requests as req

from pprint import *

#获取网页内容，经常访问不到，可以将地址拷贝到浏览器试试。

url = "http://download.geofabrik.de/"

try:

r = req.get(url)

print(r)

except Exception as ex:

print("Error:",ex)

#获取网页中的链接列表

bsr = bs(r.text,"html.parser")

alist = bsr.find_all('a')

#pprint(alist)

#另一种方式获取网页中链接的列表。

#bst = [x.extract() for x in bsr.find_all('a')]

#pprint(bst)

#查找列表

flist = open("osm-from-geofabrik.txt","w")

for a in bst:

#help(a)

a1 = a.attrs['href']

if a1.find('.bz2') > 0:

flist.writelines(a1+"\r\n")

print(a1)

flist.close()