Changeset 1123 for trunk/src/transform_field

Timestamp:

Aug 21, 2023, 5:53:05 PM (16 months ago)

Author:

sommeria

Message:

various improvements

Location:

trunk/src/transform_field

Files:

• FFT2_detrend.m (modified) (2 diffs)
• signal_spectrum.m (modified) (1 diff)

trunk/src/transform_field/FFT2_detrend.m

@@ -68 +68 @@
         VarName=DataIn.ListVarName{ivar};
         z=DataIn.(VarName);
+        z(isnan(z))=0;% set to 0 NaN values
         rhs(1) = sum(sum(x.*z)); rhs(2) = sum(sum(y.*z)); rhs(3) = sum(sum(z));
         lin_coeff = inv(coeff)*rhs';
@@ -75 +76 @@
         spec2 = spec2(Iy,Ix);
         spec2 = spec2(ss,:);
-        %DataOut.(VarName) = log(spec2);
-        DataOut.(VarName) = spec2;
+        DataOut.(VarName) = log10(spec2);
+        %DataOut.(VarName) = spec2;
         spec_sum=sum(sum(spec2));
         kx_mean=sum(sum(spec2.*kx))/spec_sum;

trunk/src/transform_field/signal_spectrum.m

@@ -170 +170 @@
     if numel(ind_good)>WindowLength
         NbPos=NbPos+1;
-        if ~isempty(ind_bad)
+        if ~isempty(ind_bad)WindowLength
             sample=sample(ind_good); % keep only  non NaN data
             NbNan=NbNan+numel(ind_bad);