I still have problem when using MCTK to batchly process mod10A2 modis data,the pro file is following:
pro test_batch_modis_conversion_grid
compile_opt idl2
FilePath = 'f:\MOD10A2input\'
FileSearch = FILE_SEARCH(FilePath, COUNT = ProductCount, 'MOD10*.HDF', /TEST_READ, /FULLY_QUALIFY_PATH)
IF(ProductCount LE 0) THEN BEGIN
Print, 'There is no MODIS MOD10 to be processed.'
RETURN
ENDIF
FOR i=0, ProductCount-1 DO BEGIN
File = FileSearch[i]
Filename = FILE_BASENAME(File)
A = STRPOS(FILENAME,'.')
output_rootname = STRMID(FILENAME,A+1,15)
;The specified output location MUST end in the appropriate path
;separator for your OS
output_location = 'f:\modoutput\'
;output_rootname = Filename
grid_name = 'MOD_Grid_Snow_500m'
sd_names = ['Eight_Day_Snow_Cover']
;Output method schema is:
;0 = Standard, 1 = Reprojected, 2 = Standard and reprojected
out_method = 1
output_projection = envi_proj_create(/geographic)
;Specify the output X and Y pixel sizes in double precision floating
;point. Sizes must be in units relevant to selected output projection
;(degrees in this example).
out_ps_x = 0.004180d
out_ps_y = 0.004180d
;The interp_method keyword uses the same numbering scheme as
;ENVI_REGISTER_DOIT's METHOD keyword
;Choosing Triangulation with Nearest Neighbor.
interpolation_method = 6
;Set reprojection background and any native fill values to -999
;Setting num_x_pts to 50 and num_y_pts to 50 results in 2500 GCPs.
;That should be more than enough....
convert_modis_data, in_file=File, $
out_path=output_location, out_root=output_rootname, $
/higher_product, /grid, gd_name=grid_name, sd_names=sd_names, $
out_method=out_method, out_proj=output_projection, $
out_ps_x=out_ps_x, out_ps_y=out_ps_y, num_x_pts=50, $
num_y_pts=50, interp_method=interpolation_method, $
background=-999, fill_replace_value=-999, $
r_fid_array=r_fid_array, r_fname_array=r_fname_array
endfor
end
the error is intermittent showed as:
Any response is appreciated
|