Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1258 1853 (N65W44 L=349)  1040    121        1 (1194/1833-1310/1879)   2.0 W32-W62 N62-N67
      2 1547 1575 (N37W48 L=353)   405     58        0 (1518/1559-1579/1590)   7.1 W44-W52 N36-N38
      3  952 1381 (N27E04 L=301)  1037    394        2 ( 891/1290-1059/1472) -59.7 E07-W03 N20-N33
      4 1551 1383 (N24W40 L=345)   338     72        0 (1526/1348-1573/1417) -78.7 W38-W41 N21-N26
      5  230 1117 (N10E60 L=246)   389     42        0 ( 217/1091- 247/1143) -53.1 E61-E57 N08-N12
      6 1747  855 (S10W51 L=356)   584     54        1 (1723/ 827-1775/ 880)  36.9 W48-W54 S11-S08
      7  793  814 (S08E15 L=290)   721    192        1 ( 717/ 776- 854/ 851)  -6.3 E20-E11 S10-S05
      8  333  650 (S19E54 L=251)   639     79        2 ( 308/ 616- 371/ 674)  18.4 E57-E48 S21-S17
      9 1546  595 (S26W37 L=342)   277     99        0 (1517/ 574-1585/ 632) -29.7 W33-W40 S27-S23
     10 1174  149 (S68W17 L=322)   650     46        0 (1159/ 120-1186/ 179) -42.9 W13-W27 S75-S62
  

    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