-
Notifications
You must be signed in to change notification settings - Fork 11
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
How to leverage other data while gridding? (Enhancement) #96
Comments
So I assume your in situ data are near surface values also. Technically you can use satellite data as if they were in situ data by adding them into the observational array using (lon, lat, val) of each pixel. This will however give too much importance to the satellite data and will make cross-validation approaches more difficult (error correlations within the satellite data). To alleviate this and also address the relative errors between the two data types you can |
Many thanks for the reply!
|
I guess you can also use the satellite observations to create a background field: as you know the DIVA analysis are performed on anomalies with respect to a background or reference field, which is, in simple cases, a uniform field with a value equal to the average value of all the observations. Here you can do this:
Doing so, you ensure that the solution, in regions where no in situ obs. are available, take the value of the background field. And you don't need to sub-sample the satellite data. |
Many thanks for chiming in! Not sure I understood fully for #2 and #3. So you assume that when the anomalies are gridded, there won't be bulls eyes in regions where obs. data are not available? I doubt that is the case. I'm thinking maybe I should do this instead: |
Yes, I think that using a background field obtained by gridding the satellite observations, then performing the analysis of the in situ observations with that background field can help avoid the bulls eyes. I'm now checking the doc and the examples to see if there is an example of how to do it. Concerning the use of CPME, that can be a possibility, though I'm not sure what is the best way to extract the error fields at the locations of the observations. |
If your satellite data are already on the same grid as the analysis you want to do, just work with anomalies with respect to your satellite data. That way you automatically will have satellite data in regions where you do not have observations. Technically, to easily calculate anomalies you can first do a dummy analysis with your in situ data, recover the structure s from it and then use newobs=DIVAnd_residualobs(s, fisat); where fisat is the griddes satellite data (needs to be exactly on the same grid as your analysis) Then you do an analysis with the anomalies of newobs and at the end sum up that analysis with your fisat. Hope it makes sense |
As another thread (#90) points out, one of the major weaknesses of the DIVA gridding is the poor handling of regions with no data. Bullseyes are often created.
I wonder if you could consider enhancing the DIVAnd.jl, so that it can automatically leverage values from another grid when there are no data and the associated gridding errors are too large?
Here is an example user case:
Many thanks for your consideration!
The text was updated successfully, but these errors were encountered: