投影坐标是否准确？

[Mo-Dabao]

AwxReader/awx/__init__.py

Line 279 in 4b1e4da

x = np.linspace(ll_x, ur_x, h2.width)

我尝试读取Lambert投影的数据，也把投影中心当作图像中心，但是反算出来的经纬度范围和二级头信息中的经纬度范围不太匹配。

推测投影中心并不是图像中心。

NSMC的经纬度对照表格点数是512*512的不确定怎么缩放：http://satellite.nsmc.org.cn/PortalSite/StaticContent/DocumentDownload.aspx?TypeID=22

方便加QQ或者微信讨论下吗？

[Mo-Dabao]

推测投影中心并不是图像中心。

如果投影中心是是图像中心，经度范围应该是关于投影中心经度对称的

[Mo-Dabao]

我用MICAPS4打开了几个AWX文件看了一下，发现MICAPS4只有一套固定投影（"D:\MICAPS4.8.0\config\Map\projections.ini"），但逻辑确实是用投影中心当作图像中心。

而AWX二级头信息中的经纬度范围确实是MICAPS4那一套固定投影产生的范围，和它本身的投影中心并不匹配。

综合文档中说二级头信息中的经纬度范围是大致范围，所以投影中心优先级应当高于经纬度范围。

所以应该可以认为这个思路没问题，不是bug。至少不是咱作为数据用户端解析的bug。

AwxReader/awx/__init__.py

Lines 273 to 280 in 4b1e4da

	cx, cy = transformer.transform(h2.clon / 100., h2.clat / 100.)
	dx, dy = h2.reso_h / 100. * 1000, h2.reso_v / 100. * 1000,
	ll_x = cx - (dx * h2.width / 2.)
	ll_y = cy - (dy * h2.height / 2.)
	ur_x = cx + (dx * h2.width / 2.)
	ur_y = cy + (dy * h2.height / 2.)
	x = np.linspace(ll_x, ur_x, h2.width)
	y = np.linspace(ur_y, ll_y, h2.height)

[wqshen]

我用MICAPS4打开了几个AWX文件看了一下，发现MICAPS4只有一套固定投影（"D:\MICAPS4.8.0\config\Map\projections.ini"），但逻辑确实是用投影中心当作图像中心。

而AWX二级头信息中的经纬度范围确实是MICAPS4那一套固定投影产生的范围，和它本身的投影中心并不匹配。

综合文档中说二级头信息中的经纬度范围是大致范围，所以投影中心优先级应当高于经纬度范围。

所以应该可以认为这个思路没问题，不是bug。至少不是咱作为数据用户端解析的bug。

AwxReader/awx/__init__.py

Lines 273 to 280 in 4b1e4da

	cx, cy = transformer.transform(h2.clon / 100., h2.clat / 100.)
	dx, dy = h2.reso_h / 100. * 1000, h2.reso_v / 100. * 1000,
	ll_x = cx - (dx * h2.width / 2.)
	ll_y = cy - (dy * h2.height / 2.)
	ur_x = cx + (dx * h2.width / 2.)
	ur_y = cy + (dy * h2.height / 2.)
	x = np.linspace(ll_x, ur_x, h2.width)
	y = np.linspace(ur_y, ll_y, h2.height)

实际测试发现，兰伯特投影和墨卡托投影的经纬度确实与实际数据不符，目前正在解决当中。
已有的结果表明，墨卡托投影下lat_ts需指定为0（默认），而不是数据头信息中的标准纬度，才能保持投影正确。
兰伯特投影也存在轻微的位置漂移，如下图，目前正在查找原因。

[sxcgc]

我研究发现，这个图的中心位置应该是对的。可能是投影的另外两个纬度参数有问题，理由如下：
图片中心附近有一水库：龙羊峡水库

水库位置投影后偏移较小
但是宝岛台湾和斯里兰卡的位置偏移方向不一致，应该是两个纬度参数有问题。

[wqshen]

最近工作任务比较重，您有精力仔细分析一下投影参数该怎么调整吗？如果解决了可以直接提交融合，感谢啊

我研究发现，这个图的中心位置应该是对的。可能是投影的另外两个纬度参数有问题，理由如下：
图片中心附近有一水库：龙羊峡水库

水库位置投影后偏移较小
但是宝岛台湾和斯里兰卡的位置偏移方向不一致，应该是两个纬度参数有问题。