Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  498 1468 (N45E31 L=128)   401    118        1 ( 433/1446- 537/1490)   5.7 E37-E27 N43-N46
      2 1242 1456 (N28W24 L=183)   607    114        2 (1191/1417-1296/1490)  12.5 W20-W29 N25-N29
      3  486 1223 (N29E32 L=127)  1865    600        5 ( 351/1049- 710/1397) -45.0 E45-E18 N14-N42
      4  188 1211 (N32E60 L= 99)   597    100        2 ( 167/1189- 218/1230)   9.5 E62-E56 N30-N32
      5  199 1117 (N26E59 L=100)   889    103        1 ( 172/1081- 234/1153)  76.0 E61-E55 N23-N27
      6 1656 1090 (S04W41 L=200)   882    164        3 (1591/1059-1710/1122) -18.4 W36-W45 S06-S01
      7 1492 1000 (S04W28 L=186)   278     58        2 (1474/ 963-1507/1030)  76.0 W26-W29 S06-S02
      8  262  983 (N17E54 L=105)   347     80        0 ( 246/ 944- 275/1025)  59.0 E56-E52 N14-N19
      9  813  903 (N04E15 L=144)  1195    248        5 ( 732/ 819- 890/1008) -63.4 E17-E11 S01-N11
     10 1517  821 (S16W26 L=185)   476     61        1 (1501/ 781-1535/ 865) -90.0 W24-W27 S18-S13
     11  804  573 (S16E23 L=136)   972    207        5 ( 716/ 508- 859/ 621) -53.1 E27-E20 S20-S11
     12  666  276 (S35E47 L=112)   994    168        1 ( 629/ 237- 694/ 316)  28.6 E53-E42 S37-S31
  

    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