Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  356 1515 (N10E58 L=132)   485     82        1 ( 330/1483- 372/1559) -78.7 E59-E55 N07-N12
      2  304 1438 (N04E59 L=131)  2241    415        9 ( 207/1377- 405/1493)  17.5 E71-E51 N00-N08
      3  903 1388 (N11E15 L=175)   606    213        2 ( 863/1284- 950/1454) -73.3 E17-E11 N05-N15
      4  461 1400 (N04E44 L=146)   304     89        2 ( 432/1365- 506/1429)  59.0 E46-E41 N01-N06
      5  126 1160 (S14E70 L=120)  3340    319        8 (  93/1059- 165/1247)  90.0 E77-E65 S19-S08
      6  345  925 (S26E43 L=147)  2178    521        6 ( 235/ 768- 481/1100)  53.5 E53-E34 S36-S13
      7 1296  885 (S09W19 L=209)   386     92        2 (1242/ 862-1337/ 907)   9.5 W15-W21 S10-S07
      8 1295  741 (S17W23 L=213)  1000    317        1 (1156/ 640-1401/ 855)  64.8 W16-W26 S25-S08
      9  943  731 (S26W02 L=192)   600    118        1 ( 882/ 685-1005/ 776)  63.4 E01-W04 S29-S21
     10  880  585 (S36W02 L=192)   241     90        0 ( 825/ 564- 930/ 601)   8.1 E02-W05 S36-S35
     11 1562  462 (S25W52 L=242)  2785    653        0 (1305/ 356-1809/ 544)  38.2 W39-W73 S38-S12
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y