Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1163 1675 (N38E02 L=244)   951    204        1 (1061/1623-1255/1720) -10.6 E11-W05 N35-N39
      2  994 1542 (N25E11 L=236)   429    150        0 ( 930/1491-1055/1574) -16.7 E15-E05 N22-N26
      3 1688 1466 (N34W42 L=289)   489     94        0 (1661/1426-1717/1511) -51.3 W39-W44 N31-N37
      4  338  907 (S22E46 L=201)   782    200        0 ( 291/ 819- 368/ 993)  65.6 E50-E42 S27-S16
      5 1069  849 (S16W05 L=252)   183     69        0 (1026/ 831-1100/ 869)  36.9 W03-W07 S17-S14
      6 1344  811 (S13W23 L=270)   671    226        1 (1260/ 749-1440/ 863)  42.0 W18-W28 S17-S08
      7  542  803 (S26E29 L=218)   212     73        0 ( 511/ 784- 565/ 837)  45.0 E30-E26 S27-S23
      8  916  667 (S29E01 L=246)  2149    658        1 ( 628/ 588-1197/ 728)  21.3 E21-W16 S36-S20
      9 1708  653 (S15W55 L=302)   393     40        1 (1687/ 631-1728/ 680)  90.0 W53-W57 S16-S12
     10 1637  561 (S22W53 L=300)   727    139        0 (1570/ 541-1699/ 588)  26.6 W47-W58 S23-S18
     11 1046  205 (S61W30 L=277)  1091     99        2 ( 996/ 158-1082/ 255) -19.3 W20-W40 S65-S57
  

    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