Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1435 1727 (N55W36 L= 79)   425     28        0 (1411/1709-1462/1741)   0.0 W32-W39 N53-N56
      2  465 1661 (N37E60 L=342)  1317    132        2 ( 406/1631- 527/1688)   2.9 E71-E50 N35-N39
      3  565 1543 (N30E41 L=  2)   370     56        1 ( 544/1520- 596/1564)  26.6 E43-E38 N28-N31
      4  892 1421 (N25E13 L= 30)  1286    320        3 ( 797/1345- 946/1544) -72.9 E18-E08 N20-N33
      5 1748 1357 (N28W55 L= 98)  4551    632        5 (1532/1171-1907/1518) -27.8 W35-W75 N17-N38
      6  335 1379 (N17E55 L=347)   355     57        1 ( 307/1356- 362/1398)   0.0 E58-E51 N15-N18
      7 1246 1321 (N22W12 L= 55)   310    121        0 (1217/1268-1293/1357)  56.3 W10-W15 N18-N24
      8 1781 1176 (N17W55 L= 98)   748    175        0 (1741/1092-1800/1234) -71.6 W50-W56 N11-N20
      9 1676 1146 (N14W44 L= 87)   506    157        0 (1623/1112-1734/1197) -16.7 W39-W49 N12-N17
     10  774  725 (S19E14 L= 29)   751    238        0 ( 741/ 638- 841/ 789)  78.7 E15-E09 S24-S14
     11 1722  641 (S17W56 L= 99)  1013    135        1 (1671/ 598-1777/ 692)  40.6 W52-W61 S19-S13
  

    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