纬度偏移

[timeroute]

我下载了GFS分辨率1°的数据，使用GDAL处理成 EPSG:3857 投影的 PNG 灰度图，渲染出来的风场与实际风场的纬度有偏差，请看一下什么原因

风场图片

uv最大小值

{"umin": -21.246559639451306, "umax": 25.762312227473963, "vmin": -23.16313204755423, "vmax": 23.55504539323816}

与nullschool的对比如下（左是我的，右是nullschool的）

[sakitam-fdd]

@timeroute 不太确定，你底图是 WGS84 和 GCJ02 的

[timeroute]

底图用的mapbox的卫星瓦片，叠加了海岸线，能匹配上，应该没问题

[timeroute]

图片预处理流程如下：

1. 读取二维矩阵数据，设置 transform = -180.0, -85.051129, 180.0, 85.051129，crs = 4326
2. 重投影到 crs = 3857
3. 获取最大最小值和归一化处理
4. 将数据写到 r、g 通道

[echohye]

图片预处理流程如下：

1. 读取二维矩阵数据，设置 transform = -180.0, -85.051129, 180.0, 85.051129，crs = 4326
2. 重投影到 crs = 3857
3. 获取最大最小值和归一化处理
4. 将数据写到 r、g 通道

@timeroute 请问这个问题解决了吗

[timeroute]

图片预处理流程如下：

1. 读取二维矩阵数据，设置 transform = -180.0, -85.051129, 180.0, 85.051129，crs = 4326
2. 重投影到 crs = 3857
3. 获取最大最小值和归一化处理
4. 将数据写到 r、g 通道

@timeroute 请问这个问题解决了吗

没有

[echohye]

图片预处理流程如下：

1. 读取二维矩阵数据，设置 transform = -180.0, -85.051129, 180.0, 85.051129，crs = 4326
2. 重投影到 crs = 3857
3. 获取最大最小值和归一化处理
4. 将数据写到 r、g 通道

@timeroute 请问这个问题解决了吗

没有

请问这个图片预处理是怎么设置crs再重投影到3857的

[sakitam-fdd]

@timeroute 我最近研究一下，之前没注意这个问题，不知道是不是跟 grib 的数据组织有关系，就是 1°的数据其实他的范围不是-180-180，好像有半个格子的偏移