Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1855 (N65E07 L= 80)  4856    429        4 (1003/1768-1348/1926)  -1.1 E41-W17 N58-N71
      2 1083 1694 (N45E09 L= 78)  2254    408        7 ( 882/1627-1231/1746)  -3.6 E28-W04 N42-N47
      3  561 1664 (N33E51 L= 36)  1290    261        2 ( 488/1581- 672/1705)  31.6 E55-E41 N26-N37
      4  353 1402 (N13E55 L= 32)  1765    341        5 ( 301/1270- 405/1541) -46.8 E63-E47 N05-N21
      5  765 1502 (N26E26 L= 61)   391    102        0 ( 731/1482- 809/1527)  36.9 E28-E22 N24-N27
      6 1539 1197 (N18W31 L=118)  1088    321        1 (1447/1068-1624/1301) -38.2 W23-W38 N11-N23
      7 1507 1147 (N15W29 L=116)   223     81        0 (1468/1123-1547/1167)  33.7 W26-W32 N12-N16
      8 1874 1050 (N14W65 L=152)   664     91        0 (1859/1014-1895/1094) -53.1 W62-W68 N12-N17
      9  387  986 (S12E42 L= 45)   364    128        0 ( 341/ 938- 437/1026)  45.0 E45-E38 S15-S08
     10  172  876 (S22E67 L= 20)  1288    107        1 ( 147/ 832- 190/ 936)  54.5 E72-E63 S24-S17
     11  282  773 (S27E53 L= 34)   402    100        0 ( 259/ 730- 307/ 817)  45.0 E56-E49 S30-S23
     12  768  748 (S21E12 L= 75)   314     62        1 ( 734/ 717- 799/ 769) -31.0 E14-E09 S22-S19
     13 1169  658 (S20W15 L=102)   759    303        0 (1048/ 557-1275/ 746)  53.7 W09-W20 S28-S13
     14  600  705 (S26E24 L= 63)   695    165        1 ( 523/ 685- 672/ 736)  20.0 E29-E18 S27-S23
     15  694  209 (S69E23 L= 64)   932     84        1 ( 662/ 178- 728/ 252) -78.7 E29-E17 S73-S63
     16  951  177 (S64W19 L=106)   655     72        0 ( 921/ 153- 987/ 210) -21.0 W12-W26 S67-S60
  

    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