Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  661 1701 (N61E07 L=245)   423     82        0 ( 640/1663- 691/1734)  35.5 E09-E01 N57-N63
      2 1275 1644 (N35W39 L=290)  1396    259        2 (1165/1583-1382/1686) -47.7 W32-W43 N28-N40
      3  447 1619 (N62E39 L=212)   385     63        1 ( 418/1600- 468/1643)  23.2 E43-E35 N59-N63
      4  320 1154 (N31E47 L=205)   811    159        2 ( 246/1127- 401/1184) -19.4 E55-E38 N27-N33
      5 1775  858 (S25W46 L=297)  3048    487        4 (1653/ 705-1890/1016)  18.4 W32-W62 S31-S18
      6  898  430 (S25E25 L=227)   636    217        1 ( 842/ 398- 953/ 475) -15.9 E28-E20 S27-S20
      7 1552  430 (S47W22 L=273)   521    115        0 (1521/ 385-1596/ 467)  32.0 W18-W26 S50-S44
      8  571  436 (S17E47 L=205)   373     97        0 ( 516/ 418- 610/ 454) -26.6 E50-E43 S18-S14
      9  530  348 (S22E57 L=194)  3063    243        1 ( 443/ 296- 636/ 416) -45.0 E64-E49 S26-S16
     10  633  386 (S22E45 L=206)   308     65        1 ( 602/ 370- 662/ 402)   0.0 E48-E42 S23-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