Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1639 (N30W39 L= 76)  1261    244        2 (1233/1584-1452/1674)  -2.9 W29-W49 N28-N31
      2 1672 1620 (N20W65 L=102)  1572    185        1 (1599/1563-1743/1673) -55.0 W60-W71 N14-N24
      3 1576 1583 (N21W53 L= 90)  2154    342        5 (1413/1487-1711/1639) -50.5 W43-W60 N11-N28
      4 1353 1492 (N21W32 L= 69)   261    108        0 (1298/1481-1404/1503) -23.2 W28-W35 N19-N22
      5 1724 1110 (S09W45 L= 82)  1258    305        2 (1661/1034-1767/1224)  70.0 W39-W49 S12-S01
      6  422 1123 (N22E36 L=  1)   701    168        1 ( 342/1067- 469/1159) -41.6 E41-E32 N17-N25
      7 1407  957 (S09W20 L= 57)   279    119        0 (1386/ 905-1425/1009)  68.2 W18-W21 S11-S06
      8  868  868 (S02E12 L= 25)  1798    441        7 ( 724/ 694- 964/ 980) -70.6 E19-E09 S13-N05
      9 1335  924 (S09W15 L= 52)   431    102        2 (1299/ 884-1358/ 975)  45.0 W12-W17 S11-S06
     10 1189  892 (S08W06 L= 43)   427    156        2 (1153/ 833-1210/ 950)  56.3 W02-W08 S10-S03
     11  348  899 (N11E45 L=351)   699    158        2 ( 287/ 865- 425/ 945) -31.0 E49-E39 N07-N14
     12 1241  916 (S08W10 L= 46)   251     36        1 (1219/ 888-1264/ 941)  45.0 W08-W11 S09-S06
     13 1442  662 (S28W17 L= 54)  1288    249        6 (1334/ 625-1544/ 713) -11.3 W09-W24 S29-S24
     14  559  669 (S06E36 L=  1)   305     38        0 ( 536/ 641- 581/ 690) -90.0 E37-E34 S08-S04
     15 1130  406 (S36E10 L= 27)  3582    665        7 ( 969/ 185-1228/ 681)  42.3 E34-W00 S49-S18
     16  542  591 (S10E40 L=357)   703    110        3 ( 486/ 554- 581/ 626) -45.0 E43-E37 S13-S08
     17  527  422 (S20E50 L=347)  2778    467        4 ( 384/ 317- 629/ 564)  14.3 E79-E36 S24-S13
     18 1045  329 (S39E19 L= 18)   292     72        0 (1031/ 290-1057/ 369)  50.2 E21-E16 S42-S36
  

    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