Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1368 1678 (N47W35 L=340)   347     75        0 (1336/1657-1401/1700)  12.5 W31-W40 N45-N48
      2 1484 1626 (N42W44 L=349)   443     62        2 (1450/1607-1519/1642)  -8.1 W39-W48 N40-N42
      3  662 1415 (N29E25 L=280)   224     45        1 ( 644/1387- 679/1437) -71.6 E26-E23 N27-N30
      4 1482 1225 (N15W31 L=336)   403    140        0 (1419/1184-1539/1263) -32.0 W26-W34 N12-N17
      5 1134 1140 (N11W07 L=312)   401     80        0 (1112/1105-1161/1189) -80.5 W05-W08 N08-N14
      6 1804  811 (S13W58 L=  3)   896    145        1 (1776/ 755-1837/ 868) -41.2 W54-W64 S16-S09
      7 1744  781 (S14W52 L=357)   387    112        0 (1723/ 722-1754/ 829)  81.9 W50-W52 S18-S11
      8  393  616 (S22E48 L=257)   444    120        0 ( 365/ 561- 405/ 664)  81.9 E50-E46 S25-S18
      9 1528  483 (S34W38 L=344)   582     91        1 (1491/ 463-1572/ 513) -18.4 W34-W43 S35-S31
  

    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