Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  974 1898 (N52E46 L=158)  1150    109        5 ( 938/1863-1009/1927) -13.2 E54-E37 N47-N55
      2 1102 1796 (N47E23 L=180)   451     97        2 (1060/1782-1152/1816)   6.3 E27-E18 N45-N48
      3  827 1526 (N19E26 L=178)   267     94        1 ( 793/1488- 862/1550)  -9.5 E28-E22 N17-N20
      4 1531 1518 (N40W21 L=224)   621    122        4 (1475/1499-1597/1540)  14.0 W15-W27 N37-N43
      5  313 1467 (N04E62 L=141)   950    110        2 ( 290/1403- 337/1521) -54.5 E66-E58 N00-N07
      6  846 1272 (N05E16 L=187)   280     97        3 ( 803/1246- 881/1302)  -9.5 E19-E13 N04-N06
      7  156 1128 (S18E69 L=135)  1518    172        4 ( 128/1075- 196/1182) -80.5 E74-E61 S21-S15
      8  252  931 (S28E53 L=151)  1978    329        3 ( 178/ 773- 290/1056)  53.6 E66-E49 S39-S20
      9  639  823 (S26E18 L=185)   307    133        0 ( 592/ 797- 698/ 860)   0.0 E22-E13 S26-S24
     10 1210  561 (S25W25 L=229)   334    114        1 (1170/ 532-1261/ 597)  74.1 W23-W26 S28-S21
     11  413  530 (S52E37 L=167)   384     66        1 ( 399/ 501- 430/ 565)  59.0 E40-E34 S54-S49
     12  492  498 (S52E27 L=177)   302     50        1 ( 463/ 482- 519/ 516)   8.1 E30-E23 S53-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