Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1438 1728 (N55W36 L= 78)   404     39        1 (1412/1714-1464/1743)   0.0 W32-W40 N53-N56
      2  468 1659 (N37E59 L=343)  1256    125        3 ( 415/1627- 533/1687)  -5.7 E70-E50 N35-N38
      3  573 1537 (N30E40 L=  2)   456     80        1 ( 550/1497- 604/1562)  90.0 E42-E37 N27-N32
      4  899 1420 (N25E12 L= 30)  1204    325        2 ( 805/1345- 949/1543) -69.0 E17-E08 N20-N33
      5 1730 1379 (N30W54 L= 96)  3414    515        2 (1537/1206-1867/1517) -32.7 W36-W67 N19-N38
      6  335 1376 (N17E55 L=347)   424     87        1 ( 306/1354- 368/1396)   0.0 E58-E51 N15-N18
      7 1255 1322 (N22W13 L= 55)   320    114        1 (1226/1271-1300/1357)  56.3 W10-W15 N18-N24
      8 1768 1205 (N19W54 L= 96)  1209    287        1 (1715/1070-1805/1308) -66.8 W49-W57 N10-N24
      9 1687 1139 (N14W45 L= 87)   442    108        0 (1656/1114-1735/1180) -23.2 W42-W49 N12-N16
     10  792  711 (S20E12 L= 30)   870    308        0 ( 750/ 609- 862/ 788) -59.0 E15-E07 S26-S14
     11 1725  639 (S17W57 L= 99)   963    132        1 (1676/ 596-1779/ 694)  40.6 W52-W61 S20-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