Skip to content. | Skip to navigation

Personal tools
You are here: GES DISC Home Services Data Tools simap example3.bat

example3.bat

;This is the batch file for simap.pro
;This batch sample will generate a single map that
;is an RGB image of bands 1,4,3.
;Substitute the paths and file names to your valids.



.full_reset_session
.compile simap

sel={name:"",chans:STRARR(38),numch:0}
sel=REPLICATE(sel,20)





inpath='/ftp/data/modis/samples/radiometric_geolocation/L1B/MOD021KM'
outpath='/var/scratch'


;===THis is your file selection============
hdfname='MOD021KM.A2000338.1810.002.2000359191516.hdf'
;==========================================


;===This is your parameter selection======
sel[0].name=['EV_250_Aggr1km_RefSB']
sel[0].chans=['1']

sel[1].name=['EV_500_Aggr1km_RefSB']
sel[1].chans=['4','3']
;==========================================




;====================FLAGS DEFINITIONS============================

intrpl=1 ;Set to 1 activate filling
;up of holes in the images. This may slow
;the process, but will clean the images.
;Set to 0 to deactivate.


extgeo=0 ;set to 1 to use EXTERNAL geolocations
;set to 0 to use INTERNAL geolocations

igeo=1 ;set to 1 to interpolate geolocations
;It must be 1 when using internal geolocations
;and L1B MODIS files

sdsprint=0 ;Flag to list SDS
;0=no list
;1=short list, SDS names only, no attributes
;2=long list, SDS names, attributes and their
;values


True_col=3 ;Flag to generate image after all procedures
;are completed.
;Set to 3 to save RGB image composite from the
;first 3 channels selected.
;Set to 1 to save grey/color scale separate
;images of all selected channels.


flcnt=0 ;set to 1 to fill continents
;set to 0 not to fill continents

automap=1 ;Set to 1 to automaticaly
;set map limits.
; Set to 0 to use the default
; map limits.


autointen=0 ;Set to 1 to automatically set image intensity.
;Set to 0 to use default image intensity
;ranges "intens"


zlog=0 ;Set to different than 0 to show image
;in log scale

gallery=1 ;Flag to generate galleries of image files.
;0=no gallery, will do global map
;1=gif or jpeg gallery of gray/color scale,
;of true color, RGB, image files. If
;True_col=1 then grey/color scale.
;If True_col=3 then RGB


debug=-1 ;debug mode: -1 will print to monitor
; 4 will print to log file
;============================================================




;============PARAMETERS DEFINITIONS===========================


;---------SUBSET-------
xstart="" ;leave as NULL string for full size
xend="" ;or give numerical value

ystart=""
yend=""

;---------Default map settings-----------------
latmin=21.
latmax=31.
lonmin=-118. ;westboundcoordinate
lonmax=-106. ;eastboundcoordinate

map_limit=[latmin,lonmin,latmax,lonmax]
center_lat=0;latmin+(latmax-latmin)/2.
center_lon=atan((sin(lonmin*!dtor)+sin(lonmax*!dtor)),$
(cos(lonmin*!dtor)+cos(lonmax*!dtor)))*!radeg


mappos=[.05,.08,.95,.98] ;map position and size,
; relative to the graphic window

proj=1 ;1 - cylindrical
;2 - stereo
;--------------------------------------------

;IMAGE INTENSITY
;default ranges for data planes:
intens=FLTARR(2,3)
intens[*,0]=[5.,80] ;red
intens[*,1]=[15.,80] ;green
intens[*,2]=[40.,100] ;blue

slice=3 ;For 3D arrays, this is the dimension
;to fix. In IDL, dimensions are
;e.g. array[x,y,z]. Set slice=1 or 2 or 3 to fix
; x or y or z dimension, correspondingly.
;For most SDS's slice=3, i.e. they are stacked
;over the z-dimension. However, for example MODIS
;Atmosphere L2 Quality_Assurance needs slice=1 !!!
;I.e. it is stacked over the x-dimension


NXM=2000 ;Request in pxl on the size of map
NYM=2000 ;to be produced.


ctab=0 ;IDL color table

backgnd=200 ;background inside map
;and ouside the granule(s)

;this will be the size of the graphic window to
;hold the map:
NX=NXM/(mappos[2]-mappos[0])
NY=NYM/(mappos[3]-mappos[1])
;================================================================





simap,hdfname

Document Actions
NASA Logo - nasa.gov
NASA Privacy Policy and Important Notices
Last updated: Sep 10, 2010 02:20 PM ET
Top