Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  619 1758 (N28E54 L=353)   585     75        1 ( 590/1719- 641/1791) -21.8 E59-E49 N25-N30
      2  457 1721 (N22E67 L=340)  2433    159        5 ( 393/1671- 518/1772)  66.0 E73-E62 N17-N26
      3  579 1605 (N17E46 L=  1)  1339    271        4 ( 512/1516- 649/1694) -21.3 E55-E37 N11-N22
      4 1667 1415 (N37W33 L= 80)  2699    469        3 (1554/1206-1783/1577) -35.3 W21-W46 N27-N46
      5  376 1434 (N02E53 L=354)   973    216        2 ( 304/1410- 481/1485)  30.3 E60-E45 N00-N07
      6  901 1307 (N07E14 L= 33)  1003    303        2 ( 783/1250-1033/1367)  47.5 E19-E08 N01-N13
      7 1575 1263 (N24W27 L= 74)   613    171        0 (1507/1209-1638/1322) -14.0 W20-W32 N22-N25
      8  243 1234 (S12E59 L=348)   443     83        0 ( 208/1208- 273/1270)  45.0 E62-E56 S13-S08
      9 1067 1009 (S05W03 L= 50)   298    125        1 (1029/ 964-1122/1057)  66.8 W01-W04 S07-N00
     10  262  997 (S25E52 L=355)  1015    198        3 ( 200/ 942- 337/1056) -11.3 E59-E43 S26-S22
     11 1836  657 (S01W70 L=117)   831    106        2 (1811/ 614-1857/ 710) -20.6 W65-W77 S02-N02
     12 1619  442 (S19W61 L=108)   560    125        3 (1581/ 404-1664/ 481)  63.4 W57-W64 S22-S16
  

    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