Combine all MODIS 1km tiles for one day and make a global matrix

Hi,

I am trying to combine all MODIS 1km tiles (there are 35*17 sinusoidal grid tiles each of size 1200x1200, see https://modis-land.gsfc.nasa.gov/MODLAND_grid.html) into a lin- spaced grid file of 1km for one day and ultimately for an year.

I am using griddata. I am running the following lines within a loop. But the griddata function to produce 1km global data for one day is taking more than 3 hours to process. I am running on a server which is very fast.

x=lon(:);y=lat(:); %where lon and lat are of size 1200x1200 each and are sinusuidal longitude and latitude of one of the 35*17 tiles. 
d=double(Data(:)); % where Data is the data extracted from the tile corresponding to lat and lon

[LAT LON]=cdtgrid(360/43200); %Global 1km mesh of lat and lon. Creaded using cdt tool box https://www.chadagreene.com/CDT/CDT_Getting_Started.html

v=griddata(x,y,LON,LAT) %% This is taking more than 3 hours to process.

Any alternatives will be very much appreciated.

[Manoj Mirge]

Hi SChow,

griddata function in MATLAB accepts minimum of 5 and maximum of 9 arguments. In your case you are giving only 4 arguments to griddata function which would generate error. Can you clarify agian wheather you are getting error from server or you are getting desired output but after very long processign time.