Skip to content

Improve estimator variance and allocation formulas #32

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.

Sign up for GitHub

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

Open
XavierCLL wants to merge 1 commit into
base: main
Choose a base branch
from
Open

Improve estimator variance and allocation formulas #32

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
22 changes: 14 additions & 8 deletions AcATaMa/gui/accuracy_assessment_results.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,6 +35,18 @@ def rf(fv, r=5):
return round(fv, r)


def post_stratified_variance_table(error_matrix):
return [
[
cell_count * (1 - cell_count / total_row) / (total_row - 1)
if total_row not in [0, 1]
else np.nan
for cell_count in row
]
for total_row, row in zip([sum(r) for r in error_matrix], error_matrix, strict=True)
]


@error_handler
def get_html(accu_asse):
###########################################################################
Expand Down Expand Up @@ -637,10 +649,7 @@ def get_html(accu_asse):

# the error for post-stratified
if accu_asse.estimator == "Simple/systematic post-stratified estimator":
ui_table = [
[(1 - i / total_row) ** 2 * i / (total_row - 1) if total_row not in [0, 1] else np.nan for i in row]
for total_row, row in zip([sum(r) for r in accu_asse.error_matrix], accu_asse.error_matrix, strict=True)
]
ui_table = post_stratified_variance_table(accu_asse.error_matrix)
wi = [accu_asse.thematic_pixels_count[value] / total_pixels_classes for value in accu_asse.values]
variance = [
((1 - total_samples / total_pixels_classes) / total_samples)
Expand Down Expand Up @@ -1038,10 +1047,7 @@ def export_to_csv(accu_asse, file_out, csv_separator, csv_decimal_separator):

# the error for post-stratified
if accu_asse.estimator == "Simple/systematic post-stratified estimator":
ui_table = [
[(1 - i / total_row) ** 2 * i / (total_row - 1) if total_row not in [0, 1] else np.nan for i in row]
for total_row, row in zip([sum(r) for r in accu_asse.error_matrix], accu_asse.error_matrix, strict=True)
]
ui_table = post_stratified_variance_table(accu_asse.error_matrix)
wi = [accu_asse.thematic_pixels_count[value] / total_pixels_classes for value in accu_asse.values]
variance = [
((1 - total_samples / total_pixels_classes) / total_samples)
Expand Down
6 changes: 4 additions & 2 deletions AcATaMa/utils/sampling_utils.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,13 +49,15 @@ def get_num_samples_by_area_based_proportion(srs_table, total_std_error):
total_pixel_count = float(sum(mask(srs_table["pixel_count"], srs_table["On"])))
ratio_pixel_count = [p_c / total_pixel_count for p_c in mask(srs_table["pixel_count"], srs_table["On"])]
Si = [(float(ui) * (1 - float(ui))) ** 0.5 for ui in mask(srs_table["ui"], srs_table["On"])]
total_num_samples = (sum(rpc * si for rpc, si in zip(ratio_pixel_count, Si, strict=True)) / total_std_error) ** 2
weighted_std_sum = sum(rpc * si for rpc, si in zip(ratio_pixel_count, Si, strict=True))
total_num_samples = (weighted_std_sum / total_std_error) ** 2

num_samples = []
idx = 0
for item_enable in srs_table["On"]:
if item_enable:
num_samples.append(str(round(ratio_pixel_count[idx] * total_num_samples)))
allocation_weight = ratio_pixel_count[idx] * Si[idx] / weighted_std_sum if weighted_std_sum > 0 else 0
num_samples.append(str(round(allocation_weight * total_num_samples)))
idx += 1
else:
num_samples.append("0")
Expand Down
Loading