Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  622 1688 (N61E10 L= 45)  1642    246        2 ( 516/1631- 724/1738)  -6.3 E25-W02 N57-N64
      2 1230 1621 (N34W34 L= 89)   287    107        0 (1182/1599-1273/1654) -71.6 W32-W36 N31-N37
      3  478 1617 (N60E31 L= 24)   347     43        0 ( 462/1595- 499/1641) -39.8 E33-E26 N57-N62
      4  400 1553 (N57E41 L= 14)   388     59        0 ( 389/1519- 415/1586) -78.7 E43-E38 N54-N59
      5  731 1495 (N42E05 L= 49)   598    159        1 ( 685/1426- 786/1552) -76.0 E07-E03 N35-N47
      6 1148 1435 (N26W20 L= 74)   482    196        0 (1058/1398-1221/1450) -41.2 W14-W23 N21-N28
      7  806 1268 (N26E06 L= 48)   964    312        6 ( 691/1219- 928/1344) -28.3 E13-W01 N21-N30
      8  820 1186 (N21E07 L= 47)   380    127        1 ( 774/1136- 845/1217)  53.1 E10-E05 N17-N22
      9  956 1179 (N17W00 L= 55)   270     73        1 ( 926/1142- 977/1212)  76.0 E01-W01 N14-N18
     10  630 1090 (N20E22 L= 33)   478    160        1 ( 565/1057- 709/1121) -36.9 E25-E17 N16-N22
     11 1802  909 (S23W49 L=103)  1398    317        0 (1696/ 828-1873/1035)  12.8 W37-W59 S27-S18
     12 1645  883 (S19W34 L= 89)   558    120        2 (1612/ 825-1682/ 926) -74.1 W31-W36 S23-S16
     13 1564  790 (S22W26 L= 81)   477    158        0 (1518/ 732-1613/ 856)  29.1 W21-W30 S24-S19
     14 1551  637 (S32W22 L= 77)   283     82        0 (1513/ 621-1587/ 654)  -9.5 W19-W25 S32-S29
     15 1089  559 (S22E09 L= 45)   547    133        1 (1043/ 501-1134/ 614)  32.0 E13-E05 S24-S19
     16  945  368 (S30E26 L= 29)   284     67        0 ( 910/ 349- 973/ 393)   9.5 E28-E22 S31-S29
     17  902  191 (S43E43 L= 12)   674    102        1 ( 868/ 152- 938/ 229)  18.4 E51-E36 S46-S40
  

    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