【QGIS】使用“创建网格”来对一个区域分割得到正方形网格

实验说明

我现在拥有的文件是北京行政区矢量地图（beijing.shp），现在我想对地图从上到下进行分割，从而得到150*150的网格四个点的坐标范围。实验得到的结果如下图：

实验步骤

1. 使用创建网格命令
2. 我们选择矩形，当然你也可以选择六边形；网格范围就选择你要进行分割的地图shp文件就可以。
   
3. 我想要得到150*150米的网格，因此我需要选择3857的坐标系。
   
4. 选择完坐标系后，输入你想得到的正方形长度，运行就可以
5. 运行成功后，我们得到了网格文件，将网格文件导出为csv文件，我们发现得到的四个点坐标不是我想要的经纬度格式，因此我导出时选择了一些参数
   
6. 导出后的第一个字段为WKT，是我们想要的经纬度格式。（当然，你也可以对之前得到的坐标进行转换）
7. 如下图，我们得到了网格的区域范围，第一个字段为WKT，第2个字段为网格的id号。我发现QGIS中的创建网格是从上到下来进行分割的。
   

mercator坐标转换

各种坐标转换blog

def mercatortowgs84(x,y):
    """
    墨卡托投影坐标转回wgs84
    :param x:
    :param y:
    :return:
    """
    lng = x / 20037508.34 * 180
    lat = 180 / math.pi * (2 * math.atan(math.exp(y / 20037508.34 * 180 * math.pi / 180)) - math.pi / 2)
    return lng,lat


def wgs84togcj02(lng, lat):
    """
    WGS84转GCJ02(火星坐标系)
    :param lng:WGS84坐标系的经度
    :param lat:WGS84坐标系的纬度
    :return:
    """
    if out_of_china(lng, lat): # 判断是否在国内
        return lng, lat
    dlat = transformlat(lng - 105.0, lat - 35.0)
    dlng = transformlng(lng - 105.0, lat - 35.0)
    radlat = lat / 180.0 * pi
    magic = math.sin(radlat)
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt(magic)
    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
    mglat = lat + dlat
    mglng = lng + dlng
    return [mglng, mglat]

近几天在实验时，发现QGIS的资料较少，希望我的文章能对你有所帮助，如果有问题，大家可以评论一起交流下~