StackOverflow2013

Note that there are some explanatory texts on larger screens.

plurals

PO
primarykey
Id
14521304
data
AcceptedAnswerId
0
AnswerCount
0
ClosedDate
CommentCount
0
CommunityOwnedDate
CreationDate
2013-01-25T11:54:33.613
FavoriteCount
0
LastActivityDate
2013-01-25T11:54:33.613
LastEditDate
LastEditorUserId
0
OwnerUserId
1257804
ParentId
9909945
PostTypeId
2
Score
0
ViewCount
0
LastEditorDisplayName
text
Body
<p>Using Visual Python I think I now have solved it:</p> <pre><code># Rotation first described for geo purposes: http://www.uwgb.edu/dutchs/mathalgo/sphere0.htm # http://stackoverflow.com/questions/6802577/python-rotation-of-3d-vector # http://vpython.org/ from visual import * from math import * import sys def ll2cart(lon,lat): # http://rbrundritt.wordpress.com/2008/10/14/conversion-between-spherical-and-cartesian-coordinates-systems/ x = cos(lat) * cos(lon) y = cos(lat) * sin(lon) z = sin(lat) return x,y,z def cart2ll(x,y,z): # http://rbrundritt.wordpress.com/2008/10/14/conversion-between-spherical-and-cartesian-coordinates-systems/ r = sqrt((x**2) + (y**2) + (z**2)) lat = asin(z/r) lon = atan2(y, x) return lon, lat def distance(lon1, lat1, lon2, lat2): # http://code.activestate.com/recipes/576779-calculating-distance-between-two-geographic-points/ # http://en.wikipedia.org/wiki/Haversine_formula dlat = lat2 - lat1 dlon = lon2 - lon1 q = sin(dlat/2)**2 + (cos(lat1) * cos(lat2) * (sin(dlon/2)**2)) return 2 * atan2(sqrt(q), sqrt(1-q)) if len(sys.argv) == 1: sys.exit() else: csv = sys.argv[1] # Points A and B defining the rotation: LonA = radians(float(sys.argv[2])) LatA = radians(float(sys.argv[3])) LonB = radians(float(sys.argv[4])) LatB = radians(float(sys.argv[5])) # A and B are both vectors # The crossproduct AxB is the rotation pole vector P: Ax, Ay, Az = ll2cart(LonA, LatA) Bx, By, Bz = ll2cart(LonB, LatB) A = vector(Ax,Ay,Az) B = vector(Bx,By,Bz) P = cross(A,B) Px,Py,Pz = P LonP, LatP = cart2ll(Px,Py,Pz) # The Rotation Angle in radians: # http://code.activestate.com/recipes/576779-calculating-distance-between-two-geographic-points/ # http://en.wikipedia.org/wiki/Haversine_formula RotAngle = distance(LonA,LatA,LonB,LatB) f = open(csv,"r") o = open(csv[:-4] + "_translated.csv","w") o.write(f.readline()) for line in f: (lon, lat) = line.strip().split(",") # Point C which will be translated: LonC = radians(float(lon)) LatC = radians(float(lat)) # Point C in Cartesian coordinates: Cx,Cy,Cz = ll2cart(LonC,LatC) C = vector(Cx,Cy,Cz) # C rotated to D: D = rotate(C,RotAngle,P) Dx,Dy,Dz = D LonD,LatD = cart2ll(Dx,Dy,Dz) o.write(str(degrees(LonD)) + "," + str(degrees(LatD)) + "\n") </code></pre>
Tags
Title
singulars
PostAcceptedAnswerId
1. This table or related slice is empty.
PostParentId
1. PORotate a sphere from coord1 to coord2, where will coord3 be?
  singulars
  PostTypePostTypeId
  PTQuestion
PostTypePostTypeId
1. PTAnswer
UserLastEditorUserId
1. This table or related slice is empty.
UserOwnerUserId
1. USjoosthoek
plurals
PostLinksPostIdRelatedPostId
1. This table or related slice is empty.
PostLinksRelatedPostIdPostId
1. This table or related slice is empty.
PostsAcceptedAnswerId
1. This table or related slice is empty.
PostsParentIdCreationDate
1. This table or related slice is empty.
VotesPostIdCreationDate
1. This table or related slice is empty.
CommentsPostId
1. This table or related slice is empty.

Querying!

Guidance

A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.