Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1173 1761 (N44E17 L=196)   575    159        0 (1110/1699-1221/1821) -20.6 E26-E10 N40-N47
      2  327 1446 (N00E59 L=154)  1008     82        2 ( 301/1399- 359/1482)  90.0 E63-E56 S02-N03
      3 1187 1371 (N18E01 L=213)   802    186        0 (1111/1295-1253/1426) -16.7 E06-W04 N14-N19
      4  363 1330 (S06E51 L=163)   349     36        0 ( 342/1304- 384/1358)  90.0 E52-E48 S07-S03
      5  214 1068 (S24E59 L=155)   444     58        1 ( 197/1038- 231/1100)  90.0 E61-E56 S25-S21
      6  774  907 (S20E11 L=203)  1245    444        0 ( 585/ 815- 962/ 964)  42.0 E21-E01 S29-S11
      7  399  802 (S36E37 L=177)   755    170        1 ( 366/ 725- 459/ 868)  61.4 E39-E31 S41-S30
      8  900  750 (S25W01 L=215)   600    182        0 ( 830/ 706- 992/ 779)  52.1 E02-W05 S29-S20
      9  973  624 (S30W09 L=223)   926    290        1 ( 837/ 606-1102/ 649)  31.0 W00-W16 S34-S24
     10 1257  499 (S28W32 L=246)   307     84        0 (1219/ 478-1298/ 524)   8.1 W28-W35 S28-S26
     11 1562  456 (S20W56 L=270)   990    138        1 (1496/ 417-1630/ 499)  56.3 W52-W59 S24-S15
     12  413  448 (S62E40 L=174)   472     48        1 ( 395/ 428- 436/ 475) -36.9 E44-E34 S63-S59
  

    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