Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1479 1740 (N32W63 L=137)  1173    101        2 (1429/1703-1526/1774) -63.4 W59-W68 N29-N35
      2  468 1461 (N48E31 L= 43)   401     73        1 ( 429/1435- 497/1486) -20.6 E35-E27 N45-N49
      3 1832 1156 (S10W58 L=132)   882    131        0 (1791/1108-1879/1211) -48.8 W54-W62 S13-S05
      4 1768  984 (S17W47 L=121)  1366    265        0 (1744/ 872-1788/1112)  66.8 W43-W49 S23-S09
      5  623 1035 (N17E24 L= 51)   270     76        2 ( 583/1011- 651/1065) -45.0 E25-E21 N15-N19
      6  374  834 (N11E46 L= 28)   760    184        1 ( 345/ 742- 408/ 901)  84.3 E48-E43 N05-N15
      7  166  761 (N09E73 L=  1)  2484    150        3 ( 120/ 699- 230/ 827) -31.6 E83-E63 N05-N13
      8 1564  429 (S48W22 L= 96)   348     54        0 (1542/ 408-1589/ 452)  11.3 W19-W25 S49-S46
      9  924  292 (S35E31 L= 43)   691    154        0 ( 871/ 244- 966/ 365)  23.2 E38-E24 S37-S30
  

    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