Re^3: Contour mapping?

And no, this isn't about plotting; its about calculating and collating all these damned angles.

Actually, since plotting isn't your goal, I'm not sure you'll need the angle-sort I suggested... or maybe even the local maxima search-connected. It really depends: do you need to keep two different "islands" separate?

              6     g f e d
       7       6            c b
     7   7     6     7 7 7     a
     7 9   7   6   7       7
       7 7    6  7     8   7
             6     7 7   7 7 
                       7
[download]

If you have a grid as pictured, would you want to keep the two rings of height=7 separate, or would they all be the same "contour", in terms of the angles you need to calculate? And for calculating angles, assuming the alphabet ridge were all at the same height, would you need to know that a was adjacent to b, and b adjacent to c, but a not adjacent to c (etc)? Or would you just not care about ordering and adjacency, because you just have to compute so many annoying angles?

If adjacency doesn't matter, and you just need bands of approximately equal height, whether or not they are part of the same ring, then you don't need the local min/max searches, and you don't need the angle sorting
If adjacency doesn't matter, but you do need to treat the two rings of 7s separately, then you'll still have to do the min/max search, but not the angle sorting
If adjacency does matter and the rings need to be separate, then the procedure I outlined before is probably a good starting thought.

Comment on Re^3: Contour mapping? Download Code

Replies are listed 'Best First'.
Re^4: Contour mapping? by BrowserUk (Patriarch) on Apr 28, 2016 at 00:10 UTC
Or would you just not care about ordering and adjacency, because you just have to compute so many annoying angles? You're further down the road thinking about this than I am :) My assumption (perhaps desire is a better term here) is that the contours won't cross those of a higher or lower value (or band); despite that there may be other contours else where that are at the same level. My first thought reading your reply was that I had a much bigger problem on my hands that I first thought. Detecting whether the contours cross other already discovered contours; or those yet to be discovered would be a exponential explosion of processing and a total nightmare. Then I thought again; and (I think) that I can get away without that. Having followed some of your links; and some of the links they contained; a word in one of them -- otherwise not related -- stuck in my brain. The word was "strata". My current plan of approach is: Discover the min and max heights whilst reading the data; divide that range into a series of strata -- I'm using 100 which is arbitrary, but divides up the data into a nice bell-curve of buckets, with ~450 in the extreme edges and 100,000 in the middle bucket. Then, within each of those buckets I find the min&max X&Y, and use that to split the area in a (say) a 10 x 10 grid of rectangular regions. I then sort the points within each region by X then Y (or vice versa) to speed searching. Then I pick a point at random within a bucket, and copy it into a 'contour array'; and then look for the nearest point in the bucket to connect it to. The 10 x 10 spacial grid means that I will only have to look in the current region; or the 8 surrounding regions for the next nearest point. Unless they are empty in which case the next one over. That will also allow me to detect if one end runs off the edge of the data space. The search completes when I find myself connecting back to the first point; or when both ends have left the building. Repeat for all strata. As algorithm definitions go; its kinda sketchy and definitely incomplete; but I have something to explore, which is a sight more than I had a couple of hours ago. So thank you for the links and questions; they've both served to prompt my brain into gear. With the rise and rise of 'Social' network sites: 'Computers are making people easier to use everyday' Examine what is said, not who speaks -- Silence betokens consent -- Love the truth but pardon error. "Science is about questioning the status quo. Questioning authority". I knew I was on the right track :) In the absence of evidence, opinion is indistinguishable from prejudice.	[reply]

Replies are listed 'Best First'.

Re^4: Contour mapping?
by BrowserUk (Patriarch) on Apr 28, 2016 at 00:10 UTC

Or would you just not care about ordering and adjacency, because you just have to compute so many annoying angles?

You're further down the road thinking about this than I am :)

My assumption (perhaps desire is a better term here) is that the contours won't cross those of a higher or lower value (or band); despite that there may be other contours else where that are at the same level.

My first thought reading your reply was that I had a much bigger problem on my hands that I first thought. Detecting whether the contours cross other already discovered contours; or those yet to be discovered would be a exponential explosion of processing and a total nightmare.

Then I thought again; and (I think) that I can get away without that.

Having followed some of your links; and some of the links they contained; a word in one of them -- otherwise not related -- stuck in my brain. The word was "strata".

My current plan of approach is:

Discover the min and max heights whilst reading the data; divide that range into a series of strata -- I'm using 100 which is arbitrary, but divides up the data into a nice bell-curve of buckets, with ~450 in the extreme edges and 100,000 in the middle bucket.
Then, within each of those buckets I find the min&max X&Y, and use that to split the area in a (say) a 10 x 10 grid of rectangular regions.
I then sort the points within each region by X then Y (or vice versa) to speed searching.
Then I pick a point at random within a bucket, and copy it into a 'contour array'; and then look for the nearest point in the bucket to connect it to.
The 10 x 10 spacial grid means that I will only have to look in the current region; or the 8 surrounding regions for the next nearest point. Unless they are empty in which case the next one over.
That will also allow me to detect if one end runs off the edge of the data space.
The search completes when I find myself connecting back to the first point; or when both ends have left the building.
Repeat for all strata.

As algorithm definitions go; its kinda sketchy and definitely incomplete; but I have something to explore, which is a sight more than I had a couple of hours ago.

So thank you for the links and questions; they've both served to prompt my brain into gear.

With the rise and rise of 'Social' network sites: 'Computers are making people easier to use everyday'

Examine what is said, not who speaks -- Silence betokens consent -- Love the truth but pardon error.

"Science is about questioning the status quo. Questioning authority". I knew I was on the right track :)

In the absence of evidence, opinion is indistinguishable from prejudice.

[reply]