Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  475 1668 (N35E55 L=360)  2107    513        6 ( 434/1591- 535/1745) -24.2 E65-E45 N30-N39
      2  918 1580 (N35E16 L= 39)   405     93        2 ( 874/1557- 963/1610) -15.9 E19-E12 N33-N36
      3 1410 1499 (N36W22 L= 77)  1361    258        5 (1348/1386-1469/1604) -60.3 W17-W26 N28-N43
      4  717 1381 (N19E25 L= 30)  1878    424        4 ( 588/1239- 811/1516) -33.7 E37-E16 N10-N26
      5 1101 1433 (N27E00 L= 55)  1006    218        2 (1030/1334-1188/1515) -39.8 E06-W06 N21-N31
      6 1692 1119 (N14W44 L= 99)  1473    348        4 (1651/ 970-1727/1300)  90.0 W41-W47 N04-N25
      7  248 1227 (N03E59 L=356)   408     91        0 ( 236/1184- 258/1279) -80.5 E60-E57 N00-N06
      8 1520 1195 (N17W30 L= 85)   762    268        0 (1450/1135-1560/1261)  48.8 W25-W33 N12-N21
      9 1319 1141 (N11W17 L= 72)   257    101        1 (1271/1105-1363/1182) -29.7 W13-W20 N09-N13
     10 1893 1070 (N13W67 L=122)  3278    355        6 (1855/ 984-1931/1155)  70.0 W61-W73 N07-N18
     11 1098 1064 (N04W04 L= 59)   566    166        2 (1015/1041-1180/1092)  26.6 W00-W08 N01-N06
     12 1811 1015 (N09W57 L=112)   473    130        1 (1778/ 957-1832/1062) -45.0 W52-W59 N05-N11
     13  466  941 (S11E35 L= 20)  1274    458        0 ( 307/ 905- 647/ 998)   0.0 E49-E22 S13-S08
     14 1192  846 (S08W12 L= 67)   283    129        0 (1155/ 782-1232/ 909)  77.5 W10-W14 S12-S03
     15  315  814 (S21E48 L=  7)   325     60        0 ( 299/ 786- 341/ 844) -45.0 E49-E45 S22-S19
     16  255  769 (S25E54 L=  1)   480    111        1 ( 225/ 749- 281/ 789)  26.6 E57-E51 S26-S23
     17 1730  708 (S10W55 L=110)   771    198        2 (1702/ 667-1782/ 764) -16.7 W50-W60 S12-S06
     18 1436  369 (S37W43 L= 98)   963    183        3 (1365/ 345-1501/ 402)  21.8 W37-W48 S39-S33
  

    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