Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  673 1328 (N01E28 L=359)   491    129        0 ( 622/1294- 732/1358)  -7.1 E31-E23 N00-N02
      2 1576 1173 (N17W28 L= 56)   325    100        0 (1539/1132-1603/1211) -21.8 W25-W31 N14-N18
      3 1482 1157 (N13W23 L= 50)   294     81        0 (1467/1110-1496/1194) -53.1 W21-W24 N10-N15
      4 1720 1016 (N13W43 L= 70)   398    113        1 (1684/ 960-1736/1059) -59.0 W39-W45 N09-N15
      5 1770  803 (N03W54 L= 81)   836    184        5 (1731/ 746-1819/ 871) -20.0 W48-W59 N00-N05
      6  512  772 (S34E27 L=  0)   759    165        3 ( 439/ 720- 574/ 809)  -9.5 E33-E21 S36-S32
      7 1649  469 (S18W62 L= 90)   383     21        0 (1630/ 449-1666/ 493) -63.4 W60-W64 S18-S16
      8  412  392 (S67E48 L=339)   699     67        1 ( 374/ 376- 453/ 406)   7.1 E61-E37 S68-S64
      9  545  361 (S65E21 L=  6)   390     43        1 ( 522/ 339- 567/ 376)   8.1 E24-E17 S67-S62
  

    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