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

实验说明

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

实验步骤

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

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的资料较少,希望我的文章能对你有所帮助,如果有问题,大家可以评论一起交流下~

  • 6
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

编程小白呀

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值