收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

海洋环境模型 -海洋环境模型图片

[复制链接]
3 T( t. K: p- b) i& d1 W

我的专栏目录:

/ ~: |& ]% n! |$ F9 s! ]8 l 7 N( j6 Y" Q2 n( M

简介:

6 R4 D' s% v" x2 r4 @' t

海洋这个要素在游戏里占比越来越多,开放世界的游戏要求角色能上天入地下水。目前游戏里做海洋的方法就那几种。(1)预烘焙法(2)Gerstner wave(3)FFT海洋。预烘焙法可以是实现烘焙好DisplacementMap或者是FFT的运算结果。Gerstner wave可以在GPU或者CPU上算。FFT的话就是拿海洋频率模型算出Displacement。

" D; D- h# s5 [/ B! }" v

【1】基础环境搭建

! Y9 q& I7 E2 D% R) y5 B7 K: |

在开始研究之前我们需要先搭建起我们的环境。我选择在ComputeShader种完成各种计算,然后在顶点着色器种直接Sample前面的ComputeShader的波形计算结果。不要把波形的计算塞到VertexShader里。把波形计算独立出来还有个好处就是我能把波形的结果储存起来拿给其它效果使用,比如制作浮力部分的时候我们就需要知道海面波形的信息,如果塞VertexShader里就拿不到这些信息了。

7 O8 X# {- j9 S) s1 O) S

搭建ComputeShader的方法前面我的文章有提到,这里我就直接贴代码了。使用的引擎版本是4.21.0,如果引擎更新了新版本可能代码有一点区别。FFT部分我会给出4.22的代码。

+ U/ W) R. j1 E* C0 k$ e3 x

如果不是在Unreal中实现或者不想做这么复杂,可以直接跳过这部分。

0 @+ g, ?0 k1 X1 G5 i; x
, E/ i% M; _3 ^0 g
! e' k* V8 m8 @: U2 h

SDHOcean.build.cs

