Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  638 1562 (N33E35 L=179)  1298    182        3 ( 569/1496- 701/1609) -20.6 E43-E27 N28-N36
      2  342 1522 (N27E63 L=151)  1565    129        3 ( 292/1462- 391/1572) -15.9 E78-E55 N23-N29
      3  149  906 (S13E71 L=143)   619     48        1 ( 132/ 878- 167/ 932)  29.7 E75-E66 S15-S11
      4  563  833 (S14E29 L=186)   499    121        1 ( 535/ 786- 593/ 895)  66.8 E30-E26 S16-S09
      5  279  783 (S20E54 L=160)   736    178        1 ( 246/ 749- 318/ 831)  32.0 E58-E49 S21-S16
      6  946  771 (S14E03 L=211)   568    203        0 ( 849/ 749-1035/ 795)   5.2 E09-W02 S15-S13
      7 1536  658 (S17W38 L=253)   497     99        0 (1493/ 636-1594/ 679)  -6.3 W34-W43 S18-S16
      8 1633  488 (S28W55 L=269)   411     48        1 (1607/ 463-1652/ 512) -63.4 W51-W57 S30-S26
      9 1533  437 (S33W47 L=262)   393     49        1 (1502/ 424-1565/ 458)   0.0 W44-W51 S33-S31
     10  603  400 (S44E31 L=183)   354     39        2 ( 574/ 379- 627/ 418)  11.3 E33-E28 S46-S42
  

    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