20190911
17511400
1.0
FALSE
MosaicToNewRaster \\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\serap_sleuth\serap_urb2100;\\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\app_sleuth\app_urb2100;\\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\gcpo_sleuth\gcpo_urb2100;\\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\gcp_sleuth\gcp_urb2100 \\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\se_sleuth.gdb sesleuth2100 PROJCS['NAD_1983_Albers',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5],PARAMETER['Latitude_Of_Origin',23.0],UNIT['Meter',1.0]] 16_BIT_SIGNED 60 1 LAST FIRST
CopyMultiple "'C:\Users\558351\Desktop\REPI\REPI_Urbanization\SLEUTH Mosaic Geodatabase.gdb\se_sleuth.gdb\sesleuth2100' RasterDataset" "C:\Users\558351\Desktop\REPI\REPI_Urbanization\New File Geodatabase.gdb" sesleuth2100 "sesleuth2100 RasterDataset sesleuth2100 #"
sesleuth2100
-567291.227000
1892168.773000
51212.078000
2396552.078000
1
file://\\IGSKFAEBFSSDMS\SciDMS\SC00E6E\D2\Working\SLEUTH\20190911\se_sleuth.gdb
Local Area Network
002
Projected
GCS_North_American_1983
Linear Unit: Meter (1.000000)
NAD_1983_Albers
<ProjectedCoordinateSystem xsi:type='typens:ProjectedCoordinateSystem' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/10.7'><WKT>PROJCS["NAD_1983_Albers",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",23.0],UNIT["Meter",1.0]]</WKT><XOrigin>-16901100</XOrigin><YOrigin>-6972200</YOrigin><XYScale>266467840.99085236</XYScale><ZOrigin>-100000</ZOrigin><ZScale>10000</ZScale><MOrigin>-100000</MOrigin><MScale>10000</MScale><XYTolerance>0.001</XYTolerance><ZTolerance>0.001</ZTolerance><MTolerance>0.001</MTolerance><HighPrecision>true</HighPrecision></ProjectedCoordinateSystem>
20190912
16014600
20190912
16014600
VAT_sesleuth2100
Table
16
OBJECTID
OBJECTID
OID
4
0
0
Internal feature number.
Esri
Sequential unique whole numbers that are automatically generated.
VALUE
VALUE
Integer
4
0
0
COUNT
COUNT
Integer
4
0
0
311B40E3-E52E-4656-9272-B8D3DFF13078
dataset
Curtis Belyea
Biodiversity and Spatial Information Center, Dept. of Applied Ecology, NCSU
Biologist, GIS and Database Specialist
919-513-7337
919-515-4454
127 David Clark Labs, Dept. of Applied Ecology, NCSU
Raleigh
NC
27695-7617
US
cbelyea@ncsu.edu
20190912
ArcGIS Metadata
1.0
Curtis Belyea
Biodiversity and Spatial Information Center, Dept. of Applied Ecology, NCSU
Biologist, GIS and Database Specialist
919-513-7337
919-515-4454
127 David Clark Labs, Dept. of Applied Ecology, NCSU
Raleigh
NC
27695-7617
US
cbelyea@ncsu.edu
File Geodatabase Raster Dataset
sesleuth2100
2016-05-23T00:00:00
Curtis Belyea
Biodiversity and Spatial Information Center, Dept. of Applied Ecology, NCSU
Biologist, GIS and Database Specialist
919-513-7337
919-515-4454
127 David Clark Labs, Dept. of Applied Ecology, NCSU
Raleigh
NC
27695-7617
US
cbelyea@ncsu.edu
This dataset represents the extent of urbanization (for the year indicated) predicted by the model SLEUTH, developed by Dr. Keith C. Clarke, at the University of California, Santa Barbara, Department of Geography and modified by David I. Donato of the United States Geological Survey (USGS) Eastern Geographic Science Center (EGSC). Further model modification and implementation was performed at the Biodiversity and Spatial Information Center at North Carolina State University.Urban growth probability extents throughout the 21st century were projected for the Southeast Regional Assessment Project (SERAP), which encompasses all or parts of the states of Alabama, Arkansas, Florida, Georgia, Illinois, Indiana, Kentucky, Louisiana, Mississippi, Missouri, North Carolina, Ohio, South Carolina, Tennessee, Virginia and West Virginia. Urban modeling for a large portion of the Appalachian, Gulf Coastal Plains & Ozarks, and Gulf Coast Prairie Landscape Conservation Cooperatives, which encompasses parts of the states of Texas, Oklahoma and Louisiana, Arkansas, Missouri and Illinois. Urban modeling for a large portion of these areas were completed earlier by the Biodiversity and Spatial Information Center as part of the Southeast Regional Assessment Project (SERAP). These data are supplementary to those projections.
Urbanization in the southeastern United States puts natural ecosystems and agro-ecosystems at risk. A picture of what future urbanization might look like, based on recent trends, may be of value to urban, conservation and transportation planners alike.
Biodiversity and Spatial Information Center; Dept of Applied Ecology at North Carolina State University, Raleigh, NC
Curtis Belyea
Biodiversity and Spatial Information Center, Dept. of Applied Ecology, NCSU
Biologist, GIS and Database Specialist
919-513-7337
919-515-4454
127 David Clark Labs, Dept. of Applied Ecology, NCSU
Raleigh
NC
27695-7617
US
cbelyea@ncsu.edu
West Virginia
Arkansas
Georgia
Pennsylvania
Indiana
Florida
North Carolina
Alabama
Mississippi
Kentucky
South Carolina
Tennessee
Ohio
Southeast
Virginia
Louisiana
Oklahoma
Texas
United States
Illinois
Missouri
Growth
SLEUTH
Urbanization
Project Gigalopolis
Prediction
Model
Development
Development
Kentucky
Indiana
Alabama
South Carolina
Georgia
Missouri
West Virginia
SLEUTH
Southeast
Florida
Growth
Virginia
Louisiana
United States
Oklahoma
Project Gigalopolis
Prediction
Tennessee
Urbanization
Model
North Carolina
Illinois
Texas
Mississippi
Arkansas
Ohio
Pennsylvania
This data set is not intended for site-specific analyses. Interpretations derived from its use are suited for regional and planning purposes only. These data are not intended to be used at scales larger than 1:100,000.
60
Version 6.2 (Build 9200) ; Esri ArcGIS 10.7.0.10450
true
-103.143808
-72.611193
21.815928
44.558369
true
-96.607701
-88.421865
29.603247
39.266971
1
-103.143808
-72.611193
44.558369
21.815928
No restrictions.
Prediction is run with the optimum set of coefficient values determined above. Output values in the Color Table were derived in order to capture probability of urbanization with a 95% confidence interval. Other important settings included: Monte Carlo Iterations = 200; Number of working grids = 6; Prediction Stop Date = 2101
Mosaicking was performed with the ArcGIS geoprocessor and the Python programming language via custom written code. For each subunit, the extent and mask environments were set to that of the study area. With these two settings, single output map algebra was used to Set null values (NoData) to 0 as such: gp.singleoutputmapalgebra_sa(“Con(Isnull(inputRaster),0,inputRaster") Cell statisctics was then performed on all subunits in the study area for each year to produce the mosaic raster for the study area for that year as such: gp.cellstatistics_sa(rasterList,studyarea_year,"MAXIMUM")
Once the initial optimal Auxillary Diffusion Multiplier was arrived at, the maximum fractional values for Area, Edges and Clusters were found. This was accomplished running calibration with the initial optimal Auxillary Diffusion Multiplier and the following coefficient values: Diffusion: 100 Breed: 100 Spread: 100 Slope: 0 Road Gravity: 100 After this step was run, the resulting log file was opened in Microsoft Excel, and the highest values for the time periods between input urban data time stamps for Fractional Area, Fractional Edges and Fractional Clusters were recorded. These were used in later calibration steps to evaluate results by weighting error indicated by Fractional Area, Fractional Edges and Fractional Clusters.
After each calibration step performed as described above, all log files produced were concatenated into a single file, which was opened in Microsoft Excel. Each of the Area, Edges and Clusters fractional measures for each time period between input urban data time stamps from the calibration were divided by their respective Maximum Fractional values determined in step 3 and summed. The absolute values of these summations are minimized through sorting, and the coefficient combinations with the lowest absolute weighted fractional error values are used as indicators for the ranges of coefficient values to be used in the next calibration step. It was observed that in initial rounds of calibration the diffusion coefficient among the best performing sets of coefficient combinations was often 1, indicating that lower values would provide better fitting results in calibration. As a result, the Auxillary Diffusion Multiplier was reduced until diffusion coefficient values other than 1 started appearing in the better performing combinations of coefficients. This step was repeated until optimal coefficient values were arrived at. The top performing combination of coefficients was then used in prediction.
After prediction was completed, the output .gif images were converted to geo-referenced .tif files. After geo-referencing was completed, projection information was defined and the .tif files were converted to ESRI format grid datasets and reclassified to remove hillshade values used as backdrop using the ArcGIS geoprocessor and the Python programming language. Values in the raster attribute table were interpreted as follows: 1 - Extent of existing urbanization at the time of the most recent input dataset for the model SLEUTH; 3 - 0-2.5% Probability of urbanization for the projected year; 4 - 2.5-5% Probability of urbanization for the projected year; 5 - 5-10% Probability of urbanization for the projected year; 6 - 10-20% Probability of urbanization for the projected year; 7 - 20-30% Probability of urbanization for the projected year; 8 - 30-40% Probability of urbanization for the projected year; 9 - 40-50% Probability of urbanization for the projected year; 10 - 50-60% Probability of urbanization for the projected year; 11 - 60-70% Probability of urbanization for the projected year; 12 - 70-80% Probability of urbanization for the projected year; 13 - 80-90% Probability of urbanization for the projected year; 14 - 90-95% Probability of urbanization for the projected year; 15 - 95-97.5% Probability of urbanization for the projected year; 16 - 97.5-100% Probability of urbanization for the projected year.
Prior to running the SLEUTH-3r model, the input ESRI format grid datasets were converted to .gif format images in order to conform to the model's requirements. Once input datasets were converted to .gif format, the first step in calibration, finding an initial optimal Auxillary Diffusion Multiplier, was performed. The scenario file for the study area was completed, with the following coefficient values: Diffusion = 100 Breed = 0 Spread = 0 Slope = 100 Road Gravity = 0 After this was run, the resulting log file was opened with Microsoft Excel, and the values for Fractional Area were examined. Ideal values for the Fractional Area were determined to be between +0.15 and 0.30 for the time periods between input urban data time stamps. If the fractional values were too low or too high the Auxillary Diffusion Multiplier was changed, according to their direct relationship, and the step was repeated.
Prior to mosaicking, grids produced from outputs of the SLEUTH model were reclassified once more. Bin values were replaced with upper-limit probability values mulitplied by 10 (to avoid the need for floating point rasters and make interpretation of values more straighforward) as described below: gp.reclassify_sa(dataset, "VALUE", "3 25;4 50;5 100;6 200;7 300;8 400;9 500;10 600;11 700;12 800;13 900;14 950;15 975;16 1000", output_dataset) where the value 3 was reclassified to 25 (upper-limit 2.5%),...16 was reclassified to 1000 (upper-limit 100%).
Calibration was initiated using the initial optimal Auxillary Multiplier, 25 Monte Carlo iterations and testing coefficient values as such, where no significant relief is present: Diffusion: 0 - 100; step interval 5 Breed: 0 - 100; step interval 5 Spread: 0 - 100; step interval 5 Slope: fixed at 25 Road Gravity: fixed at 100 due to the fact that it has been shown to have no significant relationship to fit statistics. Where significant relief is present: Diffusion: 0 - 100; step interval 10 Breed: 0 - 100; step interval 10 Spread: 0 - 100; step interval 10 Slope: 0 - 100; step interval 10 Road Gravity: fixed at 100 due to the fact that it has been shown to have no significant relationship to fit statistics.
2
11686
60.000000
17609
60.000000
true
false
-51951.227000 748292.078000
649208.773000 1804832.078000
649208.773000 748292.078000
-51951.227000 1804832.078000
298628.773000 1276562.078000
Projected urbanization
gcpo_urb2020
1000
0
Unified Code of Units of Measure
16
false
false
false
false
false
Pixel
39089
40991
60.000000
60.000000
16
1
Upper Left
FALSE
LZ77
1
pixel codes
FGDBR
TRUE
row and column
60.000000
60.000000
/9j/4AAQSkZJRgABAQEAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0a
HBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIy
MjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADIASwDASIA
AhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQA
AAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3
ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWm
p6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEA
AwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSEx
BhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElK
U1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3
uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3fevq
Khnjnk2mC48rHX5Q2anxSbccjikBUMF5gYvBxjP7sc1ajDLGodtzActjGTS5x1FG4UgFpM+1H4Uf
hQAE8c8UyJ9yYGDjiqs66hNKxt3jijAAxKuST3PHYjH5U6wjuozKLmSNySCPLXbigC5zSHODzTqQ
9DQBR/sq3acTzJ506tkSHgjHQcUwaFpYXAs0A9ATx+taNIWUdSB+NACRxpFGsaKAqjAHtSdWKjAp
wdScZGfrSD5nyOgFMAVduecn1NK2cYpcd6QHn+VIA5x1/Kmsg2kbRnFOFHc0wBenHApcc5poPyj6
U7HrQAE00cOR7CnY5FNIxIuPQ0AO60tIOpqK4LLAxWTyyP4sZxSAS6uobSBpp32RL95sE4rmtb8Z
2tpGI9PKXEzDO8cqn+Jp+p6jef2VcRJbzXTyJtAMOByQD0PpmvNyCGIIwQec9q78HQhUvKXQ8zMM
VUo2jDqbNz4r1m6g8l7wqvcxqEJ/EVjtud/MJ3MfmJJ7+9NpQQDmvRlRil7qPJhiJ3vOTYpGD70l
BYnn0pKqCfKuYipJcz5QJ711nhLw1HqKte30RMCkeUp6Oe+fYVx0kuxgoIGe9et+GrZrXw/BDKdz
qWHXI+8a4cXXXK4xZ6uBwkotVJre5ppDaxqI1hiUDoAoqKaNAknkQgOwwSg2s2BwAfxo8wg7QuPe
nurrDmLLMeV5x+NeYldXZ6t9bFHYnlqkmnSNjAHmKpP86Ibe3kdm+wrEcY+ZF5/KkSLUgy+bfRuA
fm/cYJ/HPFXK0pwad2Y1Jp6IRVVFCqoVQMAAYApc0VA7/NXQlcwbsaPmL6j8TWW02peRd7WQSjiH
IBBO48/Tbt/WtesySLU2tr9Vl2ylibZgV6Y4ByPWuM7x8El6/wBn80Bd0LGXGDtfjH9aijk1JU0/
zMsWBNydoGDxx/Opyl5/aKDCm08khmON2/selVRb6mtvZDzN0gmJuM7eY8NgDj12n86ALVi961zd
C5UCEP8AuTxkjJ9O2Nv61eqnbRXAnuTOFMZkzDgDhcdPzq1tX+6KAHUyRRkNjOODTti+lBRT2pAJ
sT+6KCif3R+VNAIkwCdoHINSd6AG7EHOxfypFADHAH4VJTcYyetAA4ypoQYRfpS9eajDAfIenY0A
P5Jx2rP1HU/7Pkj8xVET5G9mwM4+laC9Kja3ieQs6Bjx94Zx9KAMpPEdmHxMwhHIznPIxweOvNal
tcJdQLNEwZG6HHWlNvCZC3kx7yMElR0qRAEUKoAHoKAM6y1JXluoLl0ilgl24Y4yp+6ffNXp5fJi
Zz0AzXk/iS9Goa3dTKCqhti888cZpf8AhK9ZSO1jjucG3B68+Z/vetdVXDyVnFbnNh8RCfMpSSsz
0VbvWGKFdPheM/xGbaT6HGOP1pjajqpTzP7OgVUY7iboEY5744rn7P4jW5tnW6tGW6TgCM/I359K
47+0rtUvELjZdcygD73JP+NZ06Tnc2qydNq631+T2Z62NY08362gvrczP0QSAnPpVm4kgiiJuHRI
zwTIQAfbmvELd/s93DOh5Rg24dRg5rQ1XVrjV75rieQkYwqAEKAParjhnKfJcirXhCn7RbHb+IvF
LaZJBDprW8iSRbj3256Hg8cV56xLMS2cnk01SSOcde1LmvVw9ONOKt1PAxlWdSo+boFKD1pKK6Gc
aCiilApSeg4p30IZE3Nu7gcVt6P4nvtMCWcQSSAE/Jt5yevNZJAIIyKlsbr7FqMNwIxJ5R3bW6Gv
NqUI890tH+Z71HFylScZPVLT0S3/AFPTf7JLhZXu7oOx3sglIGeuPYZqWPSXbY51G9CqAAok6+9Y
Fj42haJEvomEucM8fTFbUfiXSRtYX8eAM4Oa5XhqkVZo2hi6UndSNWRdgQDcQBjJ6n60ysnT/EC6
zevHbW8gt4wSZX7nsMVrVSjKKtLcTnGbvHYZI2B71XPWpJGyTUVaxRjJ3Zs4pOaWivPPTEzS5oqt
cTzRSBYrUygjqGAwfxoAs5pCRVeCeaSUrJatGuD8xYHP+f6VZHQUAJk+lV7prgqEt/lkJB3kZAAP
IP1Gask4pAaAMktrCOJJfsmwYBVVYH36mkS11STdINSMYYlghhU7c9s98VrMMoc+lOHQUAY1yLmx
hW7utUKxRYMpMXDDPTA/pXOX/j6YXe2xgjNup+9IDlx/SrHjjXYVt20qMFpWIMjdlHXH1rgieRxX
o4bDwceaaPKxmKmp8lNnskWq2j2EV488aRSKCCzYHParPmI4BDKVIyDmvE977QhY7R0GeK7Lw7f6
WdMP2+6aKaLK4LkAqemBWdbBuEeaOppQx8ak+SWnmdusyRRnecKucmoxqljs8z7VEFPTLYziqNja
WbwbbWdpIXfdkPnHfH+fWrc/2C2tpZZRCIoR852g7fyrjSb0PQ5la9y2jrLGrxsGVhkMpzkVBqN7
HpunzXchG2NSQCcbj2H4mqt7rumaUVjnmEZMYdEVScr2xgV51r+vy65cq5QxRIgUR7sjPrXTQw0q
ju9EcmKxkKMbJ3kZl3L9qu5ZygTzHLbV6DPaqzbt6EHjoak696aNxHPHWvSkovRf1Y8uk5wXPLVe
fncYVXdkjjnOR1P+c1IDnI7VE65LMTvx29KWJSF6/l1rKEmpWtudVaEZUuZy1Vl9/wCOl7bLp2H7
QBjihcgHdxjv7Uc+2PWjI98jpW3MnqtDkcJJNPW/5/1uIjKzEqDT+tNA+XB/GnDp15q4X5bSMqvL
z3ht56h0GKKCBtzmjII7c1aa2Rg4t+8woVhg45pQB1NAXAHoKUmxwS6jVyD1pN3P3Ruz+lOxk01h
nt/SuapT/lPQoVrv95YUN8oI59/WlJyeKQDAANCkDrnr61rFyVuZnPUjTlzci6/geheC4Wj0RpGG
BLKzL9AAP5g10TcKawNDku5tLtmszaLbAAbDu3D1B9632+6a4Jy5ptnpUo8tNK5WY02lPWmN1q0Q
za5ozS0ySWOFQ0jhQTjJ9a809UdnPSlqD7Zb4B85Pmzjnr2qZWV1DKQVIyCO9ABjHSgMMkZ6UtNA
AcnAyaAKwnusjda9j/F04p8E1w8hWW38tf727NWD1oBBGR0oAR/un6UE4XJIAAzk0rfdOfSuQ1zx
np6Wc1pZP9omdTGWX7q5GDz3q4Qc5KK6k1JqEHN7I4jVLpr7VLm5faxdzg44x0H6YqltJcvnjGMU
7IHU0Y4r3eSN7LofMqpUScnsxASen606korSKsYzlfZWLVvqN7ax+XBcyxpkthGxzjFNivbiETBJ
WxMpSQHncDVelAJ6VPKlrYpTlJpXJJrie5cPM7SMFC5Y5OBwKj+lIeBxSE9h1qOdR93qbeylNKdt
AI+Ujpn0pEXC4zk0cFs5wT6GjO1cDkgVnFXlzs6Jtqn7JbvX+vIGUbgenrigABiBjjrUY3SKpIwe
uPSnH7oO4j6VKle8kjWdPlSpylrtpr/w/wDXQUcHGefSlHPpSKuCWP0FCjbnGOtXCNrXX9dDCrUu
3yvXT7+o4nGcj6Ype1NAI+8c/h0oJHOBz71cW4+9IzlFT9yC/r8BTz/Wm54AHH4UMwVeTzTHZt3B
6AYwOayqTULt7nRh6MqzSW3T/O3UmTjj0FOz0wRz2qs77FP94DimrMoY+ZnOQBuPOKwlikrLodkM
rlO8uvy8tSzxn6UZ3cfrTGP8I/D6U5TkDbzjjFaqurpM5ZYOSg2l1t5IUDtmkcHbwCacPeiulrmW
hwKThO8tbHS+DtVSzu2spmxHNypPQN/9eu9PQ144CQQR1FeneG7+TUdFilmOZFJRj647/liuXE0r
e+jswdfm9xltutRt1qV/vHFQk5NZI6GbP7r/AGRSNFBMNrKrj0PNSZ+v5UjEZU+/9K849Uj+yW+A
PJXC9OOlPWJUUKvCjgAGorm6Nvt2wSS5/uDOKhGoN/z6XJ64+SkBd2j1P51UvruGwg82UthnVFAP
UngVZhl85C2x05Iw4wayPEM1lLpl7ZzEPMLdphGOoA/i/Orpq8kmZ1ZOMG1uN1/VrWy0V5RIZPPU
pGY2zkketS+GcP4csWDE/u8Hnvk5rylp5WiWJpm8pf4S3A/CvRfD8s91pNsNNvbfy40CujJ0YZz7
jPB5rorwjShyJ63ObDTnWn7VrS1i34p1FNP0ScLNtnlGyNc8nPX9M15LGgDkYHXiuk8YPMfEMizY
3LGgyvAPyjJ/PNcy+7zVUHqegNbRhGlTjPciM516tSknyq1vuJZG2ZJG4YPSnA7gDgjPrS9RyKOM
16EU1JyPGqSi4KFtV1uFBpSec0laMxWug35ieOMU7vim5+bH589KdwOlZQlfzOmrDltpbsByR/jT
Qu3HP1oVucM3PX2p3HUg5qE1U95GjU6F6cv+H0EIJHHXtUTl+V6e9Stypxwahfcdqg5JFZ4iy7nR
l6ctNPn063/AIpQwwW5+uac0g3lOwHJpkIYFkZcYPB7U9owrZJyCMH3rGnObprudmIp0YYh9unrZ
a/d+PkPByoOeM8YoLYHrUSFVXaM8U5VL4Y8D09a1jVcoqK3ZzVMNCnN1JaRW3/DbgG3kjJA7U4DB
GeSR0oLKOBwT0qGQsjqxI3EYwKJWp+83dr+v+GCkpV3yJcqey7/Pt3H4BYs4+bOADSbxuYbcgc5o
LCJshcljyc0xmy5c9COCK46lX7K76nq0MNd+0km1ZW1/BJPy/NjWjLvtU9DkjGKmSHDMGXIPc96R
NoDLuB5+tWEPGeoq8PQhK7lv6mOPx1WCUKezW9u9n/XYjZSR8uM570oHydMjH5048jccf/XpMEjG
QK19l1XY4/rO0ZaXd3f/AIH4B1xxjHvTs5B5700E7RkEe3pW7pnha/1GETnZDE3Kl+rfhW8JqKuz
jrU3KXKtbeexT0jSZ9XuxDF8qjl3PRRXo2n6fBpFl9ngLFc7iWOST/kVBoGkjSNOETYM7EtIw7nt
+lXpm5rGrVdSVlsbYegqUeZ7kDGm0E5NFSaGlZ6hHemQJFLGUOCJE29z/hVk/eWlwKQgb1rzD1he
KOewpaMjp3oAazAIxIPArxm7lki1C6j8+RyrNHvYnLLnvXpfiDWBZ6fJLaOsksE6LKqnpz0PpnpX
FeLdP+zasblShhuh5ibTnkYyK7cLGSmn3OHGVKbptPpbYwMV1Xg3U9J0rzvtRaK4kbAk5IK/h05r
laK9KrQjUWp41DFzovTY3vGFxBd6+ZbeZJU8pRuQ5FYGB/WnUme9EKahFRQp1nUk5tiZxgetLRRW
vvIyXLLyD60ZppYggDqaUY7dqxjNSe50zpOmk2t/6vqGB170ySQYH+1/Kn4GD2pgVzjACjHTvWdS
8FaHU6MO4VHz1enfTp/XQkHTGBxQcY+lA+lB5Jz3rWMPd0OWdT39b2Aj1qMoA4IBz1zSuGIwO/HF
I4OwAHB7Gsqsk21bY6sLTcUnzfF+HqDAt9emPel5xyR74pBlV65I6imod0mccZ5OayUlFq+7Ol05
TUlH4Y636/1/WwoBfJGMdOKY0pZykakYOCSP5VINq4Gfbr0ppYRjaD83Tk5qW2tL27/5GkIxk78v
M7e6nt6vvshTHuIJODjnFRtMoOVU896TkEkPlj2J6CiKJmkYSLgZ6+1YTrucuSmtTuo4ONKHtcRK
8VstvWy662+7YETMhJ/UVKULpjaAen0p5jB2jHToaCdgH6+1bQpKEXznFVxUq04uk9V+H/D9hsS7
RjAwT1FS9O1MHYZ6dDT+gxn8q0pXiuVI5sVac/aSer/q+ont0pB05HWnDk8/yrqNB8Ki+gS8vHZY
W5SNeC3uauclHVowpxcrJPsVPCmmR6jqjNMoaKEbyp6E9h9K9FAAAAGAOwqjp2j2elNKbRGXzcbs
sT0z6/Wr9cs5Xeh2Qi0tRGOBVaU1MTk1Xk5NOKCbI6KKK0MzcAGKQqTyDyKitXuHiJuYljfcflVs
jFT15Z6wzeQRuAAPfNP4Jz+tV72/ttOtzPdSrHHnGT3PpXIRePIPNn86yOwH915Z5YZ75ranQnUV
4owq4mlSaU3Yy/GlsbLWXMRkWO6UO4z8rNmubaR2VVZ2Kr90E8D6Vr+INebW7hCI/Lgiz5aE5PPU
k/hWIZFUHkfTNexSvCmufRngVY+2rS9krpjsc4opiSK+cHGOKfWsJxkrxMatGdKXLNWYUYooqmrm
adgpAGwcnn1FL15pevepcb7mkanImkM2Y5B/OlDZHFLTQB271lycktNjo9r7WD5t+/6DwaOM5pMc
0fhWtkc3Mw/nS0lFNrqJSdrWA5PTPtTGJKkA80/BBqNwgBbHJrmndJ2O+jyuUU++nr5jQA4DM2M4
zinqm3p0/nUZHJIbB9alTcOp5PoK56cFz7HfiKslSvzWT6f8DbYryuPMI43Dnj+f+fSnEREFyBu7
05owSXHXHQ1E6MFxzz1zXPOUqLfMr9j0KUKWJjH2c3G1r627XXoNjUcBjjPHXtVlSAuAeAcUwAKm
C3IPUf596euVXBXHJHHelQg6aurDx1dYhpO9k/6tdf1cBKNpG7dyce1KpLA+oPWkKKzkMKkHuRzX
TCnKdufY8yrXp0k1SWr1flp028+i769DoByAMUZJH6cUAcdKdHGzvtjUl24AHJP0rZwcfh2ORVFP
Wa1ZJbxvNcJDEpZpCAF9TXrkaLHGqIoVVAAA7CsHw3oCabbCe4RWupADyP8AVj0HvXQE4FYVpJu0
TehCUVeQdKjLZ+lBbNJWSRs2ITgVC/Oakc4FR1pEzkR0U4r6U2rJNzPrR1NHFBORwcGvMPWPL/Fm
qXl3q01rP8kUDkJGP5n3Nc/VvUvN/tS689i8olbcx4yc1TPbrX0EEoQSR8pUbqVW5CO5XnGRUBnE
jEKM8Zp0qlvlz1HeokeJflXjBz0rzq1WU6nI3yrzPosHhadOh7WMOefleyV93+Qis27G3GOc1ajD
bQcj6VEQZHR8YX0qwowAB0xVYOlacpX06GWb4pSpQhypSe/W26+8WmkEjg4NLnGaM+nSvQbT91nh
RjONpxQCijIJxR2qotNaEzi07sPpSE880uKTr/Kpm1tfUulF7207i9uKO9FITyBTclFXZMacpuyF
opfwpKq5nYQ9Bx1pojGQe+fyqQZ6dKTBx7mudwTd5I7Y1pKPLBka4Q7ex5yafkg9se9BAYBSTjuP
WgqMdScVlFSWkEdU5wlrVlrs+vzK7SIW4GQKcCpIJzknpUpUEDjoeKAvzZJrP2E+a8n+BusbR9ny
wTXz1/L7/uIZGxKnbORnHbNSoSOqj6ClKKTk/gfTNKBzzn0qadGSm7l18VTlRio9PzuHLt0IPp60
ZIJGMAU7GTuwPakVQFA9qqUJuSV9LmUKtJQbcbu39f0xckEFeTnFegeF9COnwG5uFXz5ACoxyg9P
rXK+G9PTUdYijl5jQGRh64PSvTQMDAqZuUY8rd2xRUZy5lGyWwUxzzinE4FR1kjZhRRRTJIn600n
FObrUbHmtEZsQkmkooqhG7ntSOu9GXJXIIyOoqK5ure0jD3MyRITgFzgZqBdQt7tmhsrqKSYAN8p
3ADPfFeaou17HqOUb2bPL9f02fS9WlhmZnD/ADo7HJdT3PvWXXbeN7KXy4b24eIPxEqpnnv/AI1x
Ne7hqntKabPmsXSVKq4oQqCRSGNTkADmnelHTtTnGDequFKrWivdlZIaEAG3n86Xn1z+lOFJVeyR
H1idtdQI9OKO1FA6daEkpWsEpOUOZsCeenFFHGKKu1kZXcmHtSZ/+tQRuAwT70p9q52pTbfbY7Iu
FKKj1e9+mocE5xR3z7YopcH3xjNaTUVHVGNJ1JVLRf8AwwYpDxQOucEZ60HrRCUnuKpCEXoKf0pP
T2pQeOvbikI4yex60SlZXY6cW2kuovY0meR0oPB6jikIyKT2uOC96z6aBnOOOp9aXr1qMoqsWbIO
Kc7MOgOT6elc3tnFOUkd7wqnKMKUvnt/VlYUEjgnjpThn8KhjfMpDZIOcZGKlJ5Gcc1lh6nOubsd
GNw7oyVO2r/yFzj1wD2pQT6Dnik+tN+bdweMcit5w5I8yZxUqvtJuEl+Z1/gSHM13PkfKoTHfnn+
ldtXBeD9Rjs2uY2jkcuFYLGm49cf1rrW1eMIWFrdtjsIjXLVl712dlGFo2Xdl5zximUu7eqsAQCM
4IwaSkimFFFFAiN171CfvVaIzUT4XNWmQ0Q9KKaWy9OqyTl/FXiWz1i0itrWOQ7X3l34HTpj8a5i
C6uLbPkTyRZ67GIz+VS6jYTaZfSWk+0yJ1KnIqrXbSpwhC0djyq1WpOo5T3J7i8ubzb9puJJdgwu
9icCoDjvTVkRm2hhu9M807OR096OZOPuWK9nOM7VUxjEbjhsE84NOQYXrnvmqwKvISxKkMBycVZV
geBXLh6qnUb+49LHYaVGikl2vou3fsLRR/OmO+MAYJPSuupUjBXZ5dChOrLliiTFJSK4OcfjTgOx
pxmpLmQqlKUJcslZiUUDnoRQPei6bBRkl6CdyM4JHFL+NHQEjrikBBzWUWlPl6s6JxcqamtluL3w
aUetNyEUZI/GlznvVKTloZypxhrrbuKBgd/ak6H0+tHIPWrFsto7n7XLKi9vLQN/WhKavfYHKm7W
3/ruVgxz09uaXJIyO/Y12tp4Lsp7ZZftkzCQBkIUDg+vWtS18KaTbEHyTK3rK2f0rlnUja17ndCn
K9+W3z/M84jgmuG2pG0hI5CKSTT5rS5twGmglj93Qj+detwwQ26bIYkjX0RcU9lV1KsoZT1BGaj2
0UrKJfsajfM5fcrI8awzMQoxxyaV1DY5Ndv4k8N26Wcl7ZRski8vGgyGH07Vytppl3eXYtoYH355
3DG0ep9KunShNuUtUTWxVWmowpqzXb+vwM4JulY44AxSzLvXjO4HcMV6lpGg22l2zJ/rXkA8wuAR
nHIHHSsTxJ4ZjWJr2wjC7R88KL19xURVN3glubSq1041W/hSsvzv3T1uvM4TzGjYDkr33ccVKzbQ
Tn8zTJY3ZmUgHOcZ7YxRDtMYDEgZxz1rntOM3Te3S56MnSqUo14r3ra289fwaa+XoamjX7adqkFz
nKhsOB3U9a9VVgyhlOQRkGvGiDkYI/HvXf8AhHWRdWxsJnzND9wk/eX/AOtWtS/NyvdHBTUXD2kH
o3t2/rY6R+tNqRhkVHUobCiiigQE4qCbpmpCc01hlSKpaEvUrYoBNOHoe1LitbmZ594g1BNT1qe5
jB8snaue4HGazKKK74xUYqKPHnJzk5PdkbRpuD7RkUzeWA3Z25/hzz6VMRn61GEPQHGOwrhr02na
K3/M9vB4iMoqVSWse99Ftpv3EijyN7Dk9sYxipgAOgpAMAClrpoUY042tqefjMZOvUbvp+gHJ4FR
sBvOQOB1FSUmAckClWp3swwtbkurb6DFPzfLjFScZ6UtvbNNMscEbNI/AVecmu60PwnFboJtRjSW
Y8qmchfr6moi1Shee76Glb9/UtSei6s4iGCW4k2QxPI/XaikmnzWd1bqGmt5Y1PQuhFesx28MTFo
4kRiMZVQKc6JIpV1DKeoIyKj632Q/qF1rLU8dpNozkcH+deqPoWlSOGawgyPRcVkzeCbF95immRj
kqCQQKv6xSl8SIWFr0/gZwRA685oGa1NS0C+0yUI8ZkVvuvGCRUEWkajMyiOynO7p8hFbrl+JM5n
7S3I0UqsWNuLu9hgaRY1dgC7HAFauo+GLnTdNF3NIrNkBkQZ2/U1W0jRrnUbuIeW6Q53NIV4x7U/
aRcW0yfZTU1FrU9NhiWCCOJPuooUfQU4sBURbAAzUe7JryuW57nNbREuacG9aYvSlosCZL1pMAEk
AZPU01WxxT6llBRRRQM5PxB4ajZLvUIM+Zjf5ajA/wBo+9cMIlLe/UD1r2GeWOC3kmlOI0Usxxng
V5HIdzvIoChiWA9K3T54a9Dn/hVPc67kboRtAIGeoz2qxBLLatHPC5EqMNremKgMgUAsAcdyPWpA
4CBD0zng9utZybhbn95o2jH2l1T9xNPS3z+f5npehauNYsjKYyjoQrjtnHatAjBrkfBmpRxGTTpD
tZmLxk9/UfpXXv1pzVpabGdOTlBX3G0jHApaQjIpFjKKKO1MkhkGH470lPkHQ0yrWxDPK8c5oopy
EAnd0Ir0tjxdZOw2inhAQ2N3HcDimUKzB3CjvRSqrOwVQSxOAB3ounoFmtRK1tM8O3+pMpWIxQnr
JIMDHt610mieEYoFjudQ/eTfeEX8K/X1rqq5KuJS0gd9HBtq9TbsZWj6Da6OjGPMkrfekYc/Qela
tJketG4etcUm5O7PRjGMFaItFJketGRSKFopu8UbxRYVx1FM30hb1NOwXBucg4IPY0wkAYGKC3pT
apIhsQ9KbStSVaIY4ORTxIO/FQscUzrRy3DmsWgQehpwYiqgYqeKsK24ZqXGxalcnHNFNTpTqzLG
TRiaCSI9HUr+YryW8tns7yW3YhmicqSO9eu1zPijQ1uYzqFvH/pEeC4H8a//AFq0hNqLiupnOEXO
M5bLf0OECrjt0Oc1F5p4Uj7vfPH0z+VPfLEFFJbvUDGOJvMdSeMAehrHEVHTajJK/wCB3YGhGupT
hJ2W1vivr0v/AFr1JFkmtr6K5ibaVbP416xp2oR6np8VzHxnhl/ut3FeSebHMq7Xz6AHmus8F3jp
fS2hJ8uVN+CejD/61bRhfms01voc1ao1GmpxcXbl1Vr2bs/0O5UAnmn7R6VHT1cfxD8qhkoCgqMx
+lKZPm6cU84deDRqgsmVZBjrnNR1LcMGkIAwBUVax2Mpbnl89vNazGKeNo5B1VhiomXcuKKK9Fe/
DXqeO/3dT3ejHlsgA9qbRRVJJbEOTe4V1HgzTVubyS8lXKwYCZH8R7/h/WiissQ2qbsb4VKVVJnd
k4FRs1FFeaj2WN3e1G6iiqsRcUHNLn3oopDuNLUFjRRTsK4mTSUUUCE3UmTRRVCEooooACM1GetF
FNCCnI+00UUwLSH9akoorBm62CiiikM4LWfDN5Dd3V1bIrWwzJ97BA6kY/OuaRdj+ZIBg/d74oor
blVZx5jBVHhoz9n1RGY03MeuWLcitvwxPHba5E8sixptYFmOB0oorrdKFKk4wXQ4FiKuIrqdR3d/
1PQkuIZG2pNGzeisCan2e9FFedc9ZqzG+XmmMGjG5TRRVJ6ktaEWDNIf1oaFgcDmiiqbadiVFNXZ
/9k=