Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1881 (N60E36 L=118)  1066    113        1 ( 944/1847-1024/1925) -11.8 E54-E28 N56-N65
      2  663 1577 (N27E37 L=117)   508    194        0 ( 619/1549- 728/1635)  -5.7 E42-E32 N25-N30
      3  327 1548 (N20E68 L= 86)  1080    100        0 ( 279/1525- 377/1570)   7.6 E77-E60 N18-N20
      4  565 1462 (N18E39 L=115)   454     68        2 ( 545/1428- 590/1496) -36.9 E42-E36 N16-N20
      5 1578 1236 (N22W33 L=188)   506    165        0 (1563/1149-1592/1310) -78.7 W31-W34 N16-N26
      6 1880 1279 (N31W70 L=225)   903     95        0 (1842/1250-1916/1302)  -3.4 W62-W81 N29-N31
      7  246  849 (S24E56 L= 99)   731    142        0 ( 226/ 792- 292/ 903)  48.8 E58-E50 S28-S20
      8  386  797 (S25E41 L=113)   311     92        0 ( 353/ 781- 432/ 824)   0.0 E44-E36 S26-S23
      9 1086  622 (S24W11 L=166)  2139    559        5 ( 874/ 492-1288/ 754)  50.2 W00-W21 S36-S12
     10 1709  646 (S11W56 L=211)   465     38        2 (1688/ 623-1731/ 674) -45.0 W54-W59 S12-S09
     11 1558  617 (S16W44 L=198)   601    146        1 (1508/ 569-1624/ 657)  49.4 W41-W48 S19-S12
     12  609  255 (S67E31 L=123)   573     54        1 ( 584/ 240- 640/ 278)  16.4 E41-E24 S70-S63
  

    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