Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1386 1635 (N46W24 L=150)   297     78        0 (1364/1612-1414/1654)  31.0 W21-W26 N44-N48
      2  544 1485 (N25E39 L= 86)  2257    402        9 ( 417/1353- 634/1580) -24.6 E53-E29 N17-N30
      3 1231 1505 (N34W10 L=136)   424    128        1 (1210/1461-1263/1562) -60.3 W08-W12 N31-N38
      4 1680 1335 (N27W45 L=171)   347     98        1 (1653/1315-1724/1359) -18.4 W42-W48 N25-N28
      5 1724 1193 (N18W48 L=174)   965    152        3 (1686/1137-1757/1254) -60.3 W44-W51 N14-N21
      6  943 1125 (N07E06 L=120)   194     65        1 ( 908/1103- 969/1155)  26.6 E08-E04 N05-N08
      7  430 1094 (S00E40 L= 86)  1811    425        6 ( 327/1013- 568/1153)  15.5 E47-E29 S05-N03
      8  291  887 (S14E51 L= 75)   375    117        2 ( 273/ 842- 310/ 938) -90.0 E52-E48 S17-S11
      9  587  715 (S22E26 L=100)   262     68        1 ( 556/ 701- 619/ 733)   0.0 E28-E23 S23-S20
     10  389  663 (S28E43 L= 83)   702    238        2 ( 344/ 591- 446/ 722)  45.0 E48-E37 S33-S23
     11  729  555 (S32E16 L=110)   296     62        1 ( 697/ 536- 757/ 582)  21.8 E18-E13 S33-S29
     12  605  393 (S46E29 L= 97)   346     49        1 ( 582/ 373- 623/ 419)  38.7 E31-E26 S48-S43
  

    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