实验说明
我现在拥有的文件是北京行政区矢量地图(beijing.shp),现在我想对地图从上到下进行分割,从而得到150*150的网格四个点的坐标范围。实验得到的结果如下图:
实验步骤
- 使用创建网格命令
- 我们选择矩形,当然你也可以选择六边形;网格范围就选择你要进行分割的地图shp文件就可以。
- 我想要得到150*150米的网格,因此我需要选择3857的坐标系。
- 选择完坐标系后,输入你想得到的正方形长度,运行就可以
- 运行成功后,我们得到了网格文件,将网格文件导出为csv文件,我们发现得到的四个点坐标不是我想要的经纬度格式,因此我导出时选择了一些参数
- 导出后的第一个字段为WKT,是我们想要的经纬度格式。(当然,你也可以对之前得到的坐标进行转换)
- 如下图,我们得到了网格的区域范围,第一个字段为WKT,第2个字段为网格的id号。我发现QGIS中的创建网格是从上到下来进行分割的。
mercator坐标转换
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的资料较少,希望我的文章能对你有所帮助,如果有问题,大家可以评论一起交流下~