Opened 13 years ago

Last modified 12 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 Changed 13 years ago by hamish

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 Changed 12 years ago by neteler

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 Changed 12 years ago 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.

comment:4 Changed 12 years ago by hamish

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

Note: See TracTickets for help on using tickets.