Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1785 (N35E61 L=354)  1006     94        1 ( 514/1738- 559/1830) -38.7 E69-E55 N31-N40
      2  791 1680 (N30E30 L= 25)   301     71        0 ( 751/1663- 822/1698)  36.9 E31-E27 N28-N32
      3  814 1495 (N18E22 L= 34)  1313    336        3 ( 735/1410- 931/1634) -32.3 E31-E12 N14-N26
      4  342 1483 (N12E57 L=358)   486    121        0 ( 305/1445- 385/1514) -20.0 E62-E51 N09-N13
      5  411 1195 (S05E42 L= 13)   624    170        0 ( 354/1149- 497/1250)  37.9 E45-E36 S08-S01
      6  939 1118 (S03E06 L= 49)   340    180        0 ( 907/1065- 963/1197) -74.1 E09-E04 S06-N01
      7  421 1014 (S16E39 L= 16)   365    109        1 ( 374/ 988- 468/1037)  -7.1 E42-E34 S17-S14
      8  812  994 (S12E12 L= 43)   184     45        1 ( 778/ 979- 838/1008)  14.0 E14-E10 S12-S10
      9  554  910 (S21E28 L= 27)   227     83        0 ( 542/ 871- 566/ 954) -78.7 E28-E26 S22-S17
     10  911  745 (S25E02 L= 53)   534    290        0 ( 824/ 644- 998/ 847)  68.2 E06-W01 S32-S17
     11 1463  612 (S22W37 L= 92)   375    194        0 (1413/ 551-1496/ 660)  90.0 W33-W39 S26-S19
     12 1060  602 (S31W11 L= 66)   409    131        1 (1023/ 553-1111/ 650)  69.4 W08-W12 S34-S26
  

    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