كيف نحول الاحداثيات الجغرافية من اهليلج wgs84 الى اهليلج clarcke1880

السلام عليكم
لدي برنامج يحول الاحداثيات الجغرافية المأخوذة من كوكل ايرث ذات الاهليلج العالمي wgs84 ,,الى احداثيات جغرافية ايضا لكن باهليلج كلارك 1880 ...
السؤال هو؟؟؟:
من المعلوم ان هنالك 7 بارامترات للاهليلج كلارك كما هي لكل انماط الاهليلج وهي
Dx - Dy - Dz - Rx - Ry - Rz - Scale Factor
ولدي في البرنامج خيار لتغير هذه البارامترات ,,لكن اثناء التحويل لاحظت ان جميع قيم هذه البارامترات هي صفر (0) ,,
هل ابقى على هذه القيم الصفرية لاحول الاحداثيات الجغرافية من اهليلج wgs84 الى احداثيات جغرافية باهليلج كلارك 1880 ,,,
ولاوضح اكثر تصور ان برنامج كوكل ايرث يستخدم اهليلج كلارك 1880 ,,,هل بارامترات هذه الاهليلج السبعة جميع قيمها هي صفر ..ام اننا نحتاج لبعض التعديل في هذه البارامترات ...
ولكم جزيل الشكر
واخص بالشكر الدكتور جمعه داوود

السلام عليكم

غير منطقي أن تكون عناصر التحويل السبعة = صفر !! فلا بد من وجود عناصر تحويل بين أي أهليجين عالميين مهما كانت قيم هذه العناصر صغيرة ، والا لكان كلاهما منطبق علي الاخر ويعدا أهليج واحد !

أنصحك باستخدام برنامج هيئة المساحة العسكرية الامريكية لتحويل الاحداثيات بين أكثر من مرجع وهو برنامج مجاني وسهل الاستخدام. أنظر الرابط:
http://www.gisclub.net/vb/showthread...ighlight=wgs84

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

السلام عليكم وجزاك الله كل الخير
الدكتور جمعه ..هذا البرنامج انف الذكر استخدمته سابقا ...وحتى انني طبقت مثال عن عين العبد في السعودية ..
المشكلة هنا في سوريا اريد ان احول الى نظام الاسقاط الستيريوغرافي المحلي ذات الداتم كلارك 1880 ...في اغلب البرامج تستطيع ان تجد داتم عين العبد السعودي ..لكن لا تجد الداتم السوري كلارك 1880 ...
فكيف لي ان احول الاحداثيات الجغرافية ذات الاهليلج wgs84 الى احداثيات جغرافية بداتم كلارك المحلي 1880 ...
وكيف لي ان اصنع داتم كلارك بواسطة البرنامج انف الذكر ...
بصراحة الموضوع معقد ورايت هنا في بلدي انه لا يوجد ولا برنامج تحويل ولا يوجد اي شيء على الاطلاق ...
فكيف لنا ان نجري التحويل ,,,,

بسم الله الرحمن الرحيم

