城市建筑日照分析

现有一城市某区域建筑数据和地块单元数据。计算该地区各个单元的容积率。容积率=建筑面积/用地面积（建筑面积为各楼层建筑面积之和，用地面积为各地块单元parcel.shp的面积）。假设建筑物A与建筑物B在空间上存在阴影遮挡(即A遮挡了B)，则A在B向阳方向的前方。B的房顶会遮盖住A的阴影，给遮蔽判断带来困难，通过分析阴影与建筑物的空间叠加关系，找出不符合日照标准的建筑物。

打开工具箱中的分析工具，点击叠加分析中的标识，如下图所示。

标识完后打开buildings_identity的属性表，在ID字段右键选择汇总，选择f_area中的总和确认输出，输出后，将parce与Sum_Output连接，如下图所示。

连接后打开parce的属性表，新建双精度字段rate并计算容积率（一个单元内建筑的总面积除以地块单元的面积）如下图所示。

将buildings进行转栅格，点击工具转为栅格中的面转栅格并输出如下图所示。

将输出的buildings_PolygonToRaster进行重分类，将新值里的NoData设置为0，其他保持和旧值保持不变，点击确定如下图所示。

对坡向进行处理，点击坡向工具，如下图所示。

找到坡向后进行栅格计算器，分别把12点13点14点的阴影栅格给求出来。角度和公式方位角可参考资料，如下图所示。

选择之前重分类的Reclass_buil1数据乘以back12，back13，back14得到建筑物的阴影的高度。打开栅格计算器依此类推，并输出为dem12，dem13，dem14，如下图所示。

完后进行山体阴影的处理，选择工具里的山体阴影，如下图所示。

进行重分类，分为2类，中断值第一个改为0，将0值赋为1，其他值赋为0，完后输出，以此类推输出12.13.14。打开栅格计算器将输出的Reclass_Hill12，13，14相加并输出，如图所示。

点击重分类将新值改为0.1.1.1输出，如下图所示。

打开刚输出的Reclass_rast1的属性表选择阴影部分，点击工具中栅格转面工具并输出，最终得到会有阴影的部分，如下图所示。

选择不合适的建筑，点击选择中的按位置选择，从buildings中选择，目标图层要素的空间选择方法选择为质心在源图层要素内，应用确定，最终选中了不合格的数据，将不合格的数据导出，如下图所示。

最终导出不合格的数据，如下图所示。