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

[McIDAS #AIT-655217]: Ticket ID: FVC-112860



Hi Paul,

As you guessed, I was at the AMS convention in New Orleans, and my ability
to reply to support email was limited.

re: previous explanation
> OK...I think that makes sense although I am still looking through a weird
> myriad of product types but am doing well so far.

There are lots of NEXRAD Level III products all right!

re:
> I do have a question concerning the enhancement tables and the stretch
> utilities for these products.

OK, ready.

re:
> I am working on base reflectivity for the digital. The type is what?

Type???  I am not sure I understand what you are asking.

re:
> I know it goes from 0-255 in brightness levels but how does it range in
> dBZ's?

The range of values in a product is not fixed.  Like many other Level III
products, the range in values can change from scan to scan.  The product
itself defines the lowest value and increments, and the ADDE server uses
these values to set the display values reflected by BAR.

re:
> Also, when I run the BAR command i am getting weird intervals.
> 
> see http://climate.cod.edu/data/nexrad/SJT/N0Q/SJT.N0Q.20120125.0230.gif

I just looked at your GIF, and I don't think that the intervals are weird
at all.  Again, the lowest value and increment can change from image to image.

re:
> I have tried the key word LIN=5 but McIDAS ignores that change.

Since the range is dynamic (see above), I pretty much don't allow the user
to change the mapping.

re:
> Here is a followup to my last email. I assume you are at the AMS
> conference.

Yup, I got back yesterday (Friday) night late).

re:
> I wrote this is a batch file but am not getting results i would expect. I
> want to map the dBZ levels to certain color levels...why is this not
> right?

> EU MAKE 0 2 0 0 0 0 0 0
> EU MAKE 3 25 85 200 85 200 85 200
> EU MAKE 26 54 128 141 128 4 128 14
> EU MAKE 55 69 80 216 80 215 80 0
> EU MAKE 70 84 0 0 255 200 0 0
> EU MAKE 85 99 0 0 200 144 0 0
> EU MAKE 100 114 0 0 144 255 0 255
> EU MAKE 115 129 0 0 255 192 255 231
> EU MAKE 130 144 0 0 192 144 231 255
> EU MAKE 145 159 0 0 144 60 255 255
> EU MAKE 160 174 0 0 0 0 255 154
> EU MAKE 175 189 0 255 0 0 214 255
> EU MAKE 190 203 255 201 0 85 255 153
> EU MAKE 204 218 201 255 85 255 153 255
> EU MAKE 219 233 255 255 255 255 255 0
> EU MAKE 234 248 255 50 255 0 0 0
> EU MAKE 249 254 0 50 255 0 255 255
> EU MAKE 255 255 255 255 255 255 255 255

I see nothing wrong with this set of commands aimed at creating a
new color enhancement.

re:
> SU INI PAULBREF RAW dBZ
> SU MAKE PAULBREF -30 -28 0 2
> SU MAKE PAULBREF -27.9 -.1 3 25
> SU MAKE PAULBREF 0 4.9 26 40
> SU MAKE PAULBREF 5 9.9 41 54
> SU MAKE PAULBREF 10 14.9 55 69
> SU MAKE PAULBREF 15 19.9 70 84
> SU MAKE PAULBREF 20 24.9 85 99
> SU MAKE PAULBREF 25 29.9 100 114
> SU MAKE PAULBREF 30 34.9 115 129
> SU MAKE PAULBREF 35 39.9 130 144
> SU MAKE PAULBREF 40 44.9 145 159
> SU MAKE PAULBREF 45 49.9 160 174
> SU MAKE PAULBREF 50 54.9 175 189
> SU MAKE PAULBREF 55 59.9 190 203
> SU MAKE PAULBREF 60 64.9 204 218
> SU MAKE PAULBREF 65 69.9 219 233
> SU MAKE PAULBREF 70 74.9 234 254
> SU MAKE PAULBREF 75 79.9 249 254
> SU MAKE PAULBREF 80 90 255 255

The problem I see is in the first line:

SU INI PAULBREF RAW dBZ

'RAW' is not the 'type' of the image.  Here is the error I see when
I run your SU commands above:

BATCH CODN0QSU.BAT
SU INI PAULBREF RAW DBZ
SU: Unsupported calibration type : RAW

The HELP for SU lets us know what to do:

HELP SU
SU -- Image data stretching utility
   SU INI name type unit <keywords>
   SU MAKE name inlo inhi britlo brithi
   SU LIST string
   SU TABLE name
