20051123: MAKEAREA (makearea.pgm) mods (cont.)

This is a MIME encoded message.  Decode it with "munpack"
or any other MIME reading software.  Mpack/munpack is available
via anonymous FTP in ftp.andrew.cmu.edu:pub/mpack/

>From: Unidata Support <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200511171722.jAHHMk7s019022 McIDAS MAKNAV

Hi Valentijn,

re: My inclination is to add support to MAKNAV to extract the navigation
block out an existing image and write it into the target.

>Do you have this adjusted MAKNAV program which can do this?

I modified ~mcidas/mcidas2005/src/maknav.pgm (v2005 source was the same
as v2004, so no worries if you are running v2004) and added the ability
to extract the navigation block out of an existing AREA file and put it
into an AREA created by MAKEAREA.  I also modified makearea.pgm to
increase the size of the navigation block section to the maximum
allowed by McIDAS (this made the modification to maknav.pgm easier and
is completely general).

I have attached the source code for MAKNAV (maknav.pgm) and MAKEAREA
(makearea.pgm) below.

To check it out, you need to do the following:

<as 'mcidas'>
cd ~mcidas/mcidas2005/src
mv maknav.pgm maknav.pgm.2005
mv makearea.pgm makearea.pgm.2005

- add attached maknav.pgm source 
make maknav.k maknav.hlp
rm ~/bin/maknav.k
rm -f ~/help/maknav.hlp
ln maknav.k ~/bin
ln maknav.hlp ~/help

make makearea.k makearea.hlp
rm ~/bin/makearea.k
ln makearea.k ~/bin

The invocation syntax for MAKEAREA remains the same

The new invocation syntax of MAKNAV is:

MAKNAV -- Create navigation entry for an area
     MAKNAV area <keywords>
Parameters:
      area - area to attach navigation (def=current frame)
Keywords:
      LAMB= npl npe slat1 slat2 qlon size pole | Lambert Conformal secant
      MERC= eql eqe slat qlon size             | Mercator
      PS=   npl npe slat qlon size pole        | Polar Stereographic
      RADR= cenl cene clat clon size rot       | Radar
      RECT= lin lat ele lon deglin degele      | Rectilinear
      SIN=  cenl cene clat clon size           | Sinusoidal Equal Area
      TANC= npl npe slat qlon size             | Tangent Cone
      MOLL= eql eqe qlon size                  | Mollweide
            cenl,cene - line (cenl) and element (cenl) of radar site
            clat,clon - lat (clat) and lon (clon) of radar site
            deglin,degele - degrees/line (deglin) and element (degele)
            ele,lon - image element (ele) of longitude (lon)
            eql,eqe - image line (eql) and element (eqe) of equator
            lin,lat - image line (lin) of latitude (lat)
            npl,npe - image line (npl) and element (npe) of pole
            pole - N(orth) or S(outh) pole (def=N)
            qlon - normal longitude (center)
            rot - rotation of radar picture (def=0)
            size - size of one pixel (km) at standard latitude
                 - or map scale factor (ratio)
            slat - standard latitude (PS def=60, MERC def=0)
            slat1- standard latitude (LAMB def=30)
            slat2- standard latitude (LAMB def=60)
      NOAA= kepler day time tmsec orient flip  - NOAA polar orbiter
            kepler - name of 3-line Keplerian element file (form in Notes)
            day    - nominal DAY for image
            time   - nominal TIME for image
            tmsec  - time of first line to navigate [msec from DAY]
            orient - ascending(+1)/descending(-1) pass (def=+1)
            flip   - inverted image flag (0-normal/1-flipped) (def=0)
      AREA= iarea ullin ulele lres eres
            iarea  - AREA number containing desired navigation (no default)
            ullin  - Image line coordinate of area line 0, elem 0
                       (def= ullin from iarea)
            ulele  - Image element coordinate of area line 0, elem 0
                       (def= ulele from iarea)
            lres   - Line Resolution; spacing in image-lines between
                       consecutive area lines (def= nlin from iarea)
            eres   - Element resolution; spacing in image-elements
                       between consecutive area elements (def= nele from
                       iarea)

 ...

I tested the new versions of MAKEAREA and MAKNAV with images containing
GOES and METEOSAT-7 navigation and both worked correctly.

The sequence taking uncompressed raster image data and adding it to 
an AREA file and then navigating it is as follows:

- get a raster image from somewhere
- determine the following for the raster image:

  - number of lines in the raster
  - number of elements in the raster
  - size of the header of the raster in bytes (number of bytes at the
      beginning of the raster before the data values)
  - how many bytes in the raster represent a data value; possible values
      are 1, 2, or 4
  - endiness of the raster (big endian or little endian)
  - minimum data value represented by the raster image values
  - maximum data value represented by the raster image values

- convert the raster to a McIDAS AREA file:

  MAKEAREA - Make a McIDAS area from a flat image file of brightnesses

  MAKEAREA file area nlins neles nhead nbyte min max

  Parameters:
    file   - name of flat image file (no default)
    area   - McIDAS area number (no default)
    nlins  - number of lines of data in flat file (no default)
    neles  - number of elements per line in flat file (no default)
    nhead  - size of flat image file header in bytes (default=0)
    nbyte  - number of bytes per pixel (default=1)
    min    - minimum value of data (default=0)
    max    - maximum value of data (default=depends on nbyte)

  Keywords:
    FLIP=YES - flip bytes in 2,4 byte images (default=NO)

  Remarks:
    MAKEAREA converts a flat file of brightnesses (1, 2, or 4-byte data)
    into a 1-byte AREA scaling the data appropriately

- identify an image in an AREA file that has the navigation you want
  to add to the newly created AREA from the raster

- add the navigation using MAKNAV:

MAKNAV area AREA=iarea
        ^          ^___ AREA file number of image containing desired navigation
        |______________ AREA file created by MAKEAREA

- you may need to play around with the 'ullin' 'ulele' 'lres' and 'eres'
    values to get the navigation registration correct

Please let me know when you have questions.

Cheers,

Tom

>From address@hidden  Thu Nov 24 06:49:27 2005

Tom, thanks so much!! I'll have a look at it as soon as the current
course is finished, and my life becomes less busy. I'll probably get
back to you at that time with more questions...

Cheers, Valentijn