关于风云4A 中国区域的坐标

schliezer

一阵妖风 发表于 2019-9-27 16:44
请问raw文件是怎么读的啊，我用python和matlab读取之后结果不一样是怎么回事啊

你好，解决了吗？我也是这种，两个读出来不一样

安达充

schliezer 发表于 2020-1-14 10:58
你好，我也遇到了类似的问题。我用的风云4号16km分辨率、区域合成的大气垂直廓线产品AVP。读16000.raw的结 ...

参考上面已经解决的方法。

风云卫星官方提供的经纬度对照表是全圆盘的，包含全球。
而你下载的数据是区域的。

你需要找到类似上面的同学找到的参数。

安达充

windprofiler 发表于 2019-10-11 01:53
您好，我需要全面盘4km的对应坐标，找了很久没有找到，能方便发我一份或者链接吗？

在风云卫星的官网可以下载。

链接：
http://satellite.nsmc.org.cn/Por ... nload.aspx?TypeID=2

schliezer

安达充 发表于 2020-1-14 15:11
参考上面已经解决的方法。

风云卫星官方提供的经纬度对照表是全圆盘的，包含全球。

我用matlab读了同样的16000.raw，数组长度变成了943938。后来我找到了16000.raw的说明文档，行列数均为687，也就是总的数组长度为687*687=471969，小于943938或者7551504，全圆盘还小于中国区域？而且文档里面没有中国区域在这个圆盘中的起止点或者数组长度。

schliezer

安达充 发表于 2020-1-14 15:11
参考上面已经解决的方法。

风云卫星官方提供的经纬度对照表是全圆盘的，包含全球。

另外，我按照说明书里提示的方法，读了16000.raw的经纬度，如果latlon1是经度，latlon2是纬度，那第三个-124的纬度也不正常

安达充

schliezer 发表于 2020-1-14 16:31
另外，我按照说明书里提示的方法，读了16000.raw的经纬度，如果latlon1是经度，latlon2是纬度，那第三个- ...

请认真阅读你自己发的截图内文件格式的全部内容。

每个网格对应16字节，这16字节内，前8为经度，后8为纬度。

所以第三个-124就是第二个网格的经度，这个值没有不正常的地方。

安达充

schliezer 发表于 2020-1-14 16:22
我用matlab读了同样的16000.raw，数组长度变成了943938。后来我找到了16000.raw的说明文档，行列数均为68 ...

687*687*2（经度&纬度，两个要素）=943,938

区域的行列数请在数据文件里找到那几个line number相关的参数。

schliezer

安达充 发表于 2020-1-14 16:39
687*687*2（经度&纬度，两个要素）=943,938

区域的行列数请在数据文件里找到那几个line number相关的 ...

恕我愚钝，并没有抬杠的意思，但我确实没有理解你说的“每个网格对应16字节，这16字节内，前8为经度，后8为纬度。

所以第三个-124就是第二个网格的经度，这个值没有不正常的地方。” 我截取前100个数据，结果全是这几个数，完全没有看出坐标的规律性变化，而且，我理解的是，从第一个开始，每隔一个是经度，从第二个开始，每隔一个是纬度，也就是经度：-113，-14，127,46……，纬度是：-28，-1，-124,46……，这样确实显然不行，如果按你说的-124是经度，那就应该是-113，-28是经度，-14，-1是纬度然后以此类推吗？我后来看到了你说的687*687*2=943938的帖子，难道是只有用matlab读才能得到正确结果？多谢解答！

安达充

schliezer 发表于 2020-1-14 16:54
恕我愚钝，并没有抬杠的意思，但我确实没有理解你说的“每个网格对应16字节，这16字节内，前8为经度，后8 ...

经纬度的数据类型不是整型

schliezer

安达充 发表于 2020-1-14 17:00
经纬度的数据类型不是整型

谢谢！我试试matlab先