Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  704 1683 (N50E29 L= 56)   428     62        1 ( 678/1654- 736/1704)  53.1 E31-E25 N47-N52
      2 1621 1610 (N39W60 L=144)  1168    136        3 (1547/1586-1678/1633) -10.6 W51-W69 N36-N41
      3 1091 1440 (N30W06 L= 91)  2564    721        2 ( 796/1265-1332/1552) -29.6 E16-W21 N18-N39
      4 1599 1298 (N18W42 L=126)   811    196        1 (1575/1225-1643/1378)  78.7 W39-W45 N13-N23
      5  280 1168 (N14E55 L= 30)   371    109        2 ( 266/1116- 298/1223) -81.9 E56-E52 N10-N17
      6 1744  715 (S19W53 L=137)  1245    227        2 (1684/ 632-1805/ 817)  55.0 W48-W58 S24-S13
      7 1255  704 (S17W14 L= 98)  1099    233        2 (1185/ 599-1321/ 793)  50.7 W09-W18 S23-S10
      8  556  744 (S12E32 L= 52)   444    107        1 ( 521/ 702- 596/ 791)  45.0 E35-E28 S14-S09
      9  361  692 (S16E49 L= 35)   627    157        1 ( 334/ 611- 390/ 755)  48.0 E54-E45 S21-S11
     10 1441  463 (S35W31 L=115)  2254    394        4 (1288/ 361-1583/ 583)  33.7 W19-W40 S42-S26
     11  482  504 (S29E44 L= 40)   396     59        1 ( 451/ 487- 510/ 529)   9.5 E47-E41 S30-S26
     12 1633  440 (S39W53 L=137)  1467    129        3 (1586/ 397-1693/ 487)  14.0 W45-W60 S42-S35
  

    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