Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  720 1869 (N38E60 L=135)   613     47        1 ( 691/1849- 752/1886)   0.0 E66-E55 N36-N39
      2  653 1833 (N33E60 L=135)   747     83        1 ( 626/1796- 669/1879) -26.6 E69-E54 N30-N37
      3  634 1700 (N23E48 L=147)   332     68        0 ( 603/1678- 667/1724)  63.4 E49-E46 N21-N25
      4  800 1619 (N22E32 L=163)   578    244        1 ( 756/1504- 851/1715) -34.5 E39-E23 N16-N27
      5  369 1504 (N05E58 L=137)  3135    400       13 ( 294/1328- 473/1619)  87.1 E64-E52 S06-N14
      6 1089 1545 (N26E12 L=183)   257    101        0 (1073/1492-1106/1592) -56.3 E13-E08 N22-N28
      7 1183 1316 (N14W01 L=196)   322    130        1 (1154/1254-1199/1364) -56.3 E02-W02 N10-N16
      8  392 1074 (S19E41 L=154)  1528    311        3 ( 279/ 956- 490/1177)  54.8 E49-E35 S27-S10
      9  867 1113 (S05E11 L=184)   215    105        0 ( 834/1063- 883/1149) -45.0 E13-E09 S07-S03
     10  492  848 (S30E29 L=166)  2694    663        4 ( 208/ 731- 791/ 938)  27.0 E61-E10 S43-S17
     11 1399  809 (S07W27 L=222)   227     89        0 (1370/ 768-1437/ 838)  68.2 W26-W28 S09-S04
     12 1133  629 (S24W18 L=213)   354    146        1 (1066/ 593-1191/ 661)  54.5 W15-W21 S28-S21
     13  550  481 (S53E19 L=176)   348     91        1 ( 511/ 444- 584/ 516) -18.4 E23-E14 S54-S50
     14 1624  441 (S18W62 L=257)   740    124        4 (1583/ 395-1666/ 473)  81.9 W57-W64 S22-S15
     15 1664  404 (S19W70 L=265)   511     84        1 (1637/ 380-1685/ 442) -76.0 W66-W73 S20-S16
     16 1019  296 (S48W30 L=225)  2010    431        4 ( 838/ 203-1164/ 407)  -1.4 W09-W50 S50-S45
  

    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