Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  822 1873 (N56E51 L=203)   678     71        1 ( 782/1857- 860/1885)   4.8 E58-E45 N54-N57
      2 1152 1735 (N50E05 L=249)   374     63        0 (1132/1704-1182/1755)  78.7 E07-E02 N47-N52
      3  474 1672 (N31E61 L=193)   637     69        3 ( 447/1646- 503/1697) -10.3 E66-E55 N29-N33
      4 1303 1635 (N45W10 L=265)   428     34        0 (1275/1614-1331/1656) -26.6 W07-W12 N43-N46
      5  737 1559 (N28E30 L=224)   255     53        0 ( 706/1545- 767/1576)  33.7 E32-E28 N27-N29
      6 1792 1357 (N33W58 L=312)   536     59        1 (1756/1341-1817/1380)  16.7 W52-W62 N31-N35
      7 1580 1250 (N22W34 L=288)   616    219        0 (1501/1185-1648/1317) -26.6 W27-W39 N19-N25
      8  276  948 (S16E52 L=202)   532    126        0 ( 243/ 919- 319/1007) -20.6 E55-E47 S18-S13
      9  269  875 (S21E53 L=201)   563     78        2 ( 238/ 850- 291/ 907)  33.7 E56-E50 S22-S18
     10 1457  440 (S30W44 L=299)  1084    177        1 (1374/ 418-1541/ 471)  18.4 W38-W51 S32-S26
     11 1046  196 (S59W29 L=283)   671     69        2 (1024/ 160-1067/ 233) -36.9 W23-W35 S63-S54
  

    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