McIDAS-X Learning Guide
Version 2020

[Search Manual] [Table of Contents] [Go to Previous] [Go to Next]

Satellite Imagery - Manipulating Images

1. Remap the second latest GOES-16 Visible image in the RTGOESR/FD dataset to a mercator projection centered on Washington DC, and place it in the fourth position in the TI dataset. Set the output resolution at 4 km. IMGREMAP calculates the size needed to make the default 480x640 image.

   Type:  IMGREMAP RTGOESR/FD.-1 TI.4 STATION=DCA BAND=2 PRO=MERC RES=4

2. List out the directory information for these images.

   Type:  IMGLIST RTGOESR/FD.-1 BAND=2 FORM=BAND;IMGLIST TI.4 FORM=BAND

   Image file directory listing for:RTGOESR/FD Pos Satellite/ Date Time Center Res (km) Image_Size sensor Lat Lon Lat Lon --- ------------- ------------ -------- ---- ---- ----- ----- ------------ 192 GOES-16 17 SEP 19260 10:00:18 0 75 Band: 2 0.64 um VIS clouds fog, insol, winds 0.50 0.50 21696 x21696 191 GOES-16 17 SEP 19260 09:50:18 0 75 Band: 2 0.64 um VIS clouds fog, insol, winds 0.50 0.50 21696 x21696 Image file directory listing for:TI Pos Satellite/ Date Time Center Res (km) Image_Size sensor Lat Lon Lat Lon --- ------------- ------------ -------- ---- ---- ----- ----- ------------ 4 GOES-16 17 SEP 19260 09:50:18 39 77 Band: 2 0.64 um VIS clouds fog, insol, winds 3.12 3.12 480 x 640

   Notice that the remapped image now has a different center lat/lon and is a different size than the original image.

3. Show frame 1. Display the images in frames 1 and 2 centered over Washington DC. Draw a map of dashed lat/lon lines in color 3 on each image. You will learn more about the MAP command in the Graphics and the Cursor lesson.

   Type:  SF 1

   Type:  IMGDISP RTGOESR/FD.-1 1 BAND=2 STATION=DCA REFRESH='MAP X X LALO -3 GRA=(GRA) IMA=(GRA)'

   Type:  IMGDISP TI.4 2 STATION=DCA REFRESH='MAP LALO -3 GRA=(GRA) IMA=(GRA)'

4. Set the loop sequence, and loop between the frames.

   Type:  LS 1-2

   Press:  Alt L

   Notice that the images have different projections.


5. Stop the loop.

   Press:  Alt L

6. Create a new image using two DORIAN visible GOES-16 images in the DORIAN dataset using the formula DORIAN/M1.1 - DORIAN/M1.10. This will show the shift in the clouds from one image to the next. Put the new image into position 5 of the TI dataset.

   Type:  IMGOPER DORIAN/M1.1 DORIAN/M1.10 TI.5 COEF=1 -1 FORM=ADD MAG=-2

   IMGOPER generates a new image by applying mathematical functions to data from one or more source images. The following equation computes the data value for each line/element pair in the destination image. This operation is performed repeatedly using data from source image line/element pairs as input values (input) until the entire image is completed.

   output data value=
   FUNC[ACON+(MCON*(FORM((COEF1*(OFF1+((SIGN1)*input1)**POW1))
   	(COEF2*(OFF2+((SIGN2)*input2)**POW2)) . .
   	(COEFn*(OFFn+((SIGNn)*inputn)**POWn)))))]

   
   

   The part of the equation operating on a single input data value, (COEFn*(OFFn+((SIGNn)*inputn)**POWn)) , is referred to as a term. n represents the number of sdataset images. It may not be larger than 100. input1 . . n represents the individual data values from each of the source images. The FORM keyword determines how the terms are combined. For example, if you specify FORM=MULT, the terms are multiplied.

7. Erase the maps in graphics frames 4 through 6, and display the images that you just used in frames 3-5.

   Type:  ERASE G 4 6

   Type:  IMGDISP DORIAN/M1.1 3 MAG=-2

   Type:  IMGDISP DORIAN/M1.10 4 MAG=-2

   Type:  IMGDISP TI.5 5

8. Loop through the images.

   Type:  LS 3-5

   Press:  Alt L

9. Stop the loop.

   Press:  Alt L

10. Add a label to the image frame.

    Type:  SF 5;FRMLABEL IMA=5 "CLOUD MOVEMENT FROM 18:00:25 TO 18:09:22

[Search Manual] [Table of Contents] [Go to Previous] [Go to Next]