Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1261 1869 (N59E20 L=271)   709     40        1 (1239/1838-1279/1893) -26.6 E25-E15 N55-N63
      2 1091 1487 (N21E10 L=281)   286    108        0 (1054/1439-1114/1526) -26.6 E13-E06 N19-N22
      3  337 1412 (S02E56 L=235)   330     50        0 ( 309/1393- 362/1432)   0.0 E59-E53 S02-N00
      4 1800 1183 (N27W47 L=338)  2277    443        4 (1719/1034-1860/1365) -68.2 W38-W56 N19-N39
      5 1583 1164 (N18W29 L=320)   375    105        1 (1557/1117-1609/1220)  81.9 W27-W30 N14-N21
      6 1646 1006 (N11W37 L=328)   484     97        1 (1629/ 946-1663/1053) -56.3 W35-W39 N07-N13
      7  243  963 (S29E54 L=237)  1885    265        2 ( 125/ 941- 386/ 993)   6.7 E73-E39 S31-S26
      8 1829  748 (N04W63 L=354)   746    107        0 (1770/ 732-1879/ 768)  19.7 W56-W70 N01-N06
      9  861  717 (S28E00 L=291)   405     95        0 ( 821/ 681- 895/ 758)  81.9 E02-W00 S31-S24
     10 1300  546 (S24W32 L=323)   663    220        0 (1198/ 522-1405/ 569)  36.0 W26-W38 S27-S19
     11 1551  416 (S23W58 L=349)  2058    243        5 (1494/ 331-1619/ 504) -69.0 W54-W64 S28-S15
     12  781  305 (S59W13 L=304)   319     48        0 ( 753/ 289- 811/ 320)  31.0 W10-W15 S61-S57
  

    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