用EASE2 720x720 grid tool 读取 osisaf-450 432x432 数据

参考资料
①https://github.com/alan-turing-institute/IceNet-Project/issues/5
②http://psc.apl.washington.edu/nonwp_projects/PPobs//EASE_GRID_TOOLS.html
③EASE官网：https://nsidc.org/data/ease

数据：OSISAF-450海冰密集度产品
投影：EASE2 grid( https://nsidc.org/data/ease )
这里面grid tool提供的是720x720

ftool = 'F:/osi450/gridloc.EASE2_S25km/'
lats = np.fromfile('F:/osi450/gridloc.EASE2_S25km/EASE2_S25km.lats.720x720x1.double',
                   dtype=np.float64).reshape((720,720))
lons = np.fromfile('F:/osi450/gridloc.EASE2_S25km/EASE2_S25km.lons.720x720x1.double',
                  dtype=np.float64).reshape((720,720))

但osi-450数据是432x432，根据参考资料①②，432x432在720x720的位置居中，

对于南半球来说
下面这张图表示的是OSI-450数据覆盖范围432x432, 最多到南纬16°。

lats_osi = lats[144:576,144:576]
lons_osi = lons[144:576,144:576]