Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  551 1752 (N25E60 L=216)  1967    198        2 ( 490/1676- 601/1837) -47.1 E72-E56 N18-N32
      2  742 1753 (N29E44 L=232)   895    163        1 ( 710/1673- 771/1808) -47.7 E50-E39 N22-N33
      3  559 1672 (N19E52 L=224)  1721    239        4 ( 490/1601- 631/1736)  80.5 E55-E47 N13-N25
      4 1132 1640 (N32E12 L=264)   682    156        0 (1109/1566-1161/1713) -71.6 E14-E10 N26-N38
      5 1563 1491 (N37W22 L=298)   800    149        0 (1514/1453-1607/1546)  55.0 W18-W25 N33-N43
      6  380 1398 (S01E52 L=224)   479     73        0 ( 360/1368- 405/1439) -31.0 E54-E48 S02-N01
      7  219 1210 (S15E62 L=214)  2410    198        3 ( 166/1141- 281/1293)  54.0 E69-E56 S19-S08
      8  313 1220 (S13E52 L=224)   589    132        1 ( 275/1181- 340/1270)  63.4 E55-E49 S15-S09
      9  167 1054 (S25E66 L=210)   625    104        0 ( 131/1019- 214/1078)  18.4 E71-E59 S26-S22
     10 1795  678 (S02W63 L=339)   517     79        1 (1755/ 657-1830/ 698)  11.3 W57-W67 S02-N00
     11 1484  533 (S20W45 L=321)   484    141        1 (1429/ 493-1547/ 578)  11.3 W39-W49 S21-S17
     12 1184  538 (S29W25 L=302)   925    213        0 (1091/ 517-1299/ 561)  24.8 W19-W32 S31-S25
     13 1699  457 (S16W70 L=346)   806    108        1 (1645/ 431-1735/ 484)  32.0 W65-W74 S18-S13
  

    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