安吉玉米地适宜性评价

1、淋出拜乾大菅ZHEJIANG NORMAL UNIVERSITY数字图像处理实验报告题目：农业开发活动的适宜性分析以玉米种植业为例学院：地理与环境科学学院专 业：地理科学学生姓名：耿欢学号：12280208农业开发活动的适宜性分析一一玉米种植业为例地理122班 12280208 耿欢一、实验目的熟练使用ArcGIS软件中的空间分析工具;学会制作土地适宜性评价分级图二、实验要求及内容建立安吉县玉米种植适宜性评价体系；完成玉米种植各影响因子的空间分布等级的栅格图像;完成玉米种植适宜性分析报告三、实验原理1.建立玉米种植业适宜性评价体系适宜性分析是指土地针对某种特定开发活动的分布，这些开发活动包括农

2、业 应用，城市化选址，作物类型布局，道路选线，选着重新造林的最适宜的土地等。 因此，建立适宜性分析模型，首先确定具体的开发活动，其次选择其影响因子， 然后评判某一地域的各个因子对这种开发活动的适应程度，以作为土地利用规划 决策的依据。（1）评价对象：玉米种植用地的土地适宜性评价，通过评价将研究区分出不同 的适宜性等级S1 （最适宜），S2 （次适宜）,S3 （临界适宜），和N （不适宜）。（2）评价方法：采用基于GIS的土地质量评价法，即将玉米作物生长有关主导 生态条件与土地质量相比照，从而评定土地的适宜性等级。（3）评价过程评价对象生态条件的调查。玉米属禾本科，为一年生草本，其主导生态条 件

3、如性喜高温，需水量大，要求土壤肥沃和图层疏松，其根系生长要求防止土壤 侵蚀等。确定评价对象的影响因素和因子。根据将玉米作物生长有关的生态条件与 土地质量相比照，除了温度可以通过季节调节外，其他因子如下图所示：图1玉米作物生长有关的影响因素和因子分析图利用GIS生成影响因素数据。和表1对应的空间数据，可以用过图2所示的操作过程，生成基于栅格单元的供水条件等级数据R1,同理，可以生成供肥条件等级数据R2,供氧条件等级 数据R3和土壤侵蚀等级数据R4.计算各因素的权重和贡献函数值。由于每个影响因素或土地质量Q对评价对 象的贡献大小和适宜程度不同，因此不同因素有不同的权重W,不同级别的因素 对评价对象

4、的贡献函数值P也不同。Q.对应的权重 W(W1=0.67,W2=0.17,W3=0.08,W4=0.08)各个等级的分值 P，Pij=100，P2i =67, P引=33, P =0.表1供水影响因子及属性分级表2因素，等级评价指标表供水级别距水源距离/m土层厚度/cmQ1Q2Q3Q4150167178.08.02100-20050-20244.811.45.45.43200-50020-10322.25.62.62.645001040000通过将不同土地质量数据R的等级，切换为表2中的指标值，通过下式计算出栅 格单元的评价分值。R(T)=湍 *，/T23，ni=1根据使G(P ) = 1-1

5、 R(T)-P /1001的值最大时P所对应的J,即为所求的某个T的适宜级别S.i四、具体操作步骤1、打开2010年安吉土地利用类型图，利用Arcmap对其进行操作， 根据距离水源地进行分级。在工具包中选择“buffer”工具，分别做100,200,500米的缓冲区。分别见下图1-1,1-2,1-3：曰bufferZOOE) buffer 500图1-1: 100米缓冲区图1-2： 200米缓冲区曰乡Layers曰bufferlOObuffer200图1-3: 500米缓冲区2、选取水源（河流，水库，沟渠），使用Erase工具分别得出距离水 源（河流，水库，沟渠）小于100，100-200，2

6、00-500和大于500米 且不包含水源本身的图像。之后与旱地利用“ Intersect”工具，得 出距离水源小于100，100-200，200-500和大于500米的旱地图象。 见下图 2-1,2-2,2-3,2-4：日O _ _，.-S0 bufferl 00小于fLayers |曰 buffer200-100日 口 buffer500_200曰 buffer500大于x日乡 bufferlOOJ、于&-0 buffer200-100 buffer5OO_2OO buffer5OO?3-图2-1: 小于100米旱地图2-2： 100-200米旱地日 bufferl 00454紫色土，潮土6

7、、查阅资料，根据土壤孔隙度的等级，将土壤类型分级，见表6：表6 土壤孔隙度等级表土壤孔隙度等级土壤类型1石灰石土，初骨土2紫色土3黄壤土，红壤4水稻土，潮土7、在空间分析模块中，选择重分类“ Reclassify ”进行对象操作，然后选择栅格计算器“ Raster Calculator ”分别输入输入各因素的 权重和贡献函数值(权重表见表7,输入过程见下图7-1)，得到四大 要素的评价分值(各个等级的分值P，P1i=100，P2i =67，P3i =33, P4i =0)，评价分值结果图见下7-2,7-3,7-4,7-5.将不同土地质量数据R的等级，切换为表中的指标值，通过下式计算出栅格 单元

8、的评价分值。R(T)= 100I23，i=1根据使G(P ) = 1-1R(T)-P /1001的值最大时P所对应的J，即为所求的某个T的适宜级别SiQ1Q2Q3Q4167178.08.0244.811.45.45.4322.25.62.62.640000表7影响因素等级评价指标表图7-1：输入公式曰乡Layers56114口 170日 sy_fz0224467乡 Layers日 trfjz055114 170 日 syjzEiO2244W Layers slope_fe TOC o 1-5 h z 2654口囤B 口 trk_fz0255480日 trf_fa056114170h n sv fc图7-4： 土壤孔隙度评价分值图图7-5： 土地坡度评