I'm using NLDAS (North America Land Data Assimilation System) data which has 0.125 degree resolution with 464 X 224 grid boxes, and I'm trying to project this data to the Albers map projection using the code below. The LIMIT key used in the MAP_PROJ_IMAGE indicates the total area coverage of NLDAS. Using this code, the result returns with 464 X 224 grid boxes which are the exact same size as the original data input. How can I know the spatial resolution of the result. I believe the result should have a spaptial resolution in meters, right? Is there any way to give a specific resolution when the original data are wapred?
sMap = MAP_PROJ_INIT(103, SEMIMAJOR_AXIS=6378137.0, SEMIMINOR_AXIS=6356752.31414,$
FALSE_EASTING=0., FALSE_NORTHING=0, STANDARD_PAR1=29.5,$
STANDARD_PAR2=45.5, CENTER_LATITUDE=23., CENTER_LONGITUDE=-96.)
Result = MAP_PROJ_IMAGE(Image, [-125.,25.,-67.,53.], Map_structure=sMap, missing=0)
Thanks.
|