AKTIVE

cumulation.tcl at trunk
Login

cumulation.tcl at trunk

Annotate Artifact Blame Check-ins Using Doc Download Hex

File etc/transformer/statistics/cumulation.tcl artifact 5330713149 on branch trunk

     1
     2
     3
     4
     5
     6
     7
     8
     9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96

## -*- mode: tcl ; fill-column: 90 -*-
# # ## ### ##### ######## ############# #####################
## Transformers -- Cumulations along the various axes
##                 (rows, columns, bands), and the entire image.
##                 Tiled does not make sense.
##
## Note: Cumulation of the entire image is also known as `Sum Area Table`.
##
##       Due to the wrap around from one row to the next the SAT is __not separable__ into
##       a combination of row and column csums. It might be possible to use a veccache for
##       some perf boost.

# # ## ### ##### ######## ############# #####################

operator okind {
    op::row::cumulative    column
    op::column::cumulative row
} {
    op -> _ kind _

    example {
	aktive op sdf 2image smooth [aktive op sdf ring [aktive image sdf triangle width 32 height 32 a {10 10} b {50 80} c {80 30}] thickness 4]
	@1 | -matrix
    }

    section transform statistics

    note Returns image with the input ${kind}s transformed into cumulative sums.

    note This means that each pixel in a ${kind} is the sum of the values in the \
	${okind} before it, having the same ${kind}.

    note The result has the same geometry as the input. Only the contents change.

    # !xref-mark cumulation
    cached $kind cumulation AKTIVE_CSUM_FILL
}

# # ## ### ##### ######## ############# #####################

operator op::band::cumulative {
    section transform statistics

    input

    note Returns image with the input bands transformed into cumulative sums.

    note This means that each pixel in a band is the sum of the values in the \
	bands before it, having the same row and column.

    note The result has the same geometry as the input. Only the contents change.

    state -setup {
	aktive_geometry_copy (domain, aktive_image_get_geometry (srcs->v[0]));
    }

    blit prefixsum {
	{DH {y 0 1 up} {y 0 1 up}}
	{DW {x 0 1 up} {x 0 1 up}}
    } {raw sum-bands {
	// dstvalue = row/col start - 1-strided full band vector
	// srcvalue = row/col start - 1-strided full band vector

	aktive_cumulative_sum (dstvalue, DD, srcvalue, 1);

	// Applies only when src can be assumed to contain only positive or negative values
	// If a mix is present the result can return to the initial value
	// ASSERT (dstvalue[DD-1] > dstvalue[0], "bad cumulation is flat");
    }}

    pixels {
	// No caching is required for the bands.
	// The request is passed unchanged.
	// We can and do use a regular blitter.

	aktive_rectangle_def_as (subrequest, request);
	TRACE_RECTANGLE_M("band csum", &subrequest);
	aktive_block* src = aktive_region_fetch_area (0, &subrequest);

	#define AH (dst->height)
	#define AW (dst->width)
	#define AX (dst->x)
	#define AY (dst->y)
	@@prefixsum@@
	#undef AH
	#undef AW
	#undef AX
	#undef AY

	TRACE_DO (__aktive_block_dump ("band csum out", block));
    }
}

##
# # ## ### ##### ######## ############# #####################
::return