Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1732 (N41E43 L=296)   429     54        0 ( 558/1710- 621/1751) -18.4 E47-E38 N39-N42
      2  488 1708 (N39E52 L=287)   582     62        0 ( 468/1669- 516/1739) -80.5 E54-E50 N35-N41
      3 1558 1567 (N34W38 L= 16)   340     54        0 (1540/1534-1573/1596)  78.7 W36-W39 N31-N36
      4  333 1549 (N27E59 L=280)   591     50        1 ( 306/1525- 356/1570) -12.5 E63-E54 N25-N28
      5  320  841 (S17E49 L=290)   449     96        0 ( 306/ 783- 335/ 886)  63.4 E49-E46 S20-S14
      6  213  761 (S22E62 L=277)   958     84        1 ( 163/ 736- 261/ 787)   3.8 E70-E55 S23-S20
      7 1650  519 (S31W54 L= 32)  1228    138        1 (1582/ 480-1699/ 555)  39.8 W47-W60 S33-S28
      8  518  510 (S39E38 L=301)   981    122        2 ( 462/ 477- 574/ 550)  14.0 E44-E32 S41-S35
      9 1500  452 (S37W43 L= 21)   379     64        1 (1487/ 416-1517/ 488) -68.2 W40-W45 S39-S34
     10 1187  440 (S41W17 L=355)   347    127        0 (1136/ 415-1242/ 477)  35.5 W13-W20 S42-S37
  

    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