Wurzel 2 pi

CO_kirch_depth.m

@@ -60,12 +60,14 @@
         cosphi_up = z./sqrt(z.^2 + ((i_cmp-i_aper)*dcmp+h).^2);
         cosphi_down = z./sqrt(z.^2 + ((i_cmp-i_aper)*dcmp-h).^2);
 
-        % Compute diffraction hyperbola, /2 because data is not TWT but depth
-        r_down = sqrt(z.^2 + ((i_cmp-i_aper)*dcmp-h).^2);
+        % up and downgoing ray
         r_up   = sqrt(z.^2 + ((i_cmp-i_aper)*dcmp+h).^2);
+        r_down = sqrt(z.^2 + ((i_cmp-i_aper)*dcmp-h).^2);
+
+        % Compute diffraction hyperbola, /2 because data is not TWT but depth
         z_diff = 0.5*( r_down + r_up );
 
-        % Exit if diffraction ist out of data
+        % Zero if diffraction hyperbola ist out of data; else 1
         z_flag = (z_diff - z_max <= 0);
 
         %% Compute amplitude correction

COmig.m

@@ -114,7 +114,7 @@
         end
         if kirch_depth == 1
             [Kirchhoffdepth(:,:,i_h)] = CO_kirch_depth(filtdata(:,:,i_h), v, h(i_h), dt, dz, dcmp, aper_half, flag_interp);
-            COG = Kirchhoffdepth; %Was schlaueres ueberlegen
+            COG = Kirchhoffdepth/sqrt(2*pi); %Was schlaueres ueberlegen
         end
     end
 
@@ -125,7 +125,7 @@
 
     if v == vfinal % If loop reaches the correct velocity (estimated with constant velocity scan)
         % (estimated with constant velocity scan)
-        mig(1:z_len,1:ns) = sum(Kirchhoffdepth,3); % summing CO-Gather
+        mig(1:z_len,1:ns) = mean(Kirchhoffdepth,3); % summing CO-Gather
 
         %% Plot of the migration result
         mig_graphs('PolarPlot',mig(:,:),((1:ns)-1)*dcmp/1000,z,'CMP [km]','Depth [km]',sprintf('sum_vm%g',v))