Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  982 1446 (N25W02 L= 51)   691    173        1 ( 885/1412-1053/1471) -21.8 E04-W06 N22-N27
      2  306 1325 (N25E49 L=360)  2178    327        2 ( 192/1231- 431/1436) -25.2 E71-E37 N17-N34
      3 1712 1047 (S06W46 L= 95)   274     55        0 (1700/1014-1724/1080)  63.4 W44-W47 S08-S04
      4 1861  954 (S13W61 L=110)   888    128        2 (1826/ 899-1901/1000)  40.6 W56-W66 S16-S10
      5 1384  761 (S20W20 L= 69)   370     77        1 (1335/ 744-1425/ 775)  -8.1 W16-W23 S21-S19
      6 1272  732 (S21W13 L= 62)   214     70        1 (1259/ 694-1292/ 767)  78.7 W11-W13 S23-S18
      7 1804  663 (S31W60 L=109)  1828    162        1 (1753/ 629-1875/ 714)   0.0 W53-W72 S33-S27
      8  467  497 (S27E49 L=  1)  1508    213        2 ( 419/ 399- 520/ 581)  30.1 E59-E40 S33-S22
  

    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