Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  620 1828 (N32E66 L=265)   967     99        1 ( 592/1783- 649/1873) -39.3 E75-E60 N28-N37
      2  936 1711 (N32E29 L=301)  1232    206        1 ( 851/1657-1017/1753)  45.0 E32-E23 N26-N35
      3  784 1601 (N21E33 L=298)   362     74        1 ( 756/1564- 813/1637) -80.5 E34-E31 N17-N23
      4 1515 1555 (N42W16 L=347)   447     57        0 (1481/1533-1543/1578)  45.0 W13-W18 N39-N44
      5  577 1347 (N00E35 L=295)   393    106        0 ( 522/1334- 631/1361)   7.1 E39-E31 N00-N02
      6  303 1108 (S19E50 L=280)   358     68        0 ( 273/1085- 325/1146) -31.0 E53-E47 S20-S17
      7  637 1083 (S13E24 L=307)   592    206        0 ( 561/1030- 731/1116)   9.5 E28-E16 S15-S11
      8 1742 1009 (N15W45 L= 16)  1110    280        1 (1682/ 903-1813/1112) -28.1 W38-W54 N11-N19
      9 1274 1047 (N02W13 L=343)   583    147        1 (1234/1000-1326/1106)  76.0 W11-W14 S01-N07
     10 1535  997 (N08W30 L=  0)   269     47        0 (1520/ 969-1550/1032) -56.3 W28-W31 N06-N09
     11  266  964 (S29E52 L=279)   339     29        0 ( 251/ 937- 278/ 991) -71.6 E53-E50 S30-S26
     12  264  853 (S36E52 L=279)   628     81        1 ( 234/ 824- 315/ 882)  16.7 E55-E45 S38-S34
     13 1155  790 (S15W14 L=344)   271    103        0 (1114/ 767-1204/ 827)  63.4 W12-W15 S17-S11
     14  631  712 (S35E16 L=315)   692    202        0 ( 552/ 695- 742/ 732)  14.9 E22-E07 S37-S31
     15  425  668 (S44E33 L=297)   701    101        0 ( 381/ 645- 481/ 695)  10.3 E38-E27 S45-S41
     16 1409  442 (S26W45 L= 16)  2475    449        2 (1219/ 393-1614/ 509)  38.0 W34-W57 S34-S15
  

    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