HIMAN-359: Add 3x3 filtering to POT

fmidev · May 7, 2024 · 60399ca · 60399ca
1 parent f1e9405
commit 60399ca
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/himan-plugins/source/pot.cpp b/himan-plugins/source/pot.cpp
@@ -255,6 +255,13 @@ void pot::Calculate(shared_ptr<info<double>> myTargetInfo, unsigned short thread
 		POT = PoLift * PoThermoDyn * PoColdTop * PoMixedPhase * PoDepth * 100;
 	}
 
+	const himan::matrix<double> filter_kernel(3, 3, 1, MissingDouble(), 1 / 9.);
+
+	const auto smooth_pot =
+	    numerical_functions::Filter2D<double>(myTargetInfo->Data(), filter_kernel, itsConfiguration->UseCuda());
+
+	myTargetInfo->Base()->data = move(smooth_pot);
+
 	myThreadedLogger.Info("[" + deviceType + "] Missing values: " + to_string(myTargetInfo->Data().MissingCount()) +
 	                      "/" + to_string(myTargetInfo->Data().Size()));
 }