[ncl-talk] Connected points

Geeta Geeta geetag54 at yahoo.com
Thu Dec 1 20:00:28 MST 2016


Dear NCL Group. 
my purpose is to find out connected pixels and name them as TS1. name another set of pixels as TS2 and so on. the details are here as below
I have reflectivity data (3dim) which has been converted to 2 dim by taking the maxima in the vertical. then this data is again filtered based on another criterion and the final variable is maxZ_Conv which refers to the only convective pixels. This maxZ_Conv variable is a swath data with dimensions (9248x49). Now I have to find out the pixels which are connected. and name the group as TS1 or 1 or a. That means all the pixels which are connected/adjacent/or in the neighbourhood should be numbered as 1. 
similarly for the 2nd group, all the pixels in group2 shhould be numberd as 2.  
The criterion for connectivity is as follows:I have to search for 8 points around the reference and see if they are adjacent/or in the neighbourhood, I am calling them as connected. Similarly I have to do for all the pixels. 
I have converted the maxZ_Conv to binary (variable name is maxZ_Convb) but because of the FillValue, the code is not working.. 
Kindly suggest if i am correct logically and I have attatched my code. (0) PT4 2939 29 9.969209968386869e+36(0) PT5 2940 29 9.969209968386869e+36(0) PT6 2941 29 1(0) PT7 2939 30 9.969209968386869e+36(0) PT8 2940 30 1(0) PT9 2941 30 9.969209968386869e+36(0) PT1 2939 29 9.969209968386869e+36(0) PT2 2940 29 9.969209968386869e+36(0) PT3 2941 29 1(0) PT4 2939 30 9.969209968386869e+36(0) PT5 2940 30 1(0) PT6 2941 30 9.969209968386869e+36(0) PT7 2939 31 9.969209968386869e+36(0) PT8 2940 31 9.969209968386869e+36(0) PT9 2941 31 9.969209968386869e+36fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing functionfatal:["Execute.c":8578]:Execute: Error occurred at or near line 127 in file 2a25.080324.pix.count.ncl
I am also attaching the picture of the data interpreted by me. Kindly suggest how I should go about doing it. 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161202/dbc874bd/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2a25.080324.pix.count.ncl
Type: application/octet-stream
Size: 5782 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161202/dbc874bd/attachment-0001.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: WP_20161202_08_26_54_Pro.jpg
Type: image/jpeg
Size: 286713 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161202/dbc874bd/attachment-0001.jpg 


More information about the ncl-talk mailing list