SatyrLee
文章79
标签21
分类7

一言

文章归档

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

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

说好的今天摆烂,结果没摆成。起因是收到师姐的消息,需要处理一篇文章中所提到的栅格数据(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

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

本文作者:SatyrLee
本文链接:http://www.naive514.top/posts/671a8b9a/
版权声明:本文采用 CC BY-NC-SA 4.0 CN 协议进行许可