Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  543 1632 (N16E50 L=170)  1378    212        3 ( 521/1523- 591/1726) -47.3 E58-E44 N09-N22
      2 1531 1668 (N51W16 L=236)   819    122        1 (1473/1633-1583/1706)  -4.4 W09-W22 N49-N53
      3  173 1187 (S17E66 L=154)  1675    197        2 ( 144/1100- 228/1268)  57.5 E69-E61 S22-S11
      4  372 1064 (S20E42 L=178)  1139    205        3 ( 284/1015- 461/1114)  29.7 E49-E35 S24-S16
      5  167  897 (S34E64 L=156)  1071    222        1 ( 105/ 817- 231/ 972)  21.4 E83-E54 S39-S28
      6 1714  905 (N08W45 L=265)   310     55        0 (1688/ 883-1741/ 927)  -9.5 W42-W48 N07-N08
      7  800  792 (S25E06 L=213)   327     75        1 ( 776/ 759- 830/ 828)  78.7 E07-E05 S27-S22
      8 1520  583 (S15W44 L=264)   322     78        1 (1483/ 567-1559/ 598)  18.4 W40-W46 S16-S13
      9 1624  527 (S14W55 L=275)  1381    241        1 (1559/ 457-1678/ 582) -84.3 W49-W60 S19-S09
     10 1515  475 (S21W50 L=270)   359    104        0 (1479/ 438-1543/ 516)  90.0 W48-W51 S23-S17
     11  806  415 (S48W07 L=227)  1601    358        1 ( 642/ 369- 952/ 498)   6.1 E09-W19 S51-S44
     12 1199  189 (S48W56 L=276)  1561    223        4 (1093/ 133-1269/ 235)   3.7 W42-W74 S50-S45
  

    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