提取城市土地透水不透水信息到城市雨水栅格存储方法.docxVIP

透水不透水信息获取在划分城市雨水网格并形成面矢量图层后，图层中仅仅是的要素仅仅包含基于空间位置的信息，并没有各个栅格对应的范围的属性信息（如土地类型，透水不透水面积，人口等等），为了对各个栅格进行雨水积水数值模拟以及灾害损失评估计算，需要对上述各类属性信息进行采集并以城市雨水网格为单位进行存储，下面以提取城市土地透水不透水信息为例介绍城市雨水网格信息采集方法。利用ERDAS软件结合城市TM数据可以对城市土地利用类型进行分类，根据分类的结果将土里利用类型分为透水和不透水两种类型，并输出为矢量图层，将分类图层与城市雨水网格图层进行叠加分析，可以得到每个城市雨水网格中透水区域和不透水区域的子网格块，分别统计各个子网格块的面积并作为字段存储于城市雨水网格图层的属性表中即可完成透水不透水信息的采集。下面详细介绍具体的采集步骤：获取城市雨水网格区域内的TM遥感影像数据通过中科院网站获取武汉区域范围内的TM遥感影像数据，年份为2005年，利用ERDAS软件的Main-Interpreter-Utilities-Layer Stack模块，取1(R),4(G),5(B)三波段进行合成,并输出为IMG格式的文件（145.img），如下图所示：由于TM数据包含的局域很大，下面使用在城市雨水网格划分中曾使用过的POLYGON图层对TM数据进行掩膜处理，通过ARCGIS工具箱中的Spatial Analyst Tool-Extraction-Extract by Mask模块截取城市雨水网格试验区范围内的TM数据，并输出为img格式的文件（hust_1451.img）二、图像预处理使用ARCGIS中的Georeferencing工具结合在城市雨水网格划分中使用的路网数据对hust_1451.img进行空间几何校正，得到校正后的结果：三、进行监督分类首先确定分类模板，利用Classifier-Signature Editor来打开分类模板界面，并利用AOI获取工具来分别选取遥感影像中的特征区域，并设定其对应的透水不透水类型，首先设定AOI工具的参数，在视窗界面中AOI-Seed Properites界面，作出如下设置：接下来在地图中选取特征区域并在分类模板中选择按钮进行添加，并在分类模板中设定其透水不透水类型，在TM区域内，绿色和浅绿对应植被和土壤区域，为透水区，而紫红，白色，浅灰黄区域内的土地对应水泥铺砌，建筑道路等不透水区域，根据此进行分类模板的设置，如图：可以将分类模板导出为文件进行存储（145.sig）。确定了分类模板后可以对其进行评价，选择分类模板界面的Evaluate-Contingency功能进入可行性矩阵设置界面，进行如下设置：点击ok按钮可以得到可行性矩阵，对矩阵进行观察，判断分类的准确性，如果准确性小于85%则说明分类模型不够精确，需要修改：确定其可行性后，便可以开始进行监督分类，在Classifier-Supervised Classification打开分类界面设定如下：点击确定后进行分类，输出分类结果为classed_145.img，如下图四、数据采集预处理现在可以在ARCGIS里将分类完成后的栅格数据打开，首先观察栅格图像的属性数据，如图：很明显图像周围存在黑边，使用ARCGIS工具箱中的Data Management Tools-Raster-Raster Dataset-Copu Raster工具进行黑边去除处理，参数设置如下：去除黑边后，使用工具箱内的Conversion Tools-From Raster-Raster to Polygon方法将栅格图像转化为矢量面图层，注意Field参数应选择class_name，这样才能以透水不透水进行统计，如下图输出的数据如下图所示：显然这样的分类过于复杂，不易于进行叠加操作，因此我们需要对分类结果进行处理，去除图上的分类斑点，这里需要使用ERDAS的CLUMP和EILIMATE工具对分类结果进行加工处理，首先使用Main-Image Interpreter-GIS Analysis-Clump打开Clump界面，进行如下设置来去除栅格图上的斑点输出图如下：接下来使用Main-Image Interpreter-GIS Analysis-Eliminate打开Eliminate界面进行聚类，将相同类型的栅格块集合在一起，进行如下设置：处理完成后得到的结果在ARCGIS打开如下：由于进行聚合处理后的各个栅格属性信息没有的透水不透水字段，现在进行手工填入，对照原始分类图添加透水不透水信息，完成后的属性信息表如下图：下面使用之前的方法将栅格图像转换为矢量面涂层，注意在转换时栅格图层属性字段不能有中文，否则则会导出出错，其中FIELD字段设为class，得到的结果如下：其属性表如下：现