Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  980 1738 (N39E17 L=136)   424     91        0 ( 931/1725-1036/1754)  23.2 E20-E13 N37-N40
      2  679 1620 (N25E33 L=120)   497    123        1 ( 642/1583- 719/1660) -26.6 E37-E29 N23-N27
      3  431 1597 (N21E52 L=102)  2134    292        3 ( 355/1552- 516/1637)  -6.3 E60-E42 N17-N23
      4  292 1483 (N13E61 L= 93)  1561    332        6 ( 230/1436- 356/1539)  33.7 E67-E55 N09-N16
      5  487 1402 (N09E40 L=113)  1124    308        2 ( 468/1264- 519/1529) -63.4 E43-E34 N01-N17
      6  365 1414 (N09E50 L=103)   830    128        1 ( 318/1361- 415/1463)  54.5 E54-E47 N05-N12
      7 1084 1199 (N05W01 L=154)   390    128        2 (1026/1141-1112/1245) -31.0 E02-W03 N01-N07
      8  260 1010 (S16E52 L=101)  1065    314        1 ( 215/ 885- 349/1129)  59.0 E55-E45 S23-S08
      9 1743 1031 (N07W46 L=200)   383    114        2 (1701/ 995-1769/1062) -45.0 W43-W49 N05-N08
     10 1234  936 (S08W14 L=167)   377    165        0 (1161/ 873-1285/ 962) -14.0 W09-W17 S11-S06
     11 1798  778 (S06W58 L=211)   403     95        1 (1778/ 723-1821/ 819) -71.6 W56-W59 S09-S03
     12 1610  768 (S11W42 L=195)   257     71        2 (1576/ 736-1643/ 801)  45.0 W39-W43 S12-S08
     13 1363  766 (S16W25 L=178)   606    177        0 (1274/ 742-1450/ 799)  24.4 W19-W30 S18-S12
     14  487  709 (S33E32 L=121)   312     82        0 ( 455/ 687- 524/ 727)  31.0 E34-E29 S34-S30
     15  329  571 (S44E49 L=104)   931    167        2 ( 295/ 524- 372/ 613)  41.2 E54-E43 S47-S40
     16  487  493 (S48E35 L=118)  1079    238        3 ( 398/ 473- 568/ 528)   3.4 E45-E27 S49-S45
     17  861  204 (S67W07 L=161)   551     30        0 ( 835/ 184- 885/ 227)  18.4 W03-W10 S69-S63
  

    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