Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  462 1701 (N25E60 L=178)   717    144        4 ( 414/1666- 528/1732)  45.0 E64-E56 N22-N28
      2  705 1589 (N21E33 L=205)  1082    208        1 ( 685/1495- 741/1681) -69.0 E36-E29 N14-N27
      3 1349 1554 (N33W10 L=248)  1270    284        1 (1286/1443-1412/1663) -49.8 W03-W15 N25-N39
      4 1865 1101 (N18W58 L=296)   486     77        0 (1832/1067-1890/1119) -26.6 W53-W61 N15-N19
      5  196  946 (S25E61 L=177)  1932    252        2 ( 138/ 865- 283/1054) -39.8 E69-E50 S28-S18
      6  321  772 (S34E47 L=191)   376    103        1 ( 264/ 755- 365/ 797)   0.0 E53-E42 S34-S32
      7  343  489 (S54E57 L=181)  1540    123        1 ( 291/ 441- 403/ 539)   8.6 E78-E45 S57-S49
      8 1191  348 (S44W33 L=271)   652    133        1 (1157/ 299-1232/ 393) -15.3 W29-W40 S46-S39
      9 1528  342 (S34W61 L=299)  1379    209        2 (1455/ 316-1618/ 378)  36.0 W54-W69 S37-S28
  

    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