采用高分辨率DEM获取周围地物遮挡情况3月24日.docVIP

采用高分辨率DEM获取周围地物遮挡情况 2014年3月24日 概述 雷达站周边几十km范围内的地形对于遮挡计算至关重要。而且，离雷达较远的上百千米或更大范围的地形资料，对评估雷达的探测能力，特别是识别大气波导现象也十分重要。 因此，本文对利用高分辨率DEM获取周围地物遮挡情况，进行了研究。采用的地形高度资料是SRTM3（90m分辨率）和ASTERM V2（30m分辨率）。 SRTM3（90m，相当于3秒的分辨率） SRTM的全称是Shuttle Radar Topography Mission，即航天飞机雷达地形测绘计划。2000年2月11日上午11时44分，美国“奋进”号航天飞机在佛罗里达州卡那维拉尔角的航天发射中心发射升空，“奋进”号上搭载的SRTM系统共计进行了222小时23分钟的数据采集工作，获取北纬60度至南纬56度之间，面积超过1.19亿平方公里的 9.8万亿字节的雷达影像数据，覆盖全球陆地表面的80%以上，该计划共耗资3.64亿美元，该测量数据覆盖中国全境。 SRTM系统获取的雷达影像的数据量约9.8万亿字节，经过两年多的数据处理，制成了数字地形高程模型（DEM），即现在的SRTM地形产品数据。此数据产品2003年开始公开发布，经历多修订，目前的数据修订版本为V4.1版本。该版本由CIAT（国际热带农业中心）利用新的插值算法得到的SRTM地形数据，此方法更好的填补了SRTM 90的数据空洞。SRTM地形数据按精度可以分为SRTM1和SRTM3，分别对应的分辨率精度为30米和90米数据（目前公开数据为90米分辨率的数据）。其数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域，达到了地球陆地表面的99%。Aster Gdem数据参数覆盖范围：全球空间分辨率：1弧度秒 （约30 米）精度：垂直精度20米，水平精度30米目前共有2版，第一版V1于2009年公布，第二版V2于2011年10月公布Fun_Select_DEM_Data.m。 注意：该程序包括了将不同的DEM文件的数据进行拼接的功能 程序流程如下： 先按照0、90、180、270度，计算四个方向上，经纬度的边界 如果使用30m分辨率的ASTER GDEM的DEM数据，则采用 DEM_File_Directory_30m 目录下的数据 如果使用90m分辨率的SRTM3的DEM数据，则采用 DEM_File_Directory_90m 目录下的数据，90m分辨率的DEM文件名的命令规则比较奇怪，见网上的介绍。 再根据四个边界，选择哪几个DEM文件 然后再将这些文件中，需要的数据拷贝到以雷达经纬度为中心的数组之中去。 生成的几张DEM地图如下： 中心点经度=118.78,纬度=31.935,量程=10km（基于ASTER GDEM生成） 中心点经度=118.78,纬度=31.935,量程=10km （基于SRTM3生成） 从上面两张图对比可以看出，30m分辨率的就是好。 计算遮挡角的方法 程序见Generate_Radar_Shield_Sector.m。 主要计算过程如下： 根据用户设定的经纬度和量程，选择需要的DEM数据 得到对应DEM数据中各个点的经纬度 调用LonToPolar，计算各个点相对于雷达的仰角 用imagesc绘制各个点相对于雷达的仰角 以横轴是方位（以0.1度为分辨率），竖轴是距离，用imagesc进行显示 计算各个方位上的最大的仰角，然后绘制极坐标下的遮挡图 注意：在设定雷达高度的时候，不能直接利用GPS得到的结果，而是要从DEM数据中，找出雷达所在经纬度的那个数据点，看看这个数据点所表示的高度是多少。和GPS对比，看两者相差多少。如果相差很多，则说明这个点的GPS值是不是搞错了。此时要利用Google Earth来进行确认。 仙林X波段雷达的遮挡计算结果 下面计算仙林雷达（架设高度为60m）的遮挡情况。 仙林雷达30km量程下的DEM数据 仙林雷达四周30km范围内，各个地物相对于雷达的仰角 仙林雷达在不同方位上，各个距离点的地物相对于雷达的仰角 仙林雷达的遮挡角（直角坐标显示） 仙林雷达的遮挡角（极坐标显示） 上面的结果是利用30m分辨率的ASTER GDEM数据计算的。 下面用90m分辨率的SRTM3数据进行计算，最终结果如下： 对比30m和90m的结果，可以发现： 利用30m的DEM数据计算出的结果，在45度的方位上遮挡比90m分辨率下的计算结果更大。但由于这个方向上建了天文台，的确应该更大。 因此，30m的DEM数据更加准确。 溧水S波段雷达的遮挡计算结果 可见，这部雷达四周基本没有遮挡。 但是，C波段雷达架设在该处，发现四周有很多的遮挡。这就说明需要和采用雷达信噪比数据获取周围遮挡情况的技术相互结合，才能获取更加准确的遮挡