Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  949 1702 (N51W19 L=291)   457     74        1 ( 924/1654- 971/1737)  40.6 W15-W22 N47-N53
      2 1277 1656 (N37W39 L=312)   266     45        0 (1263/1628-1298/1679)  31.0 W36-W41 N35-N38
      3 1362 1656 (N34W45 L=318)   300     43        1 (1335/1640-1396/1669) -33.7 W42-W48 N33-N35
      4 1644 1568 (N20W62 L=335)  2012    175        5 (1574/1522-1727/1608) -25.0 W57-W72 N16-N24
      5 1490 1556 (N25W48 L=320)   415     96        0 (1469/1512-1506/1594)  63.4 W44-W50 N21-N27
      6  409 1324 (N39E38 L=234)  1211    171        0 ( 323/1276- 479/1357) -24.0 E49-E31 N34-N42
      7  877 1253 (N24E03 L=269)  1088    420        3 ( 703/1221-1070/1279) -23.2 E14-W07 N17-N27
      8 1203 1131 (N09W13 L=285)   406    118        0 (1154/1099-1243/1170)  18.4 W09-W15 N07-N10
      9 1648 1102 (S05W40 L=313)   951    155        1 (1615/1021-1675/1170)  41.2 W36-W44 S08-S01
     10  393 1051 (N22E41 L=231)  1063    184        1 ( 327/ 986- 461/1108) -42.0 E46-E36 N17-N26
     11 1545  575 (S34W24 L=296)  5360    391        8 (1329/ 331-1754/ 904)  34.3 W02-W45 S48-S18
     12 1091  832 (S05E01 L=272)   203     57        1 (1072/ 805-1113/ 864) -90.0 E01-W00 S07-S03
     13 1097  604 (S19E06 L=266)   209     58        0 (1083/ 577-1114/ 638)  78.7 E06-E05 S21-S16
     14 1074  530 (S23E09 L=263)   268     95        0 (1029/ 506-1109/ 559)   9.5 E12-E06 S23-S22
     15  500  397 (S20E55 L=217)   390     46        0 ( 486/ 361- 512/ 423)  51.3 E57-E52 S22-S17
     16  927  319 (S35E28 L=245)   846    174        1 ( 861/ 260- 988/ 382)  25.0 E35-E20 S38-S31
     17 1161  232 (S50E17 L=255)   730     84        3 (1113/ 205-1210/ 255) -12.5 E21-E11 S51-S47
  

    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