Borrar filtros
Borrar filtros

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

2 visualizaciones (últimos 30 días)
SChow
SChow el 29 de Nov. de 2022
Comentada: Manoj Mirge el 20 de Mzo. de 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 comentario
Manoj Mirge
Manoj Mirge el 20 de Mzo. de 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.

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Matrix Indexing en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by