Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1009 1811 (N42E28 L=174)   817    184        3 ( 969/1764-1045/1869) -30.3 E34-E22 N37-N46
      2  445 1625 (N15E56 L=146)   590    145        1 ( 427/1576- 461/1668) -54.5 E60-E53 N11-N18
      3  952 1356 (N09E12 L=189)   642    315        2 ( 914/1225- 996/1491) -83.3 E15-E09 N01-N18
      4 1814 1219 (N27W47 L=249)   494    124        1 (1778/1163-1843/1276) -35.5 W43-W50 N24-N30
      5 1136 1169 (N04W03 L=204)   584    211        2 (1031/1126-1229/1193)  15.3 E03-W08 N01-N05
      6 1267 1013 (S02W14 L=216)   773    238        3 (1226/ 897-1295/1124) -51.3 W09-W17 S07-N04
      7 1120  909 (S11W08 L=210)   638    212        4 (1041/ 857-1217/ 953)  48.8 W04-W12 S15-S06
      8  815  892 (S20E09 L=193)   318    126        1 ( 785/ 863- 881/ 933)   0.0 E11-E04 S21-S17
      9 1780  872 (N06W52 L=254)   498    138        1 (1736/ 849-1828/ 899)   0.0 W48-W57 N04-N07
     10  788  578 (S39E02 L=199)  1694    446        8 ( 604/ 551-1011/ 609)  17.2 E17-W12 S43-S33
     11 1372  529 (S25W36 L=238)  1419    425        3 (1191/ 500-1555/ 568)  21.8 W24-W49 S30-S18
     12 1513  465 (S24W49 L=251)   305     66        1 (1482/ 435-1536/ 492)  76.0 W47-W50 S25-S21
  

    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