Ondas en el Gran Lago Salado
Usando ImageMesh, usted puede convertir segmentos de imagen en objetos de BoundaryMeshRegion. Estas regiones de malla le permiten utilizar funciones en otras áreas, tales como métodos de elementos finitos (FEM).
Para ilustrar esta conexión FEM, determine los modos principales modos de ondas superficiales en el Gran Lago Salado en Utah.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_95.png)
img = EntityValue[Entity["Lake", "GreatSaltLake::yw8cf"], "Image"]
![](assets.es/waves-in-the-great-salt-lake/O_44.png)
Regularice la imagen con el filtro de media de desplazamiento.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_96.png)
img2 = MeanShiftFilter[img, 3, 0.1]
![](assets.es/waves-in-the-great-salt-lake/O_45.png)
Obtenga una segmentación por medio de una región en crecimiento.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_97.png)
mask = RegionBinarize[img2, \!\(\*
GraphicsBox[
TagBox[RasterBox[CompressedData["
1:eJzt1jEKwkAQQNFdK0uv4C1sLW0VD6AYxSZCFMRzCJ7XiF3SzFb7lf8ggUCK
D8MmM99f1sdJSuk67W/r3X3ZdbvHZtY/bNvr+dQ2h1V7a05Nt9h/Xnv21ytJ
kiRJkiRJkn5H7tVuGLMqzqo4q+KsiuNW8bKsimNWfbNwXVbF5Yzuqp0xZFUB
dBaui5kF/T5AsxJ1jNAsq0qgs3BlOSPDMrMLmvUBDcNnsbqov2x4FrWrdscA
NIu6SjCzoEOEZlGniF5Ta0eMmVXCqjjmCJlV5KzaDWNWxXGreFnMKvAIazeM
MUfIXEmtKsCs8hSWsCrOqjhmlSRJkiRJkvS/3tRrD1M=
"], {{0, 147}, {150,
0}}, {0, 1},
ColorFunction->GrayLevel],
BoxForm`ImageTag["Bit", ColorSpace -> Automatic, Interleaving -> None],
Selectable->False],
DefaultBaseStyle->"ImageGraphics",
ImageSizeRaw->{150, 147},
PlotRange->{{0, 150}, {0, 147}}]\), 1/5]
![](assets.es/waves-in-the-great-salt-lake/O_46.png)
Recupere un objeto de BoundaryMeshRegion de la superficie del lago.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_98.png)
\[ScriptCapitalR] = ImageMesh[mask]
![](assets.es/waves-in-the-great-salt-lake/O_47.png)
Cree una malla de la superficie del lago.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_99.png)
\[CapitalOmega] =
TriangulateMesh[\[ScriptCapitalR], MaxCellMeasure -> 8]
![](assets.es/waves-in-the-great-salt-lake/O_48.png)
Resuelva la ecuación de onda de la superficie del lago determinando las funciones propias de laplaciano dentro de la región del lago.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_100.png)
\[ScriptCapitalL] = -\!\(
\*SubsuperscriptBox[\(\[Del]\), \({x, y}\), \(2\)]\(\[CurlyPhi][x,
y]\)\);
Usando una condición de límite.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_101.png)
\[ScriptCapitalB] = DirichletCondition[\[CurlyPhi][x, y] == 0, True];
Genere una base ortonormal de funciones propias Φ con valores propios Λ.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_102.png)
{\[CapitalLambda], \[CapitalPhi]} =
NDEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]}, \[CurlyPhi][x,
y], {x, y} \[Element] \[CapitalOmega], 64];
Muestre los primeros seis modos de oscilaciones.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_103.png)
GraphicsGrid[
Partition[
Table[ContourPlot[\[CapitalPhi][[
k]], {x, y} \[Element] \[CapitalOmega], PlotRange -> All,
PlotLabel -> \[CapitalLambda][[k]], PlotTheme -> "Minimal"], {k,
6}],
3
],
ImageSize -> 512
]
![](assets.es/waves-in-the-great-salt-lake/O_49.png)
La evolución temporal de los modos de oscilación decadente es dada como sigue.
![Click for copyable input](assets.es/waves-in-the-great-salt-lake/In_104.png)
\[CapitalTheta][\[Lambda]_, \[Xi]_, t_] = FullSimplify[
DSolveValue[Join[{ \!\(
\*SubscriptBox[\(\[PartialD]\), \(t,
t\)]\(\(TraditionalForm\`\[CurlyTheta]\)[
t]\)\) + \[Xi] \!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(\(TraditionalForm\`\
\[CurlyTheta]\)[
t]\)\) == -\[Lambda] \!\(TraditionalForm\`\[CurlyTheta]\)[
t] }, {\!\(TraditionalForm\`\[CurlyTheta]\)[0] ==
1, \[CurlyTheta]'[0] ==
0} ], \!\(TraditionalForm\`\[CurlyTheta]\)[t], t],
{\[Lambda] > 0, \[Xi] > 0, \[Xi]^2 < 4 \[Lambda], t \[Element] Reals}
]
![](assets.es/waves-in-the-great-salt-lake/O_50.png)
Por medio de la expansión de una perturbación inicial en funciones propias y dejando que estas evolucionen en el tiempo, se obtiene una propagación de ondas en el lago.
![](assets.es/waves-in-the-great-salt-lake/swf_8.png)