read a grid hdf file with Matlab code

[nastakpa84]

Hello,
I have downloaded the Matlab code for a grid hdf file of SeaWifs of OBPG but it does not work.
I tried to read*chlorophyll data, that i have download from Ocean Color page but the code does not work.
Can you help me?

Anastasia.

[hyoklee]

Hi,

The data file (S20000732000080.L3m_8D_CHL_chlor_a_9km.G3.hdf) has all 0 values in the dataset.

Also, would you please post the error message you got and tell us which version of MATLAB you use?

[nastakpa84]

Hi,
I am using the Matlab7.9.0(R2009b) and the error message is:
Error using ==> convertToGeoRasterRef at 49
Function MESHGRAT expected its second input argument, R, to be either a 3-element referencing vector or a 3-by-2
referencing matrix.

Error in ==> meshgrat at 130
R = convertToGeoRasterRef(R, size(Z), 'degrees', mfilename, 'R', 2);

Error in ==> surfacem>parseInputs at 114
[lat,lon] = meshgrat(lat,lon);

Error in ==> surfacem at 77
[lat, lon, alt, grid, pvPairs] = parseInputs(varargin{:});

Error in ==> surfm at 85
surfacem(lat,lon,Z,varargin{:});

Error in ==> read_hdf_seawifs_grid_file_OBPG_initial_run at 127
surfm(lat,lon,data);


Anastasia

[johngevans]

What happens if you try

pcolor(lon,lat,data); shading flat

[nastakpa84]

Again an error has occured.
Error using ==> pcolor at 58
Matrix dimensions must agree.
Maybe the code cannot read the data from GIOVANNI page.
If i download the chlorophyll data from the ocean color page and exactly from the Level 3 Brownser, I get a file with different name (S20000812000088.L3m_8D_CHL_chlor_a_9km.hdf), and the Matlab code can read the chlorophyll data.
I cannot understand what is happened.

Thank you.

[johngevans]

Send me email at john.evans@mathworks.com and we will try to work it out with you.