ENVI软件中,手动划定地面控制点从而实现栅格图像相互间地理配准的方法;其中,所用软件版本为ENVI Classic 5.3 (64-bit)。
Link Displays”后在其中一幅图像中随机点击;此时可以看到两幅图的同一位置并不是同样的地物,而是具有一定空间位置差异,如下图所示。
Map”→“Registration”→“Select GCPs: Image to Image”;如果其中一景图像有地理信息而另一景没有(例如一景遥感影像与一幅.jpg
格式的图层),就需要选择“Select GCPs: Image to Map”。
Base Image”与“Warp Image”,亦即基准图层与需要变换的图层,在这里我们分别选中前述两景图像即可,具体二者谁是“Base Image”谁是“Warp Image”并没有强制要求;但是一定要牢记这里的设置,在后期还会用到。
GCP)选择窗口,此时就可以在图像显示区域中选择GCP了。
Link Displays”关闭后才可以选择GCP。
Add Point”即可。
Show List”,可以看到当前已经找到的GCP。
GCP的各类信息。
GCP在图中显示得不是很明显,可以通过“Set Point Colors”进行设置。
GCP的过程,从而找到更多的GCP。
Degree”设置为2
会有比较好的效果(这里“Degree”指的是用于计算RMS误差的次数或阶数,2
就指的是用二次多项式来计算误差);进一步的,RMS误差就是下图中“RMS Error”,其表示地理配准过程中,控制点原始位置与转换后控制点新位置间的像素差值,因此其越小越好。
GCP后,我们就可以用“Predict”进行辅助操作:在第一景图像中找到第一个点后,通过“Predict”就可以自动定位到第二景图像的对应位置附近,随后手动微调即可。
Auto Predict”。
GCP列表中,选中某一行GCP后,可以通过“Goto”实现直达这一GCP位置的功能。
GCP,我们可以通过“On/Off”将其暂时取消(没错,不是删除,是暂时不加入该点)。
GCP后,可以选择将GCP列表导出为文本格式:
GCP列表的过程是可选的,而接下来的操作则是必须的——我们需要保存GCP(这里就不是上面的那个GCP列表了,而是各个GCP的信息)为.pts
格式。
.pts格式的GCP信息后,之后如果我们需要再次修改对应图层的GCP,直接导入即可。
Map”→“Registration”→“Warp from GCPs: Image to Image”。
.pts格式的GCP信息文件并选中。
Input Warp Image”窗口和“Input Base Image”窗口中,要按照前述选择“Base Image”与“Warp Image”时的设置进行选择——这就是为什么前面说需要牢记“Base Image”与“Warp Image”设置的原因。
ENVI中打开(如下所示最右侧图像为地理配准后图像),用“Link Displays”进行随机选择,可以看到最右侧的图与最左侧的基准图像空间位置几乎一致,说明大功告成。