Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1173 1751 (N54W03 L= 84)   474     70        2 (1134/1729-1213/1766)   0.0 W00-W07 N51-N55
      2  437 1508 (N25E51 L= 29)  1971    377        3 ( 370/1427- 493/1634) -71.6 E57-E47 N19-N34
      3 1357 1380 (N27W19 L= 99)  1397    373        1 (1233/1250-1429/1491) -45.0 W09-W24 N19-N33
      4 1776 1392 (N32W60 L=140)   493     86        0 (1745/1360-1795/1425) -63.4 W57-W62 N29-N33
      5 1129 1230 (N15W04 L= 85)   397    165        0 (1093/1180-1183/1295) -50.2 W01-W08 N12-N18
      6 1568 1209 (N18W35 L=115)   403    150        0 (1540/1178-1595/1250) -53.1 W32-W37 N15-N20
      7  432  729 (S23E39 L= 41)   306     58        0 ( 400/ 705- 459/ 755)  38.7 E42-E36 S25-S21
      8 1249  681 (S18W18 L= 98)   611    179        1 (1181/ 627-1323/ 732)  41.2 W14-W22 S21-S13
      9 1681  608 (S18W54 L=134)  1270    195        2 (1647/ 517-1718/ 680) -79.7 W50-W58 S24-S13
     10 1612  576 (S21W49 L=129)   387     64        1 (1595/ 551-1637/ 607)  90.0 W47-W51 S22-S18
     11  414  507 (S40E47 L= 33)  1065    134        3 ( 361/ 477- 479/ 542)  18.4 E55-E37 S42-S36
  

    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