Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  714 1843 (N62E24 L=357)   583     52        1 ( 681/1815- 740/1864)   7.1 E29-E20 N58-N64
      2  822 1739 (N49E08 L= 13)   390     58        3 ( 801/1705- 843/1768) -80.5 E10-E06 N45-N51
      3 1174 1705 (N41W20 L= 42)   283     66        0 (1146/1688-1214/1717)   9.5 W17-W23 N40-N42
      4 1357 1524 (N26W28 L= 49)   979    224        2 (1240/1474-1442/1548) -26.6 W20-W33 N22-N28
      5 1702 1409 (N16W52 L= 73)  1179    240        3 (1647/1367-1748/1490)  53.1 W46-W58 N13-N21
      6 1058 1433 (N23W06 L= 27)   213     90        0 (1024/1421-1105/1450) -21.8 W04-W09 N22-N24
      7  702 1375 (N23E17 L=  4)  1000    335        3 ( 575/1296- 860/1428) -30.1 E26-E07 N16-N27
      8  352 1411 (N29E46 L=335)   375     76        0 ( 317/1391- 392/1427)   0.0 E50-E42 N28-N30
      9 1323 1276 (N11W21 L= 42)   196     60        1 (1297/1252-1350/1300) -56.3 W19-W22 N09-N12
     10 1266 1177 (N06W16 L= 37)   539    220        1 (1199/1134-1366/1237) -29.1 W12-W21 N02-N09
     11  375 1172 (N14E42 L=340)   279     95        1 ( 357/1128- 399/1222)  68.2 E43-E39 N11-N16
     12 1005 1104 (N04E00 L= 21)   189     50        0 ( 992/1077-1020/1138) -90.0 E01-W00 N02-N05
     13 1837  939 (S13W58 L= 79)   507     75        2 (1797/ 922-1869/ 964)   7.1 W53-W61 S14-S11
     14 1731  827 (S19W47 L= 69)   614     87        2 (1702/ 795-1766/ 867) -38.7 W44-W50 S21-S16
     15 1657  731 (S24W42 L= 63)   885    230        2 (1635/ 629-1687/ 838) -85.6 W39-W44 S31-S17
     16  820  579 (S26E18 L=  3)  1097    298        2 ( 710/ 477- 886/ 686)  42.9 E28-E12 S32-S19
  

    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