Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  729 1657 (N32E37 L=197)   634    126        2 ( 676/1628- 774/1680)  26.6 E40-E32 N29-N34
      2 1054 1555 (N32E10 L=224)   851    173        1 ( 997/1487-1118/1610) -28.6 E15-E04 N28-N35
      3  723 1493 (N21E30 L=204)   401    128        2 ( 660/1473- 791/1522)  26.6 E33-E25 N18-N24
      4  949 1271 (N12E09 L=224)   284    111        0 ( 903/1255-1001/1300)   0.0 E12-E05 N10-N12
      5 1610 1120 (N16W36 L=270)   630     86        2 (1587/1080-1638/1166)  78.7 W34-W38 N13-N18
      6  187 1060 (S14E64 L=170)   905    127        2 ( 156/1005- 223/1109) -80.5 E68-E58 S17-S11
      7 1627  942 (N06W40 L=274)   248     75        1 (1607/ 907-1659/ 957)  53.1 W38-W42 N03-N07
      8  885  788 (S18E04 L=230)   218     62        0 ( 863/ 767- 917/ 813)   0.0 E05-E01 S19-S16
      9 1011  761 (S17W04 L=238)   202     23        1 ( 988/ 743-1031/ 776)   0.0 W02-W05 S18-S16
     10 1426  713 (S12W32 L=265)   816    255        1 (1307/ 663-1523/ 743)  31.6 W25-W38 S17-S09
     11  805  667 (S27E07 L=227)   302    161        1 ( 732/ 642- 878/ 687)  33.7 E11-E02 S30-S24
     12 1461  592 (S18W38 L=272)   538    169        3 (1390/ 559-1534/ 627)  41.2 W33-W42 S21-S14
     13  688  332 (S56E12 L=222)   387     79        0 ( 647/ 316- 720/ 352)  20.6 E16-E08 S58-S53
  

    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