15

I am trying to create a rotated rectangle in lat/lon coordinates.

Basically, I have the coordinates of the center and rotation angle around the center.

Is it possible to do using, shapely or gdal, or something else?

Germán Carrillo
  • 36,307
  • 5
  • 123
  • 178
Oleksandr
  • 253
  • 1
  • 2
  • 4

5 Answers5

14

Shapely's affinity module supports affine transformations on any geometry, for example:

from shapely import affinity
from shapely.geometry import LineString

# Example geometry
line = LineString([(1, 3), (1, 1), (4, 1)])

# Rotate 30 degrees CCW from origin at the center of bbox
line_rot_center = affinity.rotate(line, 30, 'center')

# Rotate from origin at (1, 1)
line_rot_11 = affinity.rotate(line, 30, (1, 1))
Mike T
  • 42,095
  • 10
  • 126
  • 187
  • Unfortunately, a rotation in (lat, lon) is not performed with an affine transformation. You either have to project the data and rotate the projected coordinates or else perform a spherical rotation. – whuber Aug 11 '11 at 15:05
  • 1
    Oh right, I missed the (lat, lon) part of the question. I've never seen a affine transformation function built for geographic inputs, so the tips you provided are essential. – Mike T Aug 11 '11 at 22:27
  • Since v1.2.17, Shapely has had affine transformations, although they're not documented yet - https://github.com/Toblerity/Shapely/blob/master/shapely/affinity.py – Alex L Apr 28 '13 at 07:49
  • Also see the author's blog: http://sgillies.net/blog/1168/shapely-1-2-17 – Alex L Apr 28 '13 at 07:50
  • 1
    @AlexL type help(shapely.affinity). Actually I'm nearly done preparing the manual update for the module. – Mike T Apr 28 '13 at 10:19
  • @MikeToews Ah, I missed the bit that said you wrote it. Great work - thank you! – Alex L Apr 28 '13 at 12:22
2

Sextante has a tool that might work. It is called Transform, under the "Tools for Vector layers". To run the command, define the angle, anchor point, scale factor and translation. Sextante is an extension to gvsig, openjump, udig, others. It also runs from the command-line.

klewis
  • 7,475
  • 17
  • 19
2

I don't see any features in Python or Shapely to do so but you can do this with JTS (Java Topology Suite) http://tsusiatsoftware.net/jts/main.html (GEOS, the lib used by Shapely is a port of JTS but JTS have more functions)

You can download the soft here http://sourceforge.net/projects/jts-topo-suite/ and launch testbuilder.sh or testbuilder.bat depending on your OS to play for discovering. After you understood, the basics you will have to do the same without GUI but with Java programming only.

http://geoscript.org/ may help you to make the glue beetween python code and java (because one implementation use Jython)

ThomasG77
  • 30,725
  • 1
  • 53
  • 93
2

You can use the gdal python bindings to do affine transformations (Link).

Since it is built on top of proj/geos/gdal you can do the same things as in shapely + more.

Kadir Şahbaz
  • 76,800
  • 56
  • 247
  • 389
Ragi Yaser Burhum
  • 15,339
  • 2
  • 59
  • 76
1

Grass can do it using v.transform.

lagerratrobe
  • 981
  • 6
  • 16