创建等间距的形状文件空间网格

Question

我是 GIS 的新手，我一直在用地图上的纬度、经度坐标绘制点。我想在地图上构建一个 500 m 正方形的空间网格，以便稍后将点与这个新的空间约束相交。输出应位于 .shp 文件中。

到目前为止，我一直在使用底图包轻松阅读，并以以下形式绘制一些已经制作的格网：

Code minx miny maxx maxy

说明每个正方形（代码）的顶点坐标。我无法弄清楚如何制作一个新网格，给定一个选定的空间区域为：

# Projection CYL Cylindrical Equal Distance Projections (PlateCarree)
m = Basemap(projection='cyl',
    llcrnrlat=41.1905,urcrnrlat=41.5404,
    llcrnrlon=1.9144,urcrnrlon=2.6628,
    resolution='c')

任何帮助将不胜感激。

Answer 1

最后，

我用Geopandas和底图管理了所有的事情。首先，我找到了一个脚本，将空间部分（在UTM区域31投影中）划分为500 m的等空间网格。然后我不得不将网格坐标转换为纬度/纬度，以便使其与我的点相匹配：

import geopandas as gpd
geofile_in = 'UTMgrid.shp'
geofile_out = 'LATLONgrid.shp'
g = gpd.read_file(geofile_in)
originalcrs = {u'units': u'm', u'ellps': u'WGS84', u'datum': u'WGS84', u'proj': u'utm', u'zone': 31}
targetcrs = {u'ellps': u'WGS84', u'datum': u'WGS84', u'proj': u'longlat'}
# Set the original crs (UTM Zone 31 N)
g.crs = originalcrs
# Transform the Grid to the target crs (Lon, Lat)
g.to_crs(crs=targetcrs, inplace=True)
# Save to .shp file
g.to_file(geofile_out)

希望有帮助。