Advertisement
Guest User

Untitled

a guest
Jul 17th, 2017
58
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 2.98 KB | None | 0 0
  1. f = xmlread('/home/ivan/Downloads/sentinel1_position.xml');
  2. orbits = f.getElementsByTagName('orbit');
  3.  
  4. Lo = orbits.getLength;
  5. d = zeros(1, Lo);
  6. d = datetime(d, d, d, d, d, d);
  7. x = zeros(1, Lo);
  8. y = zeros(1, Lo);
  9. z = zeros(1, Lo);
  10. vx = zeros(1, Lo);
  11. vy = zeros(1, Lo);
  12. vz = zeros(1, Lo);
  13.  
  14. for i = 1:orbits.getLength
  15.  
  16.   orbit = orbits.item(i-1);
  17.   date = orbit.getElementsByTagName('time').item(0).item(0).getData;
  18.   d(i) = datetime(char(date), 'InputFormat', 'yyyy-MM-dd''T''hh:mm:ss.SSSSSS');
  19.   position = orbit.getElementsByTagName('position').item(0);
  20.   x(i) = str2num(position.getElementsByTagName('x').item(0).item(0).getData);
  21.   y(i) = str2num(position.getElementsByTagName('y').item(0).item(0).getData);
  22.   z(i) = str2num(position.getElementsByTagName('z').item(0).item(0).getData);
  23.   velocity = orbit.getElementsByTagName('velocity').item(0);
  24.   vx(i) = str2num(velocity.getElementsByTagName('x').item(0).item(0).getData);
  25.   vy(i) = str2num(velocity.getElementsByTagName('y').item(0).item(0).getData);
  26.   vz(i) = str2num(velocity.getElementsByTagName('z').item(0).item(0).getData);
  27.  
  28. end
  29.  
  30. attitudes = f.getElementsByTagName('attitude');
  31.  
  32. La = attitudes.getLength;
  33. t = zeros(1, La);
  34. t = datetime(t, t, t, t, t, t);
  35. q0 = zeros(1, La);
  36. q1 = zeros(1, La);
  37. q2 = zeros(1, La);
  38. q3 = zeros(1, La);
  39. wx = zeros(1, La);
  40. wy = zeros(1, La);
  41. wz = zeros(1, La);
  42. roll = zeros(1, La);
  43. pitch = zeros(1, La);
  44. yaw = zeros(1, La);
  45.  
  46.  
  47. for i = 1:attitudes.getLength
  48.  
  49.   attitude = attitudes.item(i-1);
  50.   date = attitude.getElementsByTagName('time').item(0).item(0).getData;
  51.   t(i) = datetime(char(date), 'InputFormat', 'yyyy-MM-dd''T''hh:mm:ss.SSSSSS');
  52.   q0(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  53.   q1(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  54.   q2(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  55.   q3(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  56.   wx(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  57.   wy(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  58.   wz(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  59.   roll(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  60.   pitch(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  61.   yaw(i) = str2num(attitude.getElementsByTagName('q0').item(0).item(0).getData);
  62.  
  63. end
  64.  
  65. subplot(2, 3, 1)
  66. plot(d, x, '.-')
  67. grid on
  68. title('x(datetime)')
  69.  
  70. subplot(2, 3, 2)
  71. plot(d, y, '.-')
  72. grid on
  73. title('y(datetime)')
  74.  
  75. subplot(2, 3, 3)
  76. plot(d, z, '.-')
  77. grid on
  78. title('z(datetime)')
  79.  
  80. subplot(2, 3, 4)
  81. plot(d, x, '.-')
  82. grid on
  83. title('Vx(datetime)')
  84.  
  85. subplot(2, 3, 5)
  86. plot(d, x, '.-')
  87. grid on
  88. title('Vy(datetime)')
  89.  
  90. subplot(2, 3, 6)
  91. plot(d, x, '.-')
  92. grid on
  93. title('Vz(datetime)')
  94.  
  95. dnum = datenum(d - d(1));
  96. tnum = datenum(t - t(1));
  97. dstr = datestr(d, 'hh:MM:ss');
  98. tstr = datestr(t, 'hh:MM:ss');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement