Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  964 1916 (N53E40 L=355)   530     62        0 ( 928/1906-1000/1927)  15.9 E44-E36 N51-N55
      2  864 1740 (N34E30 L=  5)  1526    207        4 ( 786/1684- 954/1807) -18.4 E39-E21 N30-N38
      3 1408 1734 (N50W10 L= 45)   562     54        2 (1379/1709-1433/1763)  63.4 W07-W13 N46-N53
      4  533 1663 (N22E50 L=344)  1455    253        4 ( 419/1622- 618/1733)  39.8 E61-E45 N18-N28
      5  567 1592 (N18E44 L=351)   706    267        0 ( 506/1505- 657/1671)  77.0 E47-E40 N11-N24
      6  892 1555 (N22E20 L= 15)   221     97        0 ( 853/1534- 939/1569)  21.8 E22-E17 N19-N22
      7 1694 1391 (N32W38 L= 73)  1457    185        1 (1626/1317-1763/1464) -21.0 W31-W45 N28-N35
      8 1300 1353 (N18W09 L= 44)   352    152        0 (1241/1303-1342/1415) -32.0 W04-W12 N15-N20
      9  835 1131 (S05E13 L= 22)   679    175        2 ( 767/1071- 897/1200)  60.9 E15-E10 S09-N00
     10 1897 1053 (N17W62 L= 97)  2659    398        4 (1851/ 851-1938/1192) -54.2 W55-W71 N06-N24
     11  669  887 (S22E19 L= 16)   305     82        0 ( 657/ 838- 685/ 929)  90.0 E19-E17 S25-S19
     12  241  824 (S33E56 L=339)  1019    228        0 ( 182/ 736- 322/ 897)  27.6 E67-E46 S39-S27
     13  483  846 (S28E32 L=  3)   258     91        1 ( 465/ 796- 497/ 890)  80.5 E33-E30 S31-S25
     14  657  592 (S41E16 L= 19)   390     83        1 ( 610/ 579- 698/ 611)  23.2 E19-E12 S42-S38
     15 1516  465 (S26W50 L= 85)   727    135        0 (1463/ 422-1573/ 515)  -4.1 W43-W57 S27-S24
     16 1533  293 (S35W69 L=104)   589     77        0 (1507/ 258-1558/ 324)  -3.2 W62-W80 S35-S34
  

    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