فهرست مطالب

نشریه علوم و فنون نقشه برداری
سال ششم شماره 2 (پیاپی 22، پاییز 1395)

  • تاریخ انتشار: 1395/11/16
  • تعداد عناوین: 10
|
  • رحیم علی عباسپور، سید حسن حسینی صفحات 1-11
    کنترل و تضمین کیفیت اطلاعات مکانی داوطلبانه یک از بزرگترین چالش های پیش روی این نوع از گردآوری اطلاعات مکانی است. روش های ارائه شده برای ارزیابی کیفیت اطلاعات مکانی داوطلبانه به صورت کلی شامل دو رویکرد مقایسه این اطلاعات با یک مرجع معتبر و یا روش های بدون نیاز به مقایسه با مرجع معتبر می شود. روش های گروه اول نیازمند دسترسی به اطلاعات مرجع هستند که دسترسی به این اطلاعات همواره امکان پذیر نمی باشد. روش های گروه دوم، یا مثل روش اعتبارسنجی اطلاعات سعی بر به دست آوردن اعتبار اطلاعات بر اساس معیارهای بدست آمده از خود اطلاعات دارند و یا با استفاده از قوانین حاکم بر داده ها مثل سازگاری منطقی داده ها به ارزیابی کیفیت داده ها می پردازند. در این تحقیق تلاش بر ارائه روشی برای ارزیابی کیفیت اطلاعات مکانی داوطلبانه با استفاده از خصوصیات قابل محاسبه از داده مکانی و بدون نیاز به مقایسه با اطلاعات مرجع گردیده است. این تحقیق با استفاده از خصوصیات قابل محاسبه از داده های موجود در شبکه راه های OpenStreetMap، به صورت خاص به تحلیل میزان دقت دسته بندی خیابان ها، به عنوان بخشی از دقت معنایی مربوط به شبکه راه ها، می پردازد. در این تحقیق از درخت تصمیم و شبکه عصبی با ساختار چند لایه برای یادگیری این خصوصیات استفاده گردید. برای تعیین میزان دقت دسته بندی بدست آمده، معیار های دقت و بازخوانی برای هر روش محاسبه شد. درخت تصمیم با دقت وزن دار بدست آمده 84٫1٪ نشان دهنده یک مدل مناسب برای این روش است. این روش بر اساس اطلاعات قابل استخراج از داده ها استوار است و می تواند برای داده های هر شهر دیگری تعمیم و استفاده گردد.
    کلیدواژگان: OpenStreetMap، کیفیت داده، داده مکانی داوطلبانه، آنالیز شبکه راه، ماشین یادگیری
  • محمدعلی برومند، عبدالرضا صفری، یحیی الله توکلی، عباس بحرودی صفحات 13-25
    تانسور گرادیان گرانی، ماتریسی مشتمل بر مشتقات مرتبه ی دوم گرانی می باشد که کاربردهای فراوانی درعلوم ژئودزی و ژئوفیزیک دارد. تا به امروز تلاش های زیادی برای اندازه گیری و برآورد تانسور گرادیان گرانی با دقت مناسب صورت گرفته است. این کمیت تانسوری از روش های متنوعی قابل برآورد بوده که یکی از این روش ها استفاده از مشاهدات ثقلی و به کارگیری روش تفاضل محدود است. روش تفاضل محدود با استفاده مستقیم از مفهوم ریاضی تانسور گرایان، بدون در نظر گرفتن فرضیات اضافی که ممکن است از واقعیت مسئله به دور باشد، می‎تواند از نقطه نظر تئوری امکان برآورد صحیح تانسور گرادیان گرانی را فراهم ‎آورد. در این مطالعه، تانسور گرادیان با استفاده از داده های جاذبی زمینی و روش تفاضل محدود فرموله و تعیین شده است. در این مقاله، روابط پیشنهادی به صورت عددی با استفاده از یک مدل جاذبی جهانی زمین مورد آزمایش قرار می‎گیرد و سپس از این روش، جهت برآورد تانسور گرادیان گرانی در 6350 ایستگاه ثقلی در منطقه ی فارس ساحلی در جنوب ایران، استفاده می‎گردد. آزمایشهای عددی این تحقیق نشان از توانایی روش در برآورد مولفه های تانسور گرادیان گرانی با دقت مناسب می‎باشد. به عنوان مثال در آزمایشهای عددی دیده می‎شود که روش پیشنهادی این تحقیق امکان برآورد مولفه های قطری تانسور گرادیان (به عبارتی مشتقات دوم گرانی در جهت‎های شرق، شمال و قائم) را با دقت‎هایی به ترتیب برابر با 12.46، 34.49 و 454.82 اتووش فرآهم نموده است.
    کلیدواژگان: تانسور گرادیان گرانی، روش تفاضل محدود، مدل ژئوپتانسیل، فارس ساحلی
  • مهدی اصفهانی، علی محمدزاده صفحات 27-42
    شناسایی عوارض شهری مهم مانند ساختمان و درخت از تصاویر هوایی و ماهواره ای به همراه داده های لایدار توجه محققان سنجش ازدور و فتوگرامتری را در چند دهه اخیر به خود متوجه ساخته است. روش های طبقه بندی و شناسایی مختلف به منظور استفاده در مناطق شهری با مشکلات پیچیده ای مانند شناسایی ساختمان و درخت های کوچک، مرز نامناسب تاج درختان، وجود سقف هایی از ساختمان ها با پوشش گیاهی، ساختمان ها محصورشده توسط درختان و پوشش های گیاهی قرارگرفته در سایه، مواجه هستند. در این تحقیق جهت بهبود مشکلات ذکر شده در مرحله اول پیش پردازش داده های لایدار، تولید مدل رقومی زمینی و مدل رقومی سطح نرمال شده انجام گردیده است؛ در مرحله دوم ویژگی های هندسی، طیفی، بافتی و ترکیبی متناسب با نوع مشکلات ذکر شده تولید و با استفاده از الگوریتم SVM_GA ویژگی های مناسب انتخاب گردیده است. از میان ویژگی ها تولیدی شاخص طیفی ترکیب یافته پوشش گیاهی با مناطق سایه و ویژگی های هندسی تولیدشده از داده های لایدار جهت شناسایی درختان نقش بسزایی بازی می کنند. در مرحله سوم طبقه بندی ماشین بردار پشتیبان در دو سطح شیءگرا و پیکسل مبنا با استفاده از ویژگی ها منتخب و مناسب با استفاده از روش SVM_GA به منظور شناسایی عوارض شهری ازجمله درخت و ساختمان های شهری استفاده گردیده است. در مرحله چهارم نتایج هریک از دو سطح طبقه بندی با استفاده از روش های پس پردازش بهبود داده شده اند. در مرحله پنجم با توجه به قدرت هر سطح از طبقه بندی در شناسایی عوارض سعی در بهبود شناسایی مرز ساختمان ها، درختان و شناسایی ساختمان ها و درختان کوچک و کم ارتفاع را با استفاده از تلفیق نتایج آن ها را داریم و درنهایت نتایج حاصله از شناسایی در هر سطح و هر مرحله با داده های مرجع مقایسه و ارزیابی شده است که شناسایی ساختمان ها و درختان در سطح شیءگرا به ترتیب دارای صحت شناسایی 961/0 و 65/0 و دقت کلی 97/0 و 943/0 هستند. در سطح پیکسل مبنا دارای صحت شناسایی 953/0 و 632/0، دقت کلی 961/0 و 94/0 و بعد از تلفیق دارای صحت شناسایی971/0 و 718/0 و دقت کلی 975/0 و 957/0 می باشند. نتایج نشان می دهد که تلفیق دو سطح از طبقه بندی شیءگرا و پیکسل مبنا باعث بهبود نتایج شناسایی به ویژه در شناسایی ساختمان ها و درختان کوچک و کم ارتفاع می گردد. نتایج این روش با توجه به بهبود شناسایی های انجام شده در مرز ساختمان ها و شناسایی عوارض کوچک تاثیر مثبتی را می تواند در استخراج عوارض داشته باشد.
    کلیدواژگان: شناسایی ساختمان و درختان شهری، طبقه بندی، تلفیق، شیءگرا، پیکسل مبنا، داده های لایدار و تصاویر هوایی
  • یحیی الله توکلی، عبدالرضا صفری، علیرضا آزموده ‎اردلان، عباس بحرودی صفحات 43-62
    هدف اصلی این تحقیق آن می‎باشد که با تلفیق مسائل معکوس گراویمتری ژئودزی و ژئوفیزیک امکان مدلسازی همزمان «میدان ثقل» و «توزیع دانسیته جرمی» برروی سطح زمین را فراهم آورد. در فیزیکال ژئودزی هدف از طرح مسئله معکوس گراویمتری، مدلسازی دقیق تابعک پتانسیل ثقل و مشتقات آن می‎باشد و در آن طرف در مباحث ژئوفیزیک اکتشافی، عموما مسئله معکوس گراویمتری با هدف تعیین توزیع دانسیته جرمی زمین طرح می‎گردد؛ در این تحقیق می‎بینیم که از دو منظر ریاضی و فیزیک، وصول به این دو هدف می‎تواند در حل یک مسئله معکوس گراویمتری تلفیقی، خلاصه و محقق گردد. بدین منظور، تئوری‎، مفاهیم و تشابه این دو نوع مسئله معکوس گراویمتری مورد بررسی قرار گرفته و ارتباط بین پارامترهای مجهول این دو نوع مسئله معکوس گراویمتری مشخص خواهد گردید. سپس روشی برای مدلسازی تابعک پتانسیل میدان جاذبه برروی سطح زمین، بر مبنای توابع پایه شعاعی غیرهارمونیک و پایدارسازی تیخونوف تعمیم‎یافته پیشنهاد می‎گردد و با استفاده از روش پیشنهادی و تلفیق معادله دیفرانسیلی پواسن و معادله انتگرالی نیوتن امکان استخراج تابعک دانسیته از مشاهدات ثقلی زمینی فراهم می‎آید. در انتها نیز در یک مطالعه موردی، روش پیشنهادشده با استفاده از 6350 ایستگاه ثقلی در منطقه فارس ساحلی مورد آزمایش عددی قرارگرفته می‎شود. نتایج عددی حکایت از آن دارد که روش پیشنهادشده امکان مدلسازی توام میدان جاذبه و توزیع دانسیته جرمی در منطقه مورد مطالعه را به ترتیب با دقتهای 0.54 میلی‎گال و 0.25 گرم‎بر‎سانتیمترمکعب فرآهم نموده است.
    کلیدواژگان: ثقل، دانسیته، مدلسازی، گرانی، چگالی، وارونسازی، مسئله معکوس
  • مهرنوش امتی، محمودرضا صاحبی صفحات 63-78
    عواملی هم چون وجود اسپکل ها در تصاویر SAR، وابستگی زیاد بین پیکسل های همسایه، احتمال رخداد تغییرات در مناطق مجاور نسبت به نقاط دور از یکدیگر و هم چنین دشواری دست یابی به نتایج موردنظر در صورت استفاده از روش های مبتنی بر پیکسل، ضرورت به کارگیری آنالیز شی مبنا را در بهبود دقت تشخیص تغییرات ایجاب می کند. هدف از این مقاله، ارائه روشی نوین در تشخیص تغییرات به هنگام پوشش های اراضی با استفاده از تصاویر پلاریمتریک SAR خواهد بود. در روش پیشنهادی ابتدا هر دو تصویر پلاریمتری اخذ شده از منطقه با استفاده از الگوریتم بهبودیافته آب پخشان به طور جداگانه قطعه بندی می شوند. سپس از میان ویژگی های مختلف پلاریمتری، ده ویژگی تفاضلی بهینه با استفاده از الگوریتم ژنتیک و معیار فاصله Jeffries-Matusita انتخاب گردیده و در مرحله پایانی با درنظرگرفتن اطلاعاتی همانند تفاوت میانگین مقادیر پیکسل های مناطق همگن قطعه بندی شده در دو تصویر، تغییر و یا عدم تغییر قطعات با استفاده از طبقه بندی کننده ماشین بردار پشتیبان (SVM) بررسی می گردد. مقایسه نتایج حاصل از اجرای این الگوریتم با داده های مرجع، مقادیر 40/92 درصد و 85/0 را برای دقت طبقه بندی و ضریب کاپا نشان می دهد. این در حالی است که دقت و ضریب کاپا در روش پیکسل مبنا به ترتیب برابر با 61/78 درصد و 58/0 به دست آمده است. هم چنین این روش قطعه بندی در مقایسه با روش پیکسل مبنا توانسته است مرز های قطعات تصویری را با درصد قابل توجهی حفظ کند.
    کلیدواژگان: تشخیص تغییرات، تصاویر پلاریمتری SAR (PolSAR)، آنالیز شی مبنا، الگوریتم بهبودیافته آب پخشان
  • مهرداد آهنگرکانی، مهدی فرنقی، محمدرضا شیرزادی صفحات 79-98
    تشخیص عوامل بیماری زا، شناسایی تجمع مکانی بیماری و یافتن الگوی انتشار آن در محیط از ضروری ترین نیازها در زمینه بهداشت عمومی و مدیریت بیماری ها به شمار می آیند. بیماری لپتوسپیروز یکی از بیماری های مشترک انسان و دام می باشد که تقریبا در تمام نقاط جهان به ویژه در مناطق گرمسیری، نیمه گرمسیری و نواحی گرم و مرطوب شیوع بیشتری دارد. شرایط آب و هوایی معتدل و مرطوب در استانهای شمالی ایران این مناطق را در معرض خطر بیشتر برای ابتلا به این بیماری قرار داده است. اهداف اصلی این تحقیق بررسی سالانه توزیع مکانی بیماری لپتوسپیروز، کشف خوشه های مکانی و مکانی- زمانی بیماری و تهیه نقشه پیش بینی توزیع مکانی- زمانی آن در استانهای شمالی ایران می باشد. در این تحقیق داده های وقوع بیماری لپتوسپیروز به همراه داده های محیطی و توپوگرافی از قبیل میانگین دمای هوا، میانگین رطوبت، مجموع بارش سالیانه، ارتفاع، شیب، جهت شیب و تعداد روزهای یخبندان در سه استان گیلان، مازندران و گلستان از ابتدای سال 1389 تا پایان سال 1393 به صورت ماهانه جمع آوری شده اند. به منظور بررسی وجود و یا عدم وجود خودهمبستگی مکانی میان موارد وقوع بیماری لپتوسپیروز از شاخص های عمومی اندازه گیری خودهمبستگی مکانی همچون Moran’s I و General G استفاده شده است. آزمون های خودهمبستگی عمومی برای هر 6 سال به صورت جداگانه فرض تصادفی بودن توزیع بیماری در منطقه را رد کرده و بیانگر این است که توزیع بیماری در منطقه در هر 6 سال خوشه ایمی باشد. به منظور کشف خوشه های مکانی و شناسایی مناطق پرخطر بیماری لپتوسپیروز از شاخص های محلی کشف خوشه های مکانی از قبیل Local Moran’s I و Gi Local استفاده شده است. در این تحقیق به منظور تهیه نقشه پیش بینی توزیع مکانی- زمانی بیماری لپتوسپیروز، از دو مدل رگرسیون وزندار جغرافیایی و شبکه عصبی پرسپترون چند لایه استفاده شده است. در نهایت معیارهای ضریب کاپا، RMSE، MAPE و R2 جهت ارزیابی عملکرد مدلهای پیشنهادی به کار گرفته شده اند.
    کلیدواژگان: بیماری لپتوسپیروز، سامانه اطلاعات مکانی، توزیع مکانی، زمانی، رگرسیون وزندار جغرافیایی، شبکه عصبی، ایران
  • سعید نادی، یوسف غیاثی، شهرام هداوند صفحات 99-118
    با گسترش مدل های ارتفاعی رقومی آزاد، کاربرد این مدل ها در زمینه های مختلف در چند سال اخیر رشد روزافزونی داشته است. عدم اطلاع کامل از دقت ارتفاعی این مدل ها در مناطق با ویژگی های مختلف یکی از معضلاتی می باشد که نتایج مطالعات در حوزه های مختلف بر پایه این مدل ها را غیر قابل ارزیابی و اطمینان می کند. نتایج حاصل از مدل های آماری رایج ارزیابی دقت نظیر خطای میانگین، انحراف استاندارد و خطای مجذور میانگین مربعات با فرض پیروی خطاها از توزیع نرمال قابل اتکا می باشند در حالیکه این شرط عمدتا با توجه به احتمال وجود خطاهای بزرگ در داده ها قابل تامین نمی باشد. در این مقاله ابتدا به بررسی رفتار آماری خطای مدل های ارتفاعی زمین GDEM و SRTM پرداخته و با استفاده از تعریف یک محدوده اطمینان برای خطای مجذور میانگین مربعات روشی برای تعیین تعداد نقاط مرجع مناسب ارائه شده است. سپس با استفاده از مدل های آماری رایج و نیز مدل های با استحکام به ارزیابی دقت این مدل ها در مقایسه با نقاط مرجع پرداخته ایم. برای این منظور ابتدا شاخصهای آماری لازم تعریف شده و سپس در چهار کلاس مناطق شهری، مسطح، تپه ماهور و کوهستان در مناطق غربی و مرکزی ایران با محاسبه این شاخصها به مقایسه و تحلیل خطای مدل های رقومی مذکور پرداخته شده است. عدم توجه به ویژگی های متفاوت مناطق مختلف و تاثیری که این ویژگی ها در مراحل مختلف تعیین ارتفاع و تهیه مدل ارتفاعی خواهند داشت یکی از نواقص مطالعات پیشین می باشد که در این مطالعه به آن پرداخته شده است. علاوه بر این بررسی تاثیر خطا بر محصولات ثانویه از طریق محاسبه انتشار خطاها در نقشه های شیب و جهت شیب حاصل از مدل های مذکور از جنبه های نوآورانه بسیار مهمی می باشد که در این تحقیق به آن پرداخته شده است. نتایج تحقیق انجام شده نشان داد که علی رغم قدرت تفکیک پایین تر، دقت SRTM در همه مناطق بهتر از GDEM است. علاوه بر این در هر دو مدل دقت در مناطق شهری به میزان قابل ملاحظه ای بهتر از دیگر مناطق می باشد.
    کلیدواژگان: مدل ارتفاعی رقومی، مقایسه دقت و صحت، مدل های آماری با استحکام، انتشار خطاها، محدوده اطمینان
  • احسان کیانا، سعید همایونی، محمدعلی شریفی، محمدرضا فریدروحانی صفحات 119-130
    در این مقاله روشی برای آشکارسازی بدون نظارت تغییرات در تصاویر سنجش از دوری قطبیده راداری با روزنه مصنوعی ارائه شده است. این روش بر پایه‏ی روش آنالیز نقطه ی تغییر است. تابع چگالی احتمال تصاویر اختلاف، که در برگیرنده ی تغییرات زمانی محیطی هستند، از توزیع‏های آماری مشخصی برای هر کلاس پیروی می کنند. مدل آمیخته گاوسی یکی از مدل های مناسب برای آنالیز نقطه ی تغییر است که توانمندی مناسبی برای برآورد پارامترهای توزیع ها از یک توزیع آمیخته دارد. این مدل می‏تواند نقطه برخورد دو توزیع را به عنوان نقطه ی تغییر دو کلاس تغییریافته و تغییرنیافته، یا همان حد آستانه را به راحتی مشخص کند. روش پیشنهادی بر روی سه داده راداری پیاده سازی و مورد ارزیابی قرار گرفت. نقشه تغییر حاصل از دو مجموعه‏ی داده در این تجزیه وتحلیل با استفاده از داده مرجع، دقت کلی برابر 82.97 و 96.34 درصد را نشان دادند. داده‏ی چند زمانه با قطبیدگی راداری کامل این پژوهش از یک منطقه کشاورزی بود. اطلاعات زمینی و اطلاعات تجزیه ماتریس پراکنش با دقت خوبی ارتباط این تغییرات را با فعالیت‏های کشاورزی مانند مراحل رشد گیاه نشان داد. همچنین روش پیشنهادی با یکی از کاراترین روش‏های برآورد حد آستانه، یعنی روش Otsu مقایسه شد که نتایج به طور ویژه دقت بالاتری را نشان داد.
    کلیدواژگان: آشکارسازی بدون نظارت تغییرات، تصاویر تمام قطبیده راداری با روزنه مصنوعی، آنالیز نقطه ی تغییر، مدل آمیخته گاوسی
  • محمد رستم نیا، مهدی آخوندزاده هنزایی صفحات 131-144
    این تحقیق برای بررسی خشکیدگی درختان جنگلی ناحیه غرب کشور در استان ایلام در یک دوره ی 15 ساله با استفاده از تصاویر ماهواره ای لندست انجام شده است. برای بارزسازی پوشش گیاهی در این تصاویر ، پنج شاخص پوشش گیاهی مورد ارزیابی قرار گرفت و با استفاده از نقاط نمونه برداری شده از روی تصاویر گوگل ارث دقت شاخص ها محاسبه گردید. شاخص گیاهی EVI با همبستگی بالا با پوشش جنگلی ، نسبت به سایر شاخص ها، به عنوان شاخص بارز کننده پوشش گیاهی در این تحقیق مورد استفاده قرار گرفت. شاخص EVI برای این دوره ی 15 ساله محاسبه شد و میزان همبستگی بین تغییرات EVI ناحیه ی جنگلی استان ایلام و تغییرات شاخص هواشناسی SPI دوره ی یک ساله و سه ماهه و نه ماهه مورد بررسی قرار گرفت که نتایج حاصل بترتیب برابر با 46/0، 09/0 و 73/0 بدست آمد. در ادامه این تحقیق، با روش PCA تغییرات هر سال نسبت به سال 1379 محاسبه شد و همچنین با بررسی رفتار تغییرات هر پیکسل در بازه ی زمانی 15 ساله، تغییرات ناشی از خشکیدگی درختان از سایر تغییرات، تفکیک شد و همچنین میزان این تغییرات برای هر سال محاسبه شد و میزان تاثیرگذاری تغییرات بارش و ریزگردها در افزایش روند خشکیدگی درختان جنگلی استان ایلام، با انجام یک رگرسیون خطی دو متغییره بدست آمد که نتایج سهم 62 درصدی تغییرات بارش و 38 درصدی گردوخاک را نشان می دهند.
    کلیدواژگان: خشکیدگی درختان، شاخص پوشش گیاهی، SPI، PCA، ایلام
  • ناصر طهانی، علی اصغر آل شیخ، مهدی فرنقی صفحات 145-160
    مکان یابی و نمایش محیط اطراف کاربر در دستگاه های همراه به طور معمول در قالب نقشه های دو بعدی امکان پذیر است که کاربر را به یک تصویر عمودی یا نمایش بالا به پایین از محیط محدود می کند. با توجه به ماهیت سه بعدی عوارض دنیای واقعی، ارائه آنها به صورت مدل های سه بعدی نقش مهمی را در بصری سازی و انجام انواع تحلیل ها ایفا می کند. در بخش بصری سازی و نمایش، مدل های سه بعدی شهری با ایجاد محیط مجازی و در نظرگرفتن خصوصیات هندسی و ظاهری عوارض، به درک بهتر کاربران کمک می کنند. علاوه بر این مدل های شهری در طیف وسیعی از کاربردها همانند تحلیل محدوده قابل رویت، ارزیابی میزان هدررفت انرژی، مدیریت تاسیسات شهری، کاداستر سه بعدی و ناوبری در محیط های داخلی به کار گرفته می شوند. پیشرفت های دستگاه های همراه منجر به استفاده آنها در کاربردهای مختلف نظیر فناوری واقعیت مجازی و مدل های سه بعدی شهری شده است. در این مطالعه به منظور رفع محدودیت های نقشه های دو بعدی و استفاده از قابلیت های نمایش عوارض به صورت سه بعدی، یک برنامه کاربردی موبایل در دو رویکرد تعاملی چشم پرنده[1] و زاویه دید اول شخص توسعه پیدا کرده است. مدل های سه بعدی به کمک ردپای دو بعدی عوارض حاصل از نقشه 1:2000 منطقه و افزودن مقادیر ارتفاع و بافت نمایشی مربوطه در محیط نرم افزار Google Sketchup تهیه شده اند. با توجه به محدودیت های رابط برنامه نویسی کاربردی[2] OpenGL ES در بارگذاری مدل های گرافیکی پیچیده، کتابخانه گرافیکی Libgdx به منظور تشکیل محیط و نمایش مدل ها استفاده شده است. فرآیند Deserialization به منظور تشکیل خودکار مدل های سه بعدی، اتصال بافت مربوط به عوارض و هم چنین جانمایی آنها در محیط برنامه صورت گرفته است. در نمای چشم پرنده کاربر به کمک دکمه های طراحی شده و لمس صفحه نمایش به گردش در محیط می پردازد. در رویکرد ناظر اول شخص یا واقعیت مجازی، موقعیت کاربر به وسیله GPS دستگاه و راستای نشانه روی با استفاده از سنجنده های مکانی شتاب سنج و مغناطیس سنج و اعمال فیلتر پایین گذر[3] به منظور کاهش خطای داده ها تعیین می شود . رویکرد تعاملی ترکیبی پیاده سازی شده در این تحقیق با فراهم نمودن قابلیت نمایش به صورت نقشه های سه بعدی و ساختار واقعیت مجازی، باعث تمایز آن از سایر مطالعات مربوط به بصری سازی عوارض شهری شده است. به منظور ارزیابی برنامه کاربردی، عوارض محوطه دانشکده مهندسی نقشه برداری دانشگاه صنعتی خواجه نصیرالدین طوسی در برنامه وارد شده و امکان مشاهده مدل های سه بعدی و تعامل کاربران با آن در دو رویکرد نمایشی فراهم شده است.
    1 Bird''s Eye View
    2 Application Programming Interface
    3 Low pass filter
    کلیدواژگان: واقعیت مجازی، مدل سه بعدی شهری، کتابخانه گرافیکی Libgdx، نمای چشم پرنده، زاویه دید اول شخص
