NDVI mean 값을 산출하기 위해 아래와 같이 코드를 짰는데

WorkDir = ‘C:\all\’
WorkDirSat = ‘C:\all\’ ;Flux result
WorkDirOut = ‘C:\meann\’

syear = 2007
fyear = 2013

sday = 137
fday = 265

ncols = 4000
nrows = 4000

;psize = 0.01

data = fltarr(ncols,nrows,17)
mdata = fltarr(ncols,nrows)
stddata = fltarr(ncols,nrows)

data[*,*,*] = -9999
mdata[*,*] = -9999
stddata[*,*] = -9999

FOR year=syear, fyear DO BEGIN
FOR doy=sday, fday, 8 DO BEGIN

print, “Now processing Water Use Efficiency data from date: “, string(doy), ‘ yr ‘

Filename = strcompress(WorkDir+string(year,format='(i04)’)+’\’+string(year,format='(i04)’)+string(doy,format='(i03)’)+’\’+string(year,format='(i04)’)+string(doy,format='(i03)’)+’_NDVI_1km.img’,/remove_all)

ENVI_OPEN_FILE, Filename, r_fid=prcp_fid
ENVI_FILE_QUERY, prcp_fid, ns=prcp_ns, nl=prcp_nl, dims=prcp_dims
NDVI = ENVI_GET_DATA(dims=prcp_dims, fid=prcp_fid, pos=0) ;unit = mm/month

if(doy eq sday) then begin
proj = ENVI_GET_PROJECTION(FID=prcp_fid)
ENDIF

NDAY = DOY/8

data[*,*,NDAY] = NDVI

ENVI_FILE_MNG, ID=prcp_fid,/remove

ENDFOR
FOR y=0, nrows-1 DO BEGIN
FOR x=0, ncols-1 DO BEGIN
no = where(data[x,y,*] ne -9999, count1)
if (count1 gt 1) then begin
mdata[x,y] = mean(data[x,y,no])*0.0001
stddata[x,y] = stddev(data[x,y,no])*0.0001
endif
ENDFOR
ENDFOR
Outfile_mean = strcompress(WorkDirOut+string(year,format='(i04)’)+’_GROW_MEAN_1km.img’,/remove_all)
Outfile_std = strcompress(WorkDirOut+string(year,format='(i04)’)+’_GROW_STD_1km.img’,/remove_all)

ENVI_WRITE_ENVI_FILE, mdata, out_name=Outfile_mean, map_info=map_info, out_dt=4, r_fid=mean_map
ENVI_WRITE_ENVI_FILE, stddata, out_name=Outfile_std, map_info=map_info, out_dt=4, r_fid=std_map

ENVI_FILE_MNG, ID=mean_map,/remove
ENVI_FILE_MNG, ID=std_map,/remove
ENDFOR ;end of yearly loop

Datum = ‘WGS-84’
ps = [1000, 1000]
mc =[0.0, 0.0, -2751788.668d,394503.516d] ;clip image to array[140,85]
units =envi_translate_projection_units(‘Meters’)
Projection_Name = ‘Lambert Conformal Conic-NE_Asia’
Params = [6378137.0D, 6356752.3D, 38.000000D, 110.000000D, 0.0, 0.0, 30.0, 60.0]
proj = ENVI_PROJ_CREATE(type=4, name=Projection_name, datum=datum, params=params)
map_info = ENVI_MAP_INFO_CREATE(type=4, name=Projection_name, proj=proj, datum=datum, params=params, mc=mc, ps=ps)

;print, data[2600,2400,*]
;print, mdata[2600,2400]

PRINT, ‘The end!’
close,/all
END

이렇게 오류가 뜨는데 혹시 잘 아시는 분 도와주세요

NDVI = ENVI_GET_DATA(dims=prcp_dims, fid=prcp_fid, pos=0) ;unit = mm/month
^
% Syntax error.
At: H:\code\MEAN_ndvi_mod13.pro, Line 37

proj = ENVI_GET_PROJECTION(FID=prcp_fid)
^
% Syntax error.
At: H:\code\MEAN_ndvi_mod13.pro, Line 40

proj = ENVI_PROJ_CREATE(type=4, name=Projection_name, datum=datum, params=params)
^
% Syntax error.
At: H:\code\MEAN_ndvi_mod13.pro, Line 76

map_info = ENVI_MAP_INFO_CREATE(type=4, name=Projection_name, proj=proj, datum=datum, params=params, mc=mc, ps=ps)
^
% Syntax error.
At: H:\code\MEAN_ndvi_mod13.pro, Line 77
% 4 Compilation error(s) in module MEAN_NDVI_MOD13.
% Attempt to call undefined procedure/function: ‘MEAN_NDVI_MOD13’.
% Execution halted at: $MAIN$
IDL>