Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1482 1742 (N32W64 L=137)  1183     98        1 (1430/1706-1527/1771) -56.3 W60-W68 N29-N35
      2  473 1464 (N48E30 L= 44)   442     76        1 ( 431/1436- 500/1492)   8.1 E34-E27 N45-N50
      3 1836 1158 (S10W58 L=132)   964    120        0 (1795/1109-1884/1212) -53.1 W54-W63 S13-S05
      4 1772  992 (S17W47 L=121)  1485    245        1 (1747/ 874-1792/1111)  70.3 W43-W50 S23-S09
      5  633 1034 (N17E23 L= 51)   271     52        1 ( 610/1013- 666/1057) -45.0 E24-E20 N15-N18
      6  369  867 (N13E46 L= 28)   481     56        2 ( 349/ 824- 390/ 902)  78.7 E47-E43 N10-N15
      7  166  765 (N10E73 L=  1)  3029    202        2 ( 115/ 699- 234/ 849) -32.7 E82-E63 N05-N14
      8 1563  425 (S48W22 L= 96)   345     54        0 (1543/ 399-1585/ 455)  45.0 W19-W24 S50-S46
      9  932  298 (S35E30 L= 44)   752    147        0 ( 876/ 248- 972/ 368)  26.6 E37-E23 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