Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1827 (N52E36 L=164)   708    100        1 ( 821/1809- 928/1842)  32.0 E40-E30 N48-N54
      2  549 1541 (N24E43 L=157)  4317    515        4 ( 398/1377- 717/1695) -17.0 E74-E25 N16-N31
      3  308 1425 (N13E60 L=140)   804    120        3 ( 290/1373- 326/1488) -49.4 E64-E56 N09-N16
      4  551 1373 (N13E36 L=163)   446    140        1 ( 517/1310- 572/1433) -40.6 E40-E33 N09-N16
      5 1658 1075 (N12W41 L=240)   561     67        1 (1633/1038-1690/1105)  51.3 W38-W43 N09-N14
      6 1637 1006 (N08W40 L=239)  1086    205        1 (1594/ 933-1695/1100) -45.0 W35-W45 N04-N13
      7 1463 1019 (N06W27 L=227)   500    159        1 (1449/ 936-1479/1091) -71.6 W25-W28 N01-N10
      8  309 1055 (S09E49 L=151)   292     62        0 ( 284/1034- 340/1082) -18.4 E51-E45 S09-S07
      9  216  985 (S15E58 L=141)   372     74        1 ( 199/ 950- 241/1007) -71.6 E60-E55 S16-S13
     10 1761  974 (N08W51 L=251)   292     40        0 (1740/ 948-1782/ 997) -26.6 W48-W53 N06-N09
     11  391  915 (S16E40 L=159)   434    132        0 ( 353/ 865- 449/ 948)  36.9 E43-E35 S19-S13
     12  673  776 (S20E18 L=181)  1059    267        1 ( 534/ 746- 757/ 827)  28.1 E28-E13 S24-S16
     13 1641  674 (S12W48 L=247)  1364    250        1 (1572/ 566-1700/ 781)  78.7 W45-W50 S19-S04
     14  713  356 (S50E13 L=186)   563     66        1 ( 680/ 335- 750/ 386)   8.1 E16-E09 S52-S48
     15  893  271 (S54W06 L=205)   903    167        2 ( 841/ 217- 943/ 329) -36.9 W00-W13 S58-S49
     16 1156  176 (S57W41 L=241)   558     54        1 (1137/ 154-1179/ 210) -22.6 W35-W47 S59-S53
  

    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