Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1832 (N54W08 L=267)   350     31        1 (1189/1808-1232/1850) -31.0 W05-W10 N51-N56
      2  426 1620 (N29E53 L=205)   503     88        2 ( 412/1567- 441/1665) -45.0 E57-E50 N25-N32
      3 1633 1555 (N35W43 L=302)   438     39        0 (1620/1525-1648/1585)  76.0 W41-W45 N33-N37
      4 1312 1509 (N28W15 L=274)   487    138        0 (1267/1460-1374/1536)  35.5 W12-W19 N24-N30
      5 1563 1210 (N11W33 L=291)   308     54        0 (1528/1187-1588/1229) -14.0 W30-W34 N09-N12
      6 1848 1072 (N07W58 L=317)   490     64        1 (1832/1034-1865/1116)  78.7 W56-W60 N04-N09
      7  197  848 (S19E62 L=197)   795    115        3 ( 149/ 802- 241/ 887) -26.6 E68-E57 S22-S17
      8 1512  740 (S17W35 L=293)   290     61        0 (1490/ 708-1533/ 773)  63.4 W33-W35 S19-S15
      9 1702  650 (S20W53 L=312)   480     57        2 (1683/ 617-1726/ 678) -45.0 W50-W55 S22-S18
     10 1446  388 (S41W43 L=302)   767     97        2 (1410/ 350-1471/ 444) -60.3 W40-W47 S43-S36
     11  528  312 (S57E49 L=209)  1081    184        1 ( 470/ 282- 618/ 336)   8.9 E66-E34 S59-S54
  

    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