كيفية التحويل من الـutm الى اسقاط محلي؟؟؟

تحية طيبة وبعد

يتوفر لدي احداثيات لمجموعة من النقاط وفق اسقاط utm واسقاط اخر محلي

ما اود معرفته هو كيفية حساب البارمترات السبعة للتحويل من الـutm الى الاسقاط المحلي اي dx,dy,dz,rx,ry,rz,k

وما علمته اني بحاجة الى برنامج لحساب هذه البارمترات (كالتي تستخدم في حساب شبكات الـGPS)، لكني حتى الآن لم استطع الحصول على واحد منها

فهل اجد لديكم اي فكرة او مساعدة؟؟؟

محمد هادي جواد كتب:ما اود معرفته هو كيفية حساب البارمترات السبعة للتحويل من الـutm الى الاسقاط المحلي اي dx,dy,dz,rx,ry,rz,k ، وما علمته اني بحاجة الى برنامج لحساب هذه البارمترات (كالتي تستخدم في حساب شبكات الـGPS)، لكني حتى الآن لم استطع الحصول على واحد منها ، فهل اجد لديكم اي فكرة او مساعدة؟؟؟

السلام عليكم

نعم أنت تحتاج لبرنامج متخصص من برامج GPS Processing مثل برنامج LGO من شركة ليكا أو برنامج Geomatics من شركة ترمبل. سياسة الشركة التي تستضيف هذا المنتدي تمنع وضع نسخ من البرامج التجارية والا يتم غلق المنتدي ! لكننا سبق أن جاوبنا علي هذا السؤال في موضوع سابق ووضعنا روابط لمواقع أخري يمكنك منها تحميل نسخة من هذه البرامج ، أبحث في المنتدي حتي تجد هذا الموضوع وبه الروابط المطلوبة.

أطيب التمنيات بالتوفيق.

ماهي المعادلات المستخدمه ف التحويل

hamadayounis كتب:ماهي المعادلات المستخدمه ف التحويل

السلام عليكم

المعادلات المستخدمة في التحويل ليست سهلة كما يعتقد البعض ! وهي موجودة في أكثر من مرجع أو كتاب من كتب المساحة الجيوديسية (مثلا كتاب د. علي شكري عن المساحة الجيوديسية وهو من أهم المراجع المساحية في مصر). لكن هناك العديد من صفحات الانترنت التي تعرض هذه المعادلات أيضا مع شرح بسيط لها ، ومنهم مثلا صفحة:

http://www.uwgb.edu/dutchs/usefuldata/utmformulas.htm

وهذه الصفحة أيضا تعرض ملف اكسل لاتمام هذه التحويلات وتجده في الرابط:
http://www.uwgb.edu/dutchs/usefuldata/UTMConversions1.xls

ملحوظة: أنا لم أجرب هذا الملف وبالتالي لا أستطيع الحكم عليه !!

__________________________

المعادلات كالاتي:

Formulas For Converting Latitude and Longitude to UTM:

These formulas are slightly modified from Army (1973). They are accurate to within less than a meter within a given grid zone. The original formulas include a now obsolete term that can be handled more simply - it merely converts radians to seconds of arc. That term is omitted here but discussed below:

Symbol

slat = latitude of point
long = longitude of point
long0 = central meridian of zone
k0 = scale along long0 = 0.9996.

Even though it's a constant, we retain it as a separate symbol to keep the numerical coefficients simpler, also to allow for systems that might use a different Mercator projection.

e = SQRT(1-b2/a2) = .08 approximately.

This is the eccentricity of the earth's elliptical cross-section.

e'2 = (ea/b)2 = e2/(1-e2) = .007 approximately.

The quantity e' only occurs in even powers so it need only be calculated as e'2.

n = (a-b)/(a+b) rho = a(1-e2)/(1-e2sin2(lat))3/2.

This is the radius of curvature of the earth in the meridian plane.

nu = a/(1-e2sin2(lat))1/2.

This is the radius of curvature of the earth perpendicular to the meridian plane.

It is also the distance from the point in question to the polar axis, measured perpendicular to the earth's surface.

p = (long-long0) in radians

(This differs from the treatment in the Army reference)

Calculate the Meridional ArcS is the meridional arc through the point in question (the distance along the earth's surface from the equator). All angles are in radians.


S = A'lat - B'sin(2lat) + C'sin(4lat) - D'sin(6lat) + E'sin(8lat)

where lat is in radians and

A' = a[1 - n + (5/4)(n2 - n3) + (81/64)(n4 - n5) ...] and

B' = (3 tanS/2)[1 - n + (7/8)(n2 - n3) + (55/64)(n4 - n5) ...] and

C' = (15 tan2S/16)[1 - n + (3/4)(n2 - n3) ...] and

D' = (35 tan3S/48)[1 - n + (11/16)(n2 - n3) ...] and

E' = (315 tan4S/512)[1 - n ...] and

The USGS gives this form, which may be more appealing to some. (They use M where the Army uses S
)


M = a[(1 - e2/4 - 3e4/64 - 5e6/256 ....)lat
- (3e2/8 + 3e4/32 + 45e6/1024...)sin(2lat)
+ (15e4/256 + 45e6/1024 + ....)sin(4lat)
- (35e6/3072 + ....) sin(6lat) + ....)] where

lat is in radians This is the hard part.

Calculating the arc length of an ellipse involves functions called elliptic integrals, which don't reduce to neat closed formulas. So they have to be represented as series.

Converting Latitude and Longitude to UTMAll angles are in radians.

y = northing = K1 + K2p2 + K3p4, where

K1 = Sk0, and
K2 = k0 nu sin(lat)cos(lat)/2 = k0 nu sin(2 lat)/4 and
K3 = [k0 nu sin(lat)cos3(lat)/24][(5 - tan2(lat) + 9e'2cos2(lat) + 4e'4cos4(lat)] SO:

x = easting = K4p + K5p3, where

K4 = k0 nu cos(lat) and
K5 = (k0 nu cos3(lat)/6)[1 - tan2(lat) + e'2cos2(lat)]

Easting x is relative to the central meridian. For conventional UTM easting add 500,000 meters to x.

________________________________

بالتوفيق للجميع بمشيئة الله

الرابط التالي يعرض ملف pdf بتاريخ يوليه 2011 و حجمه 3 ميجا و مكون من 132 صفحة (باللغة الانجليزية) من مجموعة الملفات التدريبية للمنظمة الدولية لمنتجي الغاز و البترول وبه شرح كامل لجميع المعادلات المستخدمة في تحويل الاحداثيات علي اختلاف أنواعها و أنظمتها:

http://www.epsg.org/guides/docs/G7-2.pdf