Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  418 1629 (N23E57 L=338)  1371    263        7 ( 368/1532- 475/1728) -71.6 E64-E52 N16-N31
      2  491 1616 (N23E49 L=347)   524    108        0 ( 457/1579- 522/1661) -80.5 E50-E47 N20-N26
      3  516 1558 (N19E44 L=352)   618     99        1 ( 459/1530- 552/1583)   0.0 E48-E40 N17-N21
      4 1087 1398 (N17E02 L= 33)   201     59        2 (1058/1378-1110/1423)   0.0 E03-W00 N15-N18
      5 1842 1152 (N17W56 L= 91)   661    121        1 (1822/1093-1861/1206) -74.1 W53-W58 N13-N20
      6  840  828 (S21E08 L= 27)   477    180        0 ( 762/ 799- 929/ 866)  26.6 E13-E03 S22-S17
      7  284  625 (S40E58 L=338)   772     84        0 ( 271/ 576- 299/ 679)  41.2 E62-E53 S43-S36
      8 1383  500 (S31W36 L= 71)   397    115        0 (1354/ 446-1422/ 550)  90.0 W34-W37 S34-S27
      9  558  362 (S58E34 L=  2)   676     58        2 ( 527/ 337- 585/ 387) -26.6 E39-E28 S60-S55
  

    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