Parameters:
   INI | initializes the stretch table defining the type of calibration,
         calibration units and initial breakpoint settings; must be done
         for making a new stretch table
   MAKE | enters breakpoints in the table using brightness levels
   LIST | lists the stretch tables saved on the workstation (def)
   TABLE | lists the breakpoints stored in a stretch table
   name | name of stretch table; 8 characters maximum
   string | list only stretch tables beginning with this character string
   type | calibration type for the input table (def=any calibration type);
          to find the calibration type of an image, run IMGLIST FORM=ALL
          and see the value for "type" in the output
 ...

An IMGLIST like the one called for is:

IMGLIST RTNEXRAD/N0Q ID=TYX FORM=ALL
Image file directory listing for:RTNEXRAD/N0Q
 Pos Satellite/         Date       Time      Center      Res (km)   Image_Size
     sensor                                 Lat  Lon    Lat   Lon
 --- -------------  ------------  --------  ---- ----  ----- ----- ------------
 139  RADAR         28 JAN 12028  21:35:39     TYX
   Band: 1    Base Reflectivity                         1.00  1.00   920 x  920
     proj: 1234 created: 2012028 213543  memo: BaseRef DR  0.50 DEG Precip Mode
     type:NEXR     cal type:RAW
     offsets:  data=    1964 navigation=  256 calibration=  768 auxiliary=  768
     doc length:   0   cal length:   0   lev length:   0 PREFIX=   0
     valcod:          0 zcor:  0 avg-smp: N
     start day: 2012028  start time: 213539  end day: 2012028  end time: 213539
     range:  460 KM  update interval:  5 MIN  tilt angle: 0.50 DEG
     maximum data value: 56.00 DBZ
     lcor:    1  ecor:     1  bytes per pixel: 1  ss:  7
     Resolution Factors (base=1):   Line=    1.0   Element=    1.0
IMGLIST: done

From this listing, we can see that the 'type' is NEXR (line #8 of the
listing).

Changing your SU INI invocation from:

SU INI PAULBREF RAW dBZ

to:

SU INI PAULBREF NEXR DBZ

(I also make DBZ all caps)

Allows the set of SU commands to run:

BATCH CODN0QSU.BAT
SU INI PAULBREF NEXR DBZ
SU: DONE
SU MAKE PAULBREF -30 -28 0 2
SU: DONE
SU MAKE PAULBREF -27.9 -.1 3 25
SU: DONE
SU MAKE PAULBREF 0 4.9 26 40
SU: DONE
SU MAKE PAULBREF 5 9.9 41 54
SU: DONE
SU MAKE PAULBREF 10 14.9 55 69
SU: DONE
SU MAKE PAULBREF 15 19.9 70 84
SU: DONE
SU MAKE PAULBREF 20 24.9 85 99
SU: DONE
SU MAKE PAULBREF 25 29.9 100 114
SU: DONE
SU MAKE PAULBREF 30 34.9 115 129
SU: DONE
SU MAKE PAULBREF 35 39.9 130 144
SU: DONE
SU MAKE PAULBREF 40 44.9 145 159
SU: DONE
SU MAKE PAULBREF 45 49.9 160 174
SU: DONE
SU MAKE PAULBREF 50 54.9 175 189
SU: DONE
SU MAKE PAULBREF 55 59.9 190 203
SU: DONE
SU MAKE PAULBREF 60 64.9 204 218
SU: DONE
SU MAKE PAULBREF 65 69.9 219 233
SU: DONE
SU MAKE PAULBREF 70 74.9 234 254
SU: DONE
SU MAKE PAULBREF 75 79.9 249 254
SU: DONE
SU MAKE PAULBREF 80 90 255 255
SU: DONE
BATCH: BATCH done /home/mcidas/workdata/CODN0QSU.BAT

re:
> Also, would i run su with the IMGDISP command or just the BAR command?

I would include it in the BAR invocation:

BAR SU=PAULBREF

PLEASE NOTE CAREFULLY:

- since the mapping from original value in the image to "calibrated unit"
  (dBZ in this case) is not fixed, your results should vary (i.e., not
  be consistent)

Cheers,

Tom
--
****************************************************************************
Unidata User Support                                    UCAR Unidata Program
(303) 497-8642                                                 P.O. Box 3000
address@hidden                                   Boulder, CO 80307
----------------------------------------------------------------------------
Unidata HomePage                       http://www.unidata.ucar.edu
****************************************************************************


Ticket Details
===================
Ticket ID: AIT-655217
Department: Support McIDAS
Priority: Normal
Status: Closed