风云4号AWX数据读取画图

slddbs

您有没有读取过FY4A nc的云亮温格式，读取变量为2748*1108，但风云卫星FullMask_Grid.raw的为2748*2748，这个如何解决，感谢啦

slddbs

Masterpiece 发表于 2019-3-4 12:58
怪不得...
另外建议楼主学一下pyqt，交互会好很多

您有没有读取过FY4A nc的云亮温格式，读取变量为2748*1108，但风云卫星FullMask_Grid.raw的为2748*2748，这个如何解决，感谢啦

yuyue1314

slddbs 发表于 2019-7-9 17:39
您有没有读取过FY4A nc的云亮温格式，读取变量为2748*1108，但风云卫星FullMask_Grid.raw的为2748*2748，这 ...

裁剪下，因为你读取的是中国区数据而经纬度格点是4km圆盘，所以会少点，没记错应该是切一下南北向的行列号，这个在属性文件里面都有的。

slddbs

yuyue1314 发表于 2019-7-16 17:08
裁剪下，因为你读取的是中国区数据而经纬度格点是4km圆盘，所以会少点，没记错应该是切一下南北向的行列 ...

def fy4disk(rawfile, dim):
    sz = np.fromfile(rawfile, dtype=np.float, count=2*dim*dim)
    print(sz.shape)
    latlon = np.reshape(sz, (dim, dim, 2))
    lonn = latlon[:, :, 0]
    latt = latlon[:, :, 1]

    latt = np.ma.masked_values(latt, 999999.9999)
    lonn = np.ma.masked_values(lonn, 999999.9999)

    return latt,  lonn
请问 ，我这样读取经纬度格点数据是对的吗？总感到出的结果怪怪的

非对称

感谢楼主分享。

成信之子

谢谢楼主分享，学习一下

slddbs

CARTOPY怎么画有numpy赋值缺省的经纬度，TOPOLOGYEXCEPTION side location conflict at 300 300
shapely. errors. Topological Error the operation  GEOSContains could not be performed总是提示这种错误

yuyue1314

slddbs 发表于 2019-8-5 10:15
CARTOPY怎么画有numpy赋值缺省的经纬度，TOPOLOGYEXCEPTION side location conflict at 300 300
shapely. ...

matplotlib现在不允许经纬度内有nan或者mask

HSY

为啥我们下载不到AWX的格式。。。

yuyue1314

HSY 发表于 2019-9-3 09:12
为啥我们下载不到AWX的格式。。。

这是micaps用的格式，在micaps里下或者在mdfs里下