Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1607 1459 (N28W48 L=301)  1470    243        1 (1544/1391-1686/1557) -71.6 W44-W54 N22-N34
      2  352 1420 (N32E54 L=199)   713    119        0 ( 296/1393- 399/1451)  -7.6 E62-E47 N29-N33
      3  256 1373 (N28E64 L=188)   524     65        1 ( 239/1341- 284/1406) -63.4 E67-E60 N26-N30
      4  317 1324 (N25E54 L=199)   936    137        0 ( 261/1293- 383/1346) -11.3 E61-E46 N23-N26
      5  689 1249 (N20E21 L=232)   373    121        2 ( 629/1209- 739/1286)  35.5 E25-E17 N17-N22
      6 1665  882 (S09W43 L=296)   833    175        1 (1637/ 829-1697/ 950)  60.3 W40-W45 S11-S04
      7 1749  859 (S11W51 L=304)   499     47        0 (1723/ 837-1778/ 880) -11.3 W48-W54 S12-S09
      8  585  549 (S24E34 L=219)   763    192        0 ( 559/ 458- 623/ 628)  80.5 E35-E31 S30-S18
      9  392  554 (S24E52 L=201)   535     94        2 ( 357/ 528- 432/ 596)  20.0 E56-E45 S25-S20
  

    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