Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1738 1548 (N26W73 L=138)   997    149        0 (1683/1512-1794/1595) -38.7 W68-W80 N22-N30
      2 1686 1439 (N21W57 L=122)   855    151        2 (1620/1393-1742/1480) -37.9 W51-W60 N17-N24
      3  875 1340 (N27E06 L= 60)   876    430        0 ( 735/1248- 994/1395) -17.4 E16-W01 N19-N30
      4  359 1143 (N19E46 L= 20)   422     71        3 ( 331/1110- 393/1173)  63.4 E48-E42 N16-N20
      5 1430  935 (S05W24 L= 90)   890    258        5 (1362/ 834-1475/1051) -77.9 W21-W26 S11-N03
      6 1798  976 (S08W55 L=120)   483    107        2 (1782/ 928-1818/1030) -74.1 W53-W56 S11-S04
      7 1751  855 (S15W49 L=115)   618     27        0 (1726/ 823-1776/ 885)  90.0 W46-W52 S17-S12
      8  472  810 (S02E39 L= 27)   507    114        2 ( 430/ 773- 529/ 840) -29.7 E41-E34 S04-N00
      9  917  623 (S18E12 L= 54)   424    172        3 ( 887/ 545- 946/ 682)  60.9 E14-E08 S23-S14
     10  262  538 (S20E73 L=352)  3058    166        1 ( 211/ 465- 300/ 599)  84.3 E83-E69 S25-S15
     11 1256  177 (S64W10 L= 76)   735     87        0 (1211/ 161-1293/ 196)  -4.4 W03-W16 S66-S61
  

    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