论坐标系不规范所带来的问题

说好的今天摆烂，结果没摆成。起因是收到师姐的消息，需要处理一篇文章中所提到的栅格数据（HDF 格式）。

收到的问题为在ArcGIS显示异常，且无坐标系定义。对照文档内容：

Time span: 1981-2019

Temporal resolution: daily

Spatial resolution: 0.05°×0.05°

Projection: geographic latitude/longitude projection

Geographical scope: global

Data format：HDF4 (*.hdf)

其中坐标系为地理经纬度坐标系，换句话问：您是哪国人？答：北半球。没有任何参考价值。

将HDF文件输出为ASCII文件来检测坐标系定义，得出的文件头如下所示，额外的小插曲是，转换结果会丢失小数精度：

ncols 7200
nrows 3600
xllcorner -0.5
yllcorner -3599.5
cellsize 1
NODATA_value -999900

ASCII 文件头，那能对么

考虑到并没有相应的参考，因此拖进ArcGIS自然会有问题。观察数据，其为3600*1800的数据及大概形状可知其为全球范围。因此采用WGS 1984定义，修改文件头。将cellsize 修正0.05，将左下脚的坐标修正为（-180，-90），也就是西经180度，南纬90度。

将ASCII文件头修改为如下所示：

ncols 7200
nrows 3600
xllcorner -180
yllcorner -90
cellsize 0.05
NODATA_value -9999

ASCII 修正文件头

然后定义投影为WGS 1984，与底图对照，结果一致。

整理出来的工作流程如下：

1. 将HDF转换为ASC
2. 修正文件头
3. 定义投影
4. 转为ASC
5. 提取数据
6. 输出为CSV

所以说，规范的坐标系定义十分重要。