使用QGIS将无人机航飞影像tiff转为gcj02火星坐标，并发布切片服务

一、导入数据

打开QGIS，将无人机航飞得到的影像数据和矢量范围导入，如果影像数据太大，可以给tiff建一个金字塔，在Processing Toolbox（处理工具箱）中搜索pyramids（金字塔），

二、坐标转换

一般无人机航飞的影像数据都是CGCS2000坐标系的，这里我们需要先把影像转换成web墨卡托投影坐标（EPSG:3857）

1. tiff栅格的坐标转换：

如果数据比较大，则处理速度会比较慢，可以选择开启多线程处理，会稍微快一些。

这里勾选上就可以开启多线程。

2. 矢量范围的坐标转换：

三、地理配准

1. 配准之前需要准备好控制点数据，控制点需要用到wgs84的矢量范围和gcj02的矢量范围，wgs84转gcj02的转换请自行百度，我一般用qgis的geohey插件或者91卫图助手。

2. 使用提取顶点工具，把上一步中转换成3857投影的矢量范围顶点提取出来，然后使用字段计算器把提取出来的顶点坐标计算出来：

打开属性表，字段计算器来获取顶点坐标：

3. 对转换成gcj02的矢量范围做同样的操作，但是需要先把gcj02的矢量范围同第二步中的操作一样，把4326坐标转换成3857。

4. 现在我们拥有了两份顶点数据，一个是原始的3857顶点数据，一个是加上gcj02火星坐标偏移的顶点数据：

现在我们需要把这两份顶点数据做成控制点给我们的栅格数据地理配准使用：

打开导出的两份数据，把from的数据复制到to的后面：

结果如上图

5. 然后我们做一个控制点的样本文件，到时候把这里的坐标替换掉

首先打开地理配准菜单，把刚才转换成3857的栅格数据导入

点击地图添加一个坐标点，然后随意输入或者从地图上个选择一个目标点，然后点这个保存控制点

保存之后用文本编辑器复制之前保存的cvs文件的内容，替换掉刚才保存的控制点文件：

替换之后再次添加数据的时候会自动把控制点带进来并计算出配准后的结果，

然后点这个

选择线性，选择输出位置，重点是不压缩，压缩有可能会配准失败，我试了很多次都失败了，然后等待输出即可。

四、地图裁剪

一般飞出来的影像是需要根据需求把范围裁剪一下的，我们需要用到第二步中转换成墨卡托投影的gcj02的矢量范围对栅格进行裁剪，这里我们用到的是栅格裁剪工具，栅格有两个裁剪工具，一个是裁剪出来矩形范围，一个是按矢量边界裁剪，这里我们用矢量边界裁剪：

然后又是漫长的等待。。。

五、地图切片

如果栅格数据很大，需要先构建金字塔，步骤同第一步的导入数据，然后就可以开始切片了，这里我们用生成XYZ tiles工具（生成到目录中）：

生成结束就可以查看成果了：