关于HDF文件提取数据

没错我是猫

由于研究需要由Color Ocean上下载了海水表面温度卫星数据，格式为HDF，使用HDFviewer可以打开，但是由于从来没有使用过该类软件，也不知道如何才能从下载的卫星数据中提取出我需要的经纬度处的海水表面纬度，有没有大神可以为我指点一二啊？
附上下载的卫星数据，精度显示为1km，求大神教教我如何提出某一特定地理位置的数据。

xzshzzh918

那个HDFviewer我只是用来查看HDF文件的，好像不能提取数据（或者能提取我不知道……），我用过matlab和python提取HDF数据的，都很简单。

MeteoInfo

可以用MeteoInfo脚本来提取：

1. fn = 'C:/Temp/A2014357.L3m_DAY_MZBXTNSST_sst_MZB_1km'
   
2. f = addfile(fn)
   
3. v = f['l3m_data']
   
4. data = v[:,:]
   
5. slope = v.attrvalue('Slope')[0]
   
6. intercept = v.attrvalue('Intercept')[0]
   
7. data = data * slope + intercept
   
8. lon = arange(31, 40, 0.010416667)
   
9. lat = arange(-30, -20, 0.010416667)
   
10. data.setdimvalue(0, lat)
    
11. data.setdimvalue(1, lon)
    
12. #Interpolate to station
    
13. x = 36.2
    
14. y = -23.1
    
15. v = data.tostation(x, y)
    
16. print v
    
17. #Plot
    
18. axesm()
    
19. lworld = shaperead('D:/Temp/map/country1.shp')
    
20. geoshow(lworld)
    
21. layer = imshowm(data, 20)
    
22. colorbar(layer)

fn = 'C:/Temp/A2014357.L3m_DAY_MZBXTNSST_sst_MZB_1km'<br /> f = addfile(fn)<br /> v = f['l3m_data']<br /> data = v[:,:]<br /> slope = v.attrvalue('Slope')[0]<br /> intercept = v.attrvalue('Intercept')[0]<br /> data = data * slope + intercept<br /> lon = arange(31, 40, 0.010416667)<br /> lat = arange(-30, -20, 0.010416667)<br /> data.setdimvalue(0, lat)<br /> data.setdimvalue(1, lon)<br /> #Interpolate to station<br /> x = 36.2<br /> y = -23.1<br /> v = data.tostation(x, y)<br /> print v<br /> #Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld)<br /> layer = imshowm(data, 20)<br /> colorbar(layer)

why123

我只会用C++提取，有代码

devcry

hdf文件最外面一层是有一个介绍信息的，包括数据的种类，数据的属性等等。你可以找到经纬度一栏来定位行列，用这个行列去海表温度一栏找所需信息。
以上