0 m3 E7 h; d1 F- q
// Copyright 1998-2018 Epic Games, Inc. All Rights Reserved. " ^( z" z- X/ E x) r5 B4 @' ~ , n( r* I- @" q% X5 ]* e3 ] using UnrealBuildTool; ) I! |! \1 J0 ]5 \: q! Q ; Q1 A9 Q7 K) _/ P) v2 B' g) o, X* M' X public class SDHOcean : ModuleRules& [' D" Z. O4 t$ i L {3 o9 E/ ?) m8 k* G7 a9 e/ ^, L public SDHOcean(ReadOnlyTargetRules Target) : base(Target) ' [$ _, g" ?+ c3 F1 ~8 r S {+ x0 M7 |4 T& {2 B, R. G) W m; Q PCHUsage = ModuleRules.PCHUsageMode.UseExplicitOrSharedPCHs; 3 x! f/ p3 K( \' A. `; J O 4 U! Q& ~: L3 ?, z p( I PublicIncludePaths.AddRange( / h3 K" y& l% ` w+ g) m new string[] { / i6 r; ]9 p s- F" n7 ?! W // ... add public include paths required here .... Z. ^9 K# O/ M+ ]# b- T }: F# z7 z7 i1 m* ]8 Y3 h );$ `$ E8 _9 D& [# N 8 u3 }' h$ y$ G. p% N9 e* } % z/ U) P# M5 Z" v7 n, k ]' W1 N PrivateIncludePaths.AddRange( ) w4 h& |7 g3 i6 X$ l7 E) u' ~$ v new string[] {1 d: R' t9 }% H% A3 y4 f // ... add other private include paths required here ...4 t5 T+ r' V4 W3 }" g: g& }/ ^: w } 8 \9 x4 N" E9 m2 h9 x );- {$ |* L) H0 L8 @ 1 U, N: D* y- ?/ ^$ h4 C( [ ' {- _0 k* ^- N' Y1 D* P* e% s! _) F PublicDependencyModuleNames.AddRange( * r" l. X, b. R) r; N new string[] U4 @& ~6 W# ?: ?9 U) {5 Y {4 Q9 _7 O3 J2 i7 f, t, r% t6 i" { "Core",) u' p! _- p0 ]2 ?- m4 @ "CoreUObject",( q3 U% j+ h) b "Engine", " @* [+ o8 g. U2 g2 A) m "RHI",0 O6 M" y7 \8 @+ w# g "Engine",+ V' F( F. o4 q( v+ W% F "RenderCore", ! q0 }7 p+ K/ {4 H "ShaderCore",9 g: r& s- M- p6 \0 b // ... add other public dependencies that you statically link with here ...4 b0 \7 a4 L7 H" l3 \) R+ | }4 X7 U( E: N3 H% g# r5 `3 d; B );, z5 ?" ^3 o# ^$ n ! l* _% N# U) r. p( n! D& V& e, i+ M + C& B" ~' m2 {( J/ n PrivateDependencyModuleNames.AddRange( 4 g( ~+ O4 f( H/ M- f! r2 I. { new string[]/ C/ l: M9 I1 R9 S {' a. \5 b8 j. |- K. Z3 d v# | f% [ "CoreUObject",6 m3 m+ ]" e# e% j8 u; c4 } "Engine",5 A+ ]$ a5 B0 t6 m* X, F& X1 t "Slate",6 s S1 _$ @1 _9 ~ "SlateCore", " E+ ~$ O9 ]- Y Y% B+ k* N! J$ ^ "UnrealEd",9 O, E* V" ?* |- [* L0 H "Projects",4 Y; y$ h" h& {8 c% C // ... add private dependencies that you statically link with here ... " j, `6 u4 n: O& N a! v' Y } ) I+ Z& _& b% i: q' l1 i; e d ); 3 O; A8 h, y" Y% [8 ]9 v# p. S, h0 j3 t8 G* Y! n* e 2 }/ ?8 ]: Z, V$ V. E1 e" { DynamicallyLoadedModuleNames.AddRange(+ q; \9 X# e( f9 {3 q0 ^ new string[]8 G: m& }7 O% v% A* F7 {. C { ( E5 p: _) @$ A9 H9 C7 X // ... add any modules that your module loads dynamically here ... 0 |: c. V$ i E } - e1 b9 o& D3 U6 L0 x ); & F% D" ]( z* N& `" N* J$ a6 T8 p } 1 s# a& P; }' ~" k) P) z% u5 @ }
* w9 d l( Y7 R( ^+ E; `

SDHOcean.h

3 Y/ Q' z- D8 q- M* z8 z
// Copyright 1998-2018 Epic Games, Inc. All Rights Reserved. 1 L: m+ v/ e0 ]. n' T: M5 N. }2 S( {% P% c; O* z #pragma once2 x8 T: o q5 ? 8 f/ l* t; f* z0 j+ ~9 K E6 T #include "CoreMinimal.h"5 ~0 u2 M" ~1 T& c! N' R; H8 | s- t, X #include "Modules/ModuleManager.h" W' q' U Q8 D0 i& e7 M6 C #include "Interfaces/IPluginManager.h" g" f" {& {. Q5 p; I* A% e; f #include "Misc/Paths.h"- e( n5 W1 x8 g4 ? #include "Modules/ModuleManager.h"* H2 U) M) F( E( ^3 J) K0 E * |% C b5 [! k @ N) U9 k$ i: a- f class FSDHOceanModule : public IModuleInterface + y6 q/ e! V: f( M' o) n { : p+ }' w% q+ h public: 7 r" J7 Q V0 h) K7 x$ {' C* s6 s! s% f1 M% t; o+ D' ]1 d /** IModuleInterface implementation */8 ?2 T8 t; s9 o# p' \$ O virtual void StartupModule() override;% b/ q, _- D* ~7 r& U/ ] virtual void ShutdownModule() override;; ]$ M4 k$ Z& B };
q& E s4 X' B+ N1 l) e

SDHOcean.cpp

; p' x7 x0 }8 V( q, A
// Copyright 1998-2018 Epic Games, Inc. All Rights Reserved.8 H& C! p- \: ~& V" ]- L 8 k& K5 k# u% I9 I& D0 z #include "SDHOcean.h" 1 f5 ^& p6 i( z) ~6 ?9 b8 I) E/ g/ Z- d- T #define LOCTEXT_NAMESPACE "FSDHOceanModule" * `1 g- X9 u" c4 W8 ^9 } p' R, U- |) \* |0 d6 a0 g" s; W( H void FSDHOceanModule::StartupModule()0 K( Y, n0 E' V) {. z2 W& Z! h { - {5 }, X: t y% R' k // This code will execute after your module is loaded into memory; the exact timing is specified in the .uplugin file per-module 5 \. z3 k, |! ]: u4 H' t% K: f8 I FString PluginShaderDir = FPaths::Combine(IPluginManager::Get().FindPlugin(TEXT("SDHOcean"))->GetBaseDir(), TEXT("Shaders")); 9 l, \) v) I( s) ` AddShaderSourceDirectoryMapping(TEXT("/Plugin/SDHOcean"), PluginShaderDir);1 u/ R, m' l* q } ' F3 x" S' z; [' P( g 9 \( ~2 A' J5 u( K, r" G7 d void FSDHOceanModule::ShutdownModule()3 B) m5 X+ z6 d) L. D$ S {# g; s' i9 E& T) e0 I // This function may be called during shutdown to clean up your module. For modules that support dynamic reloading, 6 z2 {; {/ T3 }7 A: { // we call this function before unloading the module. ! O+ O: v. F$ S" W7 I } % _' I Q; g: a9 G) [ X/ [/ {' ~ 6 J; z5 H+ F( f! z( p #undef LOCTEXT_NAMESPACE 6 q0 S: y& \! h6 h - @' Z& L' ]! z M+ W IMPLEMENT_MODULE(FSDHOceanModule, SDHOcean)
* ]4 z8 @; T; F8 X, z

Ocean.h

' ] W$ L4 R$ B% `% p7 j- ]
#pragma once/ Z( ^ k8 K8 h( e+ d% p2 @5 Y% m , D, O1 ~. J$ W9 y #include "CoreMinimal.h"4 s) O" z; v! j6 L #include "UObject/ObjectMacros.h" / K) a9 p4 ^! X #include "Runtime/Engine/Classes/Components/ActorComponent.h" ( k* H/ n- p" n0 V, |* A6 }" v #include "Engine/Classes/Engine/TextureRenderTarget2D.h"# w. A# F0 D! n6 V, Y #include "Ocean.generated.h"; Y$ `7 |, M; G; G& V 2 x2 s8 Q3 Y( e5 Y typedef TRefCountPtr<class FRHITexture2D> FTexture2DRHIRef;3 M1 y$ \# E! V typedef TRefCountPtr<class FRHIUnorderedAccessView> FUnorderedAccessViewRHIRef; ; M: R" C4 P- J% ]) y2 C' q" Y8 m typedef TRefCountPtr<class FRHIStructuredBuffer> FStructuredBufferRHIRef; 9 G4 `* k, {' G class FRHITexture;% Y4 U/ @) z# Q a: y& \ class FRHIUnorderedAccessView;2 }3 f/ X4 T$ \; M; W class FRHICommandListImmediate;1 u- N/ b, m0 d& C- ^4 r; P8 T * {& L( o4 m/ N: C2 k L! ^4 P b USTRUCT(BlueprintType) 4 N. _& }& q* R" [/ ^# W struct FOceanBasicStructData_GameThread( s' f0 c+ Y% f" I- ^$ ^% p% v { : o9 T" A8 Z' }+ T GENERATED_USTRUCT_BODY() - _( O* W) Z, q3 E: l" l! v4 }$ O2 z+ J5 Q8 u! v FOceanBasicStructData_GameThread(){}8 ^* Q D: f, @" K * W6 T, V( Z+ d% |, y& U! P UPROPERTY(BlueprintReadWrite, EditAnywhere)& e$ T% Q# c/ [8 @8 \ FVector4 OceanTime_GameThread;4 s: B- o H2 u, q6 n }; * K3 e3 m w$ p) w/ w' g 8 @, y3 h* J2 x5 q UCLASS(hidecategories = (Object, LOD, Physics, Collision), editinlinenew, meta = (BlueprintSpawnableComponent), ClassGroup = Rendering, DisplayName = "OceanRenderComp") - t6 G) {4 X& U) {6 i( y/ q class SDHOCEAN_API UOceanRenderComponent : public UActorComponent1 V M" l1 N& _' O- w6 h0 R; O { . l7 I. s! |$ t& @ GENERATED_BODY(), x( F- {4 P7 C3 i' ~ 2 l7 k* F5 q( [# t public:4 c9 ^; S) b% H. W1 i ?! |) q8 |2 H1 f UOceanRenderComponent(const FObjectInitializer& ObjectInitializer); . M* |9 I" I. F5 [: V: Q* o //~ Begin UActorComponent Interface.' P& F+ v; a: i3 {3 x virtual void OnRegister() override;4 I) }# G8 T5 o8 x$ V& K: z3 U virtual void TickComponent(float DeltaTime, enum ELevelTick TickType, FActorComponentTickFunction *ThisTickFunction) override; / O# O2 G6 q( K$ t" q9 m+ h( p; r/ [0 y- r1 s6 J8 w UPROPERTY(BlueprintReadWrite, EditAnywhere, Category = "OceanComponent")$ ~% T5 p w! ~. m UTextureRenderTarget2D* OutputRenderTarget2D; - l4 o# }- k7 w y0 x# Y5 o1 R/ N9 a //UniformData for Ocean render - J5 H4 V& F8 n. w4 r1 v. Y- j UPROPERTY(BlueprintReadWrite, EditAnywhere, Category = "OceanComponent") " O! X0 p' c5 }4 G4 L. l FOceanBasicStructData_GameThread OceanUniformDataBuffer; : K- R4 ^! c& y7 ]! e/ O 1 z6 @7 s# {' g( W int32 TargetSize;; A2 O# Y+ f% y# P$ X5 l6 B ETextureRenderTargetFormat RenderTargetFormat; . } y0 x) T$ {3 W' Q+ J* P" G2 V3 i0 {( f" E# l! } private:$ N ]. Y0 _: k2 Q# Z9 q $ D) l: r& a) z+ I! ~ //Render ocean render thread ; X1 D! B& _1 |" Q. E6 i+ l" u void OceanCalculating_GameThread();2 u! w. n5 C# ^" k void OceanCalculating_RenderThread& F* m( o4 Q2 X1 L7 V' F ( 7 z. ]0 y0 H/ D5 z FRHICommandListImmediate& RHICmdList,4 y- l' g4 ^! |6 F ERHIFeatureLevel::Type FeatureLevel, % q, T. h3 p; \. F% `$ r* c FRHITexture* OutputRenderTarget,3 |( u* s7 u0 A1 O int32 SurfaceSize, " V) G6 e+ D. t" p const FOceanBasicStructData_GameThread& OceanUniformData+ \, F: g/ P" ? );: S- e1 r, H7 z* j' y: o5 o( I 7 v/ k6 K) z9 [% z3 u+ H FTexture2DRHIRef OutputTexture;3 {) e4 L, T) _6 v4 S/ c2 u FUnorderedAccessViewRHIRef OutputTextureUAV; ! y2 d3 X! M. E- w7 v. L }; 7 h0 f: y4 q, a% b2 K
7 ?: A9 _2 O, J ^

Ocean.cpp

. i3 L3 Y( Y" G) v
#include "SDHOcean/Public/Ocean.h"/ @8 u2 n2 X+ {0 j8 @ $ k/ {* y2 H. | #include "ShaderCore/Public/GlobalShader.h"; h9 z9 l1 ^& C4 P# P ' @ o, V& I$ V' Y. @- l #include "Classes/Engine/World.h" ) A$ K9 l2 U0 k #include "Public/GlobalShader.h" ( U+ B; j7 l4 f7 t G) P$ M# W$ n #include "Public/PipelineStateCache.h"( a5 }( a7 [9 M4 E7 f! b #include "Public/RHIStaticStates.h" ; `: Q; j+ Z5 [0 g& F4 K7 D% m #include "Public/SceneUtils.h"6 X0 K: X5 ~7 K' O8 C. K #include "Public/SceneInterface.h"& ~: l8 n; B( n2 f #include "Public/ShaderParameterUtils.h" , |9 k$ Q" }* }0 h$ f0 j #include "Public/Logging/MessageLog.h" _3 @. o( L0 U7 V* L$ l #include "Public/Internationalization/Internationalization.h"/ X+ U& |. V% u" K" |" y) s #include "Public/StaticBoundShaderState.h" x7 o: Y* T6 r" [, X. @9 N( v #include "RHI/Public/RHICommandList.h"# s* g4 Q0 S. V# ]& Q4 A9 _ #include "RHI/Public/RHIResources.h"0 L8 o2 b2 f1 }; f1 @ #include "Engine/Classes/Kismet/KismetRenderingLibrary.h"5 W8 X' Q/ u9 r6 l, Q' A #include "Runtime/Engine/Classes/Kismet/GameplayStatics.h" ' q3 r+ e+ t6 E" M7 B; g# j, {% u% N( J #define LOCTEXT_NAMESPACE "SDHOcean"" v9 z2 D" J5 |6 t 1 y* Y0 b# _$ m7 s7 i" _ BEGIN_UNIFORM_BUFFER_STRUCT(FOceanBasicStructData, )5 Z/ V$ o/ F% {. A3 F UNIFORM_MEMBER(FVector4, OceanTime) 4 E6 ]) b, c2 b" t END_UNIFORM_BUFFER_STRUCT(FOceanBasicStructData). a, i7 Z8 z, h3 F* h J& z ; D" J$ L+ q- h2 ^ IMPLEMENT_UNIFORM_BUFFER_STRUCT(FOceanBasicStructData, TEXT("OceanBasicStructData")) 5 C9 l8 Q' q' v$ _; D* S: y- A9 a class FOceeanCSShader : public FGlobalShader , b) D9 R7 E" H$ n6 M {, ]8 O8 @6 {6 e- s 8 S1 h0 u9 U6 M0 h" K4 c DECLARE_SHADER_TYPE(FOceeanCSShader, Global)% I: o8 w- O: A/ h. G0 u) Q" z # G- V* ^6 s1 K2 s public: L& [" T3 a8 G1 N- j 0 X" K3 W# ?$ _0 ? FOceeanCSShader() {}; X$ i0 J4 Z* z FOceeanCSShader(const ShaderMetaType::CompiledShaderInitializerType& Initializer)! [; X# E3 ]+ f* i: o : FGlobalShader(Initializer) ) a, m( ^: A1 M+ F& k) P$ ?& j3 e1 { {8 z7 l/ [) c6 i k% ~5 v //TODO Bind pramerter here ! i. Q% T' s" l9 d: [: h OutputBufferSurface.Bind(Initializer.ParameterMap, TEXT("OutputBufferSurface"));/ ^/ X8 e/ I" J SurfaceClearColor.Bind(Initializer.ParameterMap, TEXT("SurfaceClearColor")); * J _* H$ H6 M }) k9 Y" F+ ]6 q0 w //----------------------------------------------------// - |; }( q2 z5 e0 J static bool ShouldCache(EShaderPlatform PlateForm)$ e4 v) y+ J0 j3 b) O {1 h* B8 Y, b: a# { return IsFeatureLevelSupported(PlateForm, ERHIFeatureLevel::SM5); + a; [( H% V8 f } 7 C) |% o3 n) s$ p, j //----------------------------------------------------//- G9 w; E9 o1 z" ?: b( i6 D+ @ static bool ShouldCompilePermutation(const FGlobalShaderPermutationParameters& Parameters) & R! H; @+ o% H, ^ { 9 n, e1 s; k/ S; B4 D5 b0 K return IsFeatureLevelSupported(Parameters.Platform, ERHIFeatureLevel::SM5); ; W+ |' I7 J: O0 g0 h }& P! o, q5 M6 `: I7 y& b2 B) E' o //----------------------------------------------------//( e/ t4 e3 R8 g: i9 ] static void ModifyCompilationEnvironment(const FGlobalShaderPermutationParameters& Parameters, FShaderCompilerEnvironment& OutEnvironment)$ K, `& j8 F6 l1 p2 Y {7 i; M0 n1 P7 b3 T, ~4 S FGlobalShader::ModifyCompilationEnvironment(Parameters, OutEnvironment); ' e* v6 T( ?' f8 q: c //Define micro here . ~4 [7 X0 ^5 T; Y4 }0 @ //OutEnvironment.SetDefine(TEXT("TEST_MICRO"), 1);7 j* \0 j" N$ _4 ~- v% Q }8 B, d" N# L2 ]) ]2 I( k //----------------------------------------------------// ) K/ V+ s+ i" v) j* }: a 4 S# _. h6 P8 j) D void SetSurface(FRHICommandList& RHICmdList,1 c7 @! O8 g9 g+ f O( B" f( j: P FUnorderedAccessViewRHIRef& OutputUAV,. Q, m) _* u, ~/ E% _8 J$ W" P2 z5 F const FLinearColor ClearColor; _& @" ~' t$ q$ z ); T( M7 r+ ?( V2 Q {9 D; ]4 K( H- r* t$ H6 j //set the UAV3 i5 u; }" ~9 ~0 i) n3 V FComputeShaderRHIParamRef ComputeShaderRHI = GetComputeShader(); 0 P' A6 }' p! l2 c if (OutputBufferSurface.IsBound()) , O o! {4 c5 |8 s RHICmdList.SetUAVParameter(ComputeShaderRHI, OutputBufferSurface.GetBaseIndex(), OutputUAV); 8 J) y* b- W0 x" B0 K, M" t& S if (SurfaceClearColor.IsBound()) 1 N4 c& x- U+ K( D( n //RHICmdList.SetShaderParameter(GetComputeShader(), SurfaceClearColor.GetBufferIndex(), SurfaceClearColor.GetBaseIndex(), SurfaceClearColor.GetNumBytes(), ClearColor); 8 C2 z; ?# n! {0 Q2 p) i; }3 R1 p# K0 | SetShaderValue(RHICmdList, GetComputeShader(), SurfaceClearColor, ClearColor);* K4 l0 i. `1 ? }2 }. O* f& P5 g% p& u 0 E( C0 e5 |2 f8 L7 V* u, |6 _, L void SetOceanUniformBuffer(FRHICommandList& RHICmdList, const FOceanBasicStructData_GameThread& OceanStructData) 6 {8 d$ G* @) @' P8 X {4 G' z2 d, o6 V1 }! k! r FOceanBasicStructData UniformData; 1 v F% j) P) c8 u5 K8 C UniformData.OceanTime = OceanStructData.OceanTime_GameThread;6 E# |0 D" e1 m5 W7 }1 e8 S SetUniformBufferParameterImmediate(RHICmdList, GetComputeShader(), GetUniformBufferParameter<FOceanBasicStructData>(), UniformData);' Y* M6 t8 d$ |2 E } ' c9 k8 u- g+ k. ]1 j- d/ A8 s/ f5 i4 y! W6 I: O. k void UnBindBuffers(FRHICommandList& RHICmdList) 1 T- A0 W/ f/ M( K {2 m. y' ^# Z# F' E, S* [ FComputeShaderRHIParamRef ComputeShaderRHI = GetComputeShader(); 1 P$ @/ ~2 N6 t* d% ? 5 n% A7 ~. ~ J) C0 R if (OutputBufferSurface.IsBound()) & R/ d0 ]& Q& c& u: W RHICmdList.SetUAVParameter(ComputeShaderRHI, OutputBufferSurface.GetBaseIndex(), FUnorderedAccessViewRHIRef()); / z' u9 l( N' ~* ^- C8 P8 X } " u* J- A1 y9 { 8 }/ |8 v G o; s: Q, Q5 P virtual bool Serialize(FArchive& Ar) override8 u ^ t+ g1 ^% C7 y3 ?1 Z0 t { ; Z# q6 w3 D/ x% R9 p U bool bShaderHasOutdatedParameters = FGlobalShader::Serialize(Ar); $ y/ v$ g9 _& c% E& b E1 z. l //Serrilize something here. J/ b- a" n+ T' K# Q2 m! r Ar << OutputBufferSurface << SurfaceClearColor; 6 j& N+ I2 R" N, H+ i, ~& P* Y return bShaderHasOutdatedParameters;4 G3 Q% g6 I) ` } ' A1 x- p$ k( }3 T0 G' m( s3 |; r7 O8 v8 x+ a- A private:6 g; T1 W2 n7 q0 i$ f ; X5 [( p5 c y1 z9 N4 `9 h FShaderResourceParameter OutputBufferSurface;7 e0 Q( ~3 ?; g! D" c# N FShaderParameter SurfaceClearColor;* \1 @3 O# x: l, z9 D" E; y/ l }; 4 Y1 w2 a5 s2 x& q: W* T# X( \% _9 S. u IMPLEMENT_SHADER_TYPE(, FOceeanCSShader, TEXT("/Plugin/SDHOcean/Ocean.usf"), TEXT("OceanMainCS"), SF_Compute) 0 y* N; e/ U/ ?& q3 @1 {4 e K$ f x5 E; n! N& T void UOceanRenderComponent::OceanCalculating_RenderThread $ e7 t5 ~! G% e* v( d (9 X0 @. j$ h! c6 G1 }. x FRHICommandListImmediate& RHICmdList, & a$ @7 J( S6 z( @- z# c" L ? ERHIFeatureLevel::Type FeatureLevel,8 ?) J1 [* {8 v1 P FRHITexture* OutputRenderTarget, 5 H6 {. v5 O3 Q7 `& D int32 SurfaceSize, ; L3 [' `4 k% T* x1 ~ const FOceanBasicStructData_GameThread& OceanUniformData 8 T* |) p" }9 _2 g6 r/ E ), ?9 g! g; {! W6 ^, o3 V { c( Z% a) n# b. ?( x% }! } check(IsInRenderingThread()); - q' Y& e2 z; ^; I, h. p check(OutputRenderTarget);7 }+ ?0 j, s- H4 i* w2 ~0 K + v o3 E. M$ w8 n1 ^ TShaderMapRef<FOceeanCSShader>OceanComputeShader(GetGlobalShaderMap(FeatureLevel)); 5 p3 y2 C6 I! D3 W( V% `* ?9 P7 O RHICmdList.SetComputeShader(OceanComputeShader->GetComputeShader());" d. }, T7 G# f9 \! M1 s " q5 x; |2 s- a- z7 R if (OutputTexture.IsValid() == false)1 Q$ t4 d+ j: d9 R+ | {7 J* p, Q8 [8 D9 Z; R7 j' k9 G9 T if (OutputTexture.IsValid()) : D u$ M2 c0 H1 R+ q OutputTexture->Release(); 5 w0 I" r8 K3 y6 Y0 y$ q if (OutputTextureUAV.IsValid())( a$ U9 j2 ], n1 `8 R" H4 N0 x OutputTextureUAV->Release(); - b8 v( c2 Y5 k$ D6 g7 G1 N: w( G) {, B FRHIResourceCreateInfo CreateInfo; 8 `% q3 A4 Z4 B OutputTexture = RHICreateTexture2D(SurfaceSize, SurfaceSize, PF_FloatRGBA, 1, 1, TexCreate_ShaderResource | TexCreate_UAV, CreateInfo);4 n4 ^. Y8 S. v" e OutputTextureUAV = RHICreateUnorderedAccessView(OutputTexture); ; [( {6 f' e$ R, l3 a" h2 y$ F } / s- O. @; R# w6 o9 n! G$ @! P/ U3 q3 p; ~! K OceanComputeShader->SetSurface(RHICmdList, OutputTextureUAV, FLinearColor(1,1,1,1));' q4 A8 Y& O, U9 V6 m! m OceanComputeShader->SetOceanUniformBuffer(RHICmdList ,OceanUniformData);5 z1 u$ P) A5 a " N; W2 w1 _9 a0 w$ O! _7 _ DispatchComputeShader(RHICmdList, *OceanComputeShader, SurfaceSize / 32, SurfaceSize / 32, 1); : u, U2 L! R; c: t& Y OceanComputeShader->UnBindBuffers(RHICmdList);$ k4 A2 I6 j7 Q 6 \; i, d0 x5 y, `7 \+ g RHICmdList.CopyToResolveTarget(OutputTexture, OutputRenderTarget, FResolveParams());, B3 i; N' D6 n: I5 W2 d9 ~ //FRHICopyTextureInfo copyinfo(SurfaceSize, SurfaceSize);6 \3 @7 z1 L$ z, q/ c //RHICmdList.CopyTexture(OutputTexture, OutputRenderTarget, copyinfo); % U8 V9 m" V+ f5 H } + u. u( A& U9 T; \/ B3 c 7 R* R! A3 {' J) u6 G UOceanRenderComponent::UOceanRenderComponent(const FObjectInitializer& ObjectInitializer), |4 S6 p) S; ^+ [0 ^- F* @* b2 c :Super(ObjectInitializer) $ ]7 x& ~& p3 [3 ]; d# `+ ] {$ h$ l" I( n2 ?2 U. Q PrimaryComponentTick.bCanEverTick = true;2 a8 o- m B& a* H( ~ bTickInEditor = true;6 R9 V# A! h" W, |! B1 Q bAutoActivate = true; , q# q+ q! B7 O7 w5 O- r/ a( |' ]* m" y/ [. X3 J3 Z0 q/ ] RenderTargetFormat = RTF_RGBA32f; 1 _5 N' I! O' y2 @$ v0 _7 o' q* D }; K) K/ |0 W; G; u0 z2 `5 C 4 I* H# v) D. W/ }) l" B void UOceanRenderComponent::OnRegister() 7 f% Q$ j! T% _& L {6 [$ V9 b8 `! L+ i8 F Super::OnRegister(); 8 D. J6 ?* S/ S8 I }7 P$ M/ ] S0 k( ` # a2 N4 ]/ s+ t N7 u% b1 n void UOceanRenderComponent::TickComponent(float DeltaTime, enum ELevelTick TickType, FActorComponentTickFunction *ThisTickFunction), z1 [. T# U; o) r {3 {: i1 a4 c: Z7 e Super::TickComponent(DeltaTime, TickType, ThisTickFunction);) Y0 Q8 r: W; Y. | / r1 Z- M j9 B- A0 y" Y //Tick render the ocean3 l: l1 l$ K4 i, F$ P! D- K OceanCalculating_GameThread(); 8 `2 l, D% ~8 I3 W$ b5 c- a! \. O }5 d6 e1 S7 j; I8 |" g& B 3 Y+ E, s7 ]0 ] void UOceanRenderComponent::OceanCalculating_GameThread() 1 B: H* g2 [1 E- b7 j: T9 s: U {' W2 g, M- h: ^5 t; y UWorld* world = GetWorld();. U l' ^) z+ G8 _ n ERHIFeatureLevel::Type FeatureLevel = world->Scene->GetFeatureLevel();4 U2 ~5 J' I: [+ @, O$ S/ D# w checkf(FeatureLevel == ERHIFeatureLevel::SM5, TEXT("Only surpport SM5"));8 R, r* y- K( S% I: K1 \ 7 X+ u4 I' S3 Z# D5 b9 u if (OutputRenderTarget2D == nullptr) return;6 o( Z% n1 A& d) C/ ? . b. ~6 @" Y& n: m% c' z9 w3 N //Using front RT to render,back buffer store last frame imformation 2 b1 E1 T: \" v. M8 O7 @ w" k, c UKismetRenderingLibrary::ClearRenderTarget2D(world, OutputRenderTarget2D); 4 B2 @' d; f. J 0 Q1 a) Z. c1 ~% o FTextureReferenceRHIRef OutputRenderTargetTextureRHI = OutputRenderTarget2D->TextureReference.TextureReferenceRHI; ! F: W3 _% n' _. v9 s- g2 |4 A2 R! b! f) B5 ]5 C0 e checkf(OutputRenderTargetTextureRHI != nullptr, TEXT("Cant get render target %d texture"));3 n5 T0 K# t6 c" Y% t2 {3 ^ ' N0 G: U$ y0 ` FRHITexture* RenderTargetTextureRef = OutputRenderTargetTextureRHI->GetTextureReference()->GetReferencedTexture(); * X( U) q' x5 V3 A5 H TargetSize = OutputRenderTarget2D->SizeX;3 _6 k: u0 f" A; v" q a# |5 I! n5 I% j* J+ W( @, t0 O. o //Update the uniform buffer& t+ D# {' `1 a% l' ~" m OceanUniformDataBuffer.OceanTime_GameThread.X = UGameplayStatics::GetRealTimeSeconds(GetWorld()); 4 P: u% _4 C# x6 X& }! y/ w; h6 a/ Z! l! ^& I% } ENQUEUE_RENDER_COMMAND(OceanRenderCommand). U% Y4 o: T3 A. J ( 7 m7 n6 c/ f8 l* m+ W# S* c$ _) d [FeatureLevel, RenderTargetTextureRef, this](FRHICommandListImmediate& RHICmdList)9 n0 K7 F# h. G3 o" w$ u) K1 l { " d1 a; w, {. `" H: b% ?, ~$ ~ OceanCalculating_RenderThread % Y% u* u! q4 a" ` ( + z. |7 v- q5 o" @. t RHICmdList,: S, E- C% |: q% T; z FeatureLevel, 1 }1 ~0 n- }& H. |$ v RenderTargetTextureRef,% v+ H9 R( _" z this->TargetSize, , v: S3 w6 W! m6 { this->OceanUniformDataBuffer: U, L. e. l o ); g4 }% f0 B9 W( P" Z r8 h' v4 a5 p }: l/ [+ I. o0 X$ Q+ ^6 a ); ) J) w$ v+ Q2 N# P/ W; C+ | }* z4 z5 g+ b a' o7 x1 E ) A1 p3 @' \6 m #undef LOCTEXT_NAMESPACE
S! b7 \* `# ^' ?

Ocean.usf

+ o+ V ? `; ~4 A
#include "/Engine/Private/Common.ush" # Y/ h9 x* j- U& J" _ # H' X/ J7 q9 `5 F RWTexture2D<float4> OutputBufferSurface; + y, C9 q7 O6 c3 U float4 SurfaceClearColor;- {4 C% y% R$ x3 p3 j9 S . Z2 L, I7 ^6 ?$ w1 O struct OceanUniform- [) V+ J- J- x; q {% ]* H9 n F/ N$ j `' e5 n float Time; ( {4 ]$ } v J }; * V7 l+ \$ M. M, w: d1 Z. n void InitOceanUniform(out OceanUniform uniformval). M/ ^" U. [, r# c7 q {" X- v( m' Q! O P1 `7 P5 | uniformval.Time = OceanBasicStructData.OceanTime.x;: p' V! V& K& u+ x9 b+ [: f! x }9 b7 p9 p- X1 p A8 y 5 G( W% l3 O: B$ E+ j/ K [numthreads(32, 32, 1)] ! t0 b2 m6 T) m) ~: S8 A ^ void OceanMainCS(uint3 ThreadId : SV_DispatchThreadID)% G2 u" P* g' E0 N: Q4 _* R {" B9 j0 v' T6 G1 X: y+ N) R //Set up some variables we are going to need ) s+ Y- O3 n) y2 Q- Y* Z6 R //The size of outputsurface and input surface is same$ O4 k h9 F3 z6 ~9 }) }; H float sizeX, sizeY;" ~. u, ^+ @& R: P( a' p* K2 p6 ~ OutputBufferSurface.GetDimensions(sizeX, sizeY); ( X9 F5 W$ p' D# m ; }' ^- a; }1 u# C; A+ ? OceanUniform OceanUniformVal;3 N: r6 R2 c9 J) }$ ^: v) l InitOceanUniform(OceanUniformVal); 4 P$ _0 o/ e+ e$ R# S. J# c4 J4 m( g2 ^% d" ]0 O float2 iResolution = float2(sizeX, sizeY); # H( `4 [+ v0 k y float2 UV = (ThreadId.xy / iResolution.xy) * 50.0f;, m( c% B, H6 Z2 h 4 N1 q6 U; S8 h+ n float4 Output = float4(1.0f, 1.0f, 1.0f, 1.0f);2 z( _( j6 k9 K# T* v! [+ ]# B: ^. c 8 R: w2 @$ J+ ^7 A+ m5 Y+ \ Output.xyz = float3(sin(UV.x + OceanUniformVal.Time), 0, 0); + u% A$ x* G( `# I; c k, B* M8 T' e, C OutputBufferSurface[ThreadId.xy] = Output; 5 z$ b$ o( f2 f1 e; c1 m }
5 ~( j; l3 \0 H7 r0 K' t9 g6 ~

我这里做了个UniformBuffer然后拿到系统的时间,把系统的时间变量塞到我的Compute shader中。

) `% K3 Z2 r: \5 O
& y) u- `) `* n" n- n( C

直接把值采出来连到VertexPositionOffset上就可以把我们的ComputeShader的结果传到顶点着色器了。

8 P9 V2 P8 w" T" y" ]! h) {3 P, W. h

【2】Gerstner Wave

; }. n3 A+ }! T6 n

在Gerstner Waves之前,先使用正玄波变形的方法模拟。对正玄波进行变形属于经验性的方法,把正弦波变形让它的形状更接近水浪。

- ?1 E' }# {1 N1 k, H$ E
+ J( ?" [% e* E q" A' C( G

F(x)=2(\frac{sin(x) + 1}{2})^{k}

% g/ b9 V9 u+ E) }

下面先来制作正玄波水面

]+ y. {( o6 m2 p, b( ]! {) L
" P8 k6 {8 z2 L h/ V, n

于是乎我们可以得到如下效果

. }8 U' C7 p& y- X( n3 h6 c) B3 T) N
# D9 X6 ^8 g9 Y7 A1 }" y
+ ?2 n; B. I; f& v
- d( t0 r6 U1 X+ J! Y' J

有了基础的波形后,剩下就是让海面变化更丰富。想让海面变丰富那就是多叠几层波

; N0 w# h/ v. E' t; N' H
1 k% B4 G7 C1 P4 l S& `# M
$ B% w( G ]6 x' s
( e5 u( R/ D6 R5 X; K; I6 K
Output.z = 0.5 * pow((sin(WaveLength * Speed + dot(direction, UV * 0.8) * WaveLength) + 1) / 2, 2.5f);7 d W7 h, y5 q! D# v, M0 f) n- i Output.z += 0.2 * pow((sin(WaveLength * Speed * 0.8f + dot(float2(0.8, 0.1), UV * 0.9) * WaveLength) + 1) / 2, 2.5f);1 p$ |" ~# w" _5 p) C9 m Output.z += 0.15 * sin(WaveLength * Speed * 1.2f + dot(float2(-0.8, -0.1), UV) * WaveLength * 1.3f); / f4 y- a9 l& _, E1 ?9 ^ P7 L8 S Output.z += 0.1 * sin(WaveLength * Speed * 1.2f + dot(float2(0.6, -0.5), UV) * WaveLength * 1.5f); 9 d+ {! A" f3 f# @9 |7 V/ X Output.z += 0.1 * sin(WaveLength * Speed * 0.5f + dot(float2(0.5, -0.1), UV) * WaveLength * 1.5f); . B- T6 ?9 A9 r: V% p1 m/ H" R 5 k" W5 {. R8 b0 T1 z Output.y = 0.5 * pow((cos(WaveLength * Speed + dot(direction, UV * 0.8) * WaveLength) + 1) / 2, 2.5f);: o2 I+ W: C5 s1 F8 E Output.y += 0.2 * pow((cos(WaveLength * Speed * 0.8f + dot(float2(0.8, 0.1), UV * 0.9) * WaveLength) + 1) / 2, 2.5f); " `$ ^$ C- s+ { v. F+ _+ n1 X Output.y += 0.15 * cos(WaveLength * Speed * 1.2f + dot(float2(-0.8, -0.1), UV) * WaveLength * 1.3f); ; e. y: V7 H8 M( i1 T$ b" v6 P Output.y += 0.1 * cos(WaveLength * Speed * 1.2f + dot(float2(0.6, -0.5), UV) * WaveLength * 1.5f);$ g. s% Z* [) b0 O" G Output.y += 0.1 * cos(WaveLength * Speed * 0.5f + dot(float2(0.5, -0.1), UV) * WaveLength * 1.5f);8 G3 G) R/ A2 k. i5 A. F3 S , Q$ W6 X- d: Q" C P$ l Output.z = 0.5 * pow((cos(WaveLength * Speed + dot(direction, UV * 0.8) * WaveLength) + 1) / 2, 2.5f); ' d, U" s* L3 f Output.z += 0.2 * pow((cos(WaveLength * Speed * 0.8f + dot(float2(0.8, 0.1), UV * 0.9) * WaveLength) + 1) / 2, 2.5f); # h8 s* S6 \' b! x; B Output.z += 0.15 * cos(WaveLength * Speed * 1.2f + dot(float2(-0.8, -0.1), UV) * WaveLength * 1.3f); r2 E$ K' J/ \! P8 _% _$ a Output.z += 0.1 * cos(WaveLength * Speed * 1.2f + dot(float2(0.6, -0.5), UV) * WaveLength * 1.5f); : {" I3 ^ E* m, X$ } ~: @. Y) V Output.z += 0.1 * cos(WaveLength * Speed * 0.5f + dot(float2(0.5, -0.1), UV) * WaveLength * 1.5f); 7 q6 D/ d' o4 @9 w+ |3 D4 [4 ?% |& v) r OutputBufferSurface[ThreadId.xy] = Output;
1 r3 \3 o0 Y+ U

可以看到正玄波水面波浪比较平,无法模拟出水波的波峰陡峭的特点。因此我们需要使用新的模拟模型:Gerstner Wave。

* `1 n8 N7 u1 u1 h
6 t8 J' v' f0 q6 ?. u
9 G( E( c' H% E0 H
# Z8 S4 Q: e* a7 T8 |# O1 T7 S8 O, g' F

PositionOffset:

- V' E) ^1 V8 a7 S
j5 Q" A& _- j6 E% h

Normal:

. v9 N, @1 P" x5 d% B) P
, U/ U$ B- |; V% `) W! O

Gerstner Wave是周期重力波欧拉方程的解, 其描述的是拥有无限深度且不可压缩的流体表面的波形 。

v+ A! ?. x: U$ G* _6 p% C
7 P9 T- ]+ B& B+ U/ D1 d+ c# |! \
% ~) [* ^$ M: Q; T3 h
* t! C8 A; S; R

代码如下:

: f G( p* g3 T$ N2 ~
; H: A/ |7 L2 v9 P' ]8 Y q) A

想要更好的效果可以优化下参数和多叠几层波,反正Computeshader里算这种东西很快的啦。

4 Q5 q! Y" b/ _3 q, S! a& c5 O7 f

【3】FFT海面理论推导

: Y) {! h. e4 d* E

FFT海面的核心思路是我们通过一系列测量得到真实海面的波的频率然后把这些频率通过FFT变换到时域然后计算出置换贴图。下面来公式推导

3 N0 f4 S p( Z s2 r/ p

设波的高度为水平方向的位置和时间的关系 h(X,t) 。在水平方向上 X= (x,z) 根据Tessendorf J.2001的论文我们可以得到如下公式

/ Y- u! O" w+ N

h(X,t)=\sum_{K}^{}{x}\tilde{h}(K,t)e^{iK\cdot X}

" ?/ t) J+ s1 G

其中 K 代表波正在运动的二维水平方向, K = (k_{x}, K_{y})

5 m% @2 @* A( B& n

k_{x} = 2\pi/L_{x} , k_{y} = 2\pi m/L_{y}

: y+ e, j3 T, v U

n 是水平方向的Domain Resolution。 m 是竖直方向的Domain Resolution

$ H3 ~4 F* W7 t1 N

所以 K 的范围是:

: a, ]# m0 e8 d3 Z' _

K=(2\pi/L_{x},2\pi m/L_{y})

7 i& y7 i5 u% I3 L. h. T+ `9 E

The fft process generates the height field at discrete points

3 }9 x" H- l* e \' f5 A2 \, }6 j

X = (nL_{x}/N, mL_{z}/N)

! T! X4 j& V% ~" u$ P8 R! o# D4 u1 F
, a3 w' Y) Z4 J

按照上面的公式可以渲染得到下面一张图:

2 o: i; P. p, c; k; g* w5 c
8 R4 f5 x) A1 G$ W& R$ Z

(2)然后需要用高斯分布生成一个Phillips spectrum

7 b8 c; J8 L' `) }6 D+ e* r
8 m, y, }& L: l- e
! @! ]' g9 p, U* c2 i3 i$ w

会得到如下效果

3 m6 `, g5 X& Z5 j# s8 y0 H& P. ]
% r1 n5 U6 X2 o4 W" G

把它和gauss分布结合后得到如下效果:

1 m. [+ e" q7 p/ v b3 G8 U9 o( O
9 v* _" l" g4 ?

然后现在我们有了海面的频谱图。下一步需要进行IFFT变换,但是在做变换前我们需要一个OmegaTexture做蝶形变换

% f! G# V2 f- k) G9 d- b# Q
8 _, k" ~) R+ o, u' p6 Y& r% j
- d; T$ U7 Y2 ~) B! k! s7 e0 t

(3)然后做IFFT变换生成Displacement,然后生成高度图和normal

7 F: |+ E3 l) S* P% a4 l$ u
9 N: i* i% ^$ l' u, E& g9 e

最后把这些生成的图弄到渲染管线种作为渲染资源渲染海面即可。下面就根据这上面的步骤来生成我们的海面。

; b- W v! `# _2 Q8 d3 F

下一卷我将给出具体FFT实现。

5 c2 J# S% m% V M0 a3 B

Enjoy it。

/ S, g! a; t' y( S; \2 B

Next:

/ Z6 J+ F6 `8 T& d9 [) i

【参考资料】

' A$ |: H8 y h2 \

【1】Ocean Shader with Gerstner Waves

7 p; ^$ S% ?8 K0 \: d! t* Y! O# ]

【2】音速键盘猫:Shader相册第6期 --- 实时水面模拟与渲染(一)

1 n! z( I' m4 F A. y

【3】https://labs.karmaninteractive.com/ocean-simulation-pt-1-introduction-df134a47150

8 K; ^( T2 U- X2 d9 z3 O' N! e

【4】Ocean simulation part one: using the discrete Fourier transform

8 L3 m% y; i: ~# x1 Z

【5】Ocean simulation part two: using the fast Fourier transform

4 G$ j0 r% M F% o0 W

【6】https://www.slideshare.net/Codemotion/an-introduction-to-realistic-ocean-rendering-through-fft-fabio-suriano-codemotion-rome-2017

& p6 v4 M6 ?6 c7 z9 l# K2 Z# V" x

【7】海洋模擬 FFT算法實現--基於GPU的基2快速傅里葉變換 2維FFT算法實現--基於GPU的基2快速二維傅里葉變換 【pbrt】使用openFrameworks調用pbrt

- N( I) Z% ~" F: f& f! Z0 [

【8】白霂凡:一小时学会快速傅里叶变换(Fast Fourier Transform)

F6 a3 V" z& V3 A8 g& L

【9】海面模拟以及渲染(计算着色器、FFT、Reflection Matrix)

4 q3 A8 T+ V$ C* i( }9 Q: R3 q

【10】wubugui/Jerry-Tessendorf-2004

. N' J& ~8 n& }2 ]7 }! Q0 u" D0 g

【11】http://evasion.imag.fr/~Fabrice.Neyret/images/fluids-nuages/waves/Jonathan/articlesCG/waterslides2001.pdf

7 w6 \; L# `& H: s

【12】https://dsqiu.iteye.com/blog/1636299

' N- M1 I/ v9 O- _

【13】https://www.bilibili.com/video/av19141078?t=1053

W% \) ]0 A# R! d' L

