ArcGIS栅格计算器con条件函数使用

不需要con函数，比如lyr是添加到ArcMap中的栅格图层，那么可使用类似如下表达式：

("lyr" = 0 ) * 1 + ("lyr" >0) * ("lyr" <0.783) * (0.6508 - 0.3436 * log10("lyr")) + ("lyr" >0.783 ) * 0

没测试，应该可用。其中表达式中("lyr" >0.783 ) * 0是多余的，("lyr" = 0 ) * 1中*1也是多余的，写上去逻辑清楚些。

先用万能地图下载器下载天地图建筑轮廓图片，然后再用ArcGIS对建筑轮廓进行分析提取。

打开水经注万能地图下载器，下载"天地图 WGS84 坐标系"电子地图，注意在新建任务对话框内地图类型选择"仅下载影像"，下载级别选择19级，其余按照默认设置即可。新建任务 回复"下载器"获取万能地图下载器，将下载的天地图电子地图加载到ArcGIS内，点击菜单栏上的"窗口"→"影像分析"。点击影像分析，在显示的"影像分析"对话框内选择上需要提取单波段的图片，点击"添加函数"。在"函数模板编辑器"对话框内，在需要处理的图片上点击右键→"插入"，选择"波段提取函数"。选择波段提取函数，在"栅格函数属性"对话框内，"波段"选择2，"组合"选择2。点击两次"确定"后可以看到处理后的效果图。在ArcToolbox内点击Spatial Analyst工具点击"重分类"→"重分类"。在显示的"重分类"对话框内，"输入栅格"选择提取单波段后的栅格，"重分类"分为三个层级，旧值分别为"100-249"、"250-251"、"252-255"，对应新值为"1"、"2"、"3"。旧值的最小和最大值不用修改，保持和原有的一致。点击"确定"后可以看到重分类的结果。同样在ArcToolbox内点击Spatial Analyst工具点击"重分类"→"重分类"，在"重分类"对话框内"输入栅格"选择之前重分类过后的图片，将旧值"3"的新值设置为"1"，这样分类后就只有两类。点击确定后可以看到再次重分类后的效果。在目录下选择预先准备好的文件夹，点击右键，选择"新建"→"Shapefile"，新建一个Shapefile文件。在显示的"新建Shapefile"对话框内"名称"命名为"折线"，"要素类型"选择"折线"，点击"编辑"，设置坐标系为"WGS 1984"。采取同样的方式，再新建一个面文件。在左侧新建的线图层上点击右键，选择"编辑要素"→"开始编辑"，将线图层变为编辑状态。在菜单栏上点击"自定义"→"工具条"，勾选上ArcScan，显示出来ArcScan工具条。在ArcScan工具条上点击"矢量化"，选择"矢量化设置"。在显示的"矢量化设置"对话框内点击"样式"，选择样式"轮廓"。应用并关闭对话框后在ArcScan工具条上点击"矢量化"，选择"生成要素"。在显示的"生成要素"对话框内一切默认。点击"确定"后会创建矢量，创建后可以看到生成的建筑轮廓。点击编辑器工具条上的"编辑器"→"保存编辑内容"，再点击"停止编辑"完成数据的保存和退出编辑状态。生成的建筑轮廓保存在之前新建的名为"面"的文件内，可以通过加载矢量的方式将其加载到万能地图下载器内，可以看到完美的叠加。

用ArcGIS提取天地图建筑轮廓的详细说明，主要包括了地图下载、波段提取、影像重分类和创建矢量等功能。

---