List of Routines By Functionality | MISR Toolkit: Main Page
-Abstract
MTK_APPLY_REGRESSION Applies regression to given data.
-Copyright
Copyright (2013), California Institute of Technology.
U.S. Government sponsorship acknowledged.
-I/O
Given:
data1 input data
valid_mask1 valid mask for input data
mapinfo map info for input data
regression_coeff regression coefficients structure
regression_coeff_map_info map info for regression coefficients
the call:
status = MTK_APPLY_REGRESSION( data1, valid_mask1, mapinfo, regression_coeff, regression_coeff_map_info, regressed, regressed_mask )
returns:
status 0 on success; otherwise failure
regressed output data
regressed_mask valid mask for output data
-Examples
;;
;; Set up input parameters
;;
data1 = RANDOMU(seed, 512, 704)
valid_mask1 = byte(replicate(1,512,704))
data2 = RANDOMU(seed, 512, 704)
data2_sigma = RANDOMU(seed, 512, 704)
valid_mask2 = byte(replicate(1,512,704))
path = 37
resolution = 275
lat = 66.0
lon = -89.0
lat_extent = 1.5
lon_extent = 1.0
status = MTK_SETREGION_BY_LATLON_EXTENT( lat, lon, lat_extent, lon_extent, "degrees", region )
status = MTK_SNAP_TO_GRID( path, resolution, region, mapinfo )
size_factor = 2
status = MTK_REGRESSION_COEFF_CALC( data1, valid_mask1, data2, data2_sigma, valid_mask2, mapinfo, size_factor, regression_coeff, regression_coeff_map_info )
;;
;; The call
;;
status = MTK_APPLY_REGRESSION( data1, valid_mask1, mapinfo, regression_coeff, regression_coeff_map_info, regressed, regressed_mask )
;;
;; Output...
;;
print, '================================================='
index = WHERE(valid_mask1 EQ 1,count)
print, count
index = WHERE(regressed_mask EQ 1, count)
print, count
IDL outputs:
=================================================
360448
33792
-Particulars
None.
-Required Reading
For important details concerning this module's function, please refer to
the MTK routine MtkApplyRegression.c.
-Version
-IDL-MTK Version 1.2.4