1 操作软件
ArcMap 10.8(3.1)
ArcGIS Pro(3.2和3.3)
2 操作任务
3 操作步骤
3.1 将影像数据进行分割(即输出每一单波段影像数据)
预处理工作:由于时序人口的栅格数据空间范围远大于待提取的深圳建筑矢量面数据的区域,因此这里可以设置一个掩膜范围,将需要的区域的栅格影像提取出来,简化后续处理步骤(这里笔者实在进行完波段的剥离后才进行的这一步,因此笔者就在后续的统计分析中反复的应用了掩膜提取,这里建议最开始就一步到位,先掩膜提取需要的区域再进行波段剥离)如下图所示
Step1:将时序人口数据中的“holiday.tif”栅格添加到图层中,这里右键查看该栅格数据,可以发现其波段共有24个(如下图),后续所谓添加的24个时段便会对应到这里的24个波段。
Step2:接下来,在图层树中选中该栅格图层,并在主菜单栏中选择【窗口】-【影响分析】,如下图。
Step3:接下来,在【影响分析】的窗口中选中栅格图层(如未选中栅格图层则下面的功能按键都会呈现灰色),然后点击【正在处理】中的【添加函数】按钮,如下图。
Step4:接下来,点击这里的“holiday.tif”栅格图层,右键其,选择“插入函数”,在弹出的一系列函数模板中选择“波段提取函数”,如下图所示。
Step5:接下来,将各波段依次提取即可,注意这里提取相应的波段所对应的【组合】也应为该波段才正确,操作如下图。
如下图为提取的band1的波段影像,就是只有波段1的情况,由此也可反向验证提取的正确性
提取出的所有波段数据都需要导出,这里笔者选择是导出到一个文件地理数据库中,便于后续操作。如下图
3.2 建立Model加快操作步骤
这里不难发现每一个波段的数据的统计都是一个反复的过程,因此这里可以尝试建立model来简化操作步骤,如下图
这里的主题思想就是先【掩膜提取】,然后【转为整型】,然后【栅格转面】,最后进行【空间连接】,这里的连接匹配选项可以设置为“最大重叠”
依次model不断输出即可得到24个波段的统计结果,这里每个字段的统计值应该都是gridcode字段才正确,如下图所示
针对该情况,显然是无法连接的(因为每个波段统计的图层字段都是这个),因此这里需要将该字段重命名。这里笔者提供的一个行之有效的方法就是将相应的gridcode重命名为holiday_n(n即为对应的波段号)
这里使用【计算字段工具】,如下图
计算完即可在对应的图层属性表中出现你想要的字段,如下图
重复此操作,将每个统计完的图层都连接一下
3.3 将波段统计图层中的时序人口数据连接到深圳建筑图层
这里可以观察一下两个图层,确定匹配的连接字段。如下图
不难发现匹配的字段就是这两个,不过总是依次相差1(从FID匹配OBJECTID)
这里笔者提供的解决方法就是在深圳建筑图层(ShenZhenBuild)再新建一个字段为objectID,这样就一致了,这也不麻烦,使用【计算字段】工具,如下图
计算完即可得到相应的需要的连接字段
接下来就是使用【连接字段】工具,将各个波段统计图层中的时序人口数据连接到深圳建筑图层,如下图
将每个字段都传输完即可得到统计结果,在深圳建筑图层中,如下图
4 附录
这里主要包括了
1.设置的掩膜区域(3.1中已描述)
2.模型.py(model的脚本)
3.ShenZhenBuild_holiday(最终操作的时序人口统计后的shp文件)