Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1691 (N36E15 L=112)  1761    408        7 ( 972/1584-1235/1794) -15.1 E29-E03 N30-N39
      2  464 1634 (N15E56 L= 71)   885    169        0 ( 383/1586- 544/1689)  48.4 E61-E51 N11-N20
      3  589 1611 (N16E44 L= 82)   475    118        3 ( 554/1581- 622/1648) -23.2 E48-E40 N14-N18
      4 1202 1506 (N25E03 L=123)   277     66        1 (1186/1477-1214/1538) -76.0 E04-E01 N23-N27
      5 1895  971 (N17W61 L=188)   462    101        1 (1881/ 931-1905/1024) -80.5 W59-W63 N13-N19
      6 1061  964 (S09W04 L=130)   404    159        0 (1032/ 899-1108/1008) -63.4 W02-W06 S12-S06
      7 1807  845 (N06W56 L=182)   423     56        1 (1790/ 817-1824/ 875) -56.3 W54-W57 N04-N07
      8 1166  712 (S20W17 L=143)   217     63        2 (1129/ 691-1197/ 730)  14.0 W14-W18 S21-S19
      9 1614  497 (S18W56 L=182)  1010    163        1 (1570/ 418-1659/ 563) -78.7 W53-W57 S22-S12
     10 1687  484 (S16W63 L=190)  1331    194        2 (1630/ 434-1730/ 538)  71.6 W58-W66 S20-S11
     11 1353  502 (S26W37 L=163)   395    144        1 (1290/ 483-1423/ 519)  24.0 W32-W41 S27-S23
     12 1432  214 (S38W62 L=189)   748    197        2 (1393/ 196-1468/ 247)  26.6 W59-W65 S40-S36
  

    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