Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1594 1711 (N29W71 L= 43)   926     85        3 (1560/1679-1627/1745)  68.2 W67-W76 N26-N31
      2 1514 1651 (N29W56 L= 28)   562     52        2 (1490/1618-1542/1684)  63.4 W53-W59 N26-N30
      3 1053 1561 (N38W18 L=350)   445     95        1 (1005/1519-1089/1583) -71.6 W15-W20 N34-N40
      4 1637 1536 (N19W58 L= 30)  2240    364        5 (1465/1496-1768/1571) -20.8 W43-W73 N13-N24
      5  287 1374 (N44E54 L=278)   885    144        3 ( 217/1343- 336/1396) -11.3 E66-E46 N41-N45
      6  470 1308 (N36E32 L=300)   333     76        1 ( 444/1279- 502/1333)  21.8 E34-E28 N34-N37
      7 1461 1255 (N09W32 L=  4)   353    161        0 (1416/1180-1489/1316)  83.7 W30-W32 N04-N13
      8 1350  828 (S12W14 L=346)  1096    284        1 (1307/ 711-1387/ 972)  52.4 W09-W19 S17-S04
      9  695  700 (S04E27 L=305)   355    100        0 ( 673/ 655- 720/ 754)  45.0 E29-E24 S06-S01
     10 1036  378 (S32E17 L=315)  1706    326        5 ( 897/ 287-1129/ 492)  22.2 E30-E08 S36-S26
     11  540  408 (S19E50 L=283)   486     58        1 ( 509/ 386- 565/ 437)  26.6 E52-E46 S21-S17
     12  747  389 (S25E35 L=297)   350     78        0 ( 718/ 356- 776/ 431)  29.7 E38-E31 S26-S22
  

    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