Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  827 1871 (N56E50 L=203)   631     63        0 ( 790/1858- 863/1881)   5.2 E56-E44 N53-N57
      2 1162 1732 (N50E04 L=249)   389     61        1 (1140/1699-1191/1754)  45.0 E06-E01 N47-N52
      3  477 1676 (N32E61 L=192)   848     92        2 ( 449/1644- 507/1716) -14.0 E71-E55 N29-N34
      4 1314 1633 (N45W11 L=265)   431     54        1 (1283/1611-1346/1655)   0.0 W08-W14 N43-N46
      5  770 1559 (N29E28 L=226)   419    154        0 ( 703/1535- 836/1579)  26.6 E31-E22 N26-N31
      6 1800 1357 (N34W59 L=313)   611     48        1 (1765/1340-1828/1381)  10.3 W53-W64 N31-N35
      7 1592 1250 (N23W35 L=288)   678    213        1 (1511/1184-1655/1317) -26.6 W28-W40 N19-N25
      8  298  935 (S17E50 L=204)   314     59        1 ( 269/ 920- 329/ 950)   0.0 E52-E46 S17-S16
      9  277  868 (S21E52 L=201)   442     59        3 ( 256/ 844- 296/ 901) -71.6 E54-E49 S22-S19
     10 1471  437 (S29W46 L=299)  1124    190        0 (1388/ 411-1560/ 459)  21.0 W39-W53 S32-S26
     11 1054  195 (S59W29 L=283)   680     72        0 (1033/ 164-1075/ 233) -36.0 W23-W34 S62-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