Skip to content

Variogram estimation from multiple fields#104

Closed
MuellerSeb wants to merge 10 commits intodirect-variofrom
multi_vario
Closed

Variogram estimation from multiple fields#104
MuellerSeb wants to merge 10 commits intodirect-variofrom
multi_vario

Conversation

@MuellerSeb
Copy link
Copy Markdown
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
Copy Markdown
Member Author

@LSchueler ready for review! 🥰

Copy link
Copy Markdown
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
Copy Markdown
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
Copy Markdown
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" ?

@MuellerSeb MuellerSeb changed the base branch from develop to direct-vario November 7, 2020 13:23
@MuellerSeb
Copy link
Copy Markdown
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