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

schliezer

墨青衫 发表于 2019-8-19 13:29
非常感谢您，我已经找到参数了。
      Begin Line Number : 179
      Begin Pixel Number : 0 ...

你好，这个参数要在哪里找呢？我用的风云4A区域合成大气廓线产品里找不到这些

ajiao7370

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

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

你好，请问一下，如果用fortran 读取raw文件，数组如果定义的话是 integer*8 data(2748,2748,2)？ 数据类型是 byte？ double？ 试了几个发现读出来都不对。请大神指教。我用的是 4km分辨率的，所以是2748

安达充

ajiao7370 发表于 2020-3-13 10:51
你好，请问一下，如果用fortran 读取raw文件，数组如果定义的话是 integer*8 data(2748,2748,2)？ 数据类 ...

首先，请先确定你的数据是整型还是浮点型，用整型去读浮点型的数据，解析出来必然是错的；

第二，integer*8代表一个数据长度为8字节。至于是用*8还是*4取决于数据文件说明文档。

ajiao7370

安达充 发表于 2020-3-13 17:23
首先，请先确定你的数据是整型还是浮点型，用整型去读浮点型的数据，解析出来必然是错的；

第二，inte ...

谢谢，应该是浮点型。用fortran 数组 real*8   geo（2,2748,2748） 就可以读出来了。

飞奔的蜗牛

请问楼主，怎么解决的呀？求教

renjifei_2005

学习了。感谢。

锦叶seven

schliezer 发表于 2020-1-14 17:04
谢谢！我试试matlab先

你好，我也遇到和你一样的问题，请问你后面用python读出来了么

schliezer

锦叶seven 发表于 2021-1-5 10:53
你好，我也遇到和你一样的问题，请问你后面用python读出来了么

后来没管那些参数了，直接：
    f=Dataset(ncpath,'r')
    lon_aim=106.47
    lat_aim=29.58
    lev=np.array(f.variables['Pressure'])
    lon=np.array(f.variables['MW_Longitude'])
    lat=np.array(f.variables['MW_Latitude'])
    flag_lat=(28<=lat)&(lat<=33)
    flag_lon=(105<=lon)&(lon<=111)
    flag_lev=(interp_lev[0]<=lev)&(lev<=interp_lev[-1])
    at=np.array(f.variables['AT_Prof'])
    at=at[flag_lev,:][:,flag_lat&flag_lon]