edu.wisc.ssec.mcidasv.util

Class GEOSgridUtil

java.lang.Object

edu.wisc.ssec.mcidasv.util.GEOSgridUtil

public class GEOSgridUtil
extends java.lang.Object

Utilities for efficiently upscaling and downscaling data referenced on the GOES-16 2km, 1km and hkm Fixed Grid Frame (FGF).

Adapted from VisAD's resample, but realizes efficiencies when target and source domain sets are both rectilinear, especially avoiding the expansion of the domain samples (getSamples) for this special case.

Follows resampling rules outlined in the GOES-16 ABI ATBD for the Cloud and Moisture Imagery Product when incoming grids are at the their full resolution. Should be thought of as a grid transfer operation, to be used in conjunction with, not necessarily a replacement for, FlatField.resample.

Field Summary

Fields

Modifier and Type	Field and Description
private static float	accumTolerance
static RealTupleType	LambdaTheta
private static double	offsetTolerance
private static int	uniqueID

Constructor Summary

Constructors

Constructor and Description
GEOSgridUtil()

Method Summary

Modifier and Type	Method and Description
static FlatField	geosResample(FlatField fld, Linear2DSet targetSet) Efficiently upscales or downscales between ABI fields with domainSets on the Fixed Grid Frame.
static FlatField	geosResample(FlatField fld, Linear2DSet targetSet, int mode) Efficiently upscales or downscales between ABI fields with domainSets on the Fixed Grid Frame.
static FlatField	geosRGB(FlatField red, FlatField green, FlatField blue) DomainSets of red, green, blue must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame.
static FlatField	geosRGB(FlatField red, FlatField green, FlatField blue, Linear2DSet targetSet)
static FlatField	geosRGB(FlatField red, FlatField green, FlatField blue, Linear2DSet targetSet, boolean copy) DomainSets of red, green, blue and targetSet must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame.
(package private) static void	geosUpsample(int[] yidxs, int[] xidxs, int lenY, int lenX, float[][] values, int targetLenY, int targetLenX, float[][] targetValues, int mode)
(package private) static void	interp(float[] yidxs, float[] xidxs, int lenY, int lenX, float[][] values, int targetLenY, int targetLenX, float[][] targetValues, int mode)
static FlatField	makeGEOSRadiansDomainField(FlatField fltFld, double scaleX, double offsetX, double scaleY, double offsetY) Transforms FlatField with DomainCoordinateSystem (fgf_x, fgf_y) <-> (Lon,Lat) based on the Geostationary projection from fixed grid coordinates to intermediate coordinates view angle coordinates in radians (lambda, theta).
static float[]	minmax(float[] values, int length, int[] indexes)
private static float[][]	valueToGrid(Linear1DSet set, float[][] value) Adapted from VisAD except NaN is returned unless: val >= First and val <= Last, First > Last val <= First and val >= Last, Last > First i.e., clamp to inside the interval inclusive

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

LambdaTheta

public static RealTupleType LambdaTheta

uniqueID

private static int uniqueID

offsetTolerance

private static double offsetTolerance

accumTolerance

private static float accumTolerance

Constructor Detail

GEOSgridUtil

public GEOSgridUtil()

Method Detail

geosRGB

public static FlatField geosRGB(FlatField red,
                                FlatField green,
                                FlatField blue)
                         throws VisADException,
                                java.rmi.RemoteException

DomainSets of red, green, blue must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame. See makeGEOSRadiansDomainField(visad.FlatField, double, double, double, double). TargetSet is taken as that of the first component field (red).

Parameters:

red -

green -

blue -

Returns:

Throws:

VisADException

java.rmi.RemoteException

See Also:

makeGEOSRadiansDomainField(FlatField, double, double, double, double)

geosRGB

public static FlatField geosRGB(FlatField red,
                                FlatField green,
                                FlatField blue,
                                Linear2DSet targetSet)
                         throws VisADException,
                                java.rmi.RemoteException

Parameters:

red -

green -

blue -

targetSet -

