Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1471 1709 (N31W54 L=105)   387     72        2 (1456/1669-1495/1742)  26.6 W50-W58 N28-N32
      2 1624 1543 (N17W55 L=105)  1851    372        2 (1482/1484-1761/1628) -41.4 W47-W66 N10-N25
      3 1299 1535 (N25W31 L= 81)  2281    611        4 (1037/1434-1557/1589) -28.7 W13-W44 N13-N32
      4 1218 1401 (N20W21 L= 71)   313    136        0 (1157/1384-1287/1420) -29.7 W17-W24 N17-N21
      5 1853 1249 (S05W62 L=113)   436     66        1 (1836/1211-1865/1281)  45.0 W60-W64 S06-S03
      6  246 1093 (N24E52 L=358)   802    127        2 ( 179/1059- 301/1127) -12.1 E60-E46 N21-N26
      7 1556 1073 (S07W32 L= 82)   242     29        1 (1531/1055-1579/1095) -18.4 W30-W33 S07-S05
      8 1564  993 (S11W31 L= 81)   416      0        0 (1499/ 969-1613/1031)  -7.1 W26-W34 S12-S09
      9 1104  837 (S09W00 L= 50)   537    119        3 (1061/ 794-1147/ 895) -80.5 E02-W02 S11-S05
     10  678  765 (S03E26 L= 25)  2040    560        5 ( 516/ 594- 866/ 894) -53.6 E34-E18 S16-N05
     11  217  841 (N10E59 L=351)   741    125        2 ( 167/ 806- 278/ 878) -24.4 E64-E53 N07-N13
     12  972  767 (S10E09 L= 42)   246     67        1 ( 935/ 749-1001/ 793)  11.3 E11-E06 S10-S08
     13 1216  594 (S26W01 L= 51)  1387    376        4 (1113/ 539-1350/ 717) -49.8 E04-W10 S29-S15
     14  412  615 (S06E49 L=  1)   470     59        3 ( 378/ 592- 447/ 637) -36.9 E51-E45 S07-S04
     15  963  337 (S36E24 L= 26)  3862    608        5 ( 822/ 152-1080/ 600)  43.2 E48-E13 S48-S18
     16  390  550 (S10E54 L=357)   654    101        3 ( 344/ 515- 423/ 586) -51.3 E57-E49 S12-S07
     17  557  551 (S13E41 L= 10)   385     44        3 ( 525/ 526- 581/ 574) -56.3 E43-E38 S14-S11
     18  459  463 (S16E53 L=358)  2034    251        2 ( 405/ 339- 505/ 566)  33.7 E69-E46 S23-S09
  

    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