Walk卫星影像定位方法.docVIP

Walk卫星影像定位方法 Walk卫星影像定位方法 1 （一）下载影像 1 （二）构建影像金字塔 2 （三）引入影像 2 （四）重定向 5 （五）工程化影像定向 9 （一）下载影像 Google影像分级为0～19，相邻级别影像分辨率相差1倍。下载18级的影像，其影像分辨率0.6米。 先在硬盘创建影像根目录，如：..\zjg1万\ 打开“稻歌Google Map截获”软件，找到紫金港校区： 进行【软件设置】，如图，将分幅大小【宽】=10000，【高】=10000： 输入 【角点一】=30.292497, 120.062112 【角点二】=30.317766, 120.104685 【影像级别】=18 点击【确定】，指定保存目录。 （二）构建影像金字塔 下载的影像为18级，像素分辨率0.6米，图像尺寸为80MB左右。为使得加载到软件的影像显示快速平滑，需要对下载的卫星影像构建影像金字塔。 在WalkField软件运行目录下的“gdal金字塔和切图”目录下双击“gomyway.wsp”脚本文件，系统启动WalkLan.exe脚本运行器： 执行脚本(F5)，选择【源图像…】： 然后执行【添加金字塔…】，脚本自动计算出合适的金字塔层数，如下图： 确认后，脚本完成“浙大紫金港.tif”的金字塔构建。 （三）引入影像 打开小组工程数据库，选择一个图层，如image图层。在工作空间将image图层置为可编，右键菜单执行【加入图象】。 弹出【引入影像】对话框。 标准影像定位过程 【选择文件…】，然后可通过【源影像定位信息…】查看影像的定位信息： 源影像是从GoogleMap下载的卫星影像，为Web Mercator投影，杭州城市坐标系采用的是西安80坐标系，卫星影像图与杭州坐标轴不平行，且纵横轴的尺度比也不同，因此需要先将Web Mercator投影到西安80，然后再将西安80转换到紫金港建筑坐标系。 Walk采用国际通用的proj4投影变换系统，在proj4中Web Mercator投影定义为： +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +no_defs 西安80以120度为中央子午线的高斯投影（不含带号）定义为： +proj=tmerc +lon_0=120 +x_0=500000 +y_0=0 +ellps=IAU76 在【引入影像】对话框中先【加入…-】卫星影像投影： 因为WebMercator投影在proj4中无选项对应，因此在【投影坐标系设置】对话框中可通过【编辑: 投影命令】直接输入投影定义。或通过【加载…】预制的投影，如上图中的“google map false”填充投影命令。注意，一定要去掉“独立平面直角坐标系”！ 然后，在【引入影像】对话框中再【加入…-】，在【投影坐标系设置】中设置西安80投影： 投影为tmerc横轴墨卡托，椭球为IAU76，中央子午线lon_0=120度，伪东x_0 源影像经WebMercator投影到西安80的120度中央子午线高斯坐标下，执行引入，然后在影像图层右键菜单的【浏览查询本层文字…】选中引入的影像： 卫星影像已投影到40带西安80高斯平面直角坐标系下。 （四）重定向 将西安80高斯坐标变换到本地坐标 本地坐标与西安80高斯坐标为同一大地基准（椭球和地心七参数相同），但平面坐标系的原点不同，且两个坐标的横轴不平行，因此需要进行平移和旋转。平移旋转需要影像和本地矢量图同名点匹配，如本地矢量图的房屋角点与影像图上同房屋地面角点匹配。 确保卫星影像所在图层为可编辑，选中影像后，点击鼠标右键菜单【图像配准(H)…】，弹出【重定向(可编层选中图象)】对话框，使用“多点法”，“3点仿射坐标变换”。 将地方坐标地形图在窗口居中，然后点击【重定向】对话框的“501. 矢量点平移到窗口中心”，将图像四角的矢量点平移到地形图周围。 比较地形图和影像，选择四个定位点。定位点【1】为中央灯塔底座中心，鼠标拖地形图上的定位点1到灯塔中心，严格居中。然后，执行“503. 图像点中心居中”，切换到图像，鼠标拖图像上的定位点1到灯塔中心，放大严格居中。 执行“502 矢量点中心居中”回到地形图，鼠标拖定位点【2】到生物实验中心的西南角，打开捕捉（s键），严格定位到地形图上。然后，执行“503. 图像点中心居中”，切换到图像，鼠标拖图像上的定位点2到生物实验中心的西南角，放大严格对中。 执行“502 矢量点中心居中”回到地形图，鼠标拖定位点【3】为生命科学院东南楼南凸之西南角，打开捕捉（s键），严格定位到地形图上。然后，执行“503. 图像点中心居中”，切换