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

7 visualizzazioni (ultimi 30 giorni)
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.
  1 Commento
Manoj Mirge
Manoj Mirge il 20 Mar 2023
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.

Accedi per commentare.

Risposte (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by