longitude / latitude on sphere
-
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 27/10/2012 at 16:17, xxxxxxxx wrote:
Below is a piece of that I used to convert longitude / latitude to sphere coordinates.
However, it is not working and I cannot see what I do wrong.Input e.g. New York lat=40.714, lon=-74.606 and radius=100
longitude = c4d.utils.Rad(lat)
latitude = c4d.utils.Rad(lon)
x = radius * math.cos(latitude)*math.cos(longitude)
y = radius * math.cos(latitude)*math.sin(longitude)
z = radius * math.sin(latitude)
sphereobj = c4d.BaseObject(c4d.Osphere) # Create new sphere
sphereobj[c4d.PRIM_SPHERE_RAD] = 20 # set sphere label radius
sphereobj.SetRelPos(c4d.Vector(x,y,z)) # Set position of sphere
doc.InsertObject(sphereobj) -
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 28/10/2012 at 13:51, xxxxxxxx wrote:
Can you be more specific about what "not working" means?
I tried it like this and it seems to be working as expected:import c4d import math def main() : lat = 40.714 lon = -74.606 #New York radius=100 longitude = c4d.utils.Rad(lat) latitude = c4d.utils.Rad(lon) x = radius * math.cos(latitude)*math.cos(longitude) y = radius * math.cos(latitude)*math.sin(longitude) z = radius * math.sin(latitude) sphereobj = c4d.BaseObject(c4d.Osphere) # Create new sphere sphereobj[c4d.PRIM_SPHERE_RAD] = 20 # set sphere label radius sphereobj.SetName("NY") #Set the name of the sphere sphereobj.SetRelPos(c4d.Vector(x,y,z)) # Set position of sphere doc.InsertObject(sphereobj) c4d.EventAdd() if __name__=='__main__': main()
-ScottA
-
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 28/10/2012 at 15:02, xxxxxxxx wrote:
Hi pgroof,
here's a function I used in a plugin I made for a client.
def llarToWorld(lat, lon, alt, rad) : # see: http://www.mathworks.de/help/toolbox/aeroblks/llatoecefposition.html f = 0 # flattening ls = atan((1 - f)**2 * tan(lat)) # lambda x = rad * cos(ls) * cos(lon) + alt * cos(lat) * cos(lon) z = rad * cos(ls) * sin(lon) + alt * cos(lat) * sin(lon) y = rad * sin(ls) + alt * sin(lat) return c4d.Vector(x, y, z)
-Nik
-
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 29/10/2012 at 03:36, xxxxxxxx wrote:
ScottA, thanks, triggered by your remark I tried again.
Changing y = -y + set the globe tp -90,90,-90 helped.
It is ok now.Hi NiklasR, yes, what you are using is the flattened version of the algorithm.
Thanks for your help.
-
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 29/10/2012 at 03:37, xxxxxxxx wrote:
Sorry,
Changing y = -y + set the globe tp -90,90,-90 helped.should be
Changing y = -y + set the globe to -90,90,-90 helped. -
THE POST BELOW IS MORE THAN 5 YEARS OLD. RELATED SUPPORT INFORMATION MIGHT BE OUTDATED OR DEPRECATED
On 29/10/2012 at 03:39, xxxxxxxx wrote:
Or even better (I'm missing the edit function here) :
Changing y to -y (y=-y) and setting the rotation of the globe to -90,90,-90 solved the issue.