Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1144 1644 (N44W06 L=263)   641    105        1 (1096/1617-1192/1668)  12.5 W01-W10 N42-N46
      2 1415 1605 (N43W30 L=287)  1725    293        2 (1268/1570-1534/1631)   0.0 W16-W42 N40-N44
      3  407 1466 (N26E50 L=207)   776    140        1 ( 354/1413- 455/1531)  63.4 E53-E46 N22-N30
      4 1869 1299 (N22W71 L=328)  1874    150        1 (1824/1237-1901/1358) -45.0 W65-W75 N18-N25
      5 1338 1142 (N11W19 L=276)  1876    410        4 (1279/ 957-1379/1324) -77.2 W15-W21 N00-N22
      6  608 1104 (N05E27 L=231)  1355    393        1 ( 557/ 956- 656/1277) -75.3 E32-E23 S04-N15
      7 1822 1255 (N19W61 L=318)   494     19        0 (1800/1234-1843/1276) -11.3 W58-W64 N18-N20
      8  816 1188 (N11E14 L=243)   736    289        1 ( 759/1096- 921/1264)  57.5 E17-E07 N05-N16
      9  259  909 (S09E54 L=203)   821    152        0 ( 235/ 845- 283/ 993) -66.0 E56-E52 S13-S04
     10  869  766 (S14E09 L=249)  1266    395        1 ( 787/ 619- 919/ 936)  71.6 E13-E05 S24-S03
     11 1497  636 (S20W35 L=292)   247     69        0 (1471/ 620-1535/ 659) -21.8 W32-W37 S20-S18
     12  804  539 (S30E13 L=244)   328     81        1 ( 785/ 501- 831/ 578) -68.2 E14-E10 S32-S27
     13 1721  460 (S32W68 L=325)   696     47        0 (1696/ 439-1748/ 480) -10.6 W61-W77 S33-S30
     14 1565  373 (S39W54 L=311)   738     84        1 (1525/ 346-1599/ 405)  33.7 W49-W57 S41-S36
     15 1223  297 (S47W23 L=281)   348     42        1 (1198/ 276-1251/ 322)  26.6 W21-W25 S48-S44
     16  695  263 (S55E29 L=228)  2023    225        4 ( 614/ 227- 805/ 292)   6.6 E43-E16 S59-S50
  

    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