Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1587 1673 (N24W65 L=276)   579     75        1 (1547/1649-1620/1696) -51.3 W62-W68 N21-N26
      2 1408 1566 (N25W42 L=252)   318    108        0 (1388/1518-1431/1614)  56.3 W40-W44 N21-N27
      3 1056 1547 (N35W19 L=229)  1747    517        4 ( 846/1503-1303/1598) -29.7 W03-W33 N24-N40
      4 1736 1481 (N09W64 L=274)   524     73        0 (1712/1451-1756/1516)   6.3 W59-W68 N07-N11
      5 1410 1342 (N13W32 L=242)   967    374        3 (1236/1262-1533/1397) -43.0 W21-W37 N04-N18
      6  909 1258 (N22E00 L=210)   756    227        1 ( 812/1226-1015/1298) -14.0 E06-W06 N20-N24
      7  398 1191 (N31E38 L=173)   302     38        0 ( 382/1159- 416/1214)  71.6 E38-E35 N29-N32
      8 1891 1174 (S12W66 L=276)  1060     98        2 (1861/1141-1928/1211) -32.0 W62-W73 S14-S09
      9  903 1141 (N16E04 L=206)   299    190        0 ( 860/1071- 977/1175)  14.0 E07-W01 N12-N16
     10  341  948 (N18E46 L=164)   592    177        2 ( 276/ 930- 425/ 972) -17.1 E52-E39 N16-N20
     11 1498  856 (S16W23 L=233)   252     42        1 (1478/ 828-1514/ 888)  45.0 W21-W24 S17-S14
     12  227  738 (N08E64 L=146)  1290    134        4 ( 201/ 686- 262/ 801)  74.1 E68-E60 N04-N11
     13  317  706 (N05E55 L=155)   446     84        2 ( 287/ 673- 347/ 745) -59.0 E56-E53 N02-N07
     14  931  536 (S19E19 L=191)   573    178        3 ( 910/ 455- 954/ 617)  63.4 E22-E16 S24-S14
     15 1786  585 (S46W50 L=260)   951     84        1 (1741/ 562-1828/ 609)  -7.6 W42-W58 S48-S43
     16 1600  463 (S47W26 L=236)   724    150        1 (1525/ 429-1668/ 499)  -3.6 W17-W34 S49-S45
     17 1257  255 (S51E12 L=198)  1031    176        1 (1176/ 222-1338/ 300) -14.9 E19-E04 S53-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