[发明专利]一种电子地形图自动几何校正方法

权利要求书

1.一种电子地形图自动几何校正方法，包括下述步骤：

步骤一、打开一个栅格格式的扫描电子地形图，以新图幅号命名，如果获取的文件名不符合新图幅号的命名规范，将跳出错误提示；如果获取的文件名符合新图幅号的命名规范，根据该新图幅号计算获得四角图廓点，由该图幅号生成四角十字交叉点的经纬度坐标；

步骤二、模板匹配：利用模板匹配的CV_TM_SQDIFF方法，对外图廓四角进行初次匹配，匹配后通过四角内角点和已知经纬度的十字丝交叉点即图廓点，相对固定的距离通过相似度计算方法计算获得图廓点图面坐标；在二次模板匹配中，设置一个以初次获得的图面坐标点为中心的正方形ROI感兴趣区；在该区域内使用十字丝模板再次匹配，获得图廓点图面坐标；

步骤三、在ENVI+IDL环境中完成配准：经过以上步骤，已经获得了四对GCP点，将地形图的路径和GCP点对传入ENVI系统，利用ENVI自有函数envi_register_doit进行配准，ENVI_PROJ_CREATE函数和创建经纬度坐标系统，最终完成电子地形图自动几何校正。

2.如权利要求1所述的方法，步骤一中所述的栅格格式的扫描电子地形图为.Jpg或.bmp或.GIFf格式文件。

3.如权利要求1所述的方法，步骤一中所述新图幅号命名规范依据的标准为：1992年12月，我国颁布了《国家基本比例尺地形图分幅和编号GB/T13989-92》文件，1993年3月开始实施。

4.如权利要求1所述的方法，步骤一中已知新图幅号计算该图幅西南图廓点的经纬度，按下式计算该图幅西南图廓点的经纬度：

B=（b-31）*6°+（d-1）*△B；

L=（a-1）*4°+（4°/△L-c）*△L；

式中：B-------图幅内某点的经度；

L-------图幅内某点的纬度；

△B-------所求比例尺地形图分幅的经差；

△L-------所求比例尺地形图分幅的纬差；

a-------1：100万地形图图幅所在纬度带字符码所对应的数字码；

b-------1：100万地形图图幅所在经度带的数字码；

c-------所求比例尺地形图在1：100万地形图图号后的行号；

d-------所求比例尺地形图在1：100万地形图图号后的列号。

5.如权利要求1所述的方法，步骤二中相似度计算方法优选平方差方法。