<div dir="ltr">As Dennis says I don't think there is a simple way to do this. Possibly in Python but I'm not sure there. <div><br></div><div>I think what you're looking for is some type of flood fill algorithm with the fill values as your boundary. You could write a flood fill in fortran and wrap it in. Depends on your desire and ability for that kind of thing. The other option I can think of would be all ncl based using predefined functions.</div><div><br></div><div>1. Draw a faux plot of your data (maybe setting the fill value to a new plottable value depending on the next step). </div><div> In the above define your plot so the contour is the distinction between regions you want and regions you don't want. </div><div>2. Use get_isolines() to retrieve the coordinates outlining your good regions. </div><div>3. Use gc_inout() to determine which region your point is enclosed by. </div><div><br></div><div>You could pretty easily make the above a function where you pass the grid and the location and return the perimeter defined by get_isolines(). </div><div><br></div><div>Good luck with it. </div><div>Alan. </div><div><br></div><br><div class="gmail_quote"><div dir="ltr">On Thu, 16 Jul 2015 at 10:20 Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">"<span style="font-size:12.8000001907349px">Is there anything coded like this in NCL (or R or Python) already?"</span><div><span style="font-size:12.8000001907349px"><br></span></div></div><div dir="ltr"><div><span style="font-size:12.8000001907349px">To my knowledge, NCL does not have a function that does what you want.</span></div><div><span style="font-size:12.8000001907349px">I doubt R does either.</span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Good Luck</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jul 10, 2015 at 6:01 AM, Erik J. Schaffernicht <span dir="ltr"><<a href="mailto:eschaffe@uni-koeln.de" target="_blank">eschaffe@uni-koeln.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I have a 2d map where certain gridpoints have a values, all the rest are set to FillValue or Missing. The spatial structure of the points with values are not a-priori known.<br>
Now, I look for a function/procedure in NCL (or maybe in R or Python, if those are better suited) that can tell me all the points that are part of a region given that at the beginning of the algorithm I provide one gridpoint that should be member as well.<br>
<br>
The closure or border of my region is defined by all the gridpoints that are FillValue or Missing and direct neighbours of the points valid values. On the total map there might be several such separated regions, so it is important to know if that I only get the indices of the region that belongs to my initial given grid point.<br>
<br>
Is there anything coded like this in NCL (or R or Python) already?<br>
<br>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div></div>