Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1775 1454 (N08W64 L=312)  2525    415       10 (1702/1362-1841/1540)  65.6 W57-W72 N01-N12
      2  123 1168 (N31E64 L=184)   884    154        1 (  95/1127- 152/1192) -45.0 E66-E60 N28-N32
      3  703 1014 (N10E19 L=229)   321     81        2 ( 654/ 983- 737/1038) -11.3 E21-E16 N07-N11
      4  871  904 (N00E11 L=236)   699    161        4 ( 786/ 877- 952/ 937) -24.0 E15-E06 S02-N03
      5  202  619 (S03E73 L=174)   991    186        2 ( 163/ 591- 233/ 646) -20.0 E80-E69 S04-N00
      6  561  482 (S17E43 L=204)   343    119        0 ( 519/ 460- 596/ 513)  -9.5 E46-E40 S18-S15
      7  321  458 (S14E67 L=180)  2641    276        4 ( 275/ 386- 360/ 512) -69.4 E73-E62 S19-S10
      8  731  477 (S21E32 L=215)   242     42        0 ( 703/ 466- 763/ 490) -26.6 E34-E30 S22-S20
      9  547  434 (S20E47 L=201)   401     63        1 ( 518/ 406- 572/ 453) -71.6 E48-E44 S21-S18
  

    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