傅里叶变换基础教程(如果对傅里叶变换完全不清楚的建议按顺序看完下面的链接)

( w" X( e. g, l3 t" D2 [" I7 n

【14】Heinrich:傅里叶分析之掐死教程(完整版)更新于2014.06.06

6 W2 D) E: Q4 g- O; T2 U

【15】https://www.bilibili.com/video/av34364399/?spm_id_from=333.788.videocard.0

, J2 J: k( D% o

【16】https://www.bilibili.com/video/av34556069/?spm_id_from=333.788.videocard.0

* w1 M$ e0 f$ |

【17】https://www.bilibili.com/video/av34845617/?spm_id_from=333.788.videocard.0

* l [: H1 B6 M7 d( K# T5 [# Y6 A+ J

【18】https://www.bilibili.com/video/av35047004/?spm_id_from=333.788.videocard.0

7 O3 U) [' b6 P. I8 ?: ` }

【19】https://www.bilibili.com/video/av35810587/?spm_id_from=333.788.videocard.0

$ ~& v! y# {) B0 @+ B( M1 \

【20】https://www.bilibili.com/video/av36343956/?spm_id_from=333.788.videocard.0

' K, x. p% x* [9 X ; H: u9 {6 s9 T8 f! @ 4 Q- p' W. W" [2 E) R9 ^6 b; C 7 u* S" Y3 G# m 3 z: r+ L' ~+ x9 j4 Q
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
深宫鲍
活跃在2026-4-6
快速回复 返回顶部 返回列表