Skip to content. | Skip to navigation

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

example4.bat

;This is the batch file for simap.pro
;Intended for stitching of L1B 1km direct broadcast granules,
;or any other L1B 1km.
;Can be reused for all other resolutions or L2. Just give the right parameters
;names and channels.



.full_reset_session
.compile simap.pro

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


;====Input and Output path for source and output data==
inpath="/var/scratch/"
outpath="/var/scratch"
;======================================================



;===This is your file selection:======
hdfname="filelist"
;=====================================




;===This is your parameters and channel selections:===
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=0 ;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=0 ;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=0 ;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=15.
latmax=50.
lonmin=-95. ;westboundcoordinate
lonmax=-70. ;eastboundcoordinate

map_limit=[latmin,lonmin,latmax,lonmax]
center_lat=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 RGB planes:
intens=FLTARR(2,3)

intens[*,0]=[5,90] ;red
intens[*,1]=[15,90] ;green
intens[*,2]=[40,110] ;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=800 ;Request in pxl on the size of map
NYM=1000 ;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