基于Mapnik创建切片的方案

最新推荐文章于 2023-06-19 16:53:04 发布

智能空间

最新推荐文章于 2023-06-19 16:53:04 发布

阅读量1.3k

点赞数

分类专栏：地图服务文章标签： mapnik 切片 python

本文链接：https://blog.csdn.net/LIJUNGANG1019473506/article/details/80083860

版权

2 篇文章 0 订阅

订阅专栏

  操作系统：win7 64bit ； python2.7.14 32bit； mapnik 2.2。 

2.1 生成格网的函数：

1）在指定范围内将坐标划分为若干等份生成格网：

  def createGrid02(envelopeParam,widthCount=1,heightCont=1): 

"""

  把输入的mapnik.Box2d按指定的行列数目进行划分，最后返回划分的格网 

"""

  reslut={} 

  print envelopeParam.minx,envelopeParam.maxx 

  s=abs(envelopeParam.maxx-envelopeParam.minx) 

  width=abs(envelopeParam.minx-envelopeParam.maxx)/(widthCount*1.0) 

  height=abs(envelopeParam.maxy-envelopeParam.miny)/(heightCont*1.0) 

  for i in range(0,widthCount): 

  for j in range(0,heightCont): 

  box=mapnik.Box2d(envelopeParam.minx+i*width,envelopeParam.miny+j*height,envelopeParam.minx+(i+1)*width,envelopeParam.miny+(j+1)*height) 

  reslut[(i,j)]=box 

  return reslut 

2）在指定范围内根据指定等级自动生成格网：

  def createGrid03(envelopeParam,zoom): 

"""

  根据输入的Box2d范围 和瓦片等级zoom 

  生成在该范围的上对应等级的格网1~zoom(不包含zoom 级) 

  每一级对应的切片横纵坐标数目是Power(2,zoom-1) 

  返回值是字典{（行号，列号，zoom）：box} 

  行、列、等级号从1开始 

"""

  reslut={} 

  for i in range(1,zoom): 

  num=2**(zoom-1) 

  tempReslut=createGrid02(envelopeParam,widthCount=num,heightCont=num) 

  for item in tempReslut: 

  reslut[(item[0][0],item[0][1],i)]=tempReslut[item] 

  return reslut 

3）从shapefile数据读取规则格网转换为mapnik的box，可以通过shapefile文件自定义生成切片的范围：（需要安装python环境下的gdal）

  def getBoxFromShp(shapeFile): 

"""

  根据shapefile的格网生成grid列表 

  返回boxs的列表(fid,Box2d) 

"""

  gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","NO") 

  gdal.SetConfigOption("SHAPE_ENCODING","") 

  #i=0 

  ogr.RegisterAll() 

  ds=ogr.Open(shapeFile,0) 

  if ds==None: 

  return 

  olayer=ds.GetLayerByIndex(0) 

  if olayer==None: 

  return 

  olayer.ResetReading() 

  oFeature=olayer.GetNextFeature() 

  if not (olayer.GetGeomType()==3): 

  return 

  boxlist=[] 

  while oFeature is not None: 

  fid=oFeature.GetFID() 

  print fid 

  oGeometry=oFeature.GetGeometryRef() 

  envelope=oGeometry.GetEnvelope() 

  box=mapnik.Box2d(envelope[0],envelope[2],envelope[1],envelope[3]) 

  boxlist.append((fid,box)) 

  oFeature=olayer.GetNextFeature() 

  return boxlist 

  根据以上的函数基本可以满足自定义格网要求 

2.2 生成切片

通过不断的更改map的zoom范围然后产生不同位置上的切片，方法如下：

map.zoom_to_box(Box b)# b是mapnik 中的Box
mapnik.save_map(map,mapXml)
mapnik.render_to_file(map,png)

关注