[GRASS-SVN] r57600 - sandbox/turek/scatter_plot
svn_grass at osgeo.org
svn_grass at osgeo.org
Thu Sep 5 08:30:54 PDT 2013
Author: turek
Date: 2013-09-05 08:30:53 -0700 (Thu, 05 Sep 2013)
New Revision: 57600
Modified:
sandbox/turek/scatter_plot/README.txt
sandbox/turek/scatter_plot/testing_patch.diff
Log:
scatter plot: bugs fixing
Modified: sandbox/turek/scatter_plot/README.txt
===================================================================
--- sandbox/turek/scatter_plot/README.txt 2013-09-05 14:16:48 UTC (rev 57599)
+++ sandbox/turek/scatter_plot/README.txt 2013-09-05 15:30:53 UTC (rev 57600)
@@ -17,27 +17,28 @@
TODO suggestions:
Scatterplot:
-- don't allow to go on when no subgroup is present:
- - error message about missing subgroup in group <%s>
- - fall back to group selection dialog + error message about missing
+- don't allow to go on when no subgroup is present:
+ - error message about missing subgroup in group <%s> #DONE
+ - fall back to group selection dialog + error message about missing #DONE
-- Dialog "Add scatter plot"
- Band 1 (x axis) <<- add
- Band 2 (y axis) <<- add
+- Dialog "Add scatter plot"
+ Band 1 (x axis) <<- add #DONE
+ Band 2 (y axis) <<- add #DONE
- Scatter plot toolbar:
- - put "Add SP" icon to the left
- - when in "Select area" mode, un-highlight zoom + pan icons
+ - put "Add SP" icon to the left #DONE
+ - when in "Select area" mode, un-highlight zoom + pan icons
+ - zoom to extend in toolbar
- Scatter plot window
- - zoom icon: also allow for zoom box
- - mouse wheel: 0.1 zoom step (not 0.5)
- - add values to axis
- - plot legend: strip @mapset off
+ - zoom icon: also allow for zoom box #DONE
+ - mouse wheel: 0.1 zoom step (not 0.5)#DONE
+ - add values to axis # it is possible but IDEA show values in statusbar
+ - plot legend: strip @mapset off #DONE
- second plot is falling out of the main window to the left
- Add class
- - rerendering zooms out, it should not
+ - rerendering zooms out, it should not #DONE
- confidence ellipse
- connect to "xx std dev" selection in toolbar
Modified: sandbox/turek/scatter_plot/testing_patch.diff
===================================================================
--- sandbox/turek/scatter_plot/testing_patch.diff 2013-09-05 14:16:48 UTC (rev 57599)
+++ sandbox/turek/scatter_plot/testing_patch.diff 2013-09-05 15:30:53 UTC (rev 57600)
@@ -1,682 +1,1285 @@
-Index: include/imagery.h
+Index: lib/vector/Vlib/open.c
===================================================================
---- include/imagery.h (revision 57581)
-+++ include/imagery.h (working copy)
-@@ -135,6 +135,56 @@
-
- } IClass_statistics;
-
-+/* Scatter Plot backend */
+--- lib/vector/Vlib/open.c (revision 57599)
++++ lib/vector/Vlib/open.c (working copy)
+@@ -240,7 +240,9 @@
+ }
+ else {
+ char file_path[GPATH_MAX];
+-
++ /* reduce to current mapset if search path was set */
++ if(strcmp(Map->mapset, "") == 0)
++ Map->mapset = G_store(G_mapset());
+ /* temporary map: reduce to current mapset if search path
+ * was set */
+ if (strcmp(Map->mapset, "") == 0)
+Index: lib/imagery/scatt.c
+===================================================================
+--- lib/imagery/scatt.c (revision 0)
++++ lib/imagery/scatt.c (working copy)
+@@ -0,0 +1,746 @@
++/*!
++ \file lib/imagery/scatt.c
+
-+#define SC_SCATT_DATA 0
-+#define SC_SCATT_CONDITIONS 1
++ \brief Imagery library - functions for wx Scatter Plot Tool.
+
++ Low level functions used by wx Scatter Plot Tool.
+
-+/*! Holds list of all catagories.
-+ It can contain selected areas for scatter plots (SC_SCATT_CONDITIONS type) or computed scatter plots (SC_SCATT_DATA type).
-+*/
-+struct scCats
++ Copyright (C) 2013 by the GRASS Development Team
++
++ This program is free software under the GNU General Public License
++ (>=v2). Read the file COPYING that comes with GRASS for details.
++
++ \author Stepan Turek <stepan.turek at seznam.cz> (Mentor: Martin Landa)
++ */
++
++#include <grass/raster.h>
++#include <grass/imagery.h>
++#include <grass/gis.h>
++#include <grass/glocale.h>
++
++#include <stdio.h>
++#include <stdlib.h>
++#include <math.h>
++#include <string.h>
++
++
++struct rast_row
+{
-+ int type; /*!< SC_SCATT_DATA -> computed scatter plots, SC_SCATT_CONDITIONS -> set conditions for scatter plots to be computed*/
++ CELL * row;
++ char * null_row;
++ struct Range rast_range; /*Range of whole raster.*/
++};
+
-+ int n_cats; /*!< number of alocated categories */
++/*!
++ \brief Create pgm header.
++
++ Scatter plot internally generates pgm files. These pgms have header in format created by this function.
+
-+ int n_bands; /*!< number of analyzed bands */
-+ int n_scatts; /*!< number of possible scattter plot which can be created from bands */
++ \param region region to be pgm header generated for
++ \param [out] header header of pgm file
++ */
++static int get_cat_rast_header(struct Cell_head * region, char * header){
++ return sprintf(header, "P5\n%d\n%d\n1\n", region->cols, region->rows);
++}
+
-+ int n_a_cats; /*!< number of used/active categories */
-+ int * cats_ids; /*!< (cat_idx->cat_id) array index is internal idx (position in cats_arr) and id is saved in it's position*/
-+ int * cats_idxs; /*!< (cat_id->cat_idx) array index is id and internal idx is saved in it's position*/
++/*!
++ \brief Create category raster conditions file.
++
++ \param cat_rast_region region to be file generated for
++ \param cat_rast path of generated category raster file
++ */
++int I_create_cat_rast(struct Cell_head * cat_rast_region, const char * cat_rast)
++{
++ FILE * f_cat_rast;
++ char cat_rast_header[1024];//TODO magic number
++ int i_row, i_col;
++ int head_nchars;
+
-+ struct scScatts ** cats_arr; /*!< array of pointers to struct scScatts */
-+};
++ unsigned char * row_data;
+
++ f_cat_rast = fopen(cat_rast, "wb");
++ if(!f_cat_rast) {
++ G_warning("Unable to create category raster condition file <%s>.", cat_rast);
++ return -1;
++ }
+
-+/*! Holds list of all scatter plots, which belongs to category.
++ head_nchars = get_cat_rast_header(cat_rast_region, cat_rast_header);
++
++ fwrite(cat_rast_header, sizeof(char), head_nchars/sizeof(char), f_cat_rast);
++ if (ferror(f_cat_rast)){
++ fclose(f_cat_rast);
++ G_warning(_("Unable to write header into category raster condition file <%s>."), cat_rast);
++ return -1;
++ }
++
++ row_data = (unsigned char *) G_malloc(cat_rast_region->cols * sizeof(unsigned char));
++ for(i_col = 0; i_col < cat_rast_region->cols; i_col++)
++ row_data[i_col] = 0 & 255;
++
++ for(i_row = 0; i_row < cat_rast_region->rows; i_row++) {
++ fwrite(row_data, sizeof(unsigned char), (cat_rast_region->cols)/sizeof(unsigned char), f_cat_rast);
++ if (ferror(f_cat_rast))
++ {
++ fclose(f_cat_rast);
++ G_warning(_("Unable to write into category raster condition file <%s>."), cat_rast);
++ return -1;
++ }
++ }
++
++ fclose(f_cat_rast);
++ return 0;
++}
++
++static int print_reg(struct Cell_head * intersec, const char * pref, int dbg_level)
++{
++ G_debug(dbg_level, "%s:\n n:%f\ns:%f\ne:%f\nw:%f\nns_res:%f\new_res:%f", pref, intersec->north, intersec->south,
++ intersec->east, intersec->west, intersec->ns_res, intersec->ew_res);
++}
++
++/*!
++ \brief Find intersection region of two regions.
++
++ \param A pointer to intersected region
++ \param B pointer to intersected region
++ \param [out] intersec pointer to intersection region of regions A B (relevant params of the region are: south, north, east, west)
++
++ \return 0 if interection exists
++ \return -1 if regions does not intersect
++ */
++static int regions_intersecion(struct Cell_head * A, struct Cell_head * B, struct Cell_head * intersec)
++{
++
++ if(B->north < A->south) return -1;
++ else if(B->north > A->north) intersec->north = A->north;
++ else intersec->north = B->north;
++
++ if(B->south > A->north) return -1;
++ else if(B->south < A->south) intersec->south = A->south;
++ else intersec->south = B->south;
++
++ if(B->east < A->west) return -1;
++ else if(B->east > A->east) intersec->east = A->east;
++ else intersec->east = B->east;
++
++ if(B->west > A->east) return -1;
++ else if(B->west < A->west) intersec->west = A->west;
++ else intersec->west = B->west;
++
++ if(intersec->north == intersec->south) return -1;
++
++ if(intersec->east == intersec->west) return -1;
++
++ return 0;
++
++}
++
++/*!
++ \brief Get rows and cols numbers, which defines intersection of the regions.
++
++ \param A pointer to intersected region
++ \param B pointer to intersected region (A and B must have same resolution)
++ \param [out] A_bounds rows and cols numbers of A stored in south, north, east, west, which defines intersection of A and B
++ \param [out] B_bounds rows and cols numbers of B stored in south, north, east, west, which defines intersection of A and B
++
++ \return 0 if interection exists
++ \return -1 if regions do not intersect
++ \return -2 resolution of regions is not same
+*/
-+struct scScatts
++static int get_rows_and_cols_bounds(struct Cell_head * A, struct Cell_head * B, struct Cell_head * A_bounds, struct Cell_head * B_bounds)
+{
-+ int n_a_scatts; /*!< number of used/active scatter plots*/
++ float ns_res, ew_res;
++
++ struct Cell_head intersec;
++
++ /* TODO is it right check? */
++ if(abs(A->ns_res - B->ns_res) > GRASS_EPSILON) {
++ G_debug(0, "'get_rows_and_cols_bounds' ns_res does not fit, A->ns_res: %f B->ns_res: %f", A->ns_res, B->ns_res);
++ return -2;
++ }
++
++ if(abs(A->ew_res - B->ew_res) > GRASS_EPSILON) {
++ G_debug(0, "'get_rows_and_cols_bounds' ew_res does not fit, A->ew_res: %f B->ew_res: %f", A->ew_res, B->ew_res);
++ return -2;
++ }
++
++ ns_res = A->ns_res;
++ ew_res = A->ew_res;
++
++ if(regions_intersecion(A, B, &intersec) == -1)
++ return -1;
++
++ A_bounds->north = ceil((A->north - intersec.north - ns_res * 0.5) / ns_res);
++ A_bounds->south = ceil((A->north - intersec.south - ns_res * 0.5) / ns_res);
++
++ A_bounds->east = ceil((intersec.east - A->west - ew_res * 0.5) / ew_res);
++ A_bounds->west = ceil((intersec.west - A->west - ew_res * 0.5) / ew_res);
++
++ B_bounds->north = ceil((B->north - intersec.north - ns_res * 0.5) / ns_res);
++ B_bounds->south = ceil((B->north - intersec.south - ns_res * 0.5) / ns_res);
++
++ B_bounds->east = ceil((intersec.east - B->west - ew_res * 0.5) / ew_res);
++ B_bounds->west = ceil((intersec.west - B->west - ew_res * 0.5) / ew_res);
++
++ return 0;
++}
++
++/*!
++ \brief Insert raster map patch into pgm file.
++
++ \param patch_rast name of raster map
++ \param cat_rast_region region of category raster file
++ \param cat_rast path to category raster file
++
++ \return 0 on success
++ \return -1 on failure
++*/
++int I_insert_patch_to_cat_rast(const char * patch_rast, struct Cell_head * cat_rast_region, const char * cat_rast)
++{
++
++ FILE * f_cat_rast;
++ struct Cell_head patch_region, patch_bounds, cat_rast_bounds;
++ char cat_rast_header[1024];//TODO magic number
++ int i_row, i_col, ncols, nrows, cat_rast_col, patch_col, val;
++ int head_nchars, ret;
++ int fd_patch_rast, init_shift, step_shift;
++ unsigned char * patch_data;
++
++ char * null_chunk_row;
++
++ const char *mapset;
++
++ struct Cell_head patch_lines, cat_rast_lines;
++
++ unsigned char * row_data;
++
++ f_cat_rast = fopen(cat_rast, "rb+");
++ if(!f_cat_rast) {
++ G_warning(_("Unable to open category raster condtions file <%s>."), cat_rast);
++ return -1;
++ }
++
++ head_nchars = get_cat_rast_header(cat_rast_region, cat_rast_header);
++ if ((mapset = G_find_raster(patch_rast,"")) == NULL) {
++ fclose(f_cat_rast);
++ G_warning(_("Unable to find patch raster <%s>."), patch_rast);
++ return -1;
++ }
++
++ Rast_get_cellhd(patch_rast, mapset, &patch_region);
++ Rast_set_window(&patch_region);
++
++ if ((fd_patch_rast = Rast_open_old(patch_rast, mapset)) < 0) {
++ fclose(f_cat_rast);
++ return -1;
++ }
++
++ ret = get_rows_and_cols_bounds(cat_rast_region, &patch_region, &cat_rast_bounds, &patch_bounds);
++ if(ret == -2) {
++ G_warning(_("Resolutions of patch <%s> and patched file <%s> are not same."), patch_rast, cat_rast);
++
++ Rast_close(fd_patch_rast);
++ fclose(f_cat_rast);
++
++ return -1;
++ }
++ else if (ret == -1){
++
++ Rast_close(fd_patch_rast);
++ fclose(f_cat_rast);
++
++ return 0;
++ }
++
++ ncols = cat_rast_bounds.east - cat_rast_bounds.west;
++ nrows = cat_rast_bounds.south - cat_rast_bounds.north;
++
++ patch_data = (unsigned char *) G_malloc(ncols * sizeof(unsigned char));
++
++ init_shift = head_nchars + cat_rast_region->cols * cat_rast_bounds.north + cat_rast_bounds.west;
++
++ if(fseek(f_cat_rast, init_shift, SEEK_SET) != 0) {
++ G_warning(_("Corrupted category raster conditions file <%s> (fseek failed)"), cat_rast);
++
++ Rast_close(fd_patch_rast);
++ G_free(null_chunk_row);
++ fclose(f_cat_rast);
++
++ return -1;
++ }
++
++ step_shift = cat_rast_region->cols - ncols;
++
++ null_chunk_row = Rast_allocate_null_buf();
+
-+ int * scatts_bands; /*!< array of bands for which represents the scatter plots,
-+ every scatter plot has assigned two bads (size of array is n_a_scatts * 2)*/
-+ int * scatt_idxs; /*!< (scatt_id->scatt_idx) internal idx of the scatter plot (position in scatts_arr)*/
++ for(i_row = 0; i_row < nrows; i_row++) {
++ Rast_get_null_value_row (fd_patch_rast, null_chunk_row, i_row + patch_bounds.north);
+
-+ struct scdScattData ** scatts_arr; /*!< array of pointers to scdScattData */
-+};
++ for(i_col = 0; i_col < ncols; i_col++) {
++ patch_col = patch_bounds.west + i_col;
+
-+/*! Holds scatter plot data.
++ if(null_chunk_row[patch_col] != 1)
++ patch_data[i_col] = 1 & 255;
++ else {
++ patch_data[i_col] = 0 & 255;
++ }
++ }
++
++ fwrite(patch_data, sizeof(unsigned char), (ncols)/sizeof(unsigned char), f_cat_rast);
++ if (ferror(f_cat_rast))
++ {
++ G_warning(_("Unable to write into category raster conditions file <%s>"), cat_rast);
++
++ Rast_close(fd_patch_rast);
++ G_free(null_chunk_row);
++ fclose(f_cat_rast);
++
++ return -1;
++ }
++ if(fseek(f_cat_rast, step_shift, SEEK_CUR) != 0) {
++ G_warning(_("Corrupted category raster conditions file <%s> (fseek failed)"), cat_rast);
++
++ Rast_close(fd_patch_rast);
++ G_free(null_chunk_row);
++ fclose(f_cat_rast);
++
++ return -1;
++ }
++ }
++
++ Rast_close(fd_patch_rast);
++ G_free(null_chunk_row);
++ fclose(f_cat_rast);
++ return 0;
++}
++
++/*!
++ \brief Updates scatter plots data in category by pixels which meets category conditions.
++
++ \param bands_rows data represents data describig one row from raster band
++ \param belongs_pix array which defines which pixels belongs to category (1 value) and which not (0 value)
++ \param [out] scatts pointer to scScatts struct of type SC_SCATT_DATA, which are modified according to values in belongs_pix
+*/
-+struct scdScattData
++static inline void update_cat_scatt_plts(struct rast_row * bands_rows, unsigned short * belongs_pix, struct scScatts * scatts)
+{
-+ int n_vals; /*!< Number of values in scatter plot. */
++ int band_axis_1, band_axis_2, i_scatt, array_idx, cat_idx, i_chunk_rows_pix, max_arr_idx;
+
-+ unsigned char * b_conds_arr; /*!< array of selected areas (used for SC_SCATT_CONDITIONS type) otherwise NULL */
-+ unsigned int * scatt_vals_arr; /*!< array of computed areas (used for SC_SCATT_DATA type) otherwise NULL */
-+};
++ CELL * b_1_row;
++ CELL * b_2_row;
++ char * b_1_null_row,* b_2_null_row;
++ struct rast_row b_1_rast_row, b_2_rast_row;
+
++ struct Range b_1_range, b_2_range;
++ int b_1_range_size;
+
- #define SIGNATURE_TYPE_MIXED 1
-
- #define GROUPFILE "CURGROUP"
-Index: include/defs/vedit.h
-===================================================================
---- include/defs/vedit.h (revision 57581)
-+++ include/defs/vedit.h (working copy)
-@@ -33,6 +33,8 @@
- int Vedit_merge_lines(struct Map_info *, struct ilist *);
-
- /* move.c */
-+int Vedit_move_areas(struct Map_info *, struct Map_info **, int,
-+ struct ilist *, double, double, double, int, double);
- int Vedit_move_lines(struct Map_info *, struct Map_info **, int,
- struct ilist *, double, double, double, int, double);
-
-Index: include/defs/imagery.h
-===================================================================
---- include/defs/imagery.h (revision 57581)
-+++ include/defs/imagery.h (working copy)
-@@ -110,6 +110,23 @@
- FILE *I_fopen_subgroup_ref_new(const char *, const char *);
- FILE *I_fopen_subgroup_ref_old(const char *, const char *);
-
-+/* scatt_plt.c */
-+void I_sc_init_cats(struct scCats *, int, int);
-+void I_sc_free_cats(struct scCats *);
-+int I_sc_add_cat(struct scCats *);
-+int I_sc_insert_scatt_data(struct scCats *, struct scdScattData *, int, int);
++ int row_size = Rast_window_cols();
+
-+void I_scd_init_scatt_data(struct scdScattData *, int, int, void *);
++ int * scatts_bands = scatts->scatts_bands;
+
-+int I_compute_scatts(struct Cell_head *, struct scCats *, const char **,
-+ const char **, int, struct scCats *, const char **);
++ for(i_scatt = 0; i_scatt < scatts->n_a_scatts; i_scatt++)
++ {
++ b_1_rast_row = bands_rows[scatts_bands[i_scatt * 2]];
++ b_2_rast_row = bands_rows[scatts_bands[i_scatt * 2 + 1]];
+
-+int I_create_cat_rast(struct Cell_head *, const char *);
-+int I_insert_patch_to_cat_rast(const char *, struct Cell_head *, const char *);
++ b_1_row = b_1_rast_row.row;
++ b_2_row = b_2_rast_row.row;
+
-+int I_id_scatt_to_bands(const int, const int, int *, int *);
-+int I_bands_to_id_scatt(const int, const int, const int, int *);
++ b_1_null_row = b_1_rast_row.null_row;
++ b_2_null_row = b_2_rast_row.null_row;
+
- /* sig.c */
- int I_init_signatures(struct Signature *, int);
- int I_new_signature(struct Signature *);
-Index: gui/wxpython/vdigit/wxdigit.py
-===================================================================
---- gui/wxpython/vdigit/wxdigit.py (revision 57581)
-+++ gui/wxpython/vdigit/wxdigit.py (working copy)
-@@ -17,7 +17,7 @@
- (and NumPy would be an excellent candidate for acceleration via
- e.g. OpenCL or CUDA; I'm surprised it hasn't happened already).
-
--(C) 2007-2011 by the GRASS Development Team
-+(C) 2007-2011, 2013 by the GRASS Development Team
-
- This program is free software under the GNU General Public License
- (>=v2). Read the file COPYING that comes with GRASS for details.
-@@ -27,6 +27,8 @@
-
- import grass.script.core as grass
-
-+from grass.pydispatch.signal import Signal
++ b_1_range = b_1_rast_row.rast_range;
++ b_2_range = b_2_rast_row.rast_range;
++
++ b_1_range_size = b_1_range.max - b_1_range.min + 1;
++ max_arr_idx = (b_1_range.max - b_1_range.min + 1) * (b_2_range.max - b_2_range.min + 1);
+
- from core.gcmd import GError
- from core.debug import Debug
- from core.settings import UserSettings
-@@ -176,7 +178,21 @@
-
- if self.poMapInfo:
- self.InitCats()
--
++ for(i_chunk_rows_pix = 0; i_chunk_rows_pix < row_size; i_chunk_rows_pix++)
++ {
++ /* pixel does not belongs to scatter plot or has null value in one of the bands */
++ if(!belongs_pix[i_chunk_rows_pix] ||
++ b_1_null_row[i_chunk_rows_pix] == 1 ||
++ b_2_null_row[i_chunk_rows_pix] == 1)
++ continue;
+
-+ self.emit_signals = False
++ /* index in scatter plot array */
++ array_idx = b_1_row[i_chunk_rows_pix] - b_1_range.min + (b_2_row[i_chunk_rows_pix] - b_2_range.min) * b_1_range_size;
+
-+ # signals which describes features changes during digitization,
-+ # activate them using EmitSignals method
-+ #TODO signal for errors?
-+ self.featureAdded = Signal('IVDigit.featureAdded')
-+ self.areasDeleted = Signal('IVDigit.areasDeleted')
-+ self.vertexMoved = Signal('IVDigit.vertexMoved')
-+ self.vertexAdded = Signal('IVDigit.vertexAdded')
-+ self.vertexRemoved = Signal('IVDigit.vertexRemoved')
-+ self.featuresDeleted = Signal('IVDigit.featuresDeleted')
-+ self.featuresMoved = Signal('IVDigit.featuresMoved')
-+ self.lineEdited = Signal('IVDigit.lineEdited')
++ if(array_idx < 0 || array_idx >= max_arr_idx) {
++ G_warning ("Data inconsistent. Value computed for scatter plot is out of initialized range.");
++ continue;
++ }
+
- def __del__(self):
- Debug.msg(1, "IVDigit.__del__()")
- Vect_destroy_line_struct(self.poPoints)
-@@ -188,7 +204,12 @@
- Vect_close(self.poBgMapInfo)
- self.poBgMapInfo = self.popoBgMapInfo = None
- del self.bgMapInfo
--
++ /* increment scatter plot value */
++ ++scatts->scatts_arr[i_scatt]->scatt_vals_arr[array_idx];
++ }
++ }
++}
++
++/*!
++ \brief Computes scatter plots data from chunk_rows.
++
++ \param scatts pointer to scScatts struct of type SC_SCATT_DATA, where are computed scatter plots stored
++ \param scatt_conds pointer to scScatts struct of type SC_SCATT_CONDITIONS, where are selected areas (condtitions)stored
++
++ \param chunk_rows data arrays of chunk_rows from analyzed raster bands (all data in chunk_rows, null_chunk_rows and belongs_pix arrays represents same region)
++ \param null_chunk_rows null data arrays of chunk_rows from analyzed raster bands
++ \param row_size size of data in chunk_rows, null_chunk_rows and belongs_pix arrays
++ \param f_cats_rasts_conds file which stores selected areas (conditions) from mapwindow see I_create_cat_rast() and I_pa
++ \param fd_cats_rasts array of openedraster maps which represents all selected pixels for category
++ \param region analysis region
++
++ \return 0 on success
++ \return -1 on failure
++*/
++static inline int compute_scatts_from_chunk_row(struct Cell_head *region, struct scCats * scatt_conds,
++ FILE ** f_cats_rasts_conds, struct rast_row * bands_rows,
++ struct scCats * scatts, int * fd_cats_rasts)
++{
++
++ int i_rows_pix, i_cat, i_scatt, n_a_scatts, n_pixs;
++ int cat_id, scatt_plts_cat_idx, array_idx, max_arr_idx;
++ char * b_1_null_row,* b_2_null_row;
++ struct rast_row b_1_rast_row, b_2_rast_row;
++ CELL * cat_rast_row;
++
++ struct scScatts * scatts_conds;
++ struct scScatts * scatts_scatt_plts;
++ struct scdScattData * conds;
++
++ struct Range b_1_range, b_2_range;
++ int b_1_range_size;
++
++ int * scatts_bands;
++ struct scdScattData ** scatts_arr;
++
++ CELL * b_1_row;
++ CELL * b_2_row;
++ unsigned char * i_scatt_conds;
++
++ int row_size = Rast_window_cols();
++
++ unsigned short * belongs_pix = (unsigned short *) G_malloc(row_size * sizeof(unsigned short));
++ unsigned char * rast_pixs = (unsigned char *) G_malloc(row_size * sizeof(unsigned char));
++ cat_rast_row = Rast_allocate_c_buf();
++
+
-+ def EmitSignals(self, emit):
-+ """!Activate/deactivate signals which describes features changes during digitization.
-+ """
-+ self.emit_signals = emit
++ for(i_cat = 0; i_cat < scatt_conds->n_a_cats; i_cat++)
++ {
++ scatts_conds = scatt_conds->cats_arr[i_cat];
+
- def CloseBackgroundMap(self):
- """!Close background vector map"""
- if not self.poBgMapInfo:
-@@ -394,7 +415,6 @@
-
- @return tuple (number of added features, feature ids)
- """
--
- layer = self._getNewFeaturesLayer()
- cat = self._getNewFeaturesCat()
-
-@@ -419,10 +439,14 @@
- return (-1, None)
-
- self.toolbar.EnableUndo()
--
-- return self._addFeature(vtype, points, layer, cat,
-- self._getSnapMode(), self._display.GetThreshold())
--
++ cat_id = scatt_conds->cats_ids[i_cat];
+
-+ ret = self._addFeature(vtype, points, layer, cat,
-+ self._getSnapMode(), self._display.GetThreshold())
-+ if ret[0] > -1 and self.emit_signals:
-+ self.featureAdded.emit(new_bboxs = [self._createBbox(points)], new_areas_cats = [[{layer : [cat]}, None]])
++ scatt_plts_cat_idx = scatts->cats_idxs[cat_id];
++ if(scatt_plts_cat_idx < 0)
++ continue;
++ scatts_scatt_plts = scatts->cats_arr[scatt_plts_cat_idx];
+
-+ return ret
++ G_zero(belongs_pix, row_size * sizeof(unsigned short));
+
- def DeleteSelectedLines(self):
- """!Delete selected features
-
-@@ -434,16 +458,27 @@
- # collect categories for deleting if requested
- deleteRec = UserSettings.Get(group = 'vdigit', key = 'delRecord', subkey = 'enabled')
- catDict = dict()
++ /* if category has no conditions defined, scatter plots without any constraint are computed (default scatter plots) */
++ if(!scatts_conds->n_a_scatts && !f_cats_rasts_conds[i_cat]) {
++ for(i_scatt = 0; i_scatt < scatts_scatt_plts->n_a_scatts; i_scatt++)
++ {
++ /* all pixels belongs */
++ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
++ belongs_pix[i_rows_pix] = 1;
++ }
++ }
++ /* compute belonging pixels for defined conditions */
++ else
++ {
++ scatts_bands = scatts_conds->scatts_bands;
+
-+ old_bboxs = []
-+ old_areas_cats = []
- if deleteRec:
- for i in self._display.selected['ids']:
-+
- if Vect_read_line(self.poMapInfo, None, self.poCats, i) < 0:
- self._error.ReadLine(i)
-
-- cats = self.poCats.contents
-- for j in range(cats.n_cats):
-- if cats.field[j] not in catDict.keys():
-- catDict[cats.field[j]] = list()
-- catDict[cats.field[j]].append(cats.cat[j])
-+ if self.emit_signals:
-+ ret = self._getLineAreaBboxCats(i)
-+ if ret:
-+ old_bboxs += ret[0]
-+ old_areas_cats += ret[1]
-+
-+ # catDict was not used -> put into comment
-+ #cats = self.poCats.contents
-+ #for j in range(cats.n_cats):
-+ # if cats.field[j] not in catDict.keys():
-+ # catDict[cats.field[j]] = list()
-+ # catDict[cats.field[j]].append(cats.cat[j])
-
- poList = self._display.GetSelectedIList()
- nlines = Vedit_delete_lines(self.poMapInfo, poList)
-@@ -456,7 +491,10 @@
- self._deleteRecords(catDict)
- self._addChangeset()
- self.toolbar.EnableUndo()
--
++ /* check conditions from category raster condtitions file */
++ if(f_cats_rasts_conds[i_cat]) {
++ n_pixs = fread(rast_pixs, sizeof(unsigned char), (row_size)/sizeof(unsigned char), f_cats_rasts_conds[i_cat]);
+
-+ if self.emit_signals:
-+ self.featuresDeleted.emit(old_bboxs = old_bboxs, old_areas_cats = old_areas_cats)
++ if (ferror(f_cats_rasts_conds[i_cat]))
++ {
++ G_free(rast_pixs);
++ G_free(belongs_pix);
++ G_warning(_("Unable to read from category raster condtition file."));
++ return -1;
++ }
++ if (n_pixs != n_pixs) {
++ G_free(rast_pixs);
++ G_free(belongs_pix);
++ G_warning(_("Invalid size of category raster conditions file."));
++ return -1;
+
- return nlines
-
- def _deleteRecords(self, cats):
-@@ -512,22 +550,173 @@
-
- @return number of deleted
- """
-+ if len(self._display.selected['ids']) < 1:
-+ return 0
-+
- poList = self._display.GetSelectedIList()
- cList = poList.contents
-
- nareas = 0
-+ old_bboxs = []
-+ old_areas_cats = []
++ }
+
- for i in range(cList.n_values):
++ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
++ {
++ if(rast_pixs[i_rows_pix] != 0 & 255)
++ belongs_pix[i_rows_pix] = 1;
++ }
++ }
+
- if Vect_get_line_type(self.poMapInfo, cList.value[i]) != GV_CENTROID:
- continue
--
++ /* check condtions defined in scatter plots*/
++ for(i_scatt = 0; i_scatt < scatts_conds->n_a_scatts; i_scatt++)
++ {
++ b_1_rast_row = bands_rows[scatts_bands[i_scatt * 2]];
++ b_2_rast_row = bands_rows[scatts_bands[i_scatt * 2 + 1]];
+
-+ if self.emit_signals:
-+ area = Vect_get_centroid_area(self.poMapInfo, cList.value[i]);
-+ if area > 0:
-+ bbox, cats = self._getaAreaBboxCats(area)
-+ old_bboxs += bbox
-+ old_areas_cats += cats
++ b_1_row = b_1_rast_row.row;
++ b_2_row = b_2_rast_row.row;
+
- nareas += Vedit_delete_area_centroid(self.poMapInfo, cList.value[i])
-
- if nareas > 0:
- self._addChangeset()
- self.toolbar.EnableUndo()
-+ if self.emit_signals:
-+ self.areasDeleted.emit(old_bboxs = old_bboxs, old_areas_cats = old_areas_cats)
++ b_1_null_row = b_1_rast_row.null_row;
++ b_2_null_row = b_2_rast_row.null_row;
+
-+ return nareas
++ b_1_range = b_1_rast_row.rast_range;
++ b_2_range = b_2_rast_row.rast_range;
++
++ b_1_range_size = b_1_range.max - b_1_range.min + 1;
++ max_arr_idx = (b_1_range.max - b_1_range.min + 1) * (b_2_range.max - b_2_range.min + 1);
++
++ i_scatt_conds = scatts_conds->scatts_arr[i_scatt]->b_conds_arr;
++
++ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
++ {
++ /* pixels already belongs to category from category raster conditions file or contains null value in one of the bands */
++ if(belongs_pix[i_rows_pix] ||
++ b_1_null_row[i_rows_pix] == 1 ||
++ b_2_null_row[i_rows_pix] == 1)
++ continue;
++
++ array_idx = b_1_row[i_rows_pix] - b_1_range.min + (b_2_row[i_rows_pix] - b_2_range.min) * b_1_range_size;
++ if(array_idx < 0 || array_idx >= max_arr_idx) {
++ G_warning ("Data inconsistent. Value computed for scatter plot is out of initialized range.");
++ continue;
++ }
++ /* pixels meets condtion defined in scatter plot */
++ if(i_scatt_conds[array_idx])
++ belongs_pix[i_rows_pix] = 1;
++ }
++ }
++ }
++
++ /* update category raster with belonging pixels */
++ if(fd_cats_rasts[i_cat] >= 0) {
++ Rast_set_null_value(cat_rast_row, Rast_window_cols(), CELL_TYPE);
++
++ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
++ if(belongs_pix[i_rows_pix])
++ cat_rast_row[i_rows_pix] = belongs_pix[i_rows_pix];
++
++ Rast_put_c_row (fd_cats_rasts[i_cat], cat_rast_row);
++ }
++
++ /* update scatter plots with belonging pixels */
++ update_cat_scatt_plts(bands_rows, belongs_pix, scatts_scatt_plts);
++ }
++
++ G_free(cat_rast_row);
++ G_free(rast_pixs);
++ G_free(belongs_pix);
++
++ return 0;
++}
++
++/*!
++ \brief Get list if bands needed to be opened for analysis from scCats struct.
++*/
++static void get_needed_bands(struct scCats * cats, int * b_needed_bands)
++{
++ // results in b_needed_bands - array of bools - if item has value 1, band (defined by item index) is needed to be opened
++ int i_cat, i_scatt, cat_id;
++
++ for(i_cat = 0; i_cat < cats->n_a_cats; i_cat++)
++ {
++ for(i_scatt = 0; i_scatt < cats->cats_arr[i_cat]->n_a_scatts; i_scatt++)
++ {
++ G_debug(3, "Active scatt %d in catt %d", i_scatt, i_cat);
++
++ b_needed_bands[cats->cats_arr[i_cat]->scatts_bands[i_scatt * 2]] = 1;
++ b_needed_bands[cats->cats_arr[i_cat]->scatts_bands[i_scatt * 2 + 1]] = 1;
++ }
++ }
++ return;
++}
++
++/*!
++ \brief Helper function for clean up.
++*/
++static void free_compute_scatts_data(int * fd_bands, struct rast_row * bands_rows, int n_a_bands, int * bands_ids,
++ int * fd_cats_rasts, FILE ** f_cats_rasts_conds, int n_a_cats)
++{
++ int i, band_id;
++
++ for(i = 0; i < n_a_bands; i++)
++ {
++ band_id = bands_ids[i];
++ if(band_id >= 0) {
++ Rast_close(fd_bands[i]);
++ G_free(bands_rows[band_id].row);
++ G_free(bands_rows[band_id].null_row);
++ }
++ }
++
++ if(f_cats_rasts_conds)
++ for(i = 0; i < n_a_cats; i++)
++ if(f_cats_rasts_conds[i])
++ fclose(f_cats_rasts_conds[i]);
++
++ if(fd_cats_rasts)
++ for(i = 0; i < n_a_cats; i++)
++ if(fd_cats_rasts[i] >= 0)
++ Rast_close(fd_cats_rasts[i]);
++
++}
++
++/*!
++ \brief Compute scatter plots data.
++
++ If category has not defined no category raster condition file and no scatter plot with consdtion,
++ default scatter plot is computed.
++
++ \param region analysis region, beaware that all input data must be prepared for this region (bands (their ranges), cats_rasts_conds rasters...)
++ \param region function calls Rast_set_window for this region
++ \param scatt_conds pointer to scScatts struct of type SC_SCATT_CONDITIONS, where are stored selected areas (conditions) in scatter plots
++ \param cats_rasts_conds paths to category raster conditions files representing selected areas (conditions) in rasters for every category
++ \param cats_rasts_conds index in array represents corresponding category id
++ \param cats_rasts_conds for manupulation with category raster conditions file see also I_id_scatt_to_bands() and I_insert_patch_to_cat_rast()
++ \param bands names of analyzed bands, order of bands is defined by their id
++ \param n_bands number of bands
++ \param [out] scatts pointer to scScatts struct of type SC_SCATT_DATA, where are computed scatter plots stored
++ \param [out] cats_rasts array of raster maps names where will be stored all selected pixels for every category
++
++ \return 0 on success
++ \return -1 on failure
++*/
++int I_compute_scatts(struct Cell_head *region, struct scCats * scatt_conds, const char ** cats_rasts_conds,
++ const char ** bands, int n_bands, struct scCats * scatts, const char ** cats_rasts)
++{
++ const char *mapset;
++ char header[1024];
++
++ int fd_cats_rasts[scatt_conds->n_a_cats];
++ FILE * f_cats_rasts_conds[scatt_conds->n_a_cats];
++
++ struct rast_row bands_rows[n_bands];
++
++ RASTER_MAP_TYPE data_type;
++
++ int nrows, i_band, n_a_bands, band_id,
++ i, i_row, head_nchars, i_cat, id_cat;
+
-+ def _getLineAreaBboxCats(self, ln_id):
-+ ltype = Vect_read_line(self.poMapInfo, None, None, ln_id)
++ int fd_bands[n_bands];
++ int bands_ids[n_bands];
++ int b_needed_bands[n_bands];
+
-+ if ltype == GV_CENTROID:
-+ #TODO centroid opttimization, can be adited also its area -> it will appear two times in new_ lists
-+ return self._getCentroidAreaBboxCats(ln_id)
-+ else:
-+ return [self._getBbox(ln_id)], [self._getLineAreasCategories(ln_id)]
++ Rast_set_window(region);
+
++ for(i_band = 0; i_band < n_bands; i_band++)
++ fd_bands[i_band] = -1;
++
++ for(i_band = 0; i_band < n_bands; i_band++)
++ bands_ids[i_band] = -1;
+
-+ def _getCentroidAreaBboxCats(self, centroid):
-+ if not Vect_line_alive(self.poMapInfo, centroid):
-+ return None
++ if (n_bands != scatts->n_bands ||
++ n_bands != scatt_conds->n_bands)
++ return -1;
+
-+ area = Vect_get_centroid_area(self.poMapInfo, centroid)
-+ if area > 0:
-+ return self._getaAreaBboxCats(area)
-+ else:
-+ return None
++ memset(b_needed_bands, 0, (size_t)n_bands * sizeof(int));
+
-+ def _getaAreaBboxCats(self, area):
++ get_needed_bands(scatt_conds, &b_needed_bands[0]);
++ get_needed_bands(scatts, &b_needed_bands[0]);
+
-+ po_b_list = Vect_new_list()
-+ Vect_get_area_boundaries(self.poMapInfo, area, po_b_list);
-+ b_list = po_b_list.contents
++ n_a_bands = 0;
+
-+ geoms = []
-+ areas_cats = []
++ /* open band rasters, which are needed for computation */
++ for(band_id = 0; band_id < n_bands; band_id++)
++ {
++ if(b_needed_bands[band_id])
++ {
++ G_debug(3, "Opening raster no. %d with name: %s", band_id, bands[band_id]);
+
-+ if b_list.n_values > 0:
-+ for i_line in range(b_list.n_values):
++ if ((mapset = G_find_raster2(bands[band_id],"")) == NULL) {
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
++ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
++ G_warning(_("Unbale to read find raster <%s>"), bands[band_id]);
++ return -1;
++ }
++
++ if ((fd_bands[n_a_bands] = Rast_open_old(bands[band_id], mapset)) < 0) {
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
++ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
++ G_warning(_("Unbale to open raster <%s>"), bands[band_id]);
++ return -1;
++ }
+
-+ line = b_list.value[i_line];
++ data_type = Rast_get_map_type(fd_bands[n_a_bands]);
++ if(data_type != CELL_TYPE) {
++ G_warning(_("Raster <%s> type is not <%s>"), bands[band_id], "CELL");
++ return -1;
++ }
+
-+ geoms.append(self._getBbox(abs(line)))
-+ areas_cats.append(self._getLineAreasCategories(abs(line)))
-
-- return nareas
-+ Vect_destroy_list(po_b_list);
++ bands_rows[band_id].row = Rast_allocate_c_buf();
++ bands_rows[band_id].null_row = Rast_allocate_null_buf();
+
-+ return geoms, areas_cats
++ if(Rast_read_range(bands[band_id], mapset, &bands_rows[band_id].rast_range) != 1){
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
++ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
++ G_warning(_("Unable to read range of raster <%s>"), bands[band_id]);
++ return -1;
++ }
+
-+ def _getLineAreasCategories(self, ln_id):
-+ if not Vect_line_alive (self.poMapInfo, ln_id):
-+ return []
++ bands_ids[n_a_bands] = band_id;
++ ++n_a_bands;
++ }
++ }
+
-+ ltype = Vect_read_line(self.poMapInfo, None, None, ln_id)
-+ if ltype != GV_BOUNDARY:
-+ return []
++ /* open category rasters condition files and category rasters */
++ for(i_cat = 0; i_cat < scatts->n_a_cats; i_cat++)
++ {
++ id_cat = scatts->cats_ids[i_cat];
++ if(cats_rasts[id_cat]) {
++ fd_cats_rasts[i_cat] = Rast_open_new(cats_rasts[id_cat], CELL_TYPE);
++ }
++ else
++ fd_cats_rasts[i_cat] = -1;
+
-+ cats = [None, None]
++ if(cats_rasts_conds[id_cat]) {
++ f_cats_rasts_conds[i_cat] = fopen(cats_rasts_conds[id_cat], "r");
++ if(!f_cats_rasts_conds[i_cat]) {
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids,
++ f_cats_rasts_conds, f_cats_rasts_conds, scatt_conds->n_a_cats);
++ G_warning(_("Unable to open category raster condtition file <%s>"), bands[band_id]);
++ return -1;
++ }
++ }
++ else
++ f_cats_rasts_conds[i_cat] = NULL;
++ }
+
-+ left = c_int()
-+ right = c_int()
++ head_nchars = get_cat_rast_header(region, header);
++ for(i_cat = 0; i_cat < scatt_conds->n_a_cats; i_cat++)
++ if(f_cats_rasts_conds[i_cat])
++ if( fseek(f_cats_rasts_conds[i_cat] , head_nchars, SEEK_SET) != 0) {
++ G_warning(_("Corrupted category raster conditions file (fseek failed)"));
++ return -1;
++ }
+
-+ if Vect_get_line_areas(self.poMapInfo, ln_id, pointer(left), pointer(right)) == 1:
-+ areas = [left.value, right.value]
++ nrows = Rast_window_rows();
+
-+ for i, a in enumerate(areas):
-+ if a > 0:
-+ centroid = Vect_get_area_centroid(self.poMapInfo, a)
-+ if centroid <= 0:
-+ continue
-+ c = self._getCategories(centroid)
-+ if c:
-+ cats[i] = c
++ /* analyze bands by rows */
++ for (i_row = 0; i_row < nrows; i_row++)
++ {
++ for(i_band = 0; i_band < n_a_bands; i_band++)
++ {
++ band_id = bands_ids[i_band];
++ Rast_get_c_row(fd_bands[i_band], bands_rows[band_id].row, i_row);
++ Rast_get_null_value_row (fd_bands[i_band], bands_rows[band_id].null_row, i_row);
++ }
++ if(compute_scatts_from_chunk_row(region, scatt_conds, f_cats_rasts_conds, bands_rows, scatts, fd_cats_rasts) == -1)
++ {
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids, fd_cats_rasts,
++ f_cats_rasts_conds, scatt_conds->n_a_cats);
++ return -1;
++ }
++
++ }
++ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids,
++ fd_cats_rasts, f_cats_rasts_conds, scatt_conds->n_a_cats);
++ return 0;
++}
+\ No newline at end of file
+Index: lib/imagery/scatt_sccats.c
+===================================================================
+--- lib/imagery/scatt_sccats.c (revision 0)
++++ lib/imagery/scatt_sccats.c (working copy)
+@@ -0,0 +1,405 @@
++/*!
++ \file lib/imagery/scatt_cat_rast.c
+
-+ return cats
++ \brief Imagery library - functions for manipulation with scatter plot structs.
+
-+ def _getCategories(self, ln_id):
-+ if not Vect_line_alive (self.poMapInfo, ln_id):
-+ return none
++ Copyright (C) 2013 by the GRASS Development Team
+
-+ poCats = Vect_new_cats_struct()
-+ if Vect_read_line(self.poMapInfo, None, poCats, ln_id) < 0:
-+ Vect_destroy_cats_struct(poCats)
-+ return None
++ This program is free software under the GNU General Public License
++ (>=v2). Read the file COPYING that comes with GRASS for details.
+
-+ cCats = poCats.contents
++ \author Stepan Turek <stepan.turek at seznam.cz> (Mentor: Martin Landa)
++ */
+
-+ cats = {}
-+ for j in range(cCats.n_cats):
-+ if cats.has_key(cCats.field[j]):
-+ cats[cCats.field[j]].append(cCats.cat[j])
-+ else:
-+ cats[cCats.field[j]] = [cCats.cat[j]]
-
-+ Vect_destroy_cats_struct(poCats)
-+ return cats
++#include <grass/raster.h>
++#include <grass/imagery.h>
++#include <grass/gis.h>
+
-+ def _getBbox(self, ln_id):
-+ if not Vect_line_alive (self.poMapInfo, ln_id):
-+ return None
++#include <stdio.h>
++#include <stdlib.h>
++#include <math.h>
++#include <string.h>
+
-+ poPoints = Vect_new_line_struct()
-+ if Vect_read_line(self.poMapInfo, poPoints, None, ln_id) < 0:
-+ Vect_destroy_line_struct(poPoints)
-+ return []
++/*!
++ \brief Compute band ids from scatter plot id.
+
-+ geom = self._convertGeom(poPoints)
-+ bbox = self._createBbox(geom)
-+ Vect_destroy_line_struct(poPoints)
-+ return bbox
++ Scatter plot id describes which bands defines the scatter plot.
+
-+ def _createBbox(self, points):
++ Let say we have 3 bands, their ids are 0, 1 and 2.
++ Scatter plot with id 0 consists of band 1 (b_1_id) 0 and band 2 (b_2_id) 1.
++ All scatter plots:
++ scatt_id b_1_id b_2_id
++ 0 0 1
++ 1 0 2
++ 2 1 2
+
-+ bbox = {}
-+ for pt in points:
-+ if not bbox.has_key('maxy'):
-+ bbox['maxy'] = pt[1]
-+ bbox['miny'] = pt[1]
-+ bbox['maxx'] = pt[0]
-+ bbox['minx'] = pt[0]
-+ continue
-+
-+ if bbox['maxy'] < pt[1]:
-+ bbox['maxy'] = pt[1]
-+ elif bbox['miny'] > pt[1]:
-+ bbox['miny'] = pt[1]
-+
-+ if bbox['maxx'] < pt[0]:
-+ bbox['maxx'] = pt[0]
-+ elif bbox['minx'] > pt[0]:
-+ bbox['minx'] = pt[0]
-+ return bbox
++ \param scatt_id scatter plot id
++ \param n_bands number of bands
++ \param [out] b_1_id id of band1
++ \param[out] b_2_id id of band2
+
-+ def _convertGeom(self, poPoints):
++ \return 0
++ */
++int I_id_scatt_to_bands(const int scatt_id, const int n_bands, int * b_1_id, int * b_2_id)
++{
++ int n_b1 = n_bands - 1;
+
-+ Points = poPoints.contents
++ * b_1_id = (int) ((2 * n_b1 + 1 - sqrt((double)((2 * n_b1 + 1) * (2 * n_b1 + 1) - 8 * scatt_id))) / 2);
+
-+ pts_geom = []
-+ for j in range(Points.n_points):
-+ pts_geom.append((Points.x[j], Points.y[j]))
++ * b_2_id = scatt_id - ((* b_1_id) * (2 * n_b1 + 1) - (* b_1_id) * (* b_1_id)) / 2 + (* b_1_id) + 1;
+
-+ return pts_geom
++ return 0;
++}
+
- def MoveSelectedLines(self, move):
- """!Move selected features
-
-@@ -536,16 +725,45 @@
- if not self._checkMap():
- return -1
-
-+ nsel = len(self._display.selected['ids'])
-+ if nsel < 1:
-+ return -1
-+
- thresh = self._display.GetThreshold()
- snap = self._getSnapMode()
-
- poList = self._display.GetSelectedIList()
+
-+ if self.emit_signals:
-+ old_bboxs = []
-+ old_areas_cats = []
-+ for sel_id in self._display.selected['ids']:
-+ ret = self._getLineAreaBboxCats(sel_id)
-+ if ret:
-+ old_bboxs += ret[0]
-+ old_areas_cats += ret[1]
-+
-+ Vect_set_updated(self.poMapInfo, 1)
-+ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
-+
- nlines = Vedit_move_lines(self.poMapInfo, self.popoBgMapInfo, int(self.poBgMapInfo is not None),
- poList,
- move[0], move[1], 0,
- snap, thresh)
++/*!
++ \brief Compute scatter plot id from band ids.
+
- Vect_destroy_list(poList)
--
++ See also I_id_scatt_to_bands().
+
-+ if nlines > 0 and self.emit_signals:
-+ new_bboxs = []
-+ new_areas_cats = []
-+ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
-+ for i in range(n_up_lines_old, n_up_lines):
-+ new_id = Vect_get_updated_line(self.poMapInfo, i)
-+ ret = self._getLineAreaBboxCats(new_id)
-+ if ret:
-+ new_bboxs += ret[0]
-+ new_areas_cats += ret[1]
++ \param n_bands number of bands
++ \param b_1_id id of band1
++ \param b_1_id id of band2
++ \param [out] scatt_id scatter plot id
+
- if nlines > 0 and self._settings['breakLines']:
- for i in range(1, nlines):
- self._breakLineAtIntersection(nlines + i, None, changeset)
-@@ -553,7 +771,13 @@
- if nlines > 0:
- self._addChangeset()
- self.toolbar.EnableUndo()
--
-+
-+ if self.emit_signals:
-+ self.featuresMoved.emit(new_bboxs = new_bboxs,
-+ old_bboxs = old_bboxs,
-+ old_areas_cats = old_areas_cats,
-+ new_areas_cats = new_areas_cats)
++ \return 0
++ */
++int I_bands_to_id_scatt(const int b_1_id, const int b_2_id, const int n_bands, int * scatt_id)
++{
++ int n_b1 = n_bands - 1;
+
- return nlines
-
- def MoveSelectedVertex(self, point, move):
-@@ -571,12 +795,21 @@
-
- if len(self._display.selected['ids']) != 1:
- return -1
--
++ * scatt_id = (b_1_id * (2 * n_b1 + 1) - b_1_id * b_1_id) / 2 + b_2_id - b_1_id - 1;
+
-+ # move only first found vertex in bbox
-+ poList = self._display.GetSelectedIList()
++ return 0;
++}
+
-+ if self.emit_signals:
-+ cList = poList.contents
-+ old_bboxs = [self._getBbox(cList.value[0])]
-+ old_areas_cats = [self._getLineAreasCategories(cList.value[0])]
++/*!
++ \brief Initialize structure for storing scatter plots data.
+
-+ Vect_set_updated(self.poMapInfo, 1)
-+ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
++ \param cats pointer to scCats struct
++ \param n_bands number of bands
++ \param type SC_SCATT_DATA - stores scatter plots
++ \param type SC_SCATT_CONDITIONS - stores selected areas in scatter plots
++ */
++void I_sc_init_cats(struct scCats * cats, int n_bands, int type)
++{
++ int i_cat;
+
- Vect_reset_line(self.poPoints)
- Vect_append_point(self.poPoints, point[0], point[1], 0.0)
--
-- # move only first found vertex in bbox
-- poList = self._display.GetSelectedIList()
++ cats->type = type;
+
- moved = Vedit_move_vertex(self.poMapInfo, self.popoBgMapInfo, int(self.poBgMapInfo is not None),
- poList, self.poPoints,
- self._display.GetThreshold(type = 'selectThresh'),
-@@ -584,7 +817,17 @@
- move[0], move[1], 0.0,
- 1, self._getSnapMode())
- Vect_destroy_list(poList)
--
++ cats->n_cats = 100;
++ cats->n_a_cats = 0;
+
-+ if moved > 0 and self.emit_signals:
-+ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
++ cats->n_bands = n_bands;
++ cats->n_scatts = (n_bands - 1) * n_bands / 2;
+
-+ new_bboxs = []
-+ new_areas_cats = []
-+ for i in range(n_up_lines_old, n_up_lines):
-+ new_id = Vect_get_updated_line(self.poMapInfo, i)
-+ new_bboxs.append(self._getBbox(new_id))
-+ new_areas_cats.append(self._getLineAreasCategories(new_id))
++ cats->cats_arr = (struct scScatts **) G_malloc(cats->n_cats * sizeof(struct scScatts *));
++ memset(cats->cats_arr, 0, cats-> n_cats * sizeof(struct scScatts *));
+
- if moved > 0 and self._settings['breakLines']:
- self._breakLineAtIntersection(Vect_get_num_lines(self.poMapInfo),
- None)
-@@ -592,7 +835,13 @@
- if moved > 0:
- self._addChangeset()
- self.toolbar.EnableUndo()
--
++ cats->cats_ids = (int *) G_malloc(cats->n_cats * sizeof(int));
++ cats->cats_idxs =(int *) G_malloc(cats->n_cats * sizeof(int));
+
-+ if self.emit_signals:
-+ self.vertexMoved.emit(new_bboxs = new_bboxs,
-+ new_areas_cats = new_areas_cats,
-+ old_areas_cats = old_areas_cats,
-+ old_bboxs = old_bboxs)
++ for(i_cat = 0; i_cat < cats->n_cats; i_cat++)
++ cats->cats_idxs[i_cat] = -1;
+
- return moved
-
- def AddVertex(self, coords):
-@@ -681,6 +930,10 @@
- self._error.ReadLine(line)
- return -1
-
-+ if self.emit_signals:
-+ old_bboxs = [self._getBbox(line)]
-+ old_areas_cats = [self._getLineAreasCategories(line)]
++ return;
++}
+
- # build feature geometry
- Vect_reset_line(self.poPoints)
- for p in coords:
-@@ -696,6 +949,9 @@
-
- newline = Vect_rewrite_line(self.poMapInfo, line, ltype,
- self.poPoints, self.poCats)
-+ if newline > 0 and self.emit_signals:
-+ new_geom = [self._getBbox(newline)]
-+ new_areas_cats = [self._getLineAreasCategories(newline)]
-
- if newline > 0 and self._settings['breakLines']:
- self._breakLineAtIntersection(newline, None)
-@@ -703,7 +959,13 @@
- if newline > 0:
- self._addChangeset()
- self.toolbar.EnableUndo()
--
++/*!
++ \brief Free data of struct scCats, the structure itself remains alocated.
++
++ \param cats pointer to existing scCats struct
++ */
++void I_sc_free_cats(struct scCats * cats)
++{
++ int i_cat;
++
++ for(i_cat = 0; i_cat < cats->n_a_cats; i_cat++)
++ {
++ if(cats->cats_arr[i_cat])
++ {
++ G_free(cats->cats_arr[i_cat]->scatt_idxs);
++ G_free(cats->cats_arr[i_cat]->scatts_bands);
++ G_free(cats->cats_arr[i_cat]->scatts_arr);
++ G_free(cats->cats_arr[i_cat]);
++ }
++ }
++
++ G_free(cats->cats_ids);
++ G_free(cats->cats_idxs);
++ G_free(cats->cats_arr);
++
++ cats->n_cats = 0;
++ cats->n_a_cats = 0;
++ cats->n_bands = 0;
++ cats->n_scatts = 0;
++ cats->type = -1;
++
++ return;
++}
++
++#if 0
++void I_sc_get_active_categories(int * a_cats_ids, int * n_a_cats, struct scCats * cats)
++{
++ a_cats_ids = cats->cats_ids;
++ * n_a_cats = cats->n_a_cats;
++}
++#endif
++
++/*!
++ \brief Add category.
+
-+ if self.emit_signals:
-+ self.lineEdited.emit(old_bboxs = old_bboxs,
-+ old_areas_cats = old_areas_cats,
-+ new_bboxs = new_bboxs,
-+ new_areas_cats = new_areas_cats)
++ Category represents group of scatter plots.
+
- return newline
++ \param cats pointer to scCats struct
++
++ \return assigned category id (starts with 0)
++ \return -1 if maximum nuber of categories was reached
++ */
++int I_sc_add_cat(struct scCats * cats)
++{
++ int i_scatt, i_cat_id, cat_id;
++ int n_a_cats = cats->n_a_cats;
++
++ if(cats->n_a_cats >= cats->n_cats)
++ return -1;
++
++ for(i_cat_id = 0; i_cat_id < cats->n_cats; i_cat_id++)
++ if(cats->cats_idxs[i_cat_id] < 0) {
++ cat_id = i_cat_id;
++ break;
++ }
++
++ cats->cats_ids[n_a_cats] = cat_id;
++ cats->cats_idxs[cat_id] = n_a_cats;
++
++ cats->cats_arr[n_a_cats] = (struct scScatts *) G_malloc(sizeof(struct scScatts));
++
++ cats->cats_arr[n_a_cats]->scatts_arr = (struct scdScattData **) G_malloc(cats->n_scatts * sizeof(struct scdScattData *));
++ memset((cats->cats_arr[n_a_cats]->scatts_arr), 0, cats->n_scatts * sizeof(struct scdScattData *));
++
++ cats->cats_arr[n_a_cats]->n_a_scatts = 0;
++
++ cats->cats_arr[n_a_cats]->scatts_bands = (int *) G_malloc(cats->n_scatts * 2 * sizeof(int));
++
++ cats->cats_arr[n_a_cats]->scatt_idxs = (int *) G_malloc(cats->n_scatts * sizeof(int));
++ for(i_scatt = 0; i_scatt < cats->n_scatts; i_scatt++)
++ cats->cats_arr[n_a_cats]->scatt_idxs[i_scatt] = -1;
++
++ ++cats->n_a_cats;
++
++ return cat_id;
++}
++
++#if 0
++int I_sc_delete_cat(struct scCats * cats, int cat_id)
++{
++ int cat_idx, i_cat;
++
++ if(cat_id < 0 || cat_id >= cats->n_cats)
++ return -1;
++
++ cat_idx = cats->cats_idxs[cat_id];
++ if(cat_idx < 0)
++ return -1;
++
++ G_free(cats->cats_arr[cat_idx]->scatt_idxs);
++ G_free(cats->cats_arr[cat_idx]->scatts_bands);
++ G_free(cats->cats_arr[cat_idx]->scatts_arr);
++ G_free(cats->cats_arr[cat_idx]);
++
++ for(i_cat = cat_idx; i_cat < cats->n_a_cats - 1; i_cat++)
++ {
++ cats->cats_arr[i_cat] = cats->cats_arr[i_cat + 1];
++ cats->cats_ids[i_cat] = cats->cats_ids[i_cat + 1];
++ }
++ cats->cats_idxs[cat_id] = -1;
++
++ --cats->n_a_cats;
++
++ return 0;
++}
++#endif
++
++/*!
++ \brief Insert scatter plot data .
++ Inserted scatt_data struct must have same type as cats struct (SC_SCATT_DATA or SC_SCATT_CONDITIONS).
++
++ \param cats pointer to scCats struct
++ \param cat_id id number of category.
++ \param scatt_id id number of scatter plot.
++
++ \return 0 on success
++ \return -1 on failure
++ */
++int I_sc_insert_scatt_data(struct scCats * cats, struct scdScattData * scatt_data, int cat_id, int scatt_id)
++{
++ int band_1, band_2, cat_idx, n_a_scatts;
++ struct scScatts * scatts;
++
++ if(cat_id < 0 || cat_id >= cats->n_cats)
++ return -1;
++
++ cat_idx = cats->cats_idxs[cat_id];
++ if(cat_idx < 0)
++ return -1;
++
++ if(scatt_id < 0 && scatt_id >= cats->n_scatts)
++ return -1;
++
++ scatts = cats->cats_arr[cat_idx];
++ if(scatts->scatt_idxs[scatt_id] >= 0)
++ return -1;
++
++ if(!scatt_data->b_conds_arr && cats->type == SC_SCATT_CONDITIONS)
++ return -1;
++
++ if(!scatt_data->scatt_vals_arr && cats->type == SC_SCATT_DATA)
++ return -1;
++
++ n_a_scatts = scatts->n_a_scatts;
++
++ scatts->scatt_idxs[scatt_id] = n_a_scatts;
++
++ I_id_scatt_to_bands(scatt_id, cats->n_bands, &band_1, &band_2);
++
++ scatts->scatts_bands[n_a_scatts * 2] = band_1;
++ scatts->scatts_bands[n_a_scatts * 2 + 1] = band_2;
++
++ scatts->scatts_arr[n_a_scatts] = scatt_data;
++ ++scatts->n_a_scatts;
++
++ return 0;
++}
++
++#if 0
++int I_sc_remove_scatt_data(struct scCats * cats, struct scdScattData * scatt_data, int cat_id, int scatt_id)
++{
++ int cat_idx, scatt_idx, n_init_scatts, i_scatt;
++ struct scScatts * scatts;
++
++ if(cat_id < 0 && cat_id >= cats->n_cats)
++ return -1;
++
++ cat_idx = cats->cats_idxs[cat_id];
++ if(cat_idx < 0)
++ return -1;
++
++ if(scatt_id < 0 || scatt_id >= cats->n_scatts)
++ return -1;
++
++ scatts = cats->cats_arr[cat_idx];
++ if(scatts->scatt_idxs[scatt_id] < 0)
++ return -1;
++
++ scatt_data = scatts->scatts_arr[scatt_idx];
++
++ for(i_scatt = scatt_idx; i_scatt < scatts->n_a_scatts - 1; i_scatt++)
++ {
++ scatts->scatts_arr[i_scatt] = scatts->scatts_arr[i_scatt + 1];
++ scatts->scatts_bands[i_scatt * 2] = scatts->scatts_bands[(i_scatt + 1)* 2];
++ scatts->scatts_bands[i_scatt * 2 + 1] = scatts->scatts_bands[(i_scatt + 1) * 2 + 1];
++ }
++ scatts->scatts_arr[scatts->n_a_scatts] = NULL;
++
++ scatts->scatt_idxs[scatt_id] = -1;
++
++ scatt_data = scatts->scatts_arr[scatt_id];
++ scatts->n_a_scatts--;
++
++ return 0;
++}
++
++int I_sc_set_value(struct scCats * cats, int cat_id, int scatt_id, int value_idx, int value)
++{
++ int n_a_scatts = cats->cats_arr[cat_id]->n_a_scatts;
++ int cat_idx, scatt_idx, ret;
++
++ cat_idx = cats->cats_idxs[cat_id];
++ if(cat_idx < 0)
++ return -1;
++
++ if(cats->cats_arr[cat_idx]->scatt_idxs[scatt_id] < 0)
++ return -1;
++
++ cat_idx = cats->cats_idxs[cat_id];
++ scatt_idx = cats->cats_arr[cat_idx]->scatt_idxs[scatt_id];
++
++ I_scd_set_value(cats->cats_arr[cat_idx]->scatts_arr[scatt_idx], value_idx, value);
++
++ return 0;
++}
++#endif
++
++/*!
++ \brief Insert scatter plot data.
++
++ \param scatt_data pointer to existing struct scdScattData
++ \param type SC_SCATT_DATA for scatter plots or SC_SCATT_CONDITIONS for selected areas in scatter plot
++ \param n_vals number of data values
++ \param data array of values (unsigned char for SC_SCATT_CONDITIONS, unsigned int for SC_SCATT_DATA)
++ */
++void I_scd_init_scatt_data(struct scdScattData * scatt_data, int type, int n_vals, void * data)
++{
++ scatt_data->n_vals = n_vals;
++
++ if(type == SC_SCATT_DATA)
++ {
++ if(data)
++ scatt_data->scatt_vals_arr = (unsigned int *) data;
++ else {
++ scatt_data->scatt_vals_arr = (unsigned int *) G_malloc(n_vals * sizeof(unsigned int));
++ memset(scatt_data->scatt_vals_arr, 0, n_vals * sizeof(unsigned int));
++ }
++ scatt_data->b_conds_arr = NULL;
++ }
++ else if(type == SC_SCATT_CONDITIONS)
++ {
++ if(data)
++ scatt_data->b_conds_arr = (unsigned char *) data;
++ else {
++ scatt_data->b_conds_arr = (unsigned char *) G_malloc(n_vals * sizeof(unsigned char));
++ memset(scatt_data->b_conds_arr, 0, n_vals * sizeof(unsigned char));
++ }
++ scatt_data->scatt_vals_arr = NULL;
++ }
++
++ return;
++}
++
++
++#if 0
++void I_scd_get_range_min_max(struct scdScattData * scatt_data, CELL * band_1_min, CELL * band_1_max, CELL * band_2_min, CELL * band_2_max)
++{
++
++ Rast_get_range_min_max(&(scatt_data->band_1_range), band_1_min, band_2_min);
++ Rast_get_range_min_max(&(scatt_data->band_2_range), band_2_min, band_2_max);
++
++ return;
++}
++s
++void * I_scd_get_data_ptr(struct scdScattData * scatt_data)
++{
++ if(!scatt_data->b_conds_arr)
++ return scatt_data->b_conds_arr;
++ else if(!scatt_data->scatt_vals_arr)
++ return scatt_data->scatt_vals_arr;
++
++ return NULL;
++}
++
++int I_scd_set_value(struct scdScattData * scatt_data, unsigned int val_idx, unsigned int val)
++{
++ if(val_idx < 0 && val_idx > scatt_data->n_vals)
++ return -1;
++
++ if(scatt_data->b_conds_arr)
++ scatt_data->b_conds_arr[val_idx] = val;
++ else if(scatt_data->scatt_vals_arr)
++ scatt_data->scatt_vals_arr[val_idx] = val;
++ else
++ return -1;
++
++ return 0;
++}
++#endif
+Index: include/imagery.h
+===================================================================
+--- include/imagery.h (revision 57599)
++++ include/imagery.h (working copy)
+@@ -135,6 +135,56 @@
+
+ } IClass_statistics;
- def FlipLine(self):
-@@ -1514,6 +1776,16 @@
- return 0
-
- poList = self._display.GetSelectedIList()
++/* Scatter Plot backend */
+
-+ if self.emit_signals:
-+ cList = poList.contents
-+
-+ old_bboxs = [self._getBbox(cList.value[0])]
-+ old_areas_cats = [self._getLineAreasCategories(cList.value[0])]
++#define SC_SCATT_DATA 0
++#define SC_SCATT_CONDITIONS 1
+
-+ Vect_set_updated(self.poMapInfo, 1)
-+ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
+
- Vect_reset_line(self.poPoints)
- Vect_append_point(self.poPoints, coords[0], coords[1], 0.0)
-
-@@ -1525,15 +1797,35 @@
- else:
- ret = Vedit_remove_vertex(self.poMapInfo, poList,
- self.poPoints, thresh)
++/*! Holds list of all catagories.
++ It can contain selected areas for scatter plots (SC_SCATT_CONDITIONS type) or computed scatter plots (SC_SCATT_DATA type).
++*/
++struct scCats
++{
++ int type; /*!< SC_SCATT_DATA -> computed scatter plots, SC_SCATT_CONDITIONS -> set conditions for scatter plots to be computed*/
+
- Vect_destroy_list(poList)
++ int n_cats; /*!< number of alocated categories */
++
++ int n_bands; /*!< number of analyzed bands */
++ int n_scatts; /*!< number of possible scattter plot which can be created from bands */
+
-+ if ret > 0 and self.emit_signals:
-+ new_bboxs = []
-+ new_areas_cats = []
++ int n_a_cats; /*!< number of used/active categories */
++ int * cats_ids; /*!< (cat_idx->cat_id) array index is internal idx (position in cats_arr) and id is saved in it's position*/
++ int * cats_idxs; /*!< (cat_id->cat_idx) array index is id and internal idx is saved in it's position*/
+
-+ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
-+ for i in range(n_up_lines_old, n_up_lines):
-+ new_id = Vect_get_updated_line(self.poMapInfo, i)
-+ new_areas_cats.append(self._getLineAreasCategories(new_id))
-+ new_bboxs.append(self._getBbox(new_id))
-
- if not add and ret > 0 and self._settings['breakLines']:
- self._breakLineAtIntersection(Vect_get_num_lines(self.poMapInfo),
- None)
--
++ struct scScatts ** cats_arr; /*!< array of pointers to struct scScatts */
++};
+
- if ret > 0:
- self._addChangeset()
--
+
-+ if ret > 0 and self.emit_signals:
-+ if add:
-+ self.vertexAdded.emit(old_bboxs = old_bboxs, new_bboxs = new_bboxs)
-+ else:
-+ self.vertexRemoved.emit(old_bboxs = old_bboxs,
-+ new_bboxs = new_bboxs,
-+ old_areas_cats = old_areas_cats,
-+ new_areas_cats = new_areas_cats)
++/*! Holds list of all scatter plots, which belongs to category.
++*/
++struct scScatts
++{
++ int n_a_scatts; /*!< number of used/active scatter plots*/
++
++ int * scatts_bands; /*!< array of bands for which represents the scatter plots,
++ every scatter plot has assigned two bads (size of array is n_a_scatts * 2)*/
++ int * scatt_idxs; /*!< (scatt_id->scatt_idx) internal idx of the scatter plot (position in scatts_arr)*/
+
- return 1
-
- def GetLineCats(self, line):
-Index: gui/wxpython/vdigit/toolbars.py
++ struct scdScattData ** scatts_arr; /*!< array of pointers to scdScattData */
++};
++
++/*! Holds scatter plot data.
++*/
++struct scdScattData
++{
++ int n_vals; /*!< Number of values in scatter plot. */
++
++ unsigned char * b_conds_arr; /*!< array of selected areas (used for SC_SCATT_CONDITIONS type) otherwise NULL */
++ unsigned int * scatt_vals_arr; /*!< array of computed areas (used for SC_SCATT_DATA type) otherwise NULL */
++};
++
++
+ #define SIGNATURE_TYPE_MIXED 1
+
+ #define GROUPFILE "CURGROUP"
+Index: include/defs/vedit.h
===================================================================
---- gui/wxpython/vdigit/toolbars.py (revision 57581)
-+++ gui/wxpython/vdigit/toolbars.py (working copy)
-@@ -17,6 +17,7 @@
- import wx
+--- include/defs/vedit.h (revision 57599)
++++ include/defs/vedit.h (working copy)
+@@ -33,6 +33,8 @@
+ int Vedit_merge_lines(struct Map_info *, struct ilist *);
- from grass.script import core as grass
-+from grass.pydispatch.signal import Signal
+ /* move.c */
++int Vedit_move_areas(struct Map_info *, struct Map_info **, int,
++ struct ilist *, double, double, double, int, double);
+ int Vedit_move_lines(struct Map_info *, struct Map_info **, int,
+ struct ilist *, double, double, double, int, double);
- from gui_core.toolbars import BaseToolbar, BaseIcons
- from gui_core.dialogs import CreateNewVector
-@@ -42,6 +43,8 @@
- self.digit = None
- self._giface = giface
-
-+ self.editingStarted = Signal("VDigitToolbar.editingStarted")
-+
- # currently selected map layer for editing (reference to MapLayer instance)
- self.mapLayer = None
- # list of vector layers from Layer Manager (only in the current mapset)
-@@ -860,6 +863,7 @@
- alpha = int(opacity * 255)
- self.digit.GetDisplay().UpdateSettings(alpha = alpha)
-
-+ self.editingStarted.emit(vectMap = mapLayer.GetName(), digit = self.digit)
- return True
+Index: include/defs/imagery.h
+===================================================================
+--- include/defs/imagery.h (revision 57599)
++++ include/defs/imagery.h (working copy)
+@@ -110,6 +110,23 @@
+ FILE *I_fopen_subgroup_ref_new(const char *, const char *);
+ FILE *I_fopen_subgroup_ref_old(const char *, const char *);
- def StopEditing(self):
-Index: gui/wxpython/iclass/dialogs.py
-===================================================================
---- gui/wxpython/iclass/dialogs.py (revision 57581)
-+++ gui/wxpython/iclass/dialogs.py (working copy)
-@@ -333,13 +333,19 @@
- toolbar.SetCategories(catNames = catNames, catIdx = cats)
- if name in catNames:
- toolbar.choice.SetStringSelection(name)
-+ cat = toolbar.GetSelectedCategoryIdx()
- elif catNames:
- toolbar.choice.SetSelection(0)
--
-+ cat = toolbar.GetSelectedCategoryIdx()
-+ else:
-+ cat = None
++/* scatt_plt.c */
++void I_sc_init_cats(struct scCats *, int, int);
++void I_sc_free_cats(struct scCats *);
++int I_sc_add_cat(struct scCats *);
++int I_sc_insert_scatt_data(struct scCats *, struct scdScattData *, int, int);
+
- if toolbar.choice.IsEmpty():
- toolbar.EnableControls(False)
- else:
- toolbar.EnableControls(True)
++void I_scd_init_scatt_data(struct scdScattData *, int, int, void *);
+
-+ self.mapWindow.CategoryChanged(cat)
- # don't forget to update maps, histo, ...
-
- def GetSelectedIndices(self, state = wx.LIST_STATE_SELECTED):
++int I_compute_scatts(struct Cell_head *, struct scCats *, const char **,
++ const char **, int, struct scCats *, const char **);
++
++int I_create_cat_rast(struct Cell_head *, const char *);
++int I_insert_patch_to_cat_rast(const char *, struct Cell_head *, const char *);
++
++int I_id_scatt_to_bands(const int, const int, int *, int *);
++int I_bands_to_id_scatt(const int, const int, const int, int *);
++
+ /* sig.c */
+ int I_init_signatures(struct Signature *, int);
+ int I_new_signature(struct Signature *);
Index: gui/wxpython/iclass/toolbars.py
===================================================================
---- gui/wxpython/iclass/toolbars.py (revision 57581)
+--- gui/wxpython/iclass/toolbars.py (revision 57599)
+++ gui/wxpython/iclass/toolbars.py (working copy)
@@ -46,9 +46,7 @@
'importAreas' : MetaIcon(img = 'layer-import',
@@ -712,7 +1315,7 @@
self.parent.OnCategoryManager),
Index: gui/wxpython/iclass/frame.py
===================================================================
---- gui/wxpython/iclass/frame.py (revision 57581)
+--- gui/wxpython/iclass/frame.py (revision 57599)
+++ gui/wxpython/iclass/frame.py (working copy)
@@ -64,6 +64,8 @@
IClassExportAreasDialog, IClassMapDialog
@@ -752,7 +1355,7 @@
def RemoveTempVector(self):
"""!Removes temporary vector map with training areas"""
ret = RunCommand(prog = 'g.remove',
-@@ -477,7 +483,7 @@
+@@ -477,21 +483,47 @@
self.Render(self.GetFirstWindow())
@@ -760,18 +1363,51 @@
+ def AddBands(self):
"""!Add imagery group"""
dlg = IClassGroupDialog(self, group = self.group)
- if dlg.ShowModal() == wx.ID_OK:
-@@ -488,7 +494,8 @@
+- if dlg.ShowModal() == wx.ID_OK:
+- self.SetGroup(dlg.GetGroup())
++
++ while True:
++ if dlg.ShowModal() == wx.ID_OK:
++ if self.SetGroup(dlg.GetGroup()):
++ break
++ else:
++ break
++
+ dlg.Destroy()
+
+ def SetGroup(self, name):
"""!Set imagery group"""
group = grass.find_file(name = name, element = 'group')
if group['name']:
-- self.group = group['name']
-+ self.group = group['name']
-+ self.groupSet.emit(group = group['name'])
++ if not self.GroupData(group['name']):
++ GError(_("No data found in group <%s>.\n" \
++ "Please create subgroup with same name as group and add the data there.") \
++ % name, parent = self)
++ return False
+ self.group = group['name']
++ self.groupSet.emit(group = group['name'])
else:
GError(_("Group <%s> not found") % name, parent = self)
-
-@@ -768,17 +775,20 @@
+-
++ return False
++
++ return True
++
++ def GroupData(self, group):
++ res = RunCommand('i.group',
++ flags = 'g',
++ group = group, subgroup = group,
++ read = True).strip()
++ bands = None
++ if res.split('\n')[0]:
++ bands = res.split('\n')
++
++ return bands
++
+ def OnImportAreas(self, event):
+ """!Import training areas"""
+ # check if we have any changes
+@@ -768,17 +800,20 @@
Updates number of stddev, histograms, layer in preview display.
"""
@@ -802,7 +1438,7 @@
def DeleteAreas(self, cats):
"""!Removes all training areas of given categories
-@@ -1105,27 +1115,6 @@
+@@ -1105,27 +1140,6 @@
self.GetFirstWindow().SetModePointer()
self.GetSecondWindow().SetModePointer()
@@ -832,7 +1468,7 @@
Index: gui/wxpython/iclass/plots.py
===================================================================
---- gui/wxpython/iclass/plots.py (revision 57581)
+--- gui/wxpython/iclass/plots.py (revision 57599)
+++ gui/wxpython/iclass/plots.py (working copy)
@@ -19,6 +19,7 @@
import wx.lib.plot as plot
@@ -979,1430 +1615,36 @@
def DrawCoincidencePlots(self):
"""!Draw coincidence plots"""
for bandIdx in range(len(self.bandList)):
-Index: gui/wxpython/mapdisp/frame.py
+Index: gui/wxpython/iclass/dialogs.py
===================================================================
---- gui/wxpython/mapdisp/frame.py (revision 57581)
-+++ gui/wxpython/mapdisp/frame.py (working copy)
-@@ -225,6 +225,7 @@
- #
- self.dialogs = {}
- self.dialogs['attributes'] = None
-+ self.dialogs['scatt_plot'] = None
- self.dialogs['category'] = None
- self.dialogs['barscale'] = None
- self.dialogs['legend'] = None
-@@ -1168,6 +1169,19 @@
- """!Returns toolbar with zooming tools"""
- return self.toolbars['map']
-
-+ def OnScatterplot2(self, event):
-+ """!Init interactive scatterplot tools
-+ """
-+ if self.dialogs['scatt_plot']:
-+ self.dialogs['scatt_plot'].Raise()
-+ return
+--- gui/wxpython/iclass/dialogs.py (revision 57599)
++++ gui/wxpython/iclass/dialogs.py (working copy)
+@@ -333,13 +333,19 @@
+ toolbar.SetCategories(catNames = catNames, catIdx = cats)
+ if name in catNames:
+ toolbar.choice.SetStringSelection(name)
++ cat = toolbar.GetSelectedCategoryIdx()
+ elif catNames:
+ toolbar.choice.SetSelection(0)
+-
++ cat = toolbar.GetSelectedCategoryIdx()
++ else:
++ cat = None
+
-+ from scatt_plot.dialogs import ScattPlotMainDialog
-+ self.dialogs['scatt_plot'] = ScattPlotMainDialog(parent=self, giface=self._giface)
-+
-+ self.dialogs['scatt_plot'].CenterOnScreen()
-+ self.dialogs['scatt_plot'].Show()
+ if toolbar.choice.IsEmpty():
+ toolbar.EnableControls(False)
+ else:
+ toolbar.EnableControls(True)
+
- def OnVNet(self, event):
- """!Dialog for v.net* modules
- """
-Index: gui/wxpython/mapdisp/toolbars.py
-===================================================================
---- gui/wxpython/mapdisp/toolbars.py (revision 57581)
-+++ gui/wxpython/mapdisp/toolbars.py (working copy)
-@@ -239,7 +239,8 @@
- (MapIcons["scatter"], self.parent.OnScatterplot),
- (MapIcons["histogram"], self.parent.OnHistogramPyPlot),
- (BaseIcons["histogramD"], self.parent.OnHistogram),
-- (MapIcons["vnet"], self.parent.OnVNet)))
-+ (MapIcons["vnet"], self.parent.OnVNet),
-+ (MapIcons["scatter"], self.parent.OnScatterplot2)))
++ self.mapWindow.CategoryChanged(cat)
+ # don't forget to update maps, histo, ...
- def OnDecoration(self, event):
- """!Decorations overlay menu
-Index: gui/wxpython/scatt_plot/plots.py
-===================================================================
---- gui/wxpython/scatt_plot/plots.py (revision 0)
-+++ gui/wxpython/scatt_plot/plots.py (working copy)
-@@ -0,0 +1,695 @@
-+"""!
-+ at package scatt_plot.dialogs
-+
-+ at brief Ploting widgets.
-+
-+Classes:
-+
-+(C) 2013 by the GRASS Development Team
-+
-+This program is free software under the GNU General Public License
-+(>=v2). Read the file COPYING that comes with GRASS for details.
-+
-+ at author Stepan Turek <stepan.turek seznam.cz> (mentor: Martin Landa)
-+"""
-+import wx
-+import numpy as np
-+
-+#TODO testing
-+import time
-+from multiprocessing import Process, Queue
-+
-+#TODO
-+#from numpy.lib.stride_tricks import as_strided
-+from copy import deepcopy
-+
-+try:
-+ import matplotlib
-+ matplotlib.use('WXAgg')
-+ from matplotlib.figure import Figure
-+ from matplotlib.backends.backend_wxagg import \
-+ FigureCanvasWxAgg as FigCanvas, \
-+ NavigationToolbar2WxAgg as NavigationToolbar
-+ from matplotlib.lines import Line2D
-+ from matplotlib.artist import Artist
-+ from matplotlib.mlab import dist_point_to_segment
-+ from matplotlib.patches import Polygon, Ellipse
-+ import matplotlib.image as mi
-+ import matplotlib.colors as mcolors
-+ import matplotlib.cbook as cbook
-+except ImportError as e:
-+ raise ImportError(_("Unable to import matplotlib (try to install it).\n%s") % e)
-+
-+
-+import grass.script as grass
-+from grass.pydispatch.signal import Signal
-+
-+#class PlotImages()?
-+
-+class ScatterPlotWidget(wx.Panel):
-+ def __init__(self, parent, scatt_id, scatt_mgr,
-+ id = wx.ID_ANY):
-+
-+ wx.Panel.__init__(self, parent, id)
-+
-+ self.parent = parent
-+ self.full_extend = None
-+
-+ self._createWidgets()
-+ self._doLayout()
-+ self.scatt_id = scatt_id
-+ self.scatt_mgr = scatt_mgr
-+ self.press_coords = None
-+
-+ self.cidpress = None
-+ self.cidrelease = None
-+
-+ self.SetSize((200, 100))
-+ self.Layout()
-+
-+ self.base_scale = 2.0
-+ self.Bind(wx.EVT_CLOSE,lambda event : self.CleanUp())
-+
-+ self.plotClosed = Signal("ScatterPlotWidget.plotClosed")
-+
-+ self.contex_menu = ScatterPlotContextMenu(plot = self)
-+
-+ self.ciddscroll = None
-+
-+ def _createWidgets(self):
-+
-+ # Create the mpl Figure and FigCanvas objects.
-+ # 5x4 inches, 100 dots-per-inch
-+ #
-+ self.dpi = 100
-+ self.fig = Figure((1.0, 1.0), dpi=self.dpi)
-+ self.canvas = FigCanvas(self, -1, self.fig)
-+
-+ self.axes = self.fig.add_axes([0,0,1,1])
-+
-+ pol = Polygon(list(zip([0], [0])), animated=True)
-+ self.axes.add_patch(pol)
-+ self.polygon_drawer = PolygonDrawer(self.axes, pol = pol, empty_pol = True)
-+
-+ self.toolbar = NavigationToolbar(self.canvas)
-+
-+ def ZoomToExtend(self):
-+ if self.full_extend:
-+ self.axes.axis(self.full_extend)
-+ self.canvas.draw()
-+
-+ def SetMode(self, mode):
-+ self._deactivateMode()
-+ if mode == 'zoom':
-+ self.ciddscroll = self.canvas.mpl_connect('scroll_event', self.zoom)
-+ elif mode == 'pan':
-+ self.toolbar.pan()
-+
-+ def GetCoords(self):
-+ return self.polygon_drawer.GetCoords()
-+
-+ def SetEmpty(self):
-+ return self.polygon_drawer.SetEmpty()
-+
-+ def SetEditingMode(self, mode):
-+ self.polygon_drawer.SetMode(mode)
-+
-+ def _deactivateMode(self):
-+
-+ #TODO do own pan
-+ if self.toolbar._active == "PAN":
-+ self.toolbar.pan()
-+
-+ if self.ciddscroll:
-+ self.canvas.mpl_disconnect(self.ciddscroll)
-+
-+ self._stopCategoryEdit()
-+
-+
-+ def OnPress(self, event):
-+ 'on button press we will see if the mouse is over us and store some data'
-+ if not event.inaxes:
-+ return
-+
-+ if event.xdata and event.ydata:
-+ self.press_coords = { 'x' : event.xdata, 'y' : event.ydata}
-+ else:
-+ self.press_coords = None
-+
-+ def OnRelease(self, event):
-+ 'on release we reset the press data'
-+
-+ if event.xdata and event.ydata and self.press_coords:
-+
-+ bbox = {}
-+ if event.ydata > self.press_coords['y']:
-+ bbox['up_y'] = event.ydata
-+ bbox['btm_y'] = self.press_coords['y']
-+ else:
-+ bbox['up_y'] = self.press_coords['y']
-+ bbox['btm_y'] = event.ydata
-+
-+ if event.xdata > self.press_coords['x']:
-+ bbox['up_x'] = event.xdata
-+ bbox['btm_x'] = self.press_coords['x']
-+ else:
-+ bbox['up_x'] = self.press_coords['x']
-+ bbox['btm_x'] = event.xdata
-+
-+ self.scatt_mgr.SetEditCatData(self.scatt_id, bbox)
-+
-+ def _stopCategoryEdit(self):
-+ 'disconnect all the stored connection ids'
-+
-+ if self.cidpress:
-+ self.canvas.mpl_disconnect(self.cidpress)
-+ if self.cidrelease:
-+ self.canvas.mpl_disconnect(self.cidrelease)
-+ #self.canvas.mpl_disconnect(self.cidmotion)
-+
-+ def _doLayout(self):
-+
-+ self.main_sizer = wx.BoxSizer(wx.VERTICAL)
-+ self.main_sizer.Add(self.canvas, 1, wx.LEFT | wx.TOP | wx.GROW)
-+ #self.main_sizer.Add(self.toolbar, 0, wx.EXPAND)
-+ self.SetSizer(self.main_sizer)
-+ self.main_sizer.Fit(self)
-+
-+ def Plot(self, scatts, ellipses, styles):
-+ """ Redraws the figure
-+ """
-+ self.axes.clear()
-+
-+ callafter_list = []
-+
-+ q = Queue()
-+ p = Process(target=MergeImg, args=(scatts, styles, q))
-+ p.start()
-+ merged_img, self.full_extend = q.get()
-+ p.join()
-+
-+ img = imshow(self.axes, merged_img,
-+ origin = 'lower',
-+ extent = self.full_extend,
-+ interpolation='nearest',
-+ aspect = "auto")
-+
-+ callafter_list.append([self.axes.draw_artist, [img]])
-+ callafter_list.append([grass.try_remove, [merged_img.filename]])
-+
-+ for cat_id, e in ellipses.iteritems():
-+ if cat_id == 0 or not e:
-+ continue
-+
-+ colors = styles[cat_id]['color'].split(":")
-+ ellip = Ellipse(xy=e['pos'], width=e['width'], height=e['height'], angle=e['theta'], edgecolor = "r", facecolor = 'None')
-+ self.axes.add_artist(ellip)
-+ callafter_list.append([self.axes.draw_artist, [ellip]])
-+
-+ callafter_list.append([self.fig.canvas.blit, []])
-+
-+ wx.CallAfter(lambda : self.CallAfter(callafter_list))
-+
-+ def CallAfter(self, funcs_list):
-+ while funcs_list:
-+ fcn, args = funcs_list.pop(0)
-+ fcn(*args)
-+
-+ self.canvas.draw()
-+
-+ def CleanUp(self):
-+ self.plotClosed.emit(scatt_id = self.scatt_id)
-+ self.Destroy()
-+
-+ def zoom(self, event):
-+ # get the current x and y limits
-+ if not event.inaxes:
-+ return
-+ # tcaswell
-+ # http://stackoverflow.com/questions/11551049/matplotlib-plot-zooming-with-scroll-wheel
-+ cur_xlim = self.axes.get_xlim()
-+ cur_ylim = self.axes.get_ylim()
-+ cur_xrange = (cur_xlim[1] - cur_xlim[0])*.5
-+ cur_yrange = (cur_ylim[1] - cur_ylim[0])*.5
-+ xdata = event.xdata
-+ ydata = event.ydata
-+ if event.button == 'up':
-+ scale_factor = 1/self.base_scale
-+ elif event.button == 'down':
-+ scale_factor = self.base_scale
-+ else:
-+ scale_factor = 1
-+
-+ self.axes.set_xlim([xdata - cur_xrange*scale_factor,
-+ xdata + cur_xrange*scale_factor])
-+ self.axes.set_ylim([ydata - cur_yrange*scale_factor,
-+ ydata + cur_yrange*scale_factor])
-+
-+ self.canvas.draw()
-+
-+class ScatterPlotContextMenu:
-+ def __init__(self, plot):
-+
-+ self.plot = plot
-+ self.canvas = plot.canvas
-+ self.cidpress = self.canvas.mpl_connect(
-+ 'button_press_event', self.ContexMenu)
-+
-+ def ContexMenu(self, event):
-+ if not event.inaxes:
-+ return
-+
-+ if event.button == 3:
-+ menu = wx.Menu()
-+ menu_items = [["zoom_to_extend", _("Zoom to scatter plot extend"), lambda event : self.plot.ZoomToExtend()]]
-+
-+ for item in menu_items:
-+ item_id = wx.ID_ANY
-+ menu.Append(item_id, text = item[1])
-+ menu.Bind(wx.EVT_MENU, item[2], id = item_id)
-+
-+ wx.CallAfter(self.ShowMenu, menu)
-+
-+ def ShowMenu(self, menu):
-+ self.plot.PopupMenu(menu)
-+ menu.Destroy()
-+ self.plot.ReleaseMouse()
-+
-+class PolygonDrawer:
-+ """
-+ An polygon editor.
-+ """
-+
-+ showverts = True
-+ epsilon = 5 #TODO settings max pixel distance to count as a vertex hit
-+ def __init__(self, ax, pol, empty_pol):
-+ if pol.figure is None:
-+ raise RuntimeError('You must first add the polygon to a figure or canvas before defining the interactor')
-+ self.ax = ax
-+ self.canvas = pol.figure.canvas
-+
-+ self.pol = pol
-+ self.empty_pol = empty_pol
-+
-+ x, y = zip(*self.pol.xy)
-+ self.line = Line2D(x, y, marker='o', markerfacecolor='r', animated=True)
-+ self.ax.add_line(self.line)
-+ #self._update_line(pol)
-+
-+ cid = self.pol.add_callback(self.poly_changed)
-+ self.moving_ver_idx = None # the active vert
-+
-+ self.mode = None
-+
-+ if self.empty_pol:
-+ self.Show(False)
-+
-+ self.canvas.mpl_connect('draw_event', self.draw_callback)
-+ self.canvas.mpl_connect('button_press_event', self.OnButtonPressed)
-+ self.canvas.mpl_connect('button_release_event', self.button_release_callback)
-+ self.canvas.mpl_connect('motion_notify_event', self.motion_notify_callback)
-+
-+ self.it = 0
-+ def SetMode(self, mode):
-+ self.mode = mode
-+
-+ def GetCoords(self):
-+ if self.empty_pol:
-+ return None
-+
-+ coords = deepcopy(self.pol.xy)
-+ return coords
-+
-+ def SetEmpty(self):
-+ self._setEmptyPol(True)
-+
-+ def _setEmptyPol(self, empty_pol):
-+ self.empty_pol = empty_pol
-+ self.Show(not empty_pol)
-+
-+ def Show(self, show):
-+
-+ self.show = show
-+
-+ self.line.set_visible(self.show)
-+ self.pol.set_visible(self.show)
-+
-+ self.Redraw()
-+
-+ def Redraw(self):
-+ if self.show:
-+ self.ax.draw_artist(self.pol)
-+ self.ax.draw_artist(self.line)
-+ self.canvas.blit(self.ax.bbox)
-+ self.canvas.draw()
-+
-+ def draw_callback(self, event):
-+ self.background = self.canvas.copy_from_bbox(self.ax.bbox)
-+ self.ax.draw_artist(self.pol)
-+ self.ax.draw_artist(self.line)
-+
-+ def poly_changed(self, pol):
-+ 'this method is called whenever the polygon object is called'
-+ # only copy the artist props to the line (except visibility)
-+ vis = self.line.get_visible()
-+ Artist.update_from(self.line, pol)
-+ self.line.set_visible(vis) # don't use the pol visibility state
-+
-+ def get_ind_under_point(self, event):
-+ 'get the index of the vertex under point if within epsilon tolerance'
-+
-+ # display coords
-+ xy = np.asarray(self.pol.xy)
-+ xyt = self.pol.get_transform().transform(xy)
-+ xt, yt = xyt[:, 0], xyt[:, 1]
-+ d = np.sqrt((xt-event.x)**2 + (yt-event.y)**2)
-+ indseq = np.nonzero(np.equal(d, np.amin(d)))[0]
-+ ind = indseq[0]
-+
-+ if d[ind]>=self.epsilon:
-+ ind = None
-+
-+ return ind
-+
-+ def OnButtonPressed(self, event):
-+ if not event.inaxes:
-+ return
-+
-+ if event.button in [2, 3]:
-+ return
-+
-+ if self.mode == "delete_vertex":
-+ self._deleteVertex(event)
-+ elif self.mode == "add_boundary_vertex":
-+ self._addVertexOnBoundary(event)
-+ elif self.mode == "add_vertex":
-+ self._addVertex(event)
-+ self.moving_ver_idx = self.get_ind_under_point(event)
-+
-+ def button_release_callback(self, event):
-+ 'whenever a mouse button is released'
-+ if not self.showverts: return
-+ if event.button != 1: return
-+ self.moving_ver_idx = None
-+
-+ def ShowVertices(self, show):
-+ self.showverts = show
-+ self.line.set_visible(self.showverts)
-+ if not self.showverts: self.moving_ver_idx = None
-+
-+ def _deleteVertex(self, event):
-+ ind = self.get_ind_under_point(event)
-+
-+ if ind is None or self.empty_pol:
-+ return
-+
-+ if len(self.pol.xy) <= 2:
-+ self.empty_pol = True
-+ self.Show(False)
-+ return
-+
-+ coords = []
-+ for i,tup in enumerate(self.pol.xy):
-+ if i == ind:
-+ continue
-+ elif i == 0 and ind == len(self.pol.xy) - 1:
-+ continue
-+ elif i == len(self.pol.xy) - 1 and ind == 0:
-+ continue
-+
-+ coords.append(tup)
-+
-+ self.pol.xy = coords
-+ self.line.set_data(zip(*self.pol.xy))
-+
-+ self.Redraw()
-+
-+ def _addVertexOnBoundary(self, event):
-+ if self.empty_pol:
-+ return
-+
-+ xys = self.pol.get_transform().transform(self.pol.xy)
-+ p = event.x, event.y # display coords
-+ for i in range(len(xys)-1):
-+ s0 = xys[i]
-+ s1 = xys[i+1]
-+ d = dist_point_to_segment(p, s0, s1)
-+
-+ if d<=self.epsilon:
-+ self.pol.xy = np.array(
-+ list(self.pol.xy[:i + 1]) +
-+ [(event.xdata, event.ydata)] +
-+ list(self.pol.xy[i + 1:]))
-+ self.line.set_data(zip(*self.pol.xy))
-+ break
-+
-+ self.Redraw()
-+
-+ def _addVertex(self, event):
-+
-+ if self.empty_pol:
-+ pt = (event.xdata, event.ydata)
-+ self.pol.xy = np.array([pt, pt])
-+ self.Show(True)
-+ self.empty_pol = False
-+ else:
-+ self.pol.xy = np.array(
-+ [(event.xdata, event.ydata)] +
-+ list(self.pol.xy[1:]) +
-+ [(event.xdata, event.ydata)])
-+
-+ self.line.set_data(zip(*self.pol.xy))
-+
-+ self.Redraw()
-+
-+ def motion_notify_callback(self, event):
-+ 'on mouse movement'
-+ if not self.mode == "move_vertex": return
-+ if not self.showverts: return
-+ if self.empty_pol: return
-+ if self.moving_ver_idx is None: return
-+ if event.inaxes is None: return
-+ if event.button != 1: return
-+
-+ self.it += 1
-+
-+ x,y = event.xdata, event.ydata
-+
-+ self.pol.xy[self.moving_ver_idx] = x,y
-+ if self.moving_ver_idx == 0:
-+ self.pol.xy[len(self.pol.xy) - 1] = x,y
-+ elif self.moving_ver_idx == len(self.pol.xy) - 1:
-+ self.pol.xy[0] = x,y
-+
-+ self.line.set_data(zip(*self.pol.xy))
-+
-+ self.canvas.restore_region(self.background)
-+
-+ self.Redraw()
-+
-+class ModestImage(mi.AxesImage):
-+ """
-+ Computationally modest image class.
-+
-+ ModestImage is an extension of the Matplotlib AxesImage class
-+ better suited for the interactive display of larger images. Before
-+ drawing, ModestImage resamples the data array based on the screen
-+ resolution and view window. This has very little affect on the
-+ appearance of the image, but can substantially cut down on
-+ computation since calculations of unresolved or clipped pixels
-+ are skipped.
-+
-+ The interface of ModestImage is the same as AxesImage. However, it
-+ does not currently support setting the 'extent' property. There
-+ may also be weird coordinate warping operations for images that
-+ I'm not aware of. Don't expect those to work either.
-+
-+ Author: Chris Beaumont <beaumont at hawaii.edu>
-+ """
-+ def __init__(self, *args, **kwargs):
-+ #why???
-+ #if 'extent' in kwargs and kwargs['extent'] is not None:
-+ # raise NotImplementedError("ModestImage does not support extents")
-+
-+ self._full_res = None
-+ self._sx, self._sy = None, None
-+ self._bounds = (None, None, None, None)
-+ super(ModestImage, self).__init__(*args, **kwargs)
-+
-+ def set_data(self, A):
-+ """
-+ Set the image array
-+
-+ ACCEPTS: numpy/PIL Image A
-+ """
-+ self._full_res = A
-+ self._A = A
-+
-+ if self._A.dtype != np.uint8 and not np.can_cast(self._A.dtype,
-+ np.float):
-+ raise TypeError("Image data can not convert to float")
-+
-+ if (self._A.ndim not in (2, 3) or
-+ (self._A.ndim == 3 and self._A.shape[-1] not in (3, 4))):
-+ raise TypeError("Invalid dimensions for image data")
-+
-+ self._imcache =None
-+ self._rgbacache = None
-+ self._oldxslice = None
-+ self._oldyslice = None
-+ self._sx, self._sy = None, None
-+
-+ def get_array(self):
-+ """Override to return the full-resolution array"""
-+ return self._full_res
-+
-+ def _scale_to_res(self):
-+ """ Change self._A and _extent to render an image whose
-+ resolution is matched to the eventual rendering."""
-+
-+ ax = self.axes
-+ ext = ax.transAxes.transform([1, 1]) - ax.transAxes.transform([0, 0])
-+ xlim, ylim = ax.get_xlim(), ax.get_ylim()
-+ dx, dy = xlim[1] - xlim[0], ylim[1] - ylim[0]
-+
-+ y0 = max(0, ylim[0] - 5)
-+ y1 = min(self._full_res.shape[0], ylim[1] + 5)
-+ x0 = max(0, xlim[0] - 5)
-+ x1 = min(self._full_res.shape[1], xlim[1] + 5)
-+ y0, y1, x0, x1 = map(int, [y0, y1, x0, x1])
-+
-+ sy = int(max(1, min((y1 - y0) / 5., np.ceil(dy / ext[1]))))
-+ sx = int(max(1, min((x1 - x0) / 5., np.ceil(dx / ext[0]))))
-+
-+ # have we already calculated what we need?
-+ if sx == self._sx and sy == self._sy and \
-+ x0 == self._bounds[0] and x1 == self._bounds[1] and \
-+ y0 == self._bounds[2] and y1 == self._bounds[3]:
-+ return
-+
-+ self._A = self._full_res[y0:y1:sy, x0:x1:sx]
-+ self._A = cbook.safe_masked_invalid(self._A)
-+ x1 = x0 + self._A.shape[1] * sx
-+ y1 = y0 + self._A.shape[0] * sy
-+
-+ self.set_extent([x0 - .5, x1 - .5, y0 - .5, y1 - .5])
-+ self._sx = sx
-+ self._sy = sy
-+ self._bounds = (x0, x1, y0, y1)
-+ self.changed()
-+
-+ def draw(self, renderer, *args, **kwargs):
-+ self._scale_to_res()
-+ super(ModestImage, self).draw(renderer, *args, **kwargs)
-+
-+
-+def imshow(axes, X, cmap=None, norm=None, aspect=None,
-+ interpolation=None, alpha=None, vmin=None, vmax=None,
-+ origin=None, extent=None, shape=None, filternorm=1,
-+ filterrad=4.0, imlim=None, resample=None, url=None, **kwargs):
-+ """Similar to matplotlib's imshow command, but produces a ModestImage
-+
-+ Unlike matplotlib version, must explicitly specify axes
-+ Author: Chris Beaumont <beaumont at hawaii.edu>
-+ """
-+
-+ if not axes._hold:
-+ axes.cla()
-+ if norm is not None:
-+ assert(isinstance(norm, mcolors.Normalize))
-+ if aspect is None:
-+ aspect = rcParams['image.aspect']
-+ axes.set_aspect(aspect)
-+ im = ModestImage(axes, cmap, norm, interpolation, origin, extent,
-+ filternorm=filternorm,
-+ filterrad=filterrad, resample=resample, **kwargs)
-+
-+ im.set_data(X)
-+ im.set_alpha(alpha)
-+ axes._set_artist_props(im)
-+
-+ if im.get_clip_path() is None:
-+ # image does not already have clipping set, clip to axes patch
-+ im.set_clip_path(axes.patch)
-+
-+ #if norm is None and shape is None:
-+ # im.set_clim(vmin, vmax)
-+ if vmin is not None or vmax is not None:
-+ im.set_clim(vmin, vmax)
-+ else:
-+ im.autoscale_None()
-+ im.set_url(url)
-+
-+ # update ax.dataLim, and, if autoscaling, set viewLim
-+ # to tightly fit the image, regardless of dataLim.
-+ im.set_extent(im.get_extent())
-+
-+ axes.images.append(im)
-+ im._remove_method = lambda h: axes.images.remove(h)
-+
-+ return im
-+
-+def MergeImg(scatts, styles, output_queue):
-+
-+ init = True
-+ merge_tmp = grass.tempfile()
-+ for cat_id, scatt in scatts.iteritems():
-+ #print "color map %d" % cat_id
-+ #TODO make more general
-+ if cat_id != 0 and (styles[cat_id]['opacity'] == 0.0 or \
-+ not styles[cat_id]['show']):
-+ continue
-+ if init:
-+ b1_i = scatt['bands_info']['b1']
-+ b2_i = scatt['bands_info']['b2']
-+
-+ full_extend = (b1_i['min'] - 0.5, b1_i['max'] + 0.5, b2_i['min'] - 0.5, b2_i['max'] + 0.5)
-+
-+ if cat_id == 0:
-+ cmap = matplotlib.cm.jet
-+ cmap.set_bad('w',1.)
-+ cmap._init()
-+ cmap._lut[len(cmap._lut) - 1, -1] = 0
-+ else:
-+ colors = styles[cat_id]['color'].split(":")
-+
-+ cmap = matplotlib.cm.jet
-+ cmap.set_bad('w',1.)
-+ cmap._init()
-+ cmap._lut[len(cmap._lut) - 1, -1] = 0
-+ cmap._lut[:, 0] = int(colors[0])/255.0
-+ cmap._lut[:, 1] = int(colors[1])/255.0
-+ cmap._lut[:, 2] = int(colors[2])/255.0
-+
-+ #if init:
-+ masked_cat = np.ma.masked_less_equal(scatt['np_vals'], 0)
-+
-+ vmax = np.amax(masked_cat)
-+ masked_cat = masked_cat / float(vmax)
-+
-+ colored_cat = np.uint8(cmap(masked_cat) * 255)
-+ del masked_cat
-+ del cmap
-+
-+ #colored_cat[...,3] = np.choose(masked_cat.mask, (255, 0))
-+
-+ if init:
-+ merged_img = np.memmap(merge_tmp, dtype='uint8', mode='w+', shape=colored_cat.shape)
-+ merged_img[:] = colored_cat[:]
-+ init = False
-+ else:
-+ #c_img_a = np.memmap(grass.tempfile(), dtype="uint16", mode='w+', shape = shape)
-+ c_img_a = colored_cat.astype('uint16')[:,:,3] * styles[cat_id]['opacity']
-+
-+ #TODO apply strides and there will be no need for loop
-+ #b = as_strided(a, strides=(0, a.strides[3], a.strides[3], a.strides[3]), shape=(3, a.shape[0], a.shape[1]))
-+
-+ for i in range(3):
-+ merged_img[:,:,i] = (merged_img[:,:,i] * (255 - c_img_a) + colored_cat[:,:,i] * c_img_a) / 255;
-+ merged_img[:,:,3] = (merged_img[:,:,3] * (255 - c_img_a) + 255 * c_img_a) / 255;
-+
-+ del c_img_a
-+
-+ del colored_cat
-+ output_queue.put((merged_img, full_extend))
-+ #return merged_img, full_extend
-\ No newline at end of file
-Index: gui/wxpython/scatt_plot/controllers.py
-===================================================================
---- gui/wxpython/scatt_plot/controllers.py (revision 0)
-+++ gui/wxpython/scatt_plot/controllers.py (working copy)
-@@ -0,0 +1,666 @@
-+"""!
-+ at package scatt_plot.controllers
-+
-+ at brief Controller layer for scatter plot tool.
-+
-+Classes:
-+
-+(C) 2013 by the GRASS Development Team
-+
-+This program is free software under the GNU General Public License
-+(>=v2). Read the file COPYING that comes with GRASS for details.
-+
-+ at author Stepan Turek <stepan.turek seznam.cz> (mentor: Martin Landa)
-+"""
-+import os
-+import sys
-+
-+#TODO just for testing
-+import time
-+
-+#TODO
-+import wx
-+
-+from core.gcmd import GException, GError, GMessage, RunCommand
-+
-+from scatt_plot.scatt_core import Core, idBandsToidScatt
-+
-+from scatt_plot.dialogs import AddScattPlotDialog
-+from scatt_plot.gthreading import gThread
-+from core.gconsole import EVT_CMD_DONE
-+from grass.pydispatch.signal import Signal
-+
-+class ScattsManager(wx.EvtHandler):
-+ def __init__(self, guiparent, giface, iclass_mapwin):
-+ #TODO remove iclass parameter
-+
-+ wx.EvtHandler.__init__(self)
-+ self.giface = giface
-+ self.mapDisp = giface.GetMapDisplay()
-+
-+ if iclass_mapwin:
-+ self.mapWin = iclass_mapwin
-+ else:
-+ self.mapWin = giface.GetMapWindow()
-+
-+ self.guiparent = guiparent
-+
-+ self.show_add_scatt_plot = False
-+
-+ self.core = Core()
-+ self.scatts_dt, self.scatt_conds_dt = self.core.GetScattsData()
-+
-+ self.cats_mgr = CategoriesManager(self, self.core)
-+
-+ self.thread = gThread(self);
-+
-+ self.plots = {}
-+ self.added_cats_rasts = {}
-+
-+ self.cats_to_update = []
-+
-+ self.plot_mode = None
-+ self.pol_sel_mode = [False, None]
-+
-+ self.data_set = False
-+
-+ if iclass_mapwin:
-+ self.mapWin_conn = MapWinConnection(self, self.mapWin, self.core.CatRastUpdater())
-+ self.iclass_conn = IClassConnection(self, iclass_mapwin.parent, self.cats_mgr)
-+ else:
-+ self.mapWin_conn = None
-+ self.iclass_conn = None
-+
-+ self.tasks_pids = {
-+ 'add_scatt' : [],
-+ 'set_data' : -1,
-+ 'set_data_add' : -1,
-+ 'set_edit_cat_data' : -1,
-+ 'mapwin_conn' : [],
-+ 'render_plots' : -1,
-+ 'render' : []
-+ }
-+
-+ self.Bind(EVT_CMD_DONE, self.OnThreadDone)
-+
-+ def CleanUp(self):
-+ self.core.CleanUp()
-+
-+ for scatt_id, scatt in self.plots.items():
-+ scatt.CleanUp()
-+
-+ def OnThreadDone(self, event):
-+
-+ if event.exception:
-+ GError(str(event.exception))
-+ return
-+
-+ if event.pid in self.tasks_pids['mapwin_conn']:
-+ self.tasks_pids['mapwin_conn'].remove(event.pid)
-+ updated_cats = event.ret
-+
-+ for cat in updated_cats:
-+ if cat not in self.cats_to_update:
-+ self.cats_to_update.append(cat)
-+
-+ if not self.tasks_pids['mapwin_conn']:
-+ self.tasks_pids['render_plots'] = self.thread.GetId()
-+ self.thread.Run(callable = self.core.ComputeCatsScatts,
-+ cats_ids = self.cats_to_update[:])
-+ del self.cats_to_update[:]
-+
-+ return
-+
-+ if self.tasks_pids['render_plots'] == event.pid:
-+ self.RenderScattPlts()
-+ return
-+
-+ if event.pid in self.tasks_pids['render']:
-+ self.tasks_pids['render'].remove(event.pid)
-+ return
-+
-+ if event.pid in self.tasks_pids['add_scatt']:
-+ self.tasks_pids['add_scatt'].remove(event.pid)
-+ self.AddScattPlotDone(event)
-+ return
-+
-+ if self.tasks_pids['set_data'] == event.pid:
-+ self.SetDataDone(event)
-+ return
-+
-+ if self.tasks_pids['set_data_add'] == event.pid:
-+ self.SetDataDone(event)
-+ self.AddScattPlot()
-+ return
-+
-+ if self.tasks_pids['set_edit_cat_data'] == event.pid:
-+ self.SetEditCatDataDone(event)
-+ return
-+
-+ def SetData(self, bands):
-+
-+ self.CleanUp()
-+
-+ self.data_set = False
-+
-+ if self.show_add_scatt_plot:
-+ self.tasks_pids['set_data_add'] = self.thread.GetId()
-+ else:
-+ self.tasks_pids['set_data'] = self.thread.GetId()
-+
-+ self.thread.Run(callable = self.core.SetData, bands = bands)
-+
-+ def SetDataDone(self, event):
-+
-+ self.data_set = True
-+ self.cats_mgr.InitCoreCats()
-+
-+ def OnOutput(self, event):
-+ """!Print thread output according to debug level.
-+ """
-+ print event.text
-+
-+ def GetBands(self):
-+ return self.core.GetBands()
-+
-+ def AddScattPlot(self):
-+ if not self.data_set and self.iclass_conn:
-+ self.show_add_scatt_plot = True
-+ self.iclass_conn.SetData()
-+ self.show_add_scatt_plot = False
-+ return
-+ if not self.data_set:
-+ GError(_('No data set.'))
-+ return
-+
-+ dlg = AddScattPlotDialog(parent = self.guiparent, bands = self.core.GetBands())
-+ if dlg.ShowModal() == wx.ID_OK:
-+ self._addScattPlot(dlg.GetScattId())
-+
-+ dlg.Destroy()
-+
-+ return
-+
-+ #for testing
-+ jj = len(self.core.GetBands()) - 1
-+
-+ max_k = 20
-+ k = 0
-+ for i in range(jj):
-+ for j in range(i, jj):
-+ print i * (len(self.core.GetBands()) - 1) + j
-+ self._addScattPlot(i * (len(self.core.GetBands()) - 1) + j)
-+ k += 1
-+ if k == max_k:
-+ break
-+ if k == max_k:
-+ break
-+ jj -= 1
-+
-+
-+ def _addScattPlot(self, scatt_id):
-+ if self.plots.has_key(scatt_id):
-+ GMessage(_("Scatter plot has been already added."))
-+ return
-+
-+ self.tasks_pids['add_scatt'].append(self.thread.GetId())
-+
-+ self.thread.Run(callable = self.core.AddScattPlot, scatt_id = scatt_id)
-+
-+ def RenderScattPlts(self, scatt_ids = None):
-+ if len(self.tasks_pids['render']) > 1:
-+ print "skip"
-+ return
-+
-+ self.tasks_pids['render'].append(self.thread.GetId())
-+ self.thread.Run(callable = self._renderscattplts, scatt_ids = scatt_ids)
-+
-+ def _renderscattplts(self, scatt_ids):
-+ cats_attrs = self.cats_mgr.GetCategoriesAttrs()
-+ for i_scatt_id, scatt in self.plots.items():
-+ if scatt_ids is not None and i_scatt_id not in scatt_ids:
-+ continue
-+
-+ scatt_dt = self.scatts_dt.GetScatt(i_scatt_id)
-+ ellipses_dt = self.scatts_dt.GetEllipses(i_scatt_id)
-+
-+ if self.pol_sel_mode[0]:
-+ self._getSelectedAreas(i_scatt_id, scatt_dt, cats_attrs)
-+
-+ scatt.Plot(scatts = scatt_dt, ellipses = ellipses_dt, styles = cats_attrs)
-+
-+ def _getSelectedAreas(self, scatt_id, scatt_dt, cats_attrs):
-+
-+ cat_id = self.cats_mgr.GetSelectedCat()
-+ if not cat_id:
-+ return
-+
-+ sel_a_cat_id = max(cats_attrs.keys()) + 1
-+
-+ s = self.scatt_conds_dt.GetScatt(scatt_id, [cat_id])
-+ if not s:
-+ return
-+ cats_attrs[sel_a_cat_id] = {'color' : "255:255:0",
-+ 'opacity' : 0.7,
-+ 'show' : True}
-+
-+ scatt_dt[sel_a_cat_id] = s[cat_id]
-+
-+ def AddScattPlotDone(self, event):
-+
-+ scatt_id = event.kwds['scatt_id']
-+
-+ #TODO guiparent - not very good
-+ self.plots[scatt_id] = self.guiparent.NewScatterPlot(scatt_id = scatt_id)
-+ self.plots[scatt_id].plotClosed.connect(self.PlotClosed)
-+
-+ if self.plot_mode:
-+ self.plots[scatt_id].SetMode(self.plot_mode)
-+
-+ self.RenderScattPlts(scatt_ids = [scatt_id])
-+
-+ def PlotClosed(self, scatt_id):
-+ del self.plots[scatt_id]
-+
-+ def SetPlotsMode(self, mode):
-+
-+ self.plot_mode = mode
-+ for scatt in self.plots.itervalues():
-+ scatt.SetMode(mode)
-+ scatt.SetEditingMode(self.pol_sel_mode[1])
-+
-+ def ActivateSelectionPolygonMode(self, activate):
-+ self.pol_sel_mode[0] = activate
-+ self.RenderScattPlts()
-+ return activate
-+
-+ def ProcessSelectionPolygons(self, process_mode):
-+ scatts_polygons = {}
-+ for scatt_id, scatt in self.plots.iteritems():
-+ coords = scatt.GetCoords()
-+ if coords is not None:
-+ scatts_polygons[scatt_id] = coords
-+
-+ if not scatts_polygons:
-+ return
-+
-+ value = 1
-+ if process_mode == 'remove':
-+ value = 0
-+
-+ sel_cat_id = self.cats_mgr.GetSelectedCat()
-+ if not sel_cat_id:
-+ dlg = wx.MessageDialog(parent = self.guiparent,
-+ message = _("In order to select arrea in scatter plot, "
-+ "you have to select class first.\n\n"
-+ "There is no class yet, "
-+ "do you want to create one?"),
-+ caption = _("No class selected"),
-+ style = wx.YES_NO)
-+ if dlg.ShowModal() == wx.ID_YES:
-+ self.iclass_conn.EmptyCategories()
-+
-+ sel_cat_id = self.cats_mgr.GetSelectedCat()
-+ if not sel_cat_id:
-+ return
-+
-+ for scatt in self.plots.itervalues():
-+ scatt.SetEmpty()
-+
-+ self.tasks_pids['set_edit_cat_data'] = self.thread.GetId()
-+ self.thread.Run(callable = self.core.UpdateCategoryWithPolygons,
-+ cat_id = sel_cat_id,
-+ scatts_pols = scatts_polygons,
-+ value = value)
-+
-+ def SetPlotsEditingMode(self, mode):
-+
-+ self.pol_sel_mode[1] = mode
-+ for scatt in self.plots.itervalues():
-+ scatt.SetEditingMode(mode)
-+
-+ def SetEditCatDataDone(self, event):
-+
-+ if event.exception:
-+ GError(_("Error occured during computation of scatter plot category:\n%s"),
-+ parent = self.guiparent, showTraceback = False)
-+
-+ cat_id = event.ret
-+
-+ self.RenderScattPlts()
-+
-+ cat_id = event.kwds["cat_id"]
-+
-+ cat_rast = self.core.GetCatRast(cat_id)
-+
-+ if cat_rast not in self.added_cats_rasts.values():
-+
-+ cats_attrs = self.cats_mgr.GetCategoryAttrs(cat_id)
-+
-+
-+ region = self.core.GetRegion()
-+ ret, err_msg = RunCommand('r.region',
-+ map = cat_rast,
-+ getErrorMsg = True,
-+ n = "%f" % region['n'],
-+ s = "%f" % region['s'],
-+ e = "%f" % region['e'],
-+ w = "%f" % region['w'],
-+ )
-+
-+ ret, err_msg = RunCommand('r.colors',
-+ map = cat_rast,
-+ rules = "-",
-+ stdin = "1 %s" % cats_attrs["color"],
-+ getErrorMsg = True)
-+
-+ if ret != 0:
-+ GError(_("r.region failed\n%s" % err_msg))
-+
-+ self.mapWin.Map.AddLayer(ltype = "raster", name = "cat_%d" % cat_id, render = True,
-+ command = ["d.rast", "map=%s" % cat_rast, "values=1"])
-+
-+
-+ if ret != 0:
-+ GError(_("r.region failed\n%s" % err_msg))
-+
-+ self.added_cats_rasts[cat_id] = cat_rast
-+
-+ self.giface.updateMap.emit()
-+
-+ def DigitDataChanged(self, vectMap, digit):
-+
-+ if self.mapWin_conn:
-+ self.mapWin_conn.DigitDataChanged(vectMap, digit)
-+ return 1
-+ else:
-+ return 0
-+
-+ def GetCategoriesManager(self):
-+ return self.cats_mgr
-+
-+
-+class CategoriesManager:
-+
-+ def __init__(self, scatt_mgr, core):
-+
-+ self.core = core
-+ self.scatt_mgr = scatt_mgr
-+
-+ self.cats = {}
-+ self.cats_ids = []
-+
-+ self.sel_cat_id = None
-+
-+ self.initialized = Signal('CategoriesManager.initialized')
-+ self.setCategoryAttrs = Signal('CategoriesManager.setCategoryAttrs')
-+ self.deletedCategory = Signal('CategoriesManager.deletedCategory')
-+ self.addedCategory = Signal('CategoriesManager.addedCategory')
-+
-+ def Clear(self):
-+
-+ self.cats.clear()
-+ del self.cats_ids[:]
-+
-+ self.sel_cat_id = None
-+
-+ def InitCoreCats(self):
-+ if self.scatt_mgr.data_set:
-+ for cat_id in self.cats_ids:
-+ self.core.AddCategory(cat_id)
-+
-+ def AddCategory(self, cat_id = None, name = None, color = None):
-+
-+ if cat_id is None:
-+ if self.cats_ids:
-+ cat_id = max(self.cats_ids) + 1
-+ else:
-+ cat_id = 1
-+
-+ if self.scatt_mgr.data_set:
-+ self.scatt_mgr.thread.Run(callable = self.core.AddCategory,
-+ cat_id = cat_id)
-+ #TODO check number of cats
-+ #if ret < 0: #TODO
-+ # return -1;
-+
-+ self.cats[cat_id] = {
-+ 'name' : _('Category %s' % cat_id ),
-+ 'color' : "0:0:0",
-+ 'opacity' : 1.0,
-+ 'show' : True
-+ }
-+
-+ self.cats_ids.append(cat_id)
-+
-+ if name is not None:
-+ self.cats[cat_id]["name"] = name
-+
-+ if color is not None:
-+ self.cats[cat_id]["color"] = color
-+
-+ self.addedCategory.emit(cat_id = cat_id,
-+ name = self.cats[cat_id]["name"],
-+ color = self.cats[cat_id]["color"] )
-+ return cat_id
-+
-+ def SetCategoryAttrs(self, cat_id, attrs_dict):
-+ render = False
-+ for k, v in attrs_dict.iteritems():
-+ if not render and k in ['name', 'color', 'opacity', 'show']:
-+ render = True
-+
-+ self.cats[cat_id][k] = v
-+
-+ #TODO optimization
-+ if render:
-+ self.scatt_mgr.RenderScattPlts()
-+
-+ self.setCategoryAttrs.emit(cat_id = cat_id, attrs_dict = attrs_dict)
-+
-+ def DeleteCategory(self, cat_id):
-+
-+ if self.scatt_mgr.data_set:
-+ self.scatt_mgr.thread.Run(callable = self.core.DeleteCategory,
-+ cat_id = cat_id)
-+ del self.cats[cat_id]
-+ self.cats_ids.remove(cat_id)
-+
-+ self.deletedCategory.emit(cat_id = cat_id)
-+
-+ #TODO emit event?
-+ def SetSelectedCat(self, cat_id):
-+ self.sel_cat_id = cat_id
-+ if self.scatt_mgr.pol_sel_mode[0]:
-+ self.scatt_mgr.RenderScattPlts()
-+
-+ def GetSelectedCat(self):
-+ return self.sel_cat_id
-+
-+ def GetCategoryAttrs(self, cat_id):
-+ #TODO is mutable
-+ return self.cats[cat_id]
-+
-+ def GetCategoriesAttrs(self):
-+ #TODO is mutable
-+ return self.cats
-+
-+ def GetCategories(self):
-+ return self.cats_ids[:]
-+
-+ def SetCategoryPosition(self):
-+ if newindex > oldindex:
-+ newindex -= 1
-+
-+ self.cats_ids.insert(newindex, self.cats_ids.pop(oldindex))
-+
-+class MapWinConnection:
-+ def __init__(self, scatt_mgr, mapWin, scatt_rast_updater):
-+ self.mapWin = mapWin
-+ self.vectMap = None
-+ self.scatt_rast_updater = scatt_rast_updater
-+ self.scatt_mgr = scatt_mgr
-+ self.cats_mgr = scatt_mgr.cats_mgr
-+
-+ self.thread = self.scatt_mgr.thread
-+
-+ #TODO
-+ self.mapWin.parent.toolbars["vdigit"].editingStarted.connect(self.DigitDataChanged)
-+
-+ #def ChangeMap(self, vectMap, layers_cats):
-+ # self.vectMap = vectMap
-+ # self.layers_cats = layers_cats
-+
-+ #ret, region, msg = RunCommand("v.to.rast",
-+ # flags = "gp",
-+ # getErrorMsg = True,
-+ # read = True)
-+
-+ def _connectSignals(self):
-+ self.digit.featureAdded.connect(self.AddFeature)
-+ self.digit.areasDeleted.connect(self.DeleteAreas)
-+ self.digit.featuresDeleted.connect(self.DeleteAreas)
-+ self.digit.vertexMoved.connect(self.EditedFeature)
-+ self.digit.vertexRemoved.connect(self.EditedFeature)
-+ self.digit.lineEdited.connect(self.EditedFeature)
-+ self.digit.featuresMoved.connect(self.EditedFeature)
-+
-+ def AddFeature(self, new_bboxs, new_areas_cats):
-+ if not self.scatt_mgr.data_set:
-+ return
-+
-+ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
-+ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
-+ new_bboxs = new_bboxs,
-+ old_bboxs = [],
-+ old_areas_cats = [],
-+ new_areas_cats = new_areas_cats)
-+
-+ def DeleteAreas(self, old_bboxs, old_areas_cats):
-+ if not self.scatt_mgr.data_set:
-+ return
-+
-+ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
-+ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
-+ new_bboxs = [],
-+ old_bboxs = old_bboxs,
-+ old_areas_cats = old_areas_cats,
-+ new_areas_cats = [])
-+
-+
-+ def EditedFeature(self, new_bboxs, new_areas_cats, old_bboxs, old_areas_cats):
-+ if not self.scatt_mgr.data_set:
-+ return
-+
-+ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
-+ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
-+ new_bboxs = new_bboxs,
-+ old_bboxs = old_bboxs,
-+ old_areas_cats = old_areas_cats,
-+ new_areas_cats = new_areas_cats)
-+
-+ def DigitDataChanged(self, vectMap, digit):
-+
-+ self.digit = digit
-+ self.vectMap = vectMap
-+
-+ self.digit.EmitSignals(emit = True)
-+
-+ self.scatt_rast_updater.SetVectMap(vectMap)
-+
-+ self._connectSignals()
-+
-+
-+class IClassConnection:
-+ def __init__(self, scatt_mgr, iclass_frame, cats_mgr):
-+ self.iclass_frame = iclass_frame
-+ self.stats_data = self.iclass_frame.stats_data
-+ self.cats_mgr = cats_mgr
-+ self.scatt_mgr = scatt_mgr
-+
-+ self.stats_data.statisticsAdded.connect(self.AddCategory)
-+ self.stats_data.statisticsDeleted.connect(self.DeleteCategory)
-+ self.stats_data.allStatisticsDeleted.connect(self.DeletAllCategories)
-+ self.stats_data.statisticsSet.connect(self.SetCategory)
-+
-+ self.iclass_frame.groupSet.connect(self.GroupSet)
-+
-+ self.cats_mgr.setCategoryAttrs.connect(self.SetStatistics)
-+ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
-+ self.cats_mgr.addedCategory.connect(self.AddStatistics)
-+
-+ self.iclass_frame.categoryChanged.connect(self.CategoryChanged)
-+
-+ self.SyncCats()
-+
-+ def SetData(self):
-+ self.iclass_frame.AddBands()
-+
-+ def EmptyCategories(self):
-+ self.iclass_frame.OnCategoryManager(None)
-+
-+ def SyncCats(self):
-+ self.cats_mgr.addedCategory.disconnect(self.AddStatistics)
-+ cats = self.stats_data.GetCategories()
-+ for c in cats:
-+ stats = self.stats_data.GetStatistics(c)
-+ self.cats_mgr.AddCategory(c, stats.name, stats.color)
-+ self.cats_mgr.addedCategory.connect(self.AddStatistics)
-+
-+ def CategoryChanged(self, cat):
-+ self.cats_mgr.SetSelectedCat(cat)
-+
-+ def AddCategory(self, cat, name, color):
-+ self.cats_mgr.addedCategory.disconnect(self.AddStatistics)
-+ self.cats_mgr.AddCategory(cat_id = cat, name = name, color = color)
-+ self.cats_mgr.addedCategory.connect(self.AddStatistics)
-+
-+ def DeleteCategory(self, cat):
-+ self.cats_mgr.deletedCategory.disconnect(self.DeleteStatistics)
-+ self.cats_mgr.DeleteCategory(cat)
-+ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
-+
-+ def DeletAllCategories(self):
-+
-+ self.cats_mgr.deletedCategory.disconnect(self.DeleteStatistics)
-+ cats = self.stats_data.GetCategories()
-+ for c in cats:
-+ self.cats_mgr.DeleteCategory(c)
-+ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
-+
-+ def SetCategory(self, cat, stats):
-+
-+ self.cats_mgr.setCategoryAttrs.disconnect(self.SetStatistics)
-+ cats_attr = {}
-+ for attr in ['name', 'color']:
-+ if stats.has_key(attr):
-+ cats_attr[attr] = stats[attr]
-+
-+ if cats_attr:
-+ self.cats_mgr.SetCategoryAttrs(cat, cats_attr)
-+ self.cats_mgr.setCategoryAttrs.connect(self.SetStatistics)
-+
-+
-+ def SetStatistics(self, cat_id, attrs_dict):
-+ self.stats_data.statisticsSet.disconnect(self.SetCategory)
-+ self.stats_data.GetStatistics(cat_id).SetStatistics(attrs_dict)
-+ self.stats_data.statisticsSet.connect(self.SetCategory)
-+
-+ def AddStatistics(self, cat_id, name, color):
-+ self.stats_data.statisticsAdded.disconnect(self.AddCategory)
-+ self.stats_data.AddStatistics(cat_id, name, color)
-+ self.stats_data.statisticsAdded.connect(self.AddCategory)
-+
-+ def DeleteStatistics(self, cat_id):
-+ self.stats_data.statisticsDeleted.disconnect(self.DeleteCategory)
-+ self.stats_data.DeleteStatistics(cat_id)
-+ self.stats_data.statisticsDeleted.connect(self.DeleteCategory)
-+
-+ def GroupSet(self, group):
-+ res = RunCommand('i.group',
-+ flags = 'g',
-+ group = group, subgroup = group,
-+ read = True).strip()
-+ if res.split('\n')[0]:
-+ bands = res.split('\n')
-+ self.scatt_mgr.SetData(bands)
-\ No newline at end of file
+ def GetSelectedIndices(self, state = wx.LIST_STATE_SELECTED):
Index: gui/wxpython/scatt_plot/gthreading.py
===================================================================
--- gui/wxpython/scatt_plot/gthreading.py (revision 0)
+++ gui/wxpython/scatt_plot/gthreading.py (working copy)
-@@ -0,0 +1,133 @@
+@@ -0,0 +1,134 @@
+"""!
+ at package scatt_plot.gthreading
+
@@ -2477,7 +1719,7 @@
+ os.environ['GRASS_MESSAGE_FORMAT'] = 'gui'
+ while True:
+ requestId, args, kwds = self.requestQ.get()
-+ for key in ('callable', 'onDone', 'onPrepare', 'userData'):
++ for key in ('callable', 'onDone', 'onPrepare', 'userdata'):
+ if key in kwds:
+ vars()[key] = kwds[key]
+ del kwds[key]
@@ -2521,10 +1763,11 @@
+ time.sleep(.1)
+
+ if self.receiver:
-+ event = wxCmdDone(kwds = kwds,
-+ args = args, #TODO expand args to kwds
++ event = wxCmdDone(kwds=kwds,
++ args=args, #TODO expand args to kwds
+ ret=ret,
+ exception=exception,
++ userdata=vars()['userdata'],
+ pid=requestId)
+
+ # send event
@@ -2541,7 +1784,7 @@
===================================================================
--- gui/wxpython/scatt_plot/dialogs.py (revision 0)
+++ gui/wxpython/scatt_plot/dialogs.py (working copy)
-@@ -0,0 +1,142 @@
+@@ -0,0 +1,149 @@
+"""!
+ at package scatt_plot.dialogs
+
@@ -2562,6 +1805,8 @@
+import wx
+from scatt_plot.scatt_core import idBandsToidScatt
+
++from core.gcmd import GError
++
+class AddScattPlotDialog(wx.Dialog):
+
+ def __init__(self, parent, bands, id = wx.ID_ANY):
@@ -2570,7 +1815,8 @@
+
+ self.bands = bands
+
-+ self.scatt_id = None
++ self.x_band = None
++ self.y_band = None
+
+ self._createWidgets()
+
@@ -2579,13 +1825,13 @@
+ self.labels = {}
+ self.params = {}
+
-+ self.band_1_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("Band 1:"))
++ self.band_1_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("%s axis:" % "x"))
+
+ self.band_1_ch = wx.ComboBox(parent = self, id = wx.ID_ANY,
+ choices = self.bands,
+ style = wx.CB_READONLY, size = (350, 30))
+
-+ self.band_2_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("Band 2:"))
++ self.band_2_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("%s axis:" % "y"))
+
+ self.band_2_ch = wx.ComboBox(parent = self, id = wx.ID_ANY,
+ choices = self.bands,
@@ -2664,31 +1910,35 @@
+ def OnOk(self, event):
+ """!
+ """
-+ band_1 = self.band_1_ch.GetSelection()
-+ band_2 = self.band_2_ch.GetSelection()
++ b_x = self.band_1_ch.GetSelection()
++ b_y = self.band_2_ch.GetSelection()
+
++ err = True
+
-+ if band_1 == band_2:
++ if b_x < 0 or b_y < 0:
++ GError(_("Select both x and y bands."))
++ elif b_y == b_x:
+ GError(_("Selected bands must be different."))
++ else:
++ err = False
++
++ if err:
++ self.band_y = None
++ self.band_x = None
+ return
-+
-+ #TODO axes selection
-+ if band_1 > band_2:
-+ tmp_band = band_2
-+ band_2 = band_1
-+ band_1 = band_2
+
-+ self.scatt_id = idBandsToidScatt(band_1, band_2, len(self.bands))
++ self.band_y = b_y
++ self.band_x = b_x
+
+ event.Skip()
+
-+ def GetScattId(self):
-+ return self.scatt_id
++ def GetBands(self):
++ return (self.band_x, self.band_y)
Index: gui/wxpython/scatt_plot/toolbars.py
===================================================================
--- gui/wxpython/scatt_plot/toolbars.py (revision 0)
+++ gui/wxpython/scatt_plot/toolbars.py (working copy)
-@@ -0,0 +1,178 @@
+@@ -0,0 +1,212 @@
+"""!
+ at package scatt_plot.toolbars
+
@@ -2724,6 +1974,7 @@
+
+ # realize the toolbar
+ self.Realize()
++ self.scatt_mgr.modeSet.connect(self.ModeSet)
+
+ def _toolbarData(self):
+
@@ -2733,7 +1984,7 @@
+ label = _('Show manual')),
+ 'add_scatt_pl' : MetaIcon(img = 'layer-raster-analyze',
+ label = _('Add scatter plot')),
-+ 'selCatPol' : MetaIcon(img = 'polygon-create',
++ 'selCatPol' : MetaIcon(img = 'polygon',
+ label = _('Select area with polygon')),
+ 'pan' : MetaIcon(img = 'pan',
+ label = _('Pan'),
@@ -2746,6 +1997,9 @@
+ }
+
+ return self._getToolbarData((
++ ('add_scatt', icons["add_scatt_pl"],
++ lambda event : self.scatt_mgr.AddScattPlot()),
++ (None, ),
+ ("cats_mgr", icons['cats_mgr'],
+ lambda event: self.parent.ShowCategoryPanel(event.Checked()), wx.ITEM_CHECK),
+ (None, ),
@@ -2757,10 +2011,7 @@
+ wx.ITEM_CHECK),
+ ('sel_pol_mode', icons['selCatPol'],
+ self.ActivateSelectionPolygonMode,
-+ wx.ITEM_CHECK),
-+ (None, ),
-+ ('add_scatt', icons["add_scatt_pl"],
-+ lambda event : self.scatt_mgr.AddScattPlot())
++ wx.ITEM_CHECK)
+ #('settings', icon["settings"],
+ # self.parent.OnSettings),
+ #('help', icons["help"],
@@ -2768,15 +2019,14 @@
+ ))
+
+ def GetToolId(self, toolName): #TODO can be useful in base
-+
+ return vars(self)[toolName]
+
+ def SetPloltsMode(self, event, tool_name):
-+
++ self.scatt_mgr.modeSet.disconnect(self.ModeSet)
+ if event.Checked() == True:
+ for i_tool_data in self._data:
+ i_tool_name = i_tool_data[0]
-+ if not i_tool_name or i_tool_name == "cats_mgr":
++ if not i_tool_name or i_tool_name in ["cats_mgr", "sel_pol_mode"]:
+ continue
+ if i_tool_name == tool_name:
+ continue
@@ -2786,6 +2036,7 @@
+ self.scatt_mgr.SetPlotsMode(mode = tool_name)
+ else:
+ self.scatt_mgr.SetPlotsMode(mode = None)
++ self.scatt_mgr.modeSet.connect(self.ModeSet)
+
+ def ActivateSelectionPolygonMode(self, event):
+
@@ -2795,6 +2046,17 @@
+ i_tool_id = vars(self)['sel_pol_mode']
+ self.ToggleTool(i_tool_id, activated)
+
++ def ModeSet(self, mode):
++ self.UnsetMode()
++
++ def UnsetMode(self):
++ for i_tool_data in self._data:
++ i_tool_name = i_tool_data[0]
++ if not i_tool_name or i_tool_name in ["cats_mgr", "sel_pol_mode"]:
++ continue
++ i_tool_id = vars(self)[i_tool_name]
++ self.ToggleTool(i_tool_id, False)
++
+class EditingToolbar(BaseToolbar):
+ """!Main toolbar
+ """
@@ -2806,6 +2068,7 @@
+
+ # realize the toolbar
+ self.Realize()
++ self.scatt_mgr.modeSet.connect(self.ModeSet)
+
+ def _toolbarData(self):
+ """!Toolbar data
@@ -2820,7 +2083,7 @@
+ 'addVertex' : MetaIcon(img = 'vertex-create',
+ label = _('Add new vertex'),
+ desc = _('Add new vertex to polygon boundary scatter plot')),
-+ 'editLine' : MetaIcon(img = 'line-edit',
++ 'editLine' : MetaIcon(img = 'polygon-create',
+ label = _('Edit line/boundary'),
+ desc = _('Add new vertex between last and first points of the boundary')),
+ 'moveVertex' : MetaIcon(img = 'vertex-move',
@@ -2829,6 +2092,9 @@
+ 'removeVertex' : MetaIcon(img = 'vertex-delete',
+ label = _('Remove vertex'),
+ desc = _('Remove boundary vertex.')),
++ 'delete' : MetaIcon(img = 'polygon-delete',
++ label = _('Edit line/boundary'),
++ desc = _('Delete polygon')),
+ }
+
+ return self._getToolbarData((
@@ -2848,10 +2114,14 @@
+ wx.ITEM_CHECK),
+ ('delete_vertex', self.icons['removeVertex'],
+ lambda event: self.SetMode(event, 'delete_vertex'),
++ wx.ITEM_CHECK),
++ ('delete', self.icons['delete'],
++ lambda event: self.SetMode(event, 'delete'),
+ wx.ITEM_CHECK)
+ ))
+
+ def SetMode(self, event, tool_name):
++ self.scatt_mgr.modeSet.disconnect(self.ModeSet)
+ if event.Checked() == True:
+ for i_tool_data in self._data:
+ i_tool_name = i_tool_data[0]
@@ -2861,10 +2131,24 @@
+ continue
+ i_tool_id = vars(self)[i_tool_name]
+ self.ToggleTool(i_tool_id, False)
-+ self.scatt_mgr.SetPlotsEditingMode(tool_name)
++ self.scatt_mgr.SetPlotsMode(tool_name)
+ else:
-+ self.scatt_mgr.SetPlotsEditingMode(None)
++ self.scatt_mgr.SetPlotsMode(None)
++ self.scatt_mgr.modeSet.connect(self.ModeSet)
+
++ def ModeSet(self, mode):
++
++ if mode in ['zoom', 'pan']:
++ self.UnsetMode()
++
++ def UnsetMode(self):
++ for i_tool_data in self._data:
++ i_tool_name = i_tool_data[0]
++ if not i_tool_name:
++ continue
++ i_tool_id = vars(self)[i_tool_name]
++ self.ToggleTool(i_tool_id, False)
++
+ def GetToolId(self, toolName):
+ return vars(self)[toolName]
Index: gui/wxpython/scatt_plot/scatt_core.py
@@ -3830,7 +3114,7 @@
===================================================================
--- gui/wxpython/scatt_plot/frame.py (revision 0)
+++ gui/wxpython/scatt_plot/frame.py (working copy)
-@@ -0,0 +1,555 @@
+@@ -0,0 +1,567 @@
+"""!
+ at package scatt_plot.dialogs
+
@@ -3903,8 +3187,8 @@
+ #self.SplitHorizontally(self.head_panel, self.plot_panel, -50)
+ self.Layout()
+
-+ def NewScatterPlot(self, scatt_id):
-+ return self.plot_panel.NewScatterPlot(scatt_id)
++ def NewScatterPlot(self, scatt_id, transpose):
++ return self.plot_panel.NewScatterPlot(scatt_id, transpose)
+
+ def ShowPlotEditingToolbar(self, show):
+ self.toolbars["editingToolbar"].Show(show)
@@ -4021,22 +3305,34 @@
+ def _getScatterPlotName(self, scatt_id):
+ return "scatter plot %d" % scatt_id
+
-+ def NewScatterPlot(self, scatt_id):
++ def NewScatterPlot(self, scatt_id, transpose):
+ #TODO needs to be resolved (should be in this class)
+
+ scatt = ScatterPlotWidget(parent = self.mainPanel,
+ scatt_mgr = self.scatt_mgr,
-+ scatt_id = scatt_id)
++ scatt_id = scatt_id,
++ transpose = transpose)
+ scatt.plotClosed.connect(self.ScatterPlotClosed)
+
+ bands = self.scatt_mgr.GetBands()
++
+ #TODO too low level
+ b1_id, b2_id = idScattToidBands(scatt_id, len(bands))
+
++ x_b = bands[b1_id].split('@')[0]
++ y_b = bands[b2_id].split('@')[0]
++
++ if transpose:
++ tmp = x_b
++ x_b = y_b
++ y_b = tmp
++
++ caption = "%s x: %s y: %s" % (_("scatter plot"), x_b, y_b)
++
+ self._mgr.AddPane(scatt,
+ aui.AuiPaneInfo().Dockable(True).Floatable(True).
+ Name(self._getScatterPlotName(scatt_id)).MinSize((-1, 300)).
-+ Caption(("%s x: %s y: %s") % (_("scatter plot"), bands[b1_id], bands[b2_id])).
++ Caption(caption).
+ Center().Position(1).MaximizeButton(True).
+ MinimizeButton(True).CaptionVisible(True).
+ CloseButton(True).Layer(0))
@@ -4281,13 +3577,13 @@
+ self.labels = {}
+ self.params = {}
+
-+ self.band_1_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("Band 1:"))
++ self.band_1_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("x axis:"))
+
+ self.band_1_ch = wx.ComboBox(parent = self, id = wx.ID_ANY,
+ choices = self.bands,
+ style = wx.CB_READONLY, size = (350, 30))
+
-+ self.band_2_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("Band 2:"))
++ self.band_2_label = wx.StaticText(parent = self, id = wx.ID_ANY, label = _("y axis:"))
+
+ self.band_2_ch = wx.ComboBox(parent = self, id = wx.ID_ANY,
+ choices = self.bands,
@@ -4402,1204 +3698,2152 @@
+ 'scatt_core',
+ 'core_c',
+ ]
-Index: gui/wxpython/Makefile
+Index: gui/wxpython/scatt_plot/plots.py
===================================================================
---- gui/wxpython/Makefile (revision 57581)
-+++ gui/wxpython/Makefile (working copy)
-@@ -13,7 +13,7 @@
- $(wildcard animation/* core/*.py dbmgr/* gcp/*.py gmodeler/* \
- gui_core/*.py iclass/* lmgr/*.py location_wizard/*.py mapwin/*.py mapdisp/*.py \
- mapswipe/* modules/*.py nviz/*.py psmap/* rlisetup/* timeline/* vdigit/* \
-- vnet/*.py web_services/*.py wxplot/*.py) \
-+ vnet/*.py web_services/*.py wxplot/*.py scatt_plot/*.py) \
- gis_set.py gis_set_error.py wxgui.py README
-
- DSTFILES := $(patsubst %,$(ETCDIR)/%,$(SRCFILES)) \
-@@ -21,8 +21,9 @@
-
- PYDSTDIRS := $(patsubst %,$(ETCDIR)/%,animation core dbmgr gcp gmodeler \
- gui_core iclass lmgr location_wizard mapwin mapdisp modules nviz psmap \
-- mapswipe vdigit wxplot web_services rlisetup vnet timeline)
-+ mapswipe vdigit wxplot web_services rlisetup vnet timeline scatt_plot)
-
+--- gui/wxpython/scatt_plot/plots.py (revision 0)
++++ gui/wxpython/scatt_plot/plots.py (working copy)
+@@ -0,0 +1,817 @@
++"""!
++ at package scatt_plot.dialogs
+
- DSTDIRS := $(patsubst %,$(ETCDIR)/%,icons scripts xml)
-
- default: $(DSTFILES)
-Index: lib/imagery/scatt.c
-===================================================================
---- lib/imagery/scatt.c (revision 0)
-+++ lib/imagery/scatt.c (working copy)
-@@ -0,0 +1,746 @@
-+/*!
-+ \file lib/imagery/scatt.c
++ at brief Ploting widgets.
+
-+ \brief Imagery library - functions for wx Scatter Plot Tool.
++Classes:
+
-+ Low level functions used by wx Scatter Plot Tool.
++(C) 2013 by the GRASS Development Team
+
-+ Copyright (C) 2013 by the GRASS Development Team
++This program is free software under the GNU General Public License
++(>=v2). Read the file COPYING that comes with GRASS for details.
+
-+ This program is free software under the GNU General Public License
-+ (>=v2). Read the file COPYING that comes with GRASS for details.
++ at author Stepan Turek <stepan.turek seznam.cz> (mentor: Martin Landa)
++"""
++import wx
++import numpy as np
+
-+ \author Stepan Turek <stepan.turek at seznam.cz> (Mentor: Martin Landa)
-+ */
++#TODO testing
++import time
++from multiprocessing import Process, Queue
+
-+#include <grass/raster.h>
-+#include <grass/imagery.h>
-+#include <grass/gis.h>
-+#include <grass/glocale.h>
++#TODO
++#from numpy.lib.stride_tricks import as_strided
++from copy import deepcopy
+
-+#include <stdio.h>
-+#include <stdlib.h>
-+#include <math.h>
-+#include <string.h>
++try:
++ import matplotlib
++ matplotlib.use('WXAgg')
++ from matplotlib.figure import Figure
++ from matplotlib.backends.backend_wxagg import \
++ FigureCanvasWxAgg as FigCanvas, \
++ NavigationToolbar2WxAgg as NavigationToolbar
++ from matplotlib.lines import Line2D
++ from matplotlib.artist import Artist
++ from matplotlib.mlab import dist_point_to_segment
++ from matplotlib.patches import Polygon, Ellipse, Rectangle
++ import matplotlib.image as mi
++ import matplotlib.colors as mcolors
++ import matplotlib.cbook as cbook
++except ImportError as e:
++ raise ImportError(_("Unable to import matplotlib (try to install it).\n%s") % e)
+
+
-+struct rast_row
-+{
-+ CELL * row;
-+ char * null_row;
-+ struct Range rast_range; /*Range of whole raster.*/
-+};
++import grass.script as grass
++from grass.pydispatch.signal import Signal
+
-+/*!
-+ \brief Create pgm header.
++#class PlotImages()?
+
-+ Scatter plot internally generates pgm files. These pgms have header in format created by this function.
++class ScatterPlotWidget(wx.Panel):
++ def __init__(self, parent, scatt_id, scatt_mgr, transpose,
++ id = wx.ID_ANY):
++
++ wx.Panel.__init__(self, parent, id)
++
++ self.parent = parent
++ self.full_extend = None
++ self.mode = None
++
++ self._createWidgets()
++ self._doLayout()
++ self.scatt_id = scatt_id
++ self.scatt_mgr = scatt_mgr
++
++ self.cidpress = None
++ self.cidrelease = None
++
++ self.transpose = transpose
++
++ self.inverse = False
++
++ self.SetSize((200, 100))
++ self.Layout()
++
++ self.base_scale = 1.2
++ self.Bind(wx.EVT_CLOSE,lambda event : self.CleanUp())
++
++ self.plotClosed = Signal("ScatterPlotWidget.plotClosed")
++
++ #self.contex_menu = ScatterPlotContextMenu(plot = self)
++
++ self.ciddscroll = None
++
++ self.canvas.mpl_connect('motion_notify_event', self.Motion)
++ self.canvas.mpl_connect('button_press_event', self.OnPress)
++ self.canvas.mpl_connect('button_release_event', self.OnRelease)
++ self.canvas.mpl_connect('draw_event', self.draw_callback)
++
++ def draw_callback(self, event):
++ self.polygon_drawer.draw_callback(event)
++ self.axes.draw_artist(self.zoom_rect)
++
++ def _createWidgets(self):
++
++ # Create the mpl Figure and FigCanvas objects.
++ # 5x4 inches, 100 dots-per-inch
++ #
++ self.dpi = 100
++ self.fig = Figure((1.0, 1.0), dpi=self.dpi)
++ self.canvas = FigCanvas(self, -1, self.fig)
++
++ self.axes = self.fig.add_axes([0,0,1,1])
++
++ pol = Polygon(list(zip([0], [0])), animated=True)
++ self.axes.add_patch(pol)
++ self.polygon_drawer = PolygonDrawer(self.axes, pol = pol, empty_pol = True)
++
++ self.zoom_wheel_coords = None
++ self.zoom_rect_coords = None
++ self.zoom_rect = Polygon(list(zip([0], [0])), facecolor = 'none')
++ self.zoom_rect.set_visible(False)
++ self.axes.add_patch(self.zoom_rect)
++
++ self.toolbar = NavigationToolbar(self.canvas)
++
++ def ZoomToExtend(self):
++ if self.full_extend:
++ self.axes.axis(self.full_extend)
++ self.canvas.draw()
++
++ def SetMode(self, mode):
++ self._deactivateMode()
++ if mode == 'zoom':
++ self.ciddscroll = self.canvas.mpl_connect('scroll_event', self.ZoomWheel)
++ self.mode = 'zoom'
++ elif mode == 'pan':
++ #self.toolbar.pan()
++ self.mode = 'pan'
++ elif mode:
++ self.polygon_drawer.SetMode(mode)
++
++ def GetCoords(self):
++
++ coords = self.polygon_drawer.GetCoords()
++ if coords is None:
++ return
++
++ if self.transpose:
++ for c in coords:
++ tmp = c[0]
++ c[0] = c[1]
++ c[1] = tmp
++
++ return coords
++
++ def SetEmpty(self):
++ return self.polygon_drawer.SetEmpty()
++
++ def _deactivateMode(self):
++
++ #TODO do own pan
++ #if self.toolbar._active == "PAN":
++ self.mode = None
++
++ if self.toolbar._active == "ZOOM":
++ self.toolbar.zoom()
++
++ if self.ciddscroll:
++ self.canvas.mpl_disconnect(self.ciddscroll)
++
++ self.zoom_rect.set_visible(False)
++ self._stopCategoryEdit()
++
++ def OnRelease(self, event):
++ if not self.mode == "zoom": return
++ self.zoom_rect.set_visible(False)
++ self.ZoomRectangle(event)
++ self.canvas.draw()
+
-+ \param region region to be pgm header generated for
-+ \param [out] header header of pgm file
-+ */
-+static int get_cat_rast_header(struct Cell_head * region, char * header){
-+ return sprintf(header, "P5\n%d\n%d\n1\n", region->cols, region->rows);
-+}
++ def OnPress(self, event):
++ 'on button press we will see if the mouse is over us and store some data'
++ if not event.inaxes:
++ return
++
++ if event.xdata and event.ydata:
++ self.zoom_wheel_coords = { 'x' : event.xdata, 'y' : event.ydata}
++ self.zoom_rect_coords = { 'x' : event.xdata, 'y' : event.ydata}
++ else:
++ self.zoom_wheel_coords = None
++ self.zoom_rect_coords = None
+
-+/*!
-+ \brief Create category raster conditions file.
++ def _stopCategoryEdit(self):
++ 'disconnect all the stored connection ids'
++
++ if self.cidpress:
++ self.canvas.mpl_disconnect(self.cidpress)
++ if self.cidrelease:
++ self.canvas.mpl_disconnect(self.cidrelease)
++ #self.canvas.mpl_disconnect(self.cidmotion)
++
++ def _doLayout(self):
++
++ self.main_sizer = wx.BoxSizer(wx.VERTICAL)
++ self.main_sizer.Add(self.canvas, 1, wx.LEFT | wx.TOP | wx.GROW)
++ #self.main_sizer.Add(self.toolbar, 0, wx.EXPAND)
++ self.SetSizer(self.main_sizer)
++ self.main_sizer.Fit(self)
+
-+ \param cat_rast_region region to be file generated for
-+ \param cat_rast path of generated category raster file
-+ */
-+int I_create_cat_rast(struct Cell_head * cat_rast_region, const char * cat_rast)
-+{
-+ FILE * f_cat_rast;
-+ char cat_rast_header[1024];//TODO magic number
-+ int i_row, i_col;
-+ int head_nchars;
++ def Plot(self, scatts, ellipses, styles):
++ """ Redraws the figure
++ """
+
-+ unsigned char * row_data;
++ callafter_list = []
+
-+ f_cat_rast = fopen(cat_rast, "wb");
-+ if(!f_cat_rast) {
-+ G_warning("Unable to create category raster condition file <%s>.", cat_rast);
-+ return -1;
-+ }
+
-+ head_nchars = get_cat_rast_header(cat_rast_region, cat_rast_header);
++ if self.full_extend:
++ cx = self.axes.get_xlim()
++ cy = self.axes.get_ylim()
++ c = cx + cy
++ else:
++ c = None
+
-+ fwrite(cat_rast_header, sizeof(char), head_nchars/sizeof(char), f_cat_rast);
-+ if (ferror(f_cat_rast)){
-+ fclose(f_cat_rast);
-+ G_warning(_("Unable to write header into category raster condition file <%s>."), cat_rast);
-+ return -1;
-+ }
++ q = Queue()
++ p = Process(target=MergeImg, args=(scatts, styles, self.transpose, q))
++ p.start()
++ merged_img, self.full_extend = q.get()
++ p.join()
+
-+ row_data = (unsigned char *) G_malloc(cat_rast_region->cols * sizeof(unsigned char));
-+ for(i_col = 0; i_col < cat_rast_region->cols; i_col++)
-+ row_data[i_col] = 0 & 255;
++ self.axes.clear()
++ img = imshow(self.axes, merged_img,
++ origin = 'lower',
++ extent = self.full_extend,
++ interpolation='nearest',
++ aspect = "auto")
+
-+ for(i_row = 0; i_row < cat_rast_region->rows; i_row++) {
-+ fwrite(row_data, sizeof(unsigned char), (cat_rast_region->cols)/sizeof(unsigned char), f_cat_rast);
-+ if (ferror(f_cat_rast))
-+ {
-+ fclose(f_cat_rast);
-+ G_warning(_("Unable to write into category raster condition file <%s>."), cat_rast);
-+ return -1;
-+ }
-+ }
++ callafter_list.append([self.axes.draw_artist, [img]])
++ callafter_list.append([grass.try_remove, [merged_img.filename]])
+
-+ fclose(f_cat_rast);
-+ return 0;
-+}
++ for cat_id, e in ellipses.iteritems():
++ if cat_id == 0 or not e:
++ continue
+
-+static int print_reg(struct Cell_head * intersec, const char * pref, int dbg_level)
-+{
-+ G_debug(dbg_level, "%s:\n n:%f\ns:%f\ne:%f\nw:%f\nns_res:%f\new_res:%f", pref, intersec->north, intersec->south,
-+ intersec->east, intersec->west, intersec->ns_res, intersec->ew_res);
-+}
++ colors = styles[cat_id]['color'].split(":")
++ if self.transpose:
++ angle = 360 - e['theta'] + 90
++ if angle >= 360:
++ angle = abs(360 - angle)
++
++ xy = [e['pos'][1], e['pos'][0]]
++ ellip = Ellipse(xy=xy,
++ width=e['width'],
++ height=e['height'],
++ angle=angle,
++ edgecolor="r",
++ facecolor = 'None')
++ else:
++ ellip = Ellipse(xy=e['pos'],
++ width=e['width'],
++ height=e['height'],
++ angle=e['theta'],
++ edgecolor="r",
++ facecolor = 'None')
+
-+/*!
-+ \brief Find intersection region of two regions.
++ self.axes.add_artist(ellip)
++ callafter_list.append([self.axes.draw_artist, [ellip]])
+
-+ \param A pointer to intersected region
-+ \param B pointer to intersected region
-+ \param [out] intersec pointer to intersection region of regions A B (relevant params of the region are: south, north, east, west)
++ callafter_list.append([self.fig.canvas.blit, []])
+
-+ \return 0 if interection exists
-+ \return -1 if regions does not intersect
-+ */
-+static int regions_intersecion(struct Cell_head * A, struct Cell_head * B, struct Cell_head * intersec)
-+{
++ if c:
++ self.axes.axis(c)
++ wx.CallAfter(lambda : self.CallAfter(callafter_list))
++
++ def CallAfter(self, funcs_list):
++ while funcs_list:
++ fcn, args = funcs_list.pop(0)
++ fcn(*args)
+
-+ if(B->north < A->south) return -1;
-+ else if(B->north > A->north) intersec->north = A->north;
-+ else intersec->north = B->north;
++ self.canvas.draw()
+
-+ if(B->south > A->north) return -1;
-+ else if(B->south < A->south) intersec->south = A->south;
-+ else intersec->south = B->south;
++ def CleanUp(self):
++ self.plotClosed.emit(scatt_id = self.scatt_id)
++ self.Destroy()
+
-+ if(B->east < A->west) return -1;
-+ else if(B->east > A->east) intersec->east = A->east;
-+ else intersec->east = B->east;
++ def ZoomWheel(self, event):
++ # get the current x and y limits
++ if not event.inaxes:
++ return
++ # tcaswell
++ # http://stackoverflow.com/questions/11551049/matplotlib-plot-zooming-with-scroll-wheel
++ cur_xlim = self.axes.get_xlim()
++ cur_ylim = self.axes.get_ylim()
++
++ xdata = event.xdata
++ ydata = event.ydata
++ if event.button == 'up':
++ scale_factor = 1/self.base_scale
++ elif event.button == 'down':
++ scale_factor = self.base_scale
++ else:
++ scale_factor = 1
+
-+ if(B->west > A->east) return -1;
-+ else if(B->west < A->west) intersec->west = A->west;
-+ else intersec->west = B->west;
++ extend = (xdata - (xdata - cur_xlim[0]) * scale_factor,
++ xdata + (cur_xlim[1] - xdata) * scale_factor,
++ ydata - (ydata - cur_ylim[0]) * scale_factor,
++ ydata + (cur_ylim[1] - ydata) * scale_factor)
+
-+ if(intersec->north == intersec->south) return -1;
++ self.axes.axis(extend)
++
++ self.canvas.draw()
+
-+ if(intersec->east == intersec->west) return -1;
+
-+ return 0;
++ def ZoomRectangle(self, event):
++ # get the current x and y limits
++ if not self.mode == "zoom": return
++ if event.inaxes is None: return
++ if event.button != 1: return
+
-+}
++ cur_xlim = self.axes.get_xlim()
++ cur_ylim = self.axes.get_ylim()
++
++ x1, y1 = event.xdata, event.ydata
++ x2 = self.zoom_rect_coords['x']
++ y2 = self.zoom_rect_coords['y']
+
-+/*!
-+ \brief Get rows and cols numbers, which defines intersection of the regions.
-+
-+ \param A pointer to intersected region
-+ \param B pointer to intersected region (A and B must have same resolution)
-+ \param [out] A_bounds rows and cols numbers of A stored in south, north, east, west, which defines intersection of A and B
-+ \param [out] B_bounds rows and cols numbers of B stored in south, north, east, west, which defines intersection of A and B
++ dx = abs(x1 - x2)
++ dy = abs(y1 - y2)
+
-+ \return 0 if interection exists
-+ \return -1 if regions do not intersect
-+ \return -2 resolution of regions is not same
-+*/
-+static int get_rows_and_cols_bounds(struct Cell_head * A, struct Cell_head * B, struct Cell_head * A_bounds, struct Cell_head * B_bounds)
-+{
-+ float ns_res, ew_res;
+
-+ struct Cell_head intersec;
++ self.axes.set_xlim(x1, x2)
++ self.axes.set_ylim(y1, y2)
++
++ def Motion(self, event):
++ self.PanMotion(event)
++ self.ZoomRectMotion(event)
+
-+ /* TODO is it right check? */
-+ if(abs(A->ns_res - B->ns_res) > GRASS_EPSILON) {
-+ G_debug(0, "'get_rows_and_cols_bounds' ns_res does not fit, A->ns_res: %f B->ns_res: %f", A->ns_res, B->ns_res);
-+ return -2;
-+ }
++ def PanMotion(self, event):
++ 'on mouse movement'
++ if not self.mode == "pan": return
++ if event.inaxes is None: return
++ if event.button != 1: return
+
-+ if(abs(A->ew_res - B->ew_res) > GRASS_EPSILON) {
-+ G_debug(0, "'get_rows_and_cols_bounds' ew_res does not fit, A->ew_res: %f B->ew_res: %f", A->ew_res, B->ew_res);
-+ return -2;
-+ }
++ cur_xlim = self.axes.get_xlim()
++ cur_ylim = self.axes.get_ylim()
+
-+ ns_res = A->ns_res;
-+ ew_res = A->ew_res;
++ x,y = event.xdata, event.ydata
++
++ mx = (x - self.zoom_wheel_coords['x']) * 0.5
++ my = (y - self.zoom_wheel_coords['y']) * 0.5
+
-+ if(regions_intersecion(A, B, &intersec) == -1)
-+ return -1;
++ extend = (cur_xlim[0] + mx, cur_xlim[1] + mx, cur_ylim[0] + my, cur_ylim[1] + my)
+
-+ A_bounds->north = ceil((A->north - intersec.north - ns_res * 0.5) / ns_res);
-+ A_bounds->south = ceil((A->north - intersec.south - ns_res * 0.5) / ns_res);
++ self.zoom_wheel_coords['x'] = x
++ self.zoom_wheel_coords['y'] = y
+
-+ A_bounds->east = ceil((intersec.east - A->west - ew_res * 0.5) / ew_res);
-+ A_bounds->west = ceil((intersec.west - A->west - ew_res * 0.5) / ew_res);
++ self.axes.axis(extend)
+
-+ B_bounds->north = ceil((B->north - intersec.north - ns_res * 0.5) / ns_res);
-+ B_bounds->south = ceil((B->north - intersec.south - ns_res * 0.5) / ns_res);
++ #self.canvas.copy_from_bbox(self.axes.bbox)
++ #self.canvas.restore_region(self.background)
++ self.canvas.draw()
++
++ def ZoomRectMotion(self, event):
++ if not self.mode == "zoom": return
++ if event.inaxes is None: return
++ if event.button != 1: return
+
-+ B_bounds->east = ceil((intersec.east - B->west - ew_res * 0.5) / ew_res);
-+ B_bounds->west = ceil((intersec.west - B->west - ew_res * 0.5) / ew_res);
++ x1, y1 = event.xdata, event.ydata
++ self.zoom_rect.set_visible(True)
++ x2 = self.zoom_rect_coords['x']
++ y2 = self.zoom_rect_coords['y']
+
-+ return 0;
-+}
++ self.zoom_rect.xy = ((x1, y1), (x1, y2), (x2, y2), (x2, y1), (x1, y1))
+
-+/*!
-+ \brief Insert raster map patch into pgm file.
-+
-+ \param patch_rast name of raster map
-+ \param cat_rast_region region of category raster file
-+ \param cat_rast path to category raster file
++ self.axes.draw_artist(self.zoom_rect)
++ self.canvas.draw()
+
-+ \return 0 on success
-+ \return -1 on failure
-+*/
-+int I_insert_patch_to_cat_rast(const char * patch_rast, struct Cell_head * cat_rast_region, const char * cat_rast)
-+{
++class ScatterPlotContextMenu:
++ def __init__(self, plot):
+
-+ FILE * f_cat_rast;
-+ struct Cell_head patch_region, patch_bounds, cat_rast_bounds;
-+ char cat_rast_header[1024];//TODO magic number
-+ int i_row, i_col, ncols, nrows, cat_rast_col, patch_col, val;
-+ int head_nchars, ret;
-+ int fd_patch_rast, init_shift, step_shift;
-+ unsigned char * patch_data;
++ self.plot = plot
++ self.canvas = plot.canvas
++ self.cidpress = self.canvas.mpl_connect(
++ 'button_press_event', self.ContexMenu)
++
++ def ContexMenu(self, event):
++ if not event.inaxes:
++ return
+
-+ char * null_chunk_row;
++ if event.button == 3:
++ menu = wx.Menu()
++ menu_items = [["zoom_to_extend", _("Zoom to scatter plot extend"), lambda event : self.plot.ZoomToExtend()]]
+
-+ const char *mapset;
++ for item in menu_items:
++ item_id = wx.ID_ANY
++ menu.Append(item_id, text = item[1])
++ menu.Bind(wx.EVT_MENU, item[2], id = item_id)
+
-+ struct Cell_head patch_lines, cat_rast_lines;
++ wx.CallAfter(self.ShowMenu, menu)
++
++ def ShowMenu(self, menu):
++ self.plot.PopupMenu(menu)
++ menu.Destroy()
++ self.plot.ReleaseMouse()
+
-+ unsigned char * row_data;
++class PolygonDrawer:
++ """
++ An polygon editor.
++ """
+
-+ f_cat_rast = fopen(cat_rast, "rb+");
-+ if(!f_cat_rast) {
-+ G_warning(_("Unable to open category raster condtions file <%s>."), cat_rast);
-+ return -1;
-+ }
++ showverts = True
++ epsilon = 5 #TODO settings max pixel distance to count as a vertex hit
++ def __init__(self, ax, pol, empty_pol):
++ if pol.figure is None:
++ raise RuntimeError('You must first add the polygon to a figure or canvas before defining the interactor')
++ self.ax = ax
++ self.canvas = pol.figure.canvas
+
-+ head_nchars = get_cat_rast_header(cat_rast_region, cat_rast_header);
-+ if ((mapset = G_find_raster(patch_rast,"")) == NULL) {
-+ fclose(f_cat_rast);
-+ G_warning(_("Unable to find patch raster <%s>."), patch_rast);
-+ return -1;
-+ }
++ self.pol = pol
++ self.empty_pol = empty_pol
+
-+ Rast_get_cellhd(patch_rast, mapset, &patch_region);
-+ Rast_set_window(&patch_region);
-+
-+ if ((fd_patch_rast = Rast_open_old(patch_rast, mapset)) < 0) {
-+ fclose(f_cat_rast);
-+ return -1;
-+ }
++ x, y = zip(*self.pol.xy)
++ self.line = Line2D(x, y, marker='o', markerfacecolor='r', animated=True)
++ self.ax.add_line(self.line)
++ #self._update_line(pol)
+
-+ ret = get_rows_and_cols_bounds(cat_rast_region, &patch_region, &cat_rast_bounds, &patch_bounds);
-+ if(ret == -2) {
-+ G_warning(_("Resolutions of patch <%s> and patched file <%s> are not same."), patch_rast, cat_rast);
++ cid = self.pol.add_callback(self.poly_changed)
++ self.moving_ver_idx = None # the active vert
+
-+ Rast_close(fd_patch_rast);
-+ fclose(f_cat_rast);
++ self.mode = None
+
-+ return -1;
-+ }
-+ else if (ret == -1){
++ if self.empty_pol:
++ self.Show(False)
+
-+ Rast_close(fd_patch_rast);
-+ fclose(f_cat_rast);
++ #self.canvas.mpl_connect('draw_event', self.draw_callback)
++ self.canvas.mpl_connect('button_press_event', self.OnButtonPressed)
++ self.canvas.mpl_connect('button_release_event', self.button_release_callback)
++ self.canvas.mpl_connect('motion_notify_event', self.motion_notify_callback)
++
++ self.it = 0
++ def SetMode(self, mode):
++ self.mode = mode
+
-+ return 0;
-+ }
++ def GetCoords(self):
++ if self.empty_pol:
++ return None
+
-+ ncols = cat_rast_bounds.east - cat_rast_bounds.west;
-+ nrows = cat_rast_bounds.south - cat_rast_bounds.north;
++ coords = deepcopy(self.pol.xy)
++ return coords
+
-+ patch_data = (unsigned char *) G_malloc(ncols * sizeof(unsigned char));
++ def SetEmpty(self):
++ self._setEmptyPol(True)
+
-+ init_shift = head_nchars + cat_rast_region->cols * cat_rast_bounds.north + cat_rast_bounds.west;
++ def _setEmptyPol(self, empty_pol):
++ self.empty_pol = empty_pol
++ self.Show(not empty_pol)
+
-+ if(fseek(f_cat_rast, init_shift, SEEK_SET) != 0) {
-+ G_warning(_("Corrupted category raster conditions file <%s> (fseek failed)"), cat_rast);
++ def Show(self, show):
+
-+ Rast_close(fd_patch_rast);
-+ G_free(null_chunk_row);
-+ fclose(f_cat_rast);
++ self.show = show
+
-+ return -1;
-+ }
++ self.line.set_visible(self.show)
++ self.pol.set_visible(self.show)
+
-+ step_shift = cat_rast_region->cols - ncols;
++ self.Redraw()
+
-+ null_chunk_row = Rast_allocate_null_buf();
++ def Redraw(self):
++ if self.show:
++ self.ax.draw_artist(self.pol)
++ self.ax.draw_artist(self.line)
++ self.canvas.blit(self.ax.bbox)
++ self.canvas.draw()
++
++ def draw_callback(self, event):
++ self.background = self.canvas.copy_from_bbox(self.ax.bbox)
++ self.ax.draw_artist(self.pol)
++ self.ax.draw_artist(self.line)
+
-+ for(i_row = 0; i_row < nrows; i_row++) {
-+ Rast_get_null_value_row (fd_patch_rast, null_chunk_row, i_row + patch_bounds.north);
++ def poly_changed(self, pol):
++ 'this method is called whenever the polygon object is called'
++ # only copy the artist props to the line (except visibility)
++ vis = self.line.get_visible()
++ Artist.update_from(self.line, pol)
++ self.line.set_visible(vis) # don't use the pol visibility state
+
-+ for(i_col = 0; i_col < ncols; i_col++) {
-+ patch_col = patch_bounds.west + i_col;
++ def get_ind_under_point(self, event):
++ 'get the index of the vertex under point if within epsilon tolerance'
+
-+ if(null_chunk_row[patch_col] != 1)
-+ patch_data[i_col] = 1 & 255;
-+ else {
-+ patch_data[i_col] = 0 & 255;
-+ }
-+ }
++ # display coords
++ xy = np.asarray(self.pol.xy)
++ xyt = self.pol.get_transform().transform(xy)
++ xt, yt = xyt[:, 0], xyt[:, 1]
++ d = np.sqrt((xt-event.x)**2 + (yt-event.y)**2)
++ indseq = np.nonzero(np.equal(d, np.amin(d)))[0]
++ ind = indseq[0]
+
-+ fwrite(patch_data, sizeof(unsigned char), (ncols)/sizeof(unsigned char), f_cat_rast);
-+ if (ferror(f_cat_rast))
-+ {
-+ G_warning(_("Unable to write into category raster conditions file <%s>"), cat_rast);
++ if d[ind]>=self.epsilon:
++ ind = None
++
++ return ind
++
++ def OnButtonPressed(self, event):
++ if not event.inaxes:
++ return
++
++ if event.button in [2, 3]:
++ return
++
++ if self.mode == "delete_vertex":
++ self._deleteVertex(event)
++ elif self.mode == "add_boundary_vertex":
++ self._addVertexOnBoundary(event)
++ elif self.mode == "add_vertex":
++ self._addVertex(event)
++ self.moving_ver_idx = self.get_ind_under_point(event)
++
++ def button_release_callback(self, event):
++ 'whenever a mouse button is released'
++ if not self.showverts: return
++ if event.button != 1: return
++ self.moving_ver_idx = None
++
++ def ShowVertices(self, show):
++ self.showverts = show
++ self.line.set_visible(self.showverts)
++ if not self.showverts: self.moving_ver_idx = None
++
++ def _deleteVertex(self, event):
++ ind = self.get_ind_under_point(event)
++
++ if ind is None or self.empty_pol:
++ return
++
++ if len(self.pol.xy) <= 2:
++ self.empty_pol = True
++ self.Show(False)
++ return
++
++ coords = []
++ for i,tup in enumerate(self.pol.xy):
++ if i == ind:
++ continue
++ elif i == 0 and ind == len(self.pol.xy) - 1:
++ continue
++ elif i == len(self.pol.xy) - 1 and ind == 0:
++ continue
++
++ coords.append(tup)
++
++ self.pol.xy = coords
++ self.line.set_data(zip(*self.pol.xy))
++
++ self.Redraw()
++
++ def _addVertexOnBoundary(self, event):
++ if self.empty_pol:
++ return
++
++ xys = self.pol.get_transform().transform(self.pol.xy)
++ p = event.x, event.y # display coords
++ for i in range(len(xys)-1):
++ s0 = xys[i]
++ s1 = xys[i+1]
++ d = dist_point_to_segment(p, s0, s1)
++
++ if d<=self.epsilon:
++ self.pol.xy = np.array(
++ list(self.pol.xy[:i + 1]) +
++ [(event.xdata, event.ydata)] +
++ list(self.pol.xy[i + 1:]))
++ self.line.set_data(zip(*self.pol.xy))
++ break
++
++ self.Redraw()
++
++ def _addVertex(self, event):
++
++ if self.empty_pol:
++ pt = (event.xdata, event.ydata)
++ self.pol.xy = np.array([pt, pt])
++ self.Show(True)
++ self.empty_pol = False
++ else:
++ self.pol.xy = np.array(
++ [(event.xdata, event.ydata)] +
++ list(self.pol.xy[1:]) +
++ [(event.xdata, event.ydata)])
++
++ self.line.set_data(zip(*self.pol.xy))
++
++ self.Redraw()
++
++ def motion_notify_callback(self, event):
++ 'on mouse movement'
++ if not self.mode == "move_vertex": return
++ if not self.showverts: return
++ if self.empty_pol: return
++ if self.moving_ver_idx is None: return
++ if event.inaxes is None: return
++ if event.button != 1: return
++
++ self.it += 1
++
++ x,y = event.xdata, event.ydata
++
++ self.pol.xy[self.moving_ver_idx] = x,y
++ if self.moving_ver_idx == 0:
++ self.pol.xy[len(self.pol.xy) - 1] = x,y
++ elif self.moving_ver_idx == len(self.pol.xy) - 1:
++ self.pol.xy[0] = x,y
++
++ self.line.set_data(zip(*self.pol.xy))
++
++ self.canvas.restore_region(self.background)
++
++ self.Redraw()
++
++class ModestImage(mi.AxesImage):
++ """
++ Computationally modest image class.
++
++ ModestImage is an extension of the Matplotlib AxesImage class
++ better suited for the interactive display of larger images. Before
++ drawing, ModestImage resamples the data array based on the screen
++ resolution and view window. This has very little affect on the
++ appearance of the image, but can substantially cut down on
++ computation since calculations of unresolved or clipped pixels
++ are skipped.
++
++ The interface of ModestImage is the same as AxesImage. However, it
++ does not currently support setting the 'extent' property. There
++ may also be weird coordinate warping operations for images that
++ I'm not aware of. Don't expect those to work either.
++
++ Author: Chris Beaumont <beaumont at hawaii.edu>
++ """
++ def __init__(self, *args, **kwargs):
++ #why???
++ #if 'extent' in kwargs and kwargs['extent'] is not None:
++ # raise NotImplementedError("ModestImage does not support extents")
++
++ self._full_res = None
++ self._sx, self._sy = None, None
++ self._bounds = (None, None, None, None)
++ super(ModestImage, self).__init__(*args, **kwargs)
++
++ def set_data(self, A):
++ """
++ Set the image array
++
++ ACCEPTS: numpy/PIL Image A
++ """
++ self._full_res = A
++ self._A = A
++
++ if self._A.dtype != np.uint8 and not np.can_cast(self._A.dtype,
++ np.float):
++ raise TypeError("Image data can not convert to float")
++
++ if (self._A.ndim not in (2, 3) or
++ (self._A.ndim == 3 and self._A.shape[-1] not in (3, 4))):
++ raise TypeError("Invalid dimensions for image data")
++
++ self._imcache =None
++ self._rgbacache = None
++ self._oldxslice = None
++ self._oldyslice = None
++ self._sx, self._sy = None, None
++
++ def get_array(self):
++ """Override to return the full-resolution array"""
++ return self._full_res
++
++ def _scale_to_res(self):
++ """ Change self._A and _extent to render an image whose
++ resolution is matched to the eventual rendering."""
++
++ ax = self.axes
++ ext = ax.transAxes.transform([1, 1]) - ax.transAxes.transform([0, 0])
++ xlim, ylim = ax.get_xlim(), ax.get_ylim()
++ dx, dy = xlim[1] - xlim[0], ylim[1] - ylim[0]
++
++ y0 = max(0, ylim[0] - 5)
++ y1 = min(self._full_res.shape[0], ylim[1] + 5)
++ x0 = max(0, xlim[0] - 5)
++ x1 = min(self._full_res.shape[1], xlim[1] + 5)
++ y0, y1, x0, x1 = map(int, [y0, y1, x0, x1])
++
++ sy = int(max(1, min((y1 - y0) / 5., np.ceil(dy / ext[1]))))
++ sx = int(max(1, min((x1 - x0) / 5., np.ceil(dx / ext[0]))))
++
++ # have we already calculated what we need?
++ if sx == self._sx and sy == self._sy and \
++ x0 == self._bounds[0] and x1 == self._bounds[1] and \
++ y0 == self._bounds[2] and y1 == self._bounds[3]:
++ return
++
++ self._A = self._full_res[y0:y1:sy, x0:x1:sx]
++ self._A = cbook.safe_masked_invalid(self._A)
++ x1 = x0 + self._A.shape[1] * sx
++ y1 = y0 + self._A.shape[0] * sy
++
++ self.set_extent([x0 - .5, x1 - .5, y0 - .5, y1 - .5])
++ self._sx = sx
++ self._sy = sy
++ self._bounds = (x0, x1, y0, y1)
++ self.changed()
++
++ def draw(self, renderer, *args, **kwargs):
++ self._scale_to_res()
++ super(ModestImage, self).draw(renderer, *args, **kwargs)
++
++
++def imshow(axes, X, cmap=None, norm=None, aspect=None,
++ interpolation=None, alpha=None, vmin=None, vmax=None,
++ origin=None, extent=None, shape=None, filternorm=1,
++ filterrad=4.0, imlim=None, resample=None, url=None, **kwargs):
++ """Similar to matplotlib's imshow command, but produces a ModestImage
++
++ Unlike matplotlib version, must explicitly specify axes
++ Author: Chris Beaumont <beaumont at hawaii.edu>
++ """
++
++ if not axes._hold:
++ axes.cla()
++ if norm is not None:
++ assert(isinstance(norm, mcolors.Normalize))
++ if aspect is None:
++ aspect = rcParams['image.aspect']
++ axes.set_aspect(aspect)
++ im = ModestImage(axes, cmap, norm, interpolation, origin, extent,
++ filternorm=filternorm,
++ filterrad=filterrad, resample=resample, **kwargs)
++
++ im.set_data(X)
++ im.set_alpha(alpha)
++ axes._set_artist_props(im)
++
++ if im.get_clip_path() is None:
++ # image does not already have clipping set, clip to axes patch
++ im.set_clip_path(axes.patch)
++
++ #if norm is None and shape is None:
++ # im.set_clim(vmin, vmax)
++ if vmin is not None or vmax is not None:
++ im.set_clim(vmin, vmax)
++ else:
++ im.autoscale_None()
++ im.set_url(url)
++
++ # update ax.dataLim, and, if autoscaling, set viewLim
++ # to tightly fit the image, regardless of dataLim.
++ im.set_extent(im.get_extent())
++
++ axes.images.append(im)
++ im._remove_method = lambda h: axes.images.remove(h)
++
++ return im
++
++def MergeImg(scatts, styles, transpose, output_queue):
++
++ init = True
++ merge_tmp = grass.tempfile()
++ for cat_id, scatt in scatts.iteritems():
++ #print "color map %d" % cat_id
++ #TODO make more general
++ if cat_id != 0 and (styles[cat_id]['opacity'] == 0.0 or \
++ not styles[cat_id]['show']):
++ continue
++ if init:
++ if transpose:
++ b1_i = scatt['bands_info']['b1']
++ b2_i = scatt['bands_info']['b2']
++ else:
++ b2_i = scatt['bands_info']['b1']
++ b1_i = scatt['bands_info']['b2']
++
++ full_extend = (b1_i['min'] - 0.5, b1_i['max'] + 0.5, b2_i['min'] - 0.5, b2_i['max'] + 0.5)
+
-+ Rast_close(fd_patch_rast);
-+ G_free(null_chunk_row);
-+ fclose(f_cat_rast);
++ if cat_id == 0:
++ cmap = matplotlib.cm.jet
++ cmap.set_bad('w',1.)
++ cmap._init()
++ cmap._lut[len(cmap._lut) - 1, -1] = 0
++ else:
++ colors = styles[cat_id]['color'].split(":")
+
-+ return -1;
-+ }
-+ if(fseek(f_cat_rast, step_shift, SEEK_CUR) != 0) {
-+ G_warning(_("Corrupted category raster conditions file <%s> (fseek failed)"), cat_rast);
++ cmap = matplotlib.cm.jet
++ cmap.set_bad('w',1.)
++ cmap._init()
++ cmap._lut[len(cmap._lut) - 1, -1] = 0
++ cmap._lut[:, 0] = int(colors[0])/255.0
++ cmap._lut[:, 1] = int(colors[1])/255.0
++ cmap._lut[:, 2] = int(colors[2])/255.0
++
++ #if init:
++ if transpose:
++ masked_cat = np.ma.masked_less_equal(scatt['np_vals'].T, 0)
++ else:
++ masked_cat = np.ma.masked_less_equal(scatt['np_vals'], 0)
++
++ vmax = np.amax(masked_cat)
++ masked_cat = masked_cat / float(vmax)
++
++ colored_cat = np.uint8(cmap(masked_cat) * 255)
++ del masked_cat
++ del cmap
+
-+ Rast_close(fd_patch_rast);
-+ G_free(null_chunk_row);
-+ fclose(f_cat_rast);
++ #colored_cat[...,3] = np.choose(masked_cat.mask, (255, 0))
+
-+ return -1;
-+ }
-+ }
++ if init:
++ merged_img = np.memmap(merge_tmp, dtype='uint8', mode='w+', shape=colored_cat.shape)
++ merged_img[:] = colored_cat[:]
++ init = False
++ else:
++ #c_img_a = np.memmap(grass.tempfile(), dtype="uint16", mode='w+', shape = shape)
++ c_img_a = colored_cat.astype('uint16')[:,:,3] * styles[cat_id]['opacity']
+
-+ Rast_close(fd_patch_rast);
-+ G_free(null_chunk_row);
-+ fclose(f_cat_rast);
-+ return 0;
-+}
++ #TODO apply strides and there will be no need for loop
++ #b = as_strided(a, strides=(0, a.strides[3], a.strides[3], a.strides[3]), shape=(3, a.shape[0], a.shape[1]))
+
-+/*!
-+ \brief Updates scatter plots data in category by pixels which meets category conditions.
-+
-+ \param bands_rows data represents data describig one row from raster band
-+ \param belongs_pix array which defines which pixels belongs to category (1 value) and which not (0 value)
-+ \param [out] scatts pointer to scScatts struct of type SC_SCATT_DATA, which are modified according to values in belongs_pix
-+*/
-+static inline void update_cat_scatt_plts(struct rast_row * bands_rows, unsigned short * belongs_pix, struct scScatts * scatts)
-+{
-+ int band_axis_1, band_axis_2, i_scatt, array_idx, cat_idx, i_chunk_rows_pix, max_arr_idx;
++ for i in range(3):
++ merged_img[:,:,i] = (merged_img[:,:,i] * (255 - c_img_a) + colored_cat[:,:,i] * c_img_a) / 255;
++ merged_img[:,:,3] = (merged_img[:,:,3] * (255 - c_img_a) + 255 * c_img_a) / 255;
++
++ del c_img_a
++
++ del colored_cat
++ output_queue.put((merged_img, full_extend))
++ #return merged_img, full_extend
+\ No newline at end of file
+Index: gui/wxpython/scatt_plot/controllers.py
+===================================================================
+--- gui/wxpython/scatt_plot/controllers.py (revision 0)
++++ gui/wxpython/scatt_plot/controllers.py (working copy)
+@@ -0,0 +1,686 @@
++"""!
++ at package scatt_plot.controllers
+
-+ CELL * b_1_row;
-+ CELL * b_2_row;
-+ char * b_1_null_row,* b_2_null_row;
-+ struct rast_row b_1_rast_row, b_2_rast_row;
++ at brief Controller layer for scatter plot tool.
+
-+ struct Range b_1_range, b_2_range;
-+ int b_1_range_size;
++Classes:
+
-+ int row_size = Rast_window_cols();
++(C) 2013 by the GRASS Development Team
+
-+ int * scatts_bands = scatts->scatts_bands;
++This program is free software under the GNU General Public License
++(>=v2). Read the file COPYING that comes with GRASS for details.
+
-+ for(i_scatt = 0; i_scatt < scatts->n_a_scatts; i_scatt++)
-+ {
-+ b_1_rast_row = bands_rows[scatts_bands[i_scatt * 2]];
-+ b_2_rast_row = bands_rows[scatts_bands[i_scatt * 2 + 1]];
++ at author Stepan Turek <stepan.turek seznam.cz> (mentor: Martin Landa)
++"""
++import os
++import sys
+
-+ b_1_row = b_1_rast_row.row;
-+ b_2_row = b_2_rast_row.row;
++#TODO just for testing
++import time
+
-+ b_1_null_row = b_1_rast_row.null_row;
-+ b_2_null_row = b_2_rast_row.null_row;
++#TODO
++import wx
+
-+ b_1_range = b_1_rast_row.rast_range;
-+ b_2_range = b_2_rast_row.rast_range;
++from core.gcmd import GException, GError, GMessage, RunCommand
++
++from scatt_plot.scatt_core import Core, idBandsToidScatt
++
++from scatt_plot.dialogs import AddScattPlotDialog
++from scatt_plot.gthreading import gThread
++from core.gconsole import EVT_CMD_DONE
++from grass.pydispatch.signal import Signal
++
++class ScattsManager(wx.EvtHandler):
++ def __init__(self, guiparent, giface, iclass_mapwin):
++ #TODO remove iclass parameter
++
++ wx.EvtHandler.__init__(self)
++ self.giface = giface
++ self.mapDisp = giface.GetMapDisplay()
++
++ if iclass_mapwin:
++ self.mapWin = iclass_mapwin
++ else:
++ self.mapWin = giface.GetMapWindow()
++
++ self.guiparent = guiparent
++
++ self.show_add_scatt_plot = False
++
++ self.core = Core()
++ self.scatts_dt, self.scatt_conds_dt = self.core.GetScattsData()
++
++ self.cats_mgr = CategoriesManager(self, self.core)
++
++ self.thread = gThread(self);
+
-+ b_1_range_size = b_1_range.max - b_1_range.min + 1;
-+ max_arr_idx = (b_1_range.max - b_1_range.min + 1) * (b_2_range.max - b_2_range.min + 1);
++ self.plots = {}
++ self.added_cats_rasts = {}
+
-+ for(i_chunk_rows_pix = 0; i_chunk_rows_pix < row_size; i_chunk_rows_pix++)
-+ {
-+ /* pixel does not belongs to scatter plot or has null value in one of the bands */
-+ if(!belongs_pix[i_chunk_rows_pix] ||
-+ b_1_null_row[i_chunk_rows_pix] == 1 ||
-+ b_2_null_row[i_chunk_rows_pix] == 1)
-+ continue;
++ self.cats_to_update = []
+
-+ /* index in scatter plot array */
-+ array_idx = b_1_row[i_chunk_rows_pix] - b_1_range.min + (b_2_row[i_chunk_rows_pix] - b_2_range.min) * b_1_range_size;
++ self.plot_mode = None
++ self.pol_sel_mode = [False, None]
+
-+ if(array_idx < 0 || array_idx >= max_arr_idx) {
-+ G_warning ("Data inconsistent. Value computed for scatter plot is out of initialized range.");
-+ continue;
-+ }
++ self.data_set = False
+
-+ /* increment scatter plot value */
-+ ++scatts->scatts_arr[i_scatt]->scatt_vals_arr[array_idx];
-+ }
-+ }
-+}
++ if iclass_mapwin:
++ self.mapWin_conn = MapWinConnection(self, self.mapWin, self.core.CatRastUpdater())
++ self.iclass_conn = IClassConnection(self, iclass_mapwin.parent, self.cats_mgr)
++ else:
++ self.mapWin_conn = None
++ self.iclass_conn = None
+
-+/*!
-+ \brief Computes scatter plots data from chunk_rows.
++ self.tasks_pids = {
++ 'add_scatt' : [],
++ 'set_data' : -1,
++ 'set_data_add' : -1,
++ 'set_edit_cat_data' : -1,
++ 'mapwin_conn' : [],
++ 'render_plots' : -1,
++ 'render' : []
++ }
+
-+ \param scatts pointer to scScatts struct of type SC_SCATT_DATA, where are computed scatter plots stored
-+ \param scatt_conds pointer to scScatts struct of type SC_SCATT_CONDITIONS, where are selected areas (condtitions)stored
++ self.Bind(EVT_CMD_DONE, self.OnThreadDone)
+
-+ \param chunk_rows data arrays of chunk_rows from analyzed raster bands (all data in chunk_rows, null_chunk_rows and belongs_pix arrays represents same region)
-+ \param null_chunk_rows null data arrays of chunk_rows from analyzed raster bands
-+ \param row_size size of data in chunk_rows, null_chunk_rows and belongs_pix arrays
-+ \param f_cats_rasts_conds file which stores selected areas (conditions) from mapwindow see I_create_cat_rast() and I_pa
-+ \param fd_cats_rasts array of openedraster maps which represents all selected pixels for category
-+ \param region analysis region
++ self.modeSet = Signal("ScattsManager.mondeSet")
+
-+ \return 0 on success
-+ \return -1 on failure
-+*/
-+static inline int compute_scatts_from_chunk_row(struct Cell_head *region, struct scCats * scatt_conds,
-+ FILE ** f_cats_rasts_conds, struct rast_row * bands_rows,
-+ struct scCats * scatts, int * fd_cats_rasts)
-+{
++ def CleanUp(self):
++ self.core.CleanUp()
+
-+ int i_rows_pix, i_cat, i_scatt, n_a_scatts, n_pixs;
-+ int cat_id, scatt_plts_cat_idx, array_idx, max_arr_idx;
-+ char * b_1_null_row,* b_2_null_row;
-+ struct rast_row b_1_rast_row, b_2_rast_row;
-+ CELL * cat_rast_row;
++ for scatt_id, scatt in self.plots.items():
++ scatt.CleanUp()
+
-+ struct scScatts * scatts_conds;
-+ struct scScatts * scatts_scatt_plts;
-+ struct scdScattData * conds;
++ def OnThreadDone(self, event):
++
++ if event.exception:
++ GError(str(event.exception))
++ return
+
-+ struct Range b_1_range, b_2_range;
-+ int b_1_range_size;
++ if event.pid in self.tasks_pids['mapwin_conn']:
++ self.tasks_pids['mapwin_conn'].remove(event.pid)
++ updated_cats = event.ret
+
-+ int * scatts_bands;
-+ struct scdScattData ** scatts_arr;
++ for cat in updated_cats:
++ if cat not in self.cats_to_update:
++ self.cats_to_update.append(cat)
+
-+ CELL * b_1_row;
-+ CELL * b_2_row;
-+ unsigned char * i_scatt_conds;
++ if not self.tasks_pids['mapwin_conn']:
++ self.tasks_pids['render_plots'] = self.thread.GetId()
++ self.thread.Run(callable = self.core.ComputeCatsScatts,
++ cats_ids = self.cats_to_update[:])
++ del self.cats_to_update[:]
++
++ return
+
-+ int row_size = Rast_window_cols();
++ if self.tasks_pids['render_plots'] == event.pid:
++ self.RenderScattPlts()
++ return
+
-+ unsigned short * belongs_pix = (unsigned short *) G_malloc(row_size * sizeof(unsigned short));
-+ unsigned char * rast_pixs = (unsigned char *) G_malloc(row_size * sizeof(unsigned char));
-+ cat_rast_row = Rast_allocate_c_buf();
++ if event.pid in self.tasks_pids['render']:
++ self.tasks_pids['render'].remove(event.pid)
++ return
++
++ if event.pid in self.tasks_pids['add_scatt']:
++ self.tasks_pids['add_scatt'].remove(event.pid)
++ self.AddScattPlotDone(event)
++ return
+
-+
-+ for(i_cat = 0; i_cat < scatt_conds->n_a_cats; i_cat++)
-+ {
-+ scatts_conds = scatt_conds->cats_arr[i_cat];
++ if self.tasks_pids['set_data'] == event.pid:
++ self.SetDataDone(event)
++ return
++
++ if self.tasks_pids['set_data_add'] == event.pid:
++ self.SetDataDone(event)
++ self.AddScattPlot()
++ return
++
++ if self.tasks_pids['set_edit_cat_data'] == event.pid:
++ self.SetEditCatDataDone(event)
++ return
+
-+ cat_id = scatt_conds->cats_ids[i_cat];
++ def SetData(self, bands):
+
-+ scatt_plts_cat_idx = scatts->cats_idxs[cat_id];
-+ if(scatt_plts_cat_idx < 0)
-+ continue;
-+ scatts_scatt_plts = scatts->cats_arr[scatt_plts_cat_idx];
++ self.CleanUp()
+
-+ G_zero(belongs_pix, row_size * sizeof(unsigned short));
++ self.data_set = False
+
-+ /* if category has no conditions defined, scatter plots without any constraint are computed (default scatter plots) */
-+ if(!scatts_conds->n_a_scatts && !f_cats_rasts_conds[i_cat]) {
-+ for(i_scatt = 0; i_scatt < scatts_scatt_plts->n_a_scatts; i_scatt++)
-+ {
-+ /* all pixels belongs */
-+ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
-+ belongs_pix[i_rows_pix] = 1;
-+ }
-+ }
-+ /* compute belonging pixels for defined conditions */
-+ else
-+ {
-+ scatts_bands = scatts_conds->scatts_bands;
++ if self.show_add_scatt_plot:
++ self.tasks_pids['set_data_add'] = self.thread.GetId()
++ else:
++ self.tasks_pids['set_data'] = self.thread.GetId()
+
-+ /* check conditions from category raster condtitions file */
-+ if(f_cats_rasts_conds[i_cat]) {
-+ n_pixs = fread(rast_pixs, sizeof(unsigned char), (row_size)/sizeof(unsigned char), f_cats_rasts_conds[i_cat]);
++ self.thread.Run(callable = self.core.SetData, bands = bands)
+
-+ if (ferror(f_cats_rasts_conds[i_cat]))
-+ {
-+ G_free(rast_pixs);
-+ G_free(belongs_pix);
-+ G_warning(_("Unable to read from category raster condtition file."));
-+ return -1;
-+ }
-+ if (n_pixs != n_pixs) {
-+ G_free(rast_pixs);
-+ G_free(belongs_pix);
-+ G_warning(_("Invalid size of category raster conditions file."));
-+ return -1;
++ def SetDataDone(self, event):
+
-+ }
++ self.data_set = True
++ self.cats_mgr.InitCoreCats()
+
-+ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
-+ {
-+ if(rast_pixs[i_rows_pix] != 0 & 255)
-+ belongs_pix[i_rows_pix] = 1;
-+ }
-+ }
++ def OnOutput(self, event):
++ """!Print thread output according to debug level.
++ """
++ print event.text
+
-+ /* check condtions defined in scatter plots*/
-+ for(i_scatt = 0; i_scatt < scatts_conds->n_a_scatts; i_scatt++)
-+ {
-+ b_1_rast_row = bands_rows[scatts_bands[i_scatt * 2]];
-+ b_2_rast_row = bands_rows[scatts_bands[i_scatt * 2 + 1]];
++ def GetBands(self):
++ return self.core.GetBands()
+
-+ b_1_row = b_1_rast_row.row;
-+ b_2_row = b_2_rast_row.row;
++ def AddScattPlot(self):
++ if not self.data_set and self.iclass_conn:
++ self.show_add_scatt_plot = True
++ self.iclass_conn.SetData()
++ self.show_add_scatt_plot = False
++ return
++ if not self.data_set:
++ GError(_('No data set.'))
++ return
+
-+ b_1_null_row = b_1_rast_row.null_row;
-+ b_2_null_row = b_2_rast_row.null_row;
++ bands = self.core.GetBands()
++ dlg = AddScattPlotDialog(parent = self.guiparent, bands = bands)
++ if dlg.ShowModal() == wx.ID_OK:
+
-+ b_1_range = b_1_rast_row.rast_range;
-+ b_2_range = b_2_rast_row.rast_range;
++ b_1, b_2 = dlg.GetBands()
++
++ #TODO axes selection
++ transpose = False
++ if b_1 > b_2:
++ transpose = True
++ tmp_band = b_2
++ b_2 = b_1
++ b_1 = tmp_band
+
-+ b_1_range_size = b_1_range.max - b_1_range.min + 1;
-+ max_arr_idx = (b_1_range.max - b_1_range.min + 1) * (b_2_range.max - b_2_range.min + 1);
++ self.scatt_id = idBandsToidScatt(b_1, b_2, len(bands))
++ self._addScattPlot(self.scatt_id, transpose)
++
++ dlg.Destroy()
+
-+ i_scatt_conds = scatts_conds->scatts_arr[i_scatt]->b_conds_arr;
++ return
+
-+ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
-+ {
-+ /* pixels already belongs to category from category raster conditions file or contains null value in one of the bands */
-+ if(belongs_pix[i_rows_pix] ||
-+ b_1_null_row[i_rows_pix] == 1 ||
-+ b_2_null_row[i_rows_pix] == 1)
-+ continue;
++ #for testing
++ jj = len(self.core.GetBands()) - 1
+
-+ array_idx = b_1_row[i_rows_pix] - b_1_range.min + (b_2_row[i_rows_pix] - b_2_range.min) * b_1_range_size;
-+ if(array_idx < 0 || array_idx >= max_arr_idx) {
-+ G_warning ("Data inconsistent. Value computed for scatter plot is out of initialized range.");
-+ continue;
-+ }
-+ /* pixels meets condtion defined in scatter plot */
-+ if(i_scatt_conds[array_idx])
-+ belongs_pix[i_rows_pix] = 1;
-+ }
-+ }
-+ }
++ max_k = 20
++ k = 0
++ for i in range(jj):
++ for j in range(i, jj):
++ print i * (len(self.core.GetBands()) - 1) + j
++ self._addScattPlot(i * (len(self.core.GetBands()) - 1) + j)
++ k += 1
++ if k == max_k:
++ break
++ if k == max_k:
++ break
++ jj -= 1
+
-+ /* update category raster with belonging pixels */
-+ if(fd_cats_rasts[i_cat] >= 0) {
-+ Rast_set_null_value(cat_rast_row, Rast_window_cols(), CELL_TYPE);
+
-+ for(i_rows_pix = 0; i_rows_pix < row_size; i_rows_pix++)
-+ if(belongs_pix[i_rows_pix])
-+ cat_rast_row[i_rows_pix] = belongs_pix[i_rows_pix];
++ def _addScattPlot(self, scatt_id, transpose):
++ if self.plots.has_key(scatt_id):
++ GMessage(_("Scatter plot has been already added."))
++ return
+
-+ Rast_put_c_row (fd_cats_rasts[i_cat], cat_rast_row);
-+ }
++ self.tasks_pids['add_scatt'].append(self.thread.GetId())
++
++ self.thread.Run(callable = self.core.AddScattPlot, scatt_id = scatt_id, userdata = {'transpose' : transpose})
+
-+ /* update scatter plots with belonging pixels */
-+ update_cat_scatt_plts(bands_rows, belongs_pix, scatts_scatt_plts);
-+ }
+
-+ G_free(cat_rast_row);
-+ G_free(rast_pixs);
-+ G_free(belongs_pix);
++ def RenderScattPlts(self, scatt_ids = None):
++ if len(self.tasks_pids['render']) > 1:
++ print "skip"
++ return
+
-+ return 0;
-+}
++ self.tasks_pids['render'].append(self.thread.GetId())
++ self.thread.Run(callable = self._renderscattplts, scatt_ids = scatt_ids)
+
-+/*!
-+ \brief Get list if bands needed to be opened for analysis from scCats struct.
-+*/
-+static void get_needed_bands(struct scCats * cats, int * b_needed_bands)
-+{
-+ // results in b_needed_bands - array of bools - if item has value 1, band (defined by item index) is needed to be opened
-+ int i_cat, i_scatt, cat_id;
++ def _renderscattplts(self, scatt_ids):
++ cats_attrs = self.cats_mgr.GetCategoriesAttrs()
++ for i_scatt_id, scatt in self.plots.items():
++ if scatt_ids is not None and i_scatt_id not in scatt_ids:
++ continue
+
-+ for(i_cat = 0; i_cat < cats->n_a_cats; i_cat++)
-+ {
-+ for(i_scatt = 0; i_scatt < cats->cats_arr[i_cat]->n_a_scatts; i_scatt++)
-+ {
-+ G_debug(3, "Active scatt %d in catt %d", i_scatt, i_cat);
++ scatt_dt = self.scatts_dt.GetScatt(i_scatt_id)
++ ellipses_dt = self.scatts_dt.GetEllipses(i_scatt_id)
+
-+ b_needed_bands[cats->cats_arr[i_cat]->scatts_bands[i_scatt * 2]] = 1;
-+ b_needed_bands[cats->cats_arr[i_cat]->scatts_bands[i_scatt * 2 + 1]] = 1;
-+ }
-+ }
-+ return;
-+}
++ if self.pol_sel_mode[0]:
++ self._getSelectedAreas(i_scatt_id, scatt_dt, cats_attrs)
+
-+/*!
-+ \brief Helper function for clean up.
-+*/
-+static void free_compute_scatts_data(int * fd_bands, struct rast_row * bands_rows, int n_a_bands, int * bands_ids,
-+ int * fd_cats_rasts, FILE ** f_cats_rasts_conds, int n_a_cats)
-+{
-+ int i, band_id;
++ scatt.Plot(scatts = scatt_dt, ellipses = ellipses_dt, styles = cats_attrs)
+
-+ for(i = 0; i < n_a_bands; i++)
-+ {
-+ band_id = bands_ids[i];
-+ if(band_id >= 0) {
-+ Rast_close(fd_bands[i]);
-+ G_free(bands_rows[band_id].row);
-+ G_free(bands_rows[band_id].null_row);
-+ }
-+ }
-+
-+ if(f_cats_rasts_conds)
-+ for(i = 0; i < n_a_cats; i++)
-+ if(f_cats_rasts_conds[i])
-+ fclose(f_cats_rasts_conds[i]);
++ def _getSelectedAreas(self, scatt_id, scatt_dt, cats_attrs):
+
-+ if(fd_cats_rasts)
-+ for(i = 0; i < n_a_cats; i++)
-+ if(fd_cats_rasts[i] >= 0)
-+ Rast_close(fd_cats_rasts[i]);
++ cat_id = self.cats_mgr.GetSelectedCat()
++ if not cat_id:
++ return
+
-+}
++ sel_a_cat_id = max(cats_attrs.keys()) + 1
+
-+/*!
-+ \brief Compute scatter plots data.
++ s = self.scatt_conds_dt.GetScatt(scatt_id, [cat_id])
++ if not s:
++ return
++ cats_attrs[sel_a_cat_id] = {'color' : "255:255:0",
++ 'opacity' : 0.7,
++ 'show' : True}
+
-+ If category has not defined no category raster condition file and no scatter plot with consdtion,
-+ default scatter plot is computed.
++ scatt_dt[sel_a_cat_id] = s[cat_id]
+
-+ \param region analysis region, beaware that all input data must be prepared for this region (bands (their ranges), cats_rasts_conds rasters...)
-+ \param region function calls Rast_set_window for this region
-+ \param scatt_conds pointer to scScatts struct of type SC_SCATT_CONDITIONS, where are stored selected areas (conditions) in scatter plots
-+ \param cats_rasts_conds paths to category raster conditions files representing selected areas (conditions) in rasters for every category
-+ \param cats_rasts_conds index in array represents corresponding category id
-+ \param cats_rasts_conds for manupulation with category raster conditions file see also I_id_scatt_to_bands() and I_insert_patch_to_cat_rast()
-+ \param bands names of analyzed bands, order of bands is defined by their id
-+ \param n_bands number of bands
-+ \param [out] scatts pointer to scScatts struct of type SC_SCATT_DATA, where are computed scatter plots stored
-+ \param [out] cats_rasts array of raster maps names where will be stored all selected pixels for every category
++ def AddScattPlotDone(self, event):
++
++ transpose = event.userdata['transpose']
++ scatt_id = event.kwds['scatt_id']
+
-+ \return 0 on success
-+ \return -1 on failure
-+*/
-+int I_compute_scatts(struct Cell_head *region, struct scCats * scatt_conds, const char ** cats_rasts_conds,
-+ const char ** bands, int n_bands, struct scCats * scatts, const char ** cats_rasts)
-+{
-+ const char *mapset;
-+ char header[1024];
++ #TODO guiparent - not very good
++ self.plots[scatt_id] = self.guiparent.NewScatterPlot(scatt_id = scatt_id, transpose = transpose)
+
-+ int fd_cats_rasts[scatt_conds->n_a_cats];
-+ FILE * f_cats_rasts_conds[scatt_conds->n_a_cats];
++ self.plots[scatt_id].plotClosed.connect(self.PlotClosed)
+
-+ struct rast_row bands_rows[n_bands];
++ if self.plot_mode:
++ self.plots[scatt_id].SetMode(self.plot_mode)
++ self.plots[scatt_id].ZoomToExtend()
+
-+ RASTER_MAP_TYPE data_type;
++ self.RenderScattPlts(scatt_ids = [scatt_id])
+
-+ int nrows, i_band, n_a_bands, band_id,
-+ i, i_row, head_nchars, i_cat, id_cat;
-+
-+ int fd_bands[n_bands];
-+ int bands_ids[n_bands];
-+ int b_needed_bands[n_bands];
++ def PlotClosed(self, scatt_id):
++ del self.plots[scatt_id]
+
-+ Rast_set_window(region);
++ def SetPlotsMode(self, mode):
+
-+ for(i_band = 0; i_band < n_bands; i_band++)
-+ fd_bands[i_band] = -1;
++ self.plot_mode = mode
++ for scatt in self.plots.itervalues():
++ scatt.SetMode(mode)
++
++ self.modeSet.emit(mode = mode)
++
++ def ActivateSelectionPolygonMode(self, activate):
++ self.pol_sel_mode[0] = activate
++ self.RenderScattPlts()
++ return activate
++
++ def ProcessSelectionPolygons(self, process_mode):
++ scatts_polygons = {}
++ for scatt_id, scatt in self.plots.iteritems():
++ coords = scatt.GetCoords()
++ if coords is not None:
++ scatts_polygons[scatt_id] = coords
++
++ if not scatts_polygons:
++ return
++
++ value = 1
++ if process_mode == 'remove':
++ value = 0
++
++ sel_cat_id = self.cats_mgr.GetSelectedCat()
++ if not sel_cat_id:
++ dlg = wx.MessageDialog(parent = self.guiparent,
++ message = _("In order to select arrea in scatter plot, "
++ "you have to select class first.\n\n"
++ "There is no class yet, "
++ "do you want to create one?"),
++ caption = _("No class selected"),
++ style = wx.YES_NO)
++ if dlg.ShowModal() == wx.ID_YES:
++ self.iclass_conn.EmptyCategories()
++
++ sel_cat_id = self.cats_mgr.GetSelectedCat()
++ if not sel_cat_id:
++ return
++
++ for scatt in self.plots.itervalues():
++ scatt.SetEmpty()
++
++ self.tasks_pids['set_edit_cat_data'] = self.thread.GetId()
++ self.thread.Run(callable = self.core.UpdateCategoryWithPolygons,
++ cat_id = sel_cat_id,
++ scatts_pols = scatts_polygons,
++ value = value)
++
++ def SetPlotsEditingMode(self, mode):
++
++ self.pol_sel_mode[1] = mode
++ for scatt in self.plots.itervalues():
++ scatt.SetEditingMode(mode)
+
-+ for(i_band = 0; i_band < n_bands; i_band++)
-+ bands_ids[i_band] = -1;
++ def SetEditCatDataDone(self, event):
+
-+ if (n_bands != scatts->n_bands ||
-+ n_bands != scatt_conds->n_bands)
-+ return -1;
++ if event.exception:
++ GError(_("Error occured during computation of scatter plot category:\n%s"),
++ parent = self.guiparent, showTraceback = False)
+
-+ memset(b_needed_bands, 0, (size_t)n_bands * sizeof(int));
++ cat_id = event.ret
+
-+ get_needed_bands(scatt_conds, &b_needed_bands[0]);
-+ get_needed_bands(scatts, &b_needed_bands[0]);
++ self.RenderScattPlts()
+
-+ n_a_bands = 0;
++ cat_id = event.kwds["cat_id"]
+
-+ /* open band rasters, which are needed for computation */
-+ for(band_id = 0; band_id < n_bands; band_id++)
-+ {
-+ if(b_needed_bands[band_id])
-+ {
-+ G_debug(3, "Opening raster no. %d with name: %s", band_id, bands[band_id]);
++ cat_rast = self.core.GetCatRast(cat_id)
+
-+ if ((mapset = G_find_raster2(bands[band_id],"")) == NULL) {
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
-+ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
-+ G_warning(_("Unbale to read find raster <%s>"), bands[band_id]);
-+ return -1;
-+ }
-+
-+ if ((fd_bands[n_a_bands] = Rast_open_old(bands[band_id], mapset)) < 0) {
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
-+ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
-+ G_warning(_("Unbale to open raster <%s>"), bands[band_id]);
-+ return -1;
-+ }
++ if cat_rast not in self.added_cats_rasts.values():
+
-+ data_type = Rast_get_map_type(fd_bands[n_a_bands]);
-+ if(data_type != CELL_TYPE) {
-+ G_warning(_("Raster <%s> type is not <%s>"), bands[band_id], "CELL");
-+ return -1;
-+ }
++ cats_attrs = self.cats_mgr.GetCategoryAttrs(cat_id)
+
-+ bands_rows[band_id].row = Rast_allocate_c_buf();
-+ bands_rows[band_id].null_row = Rast_allocate_null_buf();
++
++ region = self.core.GetRegion()
++ ret, err_msg = RunCommand('r.region',
++ map = cat_rast,
++ getErrorMsg = True,
++ n = "%f" % region['n'],
++ s = "%f" % region['s'],
++ e = "%f" % region['e'],
++ w = "%f" % region['w'],
++ )
+
-+ if(Rast_read_range(bands[band_id], mapset, &bands_rows[band_id].rast_range) != 1){
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands,
-+ bands_ids, NULL, NULL, scatt_conds->n_a_cats);
-+ G_warning(_("Unable to read range of raster <%s>"), bands[band_id]);
-+ return -1;
-+ }
++ ret, err_msg = RunCommand('r.colors',
++ map = cat_rast,
++ rules = "-",
++ stdin = "1 %s" % cats_attrs["color"],
++ getErrorMsg = True)
+
-+ bands_ids[n_a_bands] = band_id;
-+ ++n_a_bands;
-+ }
-+ }
++ if ret != 0:
++ GError(_("r.region failed\n%s" % err_msg))
+
-+ /* open category rasters condition files and category rasters */
-+ for(i_cat = 0; i_cat < scatts->n_a_cats; i_cat++)
-+ {
-+ id_cat = scatts->cats_ids[i_cat];
-+ if(cats_rasts[id_cat]) {
-+ fd_cats_rasts[i_cat] = Rast_open_new(cats_rasts[id_cat], CELL_TYPE);
-+ }
-+ else
-+ fd_cats_rasts[i_cat] = -1;
++ self.mapWin.Map.AddLayer(ltype = "raster", name = "cat_%d" % cat_id, render = True,
++ command = ["d.rast", "map=%s" % cat_rast, "values=1"])
+
-+ if(cats_rasts_conds[id_cat]) {
-+ f_cats_rasts_conds[i_cat] = fopen(cats_rasts_conds[id_cat], "r");
-+ if(!f_cats_rasts_conds[i_cat]) {
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids,
-+ f_cats_rasts_conds, f_cats_rasts_conds, scatt_conds->n_a_cats);
-+ G_warning(_("Unable to open category raster condtition file <%s>"), bands[band_id]);
-+ return -1;
-+ }
-+ }
-+ else
-+ f_cats_rasts_conds[i_cat] = NULL;
-+ }
+
-+ head_nchars = get_cat_rast_header(region, header);
-+ for(i_cat = 0; i_cat < scatt_conds->n_a_cats; i_cat++)
-+ if(f_cats_rasts_conds[i_cat])
-+ if( fseek(f_cats_rasts_conds[i_cat] , head_nchars, SEEK_SET) != 0) {
-+ G_warning(_("Corrupted category raster conditions file (fseek failed)"));
-+ return -1;
-+ }
++ if ret != 0:
++ GError(_("r.region failed\n%s" % err_msg))
+
-+ nrows = Rast_window_rows();
++ self.added_cats_rasts[cat_id] = cat_rast
+
-+ /* analyze bands by rows */
-+ for (i_row = 0; i_row < nrows; i_row++)
-+ {
-+ for(i_band = 0; i_band < n_a_bands; i_band++)
-+ {
-+ band_id = bands_ids[i_band];
-+ Rast_get_c_row(fd_bands[i_band], bands_rows[band_id].row, i_row);
-+ Rast_get_null_value_row (fd_bands[i_band], bands_rows[band_id].null_row, i_row);
-+ }
-+ if(compute_scatts_from_chunk_row(region, scatt_conds, f_cats_rasts_conds, bands_rows, scatts, fd_cats_rasts) == -1)
-+ {
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids, fd_cats_rasts,
-+ f_cats_rasts_conds, scatt_conds->n_a_cats);
-+ return -1;
-+ }
-+
-+ }
-+ free_compute_scatts_data(fd_bands, bands_rows, n_a_bands, bands_ids,
-+ fd_cats_rasts, f_cats_rasts_conds, scatt_conds->n_a_cats);
-+ return 0;
-+}
-\ No newline at end of file
-Index: lib/imagery/scatt_sccats.c
-===================================================================
---- lib/imagery/scatt_sccats.c (revision 0)
-+++ lib/imagery/scatt_sccats.c (working copy)
-@@ -0,0 +1,405 @@
-+/*!
-+ \file lib/imagery/scatt_cat_rast.c
++ self.giface.updateMap.emit()
+
-+ \brief Imagery library - functions for manipulation with scatter plot structs.
++ def DigitDataChanged(self, vectMap, digit):
++
++ if self.mapWin_conn:
++ self.mapWin_conn.DigitDataChanged(vectMap, digit)
++ return 1
++ else:
++ return 0
+
-+ Copyright (C) 2013 by the GRASS Development Team
++ def GetCategoriesManager(self):
++ return self.cats_mgr
+
-+ This program is free software under the GNU General Public License
-+ (>=v2). Read the file COPYING that comes with GRASS for details.
+
-+ \author Stepan Turek <stepan.turek at seznam.cz> (Mentor: Martin Landa)
-+ */
++class CategoriesManager:
+
-+#include <grass/raster.h>
-+#include <grass/imagery.h>
-+#include <grass/gis.h>
++ def __init__(self, scatt_mgr, core):
+
-+#include <stdio.h>
-+#include <stdlib.h>
-+#include <math.h>
-+#include <string.h>
++ self.core = core
++ self.scatt_mgr = scatt_mgr
+
-+/*!
-+ \brief Compute band ids from scatter plot id.
++ self.cats = {}
++ self.cats_ids = []
+
-+ Scatter plot id describes which bands defines the scatter plot.
++ self.sel_cat_id = None
+
-+ Let say we have 3 bands, their ids are 0, 1 and 2.
-+ Scatter plot with id 0 consists of band 1 (b_1_id) 0 and band 2 (b_2_id) 1.
-+ All scatter plots:
-+ scatt_id b_1_id b_2_id
-+ 0 0 1
-+ 1 0 2
-+ 2 1 2
++ self.initialized = Signal('CategoriesManager.initialized')
++ self.setCategoryAttrs = Signal('CategoriesManager.setCategoryAttrs')
++ self.deletedCategory = Signal('CategoriesManager.deletedCategory')
++ self.addedCategory = Signal('CategoriesManager.addedCategory')
+
-+ \param scatt_id scatter plot id
-+ \param n_bands number of bands
-+ \param [out] b_1_id id of band1
-+ \param[out] b_2_id id of band2
++ def Clear(self):
+
-+ \return 0
-+ */
-+int I_id_scatt_to_bands(const int scatt_id, const int n_bands, int * b_1_id, int * b_2_id)
-+{
-+ int n_b1 = n_bands - 1;
++ self.cats.clear()
++ del self.cats_ids[:]
+
-+ * b_1_id = (int) ((2 * n_b1 + 1 - sqrt((double)((2 * n_b1 + 1) * (2 * n_b1 + 1) - 8 * scatt_id))) / 2);
++ self.sel_cat_id = None
+
-+ * b_2_id = scatt_id - ((* b_1_id) * (2 * n_b1 + 1) - (* b_1_id) * (* b_1_id)) / 2 + (* b_1_id) + 1;
++ def InitCoreCats(self):
++ if self.scatt_mgr.data_set:
++ for cat_id in self.cats_ids:
++ self.core.AddCategory(cat_id)
+
-+ return 0;
-+}
++ def AddCategory(self, cat_id = None, name = None, color = None):
+
++ if cat_id is None:
++ if self.cats_ids:
++ cat_id = max(self.cats_ids) + 1
++ else:
++ cat_id = 1
+
-+/*!
-+ \brief Compute scatter plot id from band ids.
++ if self.scatt_mgr.data_set:
++ self.scatt_mgr.thread.Run(callable = self.core.AddCategory,
++ cat_id = cat_id)
++ #TODO check number of cats
++ #if ret < 0: #TODO
++ # return -1;
+
-+ See also I_id_scatt_to_bands().
++ self.cats[cat_id] = {
++ 'name' : _('Category %s' % cat_id ),
++ 'color' : "0:0:0",
++ 'opacity' : 1.0,
++ 'show' : True
++ }
+
-+ \param n_bands number of bands
-+ \param b_1_id id of band1
-+ \param b_1_id id of band2
-+ \param [out] scatt_id scatter plot id
++ self.cats_ids.append(cat_id)
+
-+ \return 0
-+ */
-+int I_bands_to_id_scatt(const int b_1_id, const int b_2_id, const int n_bands, int * scatt_id)
-+{
-+ int n_b1 = n_bands - 1;
++ if name is not None:
++ self.cats[cat_id]["name"] = name
++
++ if color is not None:
++ self.cats[cat_id]["color"] = color
+
-+ * scatt_id = (b_1_id * (2 * n_b1 + 1) - b_1_id * b_1_id) / 2 + b_2_id - b_1_id - 1;
++ self.addedCategory.emit(cat_id = cat_id,
++ name = self.cats[cat_id]["name"],
++ color = self.cats[cat_id]["color"] )
++ return cat_id
+
-+ return 0;
-+}
++ def SetCategoryAttrs(self, cat_id, attrs_dict):
++ render = False
++ for k, v in attrs_dict.iteritems():
++ if not render and k in ['name', 'color', 'opacity', 'show']:
++ render = True
+
-+/*!
-+ \brief Initialize structure for storing scatter plots data.
++ self.cats[cat_id][k] = v
+
-+ \param cats pointer to scCats struct
-+ \param n_bands number of bands
-+ \param type SC_SCATT_DATA - stores scatter plots
-+ \param type SC_SCATT_CONDITIONS - stores selected areas in scatter plots
-+ */
-+void I_sc_init_cats(struct scCats * cats, int n_bands, int type)
-+{
-+ int i_cat;
++ #TODO optimization
++ if render:
++ self.scatt_mgr.RenderScattPlts()
+
-+ cats->type = type;
++ self.setCategoryAttrs.emit(cat_id = cat_id, attrs_dict = attrs_dict)
+
-+ cats->n_cats = 100;
-+ cats->n_a_cats = 0;
++ def DeleteCategory(self, cat_id):
+
-+ cats->n_bands = n_bands;
-+ cats->n_scatts = (n_bands - 1) * n_bands / 2;
++ if self.scatt_mgr.data_set:
++ self.scatt_mgr.thread.Run(callable = self.core.DeleteCategory,
++ cat_id = cat_id)
++ del self.cats[cat_id]
++ self.cats_ids.remove(cat_id)
+
-+ cats->cats_arr = (struct scScatts **) G_malloc(cats->n_cats * sizeof(struct scScatts *));
-+ memset(cats->cats_arr, 0, cats-> n_cats * sizeof(struct scScatts *));
++ self.deletedCategory.emit(cat_id = cat_id)
+
-+ cats->cats_ids = (int *) G_malloc(cats->n_cats * sizeof(int));
-+ cats->cats_idxs =(int *) G_malloc(cats->n_cats * sizeof(int));
++ #TODO emit event?
++ def SetSelectedCat(self, cat_id):
++ self.sel_cat_id = cat_id
++ if self.scatt_mgr.pol_sel_mode[0]:
++ self.scatt_mgr.RenderScattPlts()
+
-+ for(i_cat = 0; i_cat < cats->n_cats; i_cat++)
-+ cats->cats_idxs[i_cat] = -1;
++ def GetSelectedCat(self):
++ return self.sel_cat_id
+
-+ return;
-+}
++ def GetCategoryAttrs(self, cat_id):
++ #TODO is mutable
++ return self.cats[cat_id]
+
-+/*!
-+ \brief Free data of struct scCats, the structure itself remains alocated.
++ def GetCategoriesAttrs(self):
++ #TODO is mutable
++ return self.cats
++
++ def GetCategories(self):
++ return self.cats_ids[:]
+
-+ \param cats pointer to existing scCats struct
-+ */
-+void I_sc_free_cats(struct scCats * cats)
-+{
-+ int i_cat;
++ def SetCategoryPosition(self):
++ if newindex > oldindex:
++ newindex -= 1
++
++ self.cats_ids.insert(newindex, self.cats_ids.pop(oldindex))
+
-+ for(i_cat = 0; i_cat < cats->n_a_cats; i_cat++)
-+ {
-+ if(cats->cats_arr[i_cat])
-+ {
-+ G_free(cats->cats_arr[i_cat]->scatt_idxs);
-+ G_free(cats->cats_arr[i_cat]->scatts_bands);
-+ G_free(cats->cats_arr[i_cat]->scatts_arr);
-+ G_free(cats->cats_arr[i_cat]);
-+ }
-+ }
++class MapWinConnection:
++ def __init__(self, scatt_mgr, mapWin, scatt_rast_updater):
++ self.mapWin = mapWin
++ self.vectMap = None
++ self.scatt_rast_updater = scatt_rast_updater
++ self.scatt_mgr = scatt_mgr
++ self.cats_mgr = scatt_mgr.cats_mgr
+
-+ G_free(cats->cats_ids);
-+ G_free(cats->cats_idxs);
-+ G_free(cats->cats_arr);
++ self.thread = self.scatt_mgr.thread
+
-+ cats->n_cats = 0;
-+ cats->n_a_cats = 0;
-+ cats->n_bands = 0;
-+ cats->n_scatts = 0;
-+ cats->type = -1;
++ #TODO
++ self.mapWin.parent.toolbars["vdigit"].editingStarted.connect(self.DigitDataChanged)
+
-+ return;
-+}
++ #def ChangeMap(self, vectMap, layers_cats):
++ # self.vectMap = vectMap
++ # self.layers_cats = layers_cats
+
-+#if 0
-+void I_sc_get_active_categories(int * a_cats_ids, int * n_a_cats, struct scCats * cats)
-+{
-+ a_cats_ids = cats->cats_ids;
-+ * n_a_cats = cats->n_a_cats;
-+}
-+#endif
++ #ret, region, msg = RunCommand("v.to.rast",
++ # flags = "gp",
++ # getErrorMsg = True,
++ # read = True)
+
-+/*!
-+ \brief Add category.
-+
-+ Category represents group of scatter plots.
++ def _connectSignals(self):
++ self.digit.featureAdded.connect(self.AddFeature)
++ self.digit.areasDeleted.connect(self.DeleteAreas)
++ self.digit.featuresDeleted.connect(self.DeleteAreas)
++ self.digit.vertexMoved.connect(self.EditedFeature)
++ self.digit.vertexRemoved.connect(self.EditedFeature)
++ self.digit.lineEdited.connect(self.EditedFeature)
++ self.digit.featuresMoved.connect(self.EditedFeature)
+
-+ \param cats pointer to scCats struct
++ def AddFeature(self, new_bboxs, new_areas_cats):
++ if not self.scatt_mgr.data_set:
++ return
+
-+ \return assigned category id (starts with 0)
-+ \return -1 if maximum nuber of categories was reached
-+ */
-+int I_sc_add_cat(struct scCats * cats)
-+{
-+ int i_scatt, i_cat_id, cat_id;
-+ int n_a_cats = cats->n_a_cats;
++ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
++ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
++ new_bboxs = new_bboxs,
++ old_bboxs = [],
++ old_areas_cats = [],
++ new_areas_cats = new_areas_cats)
+
-+ if(cats->n_a_cats >= cats->n_cats)
-+ return -1;
++ def DeleteAreas(self, old_bboxs, old_areas_cats):
++ if not self.scatt_mgr.data_set:
++ return
+
-+ for(i_cat_id = 0; i_cat_id < cats->n_cats; i_cat_id++)
-+ if(cats->cats_idxs[i_cat_id] < 0) {
-+ cat_id = i_cat_id;
-+ break;
-+ }
++ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
++ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
++ new_bboxs = [],
++ old_bboxs = old_bboxs,
++ old_areas_cats = old_areas_cats,
++ new_areas_cats = [])
+
-+ cats->cats_ids[n_a_cats] = cat_id;
-+ cats->cats_idxs[cat_id] = n_a_cats;
-+
-+ cats->cats_arr[n_a_cats] = (struct scScatts *) G_malloc(sizeof(struct scScatts));
+
-+ cats->cats_arr[n_a_cats]->scatts_arr = (struct scdScattData **) G_malloc(cats->n_scatts * sizeof(struct scdScattData *));
-+ memset((cats->cats_arr[n_a_cats]->scatts_arr), 0, cats->n_scatts * sizeof(struct scdScattData *));
-+
-+ cats->cats_arr[n_a_cats]->n_a_scatts = 0;
++ def EditedFeature(self, new_bboxs, new_areas_cats, old_bboxs, old_areas_cats):
++ if not self.scatt_mgr.data_set:
++ return
+
-+ cats->cats_arr[n_a_cats]->scatts_bands = (int *) G_malloc(cats->n_scatts * 2 * sizeof(int));
++ self.scatt_mgr.tasks_pids['mapwin_conn'].append(self.thread.GetId())
++ self.thread.Run(callable = self.scatt_rast_updater.EditedFeature,
++ new_bboxs = new_bboxs,
++ old_bboxs = old_bboxs,
++ old_areas_cats = old_areas_cats,
++ new_areas_cats = new_areas_cats)
++
++ def DigitDataChanged(self, vectMap, digit):
++
++ self.digit = digit
++ self.vectMap = vectMap
++
++ self.digit.EmitSignals(emit = True)
++
++ self.scatt_rast_updater.SetVectMap(vectMap)
++
++ self._connectSignals()
++
++
++class IClassConnection:
++ def __init__(self, scatt_mgr, iclass_frame, cats_mgr):
++ self.iclass_frame = iclass_frame
++ self.stats_data = self.iclass_frame.stats_data
++ self.cats_mgr = cats_mgr
++ self.scatt_mgr = scatt_mgr
++
++ self.stats_data.statisticsAdded.connect(self.AddCategory)
++ self.stats_data.statisticsDeleted.connect(self.DeleteCategory)
++ self.stats_data.allStatisticsDeleted.connect(self.DeletAllCategories)
++ self.stats_data.statisticsSet.connect(self.SetCategory)
++
++ self.iclass_frame.groupSet.connect(self.GroupSet)
++
++ self.cats_mgr.setCategoryAttrs.connect(self.SetStatistics)
++ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
++ self.cats_mgr.addedCategory.connect(self.AddStatistics)
++
++ self.iclass_frame.categoryChanged.connect(self.CategoryChanged)
++
++ self.SyncCats()
++
++ def SetData(self):
++ self.iclass_frame.AddBands()
++
++ def EmptyCategories(self):
++ self.iclass_frame.OnCategoryManager(None)
++
++ def SyncCats(self):
++ self.cats_mgr.addedCategory.disconnect(self.AddStatistics)
++ cats = self.stats_data.GetCategories()
++ for c in cats:
++ stats = self.stats_data.GetStatistics(c)
++ self.cats_mgr.AddCategory(c, stats.name, stats.color)
++ self.cats_mgr.addedCategory.connect(self.AddStatistics)
++
++ def CategoryChanged(self, cat):
++ self.cats_mgr.SetSelectedCat(cat)
++
++ def AddCategory(self, cat, name, color):
++ self.cats_mgr.addedCategory.disconnect(self.AddStatistics)
++ self.cats_mgr.AddCategory(cat_id = cat, name = name, color = color)
++ self.cats_mgr.addedCategory.connect(self.AddStatistics)
++
++ def DeleteCategory(self, cat):
++ self.cats_mgr.deletedCategory.disconnect(self.DeleteStatistics)
++ self.cats_mgr.DeleteCategory(cat)
++ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
++
++ def DeletAllCategories(self):
++
++ self.cats_mgr.deletedCategory.disconnect(self.DeleteStatistics)
++ cats = self.stats_data.GetCategories()
++ for c in cats:
++ self.cats_mgr.DeleteCategory(c)
++ self.cats_mgr.deletedCategory.connect(self.DeleteStatistics)
++
++ def SetCategory(self, cat, stats):
++
++ self.cats_mgr.setCategoryAttrs.disconnect(self.SetStatistics)
++ cats_attr = {}
++ for attr in ['name', 'color']:
++ if stats.has_key(attr):
++ cats_attr[attr] = stats[attr]
++
++ if cats_attr:
++ self.cats_mgr.SetCategoryAttrs(cat, cats_attr)
++ self.cats_mgr.setCategoryAttrs.connect(self.SetStatistics)
++
++
++ def SetStatistics(self, cat_id, attrs_dict):
++ self.stats_data.statisticsSet.disconnect(self.SetCategory)
++ self.stats_data.GetStatistics(cat_id).SetStatistics(attrs_dict)
++ self.stats_data.statisticsSet.connect(self.SetCategory)
++
++ def AddStatistics(self, cat_id, name, color):
++ self.stats_data.statisticsAdded.disconnect(self.AddCategory)
++ self.stats_data.AddStatistics(cat_id, name, color)
++ self.stats_data.statisticsAdded.connect(self.AddCategory)
++
++ def DeleteStatistics(self, cat_id):
++ self.stats_data.statisticsDeleted.disconnect(self.DeleteCategory)
++ self.stats_data.DeleteStatistics(cat_id)
++ self.stats_data.statisticsDeleted.connect(self.DeleteCategory)
++
++ def GroupSet(self, group):
++ res = RunCommand('i.group',
++ flags = 'g',
++ group = group, subgroup = group,
++ read = True).strip()
++ if res.split('\n')[0]:
++ bands = res.split('\n')
++ self.scatt_mgr.SetData(bands)
+\ No newline at end of file
+Index: gui/wxpython/Makefile
+===================================================================
+--- gui/wxpython/Makefile (revision 57599)
++++ gui/wxpython/Makefile (working copy)
+@@ -13,7 +13,7 @@
+ $(wildcard animation/* core/*.py dbmgr/* gcp/*.py gmodeler/* \
+ gui_core/*.py iclass/* lmgr/*.py location_wizard/*.py mapwin/*.py mapdisp/*.py \
+ mapswipe/* modules/*.py nviz/*.py psmap/* rlisetup/* timeline/* vdigit/* \
+- vnet/*.py web_services/*.py wxplot/*.py) \
++ vnet/*.py web_services/*.py wxplot/*.py scatt_plot/*.py) \
+ gis_set.py gis_set_error.py wxgui.py README
+
+ DSTFILES := $(patsubst %,$(ETCDIR)/%,$(SRCFILES)) \
+@@ -21,8 +21,9 @@
+
+ PYDSTDIRS := $(patsubst %,$(ETCDIR)/%,animation core dbmgr gcp gmodeler \
+ gui_core iclass lmgr location_wizard mapwin mapdisp modules nviz psmap \
+- mapswipe vdigit wxplot web_services rlisetup vnet timeline)
++ mapswipe vdigit wxplot web_services rlisetup vnet timeline scatt_plot)
+
++
+ DSTDIRS := $(patsubst %,$(ETCDIR)/%,icons scripts xml)
+
+ default: $(DSTFILES)
+Index: gui/wxpython/vdigit/toolbars.py
+===================================================================
+--- gui/wxpython/vdigit/toolbars.py (revision 57599)
++++ gui/wxpython/vdigit/toolbars.py (working copy)
+@@ -17,6 +17,7 @@
+ import wx
+
+ from grass.script import core as grass
++from grass.pydispatch.signal import Signal
+
+ from gui_core.toolbars import BaseToolbar, BaseIcons
+ from gui_core.dialogs import CreateNewVector
+@@ -42,6 +43,8 @@
+ self.digit = None
+ self._giface = giface
+
++ self.editingStarted = Signal("VDigitToolbar.editingStarted")
++
+ # currently selected map layer for editing (reference to MapLayer instance)
+ self.mapLayer = None
+ # list of vector layers from Layer Manager (only in the current mapset)
+@@ -860,6 +863,7 @@
+ alpha = int(opacity * 255)
+ self.digit.GetDisplay().UpdateSettings(alpha = alpha)
+
++ self.editingStarted.emit(vectMap = mapLayer.GetName(), digit = self.digit)
+ return True
+
+ def StopEditing(self):
+Index: gui/wxpython/vdigit/wxdigit.py
+===================================================================
+--- gui/wxpython/vdigit/wxdigit.py (revision 57599)
++++ gui/wxpython/vdigit/wxdigit.py (working copy)
+@@ -17,7 +17,7 @@
+ (and NumPy would be an excellent candidate for acceleration via
+ e.g. OpenCL or CUDA; I'm surprised it hasn't happened already).
+
+-(C) 2007-2011 by the GRASS Development Team
++(C) 2007-2011, 2013 by the GRASS Development Team
+
+ This program is free software under the GNU General Public License
+ (>=v2). Read the file COPYING that comes with GRASS for details.
+@@ -27,6 +27,8 @@
+
+ import grass.script.core as grass
+
++from grass.pydispatch.signal import Signal
++
+ from core.gcmd import GError
+ from core.debug import Debug
+ from core.settings import UserSettings
+@@ -176,7 +178,21 @@
+
+ if self.poMapInfo:
+ self.InitCats()
+-
++
++ self.emit_signals = False
++
++ # signals which describes features changes during digitization,
++ # activate them using EmitSignals method
++ #TODO signal for errors?
++ self.featureAdded = Signal('IVDigit.featureAdded')
++ self.areasDeleted = Signal('IVDigit.areasDeleted')
++ self.vertexMoved = Signal('IVDigit.vertexMoved')
++ self.vertexAdded = Signal('IVDigit.vertexAdded')
++ self.vertexRemoved = Signal('IVDigit.vertexRemoved')
++ self.featuresDeleted = Signal('IVDigit.featuresDeleted')
++ self.featuresMoved = Signal('IVDigit.featuresMoved')
++ self.lineEdited = Signal('IVDigit.lineEdited')
++
+ def __del__(self):
+ Debug.msg(1, "IVDigit.__del__()")
+ Vect_destroy_line_struct(self.poPoints)
+@@ -188,7 +204,12 @@
+ Vect_close(self.poBgMapInfo)
+ self.poBgMapInfo = self.popoBgMapInfo = None
+ del self.bgMapInfo
+-
+
-+ cats->cats_arr[n_a_cats]->scatt_idxs = (int *) G_malloc(cats->n_scatts * sizeof(int));
-+ for(i_scatt = 0; i_scatt < cats->n_scatts; i_scatt++)
-+ cats->cats_arr[n_a_cats]->scatt_idxs[i_scatt] = -1;
-+
-+ ++cats->n_a_cats;
++ def EmitSignals(self, emit):
++ """!Activate/deactivate signals which describes features changes during digitization.
++ """
++ self.emit_signals = emit
+
-+ return cat_id;
-+}
+ def CloseBackgroundMap(self):
+ """!Close background vector map"""
+ if not self.poBgMapInfo:
+@@ -394,7 +415,6 @@
+
+ @return tuple (number of added features, feature ids)
+ """
+-
+ layer = self._getNewFeaturesLayer()
+ cat = self._getNewFeaturesCat()
+
+@@ -419,10 +439,14 @@
+ return (-1, None)
+
+ self.toolbar.EnableUndo()
+-
+- return self._addFeature(vtype, points, layer, cat,
+- self._getSnapMode(), self._display.GetThreshold())
+-
+
-+#if 0
-+int I_sc_delete_cat(struct scCats * cats, int cat_id)
-+{
-+ int cat_idx, i_cat;
++ ret = self._addFeature(vtype, points, layer, cat,
++ self._getSnapMode(), self._display.GetThreshold())
++ if ret[0] > -1 and self.emit_signals:
++ self.featureAdded.emit(new_bboxs = [self._createBbox(points)], new_areas_cats = [[{layer : [cat]}, None]])
+
-+ if(cat_id < 0 || cat_id >= cats->n_cats)
-+ return -1;
++ return ret
+
-+ cat_idx = cats->cats_idxs[cat_id];
-+ if(cat_idx < 0)
-+ return -1;
+ def DeleteSelectedLines(self):
+ """!Delete selected features
+
+@@ -434,16 +458,27 @@
+ # collect categories for deleting if requested
+ deleteRec = UserSettings.Get(group = 'vdigit', key = 'delRecord', subkey = 'enabled')
+ catDict = dict()
+
-+ G_free(cats->cats_arr[cat_idx]->scatt_idxs);
-+ G_free(cats->cats_arr[cat_idx]->scatts_bands);
-+ G_free(cats->cats_arr[cat_idx]->scatts_arr);
-+ G_free(cats->cats_arr[cat_idx]);
++ old_bboxs = []
++ old_areas_cats = []
+ if deleteRec:
+ for i in self._display.selected['ids']:
++
+ if Vect_read_line(self.poMapInfo, None, self.poCats, i) < 0:
+ self._error.ReadLine(i)
+
+- cats = self.poCats.contents
+- for j in range(cats.n_cats):
+- if cats.field[j] not in catDict.keys():
+- catDict[cats.field[j]] = list()
+- catDict[cats.field[j]].append(cats.cat[j])
++ if self.emit_signals:
++ ret = self._getLineAreaBboxCats(i)
++ if ret:
++ old_bboxs += ret[0]
++ old_areas_cats += ret[1]
++
++ # catDict was not used -> put into comment
++ #cats = self.poCats.contents
++ #for j in range(cats.n_cats):
++ # if cats.field[j] not in catDict.keys():
++ # catDict[cats.field[j]] = list()
++ # catDict[cats.field[j]].append(cats.cat[j])
+
+ poList = self._display.GetSelectedIList()
+ nlines = Vedit_delete_lines(self.poMapInfo, poList)
+@@ -456,7 +491,10 @@
+ self._deleteRecords(catDict)
+ self._addChangeset()
+ self.toolbar.EnableUndo()
+-
+
-+ for(i_cat = cat_idx; i_cat < cats->n_a_cats - 1; i_cat++)
-+ {
-+ cats->cats_arr[i_cat] = cats->cats_arr[i_cat + 1];
-+ cats->cats_ids[i_cat] = cats->cats_ids[i_cat + 1];
-+ }
-+ cats->cats_idxs[cat_id] = -1;
++ if self.emit_signals:
++ self.featuresDeleted.emit(old_bboxs = old_bboxs, old_areas_cats = old_areas_cats)
+
-+ --cats->n_a_cats;
-+
-+ return 0;
-+}
-+#endif
+ return nlines
+
+ def _deleteRecords(self, cats):
+@@ -512,22 +550,173 @@
+
+ @return number of deleted
+ """
++ if len(self._display.selected['ids']) < 1:
++ return 0
++
+ poList = self._display.GetSelectedIList()
+ cList = poList.contents
+
+ nareas = 0
++ old_bboxs = []
++ old_areas_cats = []
+
-+/*!
-+ \brief Insert scatter plot data .
-+ Inserted scatt_data struct must have same type as cats struct (SC_SCATT_DATA or SC_SCATT_CONDITIONS).
+ for i in range(cList.n_values):
+
-+ \param cats pointer to scCats struct
-+ \param cat_id id number of category.
-+ \param scatt_id id number of scatter plot.
+ if Vect_get_line_type(self.poMapInfo, cList.value[i]) != GV_CENTROID:
+ continue
+-
+
-+ \return 0 on success
-+ \return -1 on failure
-+ */
-+int I_sc_insert_scatt_data(struct scCats * cats, struct scdScattData * scatt_data, int cat_id, int scatt_id)
-+{
-+ int band_1, band_2, cat_idx, n_a_scatts;
-+ struct scScatts * scatts;
++ if self.emit_signals:
++ area = Vect_get_centroid_area(self.poMapInfo, cList.value[i]);
++ if area > 0:
++ bbox, cats = self._getaAreaBboxCats(area)
++ old_bboxs += bbox
++ old_areas_cats += cats
+
-+ if(cat_id < 0 || cat_id >= cats->n_cats)
-+ return -1;
+ nareas += Vedit_delete_area_centroid(self.poMapInfo, cList.value[i])
+
+ if nareas > 0:
+ self._addChangeset()
+ self.toolbar.EnableUndo()
++ if self.emit_signals:
++ self.areasDeleted.emit(old_bboxs = old_bboxs, old_areas_cats = old_areas_cats)
+
-+ cat_idx = cats->cats_idxs[cat_id];
-+ if(cat_idx < 0)
-+ return -1;
++ return nareas
++
++ def _getLineAreaBboxCats(self, ln_id):
++ ltype = Vect_read_line(self.poMapInfo, None, None, ln_id)
+
-+ if(scatt_id < 0 && scatt_id >= cats->n_scatts)
-+ return -1;
++ if ltype == GV_CENTROID:
++ #TODO centroid opttimization, can be adited also its area -> it will appear two times in new_ lists
++ return self._getCentroidAreaBboxCats(ln_id)
++ else:
++ return [self._getBbox(ln_id)], [self._getLineAreasCategories(ln_id)]
+
-+ scatts = cats->cats_arr[cat_idx];
-+ if(scatts->scatt_idxs[scatt_id] >= 0)
-+ return -1;
+
-+ if(!scatt_data->b_conds_arr && cats->type == SC_SCATT_CONDITIONS)
-+ return -1;
++ def _getCentroidAreaBboxCats(self, centroid):
++ if not Vect_line_alive(self.poMapInfo, centroid):
++ return None
+
-+ if(!scatt_data->scatt_vals_arr && cats->type == SC_SCATT_DATA)
-+ return -1;
++ area = Vect_get_centroid_area(self.poMapInfo, centroid)
++ if area > 0:
++ return self._getaAreaBboxCats(area)
++ else:
++ return None
+
-+ n_a_scatts = scatts->n_a_scatts;
++ def _getaAreaBboxCats(self, area):
+
-+ scatts->scatt_idxs[scatt_id] = n_a_scatts;
++ po_b_list = Vect_new_list()
++ Vect_get_area_boundaries(self.poMapInfo, area, po_b_list);
++ b_list = po_b_list.contents
+
-+ I_id_scatt_to_bands(scatt_id, cats->n_bands, &band_1, &band_2);
++ geoms = []
++ areas_cats = []
+
-+ scatts->scatts_bands[n_a_scatts * 2] = band_1;
-+ scatts->scatts_bands[n_a_scatts * 2 + 1] = band_2;
++ if b_list.n_values > 0:
++ for i_line in range(b_list.n_values):
+
-+ scatts->scatts_arr[n_a_scatts] = scatt_data;
-+ ++scatts->n_a_scatts;
++ line = b_list.value[i_line];
+
-+ return 0;
-+}
++ geoms.append(self._getBbox(abs(line)))
++ areas_cats.append(self._getLineAreasCategories(abs(line)))
+
+- return nareas
++ Vect_destroy_list(po_b_list);
+
-+#if 0
-+int I_sc_remove_scatt_data(struct scCats * cats, struct scdScattData * scatt_data, int cat_id, int scatt_id)
-+{
-+ int cat_idx, scatt_idx, n_init_scatts, i_scatt;
-+ struct scScatts * scatts;
++ return geoms, areas_cats
+
-+ if(cat_id < 0 && cat_id >= cats->n_cats)
-+ return -1;
++ def _getLineAreasCategories(self, ln_id):
++ if not Vect_line_alive (self.poMapInfo, ln_id):
++ return []
+
-+ cat_idx = cats->cats_idxs[cat_id];
-+ if(cat_idx < 0)
-+ return -1;
++ ltype = Vect_read_line(self.poMapInfo, None, None, ln_id)
++ if ltype != GV_BOUNDARY:
++ return []
+
-+ if(scatt_id < 0 || scatt_id >= cats->n_scatts)
-+ return -1;
++ cats = [None, None]
+
-+ scatts = cats->cats_arr[cat_idx];
-+ if(scatts->scatt_idxs[scatt_id] < 0)
-+ return -1;
++ left = c_int()
++ right = c_int()
+
-+ scatt_data = scatts->scatts_arr[scatt_idx];
++ if Vect_get_line_areas(self.poMapInfo, ln_id, pointer(left), pointer(right)) == 1:
++ areas = [left.value, right.value]
+
-+ for(i_scatt = scatt_idx; i_scatt < scatts->n_a_scatts - 1; i_scatt++)
-+ {
-+ scatts->scatts_arr[i_scatt] = scatts->scatts_arr[i_scatt + 1];
-+ scatts->scatts_bands[i_scatt * 2] = scatts->scatts_bands[(i_scatt + 1)* 2];
-+ scatts->scatts_bands[i_scatt * 2 + 1] = scatts->scatts_bands[(i_scatt + 1) * 2 + 1];
-+ }
-+ scatts->scatts_arr[scatts->n_a_scatts] = NULL;
++ for i, a in enumerate(areas):
++ if a > 0:
++ centroid = Vect_get_area_centroid(self.poMapInfo, a)
++ if centroid <= 0:
++ continue
++ c = self._getCategories(centroid)
++ if c:
++ cats[i] = c
+
-+ scatts->scatt_idxs[scatt_id] = -1;
++ return cats
+
-+ scatt_data = scatts->scatts_arr[scatt_id];
-+ scatts->n_a_scatts--;
++ def _getCategories(self, ln_id):
++ if not Vect_line_alive (self.poMapInfo, ln_id):
++ return none
+
-+ return 0;
-+}
++ poCats = Vect_new_cats_struct()
++ if Vect_read_line(self.poMapInfo, None, poCats, ln_id) < 0:
++ Vect_destroy_cats_struct(poCats)
++ return None
+
-+int I_sc_set_value(struct scCats * cats, int cat_id, int scatt_id, int value_idx, int value)
-+{
-+ int n_a_scatts = cats->cats_arr[cat_id]->n_a_scatts;
-+ int cat_idx, scatt_idx, ret;
++ cCats = poCats.contents
+
-+ cat_idx = cats->cats_idxs[cat_id];
-+ if(cat_idx < 0)
-+ return -1;
-+
-+ if(cats->cats_arr[cat_idx]->scatt_idxs[scatt_id] < 0)
-+ return -1;
++ cats = {}
++ for j in range(cCats.n_cats):
++ if cats.has_key(cCats.field[j]):
++ cats[cCats.field[j]].append(cCats.cat[j])
++ else:
++ cats[cCats.field[j]] = [cCats.cat[j]]
+
++ Vect_destroy_cats_struct(poCats)
++ return cats
+
-+ cat_idx = cats->cats_idxs[cat_id];
-+ scatt_idx = cats->cats_arr[cat_idx]->scatt_idxs[scatt_id];
++ def _getBbox(self, ln_id):
++ if not Vect_line_alive (self.poMapInfo, ln_id):
++ return None
+
-+ I_scd_set_value(cats->cats_arr[cat_idx]->scatts_arr[scatt_idx], value_idx, value);
++ poPoints = Vect_new_line_struct()
++ if Vect_read_line(self.poMapInfo, poPoints, None, ln_id) < 0:
++ Vect_destroy_line_struct(poPoints)
++ return []
+
-+ return 0;
-+}
-+#endif
++ geom = self._convertGeom(poPoints)
++ bbox = self._createBbox(geom)
++ Vect_destroy_line_struct(poPoints)
++ return bbox
+
-+/*!
-+ \brief Insert scatter plot data.
++ def _createBbox(self, points):
++
++ bbox = {}
++ for pt in points:
++ if not bbox.has_key('maxy'):
++ bbox['maxy'] = pt[1]
++ bbox['miny'] = pt[1]
++ bbox['maxx'] = pt[0]
++ bbox['minx'] = pt[0]
++ continue
++
++ if bbox['maxy'] < pt[1]:
++ bbox['maxy'] = pt[1]
++ elif bbox['miny'] > pt[1]:
++ bbox['miny'] = pt[1]
++
++ if bbox['maxx'] < pt[0]:
++ bbox['maxx'] = pt[0]
++ elif bbox['minx'] > pt[0]:
++ bbox['minx'] = pt[0]
++ return bbox
++
++ def _convertGeom(self, poPoints):
++
++ Points = poPoints.contents
++
++ pts_geom = []
++ for j in range(Points.n_points):
++ pts_geom.append((Points.x[j], Points.y[j]))
++
++ return pts_geom
++
+ def MoveSelectedLines(self, move):
+ """!Move selected features
+
+@@ -536,16 +725,45 @@
+ if not self._checkMap():
+ return -1
+
++ nsel = len(self._display.selected['ids'])
++ if nsel < 1:
++ return -1
++
+ thresh = self._display.GetThreshold()
+ snap = self._getSnapMode()
+
+ poList = self._display.GetSelectedIList()
++
++ if self.emit_signals:
++ old_bboxs = []
++ old_areas_cats = []
++ for sel_id in self._display.selected['ids']:
++ ret = self._getLineAreaBboxCats(sel_id)
++ if ret:
++ old_bboxs += ret[0]
++ old_areas_cats += ret[1]
++
++ Vect_set_updated(self.poMapInfo, 1)
++ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
++
+ nlines = Vedit_move_lines(self.poMapInfo, self.popoBgMapInfo, int(self.poBgMapInfo is not None),
+ poList,
+ move[0], move[1], 0,
+ snap, thresh)
++
+ Vect_destroy_list(poList)
+-
++
++ if nlines > 0 and self.emit_signals:
++ new_bboxs = []
++ new_areas_cats = []
++ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
++ for i in range(n_up_lines_old, n_up_lines):
++ new_id = Vect_get_updated_line(self.poMapInfo, i)
++ ret = self._getLineAreaBboxCats(new_id)
++ if ret:
++ new_bboxs += ret[0]
++ new_areas_cats += ret[1]
++
+ if nlines > 0 and self._settings['breakLines']:
+ for i in range(1, nlines):
+ self._breakLineAtIntersection(nlines + i, None, changeset)
+@@ -553,7 +771,13 @@
+ if nlines > 0:
+ self._addChangeset()
+ self.toolbar.EnableUndo()
+-
++
++ if self.emit_signals:
++ self.featuresMoved.emit(new_bboxs = new_bboxs,
++ old_bboxs = old_bboxs,
++ old_areas_cats = old_areas_cats,
++ new_areas_cats = new_areas_cats)
++
+ return nlines
+
+ def MoveSelectedVertex(self, point, move):
+@@ -571,12 +795,21 @@
+
+ if len(self._display.selected['ids']) != 1:
+ return -1
+-
++
++ # move only first found vertex in bbox
++ poList = self._display.GetSelectedIList()
++
++ if self.emit_signals:
++ cList = poList.contents
++ old_bboxs = [self._getBbox(cList.value[0])]
++ old_areas_cats = [self._getLineAreasCategories(cList.value[0])]
++
++ Vect_set_updated(self.poMapInfo, 1)
++ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
++
+ Vect_reset_line(self.poPoints)
+ Vect_append_point(self.poPoints, point[0], point[1], 0.0)
+-
+- # move only first found vertex in bbox
+- poList = self._display.GetSelectedIList()
++
+ moved = Vedit_move_vertex(self.poMapInfo, self.popoBgMapInfo, int(self.poBgMapInfo is not None),
+ poList, self.poPoints,
+ self._display.GetThreshold(type = 'selectThresh'),
+@@ -584,7 +817,17 @@
+ move[0], move[1], 0.0,
+ 1, self._getSnapMode())
+ Vect_destroy_list(poList)
+-
++
++ if moved > 0 and self.emit_signals:
++ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
++
++ new_bboxs = []
++ new_areas_cats = []
++ for i in range(n_up_lines_old, n_up_lines):
++ new_id = Vect_get_updated_line(self.poMapInfo, i)
++ new_bboxs.append(self._getBbox(new_id))
++ new_areas_cats.append(self._getLineAreasCategories(new_id))
++
+ if moved > 0 and self._settings['breakLines']:
+ self._breakLineAtIntersection(Vect_get_num_lines(self.poMapInfo),
+ None)
+@@ -592,7 +835,13 @@
+ if moved > 0:
+ self._addChangeset()
+ self.toolbar.EnableUndo()
+-
++
++ if self.emit_signals:
++ self.vertexMoved.emit(new_bboxs = new_bboxs,
++ new_areas_cats = new_areas_cats,
++ old_areas_cats = old_areas_cats,
++ old_bboxs = old_bboxs)
++
+ return moved
+
+ def AddVertex(self, coords):
+@@ -681,6 +930,10 @@
+ self._error.ReadLine(line)
+ return -1
+
++ if self.emit_signals:
++ old_bboxs = [self._getBbox(line)]
++ old_areas_cats = [self._getLineAreasCategories(line)]
++
+ # build feature geometry
+ Vect_reset_line(self.poPoints)
+ for p in coords:
+@@ -696,6 +949,9 @@
+
+ newline = Vect_rewrite_line(self.poMapInfo, line, ltype,
+ self.poPoints, self.poCats)
++ if newline > 0 and self.emit_signals:
++ new_geom = [self._getBbox(newline)]
++ new_areas_cats = [self._getLineAreasCategories(newline)]
+
+ if newline > 0 and self._settings['breakLines']:
+ self._breakLineAtIntersection(newline, None)
+@@ -703,7 +959,13 @@
+ if newline > 0:
+ self._addChangeset()
+ self.toolbar.EnableUndo()
+-
+
-+ \param scatt_data pointer to existing struct scdScattData
-+ \param type SC_SCATT_DATA for scatter plots or SC_SCATT_CONDITIONS for selected areas in scatter plot
-+ \param n_vals number of data values
-+ \param data array of values (unsigned char for SC_SCATT_CONDITIONS, unsigned int for SC_SCATT_DATA)
-+ */
-+void I_scd_init_scatt_data(struct scdScattData * scatt_data, int type, int n_vals, void * data)
-+{
-+ scatt_data->n_vals = n_vals;
++ if self.emit_signals:
++ self.lineEdited.emit(old_bboxs = old_bboxs,
++ old_areas_cats = old_areas_cats,
++ new_bboxs = new_bboxs,
++ new_areas_cats = new_areas_cats)
+
-+ if(type == SC_SCATT_DATA)
-+ {
-+ if(data)
-+ scatt_data->scatt_vals_arr = (unsigned int *) data;
-+ else {
-+ scatt_data->scatt_vals_arr = (unsigned int *) G_malloc(n_vals * sizeof(unsigned int));
-+ memset(scatt_data->scatt_vals_arr, 0, n_vals * sizeof(unsigned int));
-+ }
-+ scatt_data->b_conds_arr = NULL;
-+ }
-+ else if(type == SC_SCATT_CONDITIONS)
-+ {
-+ if(data)
-+ scatt_data->b_conds_arr = (unsigned char *) data;
-+ else {
-+ scatt_data->b_conds_arr = (unsigned char *) G_malloc(n_vals * sizeof(unsigned char));
-+ memset(scatt_data->b_conds_arr, 0, n_vals * sizeof(unsigned char));
-+ }
-+ scatt_data->scatt_vals_arr = NULL;
-+ }
+ return newline
+
+ def FlipLine(self):
+@@ -1514,6 +1776,16 @@
+ return 0
+
+ poList = self._display.GetSelectedIList()
+
-+ return;
-+}
++ if self.emit_signals:
++ cList = poList.contents
++
++ old_bboxs = [self._getBbox(cList.value[0])]
++ old_areas_cats = [self._getLineAreasCategories(cList.value[0])]
+
++ Vect_set_updated(self.poMapInfo, 1)
++ n_up_lines_old = Vect_get_num_updated_lines(self.poMapInfo)
+
-+#if 0
-+void I_scd_get_range_min_max(struct scdScattData * scatt_data, CELL * band_1_min, CELL * band_1_max, CELL * band_2_min, CELL * band_2_max)
-+{
+ Vect_reset_line(self.poPoints)
+ Vect_append_point(self.poPoints, coords[0], coords[1], 0.0)
+
+@@ -1525,15 +1797,35 @@
+ else:
+ ret = Vedit_remove_vertex(self.poMapInfo, poList,
+ self.poPoints, thresh)
+
-+ Rast_get_range_min_max(&(scatt_data->band_1_range), band_1_min, band_2_min);
-+ Rast_get_range_min_max(&(scatt_data->band_2_range), band_2_min, band_2_max);
+ Vect_destroy_list(poList)
+
-+ return;
-+}
-+s
-+void * I_scd_get_data_ptr(struct scdScattData * scatt_data)
-+{
-+ if(!scatt_data->b_conds_arr)
-+ return scatt_data->b_conds_arr;
-+ else if(!scatt_data->scatt_vals_arr)
-+ return scatt_data->scatt_vals_arr;
++ if ret > 0 and self.emit_signals:
++ new_bboxs = []
++ new_areas_cats = []
+
-+ return NULL;
-+}
++ n_up_lines = Vect_get_num_updated_lines(self.poMapInfo)
++ for i in range(n_up_lines_old, n_up_lines):
++ new_id = Vect_get_updated_line(self.poMapInfo, i)
++ new_areas_cats.append(self._getLineAreasCategories(new_id))
++ new_bboxs.append(self._getBbox(new_id))
+
+ if not add and ret > 0 and self._settings['breakLines']:
+ self._breakLineAtIntersection(Vect_get_num_lines(self.poMapInfo),
+ None)
+-
+
-+int I_scd_set_value(struct scdScattData * scatt_data, unsigned int val_idx, unsigned int val)
-+{
-+ if(val_idx < 0 && val_idx > scatt_data->n_vals)
-+ return -1;
+ if ret > 0:
+ self._addChangeset()
+-
+
-+ if(scatt_data->b_conds_arr)
-+ scatt_data->b_conds_arr[val_idx] = val;
-+ else if(scatt_data->scatt_vals_arr)
-+ scatt_data->scatt_vals_arr[val_idx] = val;
-+ else
-+ return -1;
++ if ret > 0 and self.emit_signals:
++ if add:
++ self.vertexAdded.emit(old_bboxs = old_bboxs, new_bboxs = new_bboxs)
++ else:
++ self.vertexRemoved.emit(old_bboxs = old_bboxs,
++ new_bboxs = new_bboxs,
++ old_areas_cats = old_areas_cats,
++ new_areas_cats = new_areas_cats)
+
-+ return 0;
-+}
-+#endif
-Index: lib/vector/Vlib/open.c
+ return 1
+
+ def GetLineCats(self, line):
+Index: gui/wxpython/mapdisp/toolbars.py
===================================================================
---- lib/vector/Vlib/open.c (revision 57581)
-+++ lib/vector/Vlib/open.c (working copy)
-@@ -240,7 +240,9 @@
- }
- else {
- char file_path[GPATH_MAX];
--
-+ /* reduce to current mapset if search path was set */
-+ if(strcmp(Map->mapset, "") == 0)
-+ Map->mapset = G_store(G_mapset());
- /* temporary map: reduce to current mapset if search path
- * was set */
- if (strcmp(Map->mapset, "") == 0)
+--- gui/wxpython/mapdisp/toolbars.py (revision 57599)
++++ gui/wxpython/mapdisp/toolbars.py (working copy)
+@@ -239,7 +239,8 @@
+ (MapIcons["scatter"], self.parent.OnScatterplot),
+ (MapIcons["histogram"], self.parent.OnHistogramPyPlot),
+ (BaseIcons["histogramD"], self.parent.OnHistogram),
+- (MapIcons["vnet"], self.parent.OnVNet)))
++ (MapIcons["vnet"], self.parent.OnVNet),
++ (MapIcons["scatter"], self.parent.OnScatterplot2)))
+
+ def OnDecoration(self, event):
+ """!Decorations overlay menu
+Index: gui/wxpython/mapdisp/frame.py
+===================================================================
+--- gui/wxpython/mapdisp/frame.py (revision 57599)
++++ gui/wxpython/mapdisp/frame.py (working copy)
+@@ -225,6 +225,7 @@
+ #
+ self.dialogs = {}
+ self.dialogs['attributes'] = None
++ self.dialogs['scatt_plot'] = None
+ self.dialogs['category'] = None
+ self.dialogs['barscale'] = None
+ self.dialogs['legend'] = None
+@@ -1168,6 +1169,19 @@
+ """!Returns toolbar with zooming tools"""
+ return self.toolbars['map']
+
++ def OnScatterplot2(self, event):
++ """!Init interactive scatterplot tools
++ """
++ if self.dialogs['scatt_plot']:
++ self.dialogs['scatt_plot'].Raise()
++ return
++
++ from scatt_plot.dialogs import ScattPlotMainDialog
++ self.dialogs['scatt_plot'] = ScattPlotMainDialog(parent=self, giface=self._giface)
++
++ self.dialogs['scatt_plot'].CenterOnScreen()
++ self.dialogs['scatt_plot'].Show()
++
+ def OnVNet(self, event):
+ """!Dialog for v.net* modules
+ """
+Index: gui/icons/grass/polygon.png
+===================================================================
+Cannot display: file marked as a binary type.
+svn:mime-type = application/octet-stream
+Index: gui/icons/grass/polygon.png
+===================================================================
+--- gui/icons/grass/polygon.png (revision 57599)
++++ gui/icons/grass/polygon.png (working copy)
+
+Property changes on: gui/icons/grass/polygon.png
+___________________________________________________________________
+Added: svn:mime-type
+## -0,0 +1 ##
++application/octet-stream
+\ No newline at end of property
More information about the grass-commit
mailing list