Filter löschen
Filter löschen

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

4 Ansichten (letzte 30 Tage)
SChow
SChow am 29 Nov. 2022
Kommentiert: Manoj Mirge am 20 Mär. 2023
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 Kommentar
Manoj Mirge
Manoj Mirge am 20 Mär. 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.

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Kategorien

Mehr zu Matrix Indexing finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by