# Import arcpy module
import arcpy
from arcpy import env
from arcpy.sa import *

if arcpy.CheckExtension("Spatial") == "Available":
    arcpy.CheckOutExtension("Spatial")
else:
    raise LicenseError
env.workspace = "F:\\GIS_OUTPUTS"
env.overwriteOutput = 1

env.OverwriteOutput=1

# Local variables:

#Inputs
gis_osm_roads_free_1_shp = "F:\\GIS_DATA\\Ireland_uses\\gis_osm_roads_free_1.shp"
Census2011_Admin_Counties_generalised20m_shp = "F:\\GIS_DATA\\AdminCounties\\Census2011_Admin_Counties_generalised20m.shp"
CLC18_IE_ITM_shp = "F:\\GIS_DATA\\ITMData\\CLC18_IE_ITM.shp"

#Outputs
Carlowsel_shp = "F:\\GIS_OUTPUTS\\Carlowsel_1.shp"
roadclip_shp = "F:\GIS_OUTPUTS\\roadclip.shp"
Census2011_Admin_Counties_ge1 = "F:\\GIS_OUTPUTS\\Census2011_Admin_Counties_generalised20m_gel.shp"
primselect_shp = "F:\\GIS_OUTPUTS\\primselect.shp"
primselect_Buffer = "F:\\GIS_OUTPUTS\\primselect_Buffer.shp"
Output_direction_raster = ""
secsel_shp = "F:\\GIS_OUTPUTS\\secsel.shp"
secsel_Buffer = "F:\\GIS_OUTPUTS\\secsel_Buffer.shp"
Output_direction_raster__2_ = ""
MW_shp = "F:\\GIS_OUTPUTS\\MW.shp"
MW_Buffer = "F:\\GIS_OUTPUTS\\MW_Buffer.shp"
Output_direction_raster__3_ = ""
HospitalsITM_shp = "F:\\GIS_DATA\\ITMData\\HospitalsITM.shp"
HospitalsITM_Buffer = "F:\\GIS_OUTPUTS\\HospitalsITM_Buffer.shp"
Output_direction_raster__4_ = ""
pdist = "F:\\GIS_OUTPUTS\\pdist.tif"
pclass = "F:\\GIS_OUTPUTS\\pclass"
secdist = "F:\\GIS_OUTPUTS\\secdist.tif"
seclass = "F:\\GIS_OUTPUTS\\seclass"
mwdist = "F:\\GIS_OUTPUTS\\mwdist.tif"
mwclass = "F:\\GIS_OUTPUTS\\mwclass"
hospdist = "F:\\GIS_OUTPUTS\\hospdist.tif"
Reclass_hosp3 = "F:\\GIS_OUTPUTS\\Reclass_hosp3"
CLC18_IE_ITM_shp = "F:\\GIS_DATA\\ITMData\\CLC18_IE_ITM.shp"
LUcarlow = "F:\\GIS_OUTPUTS\\LUcarlow.shp"
LUrast = "F:\\GIS_OUTPUTS\\LUrast.tif"
luclass = "F:\\GIS_OUTPUTS\\LUclass"

# Process: Select County
arcpy.Select_analysis(Census2011_Admin_Counties_generalised20m_shp, Carlowsel_shp, "\"COUNTYNAME\" = 'Carlow County'")

# Process: Select Road type: Primary Roads
arcpy.Select_analysis(roadclip_shp, primselect_shp, "fclass = 'primary'")

# Process: Buffer Primary Roads
arcpy.Buffer_analysis(primselect_shp, primselect_Buffer, "100 Meters", "FULL", "ROUND", "ALL", "", "PLANAR")

# Process: Convert Primary Roads to Raster
pdist = EucDistance("primselect_Buffer.shp", 500,"pdist.tif")

# Process: Select Road type: Secondary Roads
arcpy.Select_analysis(roadclip_shp, secsel_shp, "fclass = 'secondary'")

# Process: Buffer Secondary Roads
arcpy.Buffer_analysis(secsel_shp, secsel_Buffer, "50 Meters", "FULL", "ROUND", "NONE", "", "PLANAR")

# Process: Convert Secondary Roads to Raster
secdist = EucDistance("secsel_Buffer.shp", 500,"secdist.tiff")

# Process: Select Road type: Motorways
arcpy.Select_analysis(roadclip_shp, MW_shp, "fclass = 'motorway'")

# Process: Buffer Motorways
arcpy.Buffer_analysis(MW_shp, MW_Buffer, "1.5 Kilometers", "FULL", "ROUND", "NONE", "", "PLANAR")

# Process: Convert Motorways to Raster
arcpy.gp.EucDistance_sa(MW_Buffer, mwdist, "", "1.78079161064443E-03", Output_direction_raster__3_, "PLANAR")

# Process: Buffer Hospitals
arcpy.Buffer_analysis(HospitalsITM_shp, HospitalsITM_Buffer, "10 Kilometers", "FULL", "ROUND", "ALL", "", "PLANAR")

# Process: Convert Hospitals to Raster
arcpy.gp.EucDistance_sa(HospitalsITM_Buffer, hospdist, "", "959.985741", Output_direction_raster__4_, "PLANAR")

# Process: Reclassify Primary Roads
arcpy.gp.Reclassify_sa(pdist, "VALUE", "0 0.040140 9;0.040140 0.100362 4;0.100362 0.226174 1", pclass, "DATA")

# Process: Reclassify Secondary Roads
arcpy.gp.Reclassify_sa(secdist, "VALUE", "0 0.040140 9;0.040140 0.100362 4;0.100362 0.226174 1", seclass, "DATA")

# Process: Reclassify Motorways
arcpy.gp.Reclassify_sa(mwdist, "VALUE", "0 0.011320 9;0.011320 0.030922 7;0.030922 0.052019 5;0.052019 0.075040 3;0.075040 0.113310 1", mwclass, "DATA")

# Process: Reclassify Hospitals
arcpy.gp.Reclassify_sa(hospdist, "VALUE", "0 10857.769531 1;10857.769531 25175.388672 3;25175.388672 43889.371094 5;43889.371094 69399.085938 7;69399.085938 118453.414063 9", Reclass_hosp3, "DATA")

# Process: Convert clipped Land Use to Raster
arcpy.FeatureToRaster_conversion("LUcarlow.shp", "Class_Desc", LUrast)

# Process: Reclassify Land Use
arcpy.gp.Reclassify_sa(LUrast, "CLASS_DESC", "'Discontinuous urban fabric' 5;'Industrial or commercial units' 4;'Road and rail networks and associated land' 2;'Mineral extraction sites' 0;'Dump sites' 0;'Sport and leisure facilities' 4;'Non-irrigated arable land' 8;Pastures 8;'Complex cultivation patterns' 5;'Land principally occupied by agriculture, with significant areas of natural vegetation' 6;'Broad-leaved forest' 6;'Coniferous forest' 6;'Mixed forest' 6;'Moors and heathland' 0;'Transitional woodland-shrub' 3;'Burnt areas' 3;'Peat bogs' 0;'Water courses' 0", luclass, "DATA")

# Process: Calculate best Site Selection
outRaster = Raster("luclass") + Raster("mwclass") + Raster("pclass") + Raster("seclass") + Raster("Reclass_hosp3")
outRaster.save("F:\\GIS_OUTPUTS\\Final")

#---------------------------------------------------------------------------------------------