Power spectrum of field over an arbitrarily-shaped country
The power spectrum of a scalar field gives an indication of the spatial autocorrelation of the field, i.e. the typical size of a patch of similar values.
The power spectrum of a square or rectangular field or matrix can be determined using the Fourier transform or Fourier series.
- How can the power spectrum be determined for a country which is non-rectangular?
For example, the image below of the USA can be encoded into a matrix where all cells covering ocean are filled with NaN, while all cells covering land will have values between 25 and 375.
- What is a good strategy to deal with such large patches of NaNs?
- Can the country be mapped onto a rectangle and still provide a good measure for the original's power spectrum?
- Could spherical harmonics help, by noting the USA is a subsection of a sphere?
- Could a Wavelet transform be applied, and how?
spectral-theory fourier-transform correlation wavelets
add a comment |
The power spectrum of a scalar field gives an indication of the spatial autocorrelation of the field, i.e. the typical size of a patch of similar values.
The power spectrum of a square or rectangular field or matrix can be determined using the Fourier transform or Fourier series.
- How can the power spectrum be determined for a country which is non-rectangular?
For example, the image below of the USA can be encoded into a matrix where all cells covering ocean are filled with NaN, while all cells covering land will have values between 25 and 375.
- What is a good strategy to deal with such large patches of NaNs?
- Can the country be mapped onto a rectangle and still provide a good measure for the original's power spectrum?
- Could spherical harmonics help, by noting the USA is a subsection of a sphere?
- Could a Wavelet transform be applied, and how?
spectral-theory fourier-transform correlation wavelets
add a comment |
The power spectrum of a scalar field gives an indication of the spatial autocorrelation of the field, i.e. the typical size of a patch of similar values.
The power spectrum of a square or rectangular field or matrix can be determined using the Fourier transform or Fourier series.
- How can the power spectrum be determined for a country which is non-rectangular?
For example, the image below of the USA can be encoded into a matrix where all cells covering ocean are filled with NaN, while all cells covering land will have values between 25 and 375.
- What is a good strategy to deal with such large patches of NaNs?
- Can the country be mapped onto a rectangle and still provide a good measure for the original's power spectrum?
- Could spherical harmonics help, by noting the USA is a subsection of a sphere?
- Could a Wavelet transform be applied, and how?
spectral-theory fourier-transform correlation wavelets
The power spectrum of a scalar field gives an indication of the spatial autocorrelation of the field, i.e. the typical size of a patch of similar values.
The power spectrum of a square or rectangular field or matrix can be determined using the Fourier transform or Fourier series.
- How can the power spectrum be determined for a country which is non-rectangular?
For example, the image below of the USA can be encoded into a matrix where all cells covering ocean are filled with NaN, while all cells covering land will have values between 25 and 375.
- What is a good strategy to deal with such large patches of NaNs?
- Can the country be mapped onto a rectangle and still provide a good measure for the original's power spectrum?
- Could spherical harmonics help, by noting the USA is a subsection of a sphere?
- Could a Wavelet transform be applied, and how?
spectral-theory fourier-transform correlation wavelets
spectral-theory fourier-transform correlation wavelets
edited Jan 27 '18 at 18:06
Hans
4,96821032
4,96821032
asked Jan 18 '18 at 12:38
LBogaardtLBogaardt
4316
4316
add a comment |
add a comment |
1 Answer
1
active
oldest
votes
For a given function $L^2$ function $f$ on the domain $Omega$ (your irregular country), you can define the spectral decomposition as the expansion in the basis of the set of eigenfunctions with Dirichlet boundary condition $left{(k, phi_k);Big|nabla^2 phi_k=-kphi_k, phi_kbig|_{partialOmega}=begin{cases} fbig|_{partialOmega}, &k=0 \0, &k>0end{cases};right}$. The eigenvalue equation is called the Helmholtz equation and its null boundary value problem called the Dirichlet eigenvalue problem. The set of eigenvalue function is $L^2$ norm complete, and thus forming a basis for the $L^2$ function space, and orthogonal which is guaranteed by the spectral theorem of compact self-adjoint operators. Particularly, for $k=0$, this problem is the Laplace's equation Dirichlet boundary value problem with a given function value on the boundary $fbig|_{partialOmega}$. It can be solved directly by the Green's function method through the integration of the double layer potential. So you can expand any $L^2$ function in this basis just like in the Fourier basis. Fourier series is but the same problem with $Omega$ being a rectangle.
The beauty of this method is precisely that the first $phi_0$ is essentially the average of the boundary and thus all the local minima and maxima are on the boundary, while the remaining basis $(phi_k)_{k>0}$ is only boundary geometry (shape) dependent. Incidentally, hearing the shape of a drum is an attempt at the inverse problem from the eigenvalue to the geometry of the boundary.
There is essentially the whole theory of elliptical operators and much of functional analysis behind the above solution of this problem. One has all the existence and uniqueness properties for the solution of the above eigenvalue PDE problems for the smooth boundary and boundary value. One may explore the solution properties for say fractal boundary and nowhere differentiable boundary values.
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2610505%2fpower-spectrum-of-field-over-an-arbitrarily-shaped-country%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
For a given function $L^2$ function $f$ on the domain $Omega$ (your irregular country), you can define the spectral decomposition as the expansion in the basis of the set of eigenfunctions with Dirichlet boundary condition $left{(k, phi_k);Big|nabla^2 phi_k=-kphi_k, phi_kbig|_{partialOmega}=begin{cases} fbig|_{partialOmega}, &k=0 \0, &k>0end{cases};right}$. The eigenvalue equation is called the Helmholtz equation and its null boundary value problem called the Dirichlet eigenvalue problem. The set of eigenvalue function is $L^2$ norm complete, and thus forming a basis for the $L^2$ function space, and orthogonal which is guaranteed by the spectral theorem of compact self-adjoint operators. Particularly, for $k=0$, this problem is the Laplace's equation Dirichlet boundary value problem with a given function value on the boundary $fbig|_{partialOmega}$. It can be solved directly by the Green's function method through the integration of the double layer potential. So you can expand any $L^2$ function in this basis just like in the Fourier basis. Fourier series is but the same problem with $Omega$ being a rectangle.
The beauty of this method is precisely that the first $phi_0$ is essentially the average of the boundary and thus all the local minima and maxima are on the boundary, while the remaining basis $(phi_k)_{k>0}$ is only boundary geometry (shape) dependent. Incidentally, hearing the shape of a drum is an attempt at the inverse problem from the eigenvalue to the geometry of the boundary.
There is essentially the whole theory of elliptical operators and much of functional analysis behind the above solution of this problem. One has all the existence and uniqueness properties for the solution of the above eigenvalue PDE problems for the smooth boundary and boundary value. One may explore the solution properties for say fractal boundary and nowhere differentiable boundary values.
add a comment |
For a given function $L^2$ function $f$ on the domain $Omega$ (your irregular country), you can define the spectral decomposition as the expansion in the basis of the set of eigenfunctions with Dirichlet boundary condition $left{(k, phi_k);Big|nabla^2 phi_k=-kphi_k, phi_kbig|_{partialOmega}=begin{cases} fbig|_{partialOmega}, &k=0 \0, &k>0end{cases};right}$. The eigenvalue equation is called the Helmholtz equation and its null boundary value problem called the Dirichlet eigenvalue problem. The set of eigenvalue function is $L^2$ norm complete, and thus forming a basis for the $L^2$ function space, and orthogonal which is guaranteed by the spectral theorem of compact self-adjoint operators. Particularly, for $k=0$, this problem is the Laplace's equation Dirichlet boundary value problem with a given function value on the boundary $fbig|_{partialOmega}$. It can be solved directly by the Green's function method through the integration of the double layer potential. So you can expand any $L^2$ function in this basis just like in the Fourier basis. Fourier series is but the same problem with $Omega$ being a rectangle.
The beauty of this method is precisely that the first $phi_0$ is essentially the average of the boundary and thus all the local minima and maxima are on the boundary, while the remaining basis $(phi_k)_{k>0}$ is only boundary geometry (shape) dependent. Incidentally, hearing the shape of a drum is an attempt at the inverse problem from the eigenvalue to the geometry of the boundary.
There is essentially the whole theory of elliptical operators and much of functional analysis behind the above solution of this problem. One has all the existence and uniqueness properties for the solution of the above eigenvalue PDE problems for the smooth boundary and boundary value. One may explore the solution properties for say fractal boundary and nowhere differentiable boundary values.
add a comment |
For a given function $L^2$ function $f$ on the domain $Omega$ (your irregular country), you can define the spectral decomposition as the expansion in the basis of the set of eigenfunctions with Dirichlet boundary condition $left{(k, phi_k);Big|nabla^2 phi_k=-kphi_k, phi_kbig|_{partialOmega}=begin{cases} fbig|_{partialOmega}, &k=0 \0, &k>0end{cases};right}$. The eigenvalue equation is called the Helmholtz equation and its null boundary value problem called the Dirichlet eigenvalue problem. The set of eigenvalue function is $L^2$ norm complete, and thus forming a basis for the $L^2$ function space, and orthogonal which is guaranteed by the spectral theorem of compact self-adjoint operators. Particularly, for $k=0$, this problem is the Laplace's equation Dirichlet boundary value problem with a given function value on the boundary $fbig|_{partialOmega}$. It can be solved directly by the Green's function method through the integration of the double layer potential. So you can expand any $L^2$ function in this basis just like in the Fourier basis. Fourier series is but the same problem with $Omega$ being a rectangle.
The beauty of this method is precisely that the first $phi_0$ is essentially the average of the boundary and thus all the local minima and maxima are on the boundary, while the remaining basis $(phi_k)_{k>0}$ is only boundary geometry (shape) dependent. Incidentally, hearing the shape of a drum is an attempt at the inverse problem from the eigenvalue to the geometry of the boundary.
There is essentially the whole theory of elliptical operators and much of functional analysis behind the above solution of this problem. One has all the existence and uniqueness properties for the solution of the above eigenvalue PDE problems for the smooth boundary and boundary value. One may explore the solution properties for say fractal boundary and nowhere differentiable boundary values.
For a given function $L^2$ function $f$ on the domain $Omega$ (your irregular country), you can define the spectral decomposition as the expansion in the basis of the set of eigenfunctions with Dirichlet boundary condition $left{(k, phi_k);Big|nabla^2 phi_k=-kphi_k, phi_kbig|_{partialOmega}=begin{cases} fbig|_{partialOmega}, &k=0 \0, &k>0end{cases};right}$. The eigenvalue equation is called the Helmholtz equation and its null boundary value problem called the Dirichlet eigenvalue problem. The set of eigenvalue function is $L^2$ norm complete, and thus forming a basis for the $L^2$ function space, and orthogonal which is guaranteed by the spectral theorem of compact self-adjoint operators. Particularly, for $k=0$, this problem is the Laplace's equation Dirichlet boundary value problem with a given function value on the boundary $fbig|_{partialOmega}$. It can be solved directly by the Green's function method through the integration of the double layer potential. So you can expand any $L^2$ function in this basis just like in the Fourier basis. Fourier series is but the same problem with $Omega$ being a rectangle.
The beauty of this method is precisely that the first $phi_0$ is essentially the average of the boundary and thus all the local minima and maxima are on the boundary, while the remaining basis $(phi_k)_{k>0}$ is only boundary geometry (shape) dependent. Incidentally, hearing the shape of a drum is an attempt at the inverse problem from the eigenvalue to the geometry of the boundary.
There is essentially the whole theory of elliptical operators and much of functional analysis behind the above solution of this problem. One has all the existence and uniqueness properties for the solution of the above eigenvalue PDE problems for the smooth boundary and boundary value. One may explore the solution properties for say fractal boundary and nowhere differentiable boundary values.
edited Jan 4 at 2:11
answered Jan 27 '18 at 2:16
HansHans
4,96821032
4,96821032
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2610505%2fpower-spectrum-of-field-over-an-arbitrarily-shaped-country%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown