duplicate an AreaDefinition object with coarser/finer definition

[pytroll]

duplicate an AreaDefinition object with coarser/finer definition

[mraspaud]

Sounds good. How about you can also make the grid coarser ? like half the cells in x and y ?

[TomLav]

Originally, my request also had getting the coarser grid spacing. But mixing the two in a single function meant a lot of checks and documentation. Maybe have another method like:
new_area_def = area_def.coarser_grid_spacing(3,2) ?

[mraspaud]

what about the user provides the number of cells he/she/it wants to have, like:

new_area_def = area_def.new_grid_size(x_size / 2, y_size * 3)

?

[TomLav]

Could work. But the routine has to check the validity of the input (that they are exactly multiples or fractions of what is already there).

[TomLav]

And just to think ahead: I see two applications of this:

1. motion detection. For motion detection from cross-correlation techniques you need two grids: the target grid which is the grid locations where you will have your motion vectors on, and the imaging grid, which is much finer resolution, and is where you reproject your satellite images. The motion algorithms are much easier to write if the two grids are exactly in line, but the imaging grid has finer resolution.

2. anti-aliasing in pyresample: here the idea is to instantiate a finer grid than the target area_def, do the pyresample against the finer grid, and finally collect the resampled output from the finer grid to the original target grid (simple average).

[mraspaud]

Sounds good. I was just thinking of that this could be quite generic and allow any size, but I can see the utility of having the method allowing only integer ratios.

[mraspaud]

@TomLav what's the status on this ?

[TomLav]

I have a version of that in my modules somewhere, but never pushed it to the codebase.

[mraspaud]

Please do :)