ایجاد پلیگون های تیسن در GIS

انتخاب محدوده تحت پوشش باران سنج ها (ایستگاه های باران سنجی که لایه نقطه ای هستند) کار ساده ای نیست. مقدار بارانی که توسط باران سنج اندازه گیری می شود نمایانگر بارندگی در سطح وسیعی از منطقه خود می باشد. اگر ما تعدادی ایستگاه باران سنجی (لایه نقطه ای) داشته باشیم و قصد ما آن باشد که محدوده ای را که می بایست توسط هر ایستگاه باران سنج پوشش داده شود تعیین کنیم، پلیگون های تیسن به کمک ما می آیند.

پلیگون های تیسن ، يک منطقه را بر اساس نقاطي که در آن واقع شده است، به پليگون هاي مجاور هم تقسيم مي کند به طوري که در داخل هر پليگون يک نقطه قرار مي گيرد. محدوده هاي مربوط به هر پليگون، خصوصيات نزديکي به نقطه درون آن پليگون دارند. مشخص نمودن مناطق تحت تاثير يك بانك و مناطق تحت تاثير يك ايستگاه باران سنجي از جمله كاربردهاي اين تحليل مي باشد. روش ايجاد پليگون هاي تيسن به اين نحو است که ابتدا کليه ایستگاه های باران سنجی (لایه نقطه ای) تبديل به شبکه نامنظم مثلثي مي گردد. عمود منصف اضلاع اين مثلث ها ترسيم که از بهم پیوستن آنها به يکديگر پلیگون های تیسن ايجاد مي گردد.

ایجاد پلیگون های تیسن در GIS

برای ایجاد پلیگون های تیسن کافیست در نرم افزار ArcGIS و بخش ArcToolbox ابزار Create Thiessen Polygons را از آدرس زیر فراخوانی نمایید:

حتما بخوانید تخمین بارش ماهواره ای

ArcToolbox > Analysis tools > Proximity > Create Thiessen Polygons

پلیگون های تیسن در GIS

Input Features: لایه نقطه ای (Point) خود را که میخواهیم براساس آن پلیگون های تیسن را ایجاد کنیم، در این قسمت معرفی می کنیم.

Output Feature Class : در این بخش، نام و مکان ذخیره را برای لایه پلیگون های تیسن نهایی انتخاب می کنیم.

Output Fields : دراین قسمت مشخص می کنیم که پلیگونی که به ازای هر نقطه ایجاد می شود، حاوی تمامی اطلاعات توصیفی آن نقطه باشد، و یا صرفاً ستون FID را از لایه نقطه ای دریافت کند.

تصویر خروجی به صورت زیر خواهد بود:

پلیگون های تیسن در GIS

همچنین می توانید با کدهای پایتون نیز پلیگون های تیسن در GIS را ایجاد کنید:

import arcpy
arcpy.env.workspace = "C:/data/data.gdb"
arcpy.CreateThiessenPolygons_analysis("schools", "c:/output/output.gdb/thiessen1", "ALL")

که توضیحات فوق به زبان ساده تر:

# Name: CreateThiessenPolygons_Example2.py
# Description: Creates Thiessen polygons
 
# Import system modules
#
import arcpy
from arcpy import env
 
# Set environment settings
#
env.workspace = "C:/data/data.gdb"
 
# Set local variables
#
inFeatures = "schools"
outFeatureClass = "c:/output/output.gdb/thiessen1"
outFields = "ALL"
 
# Execute CreateThiessenPolygons
#
arcpy.CreateThiessenPolygons_analysis(inFeatures, outFeatureClass, outFields)