MAPS-L Archives

Maps-L: Map Librarians, etc.

MAPS-L@LISTSERV.UGA.EDU

Options: Use Forum View

Use Monospaced Font
Show Text Part by Default
Show All Mail Headers

Message: [<< First] [< Prev] [Next >] [Last >>]
Topic: [<< First] [< Prev] [Next >] [Last >>]
Author: [<< First] [< Prev] [Next >] [Last >>]

Print Reply
Subject:
From:
Johnnie Sutherland <[log in to unmask]>
Reply To:
Maps and Air Photo Systems Forum <[log in to unmask]>
Date:
Fri, 27 Aug 1999 13:45:55 -0400
Content-Type:
TEXT/PLAIN
Parts/Attachments:
TEXT/PLAIN (32 lines)
---------- Forwarded message ----------
Date: Sun, 27 Aug 2000 02:20:41 CDT
From: Joseph McCollum <[log in to unmask]>
To: [log in to unmask]
Subject: Re: Reference Question Regarding Elevations (fwd)

Get the DEM from the National Geophysical Data Center.
Also they have a county mask.

It would require some coding, but basically it is this:

high = zeros(3200);     /* there are something like 3200 counties in
the USA */
low = 10000 * ones(3200);
all = zeros(3200);
count = zeros(3200);
for i = 1:rows,
   for j = 1:cols,
       c = county(i,j);
       e = elev(i,j);
       if e > high(c);  high(c) = e;    endif;
       if e < low(c);   low(c) = e;     endif;
       all(c) = sum(all(c),e);
        count(c) = sum(count(c),1);
   end;
end;
avg = all/count;

This would not get the actual high elevation (or actual low) elevation in
each county.    It would get the highest (or lowest) - valued DEM pixel
(itself an average).

ATOM RSS1 RSS2