[GRASS-dev] [GRASS GIS] #674: add a region option flag to v.proj / r.proj

#674: add a region option flag to v.proj / r.proj
---------------------+------------------------------------------------------
Reporter: timmie | Owner: grass-dev@lists.osgeo.org
     Type: defect | Status: new
Priority: normal | Milestone: 6.4.0
Component: default | Version: svn-develbranch6
Keywords: | Platform: Linux
      Cpu: x86-32 |
---------------------+------------------------------------------------------
* imagine a global data set (world borders in EPSG 4326) that is to be
reprojected into a samller zone (UTM N13 location).
* doing this without the region flag caused the following error:
- v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]' failed.
     * take -r flag of v.in.ogr as example

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/674&gt;
GRASS GIS <http://grass.osgeo.org>

#674: add a region option flag to v.proj / r.proj
--------------------------+-------------------------------------------------
  Reporter: timmie | Owner: grass-dev@lists.osgeo.org
      Type: enhancement | Status: new
  Priority: normal | Milestone: 6.5.0
Component: Vector | Version: svn-develbranch6
Resolution: | Keywords: v.proj
  Platform: Linux | Cpu: x86-32
--------------------------+-------------------------------------------------
Changes (by hamish):

  * keywords: => v.proj
  * type: defect => enhancement
  * component: default => Vector
  * milestone: 6.4.0 => 6.5.0

Comment:

It is not surprising it ended in error, you asked it to do something
mathematically unsuitable.

crop your vector with v.in.region + v.overlay before you try to project
it. simply selecting line features if they cross the target region is no
good as the line could extend far out of zone. simply throwing away
vertices outside the region is no good too because it breaks topology.

the error message could be clearer, beyond that I'd consider to close this
as wontfix.

Hamish

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/674#comment:1&gt;
GRASS GIS <http://grass.osgeo.org>

#674: add a region option flag to v.proj / r.proj
--------------------------+-------------------------------------------------
  Reporter: timmie | Owner: grass-dev@lists.osgeo.org
      Type: enhancement | Status: new
  Priority: normal | Milestone: 6.5.0
Component: Vector | Version: svn-develbranch6
Resolution: | Keywords: v.proj
  Platform: Linux | Cpu: x86-32
--------------------------+-------------------------------------------------
Old description:

* imagine a global data set (world borders in EPSG 4326) that is to be
reprojected into a samller zone (UTM N13 location).
* doing this without the region flag caused the following error:
- v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]' failed.
    * take -r flag of v.in.ogr as example

New description:

* imagine a global data set (world borders in EPSG 4326) that is to be
reprojected into a samller zone (UTM N13 location).
  * doing this without the region flag caused the following error:
{{{
v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]'
}}}
  failed.
  * take -r flag of v.in.ogr as example

Comment (by neteler):

See also Ticket #818 for assertion.

Do ticket #845 and changes in 6.5svn (r40559) and trunk in r40560) help?

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/674#comment:2&gt;
GRASS GIS <http://grass.osgeo.org>

#674: add a region option flag to v.proj / r.proj
--------------------------+-------------------------------------------------
  Reporter: timmie | Owner: grass-dev@lists.osgeo.org
      Type: enhancement | Status: new
  Priority: normal | Milestone: 6.5.0
Component: Vector | Version: svn-develbranch6
Resolution: | Keywords: v.proj
  Platform: Linux | Cpu: x86-32
--------------------------+-------------------------------------------------
Comment (by hamish):

to add to my previous comments to this ticket (which were perhaps a bit
harsh), a region cropping flag could be nice. even if it is not
universally workable it would probably work for most practical cases.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/674#comment:3&gt;
GRASS GIS <http://grass.osgeo.org>

#674: add a region option flag to v.proj / r.proj
--------------------------+-------------------------------------------------
  Reporter: timmie | Owner: grass-dev@lists.osgeo.org
      Type: enhancement | Status: new
  Priority: normal | Milestone: 6.5.0
Component: Vector | Version: svn-develbranch6
Resolution: | Keywords: v.proj
  Platform: Linux | Cpu: x86-32
--------------------------+-------------------------------------------------
Comment (by hamish):

(this only applies to v.proj btw, r.proj already crops)

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/674#comment:4&gt;
GRASS GIS <http://grass.osgeo.org>