تحية طيبة الى الدكتور جمعة داؤد المحترم
الاخ عبادة مالك المحترم
هنا يجب الاعتماد على البرمجة في عملية التحويل اذا لم تجد برنامج تحويل جاهز ،
ان الاسقاط الستيريو غرافي هو نوع من الاسقاطات واساسه الحفاظ على الزوايا عند التحويل من الشكل الكروي .
المهم :
يمكن اعتماد البرمجة العكسية للمعادلات في كلا النظامين (wgs84 و clark80 من قبل مبرمج محترف ، ومن ثم تطبيق مثال حقيقي يقاس من ال GPS في نقاط تثليث معروفة في سوريا ( يمكن استشارة هيئة الاستشعار عن بعد في دمشق ) ومقارنة التحويل مع القيم الحقيقية لنقاط التثليث

ارجو من الدكتور جمعة داؤد ابداء الرأي عن الفكرة


مع التقدير

أ. صباح حسين علي

[align=justify]السلام عليكم

نعم أخوتي الكرام الاساتذة عبادة و صباح: يمكننا استخدام أي لغة برمجه في اعداد برنامج للتحويل بين أي مرجعين أو datums لكن هذا يتطلب خبرة في البرمجة كما أشار الاستاذ صباح. يجب الاشارة الي أن معظم برامج نظم المعلومات الجغرافية و البرامج المساحية و برامج الجي بي اس - بصفة عامة - يكون داخلها امكانية انشاء نظام جغرافي جديد للمستخدم الذي لا يجد نظام معين يبحث عنه داخل ممكتبة نظم الاسقاط التي يضمها كل برنامج. وبالتالي فأن انشاء نظام اسقاط جديد لا يتطلب خبرة بالمعادلات الرياضية المستخدمه و لا خبرة في البرمجة ، كل ما علي هذا المستخدم أن يكون عارفا بالقيم أو العناصر المستخدمة في نوع الاسقاط المطلوب انشاؤه والتي سيقوم باستخدامها في نظامه الجديد ، وهذه القيم هي قيم الخمس عناصر التي تحدد أي نظام اسقاط وهم: خط الطول الرئيسي central meridian ودائرة العرض الرئيسية latitude of origin والاحداثيات الشرقية و الشمالية الزائفة false easting and northing بالاضافة لمعامل القياس scale factor. فاذا عرفت قيم هذه العناصر المستخدمة في خرائط دولة معينة - مثل سوريا مثلا - أستطيع انشاء نظام اسقاط جديد في حالة أني لا أجد هذا النظام داخل البرنامج.

وسنأخذ بالشرح كيفية انشاء نظام احداثيات جديد داخل برنامج Arc MAP ومن خلال 4 صور مبسطة كالاتي:


http://www.gisclub.net/vb/images/upl...aec293d386.GIF

http://www.gisclub.net/vb/images/upl...aec293f343.GIF

http://www.gisclub.net/vb/images/upl...aec3fae248.GIF

http://www.gisclub.net/vb/images/upl...aec3faea18.GIF

بهذه الطريقة أكون قد أنشأت نظام الاحداثيات السورية - مثلا - داخل برنامج Arc GIS وأستطيع بكل سهولة التعامل معه سواء بتحديده كأـساس لمشروع جديد أو طبقة جديدة ، أو لو عندي طبقة موجودة علي نظام احداثيات اخر مثل WGS84 أو UTM أستطيع تحويلها لهذا النظام الجديد من خلال استخدام أداة Project من صندوق الادوات Arc ToolBox .

هذا والله أعلي و أعلم.


[/align]

الاستاذ جمعه المحترم لدي سؤال حول هذا الموضوع
في المثال اعلاه تم عمل نظام جديد بحيث يمكن التعامل معه مستقبلا ولكن ما دور البارمترات السبعه التي تكلم عنها الاخ عباده ام يمكن اهمالها. وهل يمكن اعتماد هذا النضام في غير دول فمثلا نحن في العراق نستخدم الاهليج clarck 880 ايضا فهل يمكن استخدام النظام الذي تم عمله في المثال اعلاه للخرائط في بلدنا.

السلام عليكم أخي ياسر ورحمة الله و بركاته
بناء علي القيم التي سيتم ادخالها في الخطوة رقم 3 من الشرح السابق سيتم تحديد المنطقة التي يصلح لها هذا الاسقاط الجديد ، بمعني ادخال قيم معاملات نظام الخرائط السورية سينتج نظام لسوريا بينما ادخال القيم المستخدمة في الخرائط العراقية سينتج نظام للعرق ... وهكذا حتي و ان كانت كلا الدولتين تعتمدان علي نفس الالبسويد أو الاهليج: كلارك 1880.

ما سبق يتعلق فقط بنظم اسقاط الخرائط أما البارمترات السبعة فهي عناصر تحويل الاحداثيات من مرجع الي مرجع اخر وهي خطوة تتم قبل اسقاط الاحداثيات علي الخرائط.

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

السلام عليكم دكتور جمعه ورحمة الله وبركاته
شكرا جزيلا على ردكم ووفقكم الله لقد اصبحت الصوره واضحه الان فيما يتعلق بموضوع الاسقاط ولكن لا اخفي عليك يادكتور فان موضوع التحويل بين المراجع قد ارهقني وكلما حاولت التقدم في هذا الموضوع اعود الى نقطة الصفر لاني وبصراحه اريد فهم العلاقه بين مرجع جيودسي مركزه في مركز جذب الكره الارضيه او ما يسمى (Earth center Earth Fixed) واخر محلي (Local) لذا ساقوم بتفصيل الموضوع كالاتي.
1-نحن في العراق لدينا مرجع محلي قديم تم عمله من قبل البولونيين في السبعينيات على ما اعتقد ولدينا شبكه ضبط ارضي لهذا المرجع تسمى بالشبكه البولونيه. هذا المرجع يعتمد على الاهليج Clarck1880 وهو محلي كما ذكرت اي ان مركزه ليس متطابقا مع مركز الجذب وله نقطة اصل تقع في محافظة كربلاء علما ان هذه النقطه تمثل نقطة تقاطع سطح الاهليج مع الجيؤيد وليست مركز الاهليج هذا ما عرفته عن طريق السؤال.

2- حسب علمي عند استخدام اجهزة ال GPS فانها تعمل على نظام WGS84 العالمي علما ان نظامنا المحلي القديم غير معرف في اجهزة الGPS ولا في البرامج الملحقه بهذه الاجهزه. لذا نحن مضطرين الى استخدام نظام WGS العالمي.

3- السؤال هنا هل يمكن تحويل احداثيات النقاط المرصوده بنظام WGS84 الى نظامنا المحلي ولنسميه (كربلاء). علما اني قد سالت بعض الاصدقاء الذين سنحت لهم الفرصه بدوره تدريبيه في منظمة ال NGS وتمت مناقشة هذا الموضوع هناك فكانت اجابه خبراء منضمة ال NGS بانه لا يمكن التحويل بين المرجعين وان افضل طريقه طريقه للتحويل هي التحويل 2D حيث تقوم الجهه التي اعمل لديها الان برصد نقاط الشبكه البولونيه بنظام WGS العالمي وبهذا ستكون كل نقطه من نقاط الشبكه لديها مجموعتين من الاحداثيات واحده بالنظام البولوني والاخرى بنظام WGS84 ثم نقوم باستخدام ال Interpolation في عملية التحويل بين المرجعين فهل هذا الاسلوب صحيح علما اني قرات موضوع في الانترنت يتكلم عن نفس المشكله في الولايات المتحده وقد قامو بعمل برنامج يسمى Nadcon لاجراء نفس الغرض اي Interpolation بين النظامين. فهل من المعقول لايوجد علاقه هندسيه بين نظامين هندسيين غير متحدين المركز وفيهما اختلاف في التوجيه. شكرا واتمنى لكم التوفيق.
هذا الرابط يتكلم عن برنامج ال NadCon
http://www.ngs.noaa.gov/TOOLS/Nadcon/Nadcon.html

الدكتور جمعه اطال الله في عمرك وجزاك كل الخير ..
انت قلت ما يلي :البارمترات السبعة فهي عناصر تحويل الاحداثيات من مرجع الي مرجع اخر وهي خطوة تتم قبل اسقاط الاحداثيات علي الخرائط.
طيب على فرض انني حصلت على بارامترات الاسقاط اي (اي الستيريوغرافي السوري المحلي )
false easting and northing وcentral meridian وlatitude of origin و...الخ ..وكتبتها كما انت شرحته في هذه الصفحة ..
الان وداخل برنامج arcmap استطيع ان اغير بين انظمة الاسقاط فمثلا اختار نظام الاسقاط utm ..هنا بالتاكيد سوف تكون الاحدثيات دقيقة اليس كذلك ..اي التنقل من نظام الاسقاط الستيريوغرافي السوري الى اي نظام اخر فانه سوف يكون بشكل نظامي ودقيق
..طيب هل تؤثر بارامترات الداتم السبعة على عملية التحويل انفه الذكر ..(اي تحويل الاسقاط من الستيريوغرافي السوري الى utm ضمن واجه برنامج arcmap )
وانت تعلم ان البارامترات في البلد والدولة الواحدة تختلف من مكان لاخر ..فمثلا البارامترات جنوب سوريا غيرها وسط وشمال سوريا اليس كذلك ..
وللتوضيح هل اذا انا حصلت على بارامترات الاسقاط الخمسة ( ك central meridian ووو ) للاسقاط السيريوغرافي السوري .استطيع عندها التحويل بين جميع الاسقاطات لتعطي احداثيات دقيقة ..
وما هو تاثير البارامترات السبعة على هذه العملية ..
وما هي الحاله التي تستعدعي معرفة البارامترات السعبة واين نستخدمها

[align=justify]سأحاول توضيح الموضوع بأبسط عبارات وان كانت هذه العملية ليست سهله لأننا نتحدث عن موضوع علمي دقيق جدا ويحتاج لخبرة - متخصصة في المساحة الجيوديسية – لفهمه ويحتاج شرح معادلات و الاستعانة برسوم توضيحية ! وأدعو الله أن أصيب:

1- أولي الخطوات في العمل المساحي في منطقة كبيرة (أي فرع المساحة الجيوديسية التي لا تهمل تأثير تفلطح الأرض) أن نقوم بقياس إحداثيات المواقع في الطبيعة وبالتالي نحصل علي ثلاث قيم (خط طول ، دائرة عرض ، ارتفاع) وتكون هذه الإحداثيات منسوبة لاهليج أو البسويد معين. فمثلا لو استخدمنا تقنية الجي بي إس فأن الإحداثيات الناتجة ستكون منسوبة للاهليج العالمي WGS84.

2- الخطوة الثانية أننا نريد أن نوقع هذه الإحداثيات - للنقاط المرصودة – علي خريطة. لكن الخريطة ما هي إلا مستوي plane ثنائي الأبعاد وليست مجسم ثلاثي الأبعاد مثل الاهليج. الآن سنحتاج لمعادلات رياضية – تسمي معادلات إسقاط الخريطة – ليمكن تحويل الإحداثيات من المجسم أو الأرض الحقيقية إلي المستوي أو الخريطة. ولان هناك نظم عديدة من نظم الإسقاط فيجب الاعتماد علي النظام المستخدم في منطقة الدراسة (أو الدولة) ولا بد من معرفة عناصر هذا النظام وهم 5 عناصر كما سبق الإشارة إليهم وشرحهم.

3- إذا كنا رصدنا الإحداثيات في الموقع بالجي بي إس أي علي WGS84 وسنسقط الخريطة علي نظام يعتمد أيضا علي WGS84 فأننا لن نحتاج إلا لمعرفة عناصر الإسقاط. مثال: الإحداثيات بالجي بي إس والإسقاط سيكون بنظام UTM ويعتمد علي WGS84 أيضا ، فنحتاج فقط لمعرفة قيم العناصر الخمسة لنظام UTM.

4- النقطة السابقة هي الإجابة علي سؤال الأخ عبادة. فإذا عرفت عناصر الإسقاط الستيروغرافي السوري (وهو يعتمد علي أهليج كلارك 1880) فيمكنك تحويل الإحداثيات إلي أي نظام آخر لكن يجب أن يكون أيضا علي كلارك 1880 ، أي يمكنك التحويل من الاستيرجرافي السوري إلي UTM المبني علي كلارك 1880 أيضا.

5- لكن تأتينا المشكلة الثانية وهي: أننا نريد التحويل إلي نظام إحداثيات مبني علي أهليج آخر مختلف عن الاهليج الذي تم القياس عليه (أي WGS84 في حالة استخدامنا لأجهزة الجي بي إس).

6- في هذه الحالة يجب أن نقوم بخطوتين متتاليتين: أولا: أن نحول الإحداثيات ثلاثية الأبعاد (أي كما هي مقاسه في الطبيعة) من أهليج WGS84 إلي الاهليج المحلي (مثلا كلارك 1880) ، ثم: ثانيا: أن نسقط الإحداثيات التي تم تحويلها إلي كلارك 1880 من نوعية الإحداثيات ثلاثية الأبعاد إلي الإحداثيات المستوية لتوقيعها علي الخريطة (الخريطة تعتمد علي كلارك 1880).7-
في الخطوة الأولي نحن نقوم بما يعرف باسم التحويل بين المراجع Datum Transformation ولا بد من معرفة قيم العناصر السبعة التي تحدد العلاقة بين كلا المجسمين أو المرجعين أو الاهليجين Transformation Parameters ....... بينما في الخطوة الثانية نقوم بعملية إسقاط الخرائط Map Projection ونحتاج لمعرفة قيم العناصر الخمسة لنظام الإسقاط المطلوب التعامل معه لهذه الدولة.

8- مرة أخري سأعيد صياغة النقاط السابقة: إذا كنا نرصد نقاط علي أهليج وسنرسم الخريطة علي نفس الاهليج فلا نحتاج إلا لعناصر الإسقاط الخمسة .... وإذا كنا نرصد علي أهليج وسنرسم الخريطة علي أهليج مختلف فنحتاج لكلا من عناصر التحويل السبعة (لنقل الإحداثيات من أهليج لآخر ) وأيضا نحتاج لعناصر الإسقاط الخمسة (لإسقاط النقاط من مجسم الاهليج الثاني إلي الخريطة المستوية).

9- أعتقد أن النقطة الأخيرة تجيب علي تساؤل الأخ عبادة عن ما هي الحالة التي تستدعي معرفة البارمترات السبعة. إذا كنت سترصد بالجي بي إس علي WGS84 والخريطة التي ستستخدمها (مهما كان نظام إسقاطها: استيرجرافي أو UTM أو أي نظام آخر) لكنها تعتمد علي WGS84 أيضا فأنت لا تحتاج للبارمترات السبعة لأنك بكل بساطة لن تقوم بتغيير المرجع. وهذا الاتجاه بدأت به بعض الدول في كل أنحاء العالم – ومنهم مصر والسعودية – بأن تبدأ في إنشاء خرائط جديدة تعتمد علي أهليج أو مرجع WGS84 بدلا من المراجع المحلية القديمة وبالتالي فلن توجد حاجة لمعرفة عناصر تحويل المراجع (البارمترات السبعة) والمطلوب فقط سيكون عناصر الإسقاط الخمسة.

10- الآن نأتي لتساؤلات الأخ ياسر: حتي وان كنت أنظر إلي أهليجين عالميين وكلاهما يتخذ مركز الأرض مركزا له أي من نوع earth-fixed فتوجد علاقة رياضية تحدد أو تصف وضع أهليج بالنسبة للآخر. فإذا كان كلاهما earth-fixed فهذا معناه أن مركزهما واحد وبالتالي فأن عناصر الانتقال Translation Parameters (dX , dY , and dZ) يسووا صفر ، لكن هل محاور كل أهليج توازي محاور الاهليج الآخر ؟ من الممكن أن تكون الإجابة هي لا وبالتالي فلا بد من وجود 3 عناصر دوران Rotation Parameters بين المحاور ، أيضا هل حجم الاهليج الأول يساوي بالضبط حجم الاهليج الثاني؟ فأن كانت الإجابة هي لا فلا بد من معرفة قيمة معامل القياس Scale factor بين الاهليجين.

11- إذن لا بد من معرفة عناصر تحويل المراجع أو البارامترات السبعة Datum Transformation Parameters عند تحويل الإحداثيات من مرجع أو أهليج إلي آخر. وكمثال: من المعروف أن العلاقة بين المرجع العالمي WGS72 (الذي كان مستخدما في تقنية الدوبلر للرصد علي الأقمار الصناعية في السبعينات من القرن الميلادي الماضي) والمرجع العالمي WGS84 – ومع أن كلاهما earth-fixed – إلا أن عنصر الانتقال dZ=4.5 meter كما أتذكر. أي أن مرجع WGS84 يوازي مرجع WGS72 (3 عناصر الدوران = صفر) وأيضا له نفس الحجم (معامل القياس = صفر) ومركزهما واحد في المستوي الأفقي (dX = dY = 0) إلا أن مركز WGS84 يرتفع 4.5 متر عن مركز WGS72 . أنظر الشكل التالي:
http://www.gisclub.net/vb/images/upl...86fc1e4079.GIF

أو أحصل علي الشكل من المكتبة الرقمية المساحية في:
http://cid-0259cb4f889eaeb3.skydrive..._Transform.GIF

12- أما كيفية التحويل بين مرجعين مختلفين (حتي لو كان أحدهما عالمي و الآخر محلي) فالمعادلات موجودة بالتفصيل وبالرسم أيضا في الشكل التالي وهي موجودة في أي كتاب عن المساحة الجيوديسية أو جيوديسيا الأقمار الصناعية:
http://www.gisclub.net/vb/images/upl...86fc1ebd6a.GIF

أو أحصل علي الشكل من المكتبة الرقمية المساحية في:
http://cid-0259cb4f889eaeb3.skydrive...nsform%202.GIF

13- أما عن قول الأستاذ ياسر عن أن مركز المرجع العراقي وله نقطة أصل تقع في محافظة كربلاء ... فهنا نحتاج لبعض الإضافة: فهذه النقطة هي نقطة أصل نظام الإسقاط وليست نقطة أصل المرجع نفسه ! ما كان يحدث في الماضي (في كل الدول وليست العراق فقط) أن الدولة عندما تبدأ في إنشاء مرجع وطني لها فأنها تختار نقطة أصل - غالبا هي النقطة الرئيسية في شبكة الثوابت الأرضية للدولة – وتبدأ في بعض الافتراضات لكيفية وضع ودوران الاهليج العالمي الذي ستعتمد عليه (كلارك 1880 في حالة العراق). فمثلا: لتقليل الفرق بين الاهليج العالمي والشكل الحقيقي للأرض (الجيويد) فأننا نفترض أن المسافة بينهما (ما يسمي حيود الجيويد Geoid Undulation) = صفر عند نقطة الأصل هذه. في الحقيقة فأن هذا الفرق له قيمة وليس صفر ! لكن للسهولة ولقلة الإمكانيات التقنية في ذلك الوقت فكانت أي دولة تفترض هذا الفرض عند نقطة الأصل لها... هنا تكون الدولة قد قامت بتغيير الوضع الأصلي للاهليج (كلارك 1880 في حالة العراق) وبالتالي فلم يعد هذا الاهليج هو نفس الاهليج العالمي لكنه أصبح – بعد هذا التغيير أو هذا الفرض – مرجعا وطنيا لهذه الدولة فقط (مرجع محلي للعراق في حالتنا).

14- لنحسب قيم عناصر التحويل (البارمترات السبعة) بين أي مرجعين يجب أن نعرف إحداثيات بعض النقاط المشتركة بينهما . أي يجب أن يكون لدينا مجموعة (علي الأقل 3) من النقاط المعلوم إحداثياتها في المرجع الوطني ونقوم برصدها بأجهزة الجي بي إس لنحصل علي إحداثياتها علي المرجع العالمي WGS84. فإذا حصلنا علي هاتين المجموعتين من الإحداثيات لهذه النقاط المشتركة فنستخدم نفس المعادلة السابقة لنحسب منها قيم عناصر التحويل السبعة dX, dY, dZ, m , Ex, Ey, Ez .

15- طريقة حساب عناصر التحويل موجودة في أي برنامج حسابات GPS يأتي مع الأجهزة الهندسية (وليس الأجهزة الملاحية أو المحمولة يدويا) مثل برامج Trimble Geomatics, Lieca Office وهكذا ، فليس هناك حاجة لإنشاء برنامج يحل هذه المعادلات طالما أنك تمتلك أجهزة GPS ومعها برنامج حسابات Processing Software.

16- نأتي للجزء التالي من رد الأخ ياسر: طريقة عناصر التحويل السابق ذكرها هي العلاقة النظرية بين أي مرجعين وكانت تقريبا هي الطريقة الوحيدة في المساحة الجيوديسية قديما... لكن بدأ علماء المساحة في دراسة هذه الطريقة تفصيلا وبيان أن بها بعض العيوب ومن أهمها أن هذه الطريقة تساوي في الدقة بين إحداثيات كلا المرجعين وهذا يخالف الواقع... في الحقيقية أن إحداثيات الجي بي إس تكون دقيقة جدا (سنتيمترات أو أحيانا ملليمترات) بينما إحداثيات المرجع الوطني (المقام من عشرات السنين) لم تكن تصل إلي هذا المستوي من الدقة بسبب الأجهزة والتقنيات المستخدمة في ذلك الزمن البعيد (غالبا تكون دقة إحداثيات المرجع الوطني لأي دولة في حدود 0.5 متر أو أقل قليلا !!!). لذلك عندما أستخدم هذه الإحداثيات فأن دقة عناصر التحويل التي سأقوم بحسابها لن تكون أحسن من دقة الإحداثيات الوطنية ذاتها! أي في حدود 0.5 م وبالتالي فأننا نفقد دقة تقنية GPS عندما نحولها إلي المرجع الوطني.

17- بدأ العلماء في ابتكار وسائل جديدة لتحويل الإحداثيات بين مرجعين (بعيدا عن الطريقة التقليدية ذات العناصر السبعة) وظهرت عدة طرق جديدة منهم – علي سبيل المثال – Robust Estimation method and Surface correction method وطريقة سطح التحويل Surface correction هي ما تحدث عنه الأخ ياسر من أنها الطريقة المستخدمة في هيئة المساحة الأمريكية NGS .

18- تعتمد فكرة هذه الطريقة – ببساطة – في أننا نحسب فروق الإحداثيات (الفرق في خط الطول و الفرق في خط العرض) لمجموعة من النقاط المشتركة المعلوم لها إحداثياتها في المرجع الوطني و إحداثياتها في المرجع العالمي WGS84 . وفي الخطوة الثانية نقوم بتمثيل هذه الفروق بمعادلة انحدار regression model كدالة في الموقع. أي أننا نستنتج معادلة رياضية يمكنها حساب قيمة الفرق عند أي نقطة في منطقة الدراسة من خلال الاستنباط أو interpolation. فإذا طورنا معادلتين (أحدهما للفرق في خط الطول و الأخرى للفرق في دائرة العرض) يمكننا تمثيلهم علي سطحين تسميهما سطوح التحويل corrector surfaces ومن هنا جاء اسم الطريقة.

19- من أهم مميزات هذه الطريقة أنها سهلة التنفيذ (غالبا تكون معادلة درجة ثانية) وسهلة الحساب بأي برنامج للاستنباط interpolation. وأيضا أدق من الطريقة التقليدية للعناصر السبعة.

20- هذه الطريقة ليست جديدة تماما والحمد لله فقد قمت بعمل بحث عن تطبيقها في مصر منذ 9 سنوات (تم نشره في مؤتمر في عام 2000م) وقمت برفعه في المكتبة الرقمية المساحية في الرابط:
http://cid-0259cb4f889eaeb3.skydrive...rm%7C_2000.pdf
علما بأن فكرة هذا البحث قد تم استخدامها كموضوع رسالة ماجستير فيما بعد في مصر!

21- ردا علي تساؤل الأستاذ ياسر في نهاية رسالته: الطريقة الجديدة لا تغني عن الطريقة التقليدية لحساب عناصر التحويل بين أي مرجعين أحدهما عالمي earth-fixed والآخر محلي local datum مهما كان الاختلاف بينهما في المركز و التوجيه و الحجم. [color=#FF0000]فما زالت المعادلة الهندسية في الشكل السابق مطبقة بل هي الأصل – علميا و هندسيا – في وصف العلاقة بين أي مرجعين.[/color
]
22- أما عن قولك "أننا في العرق مضطرون لاستخدام نظام WGS84 العالمي" فليس هناك هذا الاضطرار ! إذا علمت إحداثيات مجموعة من النقاط (لا تقل عن 3 والأفضل الأكثر) في كلا المرجعين الوطني العراقي و المرجع العالمي فيمكنك حساب قيم عناصر التحويل السبعة وبالتالي يمكنك – بعد ذلك - تحويل إحداثيات أي نقطة GPS إلي المرجع العراقي الوطني باستخدام هذه العناصر.

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

هذا والله أعلي و أعلم.

[/align]

[align=justify]لمعلومات تفصيلية عن موضوع التحويل بين المراجع Datum Transformations فتوجد العديد من البحوث العلمية المنشورة علي شبكة الانترنت والتي يمكن الحصول عليها.

فعلي سبيل المثال يوجد في بحوث المؤتمر السنوي للاتحاد العالمي للمساحة FIG (والذي عقد في هونج كونج في الفترة 13-17 مايو 2007م) بحثين من مصر تناولا هذا الموضوع:

1- Remove-Restore Technique for Improving the Datum Transformation Process
لكلا من أ.د. أحمد شاكر ، أ.د. عبد الله سعد ، د. مني سعد ، د. عمرو حنفي

ويتناول فكرة تعديل الطريقة التقليدية للتحويل بين المراجع بغرض زيادة دقتها عم طريق أخذ أوزان مختلفة weights لكل نوع من الإحداثيات في الاعتبار قبل حساب قيم معاملات التحويل.
الرابط من موقع المؤتمر:
http://www.fig.net/pub/fig2007/paper..._etal_1218.pdf

2- Simple Model for Improving the Accuracy of the Egyptian Geodetic Triangulation Network
لكلا من أ.د. عبد الله سعد ، د. مني سعد.

ويتعلق بتطبيق فكرة الاستنباط interpolation لفروق الإحداثيات بين مرجعين.
الرابط من موقع المؤتمر:
http://www.fig.net/pub/fig2007/paper..._saad_1226.pdf

3- كذلك يوجد بحث آخر في نفس الموضوع منشور في مجلة البحوث المدنية CERM لكلية الهندسة بجامعة الأزهر المصرية في عدد أبريل 2005 م بعنوان:

Optimum Implementation of the Polynomials in Coordinate Transformation Process
لكلا من أ.د. أحمد شاكر ، أ.د. عبد الله سعد ، د. مني سعد

ويتعلق أيضا بتطبيق عدة طرق منهم الاستنباط interpolation لزيادة دقة عملية التحويل بين المرجعين العالمي WGS84 والمرجع الوطني المصري. وحيث أن هذه المجلة ورقية وليست رقمية فقد تم رفع هذا البحث للمكتبة الرقمية المساحية في:
http://cid-0259cb4f889eaeb3.skydrive...pt%7C_2005.pdf

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

[/align]

السلام عليكم ...بصراحة دكتور جمعه ..موضوع الاحداثيات هو امر معقد وليس بالامر السهل وكما انت قلت انه يحتاج الى خبراء وعلماء في عمليات التحويل ولكن لدي سؤال وباختصار ::
انت قلت لاخ ياسر مايلي
14- لنحسب قيم عناصر التحويل (البارمترات السبعة) بين أي مرجعين يجب أن نعرف إحداثيات بعض النقاط المشتركة بينهما . أي يجب أن يكون لدينا مجموعة (علي الأقل 3) من النقاط المعلوم إحداثياتها في المرجع الوطني ونقوم برصدها بأجهزة الجي بي إس لنحصل علي إحداثياتها علي المرجع العالمي WGS84. فإذا حصلنا علي هاتين المجموعتين من الإحداثيات لهذه النقاط المشتركة فنستخدم نفس المعادلة السابقة لنحسب منها قيم عناصر التحويل السبعة dX, dY, dZ, m , Ex, Ey, Ez .
نقطة انتهى
الدكتور جمعه :
الامور تخربطت علي قليلا ...
انا لدي مئات النقاط معلومة الاحداثيات بالاسقاط الستريوغرافي المحلي السوري ذو كلارك 1880 ..ومن خلال كوكل ايرث استطيع معرفة احداثيات هذه النقاط انفه الذكر وذلك باحداثيات جغرافية وutm ذو الاهليلج wgs84 ,,
هل انت قصدت هذا ..ام انه ينبغي علي ان احصل على احداثيات تلك النقاط (اي مئات النقاط التي امتلكها وهي ستيريوغرافي سوري ذو كلارك 1880 ) ان احصل على احداثياتها الجغرافية و utm على نفس الاهليلج كلارك 1880
..
انا لدي كل الخرائط على الاسقاط الستيريوغرافي المحلي ذو كلارك 1880 وعليها شبكة الاحداثيات الستيريوغرافي المحلية ذو كلارك 1880 ...وهي عبارة عن خرائط ورقية قديمة .وهي من نوع geotiff .وما اقوم به هو انني احفظ صور فضائية عالية الدقة لنفس تلك المناطق على شكل صورة jpeg لكي اطبقها على تلك الخرائط واستعمل لذلك برنامج global mapper وبعدها ادخل الصورة الى arcmap واجري لها georeference وبالفعل احصل على صورة فضائية تنطبق تماما وبدقة عالية جدا جدا على المخطط القديم المحلي ستيريوغرافي ذو كلارك 1880 ..علما ان جميع هذه المخططات هي مرجعه وتسقط على برنامج الاوتوكاد او arcmap باحداثياتها الحقيقة(ستيريوغرافي محلي ذو كلارك 1880 ) لانها من نوع geotiff ...وهذا ما اقوم به باختصار ,
الان دكتور جمعه على فرض انني حصلت على البارامترات السبعة الخاصة بالمرجع ...وذلك لداتم كلارك 1880 بالنسبة لااهليلج wgs84 فماذا علي ان افعل بعد ذلك ...وانا بالاصل لا املك ايضا البارامترات الخمسة ك central merdian والبقية ,الدكتور جمعه النقاش الذي انت شرحته كان في غاية الاهمية وكل ما دخلت الموقع كنت اتفاجا كيف اننا نأخذ منك الكثير من وقتك الثمين ....لقد شرحت لك كامل القصة التي اعاني منها ..ماذا يتوجب علي فعله والى اين اذهب
والله نعتذر منك كثيرا وجزاك الله كل الخير ....
وان شاء الله سوف لن يضيع الله لك هذا الجهد الذي تبذله ا لرفعه امة محمد صلى الله عليه وسلم

اقتباس:

---

المشاركة الأصلية كتبت بواسطة عبادة مالك

انا لدي مئات النقاط معلومة الاحداثيات بالاسقاط الستريوغرافي المحلي السوري ذو كلارك 1880 ..ومن خلال كوكل ايرث استطيع معرفة احداثيات هذه النقاط انفه الذكر وذلك باحداثيات جغرافية وutm ذو الاهليلج wgs84 ,,
هل انت قصدت هذا ..ام انه ينبغي علي ان احصل على احداثيات تلك النقاط (اي مئات النقاط التي امتلكها وهي ستيريوغرافي سوري ذو كلارك 1880 ) ان احصل على احداثياتها الجغرافية و utm على نفس الاهليلج كلارك 1880

انا لدي كل الخرائط على الاسقاط الستيريوغرافي المحلي ذو كلارك 1880 وعليها شبكة الاحداثيات الستيريوغرافي المحلية ذو كلارك 1880 ...وهي عبارة عن خرائط ورقية قديمة .وهي من نوع geotiff .وما اقوم به هو انني احفظ صور فضائية عالية الدقة لنفس تلك المناطق على شكل صورة jpeg لكي اطبقها على تلك الخرائط واستعمل لذلك برنامج global mapper وبعدها ادخل الصورة الى arcmap واجري لها georeference وبالفعل احصل على صورة فضائية تنطبق تماما وبدقة عالية جدا جدا على المخطط القديم المحلي ستيريوغرافي ذو كلارك 1880 ..علما ان جميع هذه المخططات هي مرجعه وتسقط على برنامج الاوتوكاد او arcmap باحداثياتها الحقيقة(ستيريوغرافي محلي ذو كلارك 1880 ) لانها من نوع geotiff ...وهذا ما اقوم به باختصار ,

الان دكتور جمعه على فرض انني حصلت على البارامترات السبعة الخاصة بالمرجع ...وذلك لداتم كلارك 1880 بالنسبة لااهليلج wgs84 فماذا علي ان افعل بعد ذلك ...وانا بالاصل لا املك ايضا البارامترات الخمسة ك central merdian والبقية ,الدكتور جمعه النقاش الذي انت شرحته كان في غاية الاهمية وكل ما دخلت الموقع كنت اتفاجا كيف اننا نأخذ منك الكثير من وقتك الثمين ....لقد شرحت لك كامل القصة التي اعاني منها ..ماذا يتوجب علي فعله والى اين اذهب

---

[align=justify]1- لحساب معاملات التحويل بين كلارك 1880 و المرجع العالمي WGS84 يحب معرفة احداثيات مجموعة من النقاط علي كلا المرجعين ، أي يجب أن تعرف احداثيات هذه النقاط (الاحداثيات الجغرافية: خط الطول ودائرة العرض و الارتفاع الجيوديسي) علي كلارك 1880 وأيضا احداثياتهم الجغرافية علي WGS84 . فان توفرت هذه البيانات يمكنك حساب عناصر التحويل - البارمترات - السبعة.

2- نعم يمكن الحصول من جوجل ايرث علي احداثيات النقاط علي مرجع WGS84 ... لكن بأي دقة ستكون هذه الاحداثيات؟ ... المعروف أن صور و احداثيات جوجل ايرث دقتها في حدود 10 أمتار .... وبالتالي فأنك ان استخدمت احداثيات WGS84 من جوجل ايرث فأن عناصر التحويل التي ستحسبها ستكون أيضا بدقة 10 أمتار !

3- جوجل ايرث حل سريع و مجاني ... لكن استخدامه يعتمد علي الدقة المطلوبة للعمل ... لا بد أن أسأل نفسي أولا: بأي دقة أريد أن أنجز عملي؟ هل يكفيني دقة 10 أمتار أم أنها غير مناسبة للعمل؟ فان كانت هذه الدقة مناسبة للعمل فالاسهل و الاسرع أن أحصل علي الاحداثيات من جوجل ايرث .... أما ان كان المشروع المطلوب يتطلب دقة أحسن (مشروع هندسي) فلن يصلح جوجل ايرث وسيكون البديل هو رصد احداثيات النقاط باستخدام أجهزة GPS هندسية (وليست محمولة يدويا hand-held receivers) حتي يمكننا حساب احداثيات WGS84 بدقة عالية وبالتالي عندما أحسب معاملات التحويل بين المرجعين سيكونوا بدقة عالية أيضا.

4- ونفس الوضع بالنسبة لاتمام عملية الارجاع الجغرافي باستخدام احداثيات من جوجل ايرث ... نفس السؤال سيتكرر مرة أخري: ما هي الدقة المطلوبة لانجاز هذا العمل أو هذا المشروع؟ .. هل 10 أمتار مناسبة أم أريد دقة أفضل؟ ... فان كانت الاجابة = نعم فما تقوم به هو الحل المطلوب ... وان كانت الاجابة = لا فلا بد من حساب معاملات تحويل دقيقية بين المرجعين ثم أستخدم أمر Project من Arc ToolBox لاتمام تغيير مرجع الطبقة من كلارك 1880 الي WGS84 أو العكس.

[/align]

السلام عليكم ...جزاك الله كل الخير الدكتور جمعه
لقد بدات الامور اكثر متعه وفهما وعسى ان تعم الفائدة لكل من سيدخل صفحة هذا النقاش المفيد
.الدكتور جمعه :
1- في برنامج arcinfo 9,2 h او اي برنامج gis هنالك عدة انواع من المساقط وايضا عدة انماط من الداتيوم .وسوف اعرض مثالا.
عند بدا تشكيل قاعدة البيانات من personal geodatabase ومن ثم اشكل طقم البيانات dataset هنا اقوم بتعريف نظام الاسقاط وليكن utm ذو الاهليلج wgs84 ثم بعدها اقوم بتشكيل feature classes ..ومن ثم ادخل واجه البرنامج arcmap 9,2 لاضيف البيانات من طبقة الشوراع مثلا والابنية وغيره ....
2- الان اريد مثلا عرض العمل بغير نظام الاسقاط هذا .فاقوم بعمل stop editing وبالزر اليمين للفارة على dataframe ومن ثم properties واقوم بتغير نظام الاسقاط وليكن مثلا utm ذو اهليلج هايفورد 1924 او مثلا جغرافي ذو اهليلج كلارك 1880 او مثلا stereographic-north-pole ذو اهليلج d-sphere-emep ,,
هنا اعتقد ان البرنامج يقوم بالتحويل بشكل نظامي لان جميع هذه الاسقاطات موجودة ضمن البرنامج اثناء التنصيب وان البرنامج يقوم بجميع الحسابات الرياضية للتحويل بين تلك الاسقاطات وبين داتيوم واخر ...اليس كذلك ام ان توقعي خاطا ,
3-وللتوضيح بفرض انني استخدم نظام الاسقاط utm ذو اهليلج wgs84 واريد ان ارى عملي ضمن واجه arcmap , وذلك بتغير نظام الاسقاط من utm ذو اهليلج wgs84 الى utm ذو اهليلج كلارك 1880 او اهليلج هايفورد 1924 هل هذا صحيح
,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
سؤال اخر .
بفرض انني استخدم نظام الاسقاط utm --wgs84 ولدي طبقات للشوارع والابنية ...الخ ..وايضا لدي صورة مخزنة على شكل raster dataset بنفس نظام الاسقاط utm-wgs84 واريد ان اصدر العمل على شكل shapefile ولكن بغير نظام الاسقاط هذا بل اريده مثلا جغرافي ذو اهليلج wgs84 هل يجب علي ان اطبق الطريقة التي انت شرحتها وهي من
أمر Project من Arc ToolBox
ام ان الطريقة التي انا شرحتها (وهي من خلال واجه arcmap وبالزر اليمين للفارة على data-frame ومن ثم تغير نوع الاسقاط ..) هي صحيحة لكي اصدر shape file بنظام اسقاط جغرافي ذو wgs84
- وهنالك طريقة اخرى وهي من واجه arccatalog استطيع تعديل نظام الاسقاط ل dataset الى نظام اخر ..ام انه حصرا كما انت شرحت وذلك من Project من Arc ToolBox
- وهل ذلك ينطبق على raster dataset اي الصورة حيث اننا قبل ادخالها الى واجه arcmap نقوم بتعرف نظام الاسقاط لها ضمن واجه arccatalog هل ايضا لتغير نظام الاسقاط لهذه الصورة يجب ان يتم عبر Project من arc toolbox
ام نستطيع تغير الاسقاط من واجه arccatalog وذلك من خلال خصائص الصورة ومن ثم تعديل نظام الاسقاط ..
,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, ,,,,,,,,,,
وجزاك الله كل الخير وان شاء الله ينفعنا الله واياكم من هذا العلم لرفعه امتنا

[align=justify]السلام عليكم

صراحة: لم أجرب الطريقة التي تتحدث عنها ولا أعرف هل هي صحيحة أم لا ! وهل مجرد تغيير نظام الاسقاط داخل خصائص الطبقة سيغير الطبقة فعلا للنظام الجديد؟ ... ومن قال لا أعلم فقد أفتي.

أما ما أعرفه فعلا وأعتقد أنه الحل الصحيح فهو: عند انشاء الطبقة في Arc Catalouge أحدد نظام الاسقاط و المرجع أو الداتم المطلوب لها ، وان لم أفعل ذلك أثناء انشاء الطبقة فيمكنني تنفيذه - لاحقا - باستخدام أمر Define Projection من صندوق الادوات Arc ToolBox ، فاذا أردت بعد ذلك تحويل هذه الطبقة من النظام الاول الي أي نظام اسقاط أو أي مرجع اخر فأقوم باستخدام أمر Project من صندوق الادوات.

نفس التسلسل للبيانات raster حيث يوجد في صنوق الادوات أيضا أمر Project Raster الخاص بتغيير نظام الاسقاط أو المرجع لهذا النوع من البيانات أو الصور. أما أمر Define Projection فهو صالح لكلا النوعين من البيانات vector data and raster data.

هذا والله أعلي و أعلم

[/align]

السلام عليكم وجزاك الله خير ا ...
الدكتور جمعه بعد هذه النقاشات الشاملة: اعتقد انه ومن خلال برنامج arcmap عند قيامنا بالتحويل بين مسقط واخر
على سبيل المثال utm ذو اهليلج كلارك 1880 الى utm ذو اهليلج هايفورد 1924 ..وحسب شرحك السابق فقدتم اعتماد اهليلج wgs84 كاساس ليتم حساب البارامترات السبعة لاي اهليلج او داتم بالنسبه له (اي بالنسبة ل wgs84 )
-الان اثناء التحويل من مسقط لاخر وحيث كلا المسقطين لهما اهليلج غير اهليلج wgs84 كما في المثال انف الذكر
(utm ذو اهليلج كلارك 1880 الى utm ذو اهليلج هايفورد 1924 ) فان برنامج arcmap سيقوم بالحساب بين هذه الداتم او الاهليلج كما في مثالنا ..اي بمعنى اوضح ان البرنامج يقوم بكافة العمليات الحسابية للبارامترات السبعة للاهليلج كلارك 1880 بالنسبة للاهليلج الرئيسي wgs84 ,,والبارامترات السبعة للاهليلج اهليلج هايفورد 1924 بالنسبة للاهليلج الرئيسي wgs84 ..وذلك للحصول على عملية التحويل بين هذه الانماط المختلفه للاهليلج ؟؟؟
- فقد لا تكون عباراتي دقيقة علميا لكنني كنت امثل العمليات التي يقوم بها برنامج arcmap ..وانت تعلم ان هذه الانماط من الاهليلج والداتم مخزنه في قاعدة بيانات برنامج arcmap ومعروفة البارامترات السبعة وبالتالي يمكنه التحويل ؟؟
..........................
هل كلامي هذا صحيح ام ان الامور تجري في غير هذا ؟؟؟؟؟؟؟؟؟
وجزاك الله كل الخير

سؤال اخر ::::؟؟
اذا كانت البارامترات السبعة في اي دولة هي غير ثابتة فهي مثلا تختلف من مكان لاخر فعلى سبيل المثال ..البارامترات السبعة في منطقة عسير في السعودية هي ليست ذاتها في منطقة تبوك شمال السعودية ...وحسب علمي فالاهليلج المستخدم في السعودية هو International 1924

Ain el Abd 1970 International 1924 -143 -236 7 Saudi Arabia 10 10 10 9
طالما ان البارامترات هي متغيرة من شمال السعودية الى جنوبها ..ولكن ارى انها ثابته وكيف تتعاملون بذلك من خلال برنامج arcmap فحسب علمي ان برنامج arcmap لديه فقط بارامترات السبعة ثابته لكل داتم او اهليلح ...
وطبعا هذا ينطبق على كل الدول
ما هو تعليق على ذلك

اقتباس:

---

المشاركة الأصلية كتبت بواسطة عبادة مالك

اثناء التحويل من مسقط لاخر وحيث كلا المسقطين لهما اهليلج غير اهليلج wgs84 (utm ذو اهليلج كلارك 1880 الى utm ذو اهليلج هايفورد 1924 ) فان برنامج arcmap سيقوم بالحساب بين هذه الداتم او الاهليلج كما في مثالنا ..اي بمعنى اوضح ان البرنامج يقوم بكافة العمليات الحسابية للبارامترات السبعة للاهليلج كلارك 1880 بالنسبة للاهليلج الرئيسي wgs84 ,,والبارامترات السبعة للاهليلج اهليلج هايفورد 1924 بالنسبة للاهليلج الرئيسي wgs84 ..وذلك للحصول على عملية التحويل بين هذه الانماط المختلفه للاهليلج ؟؟؟

---

[align=justify]لا .... برنامج Arc Map ليس برنامج مساحي في المقام الاول ، فهو لا يحسب البارمترات بين نماذج الاهليج المختلفة !! لكنه يحتوي بداخله علي قيم هذه البارمترات (لكل أنواع الاهليج علي مستوي العالم) ويقوم فقط بتطبيقها أو استخدامها عندما يريد المستخدم أن يحول طبقة من أهليج الي أهليج اخر. فقد قامت هيئة المساحة العسكرية الامريكية بمحاولة حساب البارمترات من WGS84 الي كل الانواع الاخري المستخدمة في معظم دول العالم (أكثر من 100 نظام أو مرجع محلي) ونشرت هذه القيم في تقرير لها صادر منذ عام 1988م ، ومعظم البرامج software تأخذ هذه القيم وتضعها داخل البرنامج لكي تطبقها ولا تضطر لحسابها من جديد.

مجرد ملحوظة يجب معرفتها في هذا السياق.

أما عن اختلاف قيم البارمترات من مكان لاخر (داخل نفس الدولة) فهذا مطبق فقط في الاعمال التي تتطلب دقة عالية (المشروعات الهندسية) أما بالنسبة للتطبيقات الجغرافية أو للتخطيط الاقليمي فأن الفروق ستكون بسيطة وغير مؤثرة علي دقة العمل المطلوب.

[/align]

السلام عليكم ...
اشكرك جزيل الشكر الدكتور جمعه داود ..وبالفعل هذا ما كنت اظنه ان البرنامج يختزن بداخله البارامترات ولا يقوم بالحساب ..
والفضل لله اولا وللدكتور جمعه ثانيا في فتح هذا النقاش وعسى ان يجعله الله في ميزان حسناتك ..

السلام عليكم ورحمته وبركاته
زملائي الاعزاء انا بحاجة ماسة الى معادلة التحويل من نضام احداثيات كلارك 1880 الى Wgs 84

مــــــــــــــــــــــــــــــع فائق الشكر والتقدير

المهندس صادق خليفة