1. مقدمه
اعوجاج در پیش بینی های نقشه اجتناب ناپذیر است [ 1 ، 2 ]، و همچنان یک موضوع چالش برانگیز در نقشه کشی [ 3 ] است. اعوجاج در پیش بینی های نقشه باید اندازه گیری شود [ 4 ، 5 ] برای ارزیابی [ 6 ]، مقایسه، انتخاب [ 7 ، 8 ]، بهبود، و بهینه سازی [ 9 ] پیش بینی های نقشه موجود، و همچنین برای طراحی پیش بینی های نقشه جدید [ 3 ، 9 ] , 10 , 11 ].
چندین معیار بصری [ 12 ، 13 ، 14 ] یا عددی [ 4 ، 15 ، 16 ، 17 ، 18 ، 19 ، 20 ] از اعوجاج برای پیش بینی نقشه پیشنهاد شده است. بیشتر معیارهای موجود به معادلات رو به جلو یا معکوس پیش بینی های نقشه یا وابسته به محاسبات دیفرانسیل معادلات رو به جلو یا معکوس هستند. در میان آنها، مزایای شاخص های مبتنی بر قوس دایره بزرگ کروی (GCA) [ 19]، که در مطالعه قبلی ما پیشنهاد شد، استقلال در محاسبات دیفرانسیل است. با این حال، شاخص های مبتنی بر GCA به معادلات معکوس پیش بینی نقشه ها وابسته هستند. به عنوان یک جایگزین، یک متریک مبتنی بر درون یابی بهبود یافته (متریک GCA رو به جلو یا متریک مبتنی بر FWD-GCA) [ 20 ] که مستقل از معادلات معکوس است در مطالعه قبلی ما معرفی شده است. رگرسیون مبتنی بر تابع منطقی در متریک مبتنی بر FWD-GCA [ 20 ] برای پیشبینی حداکثر اعوجاجهای زاویهای به کار گرفته شده است.
رویکردهای مبتنی بر تصویر [ 21 ، 22 ، 23 ، 24 ] به طور گسترده در حوزه های مختلف استفاده می شوند. پیش بینی های نقشه عموماً فرم های تصویری برای به تصویر کشیدن ویژگی های طبیعی یا اجتماعی هستند [ 25]. در این مطالعه، ما از ویژگی تصویر در پیشبینیهای نقشه و ویژگی استقلال تفاضلی شاخصهای مبتنی بر GCA برای ایجاد یک متریک مبتنی بر تصویر جدید و ارزیابی اعوجاجها در پیشبینیهای نقشه بهرهبرداری کردیم. تصاویر از پیش تعریف شده با الگوهای شناخته شده در رویکرد مبتنی بر تصویر ما برای نشان دادن رابطه بین رنگ ها و مختصات عرض جغرافیایی (یا طول جغرافیایی) استفاده شد. نرم افزار نقشه برداری در رویکرد ما برای تولید تصاویر تحریف شده به کار گرفته شد. استفاده از نرم افزار نقشه برداری از فرمول های صریح پیش بینی نقشه ها جلوگیری می کند. از آنجایی که خطاهای زیادی در رویکردهای مبتنی بر تصویر معمولی وجود دارد، کاهش نویز مبتنی بر اتصالات سطحی برای کاهش خطاها در معیارهای مبتنی بر تصویر معرفی شد. علاوه بر این، برخی از رویکردهای رایج فیلتر کردن تصویر (به عنوان مثال، فیلتر میانه و فیلتر گاوسی) [ 26] به معیار پیشنهادی برای بهبود بیشتر نتایج بصری یا عددی روش ما اعمال شد.
ادامه این مقاله به شرح زیر سازماندهی شده است. متریک مبتنی بر تصویر و روش کاهش نویز مبتنی بر اتصالات سطح همراه در بخش 2 معرفی شدهاند . آزمایش ها و اعتبارسنجی در بخش 3 ارائه شده است. در نهایت، بخش 4 نتایج و نتیجه گیری این مطالعه را ارائه می کند.
2. مواد و روشها
2.1. مروری بر روش ها
شکل 1 نمای کلی و مقایسه ای از معیارهای مبتنی بر دیفرانسیل سنتی، مبتنی بر GCA قبلی، مبتنی بر FWD-GCA قبلی، و معیارهای مبتنی بر تصویر پیشنهادی را ارائه می دهد.
معیارهای مبتنی بر دیفرانسیل (نگاه کنید به صفحه 24 در [ 15 ]) برای محاسبه حداکثر اعوجاج زاویه ای استفاده شد (به نماد مراجعه کنید در شکل 1 ) و اعوجاج منطقه با استفاده از مشتقات جزئی معادلات رو به جلو پیش بینی نقشه. معیارهای مبتنی بر GCA [ 19 ] (به نماد مراجعه کنید در شکل 1 ) معیارهای مستقل از دیفرانسیل هستند که برای اجتناب از محاسبات دیفرانسیل استفاده می شوند. آنها همبستگی بالایی با معیارهای مبتنی بر دیفرانسیل کلاسیک دارند [ 19 ]. اگرچه محاسبات دیفرانسیل در معیارهای مبتنی بر GCA مورد نیاز نیست، تبدیلهای معکوس (به خط آبی در شکل 1 مراجعه کنید ) پیشبینیهای نقشه مورد نیاز است. برای اجتناب از وابستگی تبدیلهای معکوس، معیارهای مبتنی بر FWD-GCA [ 20 ] (به نماد مراجعه کنید در شکل 1 ) پیشنهاد شد. این معیارها از رویکردی مبتنی بر درونیابی برداری استفاده می کنند (خط بنفش را در شکل 1 ببینید). معیارهای مبتنی بر FWD-GCA با معیارهای مبتنی بر دیفرانسیل همبستگی زیادی دارند [ 20 ].
معیارهای مبتنی بر GCA و FWD-GCA به صراحت به معادلات رو به جلو وابسته هستند (خطوط نارنجی در شکل 1)) پیش بینی های نقشه برای تبدیل نقاط نمونه. علاوه بر این، معیارهای مبتنی بر GCA به طور صریح به معادلات معکوس پیشبینیهای نقشه برای تبدیل نقاط نمونه وابسته هستند. در این مطالعه، ما یک متریک مبتنی بر تصویر را برای اندازهگیری اعوجاج پیشبینیهای نقشه پیشنهاد کردیم. متریک پیشنهادی یک روش بهبود معیارهای قبلی مبتنی بر GCA است، اما به طور صریح به معادلات رو به جلو یا معکوس پیشبینیهای نقشه وابسته نیست. روش پیشنهادی از تصاویر تحریف شده تولید شده توسط نرم افزار نقشه برداری به عنوان ورودی استفاده می کند و معیارهای اعوجاج را بر اساس تصاویر تحریف شده محاسبه می کند. برای بهبود دقت روش مبتنی بر تصویر، کاهش نویز مبتنی بر اتصال سطحی بر این اساس معرفی شد. چندین برنامه نرم افزار نقشه برداری وجود دارد [ 27]، مانند G.Projector ناسا، Geocart Mapthematics LLC و ArcGIS Esri. در این مطالعه، از نرم افزار NASA G.Projector 3 برای تولید تصاویر تقریباً 200 نقشه پیش بینی استفاده شد. این پیش بینی های نقشه با معیارهای پیشنهادی ارزیابی شدند. علاوه بر این، دقت معیارهای پیشنهادی ما مورد ارزیابی قرار گرفت.
2.2. متریک مبتنی بر کمان دایره بزرگ
معیارهای مبتنی بر GCA [ 19 ] اصل ریاضی اساسی رویکرد مبتنی بر تصویر پیشنهادی است.
برای طرح ریزی نقشه خاص، ما چهار نقطه واقع در گوشه های یک مربع بسیار کوچک در صفحه نمایش را انتخاب کردیم (نگاه کنید به در شکل 1 ). مساحت مربع بسیار کوچک در صفحه طرح ریزی در نظر گرفته شد ، جایی که یک مقدار به اندازه کافی کوچک است. همانطور که در شکل 1 نشان داده شده است ، چهار نقطه به یک کره واحد تبدیل می شوند با استفاده از معادلات معکوس (به خط آبی در شکل 1 مراجعه کنید ) پیش بینی نقشه. مختصات سه بعدی چهار نقطه روی برای محاسبه معیارهای مبتنی بر GCA استفاده شد.
نشان دادیم به عنوان مختصات دکارتی از چهار نقطه سه بعدی در ، . علاوه بر این، ما فرض کردیم که نقاط مربوطه (روی صفحه طرح ریزی ) از و در یک مورب مربع بسیار کوچک و نقاط مربوطه (روی صفحه طرح ریزی) قرار داشتند ) از و در مورب دیگر مربع قرار داشتند.
معیارهای مبتنی بر GCA به شرح زیر بدست آمد.
ابتدا بردارهای سه بعدی را محاسبه کردیم ، جایی که ، برقراری ارتباط ، جایی که و بردارهای سه بعدی نرمال شده ، جایی که . سپس طول قوس را محاسبه کردیم بر ، جایی که .
سپس، مولفه های طول و زاویه معیارهای مبتنی بر GCA محاسبه شد.
در نهایت، ما دو جزء بالا را به عنوان یک متریک اعوجاج شکل مبتنی بر GCA ترکیب کردیم با میانگین حسابی
اعوجاج زاویه ای پیش بینی نقشه از محاسبه شد با استفاده از یک تابع منطقی [ 20 ] با یک درجه برای پیش بینی نقشه های غیر استوانه ای، غیر مخروطی و غیرآزیموتالی.
جایی که ضریب با استفاده از رگرسیون مبتنی بر تابع منطقی محاسبه شد، . مقادیر نامزد برای ضریب در مطالعه قبلی به دست آمد. به این معنا که، و [ 20 ].
برای پیش بینی های نقشه استوانه ای، مخروطی یا ازیموتال، مطالعه قبلی [ 20 ] معادله ساده تری را برای محاسبه اعوجاج زاویه ای نشان داد. .
علاوه بر این، اعوجاج ناحیه مبتنی بر GCA با استفاده از آن به دست آمد
جایی که مساحت مربع بسیار کوچک در صفحه نمایش است .
از مطالعه قبلی [ 20 ] استنباط شد که تخمین خوبی از متریک زاویه ای مبتنی بر دیفرانسیل بود ، و تقریباً همان متریک متریک مساحت مبتنی بر دیفرانسیل بود.
2.3. تبدیل مختصات مبتنی بر تصویر
در متریک مبتنی بر تصویر پیشنهادی، از دو تصویر رنگی مستطیلی از پیش تعریف شده (نگاه کنید به شکل 2 ) به ترتیب برای تبدیل مختصات طول و عرض جغرافیایی استفاده کردیم. این تصاویر از پیش تعریف شده از طرح ریزی Plate Carrée برای ایجاد یک رابطه رو به جلو و معکوس شناخته شده بین رنگ ها و مختصات عرض جغرافیایی (یا طول جغرافیایی) استفاده می کنند. در شکل 2، هر پیکسل در یک ردیف در یک تصویر عرض جغرافیایی از پیش تعریف شده دارای رنگ یکسان است. به طور مشابه، هر پیکسل در همان ستون در یک تصویر طول جغرافیایی از پیش تعریف شده دارای رنگ یکسان است. پیکسل ها در ردیف های مختلف در تصویر عرض جغرافیایی و همچنین در ستون های مختلف در تصویر طول جغرافیایی دارای رنگ های متمایز هستند. بنابراین، ما رابطه شناخته شده بین رنگ ها و طول جغرافیایی (یا طول جغرافیایی) را در تصاویر از پیش تعریف شده برای به دست آوردن مختصات نقاط در تصاویر تحریف شده عرض جغرافیایی (یا طول جغرافیایی) فرض کردیم.
ما W و H را به ترتیب به عنوان عرض و ارتفاع تصاویر از پیش تعریف شده عرض جغرافیایی (یا طول جغرافیایی) نشان دادیم، جایی که . نمادهای x و y به ترتیب مختصات پیکسل افقی و عمودی را نشان میدهند ، و مبدا مختصات پیکسل در گوشه سمت چپ بالای تصاویر است. W و H برای اطمینان از دقت و کاهش خطاهای تبدیل مختصات به اندازه کافی بزرگ هستند.
تصاویر رنگی واقعی 24 بیتی از پیش تعریف شده دارای سه کانال (R، G و B) هستند. هر کانال از 8 بیت در هر پیکسل برای نمایش 256 سطح مختلف مقیاس خاکستری استفاده می کند. برای تصاویری که پارامترهای آنها H < W <256 × 256 = 65536، تنها 2 کانال (R و G در رویکرد ما) برای تصاویر از پیش تعریف شده کافی است تا رنگ های متمایز استفاده شده برای هر ردیف در تصویر عرض جغرافیایی و رنگ های متمایز استفاده شده برای هر ستون را تضمین کند. در تصویر طول جغرافیایی برای تصاویر با عرض بیشتر از 65536، هر 3 کانال در تصاویر از پیش تعریف شده مورد نیاز است.
ما فرض کردیم که از دو کانال در تصاویر از پیش تعریف شده در مطالعه ما استفاده شده است و نشان داده شده است به عنوان مختصات پیکسل خاص در تصویر از پیش تعریف شده عرض جغرافیایی (یا طول جغرافیایی)، به عنوان مقادیر رنگ کانال های R و G برای در تصویر عرض جغرافیایی از پیش تعریف شده، و ) به عنوان مقادیر رنگ کانال های R و G برای ، به ترتیب در تصویر طول جغرافیایی از پیش تعریف شده. با در نظر گرفتن پروژکشن Plate Carrée مورد استفاده در تصاویر از پیش تعریف شده، رابطه تبدیل از مختصات پیکسل در تصاویر از پیش تعریف شده به مقادیر رنگ و به عنوان ارائه می شوند
جایی که تابع کف است، عملیات مدول، P یک عدد صحیح و P < 256 است.
فرمول های بالا دو فرمول از پیش تعریف شده ( ) تصاویر با مقادیر رنگ متوالی اعداد صحیح و برای مختصات پیکسلی اعداد صحیح متوالی . تصاویر از پیش تعریف شده به عنوان تصاویر ورودی نرم افزار نقشه برداری برای ایجاد تصاویر خروجی تحریف شده با استفاده از پیش بینی های نقشه خاص با وضوح خروجی کافی برای اطمینان از صحت محاسبات بعدی استفاده شدند. نمونه ای از تصاویر خروجی تحریف شده برای طرح Bonne در شکل 3 نشان داده شده است .
سپس، تبدیلهای معکوس معادلات ( 10 )-( 13 ) را استخراج کردیم که رابطه تبدیل (به معادلات ( 14 ) و 15 را ببینید ) از مقادیر رنگ در تصاویر تحریف شده به مختصات طول و عرض جغرافیایی ایجاد کرد.
با توجه به و ، مقادیر رنگ کانال های R و G برای یک نقطه خاص در تصاویر خروجی تحریف شده طول و عرض جغرافیایی (نگاه کنید به شکل 3 ) ایجاد شده توسط نرم افزار نقشه برداری، مختصات طول و عرض جغرافیایی را می توان به دست آورد
جایی که و توسط محاسبه می شوند
2.4. کاهش نویز مبتنی بر اتصالات سطحی
با توجه به ، یعنی مختصات یک نقطه در تصاویر خروجی تحریف شده که توسط نرم افزار نقشه برداری تولید می شود، مختصات طول و عرض جغرافیایی مربوطه. از نقطه می توان به طور مستقیم با جایگزینی مقادیر رنگ محاسبه کرد و از نقطه به معادلات ( 14 )–( 17 ).
مشابه معیارهای مبتنی بر GCA، که از چهار نقطه واقع در گوشه های یک مربع بسیار کوچک در صفحه نمایش استفاده می کند. ، چهار نقطه را انتخاب کردیم که شرایط یکسان گوشه های مربع را دارند. اینها در صفحه تصاویر خروجی تحریف شده بودند، و ما از این چهار نقطه برای محاسبه یک متریک به معنای معیارهای مبتنی بر GCA استفاده کردیم.
با این حال، با توجه به گسست اعداد صحیح استفاده شده در ( 14 )-( 17 )، و وضوح نمایش محدود تصاویر ورودی از پیش تعریف شده و خروجی تحریف شده نرم افزار نقشه برداری، مختصات طول و عرض جغرافیایی محاسبه شده با تبدیل مختصات مبتنی بر تصویر ممکن است نسبتاً بزرگ باشد. خطاها این خطاها ممکن است برای معیارهای مبتنی بر تصویر از نظر معیارهای مبتنی بر GCA اعمال شوند.
در اینجا، ما یک رویکرد کاهش نویز مبتنی بر اتصالات سطح را برای کاهش خطاها و بهبود دقت محاسبه تبدیل مختصات مبتنی بر تصویر و معیارهای مبتنی بر تصویر به دست میآوریم.
توابع چند جمله ای دو متغیره برای برازش سطح در صفحه تصویر به کار گرفته شد ( شکل 3 را ببینید ). ما d را به عنوان نصف اندازه پنجره برای اتصالات سطح مشخص کردیم. در نظر گرفتن نکته و یک پنجره با اندازه ، داخل پنجره در مرکز نقطه قرار دارد در فرآیند برازش سطح استفاده شد.
دو تابع چند جمله ای دو متغیره برای ایجاد فرآیند برازش چند جمله ای برای تبدیل طول و عرض جغرافیایی به ترتیب مورد نیاز بود. ما فرآیند برازش چند جمله ای را با توابع چند جمله ای دو متغیره درجه دوم توصیف کردیم.
ما دو تابع چند جمله ای دو متغیره با دوازده ضریب مجهول را به شرح زیر فرض کردیم:
برای تبدیل طول و عرض جغرافیایی، به ترتیب.
پس از آن، مختصات طول و عرض جغرافیایی را محاسبه کردیم برای هر نقطه ، که در داخل پنجره در مرکز قرار داشت ، با استفاده از تبدیل مختصات مبتنی بر تصویر مطابق با معادلات ( 14 )–( 17 )، که در آن ، و .
بعد، مختصات دو سری را جایگزین کردیم نقاط، جایی که و ، وارد دو تابع چند جمله ای بالا شده و سپس دو سیستم معادلات خطی به دست می آید:
جایی که
ماتریس در معادله ( 22 ) مختصات سری نقاط داخل پنجره را نشان می دهد. سمت راست معادلات ( 20 ) و ( 21 ) مختصات طول و عرض جغرافیایی است که با تبدیل مختصات مبتنی بر تصویر آن سری از نقاط محاسبه می شود.
سیستم های بیش از حد تعیین شده در بالا را می توان به حداقل مربع برای به دست آوردن ضرایب عرض جغرافیایی حل کرد ، ، ، ، ، و برای معادله ( 18 )، و ضرایب طول جغرافیایی ، ، ، ، ، و برای معادله ( 19 ).
بنابراین، ما از معادلات ( 18 ) و ( 19 ) برای محاسبه مختصات طول و عرض جغرافیایی برای نقاط خاص (به عنوان مثال، نقطه) استفاده کردیم. یا نقاط مجاور آن) به ترتیب در یک صفحه نمایش (یعنی صفحه تصویر).
در معیارهای قبلی مبتنی بر GCA، چهار نقطه در گوشه های یک مربع در صفحه طرح ریزی قرار دارند. استفاده شده. از این رو، ما چهار نقطه را به عنوان انتخاب کردیم ، ، ، ، ، ، و ، ، که با معادلات ( 18 ) و ( 19 ) در صفحه تصویر تصاویر تحریف شده به دست آمد و سپس آنها را به معادلات ( 1 ) و ( 9 ) جایگزین کرد تا معیارهای اعوجاج مبتنی بر تصویر به معنای GCA- محاسبه شود. متریک مبتنی بر پیش بینی نقشه
2.5. پس از پردازش و محاسبه کانتور
چندین برجستگی نقشه از نظر شکل متقارن هستند، به عنوان مثال، تقارن مرکزی برای برخی از برجستگی های آزیموتال، و تقارن افقی و/یا عمودی برای برخی برآمدگی های شبه استوانه ای یا چند مخروطی.
به طور خاص، تقارن شکل فوق با تقارن فرمول های زیرین پیش بینی نقشه تعیین شد. تقارن فرمول های پیش بینی نقشه منجر به تقارن توزیع معیارهای اعوجاج برای پیش بینی های نقشه خاص می شود. از این رو، ما از تقارن توزیع اعوجاج برای بهبود بیشتر کیفیت معیارها استفاده کردیم. با فرض طرح ریزی نقشه با تقارن افقی، مانند طرح ریزی Bonne، تصویر حاصل از متریک اعوجاج مبتنی بر تصویر را به صورت افقی برگردانیدیم و سپس میانگین معیارهای تصاویر برگردانده شده و اصلی را محاسبه کردیم. برای پیش بینی نقشه ها با خاصیت تقارن افقی و عمودی، تصاویر اصلی را به صورت افقی، عمودی ورق زدیم و به صورت افقی و عمودی ورق زدیم. بعد از آن،
علاوه بر این، تمام برجستگیهای استوانهای به شکل مستطیل، با موازیهای مستقیم برای جنبههای عادی بودند. برای پیشبینیهای استوانهای وجه عادی، میتوان نتیجه گرفت که نقاطی در امتداد عرض جغرافیایی یکسان دارای ارزش یکسانی از معیارهای اعوجاج زاویهای و ناحیهای مبتنی بر دیفرانسیل هستند. بنابراین، ما معیارهای مبتنی بر تصویر را در جهت افقی جمعآوری و میانگین کردیم تا دقت معیارهای مبتنی بر تصویر را برای پیشبینیهای استوانهای بهبود ببخشیم.
علاوه بر این، بسیاری از تکنیکهای فیلتر کردن تصویر، مانند فیلترهای میانی و فیلترهای گاوسی، برای کاهش نویز در پردازش تصویر وجود دارد. بنابراین، ما رویکردهای فیلترینگ موجود را در فرآیند محاسبه معیارهای مبتنی بر تصویر اعمال کردیم.
همانطور که معیارهای مبتنی بر تصویر، اعوجاجها را برای نقاط نمونه در صفحه تصویر محاسبه میکنند، ما بیشتر خطوط اعوجاج را در صفحه تصویر تعیین کردیم تا نتایج بصری بصری اعوجاجها را نشان دهیم.
3. نتایج
3.1. مروری بر پیاده سازی و آزمایش ها
NASA G.Projector 3 یک برنامه جاوا و یک نرم افزار نقشه برداری است که از چندین پیش بینی نقشه پشتیبانی می کند [ 27 ]. تصاویر بیش از 200 پیش بینی نقشه در مطالعه ما توسط NASA G.Projector 3.0.2 بر روی دسکتاپ با چهار هسته 64 بیتی 3.3 گیگاهرتز CPU، 16 گیگابایت حافظه DDR3 و 64 بیت Windows 7 Ultimate و 64- تولید شد. بیت OpenJDK 11.
پیش بینی های نقشه فوق با رویکرد پیشنهادی ما ارزیابی شد. الگوریتم اصلی تبدیل مختصات مبتنی بر تصویر، کاهش نویز مبتنی بر تناسب سطح، و محاسبه متریک مبتنی بر GCA در جولیا [ 28 ]، یک زبان برنامهنویسی مدرن، سطح بالا، با کارایی بالا و پویا، برای سرعت آن پیادهسازی شد. و انعطاف پذیری [ 28]. با این حال، محاسبات، تحلیلها و نمودارهای پس پردازش و کانتور به دلیل سادگی و کاربرد آن در پایتون پیادهسازی شدند. نسخه های جولیا و پایتون مورد استفاده به ترتیب 1.6.1 و 3.8.10 بودند. چندین بسته اضافی جولیا، از جمله CSV، DataFrames، Statistics و Images در پیاده سازی ما به کار گرفته شد. چند بسته اضافی پایتون، مانند csv، numpy، scipy، scikit-image، و matplotlib نیز در آزمایشهای ما استفاده شد. خطوط در آزمایش ما توسط matplotlib محاسبه شد.
معیارهای مبتنی بر دیفرانسیل به عنوان مرجع معیارهای پیشنهادی ما استفاده شد. نتایج معیارهای مبتنی بر دیفرانسیل برای بیش از 60 پیش بینی نقشه از مطالعه قبلی [ 20 ] به دست آمد.
3.2. پارامترها و مصرف در آزمایشات
پارامتر برای تصاویر از پیش تعریف شده طول و عرض جغرافیایی استفاده شد. وضوح تصاویر از پیش تعریف شده W = 20000 پیکسل و H بود= 10000 پیکسل تصاویر از پیش تعریف شده به عنوان تصاویر ورودی برای NASA G.Projector 3 استفاده شد. سپس NASA G.Projector تصاویر خروجی تحریف شده را بدون حاشیه، بدون گرتیکول و بدون پوشش با پس زمینه سفید برای پیش بینی نقشه ارائه می کند. رزولوشن تصاویر خروجی که ممکن است دارای پسزمینه سفید غیرضروری باشد نیز 20000 × 10000 پیکسل انتخاب شده است. فرمت های تصاویر ورودی و خروجی برای NASA G.Projector، گرافیک شبکه قابل حمل است. تقریباً 11 گیگابایت حافظه برای ماشین مجازی جاوا لازم است تا برنامه NASA G.Projector 3 را برای تولید وضوح تصویر فوق اجرا کند. مصرف زمان برای تولید تصاویر به پیچیدگی محاسباتی پیش بینی نقشه بستگی دارد. معمولا، به طور متوسط حدود 130 ثانیه برای تولید تصاویر طول و عرض جغرافیایی برای هر طرح ریزی نقشه با استفاده از یک هسته واحد CPU در محیط آزمایش ما مصرف شد. اندازه متوسط تصاویر تحریف شده در طول و عرض جغرافیایی برای هر طرح ریزی نقشه 278 مگابایت بود.
تصاویر خروجی 20000 × 10000 تحریف شده از ناسا G.Projector 3 توسط الگوریتم پیشنهادی در جولیا با پشتیبانی از چندین رشته دنبال میشود. ما هر پیکسل را از بین 20000 × 10000 نمونه در تصاویر برای کارایی ارزیابی نکردیم و فقط از 10 پیکسل به عنوان فاصله پیکسل ها در هر دو جهت عمودی و افقی استفاده کردیم. به طور خاص، تصاویر خروجی نشاندهنده معیارهای مبتنی بر تصویر پیشنهادی تقریباً 20000 × 10000 پیکسل برای هر پیشبینی نقشه است.
در آزمایش خود، پارامتر را انتخاب کردیم ، و سپس استفاده می شود به عنوان اندازه پنجره برای اتصال سطحی مبتنی بر تابع چند جمله ای درجه دوم (به معادلات ( 20 ) و ( 22 ) مراجعه کنید). ضرایب در معادله ( 7 ) به عنوان مقادیر مشابه در مطالعه قبلی [ 20 ] انتخاب شدند. به این معنا که، و . دقت ممیز شناور مورد استفاده در جولیا 64 بیتی است. برای پردازش هر طرح ریزی نقشه در جولیا، تقریباً 4 گیگابایت حافظه مورد نیاز بود. زمان مصرف معمولی تقریباً 208 ثانیه برای هر طرح نقشه با استفاده از 4 هسته CPU با رشته های متعدد پشتیبانی شده در جولیا در محیط آزمایشی ما بود. زمان پردازش ممکن است با شتاب مبتنی بر GPU به دلیل شکل ساده جبر خطی، که در متریک ما استفاده میشود، کاهش یابد.
3.3. کیفیت کاهش نویز با استفاده از فیتینگ سطحی
در این بخش، نمونه ای از طرح ریزی Bonne را برای مقایسه اثر کاهش نویز مبتنی بر اتصالات سطحی ارائه می دهیم.
3.3.1. متریک مبتنی بر دیفرانسیل به عنوان متریک مرجع
شکل 4 توزیع کلاسیک حداکثر اعوجاج زاویه ای مبتنی بر دیفرانسیل را برای طرح ریزی Bonne با موازی استاندارد نشان می دهد. با نوار اسکالر ما از همان نوار اسکالر در شکل های زیر برای نمایش اعوجاج های زاویه ای پیش بینی های مختلف نقشه استفاده کردیم. نمای بزرگ شده ای از توزیع پیش بینی Bonne نیز در شکل 4 نشان داده شده است. خطوط (خطوط سیاه را ببینید) در سطوح ، ، ، ، ، و ، همچنین در شکل 4 ترسیم شده است. حداکثر اعوجاج زاویه ای مبتنی بر دیفرانسیل به عنوان متریک مرجع برای متریک مبتنی بر تصویر پیشنهادی ما برای پیش بینی های مختلف نقشه استفاده شد.
3.3.2. اعوجاج زاویه ای مبتنی بر تصویر بدون کاهش نویز
شکل 5 توزیع ها و خطوط متریک های زاویه ای رویکرد مبتنی بر تصویر را نشان می دهد که تنها از چهار نقطه نمونه در فضای تصویر برای محاسبه متریک هر نقطه استفاده می کند. بنابراین، کاهش نویز مبتنی بر اتصال سطحی در شکل 5 الف وجود ندارد. از آنجایی که نویز شدید در شکل 5 الف وجود دارد، خطا بین متریک مبتنی بر دیفرانسیل و متریک مبتنی بر تصویر بدون کاهش نویز بسیار زیاد است و همبستگی بین آنها به طور قابل توجهی کم است.
شکل 5 ب تقارن افقی برآمدگی Bonne را در نظر می گیرد و از یک فیلتر میانه با دیسکی به شعاع 5 پیکسل استفاده می کند. در شکل 5 ب، شش رنگ مختلف (قرمز، زرد، سبز، فیروزهای، آبی و سرخابی) مقادیر متفاوتی را نشان میدهند. ، ، ، ، ، و ) از خطوط. همانطور که در شکل 5 ب نشان داده شده است، فیلتر میانه نویز متریک را تا حدی کاهش می دهد، اما خطوط به دست آمده از متریک مبتنی بر تصویر و فیلتر میانه هنوز از خطوط به دست آمده از متریک مبتنی بر دیفرانسیل منحرف می شوند.
یک فیلتر گاوسی با انحراف استاندارد 3 بعد از فیلتر میانه اعمال کردیم. نتایج فیلترهای میانه و گاوسی در شکل 5 نشان داده شده است. فیلتر گاوسی نویز را بیشتر کاهش می دهد و شکل خطوط را بهبود می بخشد. با این حال، خطوط به دست آمده از معیارهای مبتنی بر تصویر و فیلترهای میانه و گاوسی هنوز از خطوط به دست آمده از معیارهای مبتنی بر دیفرانسیل منحرف می شوند. علاوه بر این، زمانی که فیلتر گاوسی اعمال می شود، برخی از مقادیر در مرز (به گوشه سمت راست بالا در شکل 5 ج) پیش بینی نقشه گم شده اند.
ما خطوط رنگی را برای نتایج در شکل 5 a ترسیم نکردیم، زیرا خطوط در سطوح مختلف با هم همپوشانی دارند. در شکل 6 ، با استفاده از رنگ های مختلف، از نمودارهای جداگانه برای شش سطح مختلف از خطوط استفاده کردیم. شکل 6 از همپوشانی خطوط اجتناب می کند، اما اثر بصری نتایج در شکل 6 ضعیف است.
از نتایج در شکل 5 و شکل 6 ، معیارهای مبتنی بر تصویر بدون کاهش نویز نمیتوانند توزیعها و خطوط خوبی ایجاد کنند، حتی اگر برخی از فیلترهای تصویر اعمال شوند.
3.3.3. نتایج کاهش نویز مبتنی بر فیتینگ
سپس، یک تابع چند جمله ای درجه دوم برای کاهش نویز مبتنی بر اتصالات سطحی اعمال شد. هنگامی که کاهش نویز اعمال شد، نتایج توزیع ها و خطوط متریک مبتنی بر تصویر به طور قابل توجهی بهبود یافت ( شکل 7 را ببینید ).
شکل 7 a نتیجه معیارهای مبتنی بر تصویر با کاهش نویز مبتنی بر اتصالات سطحی است. کاهش نویز مبتنی بر اتصالات سطحی، چندین نقطه نمونه را در داخل یک پنجره به کار میگیرد تا دقت محاسبات مختصات طول و عرض جغرافیایی را بهبود بخشد. به طور خاص، ما استفاده کردیم پیکسل ها در تصاویر برای محاسبه مختصات و معیارهای مبتنی بر تصویر ( شکل 7 a را ببینید).
علاوه بر این، ما تقارن طرح ریزی Bonne را در نظر گرفتیم و نتایج در شکل 7 ب نشان داده شده است. هنگام در نظر گرفتن تقارن، کانتورها اندکی بهبود می یابند. هنگامی که یک فیلتر میانه با دیسک شعاع 5 بیشتر اعمال می شود، می توان نتیجه بهبود بصری بیشتری به دست آورد ( شکل 7 ج را ببینید).
اگرچه شکل خطوط در شکل 7 c بهبود یافته است، خطا کمی افزایش می یابد و برخی از پیکسل های گم شده در مرز شکل 7 c وجود دارد. بنابراین، یک رویکرد عملی این است که میتوانیم متریکهای مبتنی بر تصویر را محاسبه کنیم و توزیع معیارها را با استفاده از کاهش نویز مبتنی بر اتصالات سطحی و در نظر گرفتن تقارن پیشبینیهای نقشه ترسیم کنیم. هنگامی که خطوط را تعیین می کنیم، یک فیلتر تصویر، مانند فیلتر میانی، می تواند اعمال شود.
3.4. همبستگی و تحلیل خطای معیارها
با استفاده از طرح ریزی Bonne به عنوان مثال برای نشان دادن اثربخشی متریک مبتنی بر تصویر با کاهش نویز مبتنی بر اتصالات سطحی، ما همبستگی و خطای متریک پیشنهادی را برای پیش بینی های بیشتر نقشه ارزیابی کردیم.
ما از معیارهای مبتنی بر دیفرانسیل به عنوان مرجع معیارهای پیشنهادی خود استفاده کردیم. ضریب همبستگی محصول- لحظه پیرسون (PPMCC) و ریشه میانگین مربع خطا (RMSE) برای ارزیابی همبستگی و خطا بین متریک پیشنهادی ما و متریک مبتنی بر دیفرانسیل استفاده شد.
در مطالعه قبلی، معیارهای مبتنی بر دیفرانسیل برای بیش از 80 پیش بینی محاسبه شد. در میان آنها، 64 پیش بینی نقشه تقریباً یک میلیون نقطه نمونه را در صفحه طرح ریزی برای به دست آوردن متریک مبتنی بر دیفرانسیل به کار گرفتند. در این بخش، نتایج این 64 پیش بینی نقشه را با هم مقایسه می کنیم.
توجه داشته باشید که معیارهای مبتنی بر تصویر پیشنهادی تقریباً استفاده میشوند نقاط نمونه (ممکن است حاوی پس زمینه سفید غیر ضروری باشد). در همین حال، نتایج یک مطالعه قبلی [ 20 ] از حدود یک میلیون نقطه نمونه استفاده می کند. بنابراین، ما اندازه نتایج متریک مبتنی بر دیفرانسیل را به اندازه نتایج متریک مبتنی بر تصویر برای مقایسه بصری و تجزیه و تحلیل آماری تغییر دادیم.
جدول 1 نتایج آماری 64 پیش بینی نقشه را ارائه می دهد. هیستوگرام خطاها در شکل 8 نشان داده شده است . نتایج آماری از یک فیلتر گاوسی با انحراف معیار 1 استفاده میکند. از جدول 1 و شکل 8 ، میتوان مشاهده کرد که بدترین خطا تقریباً است. برای پیش بینی های Ortelius Oval و Winkel II، در حالی که در اکثر موارد (تقریباً 89٪)، خطا کمتر از .
از آنجایی که خطوط، شاخصهای بصری همبستگی هستند و PPMCC و RMSE فقط شاخصهای آماری همبستگی یا خطا هستند، باید توزیعها و خطوط را برای برخی پیشبینیها با خطاهای بزرگ رسم کنیم تا خطا بین متریک مبتنی بر دیفرانسیل و تصویر پیشنهادی ما را تأیید کنیم. متریک مبتنی بر
شکل 9 نتایج اندازه گیری های دیفرانسیل و مبتنی بر تصویر را برای پیش بینی های Ortelius Oval و Winkel II و برخی از پیش بینی های نقشه نشان می دهد که خطاهای آنها بیشتر از . ما میتوانیم مشاهده کنیم که خطوط متریک مبتنی بر تصویر پیشنهادی تا حدی به خطوط بهدستآمده از متریک مبتنی بر دیفرانسیل منحرف میشوند، بهویژه در مناطقی که اعوجاجهای کمتری دارند (برای بیشتر پیشبینیهای نقشه، خطوط قرمز را ببینید، و برای بیکن گلوبولار طرح ریزی)؛ در حالی که خطوط متریک مبتنی بر تصویر پیشنهادی به خطوط به دست آمده از متریک مبتنی بر دیفرانسیل در مناطقی با اعوجاج بالاتر نزدیک است، که مورد توجه محققان یا نقشهسازان است.
با توجه به نتایج بصری و آماری فوق، میتوانیم از متریک مبتنی بر تصویر پیشنهادی برای ارزیابی پیشبینیهای بیشتر نقشه استفاده کنیم. شکل 10 نتایج متریک مبتنی بر تصویر را برای برخی از پیش بینی های کانترز، پیش بینی های ماورر و غیره نشان می دهد که در مطالعه قبلی ما ذکر نشده است، به عنوان کتابخانه PROJ [ 27 ، 29 ]، که در مطالعه قبلی [ 20 ] استفاده شد.]، از این پیش بینی های نقشه پشتیبانی نمی کند. از آنجایی که NASA G.Projector از بیش از 200 پیش بینی نقشه و پیش بینی های نقشه بیشتری نسبت به کتابخانه PROJ پشتیبانی می کند، می توانیم از روش مبتنی بر تصویر پیشنهادی برای ارزیابی این پیش بینی های نقشه موجود استفاده کنیم. علاوه بر این، G.Projector شامل تمام پیش بینی های نقشه نیست، بنابراین می توانیم از سایر نرم افزارهای نقشه برداری مانند Geocart به عنوان مکمل برای ارزیابی بیشتر پیش بینی های نقشه استفاده کنیم.
4. بحث و نتیجه گیری
در این مطالعه، ما یک متریک مبتنی بر تصویر جدید برای ارزیابی اعوجاج پیشبینیهای نقشه پیشنهاد کردیم و از یک رویکرد کاهش نویز مبتنی بر اتصالات سطحی جدید برای امکانپذیر ساختن متریک مبتنی بر تصویر برای ارزیابی اعوجاج استفاده کردیم. روش پیشنهادی عملی است زیرا میتواند از نرمافزار نقشهبرداری موجود برای تولید تصاویر تحریف شده برای پیشبینیهای نقشه استفاده کند، در حالی که از فرمولهای صریح پیشبینی نقشه برای محاسبه اعوجاج پیشبینیهای نقشه استفاده نمیشود. از آنجایی که NASA G.Projector که به عنوان نرمافزار نقشهبرداری استفاده میشود، بسیاری از پیشبینیهای نقشه را با چندین پارامتر قابل تنظیم ارائه میکند، و نرمافزارهای نقشهبرداری دیگری نیز در دسترس هستند، میتوانیم از متریک مبتنی بر تصویر پیشنهادی برای ارزیابی تعداد کافی از پیشبینیهای نقشه موجود استفاده کنیم.
با این حال، نرمافزار نقشهبرداری موجود نمیتواند شامل تمام پیشبینیهای نقشه، بهویژه برخی از پیشبینیهای نقشه جدید یا اخیراً طراحی شده باشد. استفاده از نرمافزار نقشهبرداری برای محاسبه متریک مبتنی بر تصویر پیشنهادی برای این پیشبینیهای نقشه جدید ممکن است در قابلیت استفاده با مشکلاتی مواجه شود. به عنوان یک جایگزین، میتوانیم به صورت دستی تصاویر طول و عرض جغرافیایی را با معادلات رو به جلو یا معکوس پیشبینیهای نقشه خاص تولید کنیم. تصاویر تولید شده را می توان برای محاسبه متریک مبتنی بر تصویر استفاده کرد. از این رو، متریک مبتنی بر تصویر پیشنهادی یک رویکرد عملی و قابل استفاده برای ارزیابی طرح ریزی نقشه است.
همچنین باید توجه داشت که ناسا G.Projector شامل برخی از پیشبینیهای نقشه موجود نیست، به عنوان مثال، Eckert I، Eckert II، Collignon و Tobler-Mercator. به عنوان یک مکمل، Geocart’s Mapthematics LLC از این پیشبینیهای نقشه و موارد دیگر پشتیبانی میکند. بنابراین، ما می توانیم بیشتر از Geocart به عنوان نرم افزار نقشه برداری استفاده کنیم و از روش پیشنهادی خود برای ارزیابی پیش بینی های بیشتر نقشه استفاده کنیم.
علاوه بر این، برای پیشبینیهای نقشه با فرمولهای شناخته شده، یکی از بهبودهای احتمالی متریک پیشنهادی، استفاده از محاسبات مستقیم بر اساس گروهی از نقاط نمونه و فرمولهای پیشبینی نقشه است، در حالی که استفاده از تصاویر تحریف شده را کنار میگذارد. این بهبود ممکن است زمان پردازش و مصرف حافظه را کاهش دهد.
الگوهای خاصی برای تصاویر طول و عرض جغرافیایی استفاده می شود. ما فکر می کنیم الگوهای دیگری که دارای نقشه برداری یک به یک هستند نیز می توانند به عنوان الگوهای از پیش تعریف شده استفاده شوند. در ابتدا، ما استفاده از یک الگوی از پیش تعریف شده را برای تبدیل طول و عرض جغرافیایی در نظر گرفتیم. با این حال، ما دریافتیم که جدا کردن طول و عرض جغرافیایی با استفاده از یک الگو دشوار است. سپس الگوهای جداگانه ای را برای تصاویر طول و عرض جغرافیایی در نظر گرفتیم. نتایج در دستنوشته ما حاکی از آن است که الگوهای اتخاذ شده برای ایجاد رابطه نگاشت اولیه امکانپذیر هستند، و عدم دقت در تبدیل مبتنی بر تصویر را میتوان با کاهش نویز مبتنی بر اتصالات سطحی تصحیح کرد. یعنی ممکن است الگوهای قابل اعتمادتری وجود داشته باشد، اما الگوهای فعلی برای متریک پیشنهادی ما کافی است.
در این مطالعه، ما از NASA G.Projector برای تولید تصاویر برای بیش از 200 پیش بینی نقشه، و محاسبه اعوجاج برای آن پیش بینی های نقشه استفاده کردیم. علاوه بر این، ما معیارهای مبتنی بر تصویر پیشنهادی خود را با معیارهای مبتنی بر دیفرانسیل مقایسه کردیم. نتایج تحلیل و مقایسه نشان داد که متریک پیشنهادی همبستگی بالا و خطای کمی با متریک مبتنی بر دیفرانسیل موجود دارد. مشاهدات دیگر در مورد متریک مبتنی بر تصویر پیشنهادی، تفاوت های قابل توجه بین متریک مبتنی بر دیفرانسیل و متریک مبتنی بر تصویر پیشنهادی در مناطقی با اعوجاج کمتر است. از آنجایی که نواحی با اعوجاج بیشتر مورد توجه محققین یا نقشه سازان هستند، ممکن است برخی از خطاهای موجود در مناطق با اعوجاج کمتر قابل قبول باشد. در آینده، بهبود معیارهای قبلی مبتنی بر GCA را در نظر خواهیم گرفت،
اعوجاج ناحیه نیز یک معیار مهم برای پیش بینی نقشه است. همانطور که در مطالعه قبلی نشان داده شد، متریکهای مساحت بسیار سادهتر از معیارهای زاویهای هستند، و متریکهای منطقه مبتنی بر GCA و FWD-GCA تقریباً همان معیارهای متریک منطقه مبتنی بر تفاضل هستند. بنابراین، ما متریک اعوجاج منطقه را با استفاده از روش مبتنی بر تصویر مورد بحث قرار ندادیم.
بدون دیدگاه