Hi, Miae!
You need to add .he2 to the HDF file name as illustrated in [1] to
retrieve lat/lon. You can find more MOIDS examples from [2].
Regards,
[1]
http://hdfeos.org/zoo/LPDAAC/MYD09A1.A2007273.h03v07.005.2007285103507_sur_refl_b02.ncl
[2] http://hdfeos.org/zoo/index_openLPDAAC_Examples.php#MOD
-- HDF: Software that Powers Science On Tue, Dec 17, 2013 at 7:44 PM, 김미애 <toa0710@naver.com> wrote: > Hello, > > > > I am Miae Kim at Ulsan National Institute of Science and Technology > (UNIST) in South Korea. > > I am using MODIS data including MOD11A2, MOD13A2, MOD15A2, MOD17A2, > MOD16A2, MOD09A1 with NCAR Command Language (NCL) for my study. > > To do weighted average from 8-day MODIS data to monthly data, I did coding > for the average with function named "dim_avg_wgt_n_Wrap". > > Although I got all the averaged monthly data, they did not have any > projection information. > > Therefore, I could not define their projection in ArcGIS because there was > no spatial information in the calculated data. > > I think.. I cannot make coding for projection information because the > MODIS data (.hdf) do not have longitude and latitude information. > > Is there any way I can add projection information in the code? The > following is my code for weighted average. > > > > Thank you so much. > > > > ========================================================================== > > fpath = "/MODIS/2009/h28v05/" > > fpaths = systemfunc("ls "+fpath+"MOD11A2.A2009*.hdf") > > f1 = addfiles(fpaths,"r") > > variable= f1[:]->LST_Day_1km > > > > dims= dimsizes(variable) > > ny=dims(0) > > nx=dims(1) > > nfile = ny/nx > > > > input=new((/nfile,nx,nx/),float) > > do i=0,nfile-1 > > input(i,:,:) = variable(i*nx:i*nx+nx-1,:)*0.02 ;t*x*y > > end do > > input@_FillValue=0 > > > > wgt1=(/1, 1, 1, 0.875/) > > wgt2=(/0.125, 1, 1, 1, 0.5/) > > wgt3=(/0.5, 1, 1, 1, 0.375/) > > wgt4=(/0.625, 1, 1, 1, 0.125/) > > wgt5=(/0.875, 1, 1, 1/) > > wgt6=(/1, 1, 1, 0.75/) > > wgt7=(/0.25, 1, 1, 1, 0.625/) > > wgt8=(/0.375, 1, 1, 1, 0.5/) > > wgt9=(/0.5, 1, 1, 1, 0.25/) > > wgt10=(/0.75, 1, 1, 1, 0.125/) > > wgt11=(/0.875, 1, 1, 0.875/) > > wgt12=(/0.125, 1, 1, 1, 0.75/) > > > > avg1=dim_avg_wgt_n_Wrap(input(0:3,:,:), wgt1, 1, 0) > > avg2=dim_avg_wgt_n_Wrap(input(3:7,:,:), wgt2, 1, 0) > > avg3=dim_avg_wgt_n_Wrap(input(7:11,:,:), wgt3, 1, 0) > > avg4=dim_avg_wgt_n_Wrap(input(11:15,:,:), wgt4, 1, 0) > > avg5=dim_avg_wgt_n_Wrap(input(15:18,:,:), wgt5, 1, 0) > > avg6=dim_avg_wgt_n_Wrap(input(19:22,:,:), wgt6, 1, 0) > > avg7=dim_avg_wgt_n_Wrap(input(22:26,:,:), wgt7, 1, 0) > > avg8=dim_avg_wgt_n_Wrap(input(26:30,:,:), wgt8, 1, 0) > > avg9=dim_avg_wgt_n_Wrap(input(30:34,:,:), wgt9, 1, 0) > > avg10=dim_avg_wgt_n_Wrap(input(34:38,:,:), wgt10, 1, 0) > > avg11=dim_avg_wgt_n_Wrap(input(38:41,:,:), wgt11, 1, 0) > > avg12=dim_avg_wgt_n_Wrap(input(41:45,:,:), wgt12, 1, 0) > > > > printVarSummary(avg1) > > fout=addfile("lst_mon1.nc","c") > > fout->lst=avg1 > > fout=addfile("lst_mon2.nc","c") > > fout->lst=avg2 > > fout=addfile("lst_mon3.nc","c") > > fout->lst=avg3 > > fout=addfile("lst_mon4.nc","c") > > fout->lst=avg4 > > fout=addfile("lst_mon5.nc","c") > > fout->lst=avg5 > > fout=addfile("lst_mon6.nc","c") > fout->lst=avg6 > fout=addfile("lst_mon7.nc","c") > fout->lst=avg7 > fout=addfile("lst_mon8.nc","c") > fout->lst=avg8 > fout=addfile("lst_mon9.nc","c") > fout->lst=avg9 > fout=addfile("lst_mon10.nc","c") > fout->lst=avg10 > fout=addfile("lst_mon11.nc","c") > fout->lst=avg11 > fout=addfile("lst_mon12.nc","c") > fout->lst=avg12 > > end > ======================================================================= > > > Thank you so much again. > > Sincerely > > > > Miae > > ------------------------------------------------- > Miae Kim > School of Urban and Environment Engineering > Ulsan National Institute of Science and Technology (UNIST) > UNIST-gil 50, Ulsan 689-798, South Korea > Mobile: +82-10-9421-8287 > E-mail: toa0710@naver.com > > _______________________________________________ > ncl-talk mailing list > List instructions, subscriber options, unsubscribe: > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > >
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 17 20:20:08 2013
This archive was generated by hypermail 2.1.8 : Wed Dec 18 2013 - 14:29:51 MST