[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20011024: conversion in AREA files

>From: "Jason J. Levit" <address@hidden>
>Organization: Center for Analysis and Prediction of Storms, University of Okla
>Keywords: 200110241902.f9OJ2n113278 McIDAS IMGCOPY


>I'm attempting to process McIDAS AREA satellite imagery files that
>contain data in the "RAW GVAR" format, and I'm trying to convert the
>data to "VISR BRIT" AREA format.  Is there an easy way to do this in

Yes.  The command you should be looking at is IMGCOPY.  You basically
"copy" the "RAW GVAR" data (which should be 2-byte) into a new dataset
specifying STYPE=VISR.  Here are some snippits from the online help

IMGCOPY -- Copies image data from a server
   IMGCOPY sdataset ddataset <keywords>
   sdataset | source ADDE dataset name and position; specify as
              alias.position or group/descriptor.position; to use default
              position, either enter "0" or omit .position portion
              (no def for alias or group/descriptor, def=0 for position)
   ddataset | destination ADDE dataset name and position; specify as
              alias.position or group/descriptor.position; only positive
              integers (absolute positions) or the default are valid for
              the position number; to use the default position, omit the
              .position portion (no def for alias or group/descriptor, for
              position def=the position following the most recent image)

*********************  Output Keywords  **********************************
   BANd=    | image band number to copy (def=existing band if image is
              single-band; no def if image is multi-band)
       =ALL | copies all bands; not valid with STYPE=VISR
   DOC=YES  | include the image documentation when copying the
              sdataset image data (def=NO)
   GAUge= | outputs text messages that indicate the percentage of the
            copy that has been completed; specify as a percentage between
            0 and 50; for example, if GAUGE=20, text messages are output
            when the copy is 20%, 40%, 60% and 80% complete
   MAG=lmag emag | line and element magnification of the data
                   (lmag def=1, emag def=lmag)
   OVERride=NO  | does not copy image if most recent image in ddataset
                  contains same date, time and band (def)
           =YES | copies image to ddataset position following the most
                  recent image even if the date, time and band are same;
                  this keyword is ignored if a position number is
                  specified in ddataset
   SIZe=line ele | number of lines and elements of the sdataset image to
                   copy (def=480 640 unless using POES data and specifying
                   C option in TIME keyword; see TIME keyword)
       =ALL | copies entire sdataset image including all bands and
              documentation; if sdataset image has multiple bands with
              different sizes and resolutions (e.g., Stretched Data Format
              files on SDI servers), all bands are copied at the lowest
              resolution and smallest size; the only other keywords valid
              with SIZE=ALL are DAY, RTIME, TIME, ALL, REPEAT and CYCLE
       =SAME  | retains the original sdataset image size; can be used to
                copy a single band from a multi-band image
   STYpe=VISR | reduces 2-byte data to 1-byte brightness data; not valid
                with BAND=ALL or SIZE=ALL
   UNIT= | calibration unit for output image file; use this keyword with
           NCDF-format destination datasets; to determine available units,
           display the source image and run the D command (def=primary
           calibration unit in source image)


>I checked out the "SU" command, but I wasn't sure if it would
>get me there, and thought I might be missing something easy.

SU is not the command you are looking for.

>Basically, the code we have to read the AREA format data and convert
>it for use in our NWP system expects the data to be in VISR BRIT
>format.  So, I need an AREA file which conforms to that file structure,
>if it's possible.

IMGCOPY is the routine for you!

>Thanks for any assistance!

You are welcome.

>ps: this is a separate issue from converting GINI data to McIDAS AREA,
>though I'm working closely with Kevin Thomas (CAPS) on that issue as
>well.  Thanks for all of your help with this!

I understand.  I sent Kevin a reply to the GINI issue a few minutes ago
(but I see that he has already responded to that one :-)


NOTE: All email exchanges with Unidata User Support are recorded in the Unidata inquiry tracking system and then made publicly available through the web. If you do not want to have your interactions made available in this way, you must let us know in each email you send to us.