arcgis 特别问题纪要

1. 10.6 以下版本，出现 投影栅格(ProjectRaster)报错：ERROR 000151: 地理变换无效。
   
   是官方未修改问题，10.6已经解决

2. 矢量操作工具：

   1. 使用合并工具：【Data Management Tools】–> 【General】–>【Merge】，将三个要素合并成同一个要素；
   2. 使用矢量编辑-union工具：实现同一矢量图层中多个要素union合并，也可以实现不同矢量图层的要素union操作；
   3. 使用溶解工具：【Data Management Tools】–> 【Generalization】–>【Dissolve】，将上述Merge操作的结果或者Union操作的结果溶解为一个连续的要素；
   4. 使用分析工具-overlay-erase：实现两个矢量图层的差集布尔运算；

3. 带洞面要素创建流程：

   1. 先创建完整外轮廓面要素a.shp；
   2. 再创建一个洞口区域面要素b.shp；
   3. 设置a.shp为编辑状态，选中b.shp的某个洞口要素，点集editor菜单项中clip工具裁剪即可，设置如下：
      

4. 小范围栅格像元值 替换 大范围栅格数据中对应范围的像元值流程
   小范围栅格a.tif 大范围栅格b.tif
   使用【栅格计算器】工具 ArcToolbox——Spatial Analyst Tools——Map Algebra——Raster Calculator
   1）设置标记值
   以b.tif作为轮廓约束，a.tif范围内的标记为0，范围外设置成1
   公式：Con(IsNull(a.tif),1,0)
   arcgis 环境参数中，处理轮廓参数按下图所示设置
   
   2）利用1)步得到的结果c.tif，和b.tif做乘法处理
   公式：“c.tif” * “b.tif”
   处理轮廓参数依然按照1）步所示设置
   3）利用2)步得到的结果d.tif，此时，在d.tif文件中a.tif范围内的值都已经标记为0，判断标记获取最终像元值
   公式：Con(“d.tif” == 0, “a.tif”, “b.tif”)
   处理轮廓参数依然按照1）步所示设置
   得到的结果就是最终替换结果

注意：栅格计算器 计算的数据最好保持分辨率、波段数一致，否则生成结果异常

5. 使用实测数据替换dem部分高程值
   1）使用dem提取等高线dgx.shp
   【3D Analyst Tools】–>【raster surface】–>【contour】
   2）根据实测纵断中线，进行缓冲区分析，生成缓冲区面buffer.shp
   【Analysis Tools】–>【Proximity】–>【Buffer】
   3）面buffer.shp转成线buffer_line.shp
   【Data Management Tools】–>【Features】–>【Polygon To Line】
   4）根据 buffer.shp + dgx.shp，切除buffer区域内的等高线，得到dgx_clip.shp
   【Analysis Tools】–>【Overlay】–>【Erase】
   5）根据 buffer_line.shp + dgx_clip.shp + 实测数据，生成tin

6. 提取栅格波段信息
   【Data Management Tools】–>【Raster】–>【Raster Processing】–>【Composite Bands】
   
   从文件夹选择栅格数据的波段即可提取所要的信息

7. arcpy批处理流程
   7.1、打开ArcCatalog，继续打开操作的文件夹目录树(比如I:/test)；
   7.2、test文件夹右键，选择菜单“New–>ToolBox”，创建工具箱(比如test.tbx)；
   7.3、test.tbx工具箱右键，选择菜单“Add–>Script”，弹出创建脚本配置对话框，根据提示一步步配置，如下图所示
   
   
   
   7.4、右键刚刚创建好的脚本，点击“Edit”，可以打开脚本编写arcpy，如下所示
   import arcpy

# 根据顺序下标获取对应的参数
dir = arcpy.GetParameterAsText(0)
shp =arcpy.GetParameterAsText(1)
outshp =arcpy.GetParameterAsText(2)
spatialRef = arcpy.GetParameter(3)

# 矢量坐标系转换
transform_method=“wgs84_cgcs2000_null” // 现在arcgis里面自己创建地理变换
arcpy.Project_management(shp, outshp, spatialRef , transform_method)

7.5、编码结束后，右键刚刚创建好的脚本，点击“OpenCV”，就可以开始运行脚本

8. arcmap矢量图斑半透明化与影像融合
   8.1 创建图斑，添加颜色属性，如下图所示
   
   8.2 图斑转tif，每一个颜色转为一个tif，如下图所示。
   Polygon to Raster
   
   最后把四个颜色分量的tif都生成出来
   
   8.3 颜色分量合并成一个RGBA四波段的tif
   Composite Bands
   
   8.4 图斑四波段栅格图片与影像融合
   
   最后结果如下所示：
   
9. RGB三波段影像扩展透明通道变为RGBA四波段影像
   9.1 提取RGB原始影像的矩形轮廓，两种方法：
   法一：3D Analyst Tools --> Conversion --> From Raster --> Raster Domain
   这种方法不能存在无效值才能提取矩形轮廓
   法二：直接创建矩形范围矢量面

9.2 给矢量轮廓面附加透明通道属性

9.3 矢量轮廓面转tif，透明通道转为一个tif，如下图所示。
Polygon to Raster

9.4 合成波段，原始RGB三波段影像与上述提取的透明通道tif合成RGBA四波段影像
Composite Bands

10. 矢量属性表自动赋予随机值
    从属性表上打开“字段计算器”，按照下图设置代码即可
    
11. 多波段影像的【栅格计算器】使用注意
    【栅格计算器】针对单波段使用，所以对于多波段影像，应该单独加载某个波段用于栅格计算器计算。
    比如针对一个RGBA四波段的a.tif影像使用【栅格计算器】完成波段1的二值化处理：

• 首先加载a.tif - band1
• 针对波段1的二值化处理的【栅格计算器】处理表达式：

Con("a.tif - band1" > 128,Con("a.tif - band1" < 200, 1, 0), 0)