PAR系统的几点记录
最近在修改PAR系统，好记性不如烂笔头，中间遇到一些问题顺便做个记录，以备后面回顾之需。
1、HDF转换到FLT功能的时候，选择的sample.hdr头文件，里面的内容是：
ncols 1471
nrows 835
xllcorner 65.673578334713
yllcorner 13.996935656021
cellsize 0.05
NODATA_value -9999
byteorder LSBFIRST
这个范围其实是整个中国的范围（外扩一圈），用这个范围裁切下载的全球aerosol数据。
2、以前的测试数据里有地理因子数据，它们的头文件为slope.hdr latitude.hdr aspect.hdr，头文件内容相同，如下所示：
ncols 315
nrows 147
xllcorner 87.833076695378
yllcorner 30.321971925468
cellsize 0.05
NODATA_value -9999
byteorder LSBFIRST
这个范围才是整个三江源的范围（外扩一圈），从中国的dem数据里提取地理因子数据所用的头文件，不用说肯定和这三个头文件一致。
3、全国的dem数据是Krasovsky_1940_Albers投影坐标系的数据，在ARCMap的空mxd文档里，首先加载这个数据是方形的（如下左图）；如果首先加载一个地理坐标系的数据，则后加的这个dem数据就会表现为上宽下窄曲线梯形（如下右图），原因是：投影坐标系上下宽度一样的情况下，顶端比底部跨越的经度范围宽，当在一个“二维的地理坐标系显示窗口”（即直角坐标系的mapcontrol里水平或垂直移动鼠标的时候，获取的水平或垂直坐标值不变）里显示的时候，兴趣区域顶端在mapcontrol里就自然表现出上宽下窄的分布。
4、有一个功能是提取地理因子的时候把所有栅格单元值设置为4，第一次的表达式是：Con( IsNull([R])|(~IsNull([R])),4) ，在arcmap下执行没有错误，可是在arcgis engine下面执行就是出现语法错误，实在不知道为什么。无奈之下就把公式改成如下： Con( IsNull([R]),4,4) ，带入IMapAlgebraOp工具，成功执行。