IMGREMAP

Remaps images into different projections.

Format

IMGREMAP sdataset ddataset [keywords]

Parameters

sdataset	source ADDE dataset name and position; 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 for position

Image Selection Keywords

DAY=	selects from sdataset images with the specified day
DISt=	maximum distance, in kilometers, from the satellite subpoint to include in the remap; use SMOOTH=YES and/or SPLINE= with a smaller value to round the edges of the output image
ID=	selects from sdataset images with the specified NEXRAD station ID, e.g., MKX
MAG=	lmag emag	line and element reduction factor to apply to the source data before it's remapped; valid range is 1 through -999 where 1, 0, and -1 mean no reduction, -2 means a 2x reduction, etc; see the Remarks (lmag default=1, emag default=lmag)
SSIze=	line element	size of the source image region to remap; the region consists of the specified number of lines and elements centered at the location specified in the LATLON or STATION keyword (default=computed; see the Remarks)
	ALL	remaps the entire source image, even if only a small subset of the data is actually saved in the destination image with the SIZE keyword
TIMe=	btime etime	selects from sdataset images within the specified time range

Output Keywords

BANd=	image band number to remap (default=existing band for single-band images; no default for multiband images)
LATlon=	clat clon	center latitude and longitude for ddataset (default=center latitude and longitude of sdataset)
MERge=	NO	sdataset data overwrites the entire ddataset (default)
	YES	sdataset data overwrites ddataset only where images overlap; valid only with PRO=DEST
NAVtype=	LALO linc einc	use lat/lon navigation instead of the regular (native) navigation; this keyword is valid only if the source image is area format; linc and einc specify the image line and element increments, e.g., NAVTYPE=LALO 4 3 requests lat/lon nav using the data from every fourth line and every third element (default=use native nav; if you specify LALO, lat/lon nav is used and the linc and einc defaults are the computed values that result in using the largest number of data points allowable without exceeding the max size of the AUX block)
PRO=	DEST	use the navigation of ddataset (default)
	LAMB	slat1 slat2 slon Lambert conformal projection, where:
		slat1	standard latitude (default=30)
		slat2	standard latitude (default=50)
		slon	standard longitude (default=center longitude)
	MERC	slat	Mercator projection; specify standard latitude (default=0)
	MOLL	slon	Mollweide projection; specify standard longitude (default=0)
	PS	slat slon polar stereographic projection, where:
		slat	standard latitude (default=60)
		slon	standard longitude (default=center longitude)
	RADAR	rot	radar projection; specify rotation angle (default=0)
	RECT	rectilinear projection
	SIN	sinusoidal equal area projection
	TANC	slat slon tangent cone projection, where:
		slat	standard latitude (default=30 if the center latitude is in the northern hemisphere; -30 if the center latitude is in the southern hemisphere)
		slon	standard longitude (default=center longitude)
RES=	resolution of ddataset in km; not valid with PRO=DEST (default=1)
SAVesrc=	YES	writes sdataset's directory information (including SS, DATE, TIME and BAND) to ddataset; valid only with PRO=DEST (default=NO)
SIZe=	line ele	number of lines and elements in ddataset; the ele value must be a multiple of 4; SIZE is not valid with PRO=DEST (default=480 640)
	ALL	creates an image large enough to include the entire geographical coverage of sdataset
SMOoth=	limb data provides limb or data smoothing
		limb	NO	no limb smoothing (default)
			YES	provides limb smoothing
		data	NO	no data smoothing (default)
			YES	provides data smoothing using a 2/3:1/3 weighting of the two closest pixels
SPLine=	spline size in image pixels; the range is from at least 1 to 100; the range minimum may be computed based on source data (default=20)
STAtion=	station ID to place at the center of the remapped image; 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
WL=	remaps band with image wavelength closest to the specified value; valid only with high spectral resolution data, e.g., Aqua AIRS
WN=	remaps band with image wave number closest to the specified value; valid only with high spectral resolution data, e.g., Aqua AIRS

Remarks

IMGREMAP subdivides the destination image into boxes. The height and width of each of the boxes is determined with the SPLINE keyword, which defaults to 20 pixels. IMGREMAP performs a navigation transformation from line/element to latitutde/longitude at each of the spline box corners. It then determines the line/element coordinates in the source image that correspond to these same latitude/longitude values. Those corner values are used in a bi-linear interpolation scheme (which is much faster than a navigation transformation) to compute the line/element coordinate of each source image pixel it will use to fill the intermediate points in the spline box. This method ensures that each pixel in the destination image has an associated pixel from the source image.

If the destination image is lower resolution than the source image, the data are sampled. If the destination image is higher resolution, the pixel values are replicated. Only the 8-bit grayscale values are stored in the destination image. By default, there is no interpolation of the 8-bit values; however, the keyword SMOOTH=X YES performs an interpolation using a 2/3:1/3 weighting of the two closest pixels.