Returns:

Throws:

VisADException

java.rmi.RemoteException

geosRGB

public static FlatField geosRGB(FlatField red,
                                FlatField green,
                                FlatField blue,
                                Linear2DSet targetSet,
                                boolean copy)
                         throws VisADException,
                                java.rmi.RemoteException

DomainSets of red, green, blue and targetSet must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame. See makeGEOSRadiansDomainField(FlatField, double, double, double, double).

Parameters:

red -

green -

blue -

targetSet - the target, or outgoing domain

copy - to copy range of resulting RGB FlatField (default=true)

Returns:

Throws:

VisADException

java.rmi.RemoteException

See Also:

makeGEOSRadiansDomainField(FlatField, double, double, double, double)

makeGEOSRadiansDomainField

public static FlatField makeGEOSRadiansDomainField(FlatField fltFld,
                                                   double scaleX,
                                                   double offsetX,
                                                   double scaleY,
                                                   double offsetY)
                                            throws VisADException,
                                                   java.rmi.RemoteException

Transforms FlatField with DomainCoordinateSystem (fgf_x, fgf_y) <-> (Lon,Lat) based on the Geostationary projection from fixed grid coordinates to intermediate coordinates view angle coordinates in radians (lambda, theta).

Parameters:

fltFld - The incoming FlatField

scaleX -

offsetX - To transform to fgf_x -> lambda (radians)

scaleY -

offsetY - To transform to fgf_y -> theta (radians)

Returns:

FlatField with transformed Domain and original (not copied) Range

Throws:

VisADException

java.rmi.RemoteException

geosResample

public static FlatField geosResample(FlatField fld,
                                     Linear2DSet targetSet)
                              throws VisADException,
                                     java.rmi.RemoteException

Efficiently upscales or downscales between ABI fields with domainSets on the Fixed Grid Frame. Can be seen as grid transfer process, to be used in conjunction with, not a replacement for, VisAD resample. DomainSets must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame. See makeGEOSRadiansDomainField(visad.FlatField, double, double, double, double).

Parameters:

fld -

targetSet -

Returns:

Result field with targetSet domain.

Throws:

VisADException

java.rmi.RemoteException

See Also:

makeGEOSRadiansDomainField(FlatField, double, double, double, double)

geosResample

public static FlatField geosResample(FlatField fld,
                                     Linear2DSet targetSet,
                                     int mode)
                              throws VisADException,
                                     java.rmi.RemoteException

Efficiently upscales or downscales between ABI fields with domainSets on the Fixed Grid Frame. Can be seen as grid transfer process, to be used in conjunction with, not a replacement for, VisAD resample. DomainSets must be (lambda, theta): the intermediate view angle coordinates (radians) for the Fixed Grid Frame. See makeGEOSRadiansDomainField(visad.FlatField, double, double, double, double).

Parameters:

fld -

targetSet -

mode - VisAD resample mode

Returns:

Result field with targetSet domain.

Throws:

VisADException

java.rmi.RemoteException

See Also:

makeGEOSRadiansDomainField(FlatField, double, double, double, double)

geosUpsample

static void geosUpsample(int[] yidxs,
                         int[] xidxs,
                         int lenY,
                         int lenX,
                         float[][] values,
                         int targetLenY,
                         int targetLenX,
                         float[][] targetValues,
                         int mode)

interp

static void interp(float[] yidxs,
                   float[] xidxs,
                   int lenY,
                   int lenX,
                   float[][] values,
                   int targetLenY,
                   int targetLenX,
                   float[][] targetValues,
                   int mode)

minmax

public static float[] minmax(float[] values,
                             int length,
                             int[] indexes)

valueToGrid

private static float[][] valueToGrid(Linear1DSet set,
                                     float[][] value)
                              throws VisADException

Adapted from VisAD except NaN is returned unless: val >= First and val <= Last, First > Last val <= First and val >= Last, Last > First i.e., clamp to inside the interval inclusive

Parameters:

set -

value -

Returns:

Throws:

VisADException