Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1165 1573 (N36W25 L=294)  1450    360        3 ( 993/1519-1308/1600) -11.8 W11-W35 N32-N39
      2 1494 1579 (N27W49 L=317)   814    170        0 (1413/1560-1575/1596) -22.6 W43-W56 N23-N29
      3 1614 1507 (N19W54 L=323)   513     92        1 (1578/1469-1650/1538) -78.7 W51-W56 N16-N21
      4 1350 1267 (N13W25 L=294)   715    269        0 (1235/1200-1454/1315) -45.0 W18-W30 N07-N18
      5  997 1189 (N17W02 L=270)   397     93        2 ( 965/1142-1033/1226)  59.0 W00-W04 N14-N19
      6 1841 1119 (S08W59 L=327)   635    109        3 (1819/1069-1874/1171) -69.4 W56-W62 S12-S04
      7  478 1061 (N20E34 L=234)   380    101        0 ( 442/1025- 526/1094) -50.2 E36-E31 N16-N22
      8 1865  860 (S25W59 L=327)  2761    424       11 (1777/ 728-1920/1030)  41.2 W47-W69 S31-S15
      9  483  876 (N09E37 L=232)   308    116        1 ( 455/ 819- 515/ 922) -66.8 E38-E35 N04-N11
     10  584  736 (S01E33 L=235)   993    317        1 ( 522/ 625- 705/ 806)  17.4 E40-E24 S06-N03
     11  469  599 (S08E46 L=223)  1529    390        6 ( 363/ 470- 561/ 687) -66.4 E52-E41 S17-N00
     12  410  565 (S09E52 L=216)   874    124        1 ( 372/ 521- 446/ 610) -90.0 E54-E49 S11-S05
     13  279  523 (S11E71 L=197)   988     76        4 ( 247/ 491- 309/ 553) -45.0 E76-E66 S12-S08
     14  478  475 (S16E51 L=217)  1608    209        5 ( 441/ 387- 521/ 544)  37.6 E60-E45 S21-S11
     15 1395  402 (S42W11 L=279)  1139    239        2 (1316/ 333-1480/ 483)  25.0 W03-W18 S46-S37
     16  851  424 (S26E26 L=242)   493    178        2 ( 783/ 397- 920/ 462)   0.0 E31-E21 S28-S23
     17  548  408 (S21E49 L=219)   481     64        2 ( 517/ 370- 569/ 438) -90.0 E51-E45 S24-S19
     18  819  268 (S37E38 L=231)   334     55        1 ( 791/ 251- 851/ 280)   0.0 E40-E34 S38-S35
  

    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