IMGREMAP uses the sdataset parameter and the DAY, TIME and ID keywords to select the source image to remap. It uses the Output Keywords to calculate the geographical region to extract from the source image for the remap. The calculation, however, sometimes underestimates or overestimates the size of the source image region needed for the remap. If that happens, use the SSIZE keyword to override the default calculation. For example, if the default calculation doesn't include the poles in the remap of a global image, specify SSIZE=ALL to ensure that they're included. The default calculation often overestimates the data region needed for a remap of polar-orbiting satellite data. To avoid this, you can use SSIZE with line and element values to choose a smaller source image region. Selecting more source image data than necessary does not affect the output image, but takes it longer to complete and adds to the cost if you pay for the data.

The MAG keyword can be used to further reduce the amount of source image data that is extracted and remapped. Use the MAG keyword only when specifying a destination image with less detail (lower resolution) than the source image. For example, if the source image is 1km resolution and you specify RES=30 for a 30km resolution destination image, you can also specify MAG=-16 to reduce the amount of data extracted and remapped by a factor of 16. It's best to pick a MAG value somewhat less than the destination image resolution because the quality of the destination image decreases as the MAG value approaches the resolution value.

When using MERGE=YES, all data in the source image (even pixels with zero data values) overwrite the collocated data in the destination image. Thus, all geographic locations in the destination image that are also in the source image will be replaced by the source image's values.

The radar projection is identical to the sinusoidal equal area projection, except that it allows you to rotate the destination image.

Display destination images with the IMGDISP command.

To create a destination dataset with a line prefix length and validity code of zero, use the IMGCOPY command with the keyword STYPE=VISR and SIZE=line ele; do not use SIZE=SAME.

For IMGREMAP to be successful using the PRO=DEST option, the following requirements must be met:

both the sdataset and ddataset must contain navigation
the ddataset must be a single band of 1-byte data (ZRES=1 and ZSIZ=1)
the sdataset may be 2- or 1-byte data
the ddataset's calibration type must be VISR
the ddataset's source type must be BRIT
the ddataset must have a line prefix length of zero
the ddataset must have a validity code of zero

You can use IMGREMAP 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.

Examples

IMGREMAP LOCAL/IMAGE.16 LOCAL/REMAP.4000

This entry remaps the image data in position 16 of dataset LOCAL/IMAGE into position 4000 of dataset LOCAL/REMAP. The source image replaces the destination image but uses its projection information.

IMGREMAP NOAA-12/HRPT NO-GOES/DATA.10 MERGE=YES PRO=DEST

This entry remaps the most recent image in dataset NOAA-12/HRPT into position 10 of dataset NO-GOES/DATA. The source image uses projection information from the destination image and replaces data only where the images overlap.

IMGREMAP G8/FDV G8/REMAP.1 RES=24 PRO=SIN

This entry remaps the most recent image in the G8/FDV dataset into the position 1 of dataset G8/REMAP. The destination image is remapped with a resolution of 24 kilometers in a sinusoidal equal area projection.

IMGREMAP G8/CONUSV G8/RADAR.2 RES=4 STATION=BOS PRO=RADAR 45

This entry remaps the most recent image in the G8/CONUSV dataset into position 2 of dataset G8/RADAR. The destination image is remapped in a radar projection, rotated 45 degrees. The image is centered at Boston, MA and has a resolution of four kilometers.

IMGREMAP WV.4 REMAP.4 PRO=LAMB 30 60 70 LATLON=34 78 RES=4

This entry remaps the image in position 4 of the dataset with the alias WV into position 4 of the dataset with the alias REMAP. The destination image is centered at 34° N and 78° W and has a resolution of 4 kilometers. The image is remapped in a Lambert conformal projection using standard latitudes of 30° and 60° N and a standard longitude of 70° W.

IMGREMAP GOES.1 A.2 PRO=PS 45 71 RES=4 SPLINE=40

This entry remaps the image in position 1 of the dataset with the alias GOES into position 2 of the dataset with the alias A. The destination image is remapped in a polar stereographic projection, with a standard latitude of 45° N and a standard longitude of 71° W. The resolution of the remapped data is four kilometers, and the spline value is 40.

IMGREMAP A/A.8016 A/A.4000 PRO=MERC 50 RES=2 SIZE=ALL

This entry remaps all of the image data in position 8016 of dataset A/A into position 4000 of dataset A/A. The destination image is as large as necessary to contain the entire remapped image: a Mercator projection, with a standard latitude of 50°N and a resolution of two kilometers.

IMGREMAP GOES/VIS LOCAL/AREA.23 BAND=1 DAY=#Y TIME=22:55 23:05 MAG=-8 RES=16 PRO=MERC LATLON=0 90

This entry remaps the image in dataset GOES/VIS with today's date and an image time between 22:55 and 23:05 into position 23 of dataset LOCAL/AREA. The destination image is in a Mercator projection centered at 0°N and 90°W with a resolution of 16 kilometers. The MAG=-8 keyword reduces the resolution of the data retrieved from the source image from 1km to 8km before remapping it to 16km resolution. This 8x reduction speeds up completion of the command with little effect on the quality of the destination image.