Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1741 1554 (N21W73 L=138)   811     52        1 (1718/1524-1761/1582)   4.8 W68-W81 N18-N21
      2  191 1327 (N37E71 L=354)  1425    114        2 ( 156/1279- 238/1366) -12.8 E83-E60 N33-N39
      3 1483 1266 (N12W34 L= 98)   303     64        2 (1442/1247-1517/1286) -36.9 W30-W35 N10-N13
      4 1822 1016 (S11W55 L=120)   456     65        0 (1808/ 980-1839/1057) -90.0 W53-W57 S13-S08
      5 1914 1003 (S15W69 L=133)   857     80        2 (1887/ 975-1936/1044)  15.9 W63-W73 S16-S12
      6 1605  978 (S08W36 L=100)   828    158        6 (1543/ 930-1657/1036) -58.0 W30-W39 S11-S03
      7 1051  885 (S02E01 L= 64)   204    128        0 ( 986/ 861-1096/ 919) -38.7 E04-W01 S03-N01
      8  251  762 (N02E60 L=  4)   399     58        0 ( 232/ 736- 273/ 791)  31.0 E62-E57 N00-N03
      9  259  659 (S05E64 L=  1)  2264    371        5 ( 174/ 565- 333/ 780) -52.4 E72-E56 S10-N03
     10 1092  643 (S17E03 L= 62)   364    112        1 (1073/ 589-1108/ 699)  63.4 E04-E01 S20-S13
     11 1737  642 (S33W46 L=111)   692     78        1 (1691/ 626-1782/ 661) -15.3 W41-W52 S34-S31
     12  808  532 (S19E23 L= 41)   645    192        2 ( 767/ 463- 833/ 624)  48.0 E28-E19 S23-S13
     13 1360  537 (S30W13 L= 77)   567    123        0 (1327/ 481-1391/ 596)  49.4 W09-W15 S33-S26
     14  377  509 (S15E58 L=  7)   458     86        0 ( 359/ 465- 401/ 545)  32.0 E62-E53 S18-S12
     15  462  467 (S19E52 L= 13)   635    121        5 ( 408/ 441- 506/ 503)   0.0 E57-E46 S20-S16
     16  694  305 (S34E43 L= 22)  1185    215        3 ( 638/ 234- 745/ 400)  28.6 E55-E33 S39-S27
     17  601  360 (S28E46 L= 18)   481    108        1 ( 544/ 337- 654/ 378)   0.0 E52-E41 S29-S27
  

    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