|
阅读：1155回复：0
ArcGIS画分幅格网并以坐标命名
问题：我需要在坐标范围X（450000m—568000m）、Y（4560000m—4770000m）的矩形区域内画出1:2000的影像图分幅格网，格网命名规范为×××. ×-×××. ×，取其右下角坐标中的一部分。
假设一个格网的右下角坐标为（557000，4618000），则此格网名称为618.0-557.0。557.0为557000/1000保留小数点后一位，618.0为（4618000-4000000）/1000保留小数点后一位，Y在前，X在后。此为我此文目的。 A．arctoolbox里面的数据管理工具——要素类——创建鱼网。  1.输出要素类，就是保存索引格网的位置，可见我把我要生成的索引格网保存成了“全市域分幅索引.shp”文件。 2.模板范围，这个是可选的，我选了“区县界线.shp”这个文件，因为我本来就是要以区县界线的范围生成分幅索引的，加载了这个文件之后，可以看到下面模板范围和鱼网原点坐标之间的左右上下里面都有了数字，这就是我的区县界线.shp文件所占的图幅坐标范围，左X坐标45204.298400——右X坐标567023.813900，上Y坐标4767334.872000——下Y坐标4562653.983400的大矩形框之间，根据这个可以参考出我要生成索引网格的坐标范围，按区县界线.shp文件的图幅坐标范围来看，我的索引网格的坐标范围可以更小一些，只要能涵盖区县界线.shp的范围就可以了。 3.鱼网原点坐标，这个比较关键，就是鱼网的左下角坐标，根据我题目中的范围，X坐标为450000，Y坐标为4560000。 4.Y轴坐标，这个很具有迷惑性，这个坐标主要是用来定向鱼网的，按照鱼网原点坐标与Y轴坐标的连线来定向鱼网的旋转角度，Y轴坐标和鱼网原点坐标不能相同，我生成的鱼网是不带旋转的，所以这里的X坐标依旧为450000，Y坐标有了一些变化，是4560010，这就是一根垂线了，而且Y轴Y坐标比鱼网原点Y坐标大，是向上的垂线，所以没有变化，如果Y轴Y坐标比鱼网原点Y坐标小，则垂线是向下的，那么就是生成鱼网要顺时针旋转180°了，这点一定要注意！ 5.像元宽度和像元高度，因为我要生成1:2000的索引网格，图面图幅大小是500mm×500mm，对应的实际图幅大小是1000m×1000m，像元宽度为1000，像元高度为1000。 6.行数和列数，行数=DY/1000m=（4770000m-4560000m）/1000m=210；列数=DX/1000m=（568000m-450000m）/1000m=118。 7.创建标注点可以不选，选的话，会生成一个点文件，我觉得没有什么用，我是忘了勾掉了。几何类型可以选线（POLYLINE），我选的是面（POLYGON）。 8.点击确定，生成最基本的索引格网。 B．为生成的索引格网命名。 1. arctoolbox里面的数据管理工具——字段——添加字段。给文件“全市域分幅索引”添加一个字段名为“name”，类型为“text”的字段。  当然，也可以鼠标放在“添加字段”上，右键选择批处理。我用批处理加了三个字段，name（text），X（LONG），Y（LONG）。  2.在arcmap窗口中打开全市域分幅索引的属性表。在表头X、Y位置鼠标右键点击，选择计算几何。X的属性是“质心X的坐标”，Y的属性是“质心Y的坐标”，坐标系选择“使用数据源的坐标系”，单位是“米”，确定就可以。这样，就把每个格网的中心坐标赋值给了字段X、Y。  3.还是属性表中，在表头name位置点击鼠标右键，选择字段计算器。解析程序用VB脚本（因为python我还不会用），网格编号应该是（（Y-40000000）/1000）.0-（X/1000）.0，转换成VB代码是“cstr(([Y]-4000000-500)/1000)&".0-"&cstr(([X]-500)/1000)&".0"”  4.用arctoolbox里面的数据管理工具——字段——删除字段，把X、Y字段删除。  5.结果如下图。  6.因为我想使这个网格适用范围广一些，所以，我在arccatalog窗口右键查看“全市域分幅网格.shp”文件的属性，在“XY坐标系”栏，点击清除（将坐标系设置为“未知”），清除了“全市域分幅网格.shp”文件的坐标系，这样的话，任何坐标系下，这个索引网格都可以应用。 结束。 难点是创建鱼网中的Y轴坐标和字段计算器中的VB脚本。 文章转载于新浪博客 版权归原作者所有，如有侵权请告知删除 |
|
|
