冰川流量计算实习报告.doc

冰川流量计算实习报告 小组成员：陈恺磊 顾亚宁 王俊达 郑智园 2014/12/13 一、实习目的 计算穿过着地线附近的冰川流量。 二、实习数据 1、grid_speed: 冰流速； 2、PCMEGA_icethi：冰厚度数据1； 3、Russian_icethi_modify：冰厚度数据2； 4、Study_area：研究区； 5、GL_Rignot：着地线位置； 6、Basin_LGB01：流域界线； 7、Flow_lines：冰流线； 8、amm1_v2_clip:卫星影像。 三、实习平台 ArcGIS 10.1 四、成员分工 数据预处理:陈恺磊 顾亚宁 着地线绘制：陈恺磊 王俊达 数据插值：郑智园 顾亚宁 实习报告整理：郑智园 PPT制作：王俊达 五、实习步骤 （一）数据预处理 数据的裁剪与掩膜 由于实习数据的量过大，而所要计算的只是研究区内的数据，因此利用分析工具→裁剪 进行裁剪。 将数据全部加载进入ArcMap； 打开ArcToolBox，选择分析工具→提取分析→裁剪 以上操作即可完成冰流线的裁剪；重复上述过程，完成着地线的裁剪。结果如图： （3）打开ArcToolBox，选择Spatial Analysist Tools→提取分析→按掩膜提取，对两个栅格数据进行掩膜。操作过程及结果如下图： 2、冰厚度数据合并 由于两个冰厚度数据的属性表中的字段不一致，不能进行合并。因此： （1）先将两个属性表中不需要的字段进行删除（留下FID、Shape、经纬度及冰层厚度字段）； （2）然后添加三个double型字段，分别命名为LATITUDE，LONGITUDE，ICE_THICK(注意：字段名称要与另一属性表的字段名完全一致，否则不能进行联结)； （3）右击新添加的字段，选择栅格计算器导入数据；删除原来的三个属性，使两个属性表保持一致性。 （4）打开ArcToolBox,选择数据管理工具→常规→合并，进行图层合并。 打开合并之后的冰厚度数据的属性表（下图），可看到两个图层已完成联结，并且属性表中有LATITUDE，LONGITUDE，ICE_THICK这些之后操作所需要的重要的数据。 3、生成规则网格（Fishnet） （1）由于研究区域是没有投影坐标的，因此要先对研究区进行投影定义。 （2）打开ArcToolBox,选择数据管理工具→投影和变换→定义投影 （3）修改图层属性，将地图单元改为米。过程为：右击Layers,打开属性，选择常规选项卡。 （4）生成Fishnet（1公里×1公里网格） 打开ArcToolBox,选择数据管理工具→要素类→创建渔网 4、冰厚度数据抽稀 （1）将图放大可以看到：一个网格内有不止一个点，这样的工作量是十分繁琐的，因此要使用一定的方法进行数据抽稀，使每个网格内至多只有一个点。 （2）由于渔网坐标是平面坐标，厚度合并数据是椭球体坐标，因此先要统一为平面X、Y坐标。 打开合并后的冰厚度数据的属性表，增加X和Y字段，类型为双精度。右击X字段，选择字段计算器，完成转换；同样的方法对Y进行转换。 （3）数据抽稀 打开ArcToolBox,选择分析工具→叠加分析→空间链接 分别右击LATITUDE，LONGITUDE，ICE_THICK，规则选择为平均值 （4）抽稀之后数据是以面图层的形式呈现，要把它转化成点要素。右击打开数据抽稀之后图层的属性表，将属性表导出为.dbf格式 （5）从File中以添加X，Y数据的方式，将抽稀之后的点添加进来。 （二）着地线绘制 1、新建一个线图层来绘制着地线； 2、在Layers中打开显示老师给的着地线位置点图层和冰流线图层 将所有的线段合并成为1条后进行分割 打开ArcToolbox，选择数据管理工具→要素→要素折点转点，生成各分段节点 栅格插值 不同栅格插值的比较并选择 普通克里金插值法 简单克里金插值法 泛克里金插值法 1）普通克里金插值法预测图 在ArcGIS中实现步骤为： （1）、在geostatistical analyst中的地统计向导命令 、一直点击下一步，直到产生如下图结果 （3）、看到标准平均值接近于0，标准均方根接近于1，符合检验拟合度的标准，单击完成。 （4）、在生成的普通克里金插值预测图上右单击数据---导出栅格。导出栅格数据 2）简单克里金插值法预测图 、在geostatistical analyst中的地统计向导命令。 （2）、一直点击下一步，直到产生如下图结果 （3）、看到标准平均值接近于0，标准均方根与1相差较大，不符合检验拟合度的标准，单击完成。 （4）、在生成的普通克里金插值预测图上右单击数据---导出栅格。导出栅格数据 3）泛克里金插值法预测图 （1）、在geostatistical analyst中的地统计向导命令。 （