|
  • R. Ali Abbaspour, S. H. Hossaini Pages 1-11
    The geospatial data quality validation is building an effort to develop quality standards for geospatial information. The growth of Volunteer Geospatial Information (VGI) raises many issues regarding quality control or quality assurance. The lack of metadata and standards combined with an unknown motivation and used credibility lead to heterogeneous datasets with unknown quality. Thus¡ two main approaches have been suggested for assessing the quality of VGI: first approach is based on comparing VGI with authorized data and second approach tries to validate VGI without using authorized data.
    The methods based on first approach need the authorized information and this information is not always accessible or in case of existence might be expensive. Moreover¡ the method with global efficiency is not reachable. The methods based on second approach generally are three types: reviewing and correction by volunteers¡ using the general rules of geospatial data and extracting information from VGI to validate itself. The methods for reviewing and correction by volunteers argue that crowdsourcing data converges on the truth if people have the opportunity to review and correct errors. But features that many people have an interest will be more accurate than features that are of interest only to a few. In second methods¡ general rules for geospatial data like logical consistency used for VGI validation. The last method tries to use information from VGI to validate itself. For example a method determines information trust just by using the data history like production¡ change or removing of data to validate VGI just with itself.
    In this article we try to present a method for validation of quality of VGI without using authorized data and by using of other VGI. In this article¡ calculable characteristics of OpenStreetMap (OSM) data¡ as one of the important sources of VGI¡ are used for validate the quality of its VGI data. Specially¡ we analysis the accuracy of classification of OSM roads network as part of semantic data of VGI. Roads¡ such as motorway and residential¡ have different design and characteristics according to their function in street network. It is tried to classify these roads by finding and learning these distinguish characteristics of each road class.
    Machine learning models with decision tree and neural network algorithms are used to learn roads characteristics from OSM street network data. Decision tree and artificial neural network with multilayer perceptron are usable for data that contains errors and we have to assume VGI data always have errors. These two Statistical analysis¡ precision and recall are used for assessing final models. In accordance with the result¡ Decision tree model have 84.1% weighted average accuracy and represent a suitable model for this method. These methods are based on extractable information from VGI and could be used for any street network to classify the streets.
    Keywords: OpenStreetMap, Data Quality, Volunteer Geospatial Information, Street Network Analysis, Learning Machine
  • M. A. Borooman, A. Safari, Y. Allahtavakoli, A. Bahroudi Pages 13-25
    Gravity gradient tensor is a matrix containing the second order derivatives of the Earth’s gravity field, which has numerous applications in geodesy and geophysics. To date, much effort has been done for estimating gravity gradient tensor with reasonable accuracy. This quantity can be estimated via using various methods, and one of these methods is applying finite-difference method to gravity observations. Finite-difference method can estimate gravity gradient tensor directly by using the mathematical concept of gradient, regardless of extra assumptions. This ability of finite-difference method, from theoretical point of view, provides the possibility of accurate estimation of gravity gradient tensor without considering additional assumptions to the problem. This study tends to introduce and evaluate Finite difference method for estimating the gradient tensor and present formulae for determining gravity gradient tensor from land-based gravity observations. In this paper, the proposed equations are numerically tested by means of using a global gravity model of the earth. Global gravity model of the earth (EGM 2008) is a geopotential model of the earth consisting of spherical harmonic coefficients up to degree 2190 and order 2159. There are numerous uses for these high degree potential coefficient models. One of these uses is modeling and estimating gravity gradient tensor.
    Finally, gravity gradient tensor is estimated by the proposed method for 6350 gravity stations located in Costal Fars region in a northern part of the Persian Gulf, between the latitudes from 26.5 N to 27.27 N and longitudes from 53.41 E to 55.58 E. The target area is about 10000 square kilometers. About 8500 square kilometers of the study area is located in moderate mountainous regions, and about 1500 square kilometers is located in flat coastal areas. The altitudinal distribution and spatial distribution of gravity in study area are shown in figure 1 and 2 respectively.
    Numerical experiments of this study demonstrate the ability of this method in gravity gradient tensor estimation with acceptable accuracy. For example, numerical experiments showed that the proposed method can estimate diagonal components of gravity gradient tensor (second order derivatives of the Earth’s gravity field in east, north and vertical directions) with the accuracy values of 12.46, 34.49 and 454.82 Eotvos respectively. The spatial distribution of the gravity gradient tensor components obtained from finite difference method in study area are shown in figure 3.
    Finally, according to the theoretical concepts discussed in this paper, It can be said because the finite difference method using from derivative and difference concepts directly for estimating gravity gradient tensor, it is expected that this method provide accurate estimation of gravity gradient tensor, As this is happen in the simulation conducted. However the accuracy of this method is very dependent on distances between sampling stations and by reducing distance between the stations, the accuracy of proposed method will be increased. The numerical results of this study also showed that the proposed method can provide accurately estimate of gravity gradient tensor components In some stations surrounded by suitable spatial distribution of gravitational observations.
    Keywords: Gravity Gradient Tensor, Finite, Difference Method, Earth Gravity Field, Land, Based Gravity Data
  • M. Esfahani, A. Mohammadzade Pages 27-42
    Photogrammetry and Remote sensing researchers pay attention to urban phenomenon detection; such as building and tree from aerial and satellite image with LiDAR data in the recent decades. Different classification and detection methods in order to use in urban area faced with complex problems, such as small trees and buildings, unfavorable boundary crown, roof covered with vegetation , buildings that are surrounded by trees and vegetation located in the shadow. In this research to improve the problem mentioned above, in the first step the following works has been carried out, preprocessing LiDAR data That in the preprocessing step, noise and outliers of the LiDAR data detected by using Frist and Last Range images and then to calculate the proper height for them used the weighted averaging method than the distance, then at this step generation digital terrain model and normalized digital surface model that for generating digital terrain model the method local base filtering is used that has two stage 1- select training data 2- applying filter with using training data that all part of the first step is implemented in the MatLab 7.12 . In the second step geometrical, spectral, textural and combined features have been produced according to the mentioned problems that to produce these features used from normalized digital surface model and first range of LiDAR data. Then these features by using the SVM_GA algorithm for detecting suitable they are used. SVM_GA algorithm implementation is done in tow rounds, that in first round same features are evaluated with each other and in the second round selected features in first round are evaluated with each other. Finally after three times running SVM_GA algorithm 27 features are selected. From the production features vegetation index combined with shadow area and geometrical features produced from LiDAR data have crucial role for trees detection. In the third step support vector machine classification has been used in order to trees and building detection in the object-based and pixel-based levels with using suitable and selected features. Training data for pixel level classification with using SVM classifier are selected semi-automatically in four classes include road, tree, building and vegetation. For this purpose, the following features are used nDSM, Laplacian and Combined_IndexNDVI . for object-based classification with SVM classifier we use training data that they were selected manually and for segmentation of image we use multi resolution segmentation in eCognition 8 software by using the flowing layers, green, red, NIR and Frist range that segmentation parameters are set manually. In order to improving the results of each levels of classification in fourth step we improved with post-processing methods that in this step for both level of classification we have used erosion and dilation morphological filter with different size, so that in object-based level classification, segments of buildings divided into two group: high and low buildings so that for each group of buildings employed morphological filter with different size Separately and then for results of object-based classification we have used a conceptual process for improving object-based level classification. Due to high spectral similarity between two group buildings-roads and vegetation-trees and miss-classification, accuracy, correctness and quality of results have reduced. To improving results have been used neighborly relations and conceptual thresholds in five steps.In fifth step due to the power of each level classification, we try to improve buildings and trees detection with use to fusion object-based and pixel-based classification results, that in this step for building and tree classes have been used different fusion algorithms. The main basis of fusion for building’s class voting results of two level classification and for trees class height pixels in segments and again voting. Results show that fusion of pixel-based and object-based classification improve buildings and trees detection accuracy especially in small and low objects of trees and building also improve crown detection. Fusion the classification results for trees class has more improvement compared to buildings class. Object-based classification level due to using the segments and similarity between trees and vegetation leads to pixels of vegetation not be detected, that in the segments of trees and with use to fusion method for both level of classification, can be detected low and little trees. Results of this method with regards to improving detection in boundary of trees and buildings have positive effect for object extraction. Also suggested method can be detected building with low height and area. After post-processing for detected trees and building classes specified pixel level can be show more capable than object level in detect small object. Finally, results of detection and each levels classification evaluated with reference data, That buildings and trees detection in object level have correctness 0.961 and 0.65 and they have overall accuracy 0.97 and 0.943 respectively, and in pixel-based level they have correctness equal to 0.953 and 0.632 and overall accuracy 0.961 and 0.94 respectively. After the fusion they have correctness 0.971 and 0.718 and they have overall accuracy 0.975 and 0.957 respectively.
    Keywords: Urban Buildings, Trees Detection, Classification, Fusion, Object, Base, Pixel, Base, Aerial Imagery, LiDAR Data
  • Y. Allahtavakoli, A. Safari, A. A. Ardalan, A. Bahrodi Pages 43-62
    This paper aims at combining the geodetic gravimetric inverse problem with the geophysical one to do ‘local gravity field modeling’ and ‘topographic mass-density anomaly determination’ simultaneously. For this purpose, it recalls the basic theories of the two types of gravimetric inverse problems and finds certain relations between their corresponding unknown parameters. The current research proposes a method, which is based on certain anharmonic Radial Basis Functions (RBFs) as well as Generalized Tikhonov regularization method, for modeling the gravity field at the Earth’s surface. Then, it shows that the proposed method can provide certain formulae, based on a combination of the Newton’s IE and the Poisson’s PDE, for extracting the mass-density anomalies from land-based gravity data. Finally, the proposed method is numerically examined by a case study at which the proposed method has been applied to 6350 gravity stations scattered over the coastal Fars of Iran in the north of the Persian Gulf. The case study numerically indicates the possibility of the main idea of the research for simultaneously solving the two types of gravimetry inverse problem.
    1.Gravimetric inverse problems: From point of view of mathematics, a geophysical gravimetric inverse problem can be defined by the following equation:where the mass-density anomaly is considered as the unknown and defined by subtracting the mass-density functional from its prior knowledge ; is the Newton integral operator; is the gravity anomaly yielded by applying the gravity reduction to land-based gravity observations. In a local coordinate system which z-axis is a tangent to the plumb line and outwards from the Earth’s surface, the above geophysical gravimetric inverse problem (i.e. Eq.) can be re-written as follows:The above equation is a type of the Newton’s IE at which is the position vector, is the topographic space of the region of the interest, and is the differential volume element at the position .
    On the other hand, the above definition is so similar to the definition of the gravimetric inverse problem associated with the Earth’s gravity field modeling in physical Geodesy. The Earth’s gravity field modeling in physical Geodesy is often involved with various types of gravimetric inverse problem such as downward continuation, at which the incremental gravity observations such as gravity anomaly are also placed at the right side of corresponding equations and the unknowns are also placed under the mathematical operators. Land-based, airborne, ship borne and satellite gravity data, as well as astronomical, leveling and tidal observations are counted as common observations of physical Geodesy, which are usually used for modeling the Earth’s gravity field. In this way, the unknowns are usually constituted by coefficients of certain algebraic expansions which enable us to model the Earth’s gravity potential as well as its derivatives such as the gravity vector and gravity gradient tensor. These algebraic expansions are commonly linear and obtained by discretizing linear integral equations achieved by Geodetic Boundary Value Problems (GBVP). Furthermore, these expansions can also be obtained by certain basis functions such as spherical harmonic functions, radial basis functions (RBFs), Slepian and wavelet basis functions. Hence, the current research considers the following equation as the general form of the gravimetric inverse problem associated with the subject of the gravity field modeling in physical Geodesy:
    where is an operator which depends on either the corresponding integral equations or the expansions of the basis functions; and the vector comprises of the unknowns which model the Earth’s gravity potential and its derivatives. Now, if the Earth’ gravity potential is able to be expanded by a set of basis functions such as , the above general form of the gravimetric inverse problem (i.e. Eq.) will have the following form in the above-mentioned local coordinate system.where and .
    This paper comprises of an idea of correlating the concepts of the abovementioned gravimetric inverse problems with each other. Hence, first it analyzes each one of the abovementioned gravimetric inverse problems separately and then combines them with each other by finding mathematical relations between their unknowns. For this purpose, the concept of the Poisson’s PDE is used for converting the unknowns of Eq. (i.e. ) into the unknown of the Newton’s IE (i.e. the mass-density anomaly ). Indeed, the theories discussed by the paper will also examined by certain numerical experiments in a real case study based on 6350 gravity stations scattered over a part of the north coasts of the Persian Gulf.
    Keywords: Gravity, Mass Density, Gravity Modeling, Inverse Problem
  • M. Omati, M. R. Sahebi Pages 63-78
    Change detection using remotely sensed data has been used in many applications, such as the detection of dynamic changes in land cover and land use, the monitoring of forestland and agricultural land, the assessment of damage from natural disasters, and the study of urban environments. Despite the numerous studies devoted to multispectral and hyperspectral imagery, however, the use of optical imaging sensors is limited to weather conditions. Synthetic aperture radar (SAR) sensors can obtain daylight, cloud coverage, and weather-independent images. Their backscattered signals are also sensitive to the form, orientation, homogeneity, and surface conditions of a target. SAR imagery can therefore serve as a useful tool for detecting land cover and land use changes. The development of SAR techniques has given rise to Polarimetric SAR (PolSAR) systems, which measure four linear polarization channels (i.e., HH, HV, VH, VV) and the phase differences among them. SAR polarimetry with the functionality of identifying different scattering mechanisms can provide more significant information than single-channel imagery. From the perspective of image analysis unit, Change detection techniques are classified into two categories, namely, pixel- and object-based approaches. Pixel-based methods rely only on the information derived from individual pixels and do not consider the spatial relationships among these pixels. Whereas, the values of neighboring pixels in an image are strongly correlated, and the probability of change occurrence in adjacent regions is more than distinct points. The use of spatial features also effectively reduces the speckle effect in PolSAR images. Given these considerations, object-based approach has been widely used for PolSAR image analysis over recent years. This paper proposes a novel image segmentation algorithm for improving the accuracy of land cover change detection. This method consists of the following three steps: 1) segmenting two PolSAR images by new segmentation technique, namely region based improved watershed; 2) selecting the optimal differential of polarimetric features based on the Genetic Algorithm (GA) and Jefferies-Matusita (JM) distance criteria; and 3) The binary classification of image objects using the differential of mean pixel values of the corresponding image objects. Despite the development of various region-based segmentation methods, watershed segmentation is appropriate for the segmentation of high resolution images based on the many advantages of this morphological algorithm. These advantages include inherent simplicity, high speed implementation, the creation of separated regions in low contrast images, and the provision of closed connected regions. Common watershed segmentation approaches, such as distance transform and the gradient method, cause over-segmentation problem given the noise or local irregularities present in a gradient image. Unlike the direct application of the watershed algorithm, using a marker-controlled approach that involves the incorporation of additional knowledge into a segmentation procedure can limit the number of segmented regions. In this method, the flooding procedure begins from a previously defined set of markers. Markers, as connected components belonging to specific areas of an image, can be defined on the basis of a set of descriptors, such as gray level value, shape, location and texture. Compared with conventional watershed and multi-resolution segmentation methods, the improved watershed reduces the speckle effect in PolSAR images and avoids the over- segmentation problem. The results of proposed change detection method on Uninhabited Aerial Vehicle Synthetic Aperture (UAVSAR) full polarimetric images achieve the overall accuracy of 92.40% and the 0.85 kappa coefficient.
    Keywords: Change Detection, Polarimetric SAR (PolSAR) Images, Object, Based Analysis, Improved Watershed Algorithm
  • M. Ahangarcani, M. Farnaghi, M. R. Shirzadi Pages 79-98
    Detection of pathogenic factors, identify the spatial accumulation of disease cases and finding its distribution pattern are of urgent need in the field of public health and disease management and control. Leptospirosis is a zoonotic disease which occurs worldwide but is most common in tropical and subtropical areas with high rainfall. Wet and mild weather conditions in the northern provinces of Iran have put these areas at high risk for Leptospirosis incidence. The main objectives of this study are to investigate the annual pattern of Leptospirosis distribution, identify the spatial and spatio-temporal clusters of the disease and generate the annual predictive map of spatio-temporal distribution of Leptospirosis at the district level in the Northern provinces of Iran. In this study, Leptospirosis incidences, census data and topographical and climate factors have been used to generate the annual predictive map of spatio-temporal distribution of Leptospirosis. The Leptospirosis incidences are continuously recorded by the Center for Disease Control and Prevention of Ministry of Health of Iran. The population census count estimates for period 2009-2014 were obtained from the Statistical Center of Iran. Topographical data were used to generate the altitude, slope and aspect maps. Climate data such as average temperature, average humidity, annual rainfall and number of freezing days were used to model other affecting parameters. Global clustering methods including Moran’s I and general G indices were applied to investigate the existence of spatial autocorrelation between cases of Leptospirosis and also analyze the annual spatial distribution of the existing patterns. Results of both Moran’s I and general G indices indicated meaningful persistent spatial autocorrelation between Leptospirosis cases and highly clustered distribution of Leptospirosis. Additionally, presence of spatial clusters of Leptospirosis and detection of high risk areas of disease were investigated using the local Moran’s I and local G indices. The results of the local Moran’s I and local G indices identified significant spatial clusters of Leptospirosis cases located in central, north-eastern and western parts of Guilan, Mazandaran and Golestan provinces, respectively. Geographically weighted regression (GWR) and multilayer perceptron neural network (MLP) models were used to generate the annual predictive map of spatio-temporal distribution of Leptospirosis and modelling the relation between the distribution of Leptospirosis cases with topographical and climate factors. Performance of GWR and MLP models were compared using Kappa coefficient, RMSE, MAPE and R2 measures. The evaluation results showed that the MLP model was able to predict the incidence rate of Leptospirosis in 2014 for each district with acceptable accuracy. MLP was able to model the relationship between Leptospirosis incidence and factors better than GWR. Additionally, results of both GWR and MLP models showed that average humidity and annual rainfall are most important affecting factors on Leptospirosis incidence in the Northern provinces of Iran, respectively. Such predictive maps can be used to provide essential guidelines for planning of effective control strategies and identification of high risk areas of Leptospirosis which should receive more preventive measures from policy makers and healthcare authorities.
    Keywords: Leptospirosis, Geographical Information Systems, Spatio, Temporal Distribution, Geographically Weighted Regression, Neural Network, Iran
  • S. Nadi, Y. Ghiasi, Sh. Hadavand Pages 99-118
    Recent developments in open source digital elevation models (DEMs) increase their applications in different disciplines. The lack of proper information about the accuracy of these models in regions with different topographic characteristics, make their analysis unreliable. Most estimates of accuracy, e.g., mean, standard deviation and RMSE depend on an implicit assumption that the errors comprise a random sample from a normal distribution. But, analytic data often depart from that assumption. Therefore, this paper first study statistical characteristics of SRTM and GDEM open source global DEMs. Then we define a confidence interval for RMSE which make it possible to investigate the appropriateness of reference points. Afterwards, measures for accuracy assessment of DEMs based on robust methods in L1 norm and Huber's method as well as conventional methods in L2 norm are discussed. In order to consider different topographic characteristics we compare SRTM and GDEM based on the defined measures in Urban, Flat, Hill and Mountain regions. Finally, we studied how the errors propagate into slope and aspect maps. It was found that although GDEM has better resolution, but SRTM performs better based on the defined measures in all the regions. It is also illustrated that the accuracy of both models in urban are better than other regions.
    Keywords: Open Source DEM, Accuracy Assessment, Robust Statistics, Error Propagation, Confidence Interval
  • E. Kiana, S. Homayouni, M. A. Sharifi, M. R. Faridrohani Pages 119-130
    Nowadays, Earth observation (EO) technology became an indispensable tool to help environmental monitoring, as well as their changes, for natural resources management, urban planning and development, water management and land use planning. In particular, radar EOs, unlike the optical ones, can be collected regardless of illumination and weather conditions. Multitemporal polarimetric synthetic aperture radar (PolSAR) images are useful source of information for detection and mapping the environmental changes, especially in wide areas, during the day and night and all weather conditions. Change detection methods can identify the change or no change conditions in land covers using the time series observations. In this paper a method is proposed for change detection in SAR remote sensing images. This method is based on the Change Point Analysis. The cumulative frequency of difference image, which contains the environmental changes, normally follows a specific class of statistical distribution. Gaussian mixture model is one of the most suitable models for Change Point Analysis. This model can efficiently estimate the parameters of mixture distribution. The intersection point of two distributions is a change point, which can be seen as a threshold. This threshold is then used to separate the change and no change classes. The proposed method is implemented and analyzed using three SAR data sets. The analytical evaluations of the final change maps from two of these data sets with reference data had the Kappa coefficients of 90% and 96% respectively. The other data set contained the multitemporal PolSAR images and had been acquired over an agricultural area. The changes in these images were enough reliable to be connected to the agricultural activities, such as crop growing stages and harvesting, based on an available crop map. Finally, the method was evaluated against the Otsu method, as one of the best threshold estimation methods, and the results showed the superiority of the proposed method, e.g. 2% better in term of kappa coefficient. . As a result, the proposed method, can be efficiently employed for land cover change detection and monition in natural resources management.
    Keywords: Unsupervised Change Detection, SAR Images, Change Point Analysis, Gaussian Mixture Model, Threshold Estimation
  • M. Rostamnia, M. Akhoondzadeh Hanzaei Pages 131-144
    Forest dieback is a complex and important phenomenon that happened in the world’s most oak forests nearly a century ago. In recent years, this has occurred in ZAGROS oak forests due to successive droughts. ZAGROS forest region, with an area of about 6 million hectares including provineces West Azerbaijan, Kurdistan, Kermanshah, Ilam, Lorestan, Khuzestan, Fars, Esfahan, Chaharmahal and Bakhtiari, Kohgiluyeh and Boyer-Ahmad, Hamedan is exposed to several threats. Forests of ILAM province are part of ZAGROS forests, which are located in the west part of the mountain chain. The phenomenon of forest dieback has been seen in this area during recent years. In this study the dieback of forests in ILAM province area has been studied using Landsat satellite imagery and Google earth images, meteorological data and information of the amount of dust in province atmosphere. This study has conducted over a period of 15 years. Also we have to mention that the dust data have been extracted from 550 nm band Modis satellites images. In order to detect vegetation and forest area of ILAM, vegetation indices have been used in Landsat satellite images. In this study, the five vegetation indices, NDVI (Normalized Difference Vegetation Index), EVI (Enhanced Vegetation Index), RVI (Relative Vigor Index), SAVI (Soil Adjusted Vegetation Index), ARVI (Atmospherically Resistant Vegetation Index) have been used. Also to determine the appropriate vegetation index for detecting the changes in vegetation coverage areas, these indicators have been compared to each other. In order to compare these indicators, Google earth imagery has been used. Also with the use of Google earth images, the vegetation coverage of 200 pixels has been estimated in Landsat satellite images. By using these 200 pixels we compared the vegetation indices and carefully calculated the accuracy of each one. In this study, the generalized correlation of the indices with the vegetation coverage changes is respectively: 0.939, 0.953, 0.945, 0.914, 0.925 and according to these results, EVI vegetation index is used as the preferable indicator for vegetation change detection index in ILAMs forests. Landsat Images were obtained for each year and the vegetation changes amount was calculated towards the previous year for each year. Using the rainfall data index SPI (Standardized Precipitation Index), in these periods of one year, nine-month and three-month. The correlation between EVI index annual changes examined using SPI index the results were respectively, 0.43, 0.76, 0.09. In this paper, we used PCA (Principal Component Analysis) method to detecte the changed areas in forest vegetation coverage. The rate of changes in pixels amount in each year is calculated respect to year 1379. Also in order to separate the forest area changes from others, the behaviour of each pixel is studied during a period of 15 years. Also with the introduction of two patterns and determination of correlation between the patterns change behaviour, ILAMs vegetation coverage changed areas were isolated from other changes, and by deleting the non–changed pixels, the average annual change was obtained for ILAM province forest. Average entered dust was calculated for each year, using meteorological and Modis satellite images data. Also by applying a two–parameter linear regression, the combined impact of two factors, rainfall and dust was determined. The study implied that precipitation is the most effective parameter on dieback. The influence of two factors, rainfall and dust are respectively 62% and 38%.
    Keywords: Trees Drought, Vegetation Index, SPI, PCA, Ilam
  • N. Tahani, A. A. Alesheikh, M. Farnaghi Pages 145-160
    Repre­­sentation of the user's environment along with functionalities like route finding and exploring point of interests (POIs), is usually provided in the form of 2D maps which limit the user to a top-down view of the nearby environment. Considering the 3D form of real-world objects, the presentation of them as 3D city model plays an important role in visualization and location-based functionalities. In the visualization part, 3D city model helps users to identify real-world objects by making a virtual environment and regarding the geometric attributes of them. Furthermore, 3D city model is utilized in variety of applications such as visibility analysis, energy demand estimation, urban facility management, 3D cadastre and indoor navigation. Recent improvements in the hardware and software specifications of smart phones along with the utilization of powerful CPUs and GPUs provided the developers with the ability of exploiting mobile devices as platforms for different location-based usages. Virtual Reality (VR) is a technology that simulates objects of the real world in a virtual environment and enables users to interact with those objects. 3D city models that contain elements like buildings, vegetation areas, roads and terrain, are the examples of VR applications in smart phones and tablets. This paper describes a mobile-based VR application to represent and explore campus objects exists in the faculty of Geodesy and Geomatics Engineering of K.N.Toosi University of technology. Generally, OpenGL ES API was used to render 2D and 3D graphic objects in mobile devices and hardware accelerated by using GPU. In spite of that, 3D city model elements like buildings are complicated objects, which produce in desktop graphic softwares such as Google Sketchup, City Engine, 3ds MAX. A variety of game engines and graphic libraries are developed to load and render graphic models based on the low-level structure of OpenGL ES. These libraries provide functionalities like loading and creating objects, supporting different graphic formats, collision detection and so on. Our mobile application is developed based on Libgdx cross-platform game engine. 3D models of campus objects are produced in Google Sketchup software form 2D footprint of buildings in 1:2000 aerial map. In the next step, models were textured and their coordinates and non-spatial attributes were stored in a KML file format. The serialization process was performed in the mobile device to construct and render models in graphic environment. In this study, we have implemented two types of representation and interaction with 3D city model: 1) Birds-eye view and 2) First-person view. In the Bird's eye view, user is elevated above the surface and control buttons and touch screen events handle user's interactions. On the other hand, the first-person view is mimic of Augmented Reality (AR) view and user can walk in the environment and explore nearby objects. In this approach, the position and orientation of the mobile device have determined by sensory data and the low-pass filter has been utilized to decrease the noise of the data. In order to represent non-spatial information of objects, we have utilized frustum-culling concept to determine the selected model by the user and retrieved corresponding attributes. For that, we created a minimum bounding box (MMB) of each model and intersected pointing vector of the observer with MMBs to select nearest models to the user. The mixed interactive approach of our research enables users to explore their nearby urban environment in two different view and gets related attribute information of models directly.
    Keywords: Virtual Reality, 3D City Model, Libgdx Graphic Library, Bird's Eye View, First Person View