Applies mathematical functions to image data.
IMGOPER sdataset1 . . sdatasetn ddataset [keywords] "comment
sdataset1 . . n |
source ADDE dataset names and positions; specify one of the following formats: |
|
|
|
group/descriptor.position alias.position |
|
a position greater than zero represents an absolute position in the dataset; a position less than or equal to zero represents a relative position in the dataset based on image time; for example, 0 is the most recent and -1 is the next most recent; to use default position, do not specify .position (no default for group/descriptor or alias; default=0 for position) |
|
ddataset |
destination ADDE dataset name and absolute position; specify one of the following formats: |
|
|
|
group/descriptor.position alias.position |
|
use only positive integers position |
|
"comment |
description placed in the image file directory's memo field; 32 characters maximum (default=no comment) |
|
LATlon= |
lat lon |
latitude and longitude of sdataset images to place at ddataset image location specified in the PLACE keyword |
||
LINele= |
line ele sys |
|||
|
|
line |
line number of the sdataset image to place at the ddataset image location specified in the PLACE keyword (default=0) |
|
|
|
ele |
element number of the sdataset image to place at the ddataset image location specified in the PLACE keyword (default=line) |
|
|
|
sys |
coordinate system of the line and ele values; specify either F for file (area) coordinates or I for image coordinates (default=F) |
|
PLAce= |
location in the ddataset image to place the point specified in the LATLON, STATION or LINELE keyword; use one of the following: |
|||
|
ULEFT |
upper-left corner of the image (default when the LINELE keyword is used) |
||
|
CENTER |
center of the image (default when the LATLON or STATION keyword is used) |
||
STAtion= |
station ID to place at the ddataset image location specified in the PLACE keyword; specify as the ID, e.g., KMSN or YSSY, optionally followed by the station type in brackets, e.g., ARX[N] or ARX[V]; see the TYPE keyword in the STNLIST command for valid station types |
|||
DAY= |
selects from sdataset images with the specified day |
|
ID= |
selects from sdataset images with the specified NEXRAD station ID, e.g., MKX |
|
RTIme= |
bmin emin |
selects from sdataset images within the specified range of minutes in the current hour; overrides the TIME keyword |
TIMe= |
btime etime |
selects from sdataset images within the specified time range |
BANd= |
b1 . . bn |
applies the function to the specified image bands (default=existing band for single-band images; no default for multiband images) |
LLMt= |
l1 . . ln |
lower limit thresholds for each input data value in the equation; if input value is less than the specified limit, the output value is zero |
MAG= |
lmag emag |
line and element magnification of the data; use a positive integer for blowup and a negative integer for blowdown (default=1 lmag) |
MISs= |
m1 . . mn |
missing data values for source datasets; if input value equals missing data value, the output value is zero (default=0) |
|
NONE |
all pixels from source datasets are included in the analysis |
SIZe= |
line ele |
applies the function to the specified number of sdataset lines and elements (default=480 640) |
|
ALL |
applies the function to the entire sdataset image; all image positioning keywords are invalid with SIZE=ALL |
TRACKing= |
YES |
tracks image if allowed by the server; tracking allows commands to receive data as it's being ingested by the server; commands complete when their data request has been fulfilled; tracking is only supported on some servers (e.g., Area, GVAR and MTSAT) and it can be disabled by the server administrator |
|
NO |
does not track image so server only sends data that is currently available (default=value of MCTRACK environment variable; if MCTRACK is not set, default=value set by server or server administrator) |
ULMt= |
u1 . . un |
upper limit thresholds for each input data value in the equation; if input value is greater than the specified limit, the output value is zero |
UNIt= |
calibration type for the data, for example, BRIT, RAD, TEMP (default=BRIT) |
|
ACOn= |
additive constant for the equation (default=0.0) |
||
COEf= |
c1 . . cn |
coefficients for each term in the equation (default=1.0) |
|
FORm= |
ADD |
addition; adds each term in the equation; to subtract terms, specify COEF=-1 or SIGN=M (default) |
|
|
CMAX max new |
conditional maximum; uses only the term in the equation with the highest value less than the specified maximum; if all terms are greater than max, the new value is used (default=255 255) |
|
|
CMIN min new |
conditional minimum; uses only the term in the equation with the lowest value greater than the specified minimum; if all terms are less than min, the new value is used (default=0 0) |
|
|
MAX |
maximum; uses only the term in the equation with the maximum value |
|
|
MIN |
minimum; uses only the term in the equation with the minimum value |
|
|
MULT |
multiplication; multiplies each term in the equation; to divide terms, specify POW=-1 |
|
FUNc= |
transcendental function used in the equation; use one of the following: |
||
|
ALOG |
antilogarithm |
|
|
EXP |
exponential function; the output pixel values are set to e raised to the power of the input pixel values (where e is the mathematical constant that is the base of the natural logarithm, approximately 2.718) |
|
|
LN |
natural log |
|
|
LOG |
logarithm |
|
|
SQRT |
square root |
|
MCOn= |
multiplicative constant for the equation (default=1.0) |
||
OFF= |
o1 . . on |
offsets for each term in the equation (default= 0.0) |
|
POWer= |
p1 . . pn |
powers for each term in the equation (default=1.0) |
|
SIGn= |
s1 . . sn |
multiplies the data values by 1 or -1 |
|
|
|
P |
plus; multiplies by 1 (default) |
|
|
M |
minus; multiplies by -1 |
CF= |
YES |
creates a CF-compliant netCDF file (default) |
|
NO |
creates a netCDF file in the previous (pre-2009.2) format, which is not CF-compliant; the CF=YES/NO keyword is valid only when writing to netCDF format destination datasets; it is ignored when writing to datasets of other formats |
DOC= |
YES |
copies the sdataset image line documentation to the ddataset image (default=NO) |
GAMma= |
gamma correction to apply to the output (default=1.0) |
|
PROd= |
user-defined product name for the ddataset image, up to four characters (default=PROD) |
|
SCAle= |
prodlo prodhi britlo brithi minimum and maximum product and brightness values for the calibration (default=0 255 0 255) |
|
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 determines how the terms are combined. For example, if you specify FORM=MULT, the terms are multiplied.
The Gamma function is then applied to the data using the following equation.
output data value =
Before the input data values are used in the equation, they are checked against the lower and upper limit thresholds set with the keywords LLMT and ULMT. If you do not specify limit thresholds, all data is used.
To use data from more than one band in a single multibanded image, you must specify that image with the sdataset parameter for each band used. For example, the following entry subtracts data values for band 10 from band 8, both in image A/A.1, and places the resulting image in A/A.5.
IMGOPER A/A.1 A/A.1 A/A.5 FORM=ADD BAND=8 10 COEF=1 -1
All source images must have the same projection, resolution and navigation. If they do not, use the command IMGCOPY or IMGREMAP before using IMGOPER. The source image can have unlimited number of lines and no more than 22,000 elements.
The destination image receives the directory information from the last source image. To change the image directory information, use the IMGCHA command.
To get better precision for operations that result in small values, use the keywords MCON and SCALE.
When using multiple source images, you must specify keyword values in the same order as the source images: value1 applies to the data in sdataset1. If a keyword allows for multiple values, the last entered value is the default for all subsequent values. The default listed for the keyword only applies if the keyword is not used. In the following command, for example, because an offset was specified only for images A/A.1 and A/A.2, image A/A.3 receives an offset of 10, the last entered value.
IMGOPER A/A.1 A/A.2 A/A.3 D/D.5 FORM=ADD OFF=1 10
Currently each data value in the output image, ddataset, can contain only 1-byte data. If the range of output data values is greater than 256, the output data loses accuracy.
You can use IMGOPER to create a netCDF image file. To do so, the destination dataset must be NCDF-format and IMAGE-type. See the Remarks section in the DSSERVE command for more information.
IMGOPER A/A.1 A/A.1 A/A.5 FORM=ADD BAND=8 10 SIGN=P M
This entry subtracts data values of band 10 from band 8, both in the image in position 1 in dataset A/A, and places the resulting image in position 5 in dataset A/A.
IMGOPER A/A.1 A/A.4 A/A.5 FORM=ADD BAND=2 3 COEF=1 -2
This entry multiplies the data values of band 3 of the image in position 4 in dataset A/A by 2 and subtracts them from band 2 of the image in position 1 in dataset A/A and places the resulting image in position 5 in dataset A/A.
IMGOPER OLD/DATA.1 OLD/DATA.1 OLD/PROD BAND=8 10 COEF=1 -1 SCALE=-10 10 PROD=DIFF UNIT=TEMP
This entry subtracts band 10 temperature data from band 8 temperature data of the image in position 1 in dataset OLD/DATA and places the resulting image in dataset OLD/PROD in the position following the most recent image. The output values are given the product name DIFF. The difference values between -10 and 10 are scaled from 0 to 255.
IMGOPER MSI.4 MSI.4 RATIO.1 FORM=MULT BAND=7 8 POW=1 -1 PROD=RATI SCALE=50 150 50 200 UNIT=BRIT MCON=100
This entry divides the brightness values from band 7 by the brightness values from band 8 of the image in position 4 in the dataset with the alias name MSI. The resulting image is placed in position 1 in the dataset with the alias name RATIO. The output values are given the product name RATI and multiplied by 100. The product values from 50 to 150 are scaled from brightness 50 to 200.
IMGOPER GOES/IR.1 GOESIR.2 PROD.4 BAND=8 8 COEF=-1 1 SCALE=-10000 10000 UNIT=RAD PROD=TEST MCON=1000
This entry subtracts band 8 radiance data in the image in position 1 in dataset GOES/IR from band 8 radiance data in the image in position 2 in dataset GOES/IR and places the resulting image in position 4 in dataset PROD. The output values are given the product name TEST and multiplied by 1000. The product values from -10000 to 10000 are scaled from brightness 0 to 255.
IMGOPER I/TEST.8 I/PROD.4 FUNC=LOG SCALE=0 20 0 255 UNIT=RAD MCON=1000
This entry multiplies the data values in the image in position 8 in dataset I/TEST by 1000 and then takes the logarithm. The results from 0 to 20 are scaled from brightness 0 to 255 and placed in position 4 in dataset I/PROD.
IMGOPER A/ALL.8344 A/ALL.4004 MISS=0 FUNC=LN PROD=NTLN SCALE=0 6 UNIT=BRIT
This entry takes the natural logarithm of the brightness image data in the image in position 8344 in dataset A/ALL and places the resulting image in position 4004 in the same dataset. The output values are given the product name NTLN. The product values from 0 to 6 are scaled from brightness 0 to 255.
IMGOPER G-7.1 MET3.1 A.1 SCALE=0 130050 0 255 PROD=PRMS POW=2 FORM=ADD "Intermediate
RMS
IMGOPER A.1 A.2 PROD=RMS
FUNC=SQRT MCON=0.5 UNIT=PRMS "Actual RMS
These entries compute the Root Mean Square (RMS) difference between the data in the images in position 1 in the datasets with the alias names G-7 and MET3. The first entry computes the sum of the squares of the data. The second entry computes the average square root of the results of the first command. Before running IMGOPER, the two source images must be remapped so that they have identical navigation.
IMGOPER GOES7/IR.30 METEOSAT/IR.3 COMPOSITE.1 FORM=CMIN SIZE=ALL
This entry takes the entire image in position 30 in dataset GOES7/IR and the entire image in position 3 in dataset METEOSAT/IR and combines the two images using the CMIN option. Before running IMGOPER, the two source images must be remapped so that they have identical navigation.
IMGOPER TEST/AREA.1 TEST/AREA.2 LLMT=108 ULMT=144 SCALE=108 144 2
38
IMGOPER TEST/AREA.1
TEST/AREA.3 LLMT=148 ULMT=196 SCALE=148 196 43 136
IMGOPER TEST/AREA.2
TEST/AREA.3 TEST/AREA.4 FORM=ADD MISS=NONE
These entries apply a series of scaling to the brightness values in the image in position 1 in dataset TEST/AREA. All brightness values between 0 and 107 are set to zero, all brightness values between 108 and 144 are scaled from 2 to 38, and all brightness values between 148 and 196 are scaled from 43 to 136. All values outside the specified LLMT and ULMT values are assigned a value of zero. The third command adds the output from the first two together so that TEST/AREA.4 has the desired results.
IMGCOPY NEXRAD/BREF1.1 DTX/RADAR.51 ID=DTX SIZE=ALL
IMGCOPY NEXRAD/BREF1.3 DTX/RADAR.52 ID=DTX SIZE=ALL
IMGCOPY NEXRAD/BREF1.5 DTX/RADAR.53 ID=DTX SIZE=ALL
IMGCOPY NEXRAD/BREF1.7 DTX/RADAR.54 ID=DTX SIZE=ALL
IMGOPER DTX/RADAR.51 DTX/RADAR.52 DTX/RADAR.53 DTX/RADAR.54 DTX/RADAR.55 FORM=ADD COEF=0.25 0.25 0.25 0.25
These four IMGCOPY entries copy the four images in positions 1, 3, 5 and 7 of dataset NEXRAD/BREF1 to positions 51-54 of dataset DTX/RADAR. The IMGOPER entry then computes the mean of the four images and writes the output to the image in position 55 of dataset DTX/RADAR.