Opened 15 years ago

Last modified 14 years ago

#674 new enhancement

add a region option flag to v.proj / r.proj

Reported by: timmie Owned by: grass-dev@…
Priority: normal Milestone: 6.5.0
Component: Vector Version: svn-develbranch6
Keywords: v.proj Cc:
CPU: x86-32 Platform: Linux

Description (last modified by neteler)

  • 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

Change History (4)

comment:1 by hamish, 15 years ago

Component: defaultVector
Keywords: v.proj added
Milestone: 6.4.06.5.0
Type: defectenhancement

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

comment:2 by neteler, 14 years ago

Description: modified (diff)

See also Ticket #818 for assertion.

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

comment:3 by hamish, 14 years ago

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.

comment:4 by hamish, 14 years ago

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

Note: See TracTickets for help on using tickets.