In this course, you will explore a variety of open-source technologies for working with geosptial data, performing spatial analysis, and undertaking general data science. The first component of the class focuses on the use of QGIS and associated technologies (GDAL, PROJ, GRASS, SAGA, and Orfeo Toolbox). The second component of the class introduces Python and associated open-source libraries and modules (NumPy, Pandas, Matplotlib, Seaborn, GeoPandas, Rasterio, WhiteboxTools, and Scikit-Learn) used by geospatial scientists and data scientists. We also provide an introduction to Structured Query Language (SQL) for performing table and spatial queries. This course is designed for individuals that have a background in GIS, such as working in the ArcGIS environment, but no prior experience using open-source software and/or coding. You will be asked to work through a series of lecture modules and videos broken into several topic areas, as outlined below. Fourteen assignments and the required data have been provided as hands-on opportunites to work with data and the discussed technologies and methods. If you have any questions or suggestions, feel free to contact us. We hope to continue to update and improve this course. This course was produced by West Virginia View (http://www.wvview.org/) with support from AmericaView (https://americaview.org/). This material is based upon work supported by the U.S. Geological Survey under Grant/Cooperative Agreement No. G18AP00077. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the opinions or policies of the U.S. Geological Survey. Mention of trade names or commercial products does not constitute their endorsement by the U.S. Geological Survey. After completing this course you will be able to: apply QGIS to visualize, query, and analyze vector and raster spatial data. use available resources to further expand your knowledge of open-source technologies. describe and use a variety of open data formats. code in Python at an intermediate-level. read, summarize, visualize, and analyze data using open Python libraries. create spatial predictive models using Python and associated libraries. use SQL to perform table and spatial queries at an intermediate-level.

Have you ever wanted to create your own maps, or integrate and visualize spatial datasets to examine changes in trends between locations and over time? Follow along with these training tutorials on QGIS, an open source geographic information system (GIS) and learn key concepts, procedures and skills for performing common GIS tasks – such as creating maps, as well as joining, overlaying and visualizing spatial datasets. These tutorials are geared towards new GIS users. We’ll start with foundational concepts, and build towards more advanced topics throughout – demonstrating how with a few relatively easy steps you can get quite a lot out of GIS. You can then extend these skills to datasets of thematic relevance to you in addressing tasks faced in your day-to-day work.

LANDISVIEW is a tool, developed at the Knowledge Engineering Laboratory at Texas A&M University, to visualize and animate 8-bit/16-bit ERDAS GIS format (e.g., LANDIS and LANDIS-II output maps). It can also convert 8-bit/16-bit ERDAS GIS format into ASCII and batch files. LANDISVIEW provides two major functions: 1) File Viewer: Files can be viewed sequentially and an output can be generated as a movie file or as an image file. 2) File converter: It will convert the loaded files for compatibility with 3rd party software, such as Fragstats, a widely used spatial analysis tool. Some available features of LANDISVIEW include: 1) Display cell coordinates and values. 2) Apply user-defined color palette to visualize files. 3) Save maps as pictures and animations as video files (*.avi). 4) Convert ERDAS files into ASCII grids for compatibility with Fragstats. (Source: http://kelab.tamu.edu/)

GIS Online Moisture Sensor Market Outlook

According to our latest research, the global GIS online moisture sensor market size reached USD 1.12 billion in 2024, reflecting robust adoption across key sectors such as agriculture, environmental monitoring, and industrial process control. The market is projected to grow at a CAGR of 8.7% from 2025 to 2033, reaching an estimated value of USD 2.43 billion by 2033. This impressive growth is primarily driven by the increasing demand for precision agriculture, advancements in sensor technologies, and the growing need for real-time environmental data to support sustainable resource management.

One of the primary growth factors fueling the GIS online moisture sensor market is the surging adoption of precision agriculture techniques worldwide. Farmers and agribusinesses are increasingly leveraging advanced moisture sensing technologies integrated with GIS platforms to monitor soil conditions, optimize irrigation schedules, and enhance crop yields. The ability to access real-time moisture data remotely has transformed traditional farming practices, allowing for data-driven decisions that conserve water and reduce operational costs. This trend is further supported by government initiatives and subsidies promoting smart farming solutions, particularly in regions facing water scarcity or climate variability. As a result, the integration of GIS and online moisture sensors has become a cornerstone in the modernization of agricultural operations, driving sustained market expansion.

Another significant driver for the GIS online moisture sensor market is the escalating focus on environmental monitoring and industrial process control. Industries such as construction, mining, and manufacturing are increasingly required to adhere to stringent environmental regulations, necessitating continuous monitoring of moisture levels in soil, air, and materials. GIS-enabled online moisture sensors provide accurate, location-based data that supports compliance, risk management, and process optimization. In addition, the proliferation of smart city initiatives and the expansion of IoT infrastructure have amplified the deployment of these sensors in urban planning, flood prediction, and infrastructure maintenance. The convergence of GIS and online sensor technologies enables seamless data visualization and analysis, making them indispensable tools for both public and private sector stakeholders.

Technological advancements in sensor design and connectivity are also playing a pivotal role in the market's growth trajectory. Innovations such as wireless and cloud-connected moisture sensors, improved accuracy through advanced materials, and miniaturization have broadened the scope of applications. These advancements have resulted in more cost-effective, durable, and easy-to-deploy solutions, fostering adoption across diverse end-user segments. Furthermore, the integration of AI and machine learning algorithms with GIS platforms is enabling predictive analytics and automated decision-making, further enhancing the value proposition of online moisture sensors. As the demand for actionable insights and real-time monitoring continues to rise, the GIS online moisture sensor market is poised for sustained innovation and expansion.

Regionally, North America and Europe are leading the market, driven by early adoption of precision agriculture, robust regulatory frameworks, and substantial investments in R&D. Asia Pacific, however, is emerging as the fastest-growing region, propelled by rapid urbanization, increasing awareness of sustainable agricultural practices, and government support for smart farming initiatives. Latin America and the Middle East & Africa are also witnessing steady growth, albeit from a smaller base, as industries in these regions recognize the benefits of GIS-enabled moisture monitoring for resource optimization and environmental management. Overall, the global market is characterized by dynamic regional trends, with each geography contributing uniquely to the market's evolution.

Prod

This resource contains data inputs and a Jupyter Notebook that is used to introduce Hydrologic Analysis using Terrain Analysis Using Digital Elevation Models (TauDEM) and Python. TauDEM is a free and open-source set of Digital Elevation Model (DEM) tools developed at Utah State University for the extraction and analysis of hydrologic information from topography. This resource is part of a HydroLearn Physical Hydrology learning module available at https://edx.hydrolearn.org/courses/course-v1:Utah_State_University+CEE6400+2019_Fall/about

In this activity, the student learns how to (1) derive hydrologically useful information from Digital Elevation Models (DEMs); (2) describe the sequence of steps involved in mapping stream networks, catchments, and watersheds; and (3) compute an approximate water balance for a watershed-based on publicly available data.

Please note that this exercise is designed for the Logan River watershed, which drains to USGS streamflow gauge 10109000 located just east of Logan, Utah. However, this Jupyter Notebook and the analysis can readily be applied to other locations of interest. If running the terrain analysis for other study sites, you need to prepare a DEM TIF file, an outlet shapefile for the area of interest, and the average annual streamflow and precipitation data. - There are several sources to obtain DEM data. In the U.S., the DEM data (with different spatial resolutions) can be obtained from the National Elevation Dataset available from the national map (http://viewer.nationalmap.gov/viewer/). Another DEM data source is the Shuttle Radar Topography Mission (https://www2.jpl.nasa.gov/srtm/), an international research effort that obtained digital elevation models on a near-global scale (search for Digital Elevation at https://www.usgs.gov/centers/eros/science/usgs-eros-archive-products-overview?qt-science_center_objects=0#qt-science_center_objects). - If not already available, you can generate the outlet shapefile by applying basic terrain analysis steps in geospatial information system models such as ArcGIS or QGIS. - You also need to obtain average annual streamflow and precipitation data for the watershed of interest to assess the annual water balance and calculate the runoff ratio in this exercise. In the U.S., the streamflow data can be obtained from the USGS NWIS website (https://waterdata.usgs.gov/nwis) and the precipitation from PRISM (https://prism.oregonstate.edu/normals/). Note that using other datasets may require preprocessing steps to make data ready to use for this exercise.

Notice: this is not the latest Heat Island Severity image service.This layer contains the relative heat severity for every pixel for every city in the United States, including Alaska, Hawaii, and Puerto Rico. Heat Severity is a reclassified version of Heat Anomalies raster which is also published on this site. This data is generated from 30-meter Landsat 8 imagery band 10 (ground-level thermal sensor) from the summer of 2023.To explore previous versions of the data, visit the links below:Heat Severity - USA 2022Heat Severity - USA 2021Heat Severity - USA 2020Heat Severity - USA 2019Federal statistics over a 30-year period show extreme heat is the leading cause of weather-related deaths in the United States. Extreme heat exacerbated by urban heat islands can lead to increased respiratory difficulties, heat exhaustion, and heat stroke. These heat impacts significantly affect the most vulnerable—children, the elderly, and those with preexisting conditions.The purpose of this layer is to show where certain areas of cities are hotter than the average temperature for that same city as a whole. Severity is measured on a scale of 1 to 5, with 1 being a relatively mild heat area (slightly above the mean for the city), and 5 being a severe heat area (significantly above the mean for the city). The absolute heat above mean values are classified into these 5 classes using the Jenks Natural Breaks classification method, which seeks to reduce the variance within classes and maximize the variance between classes. Knowing where areas of high heat are located can help a city government plan for mitigation strategies.This dataset represents a snapshot in time. It will be updated yearly, but is static between updates. It does not take into account changes in heat during a single day, for example, from building shadows moving. The thermal readings detected by the Landsat 8 sensor are surface-level, whether that surface is the ground or the top of a building. Although there is strong correlation between surface temperature and air temperature, they are not the same. We believe that this is useful at the national level, and for cities that don’t have the ability to conduct their own hyper local temperature survey. Where local data is available, it may be more accurate than this dataset. Dataset SummaryThis dataset was developed using proprietary Python code developed at Trust for Public Land, running on the Descartes Labs platform through the Descartes Labs API for Python. The Descartes Labs platform allows for extremely fast retrieval and processing of imagery, which makes it possible to produce heat island data for all cities in the United States in a relatively short amount of time.What can you do with this layer?This layer has query, identify, and export image services available. Since it is served as an image service, it is not necessary to download the data; the service itself is data that can be used directly in any Esri geoprocessing tool that accepts raster data as input.In order to click on the image service and see the raw pixel values in a map viewer, you must be signed in to ArcGIS Online, then Enable Pop-Ups and Configure Pop-Ups.Using the Urban Heat Island (UHI) Image ServicesThe data is made available as an image service. There is a processing template applied that supplies the yellow-to-red or blue-to-red color ramp, but once this processing template is removed (you can do this in ArcGIS Pro or ArcGIS Desktop, or in QGIS), the actual data values come through the service and can be used directly in a geoprocessing tool (for example, to extract an area of interest). Following are instructions for doing this in Pro.In ArcGIS Pro, in a Map view, in the Catalog window, click on Portal. In the Portal window, click on the far-right icon representing Living Atlas. Search on the acronyms “tpl” and “uhi”. The results returned will be the UHI image services. Right click on a result and select “Add to current map” from the context menu. When the image service is added to the map, right-click on it in the map view, and select Properties. In the Properties window, select Processing Templates. On the drop-down menu at the top of the window, the default Processing Template is either a yellow-to-red ramp or a blue-to-red ramp. Click the drop-down, and select “None”, then “OK”. Now you will have the actual pixel values displayed in the map, and available to any geoprocessing tool that takes a raster as input. Below is a screenshot of ArcGIS Pro with a UHI image service loaded, color ramp removed, and symbology changed back to a yellow-to-red ramp (a classified renderer can also be used): A typical operation at this point is to clip out your area of interest. To do this, add your polygon shapefile or feature class to the map view, and use the Clip Raster tool to export your area of interest as a geoTIFF raster (file extension ".tif"). In the environments tab for the Clip Raster tool, click the dropdown for "Extent" and select "Same as Layer:", and select the name of your polygon. If you then need to convert the output raster to a polygon shapefile or feature class, run the Raster to Polygon tool, and select "Value" as the field.Other Sources of Heat Island InformationPlease see these websites for valuable information on heat islands and to learn about exciting new heat island research being led by scientists across the country:EPA’s Heat Island Resource CenterDr. Ladd Keith, University of ArizonaDr. Ben McMahan, University of Arizona Dr. Jeremy Hoffman, Science Museum of Virginia Dr. Hunter Jones, NOAA Daphne Lundi, Senior Policy Advisor, NYC Mayor's Office of Recovery and ResiliencyDisclaimer/FeedbackWith nearly 14,000 cities represented, checking each city's heat island raster for quality assurance would be prohibitively time-consuming, so Trust for Public Land checked a statistically significant sample size for data quality. The sample passed all quality checks, with about 98.5% of the output cities error-free, but there could be instances where the user finds errors in the data. These errors will most likely take the form of a line of discontinuity where there is no city boundary; this type of error is caused by large temperature differences in two adjacent Landsat scenes, so the discontinuity occurs along scene boundaries (see figure below). Trust for Public Land would appreciate feedback on these errors so that version 2 of the national UHI dataset can be improved. Contact Dale.Watt@tpl.org with feedback.

This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years since 1992. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2019Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: AnnualWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

Background: Malaria continues to pose a major public health challenge in tropical regions. Despite significant efforts to control malaria in Tanzania, there are still residual transmission cases. Unfortunately, little is known about where these residual malaria transmission cases occur and how they spread. In Tanzania, for example, the transmission is heterogeneously distributed. In order to effectively control and prevent the spread of malaria, it is essential to understand the spatial distribution and transmission patterns of the disease. This study seeks to predict areas that are at high risk of malaria transmission so that intervention measures can be developed to accelerate malaria elimination efforts.

Methods: This study employs a geospatial-based model to predict and map out malaria risk area in Kilombero Valley. Environmental factors related to malaria transmission were considered and assigned valuable weights in the Analytic Hierarchy Process (AHP), an online system using a pairwise comparison technique. The malaria hazard map was generated by a weighted overlay of the altitude, slope, curvature, aspect, rainfall distribution, and distance to streams in Geographic Information Systems (GIS). Finally, the risk map was created by overlaying components of malaria risk including hazards, elements at risk, and vulnerability. Results: The study demonstrates that the majority of the study area falls under the moderate-risk level (61%), followed by the low-risk level (31%), while the high-malaria risk area covers a small area, which occupies only 8% of the total area. Conclusion: The findings of this study are crucial for developing spatially targeted interventions against malaria transmission in residual transmission settings. Predicted areas prone to malaria risk provide information that will inform decision-makers and policymakers for proper planning, monitoring, and deployment of interventions. Methods Data acquisition and description The study employed both primary and secondary data, which were collected from numerous sources based on the input required for the implementation of the predictive model. Data collected includes the locations of all public and private health centers that were downloaded free from the health portal of the United Republic of Tanzania, Ministry of Health, Community Development, Gender, Elderly, and Children, through the universal resource locator (URL) (http://moh.go.tz/hfrportal/). Human population data was collected from the 2012 population housing census (PHC) for the United Republic of Tanzania report. Rainfall data were obtained from two local offices; Kilombero Agricultural Training and Research Institute (KATRIN) and Kilombero Valley Teak Company (KVTC). These offices collect meteorological data for agricultural purposes. Monthly data from 2012 to 2017 provided from thirteen (13) weather stations. Road and stream network shapefiles were downloaded free from the MapCruzin website via URL (https://mapcruzin.com/free-tanzania-arcgis-maps-shapefiles.htm). With respect to the size of the study area, five neighboring scenes of the Landsat 8 OLI/TIRS images (path/row: 167/65, 167/66, 167/67, 168/66 and 168/67) were downloaded freely from the United States Geological Survey (USGS) website via URL: http://earthexplorer.usgs.gov. From July to November 2017, the images were selected and downloaded from the USGS Earth Explorer archive based on the lowest amount of cloud cover coverage as viewed from the archive before downloading. Finally, the digital elevation data with a spatial resolution of three arc-seconds (90m by 90m) using WGS 84 datum and the Geographic Coordinate System were downloaded free from the Shuttle Radar Topography Mission (SRTM) via URL (https://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Africa/). Only six tiles that fall in the study area were downloaded, coded tiles as S08E035, S09E035, S10E035, S08E036, S09E036, S10E036, S08E037, S09E037 and S10E037. Preparation and Creation of Model Factor Parameters Creation of Elevation Factor All six coded tiles were imported into the GIS environment for further analysis. Data management tools, with raster/raster data set/mosaic to new raster feature, were used to join the tiles and form an elevation map layer. Using the spatial analyst tool/reclassify feature, the generated elevation map was then classified into five classes as 109–358, 359–530, 531–747, 748–1017 and >1018 m.a.s.l. and new values were assigned for each class as 1, 2, 3, 4 and 5, respectively, with regards to the relationship with mosquito distribution and malaria risk. Finally, the elevation map based on malaria risk level is levelled as very high, high, moderate, low and very low respectively. Creation of Slope Factor A slope map was created from the generated elevation map layer, using a spatial analysis tool/surface/slope feature. Also, the slope raster layer was further reclassified into five subgroups based on predefined slope classes using standard classification schemes, namely quantiles as 0–0.58, 0.59–2.90, 2.91–6.40, 6.41–14.54 and >14.54. This classification scheme divides the range of attribute values into equal-sized sub-ranges, which allow specifying the number of the intervals while the system determines where the breaks should be. The reclassified slope raster layer subgroups were ranked 1, 2, 3, 4 and 5 according to the degree of suitability for malaria incidence in the locality. To elaborate, the steeper slope values are related to lesser malaria hazards, and the gentler slopes are highly susceptible to malaria incidences. Finally, the slope map based on malaria risk level is leveled as very high, high, moderate, low and very low respectively. Creation of Curvature Factor Curvature is another topographical factor that was created from the generated elevation map using the spatial analysis tool/surface/curvature feature. The curvature raster layer was further reclassified into five subgroups based on predefined curvature class. The reclassified curvature raster layer subgroups were ranked to 1, 2, 3, 4 and 5 according to their degree of suitability for malaria occurrence. To explain, this affects the acceleration and deceleration of flow across the surface. A negative value indicates that the surface is upwardly convex, and flow will be decelerated, which is related to being highly susceptible to malaria incidences. A positive profile indicates that the surface is upwardly concave and the flow will be accelerated which is related to a lesser malaria hazard, while a value of zero indicates that the surface is linear and related to a moderate malaria hazard. Lastly, the curvature map based on malaria risk level is leveled as very high, high, moderate, low, and very low respectively.
Creation of Aspect Factor As a topographic factor associated with mosquito larval habitat formation, aspect determines the amount of sunlight an area receives. The more sunlight received the stronger the influence on temperature, which may affect mosquito larval survival. The aspect of the study area also was generated from the elevation map using spatial analyst tools/ raster /surface /aspect feature. The aspect raster layer was further reclassified into five subgroups based on predefined aspect class. The reclassified aspect raster layer subgroups were ranked as 1, 2, 3, 4 and 5 according to the degree of suitability for malaria incidence, and new values were re-assigned in order of malaria hazard rating. Finally, the aspect map based on malaria risk level is leveled as very high, high, moderate, low, and very low, respectively. Creation of Human Population Distribution Factor Human population data was used to generate a population distribution map related to malaria occurrence. Kilombero Valley has a total of 42 wards, the data was organized in Ms excel 2016 and imported into the GIS environment for the analysis, Inverse Distance Weighted (IDW) interpolation in the spatial analyst tool was applied to interpolate the population distribution map. The population distribution map was further reclassified into five subgroups based on potential to malaria risk. The reclassified map layer subgroups were ranked according to the vulnerability to malaria incidence in the locality such as areas having high population having the highest vulnerability and the less population having less vulnerable, and the new value was assigned as 1, 2, 3, 4 and 5, and then leveled as very high, high, moderate, low and very low malaria risk level, respectively. Creation of Proximity to Health Facilities Factor The distribution of health facilities has a significant impact on the malaria vulnerability of the population dwellings in the Kilombero Valley. The health facility layer was created by computing distance analysis using proximity multiple ring buffer features in spatial analyst tool/multiple ring buffer. Then the map layer was reclassified into five sub-layers such as within (0–5) km, (5.1–10) km, (10.1–20) km, (20.1–50) km and >50km. According to a WHO report, it is indicated that the human population who live nearby or easily accessible to health facilities is less vulnerable to malaria incidence than the ones who are very far from the health facilities due to the distance limitation for the health services. Later on, the new values were assigned as 1, 2, 3, 4 and 5, and then reclassified as very high, high, moderate, low and very low malaria risk levels, respectively. Creation of Proximity to Road Network Factor The distance to the road network is also a significant factor, as it can be used as an estimation of the access to present healthcare facilities in the area. Buffer zones were calculated on the path of the road to determine the effect of the road on malaria prevalence. The road shapefile of the study area was inputted into GIS environment and spatial analyst tools / multiple ring buffer feature were used to generate five buffer zones with the

Chicago sites that offer free or affordable technology resources and services, like computers with Internet access, Wi-Fi hotspots and technology training. Call or visit the organization's website before going to the location. For more information, visit http://locations.weconnectchicago.org/.

This web map presents a vector basemap of OpenStreetMap (OSM) data hosted by Esri. Esri created this vector tile basemap from the Daylight map distribution of OSM data, which is supported by Facebook and supplemented with additional data from Microsoft. This version of the map is rendered using OSM cartography. The OSM Daylight map will be updated every month with the latest version of OSM Daylight data.OpenStreetMap is an open collaborative project to create a free editable map of the world. Volunteers gather location data using GPS, local knowledge, and other free sources of information and upload it. The resulting free map can be viewed and downloaded from the OpenStreetMap site:www.OpenStreetMap.org. Esri is a supporter of the OSM project and is excited to make this enhanced vector basemap available to the ArcGIS user and developer communities.

Please note that this dataset is not an official City of Toronto land use dataset. It was created for personal and academic use using City of Toronto Land Use Maps (2019) found on the City of Toronto Official Plan website at https://www.toronto.ca/city-government/planning-development/official-plan-guidelines/official-plan/official-plan-maps-copy, along with the City of Toronto parcel fabric (Property Boundaries) found at https://open.toronto.ca/dataset/property-boundaries/ and Statistics Canada Census Dissemination Blocks level boundary files (2016). The property boundaries used were dated November 11, 2021. Further detail about the City of Toronto's Official Plan, consolidation of the information presented in its online form, and considerations for its interpretation can be found at https://www.toronto.ca/city-government/planning-development/official-plan-guidelines/official-plan/ Data Creation Documentation and Procedures Software Used The spatial vector data were created using ArcGIS Pro 2.9.0 in December 2021. PDF File Conversions Using Adobe Acrobat Pro DC software, the following downloaded PDF map images were converted to TIF format. 9028-cp-official-plan-Map-14_LandUse_AODA.pdf 9042-cp-official-plan-Map-22_LandUse_AODA.pdf 9070-cp-official-plan-Map-20_LandUse_AODA.pdf 908a-cp-official-plan-Map-13_LandUse_AODA.pdf 978e-cp-official-plan-Map-17_LandUse_AODA.pdf 97cc-cp-official-plan-Map-15_LandUse_AODA.pdf 97d4-cp-official-plan-Map-23_LandUse_AODA.pdf 97f2-cp-official-plan-Map-19_LandUse_AODA.pdf 97fe-cp-official-plan-Map-18_LandUse_AODA.pdf 9811-cp-official-plan-Map-16_LandUse_AODA.pdf 982d-cp-official-plan-Map-21_LandUse_AODA.pdf Georeferencing and Reprojecting Data Files The original projection of the PDF maps is unknown but were most likely published using MTM Zone 10 EPSG 2019 as per many of the City of Toronto's many datasets. They could also have possibly been published in UTM Zone 17 EPSG 26917 The TIF images were georeferenced in ArcGIS Pro using this projection with very good results. The images were matched against the City of Toronto's Centreline dataset found here The resulting TIF files and their supporting spatial files include: TOLandUseMap13.tfwx TOLandUseMap13.tif TOLandUseMap13.tif.aux.xml TOLandUseMap13.tif.ovr TOLandUseMap14.tfwx TOLandUseMap14.tif TOLandUseMap14.tif.aux.xml TOLandUseMap14.tif.ovr TOLandUseMap15.tfwx TOLandUseMap15.tif TOLandUseMap15.tif.aux.xml TOLandUseMap15.tif.ovr TOLandUseMap16.tfwx TOLandUseMap16.tif TOLandUseMap16.tif.aux.xml TOLandUseMap16.tif.ovr TOLandUseMap17.tfwx TOLandUseMap17.tif TOLandUseMap17.tif.aux.xml TOLandUseMap17.tif.ovr TOLandUseMap18.tfwx TOLandUseMap18.tif TOLandUseMap18.tif.aux.xml TOLandUseMap18.tif.ovr TOLandUseMap19.tif TOLandUseMap19.tif.aux.xml TOLandUseMap19.tif.ovr TOLandUseMap20.tfwx TOLandUseMap20.tif TOLandUseMap20.tif.aux.xml TOLandUseMap20.tif.ovr TOLandUseMap21.tfwx TOLandUseMap21.tif TOLandUseMap21.tif.aux.xml TOLandUseMap21.tif.ovr TOLandUseMap22.tfwx TOLandUseMap22.tif TOLandUseMap22.tif.aux.xml TOLandUseMap22.tif.ovr TOLandUseMap23.tfwx TOLandUseMap23.tif TOLandUseMap23.tif.aux.xml TOLandUseMap23.tif.ov Ground control points were saved for all georeferenced images. The files are the following: map13.txt map14.txt map15.txt map16.txt map17.txt map18.txt map19.txt map21.txt map22.txt map23.txt The City of Toronto's Property Boundaries shapefile, "property_bnds_gcc_wgs84.zip" were unzipped and also reprojected to EPSG 26917 (UTM Zone 17) into a new shapefile, "Property_Boundaries_UTM.shp" Mosaicing Images Once georeferenced, all images were then mosaiced into one image file, "LandUseMosaic20211220v01", within the project-generated Geodatabase, "Landuse.gdb" and exported TIF, "LandUseMosaic20211220.tif" Reclassifying Images Because the original images were of low quality and the conversion to TIF made the image colours even more inconsistent, a method was required to reclassify the images so that different land use classes could be identified. Using Deep learning Objects, the images were re-classified into useful consistent colours. Deep Learning Objects and Training The resulting mosaic was then prepared for reclassification using the Label Objects for Deep Learning tool in ArcGIS Pro. A training sample, "LandUseTrainingSamples20211220", was created in the geodatabase for all land use types as follows: Neighbourhoods Insitutional Natural Areas Core Employment Areas Mixed Use Areas Apartment Neighbourhoods Parks Roads Utility Corridors Other Open Spaces General Employment Areas Regeneration Areas Lettering (not a land use type, but an image colour (black), used to label streets). By identifying the letters, it then made the reclassification and vectorization results easier to clean up of unnecessary clutter caused by the labels of streets. Reclassification Once the training samples were created and saved, the raster was then reclassified using the Image Classification Wizard tool in ArcGIS Pro, using the Support...

This land cover data set is derived from the original raster based Globcover regional (Africa) archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (46 classes). This database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code)

You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend tables (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based Globcover database (regional version). Globcover has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.1 skim (10 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (46 classes); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Download - Land cover of United Republic of Tanzania - Shape file format

GLOBCOVER on the ESA Web site

Global Land Cover Network - GLCN

This land cover data set is derived from the original raster based Globcover regional (Africa) archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (46 classes). This database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code)

You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend tables (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based Globcover database (regional version). Globcover has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.1 skim (10 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (46 classes); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Download - Land cover of Mozambique - Shape file format

GLOBCOVER on the ESA Web site

Global Land Cover Network - GLCN

This land cover data set is derived from the original raster based Globcover regional (Africa) archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (46 classes). This database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code)

You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend tables (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based Globcover database (regional version). Globcover has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.1 skim (10 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (46 classes); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Download - Land cover of Seychelles - Shape file format

GLOBCOVER on the ESA Web site

Global Land Cover Network - GLCN

DamaGIS is a GIS database which aims to collect and assess flood-related damage data at the local scale. What sets this database apart is the type of sources it uses. Indeed, all types of sources are considered such as on-site observations, online media or even social networks. Moreover, the database aims to be filled out by a large number of contributors for the sake of data completeness and accuracy.

The reason for creating this database was the lack of precise damage data available to calibrate and validate flood risk assessment models. To this end, DamaGIS offers highly precise and easily accessible flood-related damage data. A simplified method to easily assess damage severity is also proposed to enable data comparison across time and location.

Since 2011, 729 damage entries caused by 23 flood events in the South of France have been reported within the database. The geodatabase contains three feature classes and two relationship classes which join those feature classes.

The first feature class is called “EVENT”. It is a table with a shape field containing polygon geometries for geographic features. The EVENT feature class identifies flood events that have caused damage in the South of France since 2011.

The second feature class is called “BASIN”. It is a table with a shape field containing polygon geometries for geographic features. It is related to French river basins which were extracted from the French Carthage database . Carthage joins a set of information layers representing geographic objects related to hydrology. It is available online as a free download (through SANDRE http://www.sandre.eaufrance.fr/).

The third and main feature class of the structure is the DAMAGE database, which catalogues flood-related damage. It is a table with a shape field containing point geometries for geographic features.

This land cover data set is derived from the original raster based Globcover regional (Africa) archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (46 classes). This database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code)

You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend tables (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based Globcover database (regional version). Globcover has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.1 skim (10 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (46 classes); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Download - Land cover of Senegal - Shape file format

GLOBCOVER on the ESA Web site

This land cover data set is derived from the original raster based Globcover global archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (22 classes worldwide). The database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -Perimeter (m) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code)

You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend table (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based global archive, Globcover. Globcover database has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.1 skim (10 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (24 classes globally); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Download Land cover of Viet Nam - Shape file format

GLOBCOVER on the ESA Web site

Global Land Cover Network - GLCN

This land cover data set is derived from the original raster based Globcover global archive. It has been post-processed to generate a vector version at national extent with the LCCS regional legend (22 classes worldwide). The database can be analyzed in the GLCN software Advanced Database Gateway (ADG), which provides a user-friendly interface and advanced functionalities to breakdown the LCCS classes in their classifiers for further aggregations and analysis.

The data set is intended for free public access.

The shape file's attributes contain the following fields: -Area (sqm) -Perimeter (m) -ID -Gridcode (Globcover cell value) -LCCCode (unique LCCS code) You can download a zip archive containing: -the shape file (.shp) -the ArcGis layer file with global legend (.lyr) -the ArcView 3 legend file (.avl) -the LCCS legend table (.xls)

Supplemental Information:

This land cover product is a vector version (ESRI shape) of the Globcover archive that was published in 2008 as result of an initiative launched in 2004 by the European Space Agency (ESA). Globcover is currently the most recent (2005) and resoluted (300 m) datasets on land cover globally. Given the need of this valuable information for environmental studies, natural resources management and policy formulation, through activities of the Global Land Cover Network (GLCN) programme, the Globcover has been reprocessed to generate databases at national extent that can be analyzed through the Advanced Database Gateway software (ADG) by GLCN. ADG is a cross-cutting interrogation software that allows the easy and fast recombination of land cover polygons according to the individual end-user requirements. Aggregated land cover classes can be generated not only by name, but also using the set of existing classifiers. ADG uses land cover data with a Land Cover Classification System (LCCS) legend. The ADG software is available for download on the GLCN web site at http://www.glcn.org/sof_7_en.jsp

Contact points:

Metadata Contact: FAO-Data

Resource Contact: Antonio Martucci

Data lineage:

This land cover database is provided as ESRI shape file (vector format) and derives from reprocessing the raster based global archive, Globcover. Globcover database has undergone the following process: a) vectoralization at the national extent using ESRI ArcGis (arcinfo) 9.3; b) topological reconstruction (custom AML scripts launched inside ArcGis-arcinfo 9.3); c) simplification of areas according to a minimum mapping unit of 0.2 skim (20 ha) (custom AML scripts launched inside ArcGis-arcinfo 9.3); application of the FAO/UNEP Land Cover Classification System (LCCS) legend (24 classes globally); final processing to assure full compatibility with the GLCN software Advanced Database Gateway (ADG).

Online resources:

Land cover of Democratic Republic of the Congo - Shape file format

GLOBCOVER on the ESA Web site

Global Land Cover Network - GLCN

CanVec contains more than 60 topographic features classes organized into 8 themes: Transport Features, Administrative Features, Hydro Features, Land Features, Manmade Features, Elevation Features, Resource Management Features and Toponymic Features. This multiscale product originates from the best available geospatial data sources covering Canadian territory. It offers quality topographic information in vector format complying with international geomatics standards. CanVec can be used in Web Map Services (WMS) and geographic information systems (GIS) applications and used to produce thematic maps. Because of its many attributes, CanVec allows for extensive spatial analysis. Related Products: Constructions and Land Use in Canada - CanVec Series - Manmade Features Lakes, Rivers and Glaciers in Canada - CanVec Series - Hydrographic Features Administrative Boundaries in Canada - CanVec Series - Administrative Features Mines, Energy and Communication Networks in Canada - CanVec Series - Resources Management Features Wooded Areas, Saturated Soils and Landscape in Canada - CanVec Series - Land Features Transport Networks in Canada - CanVec Series - Transport Features Elevation in Canada - CanVec Series - Elevation Features Map Labels - CanVec Series - Toponymic Features

Benchmark set at 77.1% O.A at: https://doi.org/10.1117/1.JRS.14.048503

The dataset consists of 60,000 images, corresponding to Landsat patches of 33x33 pixels with 102 bands. Randomly selected from Mexico (country). Each patch is labeled with one of 12 Land Use and Vegetation classes according to the classification described at https://doi.org/10.3390/rs6053923.

The zip file contains 12 folders numbered 1-12 and each contains 5,000 .npy python files (can be loaded with the NumPy library).

The labeled classes correspond to the following identifier.

1, Temperate Coniferous forest
2, Temperate Decidius Forest
3, Temperate Mixed Forest
4, Tropical Evergreen Forest
5, Tropical Deciduous Forest
6, Scrubland
7, Wetland Vegetation
8, Agriculture
9, Grassland
10, Water body
11, Barren Land
12, Urban Area

To build that dataset, we take the information of the National Continuum of Land Use and Vegetation series number 5 generated by the National Institute of Statistics and Geography from Mexico (INEGI) from The National Commission for the Knowledge and Use of Biodiversity (CONABIO) web page (http://geoportal.conabio.gob.mx/metadatos/doc/html/usv250s5ugw.html).

The file used for this dataset construction is the shape format file with geographic coordinates located in http://www.conabio.gob.mx/informacion/gis/maps/geo/usv250s5ugw.zip.
Later, a transformation to Albers equal-area conic projection was done with the followings parameters:

Fake east: 2500000.0
Fake North: 0.0
Origin longitude: -102.0º
Origin latitude: 12.0º
First standard parallel: 17.5º
Second standard parallel: 29.5º
Linear unit: Meter (1.0)
Reference ellipsoid: GRS80

Once the data was projected, using the classes identified in the National Continuum of Land Use and Vegetation, correspondence was applied to the classes identified in https://doi.org/10.3390/rs6053923, these classes being: Agriculture, Barren land, Grassland, Scrubland, Temperate coniferous forest, Temperate deciduous forest, Temperate mixed forest, Tropical deciduous forest, Tropical evergreen forest, Urban area, Waterbody and Wetland vegetation.

Once the information layer was generated with the 12 classes indicated above, the reference layer was rasterized.
Thus, a national grid of 1,975,940 regions of 1 x 1 kilometers was generated and the percentage of pixels of the dominant class in each corresponding 1 km region was associated.

A total of cells with 70% or more pixels from one dominant class corresponds to 1,640,827 which represents a total of 83% of the Mexican territory. That means, only 17% of cells have less than 70% of their pixels from one dominant class.
Then, 5000 regions were randomly selected from each land cover class at the national level. For this random selection only were selected the regions in which cells have 70% or more of their pixels from one dominant class. The above, for looking to have consistent and reliable data for the automatic classification task. This random selection generates a total of 60,000 regions selected.

Image patches were extracted from the selected regions in the sample.

The image used is the result of the application of multiple time series analysis algorithms on a cube of image data with mainly Tier 1 (T1) quality and a few Tier 2 (T2) as described in https: // www. usgs.gov/land-resources/nli/landsat/landsat-collection-1. An Open Data Cube (ODC, https://www.opendatacube.org/) was constructed from 3,515 Landsat 5 and 7 images corresponding to the year 2011, which is the same reference year of the National Continuum of Land Use and Vegetation Series 5.

From the analysis of the ODC images, the Geomedian (https://doi.org/10.1109/TGRS.2017.2723896) was calculated, which generated a national cloud-free mosaic from 2011, pixels at 30 meters resolution and 6 spectral bands (blue, green, red, nir, swir 1, swir 2). Finally, 15 spectral indices were calculated for each pixel in the image. This resulted in 15 national mosaics from the analysis of the time series of each pixel available for the year 2011 using all the combinations of normalized difference indices, which were possible with the 6 bands that were incorporated into the data cube, with which resulted in 102 information channels. Since Landsat images have a resolution of 30 meters, we have images of 33 pixels x 33 pixels for each region of 1 km x 1 km.

The 102 channels in the patches correspond to:

Geomedian Bands (6): blue, green, red, nir, swir 1, swir 2
Geomedian Based Indexes (15): evi, bu, sr, arvi, ui, ndbi, ibi, ndvi, ndwi, mndwi, nbi, brba, nbai, baei, bi
Geomedian Based Tasseled cap transformation (6): brightness, greenness, wetness, fourth, fifth, sixth

2011 Landsat Time Analysis Series by Pixel

(red-swir 1)/(red+swir 1); (5): min, mean, max, std, median
(red-nir)/( red+nir); (5): min, mean, max, std, median
(swir 1-swir 2)/( swir 1+swir 2); (5): min, mean, max, std, median
(nir-swir 2)/(nir+swir 2); (5): min, mean, max, std, median
(nir-swir 1)/( nir+swir 1); (5): min, mean, max, std, median
(red-swir 2)/( red+swir 2); (5): min, mean, max, std, median
(green-swir 2)/(green+swir 2); (5): min, mean, max, std, median
(green-swir 1)/(green+swir 1); (5): min, mean, max, std, median
(green-red)/(green+red); (5): min, mean, max, std, median
(green-nir)/(green+nir); (5): min, mean, max, std, median
(blue-swir 2)/(blue+swir 2); (5): min, mean, max, std, median
(blue-swir 1)/(blue+swir 1); (5): min, mean, max, std, median
(blue-red)/(blue+red); (5): min, mean, max, std, median
(blue-nir)/(blue+nir); (5): min, mean, max, std, median
(blue-green)/( blue+green); (5): min, mean, max, std, median