[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[McIDAS #ZNQ-188340]: Displaying SST via L2 GOES-16 Products
- Subject: [McIDAS #ZNQ-188340]: Displaying SST via L2 GOES-16 Products
- Date: Tue, 13 Apr 2021 12:06:03 -0600
Hi Carol,
re:
> I'm trying to display the GOES-16 L2 SST product over the Carribean. See
> the attached image to witness my first go at this.
OK. I just created a display like yours.
re:
> First, I created the various L2 datasets via desserve using ABIN as the
> format. It seems like that was the correct format to use based on the User
> Guide
> <https://www.ssec.wisc.edu/mcidas/doc/users_guide/current/dsserve.html>
> section
> about the GOES-R series.
Depending on the version of McIDAS you are running, ABIN would be correct.
re:
> So I'm getting confused based on what to use for the EU and SU commands. I
> see that there is a L2-SST.ET and SST.ST which I figured I could use.
The enhancement L2-SST.ET is correct. The stretch table SST.ST might
be wrong since it was crated for the old SST product that was distributed
before the GOES-R series of satellites.
re:
> Unfortunately, it seems like something isn't getting translated properly in
> the units of the BAR. I don't think the SST values range from 3-255 K.......
>
> Here are the commands I'm running to get this plot. I couldn't get the SU
> to work....
>
> imgremap.k RTGOESR/L2FDSST MYDATA/AREAS.10 TIME=15:00:19 LATLON=20 75
> PRO=MERC RES=2.5 SIZE=1024 1600 SSIZE=ALL BAND=25
> imgdisp.k MYDATA/AREAS.10 MAG= EU=L2-SST.ET SU=X
> map.k X 1 FILE=OUTLUSAM WID=2
> map.k X 1 FILE=OUTLHPOL WID=2
> bar.k X COLOR=2 SU=X
>
> When running with the SU=SST, it can't find the correct units.
>
> imgdisp.k MYDATA/AREAS.10 MAG= EU=L2-SST.ET SU=SST
> imgdisp.k: Units requested are not available for this image
> imgdisp.k: done
This is undoubtedly being caused by using a stretch table that
was created for a previous incarnation of SST images where the
calibrated unit was degK.
You can see what the calibrated unit of the image you displayed
is by:
- IMGDISP the image
- move the cursor over an interesting feature
- run the D command (e.g., D<return> or CTRL-D)
I just did this and got the following output:
Image Name Day Nominal Time Scan Time Band
---------------- ------- ------------ --------- ----
NPGOESR/SSTFD.96 13 Apr 21103 16:00:19 16:03:24 25
File Nominal Image RAW SST BRIT
Lat/Lon Line/Element Line/Element K
16:26:18/ 74:08:29 1824/ 2757 7297/11029 47918 297.00
185
IMGPROBE: Done
So, the unit of the GOES-16 L2 SST product is 'K', not 'degK'.
re:
> I really haven't worked directly with stretch tables very much, so there
> could be something that I'm totally missing.
A new stretch table should be made and used with the GOES-16 L2 SST
images.
Here is what to try:
IMGLIST NPGOESR/SSTFD FORM=ALL
Image file directory listing for:NPGOESR/SSTFD
Pos Satellite/ Date Time Center Res (km) Image_Size
sensor Lat Lon Lat Lon
--- ------------- ------------ -------- ---- ---- ----- ----- ------------
96 GOES-16 L2 13 APR 21103 16:00:19 0 75
Band: 25 L2: SST - Sea Surface Temperature 2.02 2.00 5424 x 5424
proj: 0 created: 2021103 160019 memo: GOES Full Disk - Mode 6
type:ABIN cal type:RAW
offsets: data= 768 navigation= 256 calibration= 768 auxiliary= 768
doc length: 0 cal length: 0 lev length: 0 PREFIX= 0
valcod: 0 zcor: 0 avg-smp: N
lcor: 1 ecor: 1 bytes per pixel: 2 ss:187
Resolution Factors (base=1): Line= 4.0 Element= 4.0
IMGLIST: done
This shows that 'type' is ABIN.
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
unit | calibration units for the input table; valid units are RAW,
BRIT, RAD, TEMP, ALB or user-defined units; use command IMGPROBE
or D to list available units for an image (def=BRIT)
inlo | low value of the breakpoint for input values; specify in units
input with the 'unit' parameter
inhi | high value of the breakpoint for input values; specify in units
input with the 'unit' parameter
britlo | low values of the breakpoints for output brightness values;
the range is 0-255
brithi | high values of the breakpoints for output brightness values;
the range is 0-255
...
1) create a new, empty stretch table
SU INI L2-SST ABIN SST
2) define a mapping between calibrated unit and brightnesses
SU MAKE LW-SST 270 307 3 251
NB: you can create any mapping that you want. The most logical
thing (to me, at least) is to sample values in the full range of
brightnesses and use those values for the mapping. The values
I used here simply mimic the values in the existing SST.ST
stretch table.
3) list out the values in the stretch table just created:
SU TABLE L2-SST
BREAKPOINTS STORED IN TABLE : L2-SST.ST
INPUT OUTPUT
----- ------
270 3
307 251
CALIBRATION TYPE : ABIN
CALIBRATION UNITS : K
BAND NUMBER : -1
INTERPOLATION TYPE: LIN
SU: DONE
4) specify the use of the L2-SST.ST stretch in a BAR command and/or
in the IMGDISP you are using to display the image
Please let me know if/when you have questions.
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: ZNQ-188340
Department: Support McIDAS
Priority: Normal
Status: Closed
===================
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.