X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 14 Oct 2016 05:58 PM by  anon
MODIS Conversion Toolkit - IDL Execution Error for MOD11A1 products
 0 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
14 Oct 2016 05:58 PM
    Dear All, I running this code but i have this message error: % Keyword GRID not allowed in call to: CONVERT_MODIS_DATA % Execution halted at: REPROJECT_MOD11A1 57 C:\Users\User\IDLWorkspace80\Default\reproject_mod11a1.% $MAIN$ Here is my code. ;batch det LST_day_time from MOD11A1 PRO reproject_mod11a1 COMPILE_OPT IDL2 ENVI, /RESTORE_BASE_SAVE_FILES ENVI_BATCH_INIT; ;choice your data floder input_directory=dialog_pickfile(title='chioce data floder',/directory) output_directory=input_directory IF ~STRCMP(STRMID(input_directory, 0, 1, /REVERSE_OFFSET), '\') THEN BEGIN input_directory = input_directory + '\' ENDIF IF ~STRCMP(STRMID(output_directory, 0, 1, /REVERSE_OFFSET), '\') THEN BEGIN output_directory = output_directory + '\' ENDIF input_filenames = FILE_SEARCH(input_directory + 'MOD11A1' + '*.hdf') IF SIZE(input_filenames, /N_ELEMENT) GE 1 THEN BEGIN ;Output method schema is: ;0 = Standard, 1 = Projected, 2 = Standard and Projected out_method = 2 output_projection = envi_proj_create(/geographic) out_ps_x = 0.05d out_ps_y = 0.05d interpolation_method = 6 ; swath_name = 'mod04' ; sd_names = ['Image_Optical_Depth_Land_And_Ocean'] grid_name = 'MODIS_Grid_Daily_1km_LST' sd_names = ['LST_Day_1km'] filecount = SIZE(input_filenames, /DIMENSIONS) print,filecount FOR i = 0, filecount[0] - 1 DO BEGIN tmp_filename = input_filenames[i] ; fname = getfilename(tmp_filename) fname = file_basename(tmp_filename) ; out_rootname = STRMID(fname, 0, STRLEN(fname) - 22) ff=strsplit(fname,'.',/extract) ; print,ff[2] data=STRMID(ff[1], 1, 7) out_rootname =strtrim(string(data),2) CONVERT_MODIS_DATA, IN_FILE = tmp_filename, OUT_PATH = output_directory, $ OUT_ROOT = out_rootname, /HIGHER_PRODUCT, /GRID, $ Gd_NAME = grid_name, OUT_METHOD = out_method,out_ps_x=out_ps_x, out_ps_y=out_ps_y, $ INTERP_METHOD = interpolation_method, SD_NAMES = sd_names, $ NUM_X_PTS = 50, NUM_Y_PTS = 50, /NO_MSG, OUT_PROJ = output_projection, $ BACKGROUND = -999, FILL_REPLACE_VALUE = -999,r_fid_array=r_fid_array, r_fname_array=r_fname_array ENDFOR ENDIF PRINT, 'Job Done! END anyone can help me Please
    You are not authorized to post a reply.