logo资料库

IDRISI操作全过程.docx

第1页 / 共9页
第2页 / 共9页
第3页 / 共9页
第4页 / 共9页
第5页 / 共9页
第6页 / 共9页
第7页 / 共9页
第8页 / 共9页
资料共9页,剩余部分请下载后查看
三、IDRISI操作
1、首先在projects中创建一个工程目录(yuce)。
1、获取马尔科夫矩阵
四、实现CA_Markov模型预测
五、精度评价
IDRISI 软件操作步骤 一、软件获取 IDRISI 软 件 主 要 用 于 预 测 , 网 上 有 很 多 版 本 可 以 下 载 , 目 前 主 要 用 IDRISI17.0。以 2005、2010、2015 年三期数据为例,解释预测过程。 二、转化数据 1、先在 Arcgis 中 将 raster 文件进行分类,参考分类体系分成为 1、2、3、4、 5、6、7 类。 2、Arcgis 中 raster 文件转 ASCII 文件 Arcgis 的 raster 格式转为 IDRISI 的 raster 格式时,需要用到 Arcgis 的 ASCII 文件进行转换,因此在 Arcgis 中先将分类好的 raster 文件转为 ASCII 文 件。此项工作是在 arctoolbox 中进行的。 三、IDRISI 操作 1、首先在 projects 中创建一个工程目录(yuce)。
2、IDRISI 中数据格式的转换 在 IDRISI 中 将 ASCII 文 件转 为 IDRISI 中 的 raster 文 件。 Import → Software-Specific Formates→ESRI Formate→ARCRASTER。 在 Input file 中输入在 Arcgis 的 ASCII 中转换的文件(2005.txt), Output file 中输出文件。 点击 output reference information,选择投影参数,比如,原来文件投影 是 utm zone 50n 的,就要选这样的参数文件,其他不变,点击 OK 后运行得到 2005.rst 文件,即 IDRISI 中的栅格图像。
3、重分类 直接在右上角搜索栏中输入 RECLASS 做重分类 也可以选择 GIS Analysis→ Database query→reclass,输入要做的重分类的 rst 文件(2005),选择保存路 径。IDRISI 中转化的数据属性值是土地利用类型的编码,没有属性值的默认 -9999,这在下一步的运行中会出错,要进行重分类,将-9999 转换为 0,其他编 码不变,第一列是赋给的新值,第二列和第三列是旧地类编码的范围。
三、Markov,生成适宜性图集 1、获取马尔科夫矩阵 方法:Modeling→Environmental/Simulation models→MARKOV
1 表示重分类后的转换矩阵前一期影像,为 2005 年遥感影像; 2 表示重分类后的转换矩阵的后一期影像,为 2010 年遥感影像; 3 是这个模型中输出条件概率的前缀,表示的是从 05 到 10 变化的一些信息,自己命名, 比如说 0510; 4 表示第一个与第二个影像之间的时间间隔,这里为 5 年; 5 表示向前预测的时间周期,这里也设置为 5 年,即模拟 2015 年的土地利用情况; 6 是比例误差(我看的资料里面一般都设置的是 0.1)。 获取的马尔科夫矩阵记录了在下一个时期,从每个土地利用类型转换为其他土地利用类 型的概率,得到面积转移矩阵。 四、实现 CA_Markov 模型预测 土地利用变化模拟使用的是 IDRISI 软件中的 CA-Markov 模型, 选择 Modeling →Environmental/Simulation models→CA_Markov。 1 表示模拟影像需要依据的影像,即为我们的 10 年遥感影像; 2 表示马尔科夫转换矩阵面积文件,选择第一个,这里选择的是马尔科夫转换面积矩阵 0510transition_areas;
3 即为转换适宜性图集(把从 05 转换为 10 年影像中产生的那个 0510 文件作为适宜性 图集,一般都是自己重新做一个这种图集,需要道路、河流、坡度等信息,此步骤还没有做, 暂时精度上不知道哪个更高点); 4 表示输出的土地利用变化数据,命名为 15; 5 表示元胞自动机循环次数,一般为两个年份之间间隔的整数倍,这里可以取 5、10、 15 等,但是数字越大,需要的模拟时间越多; 6 是 CA 模型中邻域结构的设定,系统默认的是 5*5 型 五、精度评价 将预测的 2015 年数据与实际 2015 年数据进行精度判断,用来分析以这种方法预测其他 年份的数据可靠性。精度越高表示该种方法预测越准确,效果越好。 GIS Analysis→Database Query→CROSSTAB
1 表示第一幅影像数据,这里的 yuce2015 为模拟产生的数据; 2 表示第二幅影像数据,这里的 rec2015.rst 为实际的数据。 二者的顺序是可以调换的,里面有不同的评价指标值。模型运行后得出的结果 为 Kappa=0.9697。
六、IDRISI 导入 ARCGIS 文件的方法 (1)File -> Import ->Desktop Publishing Formats -> GEOTIFF/TIFF (2) 然后在弹出的模块中依次选择要导入的 ASCII 文件、输出的文件名、行列 数、起止坐标、投影信息等,Pelatte to export with image 随便选一个都行。 (3) IDRISI 系统自带几百个常用的投影文件(所在目录:IDRISI 安装目录 \Idrisi32\Georef),文件名后缀为“.ref”,可用 IDRISI 的 Edit 工具打开和
分享到:
收藏