Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Variogram estimation from multiple fields #104

Closed
wants to merge 10 commits into from
Closed

Conversation

MuellerSeb
Copy link
Member

With this PR a new option for unstructured variogram estimation is added:
You can now estimate one variogram from multiple fields on the same point-set. This could be useful to estimate a single variogram to describe multiple realizations of the same process.

Example:

  • time series of precipitation data -> one variogram for all days needed

In addition, these fields can have missing values, denoted by a new no_data argument.

Checklist:

  • blackened
  • test written
  • examples added

@MuellerSeb MuellerSeb added the enhancement New feature or request label Sep 21, 2020
@MuellerSeb MuellerSeb added this to the 1.3 milestone Sep 21, 2020
@MuellerSeb MuellerSeb self-assigned this Sep 21, 2020
@MuellerSeb
Copy link
Member Author

@LSchueler ready for review! 🥰

Copy link
Member

@LSchueler LSchueler left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

That is actually quite nice syntactic sugar. I think I like it. Although it makes it even more complicated dealing with all the different field shapes internally. But that will hopefully one day change with version 2.0.

Shouldn't this be extended to the structured variogram estimation? - I imagine this could lead to confusions, when someone changes from unstructured to structured fields.

And some minor requests in the code comments.

Cheers!

fields = [field1, field2]

###############################################################################
# Now we estimate the variograms for both fields individual and then again
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

individual => individually


plt.plot(bin_center, gamma1, label="field 1")
plt.plot(bin_center, gamma2, label="field 2")
plt.plot(bin_center, gamma, label="field joint")
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

"field joint" => "joined fields" ?

gstools/variogram/estimator.pyx Show resolved Hide resolved
@MuellerSeb MuellerSeb changed the base branch from develop to direct-vario November 7, 2020 13:23
@MuellerSeb
Copy link
Member Author

This was added in the direct-vario branch with c044b62

Closing this PR in favor of the direct-vario branch.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants