Ken Mankoff writes:
> I've read the previous posts on inverse map projections and the
> lengthy tutorial by David Fanning, but still cannot get things to line
> up quite right. So I'm posting here for help...
>
> I have a data set (BEDMAP) with this information in the header:
>
> ncols 1371
> nrows 1371
> xllcorner -3426225.75
> yllcorner -3426225.75
> cellsize 5000
> NODATA_value -9999
>
> And this information on the website:
>
> Polar Stereographic projection with 71=B0S as the latitude of true scale
> and 0=B0E as the central meridian.
>
> I've managed to load the data, and inverse project it approximately
> such that things roughly line up. But I cannot get it accurate where
> my reference for 'accurate' is the /MAP_CONTINENTS, /HIRES keywords.
Oh, oh. This leads me to think you are using MAP_SET
to set up your map projections, instead of the more
accurate MAP_PROJ_INIT. After my mapping "epiphany" of
a couple of weeks ago, I have given up on the MAP_SET
projections entirely, except--possibly--for cartoon
maps.
And, of course, there is apparently a bug in the
"more accurate" MAP_PROJ_INIT routines, in that if
you use the UVBOX information in the map structure coming
directly from MAP_PROJ_INIT to set up your "data
coordinate space" for map overlays, you will still
be "slightly off". You need to use the UVBOX information
coming from MAP_PROJ_IMAGE for completely accurate
results. I've had a call into ITTVIS for three weeks
about this, but so far without results.
http://www.dfanning.com/map_tips/tiffoverlay.html
There could also be some confusion about whether the
re****ted corner pixel coordinates are in the center
of the pixel (likely) or on the edge of the pixel.
If it is the center, then you are going to have to
move the coordinates to the edge of the pixel, which
is what IDL needs.
http://www.dfanning.com/map_tips/precipmap.html
Let's see what you are doing. And can your provide
a link to an image?
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming (www.dfanning.com)
Se****e ma de ni thui. ("Perhaps thou speakest truth.")
|
