Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1082 1864 (N55W50 L=272)   889     97        1 (1049/1826-1123/1900)   6.0 W41-W60 N53-N56
      2 1707 1582 (N15W70 L=292)  1982    197        2 (1632/1503-1769/1640) -79.7 W64-W75 N09-N20
      3  908 1468 (N34W06 L=229)  1420    451        1 ( 728/1425-1134/1502) -30.3 E06-W18 N25-N40
      4 1639 1417 (N09W51 L=273)   564    133        0 (1594/1377-1668/1466)  23.2 W46-W55 N07-N11
      5 1211 1243 (N13W16 L=239)  1296    505        2 ( 976/1139-1374/1306) -41.8 W03-W23 N02-N20
      6  744 1178 (N22E12 L=210)   749    192        2 ( 656/1146- 842/1209) -18.4 E18-E06 N19-N23
      7  276 1132 (N30E50 L=172)   333     38        1 ( 256/1102- 291/1158) -71.6 E52-E48 N28-N31
      8 1811 1115 (S12W54 L=277)   579     98        2 (1774/1083-1855/1151) -24.0 W50-W59 S15-S09
      9  760 1059 (N15E14 L=208)   589    239        1 ( 680/1003- 852/1087) -10.3 E20-E09 N11-N16
     10  248  895 (N17E57 L=166)   937    220        1 ( 163/ 842- 339/ 939) -24.0 E67-E48 N12-N20
     11 1580  823 (S21W28 L=250)   384    115        0 (1545/ 777-1618/ 867)  23.2 W24-W31 S22-S19
     12  156  723 (N08E77 L=145)  1411    100        0 ( 139/ 678- 180/ 778)  66.8 E83-E73 N04-N11
     13 1754  589 (S44W44 L=266)  2786    316        3 (1650/ 505-1839/ 725)  23.2 W31-W54 S47-S37
     14  766  441 (S21E33 L=189)   660    152        1 ( 743/ 359- 784/ 507)  55.0 E37-E29 S26-S16
     15 1491  406 (S47W14 L=237)   738    190        1 (1421/ 359-1563/ 462)  12.1 W07-W21 S49-S45
     16 1013  241 (S43E30 L=193)   399     50        0 ( 982/ 222-1040/ 262)   9.5 E33-E26 S43-S41
     17 1175  217 (S51E21 L=202)  1403    254        2 (1064/ 175-1286/ 248)  -7.1 E33-E09 S54-S48
  

    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